Genome-wide Gene Regulatory Networks Research Articles

A systems genomics and genetics approach to identify the genetic regulatory network for lignin content in Brassica napus seeds.

Seed quality traits of oilseed rape, Brassica napus (B. napus), exhibit quantitative inheritance determined by its genetic makeup and the environment via the mediation of a complex genetic architecture of hundreds to thousands of genes. Thus, instead of single gene analysis, network-based systems genomics and genetics approaches that combine genotype, phenotype, and molecular phenotypes offer a promising alternative to uncover this complex genetic architecture. In the current study, systems genetics approaches were used to explore the genetic regulation of lignin traits in B. napus seeds. Four QTL (qLignin_A09_1, qLignin_A09_2, qLignin_A09_3, and qLignin_C08) distributed on two chromosomes were identified for lignin content. The qLignin_A09_2 and qLignin_C08 loci were homologous QTL from the A and C subgenomes, respectively. Genome-wide gene regulatory network analysis identified eighty-three subnetworks (or modules); and three modules with 910 genes in total, were associated with lignin content, which was confirmed by network QTL analysis. eQTL (expression quantitative trait loci) analysis revealed four cis-eQTL genes including lignin and flavonoid pathway genes, cinnamoyl-CoA-reductase (CCR1), and TRANSPARENT TESTA genes TT4, TT6, TT8, as causal genes. The findings validated the power of systems genetics to identify causal regulatory networks and genes underlying complex traits. Moreover, this information may enable the research community to explore new breeding strategies, such as network selection or gene engineering, to rewire networks to develop climate resilience crops with better seed quality.

Abstract 6217: Differential genome-wide gene regulatory networks for lung cancer

