Predicting gene regulatory networks from cell atlases

Gene regulatory networks in disease and ageing

The precise control of gene expression is required for the maintenance of cellular homeostasis and proper cellular function, and the declining control of gene expression with age is considered a major contributor to age-associated changes in cellular physiology and disease. The coordination of gene expression can be represented through models of the molecular interactions that govern gene expression levels, so-called gene regulatory networks. Gene regulatory networks can represent interactions that occur through signal transduction, those that involve regulatory transcription factors, or statistical models of gene-gene relationships based on the premise that certain sets of genes tend to be coexpressed across a range of conditions and cell types. Advances in experimental and computational technologies have enabled the inference of these networks on an unprecedented scale and at unprecedented precision. Here, we delineate different types of gene regulatory networks and their cell-biological interpretation. We describe methods for inferring such networks from large-scale, multi-omics datasets and present applications that have aided our understanding of cellular ageing and disease mechanisms.

Comprehensive characterization of IFNγ signaling in acute myeloid leukemia reveals prognostic and therapeutic strategies

Interferon gamma (IFNγ) is a critical cytokine known for its diverse roles in immune regulation, inflammation, and tumor surveillance. However, while IFNγ levels were elevated in sera of most newly diagnosed acute myeloid leukemia (AML) patients, its complex interplay in AML remains insufficiently understood. We aim to characterize these complex interactions through comprehensive bulk and single-cell approaches in bone marrow of newly diagnosed AML patients. We identify monocytic AML as having a unique microenvironment characterized by IFNγ producing T and NK cells, high IFNγ signaling, and immunosuppressive features. IFNγ signaling score strongly correlates with venetoclax resistance in primary AML patient cells. Additionally, IFNγ treatment of primary AML patient cells increased venetoclax resistance. Lastly, a parsimonious 47-gene IFNγ score demonstrates robust prognostic value. In summary, our findings suggest that inhibiting IFNγ is a potential treatment strategy to overcoming venetoclax resistance and immune evasion in AML patients.

Benchmarking causal reasoning algorithms for gene expression-based compound mechanism of action analysis

BACKGROUND

Elucidating compound mechanism of action (MoA) is beneficial to drug discovery, but in practice often represents a significant challenge. Causal Reasoning approaches aim to address this situation by inferring dysregulated signalling proteins using transcriptomics data and biological networks; however, a comprehensive benchmarking of such approaches has not yet been reported. Here we benchmarked four causal reasoning algorithms (SigNet, CausalR, CausalR ScanR and CARNIVAL) with four networks (the smaller Omnipath network vs. 3 larger MetaBase™ networks), using LINCS L1000 and CMap microarray data, and assessed to what extent each factor dictated the successful recovery of direct targets and compound-associated signalling pathways in a benchmark dataset comprising 269 compounds. We additionally examined impact on performance in terms of the functions and roles of protein targets and their connectivity bias in the prior knowledge networks.

RESULTS

According to statistical analysis (negative binomial model), the combination of algorithm and network most significantly dictated the performance of causal reasoning algorithms, with the SigNet recovering the greatest number of direct targets. With respect to the recovery of signalling pathways, CARNIVAL with the Omnipath network was able to recover the most informative pathways containing compound targets, based on the Reactome pathway hierarchy. Additionally, CARNIVAL, SigNet and CausalR ScanR all outperformed baseline gene expression pathway enrichment results. We found no significant difference in performance between L1000 data or microarray data, even when limited to just 978 'landmark' genes. Notably, all causal reasoning algorithms also outperformed pathway recovery based on input DEGs, despite these often being used for pathway enrichment. Causal reasoning methods performance was somewhat correlated with connectivity and biological role of the targets.

CONCLUSIONS

Overall, we conclude that causal reasoning performs well at recovering signalling proteins related to compound MoA upstream from gene expression changes by leveraging prior knowledge networks, and that the choice of network and algorithm has a profound impact on the performance of causal reasoning algorithms. Based on the analyses presented here this is true for both microarray-based gene expression data as well as those based on the L1000 platform.

