scholarly journals Deciphering the genetic links between NAFLD and co-occurring conditions using a liver gene regulatory network

2021 ◽  
Author(s):  
Sreemol Gokuladhas ◽  
William Schierding ◽  
Roan Eltigani Zaied ◽  
Tayaza Fadason ◽  
Murim Choi ◽  
...  

Background & Aims: Non-alcoholic fatty liver disease (NAFLD) is a multi-system metabolic disease that co-occurs with various hepatic and extra-hepatic diseases. The phenotypic manifestation of NAFLD is primarily observed in the liver. Therefore, identifying liver-specific gene regulatory interactions between variants associated with NAFLD and multimorbid conditions may help to improve our understanding of underlying shared aetiology. Methods: Here, we constructed a liver-specific gene regulatory network (LGRN) consisting of genome-wide spatially constrained expression quantitative trait loci (eQTLs) and their target genes. The LGRN was used to identify regulatory interactions involving NAFLD-associated genetic modifiers and their inter-relationships to other complex traits. Results and Conclusions: We demonstrate that MBOAT7 and IL32, which are associated with NAFLD progression, are regulated by spatially constrained eQTLs that are enriched for an association with liver enzyme levels. MBOAT7 transcript levels are also linked to eQTLs associated with cirrhosis, and other traits that commonly co-occur with NAFLD. In addition, genes that encode interacting partners of NAFLD-candidate genes within the liver-specific protein-protein interaction network were affected by eQTLs enriched for phenotypes relevant to NAFLD (e.g. IgG glycosylation patterns, OSA). Furthermore, we identified distinct gene regulatory networks formed by the NAFLD-associated eQTLs in normal versus diseased liver, consistent with the context-specificity of the eQTLs effects. Interestingly, genes targeted by NAFLD-associated eQTLs within the LGRN were also affected by eQTLs associated with NAFLD-related traits (e.g. obesity and body fat percentage). Overall, the genetic links identified between these traits expand our understanding of shared regulatory mechanisms underlying NAFLD multimorbidities.

2021 ◽  
Author(s):  
Xiangyu Pan ◽  
Zhaoxia Ma ◽  
Xinqi Sun ◽  
Hui Li ◽  
Tingting Zhang ◽  
...  

Biologists long recognized that the genetic information encoded in DNA leads to trait innovation via gene regulatory network (GRN) in development. Here, we generated paired expression and chromatin accessibility data during rumen and esophagus development in sheep and revealed 1,601 active ruminant-specific conserved non-coding elements (active-RSCNEs). To interpret the function of these active-RSCNEs, we developed a Conserved Non-coding Element interpretation method by gene Regulatory network (CNEReg) to define toolkit transcription factors (TTF) and model its regulation on rumen specific gene via batteries of active-RSCNEs during development. Our developmental GRN reveals 18 TTFs and 313 active-RSCNEs regulating the functional modules of the rumen and identifies OTX1, SOX21, HOXC8, SOX2, TP63, PPARG and 16 active-RSCNEs that functionally distinguish the rumen from the esophagus. We argue that CNEReg is an attractive systematic approach to integrate evo-devo concepts with omics data to understand how gene regulation evolves and shapes complex traits.


2021 ◽  
Author(s):  
Marouen Ben Guebila ◽  
Camila M Lopes-Ramos ◽  
Deborah Weighill ◽  
Abhijeet Rajendra Sonawane ◽  
Rebekka Burkholz ◽  
...  

