Core liver homeostatic co-expression networks are preserved but respond to perturbations in an organism- and disease-specific manner

Critical reasoning on the co-expression module QTL in the dorsolateral prefrontal cortex

Expression quantitative trait locus (eQTL) analysis is a popular method of gaining insight into the function of regulatory variation. While cis-eQTL resources have been instrumental in linking genome-wide association study variants to gene function, complex trait heritability may be additionally mediated by other forms of gene regulation. Toward this end, novel eQTL methods leverage gene co-expression (module-QTL) to investigate joint regulation of gene modules by single genetic variants. Here we broadly define a "module-QTL" as the association of a genetic variant with a summary measure of gene co-expression. This approach aims to reduce the multiple testing burden of a trans-eQTL search through the consolidation of gene-based testing and provide biological context to eQTLs shared between genes. In this article we provide an in-depth examination of the co-expression module eQTL (module-QTL) through literature review, theoretical investigation, and real-data application of the module-QTL to three large prefrontal cortex genotype-RNA sequencing datasets. We find module-QTLs in our study that are disease associated and reproducible are not additionally informative beyond cis- or trans-eQTLs for module genes. Through comparison to prior studies, we highlight promises and limitations of the module-QTL across study designs and provide recommendations for further investigation of the module-QTL framework.

Exploring the role of transcriptional and post-transcriptional processes in mRNA co-expression

Co-expression of two or more genes at the single-cell level is usually associated with functional co-regulation. While mRNA co-expression-measured as the correlation in mRNA levels-can be influenced by both transcriptional and post-transcriptional events, transcriptional regulation is typically considered dominant. We review and connect the literature describing transcriptional and post-transcriptional regulation of co-expression. To enhance our understanding, we integrate four datasets spanning single-cell gene expression data, single-cell promoter activity data and individual transcript half-lives. Confirming expectations, we find that positive co-expression necessitates promoter coordination and similar mRNA half-lives. Surprisingly, negative co-expression is favored by differences in mRNA half-lives, contrary to initial predictions from stochastic simulations. Notably, this association manifests specifically within clusters of genes. We further observe a striking compensation between promoter coordination and mRNA half-lives, which additional stochastic simulations suggest might give rise to the observed co-expression patterns. These findings raise intriguing questions about the functional advantages conferred by this compensation between distal kinetic steps.

Polo-like kinase-1 mediates hepatitis C virus-induced cell migration, a drug target for liver cancer

Polo-like kinase 1 (PLK1) is a regulator of cell mitosis and cytoskeletal dynamics. PLK1 overexpression in liver cancer is associated with tumour progression, metastasis, and vascular invasion. Hepatitis C virus (HCV) NS5A protein stimulates PLK1-mediated phosphorylation of host proteins, so we hypothesised that HCV-PLK1 interactions might be a mechanism for HCV-induced liver cancer. We used a HCV cell-culture model (Jc1) to investigate the effects of virus infection on the cytoskeleton. In HCV-infected cells, a novel posttranslational modification in β-actin was observed with phosphorylation at Ser239. Using in silico and in vitro approaches, we identified PLK1 as the mediating kinase. In functional experiments with a phosphomimetic mutant form of β-actin, Ser239 phosphorylation influences β-actin polymerization and distribution, resulting in increased cell motility. The changes were prevented by treating cells with the PLK1 inhibitor volasertib. In HCV-infected hepatocytes, increased cell motility contributes to cancer cell migration, invasion, and metastasis. PLK1 is an important mediator of these effects and early treatment with PLK1 inhibitors may prevent or reduce HCC progression, particularly in people with HCV-induced HCC.

Identifying multiscale translational safety biomarkers using a network-based systems approach

Animal testing is the current standard for drug and chemicals safety assessment, but hazards translation to human is uncertain. Human in vitro models can address the species translation but might not replicate in vivo complexity. Herein, we propose a network-based method addressing these translational multiscale problems that derives in vivo liver injury biomarkers applicable to in vitro human early safety screening. We applied weighted correlation network analysis (WGCNA) to a large rat liver transcriptomic dataset to obtain co-regulated gene clusters (modules). We identified modules statistically associated with liver pathologies, including a module enriched for ATF4-regulated genes as associated with the occurrence of hepatocellular single-cell necrosis, and as preserved in human liver in vitro models. Within the module, we identified TRIB3 and MTHFD2 as a novel candidate stress biomarkers, and developed and used BAC-eGFPHepG2 reporters in a compound screening, identifying compounds showing ATF4-dependent stress response and potential early safety signals.

