Identification of core gene networks and hub genes associated with progression of non-alcoholic fatty liver disease by RNA sequencing

Bioinformatics analysis and validation of genes related to paclitaxel's anti-breast cancer effect through immunogenic cell death

Research indicated that Paclitaxel (PTX) can induce immunogenic cell death (ICD) through immunogenic modulation. However, the combination of PTX and ICD has not been extensively studied in breast cancer (BRCA). The TCGA-BRCA and GSE20685 datasets were enrolled in this study. Samples from the TCGA-BRCA dataset were consistently clustered based on selected immunogenic cell death-related genes (ICD-RGs). Next, candidate genes were obtained by overlapping differentially expressed genes (DEGs) between BRCA and normal groups, intersecting genes common to DEGs between cluster1 and cluster2 and hub module genes, and target genes of PTX from five databases. The univariate Cox algorithm and the least absolute shrinkage and selection operator (LASSO) were performed to obtain biomarkers and build a risk model. Following observing the immune microenvironment in differential risk subgroups, single-gene gene set enrichment analysis (GSEA) was carried out in all biomarkers. Finally, the expression of biomarkers was analyzed. Enrichment analysis showed that 626 intersecting genes were linked with inflammatory response. Further five biomarkers (CHI3L1, IL18, PAPLN, SH2D2A, and UBE2L6) were identified and a risk model was built. The model's performance was validated using GSE20685 dataset. Furthermore, the biomarkers were enriched with adaptive immune response. Lastly, the experimental results indicated that the alterations in IL18, SH2D2A, and CHI3L1 expression after treatment matched those in the public database. In this study, Five PTX-ICD-related biomarkers (CHI3L1, IL18, PAPLN, SH2D2A, and UBE2L6) were identified to aid in predicting BRCA treatment outcomes.

Emerging roles for the ADAMTS-like family of matricellular proteins in cardiovascular disease through regulation of the extracellular microenvironment

Dysregulation of the extracellular matrix (ECM) occurs widely across cardiovascular pathologies. Recent work has revealed important roles for the «a disintegrin-like and metalloprotease domain with thrombospondin-type 1 motifs like" (ADAMTSL) family of secreted glycoproteins in cardiovascular tissues during development and disease. Key insights in this regard have come from naturally occurring gene mutations in humans and animals that result in severe diseases with cardiovascular manifestations or aortopathies. Expression of ADAMTSL genes is greatly increased in the myocardium during heart failure. Genetically modified mice recapitulate phenotypes of patients with ADAMTSL mutations and demonstrate important functions in the ECM. The novel functions thus disclosed are intriguing because, while these proteins are neither structural, nor proteases like the related ADAMTS proteases, they appear to act as regulatory, i.e., matricellular proteins. Evidence from genetic variants, genetically engineered mouse mutants, and in vitro investigations have revealed regulatory functions of ADAMTSLs related to fibrillin microfibrils and growth factor signaling. Interestingly, the ability to regulate transforming growth factor (TGF)β signaling may be a shared characteristic of some ADAMTSLs. TGFβ signaling is important in cardiovascular development, health and disease and a central driver of ECM remodeling and cardiac fibrosis. New strategies to target dysregulated TGFβ signaling are warranted in aortopathies and cardiac fibrosis. With their emerging roles in cardiovascular tissues, the ADAMTSL proteins may provide causative genes, diagnostic biomarkers and novel treatment targets in cardiovascular disease. Here, we discuss the relevance of ADAMTSLs to cardiovascular medicine.

Identification of differentially methylated regions associated with both liver fibrosis and hepatocellular carcinoma

BACKGROUND

Liver fibrosis is a major risk factor for hepatocellular carcinoma (HCC). We have previously reported that differentially methylated regions (DMRs) are correlated with the fibrosis stages of metabolic dysfunction-associated steatotic liver disease (MASLD). In this study, the methylation levels of those DMRs in liver fibrosis and subsequent HCC were examined.

METHODS

The methylation levels of DMRs were investigated using alcoholic cirrhosis and HCC (GSE60753). The data of hepatitis C virus-infected cirrhosis and HCC (GSE60753), and two datasets (GSE56588 and GSE89852) were used for replication analyses. The transcriptional analyses were performed using GSE114564, GSE94660, and GSE142530.

RESULTS