Abstract Gene regulation plays a fundamental role in shaping tissue identity, function, and response to perturbation. Regulatory processes are controlled by complex networks of interacting elements, including transcription factors, miRNAs and their target genes. The structure of these networks helps to determine phenotypes and can ultimately influence the development of disease or response to therapy. We developed GRAND (https://grand.networkmedicine.org) as a database for computationally-inferred, context-specific gene regulatory network models that can be compared between biological states, or used to predict which drugs produce changes in regulatory network structure. The database includes 12 468 genome-scale networks covering 36 human tissues, 28 cancers, 1378 unperturbed cell lines, as well as 173 013 TF and gene targeting scores for 2858 small molecule-induced cell line perturbation paired with phenotypic information. GRAND allows the networks to be queried using phenotypic information and visualized using a variety of interactive tools. In addition, it includes a web application that matches disease states to potentially therapeutic small molecule drugs using regulatory network properties.


2019 ◽  
Author(s):  
Dan Ramirez ◽  
Vivek Kohar ◽  
Ataur Katebi ◽  
Mingyang Lu

AbstractEpithelial-mesenchymal transition (EMT) plays a crucial role in embryonic development and tumorigenesis. Although EMT has been extensively studied with both computational and experimental methods, the gene regulatory mechanisms governing the transition are not yet well understood. Recent investigations have begun to better characterize the complex phenotypic plasticity underlying EMT using a computational systems biology approach. Here, we analyzed recently published single-cell RNA sequencing data from E9.5 to E11.5 mouse embryonic skin cells and identified the gene expression patterns of both epithelial and mesenchymal phenotypes, as well as a clear hybrid state. By integrating the scRNA-seq data and gene regulatory interactions from the literature, we constructed a gene regulatory network model governing the decision-making of EMT in the context of the developing mouse embryo. We simulated the network using a recently developed mathematical modeling method, named RACIPE, and observed three distinct phenotypic states whose gene expression patterns can be associated with the epithelial, hybrid, and mesenchymal states in the scRNA-seq data. Additionally, the model is in agreement with published results on the composition of EMT phenotypes and regulatory networks. We identified Wnt signaling as a major pathway in inducing the EMT and its role in driving cellular state transitions during embryonic development. Our findings demonstrate a new method of identifying and incorporating tissue-specific regulatory interactions into gene regulatory network modeling.Author SummaryEpithelial-mesenchymal transition (EMT) is a cellular process wherein cells become disconnected from their surroundings and acquire the ability to migrate through the body. EMT has been observed in biological contexts including development, wound healing, and cancer, yet the regulatory mechanisms underlying it are not well understood. Of particular interest is a purported hybrid state, in which cells can retain some adhesion to their surroundings but also show mesenchymal traits. Here, we examine the prevalence and composition of the hybrid state in the context of the embryonic mouse, integrating gene regulatory interactions from published experimental results as well as from the specific single cell RNA sequencing dataset of interest. Using mathematical modeling, we simulated a regulatory network based on these sources and aligned the simulated phenotypes with those in the data. We identified a hybrid EMT phenotype and revealed the inducing effect of Wnt signaling on EMT in this context. Our regulatory network construction process can be applied beyond EMT to illuminate the behavior of any biological phenomenon occurring in a specific context, allowing better identification of therapeutic targets and further research directions.


eLife ◽  
2016 ◽  
Vol 5 ◽  
Author(s):  
Erik Clark ◽  
Michael Akam

The Drosophila embryo transiently exhibits a double-segment periodicity, defined by the expression of seven 'pair-rule' genes, each in a pattern of seven stripes. At gastrulation, interactions between the pair-rule genes lead to frequency doubling and the patterning of 14 parasegment boundaries. In contrast to earlier stages of Drosophila anteroposterior patterning, this transition is not well understood. By carefully analysing the spatiotemporal dynamics of pair-rule gene expression, we demonstrate that frequency-doubling is precipitated by multiple coordinated changes to the network of regulatory interactions between the pair-rule genes. We identify the broadly expressed but temporally patterned transcription factor, Odd-paired (Opa/Zic), as the cause of these changes, and show that the patterning of the even-numbered parasegment boundaries relies on Opa-dependent regulatory interactions. Our findings indicate that the pair-rule gene regulatory network has a temporally modulated topology, permitting the pair-rule genes to play stage-specific patterning roles.


2019 ◽  
Author(s):  
Daniel Morgan ◽  
Matthew Studham ◽  
Andreas Tjärnberg ◽  
Holger Weishaupt ◽  
Fredrik J. Swartling ◽  
...  

AbstractThe gene regulatory network (GRN) of human cells encodes mechanisms to ensure proper functioning. However, if this GRN is dysregulated, the cell may enter into a disease state such as cancer. Understanding the GRN as a system can therefore help identify novel mechanisms underlying disease, which can lead to new therapies. Reliable inference of GRNs is however still a major challenge in systems biology.To deduce regulatory interactions relevant to cancer, we applied a recent computational inference framework to data from perturbation experiments in squamous carcinoma cell line A431. GRNs were inferred using several methods, and the false discovery rate was controlled by the NestBoot framework. We developed a novel approach to assess the predictiveness of inferred GRNs against validation data, despite the lack of a gold standard. The best GRN was significantly more predictive than the null model, both in crossvalidated benchmarks and for an independent dataset of the same genes under a different perturbation design. It agrees with many known links, in addition to predicting a large number of novel interactions from which a subset was experimentally validated. The inferred GRN captures regulatory interactions central to cancer-relevant processes and thus provides mechanistic insights that are useful for future cancer research.Data available at GSE125958Inferred GRNs and inference statistics available at https://dcolin.shinyapps.io/CancerGRN/ Software available at https://bitbucket.org/sonnhammergrni/genespider/src/BFECV/Author SummaryCancer is the second most common cause of death globally, and although cancer treatments have improved in recent years, we need to understand how regulatory mechanisms are altered in cancer to combat the disease efficiently. By applying gene perturbations and inference of gene regulatory networks to 40 genes known or suspected to have a role in cancer due to interactions with the oncogene MYC, we deduce their underlying regulatory interactions. Using a recent computational framework for inference together with a novel method for cross validation, we infer a reliable regulatory model of this system in a completely data driven manner, not reliant on literature or priors. The novel interactions add to the understanding of the progressive oncogenic regulatory process and may provide new targets for therapy.


Author(s):  
Xin Zhou ◽  
Xiaodong Cai

Abstract Motivation Genetic variations of expression quantitative trait loci (eQTLs) play a critical role in influencing complex traits and diseases development. Two main factors that affect the statistical power of detecting eQTLs are: 1) relatively small size of samples available, and 2) heavy burden of multiple testing due to a very large number of variants to be tested. The later issue is particularly severe when one tries to identify trans-eQTLs that are far away from the genes they influence. If one can exploit co-expressed genes jointly in eQTL-mapping, effective sample size can be increased. Furthermore, using the structure of the gene regulatory network (GRN) may help to identify trans-eQTLs without increasing multiple testing burden. Results In this paper, we employ the structure equation model (SEM) to model both GRN and effect of eQTLs on gene expression, and then develop a novel algorithm, named sparse SEM for eQTL mapping (SSEMQ), to conduct joint eQTL mapping and GRN inference. The SEM can exploit co-expressed genes jointly in eQTL mapping and also use GRN to determine trans-eQTLs. Computer simulations demonstrate that our SSEMQ significantly outperforms nine existing eQTL mapping methods. SSEMQ is further employed to analyze two real datasets of human breast and whole blood tissues, yielding a number of cis- and trans-eQTLs. Availability R package ssemQr is available at https://github.com/Ivis4ml/ssemQr.git. Supplementary information Supplementary data are available at Bioinformatics online.