Integrative network-based analysis on multiple Gene Expression Omnibus datasets identifies novel immune molecular markers implicated in non-alcoholic steatohepatitis

Introduction

Non-alcoholic steatohepatitis (NASH), an advanced subtype of non-alcoholic fatty liver disease (NAFLD), has becoming the most important aetiology for end-stage liver disease, such as cirrhosis and hepatocellular carcinoma. This study were designed to explore novel genes associated with NASH.

Methods

Here, five independent Gene Expression Omnibus (GEO) datasets were combined into a single cohort and analyzed using network biology approaches.

Results

11 modules identified by weighted gene co-expression network analysis (WGCNA) showed significant association with the status of NASH. Further characterization of four gene modules of interest demonstrated that molecular pathology of NASH involves the upregulation of hub genes related to immune response, cholesterol and lipid metabolic process, extracellular matrix organization, and the downregulation of hub genes related to cellular amino acid catabolic, respectively. After DEGs enrichment analysis and module preservation analysis, the Turquoise module associated with immune response displayed a remarkably correlation with NASH status. Hub genes with high degree of connectivity in the module, including CD53, LCP1, LAPTM5, NCKAP1L, C3AR1, PLEK, FCER1G, HLA-DRA and SRGN were further verified in clinical samples and mouse model of NASH. Moreover, single-cell RNA-seq analysis showed that those key genes were expressed by distinct immune cells such as microphages, natural killer, dendritic, T and B cells. Finally, the potential transcription factors of Turquoise module were characterized, including NFKB1, STAT3, RFX5, ILF3, ELF1, SPI1, ETS1 and CEBPA, the expression of which increased with NASH progression.

Discussion

In conclusion, our integrative analysis will contribute to the understanding of NASH and may enable the development of potential biomarkers for NASH therapy.

Loss of FOCAD, operating via the SKI messenger RNA surveillance pathway, causes a pediatric syndrome with liver cirrhosis

Cirrhosis is usually a late-onset and life-threatening disease characterized by fibrotic scarring and inflammation that disrupts liver architecture and function. While it is typically the result of alcoholism or hepatitis viral infection in adults, its etiology in infants is much less understood. In this study, we report 14 children from ten unrelated families presenting with a syndromic form of pediatric liver cirrhosis. By genome/exome sequencing, we found recessive variants in FOCAD segregating with the disease. Zebrafish lacking focad phenocopied the human disease, revealing a signature of altered messenger RNA (mRNA) degradation processes in the liver. Using patient's primary cells and CRISPR-Cas9-mediated inactivation in human hepatic cell lines, we found that FOCAD deficiency compromises the SKI mRNA surveillance pathway by reducing the levels of the RNA helicase SKIC2 and its cofactor SKIC3. FOCAD knockout hepatocytes exhibited lowered albumin expression and signs of persistent injury accompanied by CCL2 overproduction. Our results reveal the importance of FOCAD in maintaining liver homeostasis and disclose a possible therapeutic intervention point via inhibition of the CCL2/CCR2 signaling axis.

A network-based approach reveals the dysregulated transcriptional regulation in non-alcoholic fatty liver disease

Non-alcoholic fatty liver disease (NAFLD) is a leading cause of chronic liver disease worldwide. We performed network analysis to investigate the dysregulated biological processes in the disease progression and revealed the molecular mechanism underlying NAFLD. Based on network analysis, we identified a highly conserved disease-associated gene module across three different NAFLD cohorts and highlighted the predominant role of key transcriptional regulators associated with lipid and cholesterol metabolism. In addition, we revealed the detailed metabolic differences between heterogeneous NAFLD patients through integrative systems analysis of transcriptomic data and liver-specific genome-scale metabolic model. Furthermore, we identified transcription factors (TFs), including SREBF2, HNF4A, SREBF1, YY1, and KLF13, showing regulation of hepatic expression of genes in the NAFLD-associated modules and validated the TFs using data generated from a mouse NAFLD model. In conclusion, our integrative analysis facilitates the understanding of the regulatory mechanism of these perturbed TFs and their associated biological processes.

•

Disease-associated gene modules are conserved across multiple NAFLD cohorts

•

The central genes in disease-associated modules are key enzymes in cholesterol synthesis

•

YY1 and KLF13 are potential key transcriptional regulators of NAFLD development

Found in translation-core network preservation across liver diseases and species

Esmaili et al.¹ conducted co-expression network analysis to uncover core homeostatic modules critical to a broad spectrum of liver diseases in mouse and human. Perturbation state of core modules may underlie disease stages across species and serve as therapeutic targets.