Susceptibility identification for seasonal influenza A/H3N2 based on baseline blood transcriptome

Introduction

Influenza susceptibility difference is a widely existing trait that has great practical significance for the accurate prevention and control of influenza.

Methods

Here, we focused on the human susceptibility to the seasonal influenza A/H3N2 of healthy adults at baseline level. Whole blood expression data for influenza A/H3N2 susceptibility from GEO were collected firstly (30 symptomatic and 19 asymptomatic). Then to explore the differences at baseline, a suite of systems biology approaches - the differential expression analysis, co-expression network analysis, and immune cell frequencies analysis were utilized.

Results

We found the baseline condition, especially immune condition between symptomatic and asymptomatic, was different. Co-expression module that is positively related to asymptomatic is also related to immune cell type of naïve B cell. Function enrichment analysis showed significantly correlation with "B cell receptor signaling pathway", "immune response-activating cell surface receptor signaling pathway" and so on. Also, modules that are positively related to symptomatic are also correlated to immune cell type of neutrophils, with function enrichment analysis showing significantly correlations with "response to bacterium", "inflammatory response", "cAMP-dependent protein kinase complex" and so on. Responses of symptomatic and asymptomatic hosts after virus exposure show differences on resisting the virus, with more effective frontline defense for asymptomatic hosts. A prediction model was also built based on only baseline transcription information to differentiate symptomatic and asymptomatic population with accuracy of 0.79.

Discussion

The results not only improve our understanding of the immune system and influenza susceptibility, but also provide a new direction for precise and targeted prevention and therapy of influenza.

Dynamics of hepatocyte-cholangiocyte cell-fate decisions during liver development and regeneration

The immense regenerative potential of the liver is attributed to the ability of its two key cell types – hepatocytes and cholangiocytes – to trans-differentiate to one another either directly or through intermediate progenitor states. However, the dynamic features of decision-making between these cell-fates during liver development and regeneration remains elusive. Here, we identify a core gene regulatory network comprising c/EBPα, TGFBR2, and SOX9 which is multistable in nature, enabling three distinct cell states – hepatocytes, cholangiocytes, and liver progenitor cells (hepatoblasts/oval cells) – and stochastic switching among them. Predicted expression signature for these three states are validated through multiple bulk and single-cell transcriptomic datasets collected across developmental stages and injury-induced liver repair. This network can also explain the experimentally observed spatial organization of phenotypes in liver parenchyma and predict strategies for efficient cellular reprogramming. Our analysis elucidates how the emergent dynamics of underlying regulatory networks drive diverse cell-fate decisions in liver development and regeneration.

•

Identified minimal regulatory network to model liver development and regeneration

•

Changes in phenotypic landscapes by in-silico perturbations of regulatory networks

•

Ability to explain physiological spatial patterning of liver cell types

•

Decoded strategies for efficient reprogramming among liver cell phenotypes

TENET: gene network reconstruction using transfer entropy reveals key regulatory factors from single cell transcriptomic data

Accurate prediction of gene regulatory rules is important towards understanding of cellular processes. Existing computational algorithms devised for bulk transcriptomics typically require a large number of time points to infer gene regulatory networks (GRNs), are applicable for a small number of genes and fail to detect potential causal relationships effectively. Here, we propose a novel approach 'TENET' to reconstruct GRNs from single cell RNA sequencing (scRNAseq) datasets. Employing transfer entropy (TE) to measure the amount of causal relationships between genes, TENET predicts large-scale gene regulatory cascades/relationships from scRNAseq data. TENET showed better performance than other GRN reconstructors, in identifying key regulators from public datasets. Specifically from scRNAseq, TENET identified key transcriptional factors in embryonic stem cells (ESCs) and during direct cardiomyocytes reprogramming, where other predictors failed. We further demonstrate that known target genes have significantly higher TE values, and TENET predicted higher TE genes were more influenced by the perturbation of their regulator. Using TENET, we identified and validated that Nme2 is a culture condition specific stem cell factor. These results indicate that TENET is uniquely capable of identifying key regulators from scRNAseq data.