Time-course RNA sequencing reveals high similarity in mRNAome between hepatic stellate cells activated by agalactosyl IgG and TGF-β1

Previous studies have demonstrated the clinical relevance of aberrant serum immunoglobulin G (IgG) N-glycomic profiles in liver fibrosis and the pathogenic effects of agalactosyl IgG on activating hepatic stellate cells (HSCs). However, the dynamics of gene expression changes during HSC activation by agalactosyl IgG remain poorly understood. We performed RNA sequencing to analyze the mRNAome of human LX-2 HSCs at multiple time points after treatment with agalactosyl IgG and then compared these results with those obtained after normal IgG and transforming growth factor (TGF)-β1 treatments. Gene expression changes were significantly pronounced on day 5 and subsided by day 11 after HSC activation. A high degree of similarity in gene expression patterns between HSCs treated with agalactosyl IgG and TGF-β1 was observed, of which 1796 and 1785 differentially expressed genes (DEGs) were identified, respectively. Disease ontology analyses revealed that 114 and 105 DEGs in activated HSCs following agalactosyl IgG and TGF-β1 treatments, respectively, were linked to liver cirrhosis, hepatitis, fatty liver disease, hepatitis B, and alcoholic hepatitis, with CCL5 and FAS being the most commonly affected genes. DEGs associated with liver fibrosis or aforementioned liver diseases involved in gene annotation, physiological functions, and signaling pathways regarding secretion of cytokines and chemokines, expression of fibrosis-related growth factors and their receptors, modification of extracellular matrices, and regulation of cell viability in activated HSCs. In conclusion, this study characterized the dynamics of mRNAome and gene networks and identified the liver fibrosis-related DEGs during HSC activation by agalactosyl IgG and TGF-β1.