2021 ◽  
Author(s):  
Deborah Weighill ◽  
Marouen Ben Guebila ◽  
Kimberly Glass ◽  
John Quackenbush ◽  
John Platig

AbstractThe majority of disease-associated genetic variants are thought to have regulatory effects, including the disruption of transcription factor (TF) binding and the alteration of downstream gene expression. Identifying how a person’s genotype affects their individual gene regulatory network has the potential to provide important insights into disease etiology and to enable improved genotype-specific disease risk assessments and treatments. However, the impact of genetic variants is generally not considered when constructing gene regulatory networks. To address this unmet need, we developed EGRET (Estimating the Genetic Regulatory Effect on TFs), which infers a genotype-specific gene regulatory network (GRN) for each individual in a study population by using message passing to integrate genotype-informed TF motif predictions - derived from individual genotype data, the predicted effects of variants on TF binding and gene expression, and TF motif predictions - with TF protein-protein interactions and gene expression. Comparing EGRET networks for two blood-derived cell lines identified genotype-associated cell-line specific regulatory differences which were subsequently validated using allele-specific expression, chromatin accessibility QTLs, and differential TF binding from ChIP-seq. In addition, EGRET GRNs for three cell types across 119 individuals captured regulatory differences associated with disease in a cell-type-specific manner. Our analyses demonstrate that EGRET networks can capture the impact of genetic variants on complex phenotypes, supporting a novel fine-scale stratification of individuals based on their genetic background. EGRET is available through the Network Zoo R package (netZooR v0.9; netzoo.github.io).