Hypomethylated DMR and increased transcriptional level of zinc finger and BTB domain containing 38 (ZBTB38) were observed in HCC. Hypermethylated DMRs, and increased transcriptional levels of forkhead box K1 (FOXK1) and zinc finger CCCH-type containing 3 (ZC3H3) were observed in HCC. The methylation levels of DMR of kazrin, periplakin interacting protein (KAZN) and its expression levels were gradually decreased as cirrhosis progressed to HCC.

CONCLUSIONS

Changes in the methylation and transcriptional levels of ZBTB38, ZC3H3, FOXK1, and KAZN are important for the development of fibrosis and HCC; and are therefore potential therapeutic targets and diagnostic tools for cirrhosis and HCC.

Dietary Long-Chain Fatty Acids Accelerate Metabolic Dysfunction in Guinea Pigs with Non-Alcoholic Steatohepatitis

The composition of dietary fatty acids may be important for the development and progression of metabolic syndrome and non-alcoholic steatohepatitis (NASH). This study investigated the effect of two high-fat diets based on coconut oil, containing predominantly medium-chain fatty acids (MCFA), or cocoa butter, containing mainly long-chain fatty acids (LCFA), on glucose homeostasis and NASH in guinea pigs following 16 and 32 weeks of diet. At week 16, glucose intolerance was increased in the LCFA animals compared to the MCFA animals (p < 0.001), with both groups differing from the controls by week 32 (p < 0.0001), supported by increased hemoglobin A1c (p < 0.05). NASH was present in both high-fat groups from week 16, with advancing fibrosis appearing more progressive in the LCFA animals at week 16. In agreement, gene expression showed overall increased expression of NASH target genes in the LCFA animals compared to the MCFA animals at weeks 16 and 32 (p < 0.05 and p < 0.0001, respectively). The LCFA animals also displayed increased plasma uric acid at both time points (p < 0.05), a phenomenon linked to NASH in humans. In conclusion, this study reports that a diet high in LCFA promotes metabolic imbalance and may accelerate NASH-associated hepatic fibrosis. This highlights the importance of a critical evaluation of fatty acid composition when investigating NASH-associated endpoints.

Bioinformatics analysis reveals molecular connections between non-alcoholic fatty liver disease (NAFLD) and COVID-19

The ongoing coronavirus disease 2019 (COVID-19) pandemic caused by SARS-CoV-2 has devastatingly impacted people's lives. Non-alcoholic fatty liver disease (NAFLD) is fatal comorbidity of COVID-19 seen with potential risk factors to develop severe symptoms. This research focuses on determining and elucidating the molecular factors and connections that might contribute to the severity of SARS-CoV-2 infection in NAFLD patients. Here, we comprehensively inspected the genes involved in NAFLD and SARS-CoV-2 entry factors (SCEFs) found by searching through the DisGeNet database and literature review, respectively. Further, we identified the SCEFs-related proteins through protein-protein interaction (PPI) network construction, MCODE, and Cytohubba. Next, the shared genes involved in NAFLD and SARS-CoV-2 entry, and hub gene were determined, followed by the GO and KEGG pathways analysis. X2K database was used to construct the upstream regulatory network of hub genes, as well as to identify the top ten candidates of transcription factors (TFs) and protein kinases (PKs). PPI analysis identified connections between 4 top SCEFs, including ACE, ADAM17, DPP4, and TMPRSS2 and NAFLD-related genes such as ACE, DPP4, IL-10, TNF, and AKT1. GO and KEGG analysis revealed the top ten biological processes and pathways, including cytokine-mediated signaling, PI3K-Akt, AMPK, and mTOR signaling pathways. The upstream regulatory network revealed that AKT1 and MAPK14 as important PKs and HIF1A and SP1 as important TFs associated with AKT1, IL-10, and TNF. The molecular connections identified between COVID-19 and NAFLD may shed light on discovering the causes of the severity of SARS-CoV-2 infected NAFLD patients.

Two differentially methylated region networks in nonalcoholic fatty liver disease, viral hepatitis, and hepatocellular carcinoma

Background

We previously reported that two differentially methylated region (DMR) networks identified by DMR and co-methylation analyses are strongly correlated with the fibrosis stages of nonalcoholic fatty liver disease (NAFLD). In the current study, we examined these DMR networks in viral hepatitis and hepatocellular carcinoma (HCC).

Methods

