An approach for determining and measuring network hierarchy applied to comparing the phosphorylome and the regulome

Representing core gene expression activity relationships using the latent structure implicit in Bayesian networks

MOTIVATION

Many types of networks, such as co-expression or ChIP-seq-based gene-regulatory networks, provide useful information for biomedical studies. However, they are often too full of connections and difficult to interpret, forming "indecipherable hairballs."

RESULTS

To address this issue, we propose that a Bayesian network can summarize the core relationships between gene expression activities. This network, which we call the LatentDAG, is substantially simpler than conventional co-expression network and ChIP-seq networks (by two orders of magnitude). It provides clearer clusters, without extraneous cross-cluster connections, and clear separators between modules. Moreover, one can find a number of clear examples showing how it bridges the connection between steps in the transcriptional regulatory network and other networks (e.g. RNA-binding protein). In conjunction with a graph neural network, the LatentDAG works better than other biological networks in a variety of tasks, including prediction of gene conservation and clustering genes.

AVAILABILITY AND IMPLEMENTATION

Code is available at https://github.com/gersteinlab/LatentDAG.

Single-cell network biology characterizes cell type gene regulation for drug repurposing and phenotype prediction in Alzheimer's disease

Dysregulation of gene expression in Alzheimer's disease (AD) remains elusive, especially at the cell type level. Gene regulatory network, a key molecular mechanism linking transcription factors (TFs) and regulatory elements to govern gene expression, can change across cell types in the human brain and thus serve as a model for studying gene dysregulation in AD. However, AD-induced regulatory changes across brain cell types remains uncharted. To address this, we integrated single-cell multi-omics datasets to predict the gene regulatory networks of four major cell types, excitatory and inhibitory neurons, microglia and oligodendrocytes, in control and AD brains. Importantly, we analyzed and compared the structural and topological features of networks across cell types and examined changes in AD. Our analysis shows that hub TFs are largely common across cell types and AD-related changes are relatively more prominent in some cell types (e.g., microglia). The regulatory logics of enriched network motifs (e.g., feed-forward loops) further uncover cell type-specific TF-TF cooperativities in gene regulation. The cell type networks are also highly modular and several network modules with cell-type-specific expression changes in AD pathology are enriched with AD-risk genes. The further disease-module-drug association analysis suggests cell-type candidate drugs and their potential target genes. Finally, our network-based machine learning analysis systematically prioritized cell type risk genes likely involved in AD. Our strategy is validated using an independent dataset which showed that top ranked genes can predict clinical phenotypes (e.g., cognitive impairment) of AD with reasonable accuracy. Overall, this single-cell network biology analysis provides a comprehensive map linking genes, regulatory networks, cell types and drug targets and reveals cell-type gene dysregulation in AD.

Graph hierarchy: a novel framework to analyse hierarchical structures in complex networks

Trophic coherence, a measure of a graph's hierarchical organisation, has been shown to be linked to a graph's structural and dynamical aspects such as cyclicity, stability and normality. Trophic levels of vertices can reveal their functional properties, partition and rank the vertices accordingly. Trophic levels and hence trophic coherence can only be defined on graphs with basal vertices, i.e. vertices with zero in-degree. Consequently, trophic analysis of graphs had been restricted until now. In this paper we introduce a hierarchical framework which can be defined on any simple graph. Within this general framework, we develop several metrics: hierarchical levels, a generalisation of the notion of trophic levels, influence centrality, a measure of a vertex's ability to influence dynamics, and democracy coefficient, a measure of overall feedback in the system. We discuss how our generalisation relates to previous attempts and what new insights are illuminated on the topological and dynamical aspects of graphs. Finally, we show how the hierarchical structure of a network relates to the incidence rate in a SIS epidemic model and the economic insights we can gain through it.

Identification of human glucocorticoid response markers using integrated multi-omic analysis from a randomized crossover trial

Background:

