A Bayesian Approach for Learning Gene Networks Underlying Disease Severity in COPD.
Ontology highlight
ABSTRACT:
SUBMITTER: Shaddox E
PROVIDER: S-EPMC8078135 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:There is a lack of predictive markers for early and rapid identification of disease progression in COVID-19 patients. Our study aims at identifying microRNAs (miRNAs)/small nucleolar RNAs (snoRNAs) as potential biomarkers of COVID-19 severity. Using differential expression analysis of microarray data (n = 29), we identified hsa-miR-1246, ACA40, hsa-miR-4532, hsa-miR-145-5p, and ACA18 as the top five differentially expressed transcripts in severe versus asymptomatic, and ACA40, hsa-miR-3609, ENSG00000212378 (SNORD78), hsa-miR-1231, hsa-miR-885-3p as the most significant five in severe versus mild cases. Moreover, we found that white blood cell (WBC) count, absolute neutrophil count (ANC), neutrophil (%), lymphocyte (%), red blood cell (RBC) count, hemoglobin, hematocrit, D-Dimer, and albumin are significantly correlated with the identified differentially expressed miRNAs and snoRNAs. We report a unique miRNA and snoRNA profile that is associated with a higher risk of severity in a cohort of SARS-CoV-2 infected patients. Altogether, we present a differential expression analysis of COVID-19-associated microRNA (miRNA)/small nucleolar RNA (snoRNA) signature, highlighting their importance in SARS-CoV-2 infection.
| S-EPMC8538251 | biostudies-literature
Project description:MotivationOscillations lie at the core of many biological processes, from the cell cycle, to circadian oscillations and developmental processes. Time-keeping mechanisms are essential to enable organisms to adapt to varying conditions in environmental cycles, from day/night to seasonal. Transcriptional regulatory networks are one of the mechanisms behind these biological oscillations. However, while identifying cyclically expressed genes from time series measurements is relatively easy, determining the structure of the interaction network underpinning the oscillation is a far more challenging problem.ResultsHere, we explicitly leverage the oscillatory nature of the transcriptional signals and present a method for reconstructing network interactions tailored to this special but important class of genetic circuits. Our method is based on projecting the signal onto a set of oscillatory basis functions using a Discrete Fourier Transform. We build a Bayesian Hierarchical model within a frequency domain linear model in order to enforce sparsity and incorporate prior knowledge about the network structure. Experiments on real and simulated data show that the method can lead to substantial improvements over competing approaches if the oscillatory assumption is met, and remains competitive also in cases it is not.AvailabilityDSS, experiment scripts and data are available at http://homepages.inf.ed.ac.uk/gsanguin/DSS.zip.Contactd.trejo-banos@sms.ed.ac.ukSupplementary informationSupplementary data are available at Bioinformatics online.
| S-EPMC4817140 | biostudies-literature
Project description:Availability of genome sequence, molecular, and clinical phenotype data for large patient cohorts generated by recent technological advances provides an opportunity to dissect the genetic architecture of complex diseases at system level. However, previous analyses of such data have largely focused on the co-localization of SNPs associated with clinical and expression traits, each identified from genome-wide association studies and expression quantitative trait locus mapping. Thus, their description of the molecular mechanisms behind the SNPs influencing clinical phenotypes was limited to the single gene linked to the co-localized SNP. Here we introduce PerturbNet, a statistical framework for learning gene networks that modulate the influence of genetic variants on phenotypes, using genetic variants as naturally occurring perturbation of a biological system. PerturbNet uses a probabilistic graphical model to directly model the cascade of perturbation from genetic variants to the gene network to the phenotype network along with the networks at each layer of the biological system. PerturbNet learns the entire model by solving a single optimization problem with an efficient algorithm that can analyze human genome-wide data within a few hours. PerturbNet inference procedures extract a detailed description of how the gene network modulates the genetic effects on phenotypes. Using simulated and asthma data, we demonstrate that PerturbNet improves statistical power for detecting disease-linked SNPs and identifies gene networks and network modules mediating the SNP effects on traits, providing deeper insights into the underlying molecular mechanisms.
| S-EPMC7584257 | biostudies-literature
Project description:We present the Fast Greedy Equivalence Search (FGES)-Merge, a new method for learning the structure of gene regulatory networks via merging locally learned Bayesian networks, based on the fast greedy equivalent search algorithm. The method is competitive with the state of the art in terms of the Matthews correlation coefficient, which takes into account both precision and recall, while also improving upon it in terms of speed, scaling up to tens of thousands of variables and being able to use empirical knowledge about the topological structure of gene regulatory networks. To showcase the ability of our method to scale to massive networks, we apply it to learning the gene regulatory network for the full human genome using data from samples of different brain structures (from the Allen Human Brain Atlas). Furthermore, this Bayesian network model should predict interactions between genes in a way that is clear to experts, following the current trends in explainable artificial intelligence. To achieve this, we also present a new open-access visualization tool that facilitates the exploration of massive networks and can aid in finding nodes of interest for experimental tests.
| S-EPMC10745139 | biostudies-literature
Project description:BACKGROUND:Reverse engineering of transcriptional regulatory networks (TRN) from genomics data has always represented a computational challenge in System Biology. The major issue is modeling the complex crosstalk among transcription factors (TFs) and their target genes, with a method able to handle both the high number of interacting variables and the noise in the available heterogeneous experimental sources of information. RESULTS:In this work, we propose a data fusion approach that exploits the integration of complementary omics-data as prior knowledge within a Bayesian framework, in order to learn and model large-scale transcriptional networks. We develop a hybrid structure-learning algorithm able to jointly combine TFs ChIP-Sequencing data and gene expression compendia to reconstruct TRNs in a genome-wide perspective. Applying our method to high-throughput data, we verified its ability to deal with the complexity of a genomic TRN, providing a snapshot of the synergistic TFs regulatory activity. Given the noisy nature of data-driven prior knowledge, which potentially contains incorrect information, we also tested the method's robustness to false priors on a benchmark dataset, comparing the proposed approach to other regulatory network reconstruction algorithms. We demonstrated the effectiveness of our framework by evaluating structural commonalities of our learned genomic network with other existing networks inferred by different DNA binding information-based methods. CONCLUSIONS:This Bayesian omics-data fusion based methodology allows to gain a genome-wide picture of the transcriptional interplay, helping to unravel key hierarchical transcriptional interactions, which could be subsequently investigated, and it represents a promising learning approach suitable for multi-layered genomic data integration, given its robustness to noisy sources and its tailored framework for handling high dimensional data.
| S-EPMC7257163 | biostudies-literature
Project description:Bayesian networks are probabilistic models that represent complex distributions in a modular way and have become very popular in many fields. There are many methods to build Bayesian networks from a random sample of independent and identically distributed observations. However, many observational studies are designed using some form of clustered sampling that introduces correlations between observations within the same cluster and ignoring this correlation typically inflates the rate of false positive associations. We describe a novel parameterization of Bayesian networks that uses random effects to model the correlation within sample units and can be used for structure and parameter learning from correlated data without inflating the Type I error rate. We compare different learning metrics using simulations and illustrate the method in two real examples: an analysis of genetic and non-genetic factors associated with human longevity from a family-based study, and an example of risk factors for complications of sickle cell anemia from a longitudinal study with repeated measures.
| S-EPMC4857179 | biostudies-literature
Project description:Transcriptional changes occur presymptomatically and throughout Huntington's disease (HD), motivating the study of transcriptional regulatory networks (TRNs) in HD We reconstructed a genome-scale model for the target genes of 718 transcription factors (TFs) in the mouse striatum by integrating a model of genomic binding sites with transcriptome profiling of striatal tissue from HD mouse models. We identified 48 differentially expressed TF-target gene modules associated with age- and CAG repeat length-dependent gene expression changes in Htt CAG knock-in mouse striatum and replicated many of these associations in independent transcriptomic and proteomic datasets. Thirteen of 48 of these predicted TF-target gene modules were also differentially expressed in striatal tissue from human disease. We experimentally validated a specific model prediction that SMAD3 regulates HD-related gene expression changes using chromatin immunoprecipitation and deep sequencing (ChIP-seq) of mouse striatum. We found CAG repeat length-dependent changes in the genomic occupancy of SMAD3 and confirmed our model's prediction that many SMAD3 target genes are downregulated early in HD.
| S-EPMC5868199 | biostudies-literature
Project description:Current algorithms for gene regulatory network construction based on Gaussian graphical models focuses on the deterministic decision of whether an edge exists. Both the probabilistic inference of edge existence and the relative strength of edges are often overlooked, either because the computational algorithms cannot account for this uncertainty or because it is not straightforward in implementation. In this study, we combine the Bayesian Markov random field and the conditional autoregressive (CAR) model to tackle simultaneously these two tasks. The uncertainty of edge existence and the relative strength of edges can be measured and quantified based on a Bayesian model such as the CAR model and the spike-and-slab lasso prior. In addition, the strength of the edges can be utilized to prioritize the importance of the edges in a network graph. Simulations and a glioblastoma cancer study were carried out to assess the proposed model's performance and to compare it with existing methods when a binary decision is of interest. The proposed approach shows stable performance and may provide novel structures with biological insights.
| S-EPMC9684306 | biostudies-literature
Project description:Motivated by the study of the molecular mechanism underlying type 1 diabetes with gene expression data collected from both patients and healthy controls at multiple time points, we propose a hybrid Bayesian method for jointly estimating multiple dependent Gaussian graphical models with data observed under distinct conditions, which avoids inversion of high-dimensional covariance matrices and thus can be executed very fast. We prove the consistency of the proposed method under mild conditions. The numerical results indicate the superiority of the proposed method over existing ones in both estimation accuracy and computational efficiency. Extension of the proposed method to joint estimation of multiple mixed graphical models is straightforward.
| S-EPMC8035990 | biostudies-literature
Project description:The causative organism, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), exhibits a wide spectrum of clinical manifestations in disease-ridden patients. Differences in the severity of COVID-19 ranges from asymptomatic infections and mild cases to the severe form, leading to acute respiratory distress syndrome (ARDS) and multiorgan failure with poor survival. MiRNAs can regulate various cellular processes, including proliferation, apoptosis, and differentiation, by binding to the 3′UTR of target mRNAs inducing their degradation, thus serving a fundamental role in post-transcriptional repression. Alterations of miRNA levels in the blood have been described in multiple inflammatory and infectious diseases, including SARS-related coronaviruses. We used microarrays to delineate the miRNAs and snoRNAs signature in the peripheral blood of severe COVID-19 cases (n=9), as compared to mild (n=10) and asymptomatic (n=10) patients, and identified differentially expressed transcripts in severe versus asymptomatic, and others in severe versus mild COVID-19 cases. A cohort of 29 male age-matched patients were selected. All patients were previously diagnosed with COVID-19 using TaqPath COVID-19 Combo Kit (Thermo Fisher Scientific, Waltham, Massachusetts), or Cobas SARS-CoV-2 Test (Roche Diagnostics, Rotkreuz, Switzerland), with a CT value < 30. Additional criterion for selection was age between 35 and 75 years. Participants were grouped into severe, mild and asymptomatic. Classifying severe cases was based on requirement of high-flow oxygen support and ICU admission (n=9). Whereas mild patients were identified based on symptoms and positive radiographic findings with pulmonary involvement (n=10). Patients with no clinical presentation were labelled as asymptomatic cases (n=10).
2022-03-01 | E-MTAB-10970 | biostudies-arrayexpress