We performed co-methylation analysis of DMRs using a normal dataset (GSE48325), two NAFLD datasets (JGAS000059 and GSE31803), and two HCC datasets (GSE89852 and GSE56588). The dataset GSE60753 was used for validation.

Results

One DMR network was clearly observed in viral hepatitis and two HCC populations. Methylation levels of genes in this network were higher in viral hepatitis and cirrhosis, and lower in HCC. Fatty acid binding protein 1 (FABP1), serum/glucocorticoid regulated kinase 2 (SGK2), and hepatocyte nuclear factor 4 α (HNF4A) were potential hub genes in this network. Increased methylation levels of the FABP1 gene may be correlated with reduced protection of hepatocytes from oxidative metabolites in NAFLD and viral hepatitis. The decreased methylation levels of SGK2 may facilitate the growth and proliferation of HCC cells. Decreased methylation levels of HNF4A in HCC may be associated with tumorigenesis. The other DMR network was observed in NAFLD, but not in viral hepatitis or HCC. This second network included genes involved in transcriptional regulation, cytoskeleton organization, and cellular proliferation, which are specifically related to fibrosis and/or tumorigenesis in NAFLD.

Conclusions

Our results suggest that one DMR network was associated with fibrosis and tumorigenesis in both NAFLD and viral hepatitis, while the other network was specifically associated with NAFLD progression. Furthermore, FABP1, SGK2, and HNF4A are potential candidate targets for the prevention and treatment of HCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12876-022-02360-4.

Ficolin-2 Plasma Level Assesses Liver Fibrosis in Non-Alcoholic Fatty Liver Disease

Fibrosis is the strongest predictor for disease-specific mortality in non-alcoholic fatty liver diseases (NAFLD), but the need for liver biopsy limits its diagnosis. We assessed the performance of plasma ficolin-2 (FCN-2) as a biomarker of fibrosis identified by an in silico discovery strategy. Two hundred and thirty-five morbidly obese (MO) subjects with biopsy-proven NAFLD stratified by fibrosis stage (F0, n = 44; F1, n = 134; F2, n = 46; F3/F4, n = 11) and 40 cirrhotic patients were enrolled. The cohort was subdivided into discovery (n = 76) and validation groups (n = 159). The plasma level of FCN-2 and other candidate markers was determined. FCN-2 was inversely correlated with the stage of liver fibrosis (ρ = −0.49, p < 0.001) independently of steatosis (p = 0.90), inflammation (p = 0.57), and ballooning (p = 0.59). In the global cohort, FCN-2 level decreased significantly in a stepwise fashion from F0/F1 (median 4753 ng/mL) to F2−F3−F4 (2760 ng/mL) and in cirrhotic subjects (1418 ng/mL). The diagnostic performance of FCN-2 in detecting F ≥ 2 was higher than other indexes (APRI, FIB-4) (AUROC 0.82, 0.68, and 0.6, respectively). The accuracy improved when combined with APRI score and HDL values (FCNscore, AUROC 0.85). Overall, the FCN-2 plasma level can accurately discriminate liver fibrosis status (minimal vs. moderate/advanced) significantly improving the fibrosis diagnostic algorithms.

The molecular phenotypes of injury, steatohepatitis, and fibrosis in liver transplant biopsies in the INTERLIVER study