Glucocorticoids are among the most commonly prescribed drugs, but there is no biomarker that can quantify their action. The aim of the study was to identify and validate circulating biomarkers of glucocorticoid action.

Methods:

In a randomized, crossover, single-blind, discovery study, 10 subjects with primary adrenal insufficiency (and no other endocrinopathies) were admitted at the in-patient clinic and studied during physiological glucocorticoid exposure and withdrawal. A randomization plan before the first intervention was used. Besides mild physical and/or mental fatigue and salt craving, no serious adverse events were observed. The transcriptome in peripheral blood mononuclear cells and adipose tissue, plasma miRNAomic, and serum metabolomics were compared between the interventions using integrated multi-omic analysis.

Results:

We identified a transcriptomic profile derived from two tissues and a multi-omic cluster, both predictive of glucocorticoid exposure. A microRNA (miR-122-5p) that was correlated with genes and metabolites regulated by glucocorticoid exposure was identified (p=0.009) and replicated in independent studies with varying glucocorticoid exposure (0.01 ≤ p≤0.05).

Conclusions:

We have generated results that construct the basis for successful discovery of biomarker(s) to measure effects of glucocorticoids, allowing strategies to individualize and optimize glucocorticoid therapy, and shedding light on disease etiology related to unphysiological glucocorticoid exposure, such as in cardiovascular disease and obesity.

Funding:

The Swedish Research Council (Grant 2015-02561 and 2019-01112); The Swedish federal government under the LUA/ALF agreement (Grant ALFGBG-719531); The Swedish Endocrinology Association; The Gothenburg Medical Society; Wellcome Trust; The Medical Research Council, UK; The Chief Scientist Office, UK; The Eva Madura’s Foundation; The Research Foundation of Copenhagen University Hospital; and The Danish Rheumatism Association.

Clinical trial number:

NCT02152553.

Several diseases, including asthma, arthritis, some skin conditions, and cancer, are treated with medications called glucocorticoids, which are synthetic versions of human hormones. These drugs are also used to treat people with a condition call adrenal insufficiency who do not produce enough of an important hormone called cortisol. Use of glucocorticoids is very common, the proportion of people in a given country taking them can range from 0.5% to 21% of the population depending on the duration of the treatment. But, like any medication, glucocorticoids have both benefits and risks: people who take glucocorticoids for a long time have an increased risk of diabetes, obesity, cardiovascular disease, and death.

Because of the risks associated with taking glucocorticoids, it is very important for physicians to tailor the dose to each patient’s needs. Doing this can be tricky, because the levels of glucocorticoids in a patient’s blood are not a good indicator of the medication’s activity in the body. A test that can accurately measure the glucocorticoid activity could help physicians personalize treatment and reduce harmful side effects.

As a first step towards developing such a test, Chantzichristos et al. identified a potential way to measure glucocorticoid activity in patient’s blood. In the experiments, blood samples were collected from ten patients with adrenal insufficiency both when they were on no medication, and when they were taking a glucocorticoid to replace their missing hormones. Next, the blood samples were analyzed to determine which genes were turned on and off in each patient with and without the medication. They also compared small molecules in the blood called metabolites and tiny pieces of genetic material called microRNAs that turn genes on and off.

The experiments revealed networks of genes, metabolites, and microRNAs that are associated with glucocorticoid activity, and one microRNA called miR-122-5p stood out as a potential way to measure glucocorticoid activity. To verify this microRNA’s usefulness, Chantzichristos et al. looked at levels of miR-122-5p in people participating in three other studies and confirmed that it was a good indicator of the glucocorticoid activity.

More research is needed to confirm Chantzichristos et al.’s findings and to develop a test that can be used by physicians to measure glucocorticoid activity. The microRNA identified, miR-122-5p, has been previously linked to diabetes, so studying it further may also help scientists understand how taking glucocorticoids may increase the risk of developing diabetes and related diseases.

Integrative inference of transcriptional networks in Arabidopsis yields novel ROS signalling regulators