Abstract Background: Lung cancer remains the leading cause of cancer death worldwide due to the fact that most cases are diagnosed at distant stages, and treatment is less effective. During the past two decades, much effort has been devoted to understanding tumor biology and developing targeted therapeutics, which have significantly improved patient survival. With the increasingly available multi-omics data from healthy controls and lung cancer patients, we constructed gene regulatory networks for the two main subtypes of non-small-cell lung cancer (NSCLC), the most common type of lung cancer. In addition, differential networks were constructed by integrative analysis of data from lung cancer tissue and tissues from healthy lungs. Methods: Transcriptomics and genomic data were obtained from the patients with lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC) tissue available in The Cancer Genome Atlas (TCGA) project, and both types of omics data of normal lung tissue were downloaded from The Genotype-Tissue Expression (GTEx) consortium. The networks for each subtype and the differential networks were constructed by utilizing our newly developed machine learning tools, two-stage penalized least square (2SPLS) approach, and analysis of variance of directed networks (NetANOVA), respectively. Results: For the top three identified subnetworks of LUAD, we identified some well-known genes associated with lung cancer, such as EGFR, IRF1, KRAS, ERBB2, PIK3CA, MYC, at bootstrap frequency >90%. Based on the top three subnetworks with a bootstrap frequency of 100%, results from the Ingenuity Pathway Analysis (IPA) indicated several significant pathways, primary immunodeficiency signaling and communication between innate and adaptive immune cells. GO enrichment analysis showed statistically significant biological processes, including immune response, extracellular matrix organization, and significant KEGG pathways such as Hematopoietic cell lineage. The differential networks between LUAD, LUSC, and healthy lungs showed the same regulations across all three groups, such as the regulation between EEF2 and RACK1, the perturbations between LUAD and healthy (two-way regulations between HIGD2A and ARL10), and LUSC and healthy (GRB14 regulates EIF3E) respectively. In addition, we identified the differences between LUAD and LUSC, such as bi-directional regulations between CLTB and ARL10. Conclusion: Although all groups share most regulations, we identified differential gene regulations between healthy lung tissue and lung cancer tissue and between LUAD and LUSC of NSCLC. Using multi-omics data generated from lung tissues, coupled with advanced machine learning methods, the causal relationships between genome-wide genes from the analysis will help us understand the molecular mechanism of lung cancer and facilitate the development of personalized treatment strategies. Citation Format: Min Zhang, Zhongli Jiang, Dabao Zhang. Differential genome-wide gene regulatory networks for lung cancer [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6217.

In silico identification of papaya genome-encoded microRNAs to target begomovirus genes in papaya leaf curl disease.

Papaya leaf curl disease (PaLCuD) is widespread and classified in the genus begomovirus (Geminiviridae), disseminated by the vector whitefly Bemisia tabaci. RNA interference (RNAi)-based antiviral innate immunity stands as a pivotal defense mechanism and biological process in limiting viral genomes to manage plant diseases. The current study aims to identify and analyze Carica Papaya locus-derived capa-microRNAs with predicted potential for targeting divergent begomovirus species-encoded mRNAs using a 'four integrative in silico algorithms' approach. This research aims to experimentally activate the RNAi catalytic pathway using in silico-predicted endogenous capa-miRNAs and create papaya varieties capable of assessing potential resistance against begomovirus species and monitoring antiviral capabilities. This study identified 48 predicted papaya locus-derived candidates from 23 miRNA families, which were further investigated for targeting begomovirus genes. Premised all the four algorithms combined, capa-miR5021 was the most anticipated miRNA followed by capa-miR482, capa-miR5658, capa-miR530b, capa-miR3441.2, and capa-miR414 'effective' papaya locus-derived candidate capa-miRNA and respected putative binding sites for targets at the consensus nucleotide position. It was predicted to bind and target mostly to AC1 gene of the complementary strand and the AV1 gene of the virion strand of different begomovirus isolates, which were associated with replication-associated protein and encapsidation, respectively, during PaLCuD. These miRNAs were also found targeting betaC1 gene of betasatellite which were associated with retardation in leaf growth and developmental abnormalities with severe symptoms during begomovirus infection. To validate target prediction accuracy, we created an integrated Circos plot for comprehensive visualization of host-virus interaction. In silico-predicted papaya genome-wide miRNA-mediated begomovirus target gene regulatory network corroborated interactions that permit in vivo analysis, which could provide biological material and valuable evidence, leading to the development of begomovirus-resistant papaya plants. The integrative nature of our research positions it at the forefront of efforts to ensure the sustainable cultivation of papaya, particularly in the face of evolving pathogenic threats. As we move forward, the knowledge gained from this study provides a solid foundation for continued exploration and innovation in the field of papaya virology, and to the best of our knowledge, this study represents a groundbreaking endeavor, undertaken for the first time in the context of PaLCuD research.

Predicting proximal tubule failed repair drivers through regularized regression analysis of single cell multiomic sequencing

Renal proximal tubule epithelial cells have considerable intrinsic repair capacity following injury. However, a fraction of injured proximal tubule cells fails to undergo normal repair and assumes a proinflammatory and profibrotic phenotype that may promote fibrosis and chronic kidney disease. The healthy to failed repair change is marked by cell state-specific transcriptomic and epigenomic changes. Single nucleus joint RNA- and ATAC-seq sequencing offers an opportunity to study the gene regulatory networks underpinning these changes in order to identify key regulatory drivers. We develop a regularized regression approach to construct genome-wide parametric gene regulatory networks using multiomic datasets. We generate a single nucleus multiomic dataset from seven adult human kidney samples and apply our method to study drivers of a failed injury response associated with kidney disease. We demonstrate that our approach is a highly effective tool for predicting key cis- and trans-regulatory elements underpinning the healthy to failed repair transition and use it to identify NFAT5 as a driver of the maladaptive proximal tubule state.

Heterogeneity in the gene regulatory landscape of leiomyosarcoma.

Characterizing inter-tumor heterogeneity is crucial for selecting suitable cancer therapy, as the presence of diverse molecular subgroups of patients can be associated with disease outcome or response to treatment. While cancer subtypes are often characterized by differences in gene expression, the mechanisms driving these differences are generally unknown. We set out to model the regulatory mechanisms driving sarcoma heterogeneity based on patient-specific, genome-wide gene regulatory networks. We developed a new computational framework, PORCUPINE, which combines knowledge on biological pathways with permutation-based network analysis to identify pathways that exhibit significant regulatory heterogeneity across a patient population. We applied PORCUPINE to patient-specific leiomyosarcoma networks modeled on data from The Cancer Genome Atlas and validated our results in an independent dataset from the German Cancer Research Center. PORCUPINE identified 37 heterogeneously regulated pathways, including pathways representing potential targets for treatment of subgroups of leiomyosarcoma patients, such as FGFR and CTLA4 inhibitory signaling. We validated the detected regulatory heterogeneity through analysis of networks and chromatin states in leiomyosarcoma cell lines. We showed that the heterogeneity identified with PORCUPINE is not associated with methylation profiles or clinical features, thereby suggesting an independent mechanism of patient heterogeneity driven by the complex landscape of gene regulatory interactions.

Genetic variation and domestication of horses revealed by 10 chromosome-level genomes and whole-genome resequencing.

Understanding the genetic variations of the horse (Equus caballus) genome will improve breeding conservation and welfare. However, genetic variations in long segments, such as structural variants (SVs), remain understudied. We de novo assembled 10 chromosome-level three-dimensional horse genomes, each representing a distinct breed, and analysed horse SVs using a multi-assembly approach. Our findings suggest that SVs with the accumulation of mammalian-wide interspersed repeats related to long interspersed nuclear elements might be a horse-specific mechanism to modulate genome-wide gene regulatory networks. We found that olfactory receptors were commonly loss and accumulated deleterious mutations, but no purge of deleterious mutations occurred during horse domestication. We examined the potential effects of SVs on the spatial structure of chromatin via topologically associating domains (TADs). Breed-specific TADs were significantly enriched by breed-specific SVs. We identified 4199 unique breakpoint-resolved novel insertions across all chromosomes that account for 2.84 Mb sequences missing from the reference genome. Several novel insertions might have potential functional consequences, as 519 appeared to reside within 449 gene bodies. These genes are primarily involved in pathogen recognition, innate immune responses and drug metabolism. Moreover, 37 diverse horses were resequenced. Combining this with public data, we analysed 97 horses through a comparative population genomics approach to identify the genetic basis underlying breed characteristics using Thoroughbreds as a case study. We provide new scientific evidence for horse domestication, an understanding of the genetic mechanism underlying the phenotypic evolution of horses, and a comprehensive genetic variation resource for further genetic studies of horses.

Genome-Wide Transcription Factor DNA Binding Sites and Gene Regulatory Networks in Clostridium thermocellum.

Clostridium thermocellum is a thermophilic bacterium recognized for its natural ability to effectively deconstruct cellulosic biomass. While there is a large body of studies on the genetic engineering of this bacterium and its physiology to-date, there is limited knowledge in the transcriptional regulation in this organism and thermophilic bacteria in general. The study herein is the first report of a large-scale application of DNA-affinity purification sequencing (DAP-seq) to transcription factors (TFs) from a bacterium. We applied DAP-seq to > 90 TFs in C. thermocellum and detected genome-wide binding sites for 11 of them. We then compiled and aligned DNA binding sequences from these TFs to deduce the primary DNA-binding sequence motifs for each TF. These binding motifs are further validated with electrophoretic mobility shift assay (EMSA) and are used to identify individual TFs’ regulatory targets in C. thermocellum. Our results led to the discovery of novel, uncharacterized TFs as well as homologues of previously studied TFs including RexA-, LexA-, and LacI-type TFs. We then used these data to reconstruct gene regulatory networks for the 11 TFs individually, which resulted in a global network encompassing the TFs with some interconnections. As gene regulation governs and constrains how bacteria behave, our findings shed light on the roles of TFs delineated by their regulons, and potentially provides a means to enable rational, advanced genetic engineering of C. thermocellum and other organisms alike toward a desired phenotype.

Genome-environment association methods comparison supports omnigenic adaptation to ecological niche in malaria vector mosquitoes.

The concept of a fundamental ecological niche is central to questions of geographic distribution, population demography, species conservation, and evolutionary potential. However, robust inference of genomic regions associated with evolutionary adaptation to particular environmental conditions remains difficult due to the myriad of potential confounding processes that can generate heterogeneous patterns of variation across the genome. Here, we interrogate the potential role of genome environment association (GEA) testing as an initial step in building an understanding of the genetic basis of ecological niche. We leverage publicly available genomic data from the Anopheles gambiae 1000 Genomes (Ag1000g) Consortium to test the ability of multiple analytically unique GEA methods to handle confounding patterns of genetic variation, control false positive rates, and discern associations with broadly relevant climate variables from random allele frequency patterns throughout the genome. We found evidence supporting the ability of commonly implemented GEA methods to account for confounding patterns of spatial and genetic variation, and control false positive rates. However, we fail to find evidence supporting the ability of GEA tests to reject signals of adaptation to randomly simulated environmental variables, indicating that discerning between true signals of genome environment adaptation and genome environment correlations resulting from alternative evolutionary processes, remains challenging. Because signals of environmental adaptation are so diffuse and confounded throughout the genome, we argue that genomic adaptation to ecological niche is likely best understood under an omnigenic model wherein highly interconnected, genome-wide gene regulatory networks shape genomic adaptation to key environmental conditions.

Inference of Bacterial Small RNA Regulatory Networks and Integration with Transcription Factor-Driven Regulatory Networks.

ABSTRACTSmall noncoding RNAs (sRNAs) are key regulators of bacterial gene expression. Through complementary base pairing, sRNAs affect mRNA stability and translation efficiency. Here, we describe a network inference approach designed to identify sRNA-mediated regulation of transcript levels. We use existing transcriptional data sets and prior knowledge to infer sRNA regulons using our network inference tool, the Inferelator. This approach produces genome-wide gene regulatory networks that include contributions by both transcription factors and sRNAs. We show the benefits of estimating and incorporating sRNA activities into network inference pipelines using available experimental data. We also demonstrate how these estimated sRNA regulatory activities can be mined to identify the experimental conditions where sRNAs are most active. We uncover 45 novel experimentally supported sRNA-mRNA interactions in Escherichia coli, outperforming previous network-based efforts. Additionally, our pipeline complements sequence-based sRNA-mRNA interaction prediction methods by adding a data-driven filtering step. Finally, we show the general applicability of our approach by identifying 24 novel, experimentally supported, sRNA-mRNA interactions in Pseudomonas aeruginosa, Staphylococcus aureus, and Bacillus subtilis. Overall, our strategy generates novel insights into the functional context of sRNA regulation in multiple bacterial species.IMPORTANCE Individual bacterial genomes can have dozens of small noncoding RNAs with largely unexplored regulatory functions. Although bacterial sRNAs influence a wide range of biological processes, including antibiotic resistance and pathogenicity, our current understanding of sRNA-mediated regulation is far from complete. Most of the available information is restricted to a few well-studied bacterial species; and even in those species, only partial sets of sRNA targets have been characterized in detail. To close this information gap, we developed a computational strategy that takes advantage of available transcriptional data and knowledge about validated and putative sRNA-mRNA interactions for inferring expanded sRNA regulons. Our approach facilitates the identification of experimentally supported novel interactions while filtering out false-positive results. Due to its data-driven nature, our method prioritizes biologically relevant interactions among lists of candidate sRNA-target pairs predicted in silico from sequence analysis or derived from sRNA-mRNA binding experiments.

Supervised Learning of Gene Regulatory Networks.

Identifying the entirety of gene regulatory interactions in a biological system offers the possibility to determine the key molecular factors that affect important traits on the level of cells, tissues, and whole organisms. Despite the development of experimental approaches and technologies for identification of direct binding of transcription factors (TFs) to promoter regions of downstream target genes, computational approaches that utilize large compendia of transcriptomics data are still the predominant methods used to predict direct downstream targets of TFs, and thus reconstruct genome-wide gene-regulatory networks (GRNs). These approaches can broadly be categorized into unsupervised and supervised, based on whether data about known, experimentally verified gene-regulatory interactions are used in the process of reconstructing the underlying GRN. Here, we first describe the generic steps of supervised approaches for GRN reconstruction, since they have been recently shown to result in improved accuracy of the resulting networks? We also illustrate how they can be used with data from model organisms to obtain more accurate prediction of gene regulatory interactions. © 2020 The Authors. Basic Protocol 1: Construction of features used in supervised learning of gene regulatory interactions Basic Protocol 2: Learning the non-interacting TF-gene pairs Basic Protocol 3: Learning a classifier for gene regulatory interactions.

Network-based analysis of prostate cancer cell lines reveals novel marker gene candidates associated with radioresistance and patient relapse.

Radiation therapy is an important and effective treatment option for prostate cancer, but high-risk patients are prone to relapse due to radioresistance of cancer cells. Molecular mechanisms that contribute to radioresistance are not fully understood. Novel computational strategies are needed to identify radioresistance driver genes from hundreds of gene copy number alterations. We developed a network-based approach based on lasso regression in combination with network propagation for the analysis of prostate cancer cell lines with acquired radioresistance to identify clinically relevant marker genes associated with radioresistance in prostate cancer patients. We analyzed established radioresistant cell lines of the prostate cancer cell lines DU145 and LNCaP and compared their gene copy number and expression profiles to their radiosensitive parental cells. We found that radioresistant DU145 showed much more gene copy number alterations than LNCaP and their gene expression profiles were highly cell line specific. We learned a genome-wide prostate cancer-specific gene regulatory network and quantified impacts of differentially expressed genes with directly underlying copy number alterations on known radioresistance marker genes. This revealed several potential driver candidates involved in the regulation of cancer-relevant processes. Importantly, we found that ten driver candidates from DU145 (ADAMTS9, AKR1B10, CXXC5, FST, FOXL1, GRPR, ITGA2, SOX17, STARD4, VGF) and four from LNCaP (FHL5, LYPLAL1, PAK7, TDRD6) were able to distinguish irradiated prostate cancer patients into early and late relapse groups. Moreover, in-depth in vitro validations for VGF (Neurosecretory protein VGF) showed that siRNA-mediated gene silencing increased the radiosensitivity of DU145 and LNCaP cells. Our computational approach enabled to predict novel radioresistance driver gene candidates. Additional preclinical and clinical studies are required to further validate the role of VGF and other candidate genes as potential biomarkers for the prediction of radiotherapy responses and as potential targets for radiosensitization of prostate cancer.

Hierarchical graphical model reveals HFR1 bridging circadian rhythm and flower development in Arabidopsis thaliana

sTo study systems-level properties of the cell, it is necessary to go beyond individual regulators and target genes to study the regulatory network among transcription factors (TFs). However, it is difficult to directly dissect the TFs mediated genome-wide gene regulatory network (GRN) by experiment. Here, we proposed a hierarchical graphical model to estimate TF activity from mRNA expression by building TF complexes with protein cofactors and inferring TF’s downstream regulatory network simultaneously. Then we applied our model on flower development and circadian rhythm processes in Arabidopsis thaliana. The computational results show that the sequence specific bHLH family TF HFR1 recruits the chromatin regulator HAC1 to flower development master regulator TF AG and further activates AG’s expression by histone acetylation. Both independent data and experimental results supported this discovery. We also found a flower tissue specific H3K27ac ChIP-seq peak at AG gene body and a HFR1 motif in the center of this H3K27ac peak. Furthermore, we verified that HFR1 physically interacts with HAC1 by yeast two-hybrid experiment. This HFR1–HAC1–AG triplet relationship may imply that flower development and circadian rhythm are bridged by epigenetic regulation and enrich the classical ABC model in flower development. In addition, our TF activity network can serve as a general method to elucidate molecular mechanisms on other complex biological regulatory processes.

A scalable algorithm for structure identification of complex gene regulatory network from temporal expression data

BackgroundGene regulatory interactions are of fundamental importance to various biological functions and processes. However, only a few previous computational studies have claimed success in revealing genome-wide regulatory landscapes from temporal gene expression data, especially for complex eukaryotes like human. Moreover, recent work suggests that these methods still suffer from the curse of dimensionality if a network size increases to 100 or higher.ResultsHere we present a novel scalable algorithm for identifying genome-wide gene regulatory network (GRN) structures, and we have verified the algorithm performances by extensive simulation studies based on the DREAM challenge benchmark data. The highlight of our method is that its superior performance does not degenerate even for a network size on the order of 104, and is thus readily applicable to large-scale complex networks. Such a breakthrough is achieved by considering both prior biological knowledge and multiple topological properties (i.e., sparsity and hub gene structure) of complex networks in the regularized formulation. We also validate and illustrate the application of our algorithm in practice using the time-course gene expression data from a study on human respiratory epithelial cells in response to influenza A virus (IAV) infection, as well as the CHIP-seq data from ENCODE on transcription factor (TF) and target gene interactions. An interesting finding, owing to the proposed algorithm, is that the biggest hub structures (e.g., top ten) in the GRN all center at some transcription factors in the context of epithelial cell infection by IAV.ConclusionsThe proposed algorithm is the first scalable method for large complex network structure identification. The GRN structure identified by our algorithm could reveal possible biological links and help researchers to choose which gene functions to investigate in a biological event. The algorithm described in this article is implemented in MATLAB Ⓡ, and the source code is freely available from https://github.com/Hongyu-Miao/DMI.git.

Model-based transcriptome engineering promotes a fermentative transcriptional state in yeast

The ability to rationally manipulate the transcriptional states of cells would be of great use in medicine and bioengineering. We have developed an algorithm, NetSurgeon, which uses genome-wide gene-regulatory networks to identify interventions that force a cell toward a desired expression state. We first validated NetSurgeon extensively on existing datasets. Next, we used NetSurgeon to select transcription factor deletions aimed at improving ethanol production in Saccharomyces cerevisiae cultures that are catabolizing xylose. We reasoned that interventions that move the transcriptional state of cells using xylose toward that of cells producing large amounts of ethanol from glucose might improve xylose fermentation. Some of the interventions selected by NetSurgeon successfully promoted a fermentative transcriptional state in the absence of glucose, resulting in strains with a 2.7-fold increase in xylose import rates, a 4-fold improvement in xylose integration into central carbon metabolism, or a 1.3-fold increase in ethanol production rate. We conclude by presenting an integrated model of transcriptional regulation and metabolic flux that will enable future efforts aimed at improving xylose fermentation to prioritize functional regulators of central carbon metabolism.

ARMADA: Using motif activity dynamics to infer gene regulatory networks from gene expression data

Analysis of gene expression data remains one of the most promising avenues toward reconstructing genome-wide gene regulatory networks. However, the large dimensionality of the problem prohibits the fitting of explicit dynamical models of gene regulatory networks, whereas machine learning methods for dimensionality reduction such as clustering or principal component analysis typically fail to provide mechanistic interpretations of the reduced descriptions. To address this, we recently developed a general methodology called motif activity response analysis (MARA) that, by modeling gene expression patterns in terms of the activities of concrete regulators, accomplishes dramatic dimensionality reduction while retaining mechanistic biological interpretations of its predictions (Balwierz, 2014).Here we extend MARA by presenting ARMADA, which models the activity dynamics of regulators across a time course, and infers the causal interactions between the regulators that drive the dynamics of their activities across time. We have implemented ARMADA as part of our ISMARA webserver, ismara.unibas.ch, allowing any researcher to automatically apply it to any gene expression time course. To illustrate the method, we apply ARMADA to a time course of human umbilical vein endothelial cells treated with TNF. Remarkably, ARMADA is able to reproduce the complex observed motif activity dynamics using a relatively small set of interactions between the key regulators in this system. In addition, we show that ARMADA successfully infers many of the key regulatory interactions known to drive this inflammatory response and discuss several novel interactions that ARMADA predicts. In combination with ISMARA, ARMADA provides a powerful approach to generating plausible hypotheses for the key interactions between regulators that control gene expression in any system for which time course measurements are available.

Urothelial cancer gene regulatory networks inferred from large-scale RNAseq, Bead and Oligo gene expression data.

BackgroundUrothelial pathogenesis is a complex process driven by an underlying network of interconnected genes. The identification of novel genomic target regions and gene targets that drive urothelial carcinogenesis is crucial in order to improve our current limited understanding of urothelial cancer (UC) on the molecular level. The inference of genome-wide gene regulatory networks (GRN) from large-scale gene expression data provides a promising approach for a detailed investigation of the underlying network structure associated to urothelial carcinogenesis.MethodsIn our study we inferred and compared three GRNs by the application of the BC3Net inference algorithm to large-scale transitional cell carcinoma gene expression data sets from Illumina RNAseq (179 samples), Illumina Bead arrays (165 samples) and Affymetrix Oligo microarrays (188 samples). We investigated the structural and functional properties of GRNs for the identification of molecular targets associated to urothelial cancer.ResultsWe found that the urothelial cancer (UC) GRNs show a significant enrichment of subnetworks that are associated with known cancer hallmarks including cell cycle, immune response, signaling, differentiation and translation. Interestingly, the most prominent subnetworks of co-located genes were found on chromosome regions 5q31.3 (RNAseq), 8q24.3 (Oligo) and 1q23.3 (Bead), which all represent known genomic regions frequently deregulated or aberated in urothelial cancer and other cancer types. Furthermore, the identified hub genes of the individual GRNs, e.g., HID1/DMC1 (tumor development), RNF17/TDRD4 (cancer antigen) and CYP4A11 (angiogenesis/ metastasis) are known cancer associated markers. The GRNs were highly dataset specific on the interaction level between individual genes, but showed large similarities on the biological function level represented by subnetworks. Remarkably, the RNAseq UC GRN showed twice the proportion of significant functional subnetworks. Based on our analysis of inferential and experimental networks the Bead UC GRN showed the lowest performance compared to the RNAseq and Oligo UC GRNs.ConclusionTo our knowledge, this is the first study investigating genome-scale UC GRNs. RNAseq based gene expression data is the data platform of choice for a GRN inference. Our study offers new avenues for the identification of novel putative diagnostic targets for subsequent studies in bladder tumors.Electronic supplementary materialThe online version of this article (doi:10.1186/s12918-015-0165-z) contains supplementary material, which is available to authorized users.

Reverse Engineering of Genome-wide Gene Regulatory Networks from Gene Expression Data.

Transcriptional regulation plays vital roles in many fundamental biological processes. Reverse engineering of genome-wide regulatory networks from high-throughput transcriptomic data provides a promising way to characterize the global scenario of regulatory relationships between regulators and their targets. In this review, we summarize and categorize the main frameworks and methods currently available for inferring transcriptional regulatory networks from microarray gene expression profiling data. We overview each of strategies and introduce representative methods respectively. Their assumptions, advantages, shortcomings, and possible improvements and extensions are also clarified and commented.

Genome-wide enhanced gene regulatory networks

Genome-wide enhanced gene regulatory networks

Reverse Engineering and Analysis of Genome-Wide Gene Regulatory Networks from Gene Expression Profiles Using High-Performance Computing

Regulation of gene expression is a carefully regulated phenomenon in the cell. “Reverse-engineering” algorithms try to reconstruct the regulatory interactions among genes from genome-scale measurements of gene expression profiles (microarrays). Mammalian cells express tens of thousands of genes; hence, hundreds of gene expression profiles are necessary in order to have acceptable statistical evidence of interactions between genes. As the number of profiles to be analyzed increases, so do computational costs and memory requirements. In this work, we designed and developed a parallel computing algorithm to reverse-engineer genome-scale gene regulatory networks from thousands of gene expression profiles. The algorithm is based on computing pairwise Mutual Information between each gene-pair. We successfully tested it to reverse engineer the Mus Musculus (mouse) gene regulatory network in liver from gene expression profiles collected from a public repository. A parallel hierarchical clustering algorithm was implemented to discover “communities” within the gene network. Network communities are enriched for genes involved in the same biological functions. The inferred network was used to identify two mitochondrial proteins.

An orphan sensor kinase controls quinolone signal production via MexT in Pseudomonas aeruginosa

Pseudomonas aeruginosa employs both N-acylhomoserine lactone and 2-alkyl-4(1H)-quinolone (AQ)-mediated interbacterial signalling for the orchestration of a genome-wide gene regulatory network. Despite the many advances that have been made in understanding the target genes of quorum sensing regulation, little is known on how quorum sensing systems are influenced by environmental cues. In this study, we show that AQ production is modulated by an orphan P. aeruginosa sensor kinase. Transcriptional studies of the sensor kinase (MxtR) mutant demonstrated that an induced expression of MexT, a LysR-type transcriptional regulator, largely determined the global transcriptional profile. Thereby, overexpression of the MexT-regulated MexEF-OprN efflux pump led to a delayed expression of the AQ biosynthetic genes and of AQ-dependent virulence factors. Furthermore, we demonstrated that autophosphorylation of MxtR was inhibited by ubiquinone, the central electron carrier of respiration in in vitro experiments. Our results elucidate on a mechanism by which P. aeruginosa senses environmental conditions and adapts by controlling the production of interbacterial AQ signal molecules. A regulatory function of a sensor kinase may indicate that there is a pre-emptive role of adaptation mechanisms that are turned on under distinct environmental conditions and that are important for efficient colonization and pathogenesis.