The architecture of the gene regulatory networks of different tissues

Hierarchical network structure as the source of hierarchical dynamics (power-law frequency spectra) in living and non-living systems: How state-trait continua (body plans, personalities) emerge from first principles in biophysics

Living systems are hierarchical control systems that display a small world network structure. In such structures, many smaller clusters are nested within fewer larger ones, producing a fractal-like structure with a 'power-law' cluster size distribution (a mereology). Just like their structure, the dynamics of living systems shows fractal-like qualities: the timeseries of inner message passing and overt behavior contain high frequencies or 'states' (treble) that are nested within lower frequencies or 'traits' (bass), producing a power-law frequency spectrum that is known as a 'state-trait continuum' in the behavioral sciences. Here, we argue that the power-law dynamics of living systems results from their power-law network structure: organisms 'vertically encode' the deep spatiotemporal structure of their (anticipated) environments, to the effect that many small clusters near the base of the hierarchy produce high frequency signal changes and fewer larger clusters at its top produce ultra-low frequencies. Such ultra-low frequencies exert a tonic regulatory pressure that produces morphological as well as behavioral traits (i.e., body plans and personalities). Nested-modular structure causes higher frequencies to be embedded within lower frequencies, producing a power-law state-trait continuum. At the heart of such dynamics lies the need for efficient energy dissipation through networks of coupled oscillators, which also governs the dynamics of non-living systems (e.q., earthquakes, stock market fluctuations). Since hierarchical structure produces hierarchical dynamics, the development and collapse of hierarchical structure (e.g., during maturation and disease) should leave specific traces in system dynamics (shifts in lower frequencies, i.e. morphological and behavioral traits) that may serve as early warning signs to system failure. The applications of this idea range from (bio)physics and phylogenesis to ontogenesis and clinical medicine.

Identification of biomarkers complementary to homologous recombination deficiency for improving the clinical outcome of ovarian serous cystadenocarcinoma

Ovarian cancer patients with homologous recombination deficiency (HRD) tumors would benefit from PARP inhibitor (PARPi) therapy. However, patients with HRD tumors account for less than 50% of the whole cohort, so new biomarkers still need to be developed. Based on the data from the SNP array and somatic mutation profiles in the ovarian cancer genome, we found that high frequency of actionable mutations existed in patients with non-HRD tumors. Through transcriptome analysis, we identified that a downstream target of the cGAS-STING pathway, CXCL11, was upregulated in HRD tumors and could be used as a predictor of survival outcome. Further comprehensive analysis of the tumor immune microenvironment (TIME) revealed that CXCL11 expression signature was closely correlated with cytotoxic cells, neoantigen load and immune checkpoint blockade (ICB). Clinical trial data confirmed that the expression of CXCL11 could be used as a biomarker for anti-PD-1/PD-L1 therapy. Finally, in vivo and in vitro experiments showed that cancer cells with PARPi treatment increased the expression of CXCL11. Collectively, our study not only provides biomarkers of ovarian cancer complementary to the HRD score but also introduces a potential new perspective for identifying prognostic biomarkers of immunotherapy.

How higher goals are constructed and collapse under stress: A hierarchical Bayesian control systems perspective

In this paper, we show that organisms can be modeled as hierarchical Bayesian control systems with small world and information bottleneck (bow-tie) network structure. Such systems combine hierarchical perception with hierarchical goal setting and hierarchical action control. We argue that hierarchical Bayesian control systems produce deep hierarchies of goal states, from which it follows that organisms must have some form of 'highest goals'. For all organisms, these involve internal (self) models, external (social) models and overarching (normative) models. We show that goal hierarchies tend to decompose in a top-down manner under severe and prolonged levels of stress. This produces behavior that favors short-term and self-referential goals over long term, social and/or normative goals. The collapse of goal hierarchies is universally accompanied by an increase in entropy (disorder) in control systems that can serve as an early warning sign for tipping points (disease or death of the organism). In humans, learning goal hierarchies corresponds to personality development (maturation). The failure of goal hierarchies to mature properly corresponds to personality deficits. A top-down collapse of such hierarchies under stress is identified as a common factor in all forms of episodic mental disorders (psychopathology). The paper concludes by discussing ways of testing these hypotheses empirically.