Gene regulation is a dynamic process in which transcription factors (TFs) play an important role in controlling spatiotemporal gene expression. To enhance our global understanding of regulatory interactions in Arabidopsis thaliana, different regulatory input networks capturing complementary information about DNA motifs, open chromatin, TF-binding and expression-based regulatory interactions were combined using a supervised learning approach, resulting in an integrated gene regulatory network (iGRN) covering 1,491 TFs and 31,393 target genes (1.7 million interactions). This iGRN outperforms the different input networks to predict known regulatory interactions and has a similar performance to recover functional interactions compared to state-of-the-art experimental methods. The iGRN correctly inferred known functions for 681 TFs and predicted new gene functions for hundreds of unknown TFs. For regulators predicted to be involved in reactive oxygen species (ROS) stress regulation, we confirmed in total 75% of TFs with a function in ROS and/or physiological stress responses. This includes 13 ROS regulators, previously not connected to any ROS or stress function, that were experimentally validated in our ROS-specific phenotypic assays of loss- or gain-of-function lines. In conclusion, the presented iGRN offers a high-quality starting point to enhance our understanding of gene regulation in plants by integrating different experimental data types.

TGF-beta: a master immune regulator

Introduction: Transforming Growth Factor-Beta (TGF-β) is a master regulator of numerous cellular functions including cellular immunity. In cancer, TGF-β can function as a tumor promoter via several mechanisms including immunosuppression. Since the immune checkpoint pathways are co-opted in cancer to induce T cell tolerance, this review posits that TGF-β is a master checkpoint in cancer, whose negative regulatory influence overrides and controls that of other immune checkpoints.Areas Covered: This review examines therapeutic agents that target TGF-β and its signaling pathways for the treatment of cancer which may be classifiable as checkpoint inhibitors in the broadest sense. This concept is supported by the observations that 1) only a subset of patients benefit from current checkpoint inhibitor therapies, 2) the presence of TGF-β in the tumor microenvironment is associated with excluded or cold tumors, and resistance to checkpoint inhibitors, and 3) existing biomarkers such as PD-1, PD-L1, microsatellite instability and tumor mutational burden are inadequate to reliably and adequately identify immuno-responsive patients. By contrast, TGF-β overexpression is a widespread and profoundly negative molecular hallmark in multiple tumor types.Expert Opinion: TGF-β status may serve as a biomarker to predict responsiveness and as a therapeutic target to increase the activity of immunotherapies.

Plants are intelligent, here's how

HYPOTHESES

The drive to survive is a biological universal. Intelligent behaviour is usually recognized when individual organisms including plants, in the face of fiercely competitive or adverse, real-world circumstances, change their behaviour to improve their probability of survival.

SCOPE

This article explains the potential relationship of intelligence to adaptability and emphasizes the need to recognize individual variation in intelligence showing it to be goal directed and thus being purposeful. Intelligent behaviour in single cells and microbes is frequently reported. Individual variation might be underpinned by a novel learning mechanism, described here in detail. The requirements for real-world circumstances are outlined, and the relationship to organic selection is indicated together with niche construction as a good example of intentional behaviour that should improve survival. Adaptability is important in crop development but the term may be complex incorporating numerous behavioural traits some of which are indicated.

CONCLUSION

There is real biological benefit to regarding plants as intelligent both from the fundamental issue of understanding plant life but also from providing a direction for fundamental future research and in crop breeding.

The Architecture of the Human RNA-Binding Protein Regulatory Network

RNA-binding proteins (RBPs) are key players of post-transcriptional regulation of gene expression, relying on competitive and cooperative interactions to fine-tune their action. Several studies have described individual interactions of RBPs with RBP mRNAs. Here we present a systematic network investigation of fifty thousand interactions between RBPs and the UTRs of RBP mRNAs. We identified two structural features in this network. RBP clusters are groups of densely interconnected RBPs co-binding their targets, suggesting a tight control of cooperative and competitive behaviors. RBP chains are hierarchical structures connecting RBP clusters and driven by evolutionarily ancient RBPs. These features suggest that RBP chains may coordinate the different cell programs controlled by RBP clusters. Under this model, the regulatory signal flows through chains from one cluster to another, implementing elaborate regulatory plans. This work thus suggests RBP-RBP interactions as a backbone driving post-transcriptional regulation of gene expression to control RBPs action on their targets.