2021 ◽  
Author(s):  
Basavaraj Mallikarjunayya Vastrad ◽  
Chanabasayya Mallikarjunayya Vastrad

Non alcoholic fatty liver disease (NAFLD) is the most common metabolic disease in humans, affecting the majority of individuals. In the current investigation, we aim to identify potential key genes linked with NAFLD through bioinformatics analyses of next generation sequencing (NGS) dataset. NGS dataset of GSE135251 from the Gene Expression Omnibus (GEO) database were retrieved. Differentially expressed genes (DEGs) were obtained by DESeq2 package. g:Profiler database was further used to identify the potential gene ontology (GO) and REACTOME pathways. Protein-protein interaction (PPI) network was constructed using the Hippie interactome database. miRNet and NetworkAnalyst databases were used to establish a miRNA-hub gene regulatory network and TF-hub gene regulatory network for the hub genes. Hub genes were verified based on receiver operating characteristic curve (ROC) analysis. Totally, 951 DEGs were identified including 476 up regulated genes and 475 down regulated genes screened in NAFLD and normal control. GO showed that DEGs were significantly enhanced for signaling and regulation of biological quality. REACTOME pathway analysis revealed that DEGs were enriched in signaling by interleukins and extracellular matrix organization. ESR2, JUN, PTN, PTGER3, CEBPB, IKBKG, HSPA8, SFN, CDKN1A and E2F1 were indicated as hub genes from PPI network, miRNA-hub gene regulatory network and TF-hub gene regulatory network. Furthermore, ROC analysis revealed that ESR2, JUN, PTN, PTGER3, CEBPB, IKBKG, HSPA8, SFN, CDKN1A and E2F1 might serve as diagnostic biomarkers in NAFLD. The current investigation provided insights into the molecular mechanism of NAFLD that might be useful in further investigations.


2020 ◽  
Author(s):  
Maud Fagny ◽  
Marieke Lydia Kuijjer ◽  
Maike Stam ◽  
Johann Joets ◽  
Olivier Turc ◽  
...  

AbstractEnhancers are important regulators of gene expression during numerous crucial processes including tissue differentiation across development. In plants, their recent molecular characterization revealed their capacity to activate the expression of several target genes through the binding of transcription factors. Nevertheless, identifying these target genes at a genome-wide level remains a challenge, in particular in species with large genomes, where enhancers and target genes can be hundreds of kilobases away. Therefore, the contribution of enhancers to regulatory network is still poorly understood in plants. In this study, we investigate the enhancer-driven regulatory network of two maize tissues at different stages: leaves at seedling stage and husks (bracts) at flowering. Using a systems biology approach, we integrate genomic, epigenomic and transcriptomic data to model the regulatory relationship between transcription factors and their potential target genes. We identify regulatory modules specific to husk and V2-IST, and show that they are involved in distinct functions related to the biology of each tissue. We evidence enhancers exhibiting binding sites for two distinct transcription factor families (DOF and AP2/ERF) that drive the tissue-specificity of gene expression in seedling immature leaf and husk. Analysis of the corresponding enhancer sequences reveals that two different transposable element families (TIR transposon Mutator and MITE Pif/Harbinger) have shaped the regulatory network in each tissue, and that MITEs have provided new transcription factor binding sites that are involved in husk tissue-specificity.SignificanceEnhancers play a major role in regulating tissue-specific gene expression in higher eukaryotes, including angiosperms. While molecular characterization of enhancers has improved over the past years, identifying their target genes at the genome-wide scale remains challenging. Here, we integrate genomic, epigenomic and transcriptomic data to decipher the tissue-specific gene regulatory network controlled by enhancers at two different stages of maize leaf development. Using a systems biology approach, we identify transcription factor families regulating gene tissue-specific expression in husk and seedling leaves, and characterize the enhancers likely to be involved. We show that a large part of maize enhancers is derived from transposable elements, which can provide novel transcription factor binding sites crucial to the regulation of tissue-specific biological functions.


Sign in / Sign up

Export Citation Format

Share Document