To extend previous molecular analyses of rejection in liver transplant biopsies in the INTERLIVER study (ClinicalTrials.gov #NCT03193151), the present study aimed to define the gene expression selective for parenchymal injury, fibrosis, and steatohepatitis. We analyzed genome-wide microarray measurements from 337 liver transplant biopsies from 13 centers. We examined expression of genes previously annotated as increased in injury and fibrosis using principal component analysis (PCA). PC1 reflected parenchymal injury and related inflammation in the early posttransplant period, slowly regressing over many months. PC2 separated early injury from late fibrosis. Positive PC3 identified a distinct mildly inflamed state correlating with histologic steatohepatitis. Injury PCs correlated with liver function and histologic abnormalities. A classifier trained on histologic steatohepatitis predicted histologic steatohepatitis with cross-validated AUC = 0.83, and was associated with pathways reflecting metabolic abnormalities distinct from fibrosis. PC2 predicted histologic fibrosis (AUC = 0.80), as did a molecular fibrosis classifier (AUC = 0.74). The fibrosis classifier correlated with matrix remodeling pathways with minimal overlap with those selective for steatohepatitis, although some biopsies had both. Genome-wide assessment of liver transplant biopsies can not only detect molecular changes induced by rejection but also those correlating with parenchymal injury, steatohepatitis, and fibrosis, offering potential insights into disease mechanisms for primary diseases.

Bile acid dysmetabolism in the gut-microbiota-liver axis under hepatitis C virus infection

BACKGROUND & AIMS

We recently analysed and reported the features of the micro biome under hepatitis C virus (HCV) infection, but the effect of HCV infection on bile acid (BA) metabolism in the gut-liver axis remains poorly understood. The aim of this study was to clarify the characteristics of the gut-liver axis in HCV-infected patients.

METHODS

The faecal BAs composition and gut microbiota from 100 chronic hepatitis C (CHC) patients were compared with those from 23 healthy individuals. For transcriptional analysis of the liver, 22 mild CHC (fibrosis stages [F] 0-2) and 42 advanced CHC (F3-4) cases were compared with 12 healthy individuals. The findings were confirmed using chimeric mice with human hepatocytes infected with HCV HCR6.

RESULTS

Chronic hepatitis C patients, even at earlier disease stages, showed BA profiles distinct from healthy individuals, in which faecal deoxycholic acid (DCA) was significantly reduced and lithocholic acid or ursodeoxycholic acid became dominant. The decrease in faecal DCA was correlated with reduction in commensal Clostridiales and increase in oral Lactobacillales. Impaired biosynthesis of cholic acid (CA) was observed as a reduction in the transcription level of cytochrome P450 8B1 (CYP8B1), a key enzyme in CA biosynthesis. The reductions in faecal DCA and liver CYP8B1 were also observed in HCV-infected chimeric mice.

CONCLUSIONS

Chronic hepatitis C alters the intestinal BA profile, in association with the imbalance of BA biosynthesis, which differs from the pattern in NAFLD. These imbalances appear to drive disease progression through the gut-microbiome-liver axis.

A joint analysis using exome and transcriptome data identifiescandidate polymorphisms and genes involved with umbilical hernia in pigs

BACKGROUND

Umbilical Hernia (UH) is characterized by the passage of part of the intestine through the umbilical canal forming the herniary sac. There are several potential causes that can lead to the umbilical hernia such as bacterial infections, management conditions and genetic factors. Since the genetic components involved with UH are poorly understood, this study aimed to identify polymorphisms and genes associated with the manifestation of umbilical hernia in pigs using exome and transcriptome sequencing in a case and control design.

RESULTS

In the exome sequencing, 119 variants located in 58 genes were identified differing between normal and UH-affected pigs, and in the umbilical ring transcriptome, 46 variants were identified, located in 27 genes. Comparing the two methodologies, we obtained 34 concordant variants between the exome and transcriptome analyses, which were located in 17 genes, distributed in 64 biological processes (BP). Among the BP involved with UH it is possible to highlight cell adhesion, cell junction regulation, embryonic morphogenesis, ion transport, muscle contraction, within others.

CONCLUSIONS

We have generated the first exome sequencing related to normal and umbilical hernia-affected pigs, which allowed us to identify several variants possibly involved with this disorder. Many of those variants present in the DNA were confirmed with the RNA-Seq results. The combination of both exome and transcriptome sequencing approaches allowed us to better understand the complex molecular mechanisms underlying UH in pigs and possibly in other mammals, including humans. Some variants found in genes and other regulatory regions are highlighted as strong candidates to the development of UH in pigs and should be further investigated.

Deciphering the role of epigenetic modifications in fatty liver disease: A systematic review

BACKGROUND

Fatty liver disease (FLD), primarily nonalcoholic fatty liver disease (NAFLD), is the most common liver disorder that affects a quarter of the global population. NAFLD is a spectrum of disease ranging from simple steatosis to nonalcoholic steatohepatitis, which is associated with increased risk of developing liver cancer. Given that the pathogenic mechanisms of fatty liver remain largely elusive, it is important to further investigate potential underlying mechanisms including epigenetic modifications. Here, we performed a systematic review of human epigenetic studies on FLD presence.

METHODS

Five bibliographic databases were screened until 28 August 2020. We included cross-sectional, case-control and cohort studies in humans that examined the association of epigenetic modifications including global, candidate or epigenome-wide methylation of DNA, noncoding RNAs and histone modifications with FLD.

RESULTS

In total 36 articles, based on 33 unique studies, consisting of 12 112 participants met the inclusion criteria. Among these, two recent epigenome-wide association studies conducted among large population-based cohorts have reported the association between cg06690548 (SLC7A11) and FLD. Moreover, several studies have demonstrated the association between microRNAs (miRNAs) and FLD, in which miR-122, miR-34a and miR-192 were recognized as the most relevant miRNAs as biomarkers for FLD. We did not find any studies examining histone modifications in relation to FLD.

CONCLUSIONS

Cumulative evidence suggests a link between epigenetic mechanisms, specifically DNA methylation and miRNAs, and FLD. Further efforts should investigate the molecular pathways by which these epigenetic markers may regulate FLD and also the potential role of histone modifications in FLD.

Serum thrombospondin 2 is a novel predictor for the severity in the patients with NAFLD

AIM

Thrombospondins are a family of multidomain and secretory glycoproteins. Among them, thrombospondin 2 (TSP2) encoded by TSP2 gene has been reported to be involved in various functions such as collagen/fibrin formation, maintenance of normal blood vessel density and cell adhesion properties. Microarray analyses ranked TSP2 as one of the most highly up-regulated genes in the fibrotic liver in patients with non-alcoholic fatty liver disease (NAFLD). Since TSP2 possesses unique properties as a secretory protein, we hypothesized that hepatic TSP2 gene expression levels would be reflected in serum TSP2 levels. In this study, we examined the relationship between serum TSP2 concentrations and clinicopathological findings in NAFLD patients.

METHODS

One hundred and thirty NAFLD patients who had undergone liver biopsy between 2009 and 2015 were retrospectively enrolled. Serum samples were collected at the time of biopsy, and TSP2 was measured by enzyme immunoassays.

RESULTS

Serum TSP2 levels moderately correlated with ballooning (r = 0.56, P < .001) and fibrosis stage (r = 0.53, P < .001). The AUC values of TSP2 for predicting mild fibrosis (≧F1), moderate fibrosis (≧F2) and severe fibrosis (≧F3) were 0.73, 0.76 and 0.82 respectively. Additionally, NAFLD activity score (NAS) correlated best with TSP2 (r = 0.52, P < .001) compared to conventional NAFLD-related biomarkers, such as cytokeratin 18 M30, hyaluronic acid, type IV collagen 7S, APRI and FIB-4 index.

CONCLUSION

Serum TSP2 levels reflected hepatocyte ballooning, fibrosis and NAS in NAFLD patients. For clinical application of serum TSP2 as a predictor of NAFLD histological activity, additional validation and mechanistic investigations are required.

Weighted gene coexpression network analysis reveals hub genes involved in cholangiocarcinoma progression and prognosis

AIM

Cholangiocarcinoma (CCA) is a highly malignant tumor found in the bile duct epithelial cells, and the second most common primary tumor of the liver. However, the pivotal roles of molecular biomarkers in oncogenesis of CCA are unclear. Therefore, we aim to explore the underlying mechanisms of progression and screen for novel prognostic biomarkers and treatment targets.

METHOD

The data of mRNA sequencing and clinical information of CCA patients in The Cancer Genome Atlas was analyzed by weighted gene coexpression network analysis (WGCNA). Modules and clinical traits were constructed according to Pearson's correlation analysis, and Gene Ontology and pathway enrichment analysis were applied. Hub genes of these modules were screened by intramodule analysis; Cytoscape with Search Tool for the Retrieval of Interacting Genes was utilized to visualize protein-protein interaction of these modules; hub genes of these modules were validated afterwards. Furthermore, the significance of these genes was confirmed by survival analysis.

RESULTS

Genes MRPS18A, CST1, and SCP2 were identified as candidate genes in the module, which was associated with clinical traits including pathological stage, histological grade, and liver function and which also affected overall survival of CCA patients. Nineteen hub genes were analyzed together and were associated with progression and prognosis of CCA. Survival analyses found that several of the multiple genes could serve as biomarkers to stratify CCA patients into low- and high-risk groups.

CONCLUSION

These candidate genes could be involved in progression of CCA, which could serve as novel prognostic markers and treatment targets. Moreover, most of them were first reported in CCA and deserve further research.

Multitranscriptome analyses reveal prioritized genes specifically associated with liver fibrosis progression independent of etiology

Elimination or suppression of causative factors can raise the possibility of liver fibrosis regression. However, different injurious stimuli will give fibrosis from somewhat different etiologies, which, in turn, may hamper the discovery of liver fibrosis-specific therapeutic drugs. Therefore, the analogical cellular and molecular events shared by various etiology-evoked liver fibrosis should be clarified. Our present study systematically integrated five publicly available transcriptomic data sets regarding liver fibrosis with different etiologies from the Gene Expression Omnibus database and performed a series of bioinformatics analyses and experimental verifications. A total of 111 significantly upregulated and 16 downregulated genes were identified specific to liver fibrosis independent of any etiology. These genes were predominately enriched in some Kyoto Encyclopedia of Genes and Genomes pathways, including the "PI3K-AKT signaling pathway," "Focal adhesion," and "ECM-receptor interaction." Subsequently, five prioritized liver fibrosis-specific genes, including COL4A2, THBS2, ITGAV, LAMB1, and PDGFRA, were screened. These genes were positively associated with each other and liver fibrosis progression. In addition, they could robustly separate all stages of samples in both training and validation data sets with diverse etiologies when they were regarded as observed variables applied to principal component analysis plots. Expressions of all five genes were confirmed in activated primary mouse hepatic stellate cells (HSCs) and transforming growth factor β1-treated LX-2 cells. Moreover, THBS2 protein was enhanced in liver fibrosis rodent models, which could promote HSC activation and proliferation and facilitate NOTCH1/JAG1 expression in HSCs. Overall, our current study may provide potential targets for liver fibrosis therapy and aid to a deeper understanding of the molecular underpinnings of liver fibrosis. NEW & NOTEWORTHY Prioritized liver fibrosis-specific genes THBS2, COL4A2, ITGAV, LAMB1, and PDGFRA were identified and significantly associated with liver fibrosis progression and could be combined to discriminate liver fibrosis stages regardless of any etiology. Among the identified prioritized liver fibrosis-specific targets, THBS2 protein was confirmed to be enhanced in liver fibrosis rodent models, which could promote hepatic stellate cell (HSC) activation and proliferation and facilitate NOTCH1/JAG1 expression in HSCs.

Hepatic transcriptome signatures in patients with varying degrees of nonalcoholic fatty liver disease compared with healthy normal-weight individuals

Nonalcoholic fatty liver disease (NAFLD) represents a spectrum of conditions ranging from simple steatosis (NAFL), over nonalcoholic steatohepatitis (NASH) with or without fibrosis, to cirrhosis with end-stage disease. The hepatic molecular events underlying the development of NAFLD and transition to NASH are poorly understood. The present study aimed to determine hepatic transcriptome dynamics in patients with NAFL or NASH compared with healthy normal-weight and obese individuals. RNA sequencing and quantitative histomorphometry of liver fat, inflammation and fibrosis were performed on liver biopsies obtained from healthy normal-weight ( n = 14) and obese ( n = 12) individuals, NAFL ( n = 15) and NASH ( n = 16) patients. Normal-weight and obese subjects showed normal liver histology and comparable gene expression profiles. Liver transcriptome signatures were largely overlapping in NAFL and NASH patients, however, clearly separated from healthy normal-weight and obese controls. Most marked pathway perturbations identified in both NAFL and NASH were associated with markers of lipid metabolism, immunomodulation, extracellular matrix remodeling, and cell cycle control. Interestingly, NASH patients with positive Sonic hedgehog hepatocyte staining showed distinct transcriptome and histomorphometric changes compared with NAFL. In conclusion, application of immunohistochemical markers of hepatocyte injury may serve as a more objective tool for distinguishing NASH from NAFL, facilitating improved resolution of hepatic molecular changes associated with progression of NAFLD. NEW & NOTEWORTHY Nonalcoholic fatty liver disease (NAFLD) is the most common liver disease in Western countries. NAFLD is associated with the metabolic syndrome and can progress to the more serious form, nonalcoholic steatohepatitis (NASH), and ultimately lead to irreversible liver damage. Using gold standard molecular and histological techniques, this study demonstrates that the currently used diagnostic tools are problematic for differentiating mild NAFLD from NASH and emphasizes the marked need for developing improved histological markers of NAFLD progression.

Identification of differentially methylated region (DMR) networks associated with progression of nonalcoholic fatty liver disease

The progression of nonalcoholic fatty liver disease (NAFLD) is affected by epigenetics. We performed differentially methylated region (DMR) and co-methylation analyses to identify DMR networks associated with the progression of NAFLD. DMRs displaying differences in multiple consecutive differentially methylated CpGs between mild and advanced NAFLD were extracted. The average values of topological overlap measures for the CpG matrix combining two different DMRs were calculated and two DMR networks that strongly correlated with the stages of fibrosis were identified. The annotated genes of one network included genes involved in transcriptional regulation, cytoskeleton organization, and cellular proliferation. The annotated genes of the second network were primarily associated with metabolic pathways. The CpG methylation levels in these networks were strongly affected by age and fasting plasma glucose levels, which may be important co-regulatory factors. The methylation status of five DMRs in the second network was reversible following weight loss. Our results suggest that CpG methylation in DMR networks is regulated concomitantly via aging and hyperglycemia and plays important roles in hepatic metabolic dysfunction, fibrosis, and potential tumorigenesis, which occur during the progression of NAFLD. By controlling weight and blood glucose levels, the methylation of DMRs in the second network may be reduced.

Hepatic expression profiling identifies steatosis-independent and steatosis-driven advanced fibrosis genes

Chronic liver disease (CLD) is associated with tissue-destructive fibrosis. Considering that common mechanisms drive fibrosis across etiologies, and that steatosis is an important cofactor for pathology, we performed RNA sequencing on liver biopsies of patients with different fibrosis stages, resulting from infection with hepatitis C virus (HCV) (with or without steatosis) or fatty liver disease. In combination with enhanced liver fibrosis score correlation analysis, we reveal a common set of genes associated with advanced fibrosis, as exemplified by those encoding the transcription factor ETS-homologous factor (EHF) and the extracellular matrix protein versican (VCAN). We identified 17 fibrosis-associated genes as candidate EHF targets and demonstrated that EHF regulates multiple fibrosis-associated genes, including VCAN, in hepatic stellate cells. Serum VCAN levels were also elevated in advanced fibrosis patients. Comparing biopsies from patients with HCV with or without steatosis, we identified a steatosis-enriched gene set associated with advanced fibrosis, validating follistatin-like protein 1 (FSTL1) as an exemplar of this profile. In patients with advanced fibrosis, serum FSTL1 levels were elevated in those with steatosis (versus those without). Liver Fstl1 mRNA levels were also elevated in murine CLD models. We thus reveal a common gene signature for CLD-associated liver fibrosis and potential biomarkers and/or targets for steatosis-associated liver fibrosis.

Identification of the genomic region under epigenetic regulation during non-alcoholic fatty liver disease progression

AIM

The progression of non-alcoholic fatty liver disease (NAFLD) is affected by epigenetics. We undertook co-methylation and differentially methylated region (DMR) analyses to identify the genomic region that is under epigenetic regulation during NAFLD progression.

METHODS

We collected liver biopsy specimens from 60 Japanese patients with NAFLD and classified these into mild (fibrosis stages 0-2) or advanced (fibrosis stages 3-4) NAFLD. We carried out a genome-wide DNA methylation analysis and identified the differentially methylated CpGs between mild and advanced NAFLD. Differentially methylated regions with multiple consecutive differentially methylated CpGs between mild and advanced NAFLD were extracted.

RESULTS

Co-methylation analysis showed that individual differentially methylated CpG sites were clustered into three modules. The CpG sites clustered in one module were hypomethylated in advanced NAFLD and their annotated genes were enriched for "immune system" function. The CpG sites in another module were hypermethylated and their annotated genes were enriched for "mitochondria" or "lipid particle", and "lipid metabolism" or "oxidoreductase activity". Hypomethylated DMRs included tumorigenesis-related genes (FGFR2, PTGFRN, and ZBTB38), the expressions of which are upregulated in advanced NAFLD. Tumor suppressor MGMT had two DMRs and was downregulated. Conversely, FBLIM1 and CYR61, encoding proteins that reduce cell proliferation, showed hypomethylated DMRs and were upregulated. Expression of the antioxidant gene NQO1 was upregulated, with a hypomethylated DMR. The DMR containing cancer-related MIR21 was hypomethylated in advanced NAFLD.

CONCLUSIONS

Co-methylation and DMR analyses suggest that the NAFLD liver undergoes mitochondrial dysfunction, decreased lipid metabolism, and impaired oxidoreductase activity, and acquires tumorigenic potential at the epigenetic level.