•

The RBP-RBP network is a robust and efficient hierarchical structure

•

The RBP-RBP network is formed by RBP clusters and chains

•

RBP-RBP interactions cluster to cooperate and compete on common target mRNAs

•

RBP chains are master regulatory units of the cell led by ancient RBPs

Hierarchical cooperation of transcription factors from integration analysis of DNA sequences, ChIP-Seq and ChIA-PET data

Background

Chromosomal architecture, which is constituted by chromatin loops, plays an important role in cellular functions. Gene expression and cell identity can be regulated by the chromatin loop, which is formed by proximal or distal enhancers and promoters in linear DNA (1D). Enhancers and promoters are fundamental non-coding elements enriched with transcription factors (TFs) to form chromatin loops. However, the specific cooperation of TFs involved in forming chromatin loops is not fully understood.

Results

Here, we proposed a method for investigating the cooperation of TFs in four cell lines by the integrative analysis of DNA sequences, ChIP-Seq and ChIA-PET data. Results demonstrate that the interaction of enhancers and promoters is a hierarchical and dynamic complex process with cooperative interactions of different TFs synergistically regulating gene expression and chromatin structure. The TF cooperation involved in maintaining and regulating the chromatin loop of cells can be regulated by epigenetic factors, such as other TFs and DNA methylation.

Conclusions

Such cooperation among TFs provides the potential features that can affect chromatin’s 3D architecture in cells. The regulation of chromatin 3D organization and gene expression is a complex process associated with the hierarchical and dynamic prosperities of TFs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5535-2) contains supplementary material, which is available to authorized users.

Novel putative drugs and key initiating genes for neurodegenerative disease determined using network-based genetic integrative analysis

Understanding the genetic causes of neurodegenerative disease (ND) can be useful for their prevention and treatment. Among the genetic variations responsible for ND, heritable germline variants have been discovered in genome-wide association studies (GWAS), and nonheritable somatic mutations have been discovered in sequencing projects. Distinguishing the important initiating genes in ND and comparing the importance of heritable and nonheritable genetic variants for treating ND are important challenges. In this study, we analysed GWAS results, somatic mutations and drug targets of ND from large databanks by performing directed network-based analysis considering a randomised network hypothesis testing procedure. A disease-associated biological network was created in the context of the functional interactome, and the nonrandom topological characteristics of directed-edge classes were interpreted. Hierarchical network analysis indicated that drug targets tend to lie upstream of somatic mutations and germline variants. Furthermore, using directed path length information and biological explanations, we provide information on the most important genes in these created node classes and their associated drugs. Finally, we identified nine germline variants overlapping with drug targets for ND, seven somatic mutations close to drug targets from the hierarchical network analysis and six crucial genes in controlling other genes from the network analysis. Based on these findings, some drugs have been proposed for treating ND via drug repurposing. Our results provide new insights into the therapeutic actionability of GWAS results and somatic mutations for ND. The interesting properties of each node class and the existing relationships between them can broaden our knowledge of ND.

A network of epigenomic and transcriptional cooperation encompassing an epigenomic master regulator in cancer