MyomirDB: A unified database and server platform for muscle atrophy myomiRs, coregulatory networks and regulons

Muscular atrophy or muscle loss is a multifactorial clinical condition during many critical illnesses like cancer, cardiovascular diseases, diabetes, pulmonary diseases etc. leading to fatigue and weakness and contributes towards a decreased quality of life. The proportion of older adults (>65 y) in the overall population is also growing and aging is another important factor causing muscle loss. Some muscle miRNAs (myomiRs) and their target genes have even been proposed as potential diagnostic, therapeutic and predictive markers for muscular atrophy. MyomirDB (http://www.myomirdb.in/) is a unique resource that provides a comprehensive, curated, user- friendly and detailed compilation of various miRNA bio-molecular interactions; miRNA-Transcription Factor-Target Gene co-regulatory networks and ~8000 tripartite regulons associated with 247 myomiRs which have been experimentally validated to be associated with various muscular atrophy conditions. For each database entry, MyomirDB compiles source organism, muscle atrophic condition, experiment duration, its level of expression, fold change, tissue of expression, experimental validation, disease and drug association, tissue-specific expression level, Gene Ontology and KEGG pathway associations. The web resource is a unique server platform which uses in-house scripts to construct miRNA-Transcription Factor-Target Gene co-regulatory networks and extract tri-partite regulons also called Feed Forward Loops. These unique features helps to offer mechanistic insights in disease pathology. Hence, MyomirDB is a unique platform for researchers working in this area to explore, fetch, compare and analyse atrophy associated miRNAs, their co-regulatory networks and FFL regulons.

Function, dynamics and evolution of network motif modules in integrated gene regulatory networks of worm and plant

Gene regulatory networks (GRNs) consist of different molecular interactions that closely work together to establish proper gene expression in time and space. Especially in higher eukaryotes, many questions remain on how these interactions collectively coordinate gene regulation. We study high quality GRNs consisting of undirected protein–protein, genetic and homologous interactions, and directed protein–DNA, regulatory and miRNA–mRNA interactions in the worm Caenorhabditis elegans and the plant Arabidopsis thaliana. Our data-integration framework integrates interactions in composite network motifs, clusters these in biologically relevant, higher-order topological network motif modules, overlays these with gene expression profiles and discovers novel connections between modules and regulators. Similar modules exist in the integrated GRNs of worm and plant. We show how experimental or computational methodologies underlying a certain data type impact network topology. Through phylogenetic decomposition, we found that proteins of worm and plant tend to functionally interact with proteins of a similar age, while at the regulatory level TFs favor same age, but also older target genes. Despite some influence of the duplication mode difference, we also observe at the motif and module level for both species a preference for age homogeneity for undirected and age heterogeneity for directed interactions. This leads to a model where novel genes are added together to the GRNs in a specific biological functional context, regulated by one or more TFs that also target older genes in the GRNs. Overall, we detected topological, functional and evolutionary properties of GRNs that are potentially universal in all species.

MicroRNA Profiling of Salivary Duct Carcinoma Versus Her2/Neu Overexpressing Breast Carcinoma Identify miR-10a as a Putative Breast Related Oncogene

Salivary duct carcinomas (SDC) and Her2/Neu3-overexpressing invasive breast carcinomas (HNPIBC/IBC) are histologically indistinguishable. We investigated whether common histopathologic and immunophenotypic features of SDC and IBC are mirrored by a similar microRNA (miRNA) profile. MiRNA profiling of 5 SDCs, 6 IBCs Her2/Neu3+, and 5 high-grade ductal breast carcinoma in situ (DCIS) was performed by NanoString platform. Selected miRNAs and HOXA1 gene were validated by RT-PCR. We observed similar miRNA expression profiles between IBC and SDC with the exception of 2 miRNAs, miR-10a and miR-142-3p, which were higher in IBC tumors. DCIS tumors displayed increased expression of miR-10a, miR-99a, miR-331-3p and miR-335, and decreased expression of miR-15a, miR-16 and miR-19b compared to SDC. The normal salivary gland and breast tissues also showed similar expression profiles. Interestingly, miR-10a was selectively increased in both IBC and normal breast tissue compared to SDC and normal salivary gland tissue. Moreover, our NanoString and RT-PCR data confirmed that miR-10a was upregulated in IBC and DCIS compared to SDC. Finally, we show downregulation of HOXA1, a miR-10 target, in IBC tumors compared to normal breast tissue. Taken together, our data demonstrates that, based on miRNA profiling, SDC is closely related to HNPIBC. Our results also suggest that miR-10a is differentially expressed in IBC compared to SDC and may have potential utility as a diagnostic biomarker in synchronous or metachronous malignant epithelial malignancies involving both organs. In addition, miR-10a could be playing an important role as a mammary-specific oncogene, involved in breast cancer initiation (DCIS) and progression (IBC), through mechanisms that include modulation of HOXA1 gene expression.

The core regulatory network in human cells

In order to discover the common characteristics of various cell types in the human body, many researches have been conducted to find the set of genes commonly expressed in various cell types and tissues. However, the functional characteristics of a cell is determined by the complex regulatory relationships among the genes rather than by expressed genes themselves. Therefore, it is more important to identify and analyze a core regulatory network where all regulatory relationship between genes are active across all cell types to uncover the common features of various cell types. Here, based on hundreds of tissue-specific gene regulatory networks constructed by recent genome-wide experimental data, we constructed the core regulatory network. Interestingly, we found that the core regulatory network is organized by simple cascade and has few complex regulations such as feedback or feed-forward loops. Moreover, we discovered that the regulatory links from genes in the core regulatory network to genes in the peripheral regulatory network are much more abundant than the reverse direction links. These results suggest that the core regulatory network locates at the top of regulatory network and plays a role as a 'hub' in terms of information flow, and the information that is common to all cells can be modified to achieve the tissue-specific characteristics through various types of feedback and feed-forward loops in the peripheral regulatory networks. We also found that the genes in the core regulatory network are evolutionary conserved, essential and non-disease, non-druggable genes compared to the peripheral genes. Overall, our study provides an insight into how all human cells share a common function and generate tissue-specific functional traits by transmitting and processing information through regulatory network.

Chromatin states modify network motifs contributing to cell-specific functions

Epigenetic modification can affect many important biological processes, such as cell proliferation and apoptosis. It can alter chromatin conformation and contribute to gene regulation. To investigate how chromatin states associated with network motifs, we assembled chromatin state-modified regulatory networks by combining 269 ChIP-seq data and chromatin states in four cell types. We found that many chromatin states were significantly associated with network motifs, especially for feedforward loops (FFLs). These distinct chromatin state compositions contribute to different expression levels and translational control of targets in FFLs. Strikingly, the chromatin state-modified FFLs were highly cell-specific and, to a large extent, determined cell-selective functions, such as the embryonic stem cell-specific bivalent modification-related FFL with an important role in poising developmentally important genes for expression. Besides, comparisons of chromatin state-modified FFLs between cancerous/stem and primary cell lines revealed specific type of chromatin state alterations that may act together with motif structural changes cooperatively contribute to cell-to-cell functional differences. Combination of these alterations could be helpful in prioritizing candidate genes. Together, this work highlights that a dynamic epigenetic dimension can help network motifs to control cell-specific functions.

Regulatory networks in retinal ischemia-reperfusion injury

BACKGROUND

Retinal function is ordered by interactions between transcriptional and posttranscriptional regulators at the molecular level. These regulators include transcription factors (TFs) and posttranscriptional factors such as microRNAs (miRs). Some studies propose that miRs predominantly target the TFs rather than other types of protein coding genes and such studies suggest a possible interconnection of these two regulators in co-regulatory networks.

RESULTS

Our lab has generated mRNA and miRNA microarray expression data to investigate time-dependent changes in gene expression, following induction of ischemia-reperfusion (IR) injury in the rat retina. Data from different reperfusion time points following retinal IR-injury were analyzed. Paired expression data for miRNA-target gene (TG), TF-TG, miRNA-TF were used to identify regulatory loop motifs whose expressions were altered by the IR injury paradigm. These loops were subsequently integrated into larger regulatory networks and biological functions were assayed. Systematic analyses of the networks have provided new insights into retinal gene regulation in the early and late periods of IR. We found both overlapping and unique patterns of molecular expression at the two time points. These patterns can be defined by their characteristic molecular motifs as well as their associated biological processes. We highlighted the regulatory elements of miRs and TFs associated with biological processes in the early and late phases of ischemia-reperfusion injury.

CONCLUSIONS

The etiology of retinal ischemia-reperfusion injury is orchestrated by complex and still not well understood gene networks. This work represents the first large network analysis to integrate miRNA and mRNA expression profiles in context of retinal ischemia. It is likely that an appreciation of such regulatory networks will have prognostic potential. In addition, the computational framework described in this study can be used to construct miRNA-TF interactive systems networks for various diseases/disorders of the retina and other tissues.

Large differences in global transcriptional regulatory programs of normal and tumor colon cells

BACKGROUND

Dysregulation of transcriptional programs leads to cell malfunctioning and can have an impact in cancer development. Our study aims to characterize global differences between transcriptional regulatory programs of normal and tumor cells of the colon.

METHODS

Affymetrix Human Genome U219 expression arrays were used to assess gene expression in 100 samples of colon tumor and their paired adjacent normal mucosa. Transcriptional networks were reconstructed using ARACNe algorithm using 1,000 bootstrap replicates consolidated into a consensus network. Networks were compared regarding topology parameters and identified well-connected clusters. Functional enrichment was performed with SIGORA method. ENCODE ChIP-Seq data curated in the hmChIP database was used for in silico validation of the most prominent transcription factors.

RESULTS

The normal network contained 1,177 transcription factors, 5,466 target genes and 61,226 transcriptional interactions. A large loss of transcriptional interactions in the tumor network was observed (11,585; 81% reduction), which also contained fewer transcription factors (621; 47% reduction) and target genes (2,190; 60% reduction) than the normal network. Gene silencing was not a main determinant of this loss of regulatory activity, since the average gene expression was essentially conserved. Also, 91 transcription factors increased their connectivity in the tumor network. These genes revealed a tumor-specific emergent transcriptional regulatory program with significant functional enrichment related to colorectal cancer pathway. In addition, the analysis of clusters again identified subnetworks in the tumors enriched for cancer related pathways (immune response, Wnt signaling, DNA replication, cell adherence, apoptosis, DNA repair, among others). Also multiple metabolism pathways show differential clustering between the tumor and normal network.

CONCLUSIONS

These findings will allow a better understanding of the transcriptional regulatory programs altered in colon cancer and could be an invaluable methodology to identify potential hubs with a relevant role in the field of cancer diagnosis, prognosis and therapy.

FunMod: a Cytoscape plugin for identifying functional modules in undirected protein-protein networks

The characterization of the interacting behaviors of complex biological systems is a primary objective in protein-protein network analysis and computational biology. In this paper we present FunMod, an innovative Cytoscape version 2.8 plugin that is able to mine undirected protein-protein networks and to infer sub-networks of interacting proteins intimately correlated with relevant biological pathways. This plugin may enable the discovery of new pathways involved in diseases. In order to describe the role of each protein within the relevant biological pathways, FunMod computes and scores three topological features of the identified sub-networks. By integrating the results from biological pathway clustering and topological network analysis, FunMod proved to be useful for the data interpretation and the generation of new hypotheses in two case studies.

Integrative identification of deregulated miRNA/TF-mediated gene regulatory loops and networks in prostate cancer

MicroRNAs (miRNAs) have attracted a great deal of attention in biology and medicine. It has been hypothesized that miRNAs interact with transcription factors (TFs) in a coordinated fashion to play key roles in regulating signaling and transcriptional pathways and in achieving robust gene regulation. Here, we propose a novel integrative computational method to infer certain types of deregulated miRNA-mediated regulatory circuits at the transcriptional, post-transcriptional and signaling levels. To reliably predict miRNA-target interactions from mRNA/miRNA expression data, our method collectively utilizes sequence-based miRNA-target predictions obtained from several algorithms, known information about mRNA and miRNA targets of TFs available in existing databases, certain molecular structures identified to be statistically over-represented in gene regulatory networks, available molecular subtyping information, and state-of-the-art statistical techniques to appropriately constrain the underlying analysis. In this way, the method exploits almost every aspect of extractable information in the expression data. We apply our procedure on mRNA/miRNA expression data from prostate tumor and normal samples and detect numerous known and novel miRNA-mediated deregulated loops and networks in prostate cancer. We also demonstrate instances of the results in a number of distinct biological settings, which are known to play crucial roles in prostate and other types of cancer. Our findings show that the proposed computational method can be used to effectively achieve notable insights into the poorly understood molecular mechanisms of miRNA-mediated interactions and dissect their functional roles in cancer in an effort to pave the way for miRNA-based therapeutics in clinical settings.

Human-specific microRNA regulation of FOXO1: implications for microRNA recognition element evolution

MicroRNAs (miRNAs) have been established as important negative post-transcriptional regulators for gene expression. Within the past decade, miRNAs targeting transcription factors (TFs) has emerged as an important mechanism for gene expression regulation. Here, we tested the hypothesis that in TF 3'UTRs, human-specific single nucleotide change(s) that create novel miRNA recognition elements (MREs) contribute to species-specific differences in TF expression. From several potential human-specific TF MREs, one candidate, a member of the Forkhead Box O (FOXO) subclass in the Forkhead family known as Forkhead Box O1 (FOXO1; FKHR; NM_002015) was tested further. Human FOXO1 contains two sites predicted to confer miR-183-mediated post-transcriptional regulation: one specific to humans and the other conserved. Utilizing dual luciferase expression reporters, we show that only the human FOXO1 3'UTR contains a functional miR-183 site, not found in chimpanzee or mouse 3'untranslated regions (UTRs). Site-directed mutagenesis supports functionality of the human-specific miR-183 site, but not the conserved miR-183 site. Via overexpression and target site protection assays, we show that human FOXO1 is regulated by miR-183, but mouse FOXO1 is not. Finally, FOXO1-regulated cellular phenotypes, including cell invasion and proliferation, are impacted by miR-183 targeting only in human cells. These results provide strong evidence for human-specific gain of TF MREs, a process that may underlie evolutionary differences between phylogenic groups.

Characterization of regulatory features of housekeeping and tissue-specific regulators within tissue regulatory networks

Background

Transcription factors (TFs) and miRNAs are essential for the regulation of gene expression; however, the global view of human gene regulatory networks remains poorly understood. For example, how is the expression of so many genes regulated by limited cohorts of regulators and how are genes differentially expressed in different tissues despite the genetic code being the same in all tissues?

Results

We analyzed the network properties of housekeeping and tissue-specific genes in gene regulatory networks from seven human tissues. Our results show that different classes of genes behave quite differently in these networks. Tissue-specific miRNAs show a higher average target number compared with non-tissue specific miRNAs, which indicates that tissue-specific miRNAs tend to regulate different sets of targets. Tissue-specific TFs exhibit higher in-degree, out-degree, cluster coefficient and betweenness values, indicating that they occupy central positions in the regulatory network and that they transfer genetic information from upstream genes to downstream genes more quickly than other TFs. Housekeeping TFs tend to have higher cluster coefficients compared with other genes that are neither housekeeping nor tissue specific, indicating that housekeeping TFs tend to regulate their targets synergistically. Several topological properties of disease-associated miRNAs and genes were found to be significantly different from those of non-disease-associated miRNAs and genes.

Conclusions

Tissue-specific miRNAs, TFs and disease genes have particular topological properties within the transcriptional regulatory networks of the seven human tissues examined. The tendency of tissue-specific miRNAs to regulate different sets of genes shows that a particular tissue-specific miRNA and its target gene set may form a regulatory module to execute particular functions in the process of tissue differentiation. The regulatory patterns of tissue-specific TFs reflect their vital role in regulatory networks and their importance to biological functions in their respective tissues. The topological differences between disease and non-disease genes may aid the discovery of new disease genes or drug targets. Determining the network properties of these regulatory factors will help define the basic principles of human gene regulation and the molecular mechanisms of disease.

Reducing the complexity of complex gene coexpression networks by coupling multiweighted labeling with topological analysis

Undirected gene coexpression networks obtained from experimental expression data coupled with efficient computational procedures are increasingly used to identify potentially relevant biological information (e.g., biomarkers) for a particular disease. However, coexpression networks built from experimental expression data are in general large highly connected networks with an elevated number of false-positive interactions (nodes and edges). In order to infer relevant information, the network must be properly filtered and its complexity reduced. Given the complexity and the multivariate nature of the information contained in the network, this requires the development and application of efficient feature selection algorithms to be able to exploit the topological characteristics of the network to identify relevant nodes and edges. This paper proposes an efficient multivariate filtering designed to analyze the topological properties of a coexpression network in order to identify potential relevant genes for a given disease. The algorithm has been tested on three datasets for three well known and studied diseases: acute myeloid leukemia, breast cancer, and diffuse large B-cell lymphoma. Results have been validated resorting to bibliographic data automatically mined using the ProteinQuest literature mining tool.

TFClass: an expandable hierarchical classification of human transcription factors

TFClass (http://tfclass.bioinf.med.uni-goettingen.de/) provides a comprehensive classification of human transcription factors based on their DNA-binding domains. Transcription factors constitute a large functional family of proteins directly regulating the activity of genes. Most of them are sequence-specific DNA-binding proteins, thus reading out the information encoded in cis-regulatory DNA elements of promoters, enhancers and other regulatory regions of a genome. TFClass is a database that classifies human transcription factors by a six-level classification schema, four of which are abstractions according to different criteria, while the fifth level represents TF genes and the sixth individual gene products. Altogether, nine superclasses have been identified, comprising 40 classes and 111 families. Counted by genes, 1558 human TFs have been classified so far or >2900 different TFs when including their isoforms generated by alternative splicing or protein processing events. With this classification, we hope to provide a basis for deciphering protein-DNA recognition codes; moreover, it can be used for constructing expanded transcriptional networks by inferring additional TF-target gene relations.

Using potential master regulator sites and paralogous expansion to construct tissue-specific transcriptional networks

Background

Transcriptional networks of higher eukaryotes are difficult to obtain. Available experimental data from conventional approaches are sporadic, while those generated with modern high-throughput technologies are biased. Computational predictions are generally perceived as being flooded with high rates of false positives. New concepts about the structure of regulatory regions and the function of master regulator sites may provide a way out of this dilemma.

Methods

We combined promoter scanning with positional weight matrices with a 4-genome conservativity analysis to predict high-affinity, highly conserved transcription factor (TF) binding sites and to infer TF-target gene relations. They were expanded to paralogous TFs and filtered for tissue-specific expression patterns to obtain a reference transcriptional network (RTN) as well as tissue-specific transcriptional networks (TTNs).

Results

When validated with experimental data sets, the predictions done showed the expected trends of true positive and true negative predictions, resulting in satisfying sensitivity and specificity characteristics. This also proved that confining the network reconstruction to the 1% top-ranking TF-target predictions gives rise to networks with expected degree distributions. Their expansion to paralogous TFs enriches them by tissue-specific regulators, providing a reasonable basis to reconstruct tissue-specific transcriptional networks.

Conclusions

The concept of master regulator or seed sites provides a reasonable starting point to select predicted TF-target relations, which, together with a paralogous expansion, allow for reconstruction of tissue-specific transcriptional networks.