Coordinated experiments focused on transcriptional responses and chromatin states are well-equipped to capture different epigenomic and transcriptomic levels governing the circuitry of a regulatory network. We propose a workflow for the genome-wide identification of epigenomic and transcriptional cooperation to elucidate transcriptional networks in cancer. Gene promoter annotation in combination with network analysis and sequence-resolution of enriched transcriptional motifs in epigenomic data reveals transcription factor families that act synergistically with epigenomic master regulators. By investigating complementary omics levels, a close teamwork of the transcriptional and epigenomic machinery was discovered. The discovered network is tightly connected and surrounds the histone lysine demethylase KDM3A, basic helix-loop-helix factors MYC, HIF1A, and SREBF1, as well as differentiation factors AP1, MYOD1, SP1, MEIS1, ZEB1, and ELK1. In such a cooperative network, one component opens the chromatin, another one recognizes gene-specific DNA motifs, others scaffold between histones, cofactors, and the transcriptional complex. In cancer, due to the ability to team up with transcription factors, epigenetic factors concert mitogenic and metabolic gene networks, claiming the role of a cancer master regulators or epioncogenes. Significantly, specific histone modification patterns are commonly associated with open or closed chromatin states, and are linked to distinct biological outcomes by transcriptional activation or repression. Disruption of patterns of histone modifications is associated with the loss of proliferative control and cancer. There is tremendous therapeutic potential in understanding and targeting histone modification pathways. Thus, investigating cooperation of chromatin remodelers and the transcriptional machinery is not only important for elucidating fundamental mechanisms of chromatin regulation, but also necessary for the design of targeted therapeutics.

Similar bowtie structures and distinct largest strong components are identified in the transcriptional regulatory networks of Arabidopsis thaliana during photomorphogenesis and heat shock

Photomorphogenesis and heat shock are critical biological processes of plants. A recent research constructed the transcriptional regulatory networks (TRNs) of Arabidopsis thaliana during these processes using DNase-seq. In this study, by strong decomposition, we revealed that each of these TRNs can be represented as a similar bowtie structure with only one non-trivial and distinct strong component. We further identified distinct patterns of variation of a few light-related genes in these bowtie structures during photomorphogenesis. These results suggest that bowtie structure may be a common property of TRNs of plants, and distinct variation patterns of genes in bowtie structures of TRNs during biological processes may reflect distinct functions. Overall, our study provides an insight into the molecular mechanisms underlying photomorphogenesis and heat shock, and emphasizes the necessity to investigate the strong connectivity structures while studying TRNs.

Delineating functional principles of the bow tie structure of a kinase-phosphatase network in the budding yeast

Background

Kinases and phosphatases (KP) form complex self-regulating networks essential for cellular signal processing. In spite of having a wealth of data about interactions among KPs and their substrates, we have very limited models of the structures of the directed networks they form and consequently our ability to formulate hypotheses about how their structure determines the flow of information in these networks is restricted.

Results

We assembled and studied the largest bona fide kinase-phosphatase network (KP-Net) known to date for the yeast Saccharomyces cerevisiae. Application of the vertex sort (VS) algorithm on the KP-Net allowed us to elucidate its hierarchical structure in which nodes are sorted into top, core and bottom layers, forming a bow tie structure with a strongly connected core layer. Surprisingly, phosphatases tend to sort into the top layer, implying they are less regulated by phosphorylation than kinases. Superposition of the widest range of KP biological properties over the KP-Net hierarchy shows that core layer KPs: (i), receive the largest number of inputs; (ii), form bottlenecks implicated in multiple pathways and in decision-making; (iii), and are among the most regulated KPs both temporally and spatially. Moreover, top layer KPs are more abundant and less noisy than those in the bottom layer. Finally, we showed that the VS algorithm depends on node degrees without biasing the biological results of the sorted network. The VS algorithm is available as an R package (https://cran.r-project.org/web/packages/VertexSort/index.html).

Conclusions

The KP-Net model we propose possesses a bow tie hierarchical structure in which the top layer appears to ensure highest fidelity and the core layer appears to mediate signal integration and cell state-dependent signal interpretation. Our model of the yeast KP-Net provides both functional insight into its organization as we understand today and a framework for future investigation of information processing in yeast and eukaryotes in general.

Electronic supplementary material

The online version of this article (doi:10.1186/s12918-017-0418-0) contains supplementary material, which is available to authorized users.

Understanding Physiology in the Continuum: Integration of Information from Multiple -Omics Levels

In this paper, we discuss approaches for integrating biological information reflecting diverse physiologic levels. In particular, we explore statistical and model-based methods for integrating transcriptomic, proteomic and metabolomics data. Our case studies reflect responses to a systemic inflammatory stimulus and in response to an anti-inflammatory treatment. Our paper serves partly as a review of existing methods and partly as a means to demonstrate, using case studies related to human endotoxemia and response to methylprednisolone (MPL) treatment, how specific questions may require specific methods, thus emphasizing the non-uniqueness of the approaches. Finally, we explore novel ways for integrating -omics information with PKPD models, toward the development of more integrated pharmacology models.

Integrative analysis of cancer genes in a functional interactome

The post-genomic era has resulted in the accumulation of high-throughput cancer data from a vast array of genomic technologies including next-generation sequencing and microarray. As such, the large amounts of germline variant and somatic mutation data that have been generated from GWAS and sequencing projects, respectively, show great promise in providing a systems-level view of these genetic aberrations. In this study, we analyze publicly available GWAS, somatic mutation, and drug target data derived from large databanks using a network-based approach that incorporates directed edge information under a randomized network hypothesis testing procedure. We show that these three classes of disease-associated nodes exhibit non-random topological characteristics in the context of a functional interactome. Specifically, we show that drug targets tend to lie upstream of somatic mutations and disease susceptibility germline variants. In addition, we introduce a new approach to measuring hierarchy between drug targets, somatic mutants, and disease susceptibility genes by utilizing directionality and path length information. Overall, our results provide new insight into the intrinsic relationships between these node classes that broaden our understanding of cancer. In addition, our results align with current knowledge on the therapeutic actionability of GWAS and somatic mutant nodes, while demonstrating relationships between node classes from a global network perspective.

Reconstruction and topological characterization of the sigma factor regulatory network of Mycobacterium tuberculosis

Accessory sigma factors, which reprogram RNA polymerase to transcribe specific gene sets, activate bacterial adaptive responses to noxious environments. Here we reconstruct the complete sigma factor regulatory network of the human pathogen Mycobacterium tuberculosis by an integrated approach. The approach combines identification of direct regulatory interactions between M. tuberculosis sigma factors in an E. coli model system, validation of selected links in M. tuberculosis, and extensive literature review. The resulting network comprises 41 direct interactions among all 13 sigma factors. Analysis of network topology reveals (i) a three-tiered hierarchy initiating at master regulators, (ii) high connectivity and (iii) distinct communities containing multiple sigma factors. These topological features are likely associated with multi-layer signal processing and specialized stress responses involving multiple sigma factors. Moreover, the identification of overrepresented network motifs, such as autoregulation and coregulation of sigma and anti-sigma factor pairs, provides structural information that is relevant for studies of network dynamics.

Sigma factors are regulatory proteins that reprogram the bacterial RNA polymerase in response to stress conditions to transcribe certain genes, including those for other sigma factors. Here, Chauhan et al. describe the complete sigma factor regulatory network of the pathogen Mycobacterium tuberculosis.

Multi-hierarchical profiling: an emerging and quantitative approach to characterizing diverse biological networks

Multi-hierarchical profiling may offer valuable insights into the structural stability and functional direction of biological networks in cellular development, pathological process and disease variation. Owing to the emergence of several new techniques, such as bioinformatics for omics data, structural biology and structural bioinformatics, the pace of network hierarchical research has accelerated a lot in recent years. Here, we discuss and compare the techniques available for quantifying multilevel hierarchies, with a focus on their features, capabilities and drawbacks when used for different applications. Then, we classify these methods into three types: topological spatial-scales, multilevel hierarchical control and feature ordering. We observe that challenges and limitations do exist in functional hierarchical identification. And, we also provide useful suggestions on how to analyze the dynamic data of complex network studies.