1
|
Labani M, Beheshti A, O’Brien TA. GENet: A Graph-Based Model Leveraging Histone Marks and Transcription Factors for Enhanced Gene Expression Prediction. Genes (Basel) 2024; 15:938. [PMID: 39062717 PMCID: PMC11275947 DOI: 10.3390/genes15070938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2024] [Revised: 07/16/2024] [Accepted: 07/17/2024] [Indexed: 07/28/2024] Open
Abstract
Understanding the regulatory mechanisms of gene expression is a crucial objective in genomics. Although the DNA sequence near the transcription start site (TSS) offers valuable insights, recent methods suggest that analyzing only the surrounding DNA may not suffice to accurately predict gene expression levels. We developed GENet (Gene Expression Network from Histone and Transcription Factor Integration), a novel approach that integrates essential regulatory signals from transcription factors and histone modifications into a graph-based model. GENet extends beyond simple DNA sequence analysis by incorporating additional layers of genetic control, which are vital for determining gene expression. Our method markedly enhances the prediction of mRNA levels compared to previous models that depend solely on DNA sequence data. The results underscore the significance of including comprehensive regulatory information in gene expression studies. GENet emerges as a promising tool for researchers, with potential applications extending from fundamental biological research to the development of medical therapies.
Collapse
Affiliation(s)
- Mahdieh Labani
- School of Computing, Macquarie University, Sydney 2109, Australia; (M.L.); (T.A.O.)
| | - Amin Beheshti
- School of Computing, Macquarie University, Sydney 2109, Australia; (M.L.); (T.A.O.)
| | - Tracey A. O’Brien
- School of Computing, Macquarie University, Sydney 2109, Australia; (M.L.); (T.A.O.)
- Cancer Institute NSW, Sydney 2065, Australia
- School of Clinical Medicine, Medicine & Health, University of New South Wales (UNSW), Sydney 2052, Australia
| |
Collapse
|
2
|
Singh R, Wu AP, Mudide A, Berger B. Causal gene regulatory analysis with RNA velocity reveals an interplay between slow and fast transcription factors. Cell Syst 2024; 15:462-474.e5. [PMID: 38754366 DOI: 10.1016/j.cels.2024.04.005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Revised: 08/25/2023] [Accepted: 04/18/2024] [Indexed: 05/18/2024]
Abstract
Single-cell expression dynamics, from differentiation trajectories or RNA velocity, have the potential to reveal causal links between transcription factors (TFs) and their target genes in gene regulatory networks (GRNs). However, existing methods either overlook these expression dynamics or necessitate that cells be ordered along a linear pseudotemporal axis, which is incompatible with branching trajectories. We introduce Velorama, an approach to causal GRN inference that represents single-cell differentiation dynamics as a directed acyclic graph of cells, constructed from pseudotime or RNA velocity measurements. Additionally, Velorama enables the estimation of the speed at which TFs influence target genes. Applying Velorama, we uncover evidence that the speed of a TF's interactions is tied to its regulatory function. For human corticogenesis, we find that slow TFs are linked to gliomas, while fast TFs are associated with neuropsychiatric diseases. We expect Velorama to become a critical part of the RNA velocity toolkit for investigating the causal drivers of differentiation and disease.
Collapse
Affiliation(s)
- Rohit Singh
- Computer Science and Artificial Intelligence Laboratory, MIT, Cambridge, MA 02139, USA.
| | - Alexander P Wu
- Computer Science and Artificial Intelligence Laboratory, MIT, Cambridge, MA 02139, USA
| | - Anish Mudide
- Phillips Exeter Academy, Exeter, NH 03883, USA; Computer Science and Artificial Intelligence Laboratory and Department of Mathematics, MIT, Cambridge, MA 02139, USA
| | - Bonnie Berger
- Computer Science and Artificial Intelligence Laboratory and Department of Mathematics, MIT, Cambridge, MA 02139, USA.
| |
Collapse
|
3
|
Swindell WR. Meta-analysis of differential gene expression in lower motor neurons isolated by laser capture microdissection from post-mortem ALS spinal cords. Front Genet 2024; 15:1385114. [PMID: 38689650 PMCID: PMC11059082 DOI: 10.3389/fgene.2024.1385114] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2024] [Accepted: 04/03/2024] [Indexed: 05/02/2024] Open
Abstract
Introduction ALS is a fatal neurodegenerative disease for which underlying mechanisms are incompletely understood. The motor neuron is a central player in ALS pathogenesis but different transcriptome signatures have been derived from bulk analysis of post-mortem tissue and iPSC-derived motor neurons (iPSC-MNs). Methods This study performed a meta-analysis of six gene expression studies (microarray and RNA-seq) in which laser capture microdissection (LCM) was used to isolate lower motor neurons from post-mortem spinal cords of ALS and control (CTL) subjects. Differentially expressed genes (DEGs) with consistent ALS versus CTL expression differences across studies were identified. Results The analysis identified 222 ALS-increased DEGs (FDR <0.10, SMD >0.80) and 278 ALS-decreased DEGs (FDR <0.10, SMD < -0.80). ALS-increased DEGs were linked to PI3K-AKT signaling, innate immunity, inflammation, motor neuron differentiation and extracellular matrix. ALS-decreased DEGs were associated with the ubiquitin-proteosome system, microtubules, axon growth, RNA-binding proteins and synaptic membrane. ALS-decreased DEG mRNAs frequently interacted with RNA-binding proteins (e.g., FUS, HuR). The complete set of DEGs (increased and decreased) overlapped significantly with genes near ALS-associated SNP loci (p < 0.01). Transcription factor target motifs with increased proximity to ALS-increased DEGs were identified, most notably DNA elements predicted to interact with forkhead transcription factors (e.g., FOXP1) and motor neuron and pancreas homeobox 1 (MNX1). Some of these DNA elements overlie ALS-associated SNPs within known enhancers and are predicted to have genotype-dependent MNX1 interactions. DEGs were compared to those identified from SOD1-G93A mice and bulk spinal cord segments or iPSC-MNs from ALS patients. There was good correspondence with transcriptome changes from SOD1-G93A mice (r ≤ 0.408) but most DEGs were not differentially expressed in bulk spinal cords or iPSC-MNs and transcriptome-wide effect size correlations were weak (bulk tissue: r ≤ 0.207, iPSC-MN: r ≤ 0.037). Conclusion This study defines a robust transcriptome signature from LCM-based motor neuron studies of post-mortem tissue from ALS and CTL subjects. This signature differs from those obtained from analysis of bulk spinal cord segments and iPSC-MNs. Results provide insight into mechanisms underlying gene dysregulation in ALS and highlight connections between these mechanisms, ALS genetics, and motor neuron biology.
Collapse
Affiliation(s)
- William R. Swindell
- Department of Internal Medicine, Division of Hospital Medicine, University of Texas Southwestern Medical Center, Dallas, TX, United States
| |
Collapse
|
4
|
Liao Z, Tang S, Nozawa K, Shimada K, Ikawa M, Monsivais D, Matzuk M. Affinity-tagged SMAD1 and SMAD5 mouse lines reveal transcriptional reprogramming mechanisms during early pregnancy. eLife 2024; 12:RP91434. [PMID: 38536963 PMCID: PMC10972565 DOI: 10.7554/elife.91434] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/04/2024] Open
Abstract
Endometrial decidualization, a prerequisite for successful pregnancies, relies on transcriptional reprogramming driven by progesterone receptor (PR) and bone morphogenetic protein (BMP)-SMAD1/SMAD5 signaling pathways. Despite their critical roles in early pregnancy, how these pathways intersect in reprogramming the endometrium into a receptive state remains unclear. To define how SMAD1 and/or SMAD5 integrate BMP signaling in the uterus during early pregnancy, we generated two novel transgenic mouse lines with affinity tags inserted into the endogenous SMAD1 and SMAD5 loci (Smad1HA/HA and Smad5PA/PA). By profiling the genome-wide distribution of SMAD1, SMAD5, and PR in the mouse uterus, we demonstrated the unique and shared roles of SMAD1 and SMAD5 during the window of implantation. We also showed the presence of a conserved SMAD1, SMAD5, and PR genomic binding signature in the uterus during early pregnancy. To functionally characterize the translational aspects of our findings, we demonstrated that SMAD1/5 knockdown in human endometrial stromal cells suppressed expressions of canonical decidual markers (IGFBP1, PRL, FOXO1) and PR-responsive genes (RORB, KLF15). Here, our studies provide novel tools to study BMP signaling pathways and highlight the fundamental roles of SMAD1/5 in mediating both BMP signaling pathways and the transcriptional response to progesterone (P4) during early pregnancy.
Collapse
Affiliation(s)
- Zian Liao
- Department of Pathology & Immunology, Baylor College of MedicineHoustonUnited States
- Graduate Program of Genetics and Genomics, Baylor College of MedicineHoustonUnited States
- Department of Molecular and Human Genetics, Baylor College of MedicineHoustonUnited States
- Center for Drug Discovery, Baylor College of MedicineHoustonUnited States
| | - Suni Tang
- Department of Pathology & Immunology, Baylor College of MedicineHoustonUnited States
- Department of Molecular and Human Genetics, Baylor College of MedicineHoustonUnited States
| | - Kaori Nozawa
- Department of Pathology & Immunology, Baylor College of MedicineHoustonUnited States
- Center for Drug Discovery, Baylor College of MedicineHoustonUnited States
| | - Keisuke Shimada
- Research Institute for Microbial Diseases, Osaka UniversityOsakaJapan
| | - Masahito Ikawa
- Research Institute for Microbial Diseases, Osaka UniversityOsakaJapan
| | - Diana Monsivais
- Department of Pathology & Immunology, Baylor College of MedicineHoustonUnited States
- Center for Drug Discovery, Baylor College of MedicineHoustonUnited States
| | - Martin Matzuk
- Department of Pathology & Immunology, Baylor College of MedicineHoustonUnited States
- Graduate Program of Genetics and Genomics, Baylor College of MedicineHoustonUnited States
- Department of Molecular and Human Genetics, Baylor College of MedicineHoustonUnited States
- Center for Drug Discovery, Baylor College of MedicineHoustonUnited States
| |
Collapse
|
5
|
Liao Z, Tang S, Nozawa K, Shimada K, Ikawa M, Monsivais D, Matzuk MM. Affinity-tagged SMAD1 and SMAD5 mouse lines reveal transcriptional reprogramming mechanisms during early pregnancy. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.09.25.559321. [PMID: 38106095 PMCID: PMC10723262 DOI: 10.1101/2023.09.25.559321] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Endometrial decidualization, a prerequisite for successful pregnancies, relies on transcriptional reprogramming driven by progesterone receptor (PR) and bone morphogenetic protein (BMP)-SMAD1/SMAD5 signaling pathways. Despite their critical roles in early pregnancy, how these pathways intersect in reprogramming the endometrium into a receptive state remains unclear. To define how SMAD1 and/or SMAD5 integrate BMP signaling in the uterus during early pregnancy, we generated two novel transgenic mouse lines with affinity tags inserted into the endogenous SMAD1 and SMAD5 loci (Smad1HA/HA and Smad5PA/PA). By profiling the genome-wide distribution of SMAD1, SMAD5, and PR in the mouse uterus, we demonstrated the unique and shared roles of SMAD1 and SMAD5 during the window of implantation. We also showed the presence of a conserved SMAD1, SMAD5, and PR genomic binding signature in the uterus during early pregnancy. To functionally characterize the translational aspects of our findings, we demonstrated that SMAD1/5 knockdown in human endometrial stromal cells suppressed expressions of canonical decidual markers (IGFBP1, PRL, FOXO1) and PR-responsive genes (RORB, KLF15). Here, our studies provide novel tools to study BMP signaling pathways and highlight the fundamental roles of SMAD1/5 in mediating both BMP signaling pathways and the transcriptional response to progesterone (P4) during early pregnancy.
Collapse
Affiliation(s)
- Zian Liao
- Department of Pathology & Immunology, Baylor College of Medicine, Houston, TX, 77030, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
- Graduate Program of Genetics and Genomics, Baylor College of Medicine, Houston, TX, 77030, USA
- Center for Drug Discovery, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Suni Tang
- Department of Pathology & Immunology, Baylor College of Medicine, Houston, TX, 77030, USA
- Center for Drug Discovery, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Kaori Nozawa
- Department of Pathology & Immunology, Baylor College of Medicine, Houston, TX, 77030, USA
- Center for Drug Discovery, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Keisuke Shimada
- Research Institute for Microbial Diseases, Osaka University, Osaka, 565-0871, Japan
| | - Masahito Ikawa
- Research Institute for Microbial Diseases, Osaka University, Osaka, 565-0871, Japan
| | - Diana Monsivais
- Department of Pathology & Immunology, Baylor College of Medicine, Houston, TX, 77030, USA
- Center for Drug Discovery, Baylor College of Medicine, Houston, TX, 77030, USA
| | - Martin M. Matzuk
- Department of Pathology & Immunology, Baylor College of Medicine, Houston, TX, 77030, USA
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX, 77030, USA
- Graduate Program of Genetics and Genomics, Baylor College of Medicine, Houston, TX, 77030, USA
- Center for Drug Discovery, Baylor College of Medicine, Houston, TX, 77030, USA
| |
Collapse
|
6
|
Tan W, Shen Y. Multimodal learning of noncoding variant effects using genome sequence and chromatin structure. Bioinformatics 2023; 39:btad541. [PMID: 37669132 PMCID: PMC10502240 DOI: 10.1093/bioinformatics/btad541] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Revised: 08/28/2023] [Accepted: 09/04/2023] [Indexed: 09/07/2023] Open
Abstract
MOTIVATION A growing amount of noncoding genetic variants, including single-nucleotide polymorphisms, are found to be associated with complex human traits and diseases. Their mechanistic interpretation is relatively limited and can use the help from computational prediction of their effects on epigenetic profiles. However, current models often focus on local, 1D genome sequence determinants and disregard global, 3D chromatin structure that critically affects epigenetic events. RESULTS We find that noncoding variants of unexpected high similarity in epigenetic profiles, with regards to their relatively low similarity in local sequences, can be largely attributed to their proximity in chromatin structure. Accordingly, we have developed a multimodal deep learning scheme that incorporates both data of 1D genome sequence and 3D chromatin structure for predicting noncoding variant effects. Specifically, we have integrated convolutional and recurrent neural networks for sequence embedding and graph neural networks for structure embedding despite the resolution gap between the two types of data, while utilizing recent DNA language models. Numerical results show that our models outperform competing sequence-only models in predicting epigenetic profiles and their use of long-range interactions complement sequence-only models in extracting regulatory motifs. They prove to be excellent predictors for noncoding variant effects in gene expression and pathogenicity, whether in unsupervised "zero-shot" learning or supervised "few-shot" learning. AVAILABILITY AND IMPLEMENTATION Codes and data can be accessed at https://github.com/Shen-Lab/ncVarPred-1D3D and https://zenodo.org/record/7975777.
Collapse
Affiliation(s)
- Wuwei Tan
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, TX 77843, United States
| | - Yang Shen
- Department of Electrical and Computer Engineering, Texas A&M University, College Station, TX 77843, United States
- Department of Computer Science and Engineering, Texas A&M University, College Station, TX 77843, United States
- Institute of Biosciences and Technology and Department of Translational Medical Sciences, College of Medicine, Texas A&M University, Houston, TX 77030, United States
| |
Collapse
|
7
|
Berres S, Gromoll J, Wöste M, Sandmann S, Laurentino S. OGRE: calculate, visualize, and analyze overlap between genomic input regions and public annotations. BMC Bioinformatics 2023; 24:300. [PMID: 37496002 PMCID: PMC10369718 DOI: 10.1186/s12859-023-05422-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2022] [Accepted: 07/18/2023] [Indexed: 07/28/2023] Open
Abstract
BACKGROUND Modern genome sequencing leads to an ever-growing collection of genomic annotations. Combining these elements with a set of input regions (e.g. genes) would yield new insights in genomic associations, such as those involved in gene regulation. The required data are scattered across different databases making a manual approach tiresome, unpractical, and prone to error. Semi-automatic approaches require programming skills in data parsing, processing, overlap calculation, and visualization, which most biomedical researchers lack. Our aim was to develop an automated tool providing all necessary algorithms, benefiting both bioinformaticians and researchers without bioinformatic training. RESULTS We developed overlapping annotated genomic regions (OGRE) as a comprehensive tool to associate and visualize input regions with genomic annotations. It does so by parsing regions of interest, mining publicly available annotations, and calculating possible overlaps between them. The user can thus identify location, type, and number of associated regulatory elements. Results are presented as easy to understand visualizations and result tables. We applied OGRE to recent studies and could show high reproducibility and potential new insights. To demonstrate OGRE's performance in terms of running time and output, we have conducted a benchmark and compared its features with similar tools. CONCLUSIONS OGRE's functions and built-in annotations can be applied as a downstream overlap association step, which is compatible with most genomic sequencing outputs, and can thus enrich pre-existing analyses pipelines. Compared to similar tools, OGRE shows competitive performance, offers additional features, and has been successfully applied to two recent studies. Overall, OGRE addresses the lack of tools for automatic analysis, local genomic overlap calculation, and visualization by providing an easy to use, end-to-end solution for both biologists and computational scientists.
Collapse
Affiliation(s)
- Sven Berres
- Centre of Reproductive Medicine and Andrology, University of Münster, Albert-Schweitzer-Campus 1 Building D11, 48149, Munster, Germany
| | - Jörg Gromoll
- Centre of Reproductive Medicine and Andrology, University of Münster, Albert-Schweitzer-Campus 1 Building D11, 48149, Munster, Germany
| | - Marius Wöste
- Institute of Medical Informatics, University of Münster, Albert-Schweitzer-Campus 1 Building A11, 48149, Munster, Germany
| | - Sarah Sandmann
- Institute of Medical Informatics, University of Münster, Albert-Schweitzer-Campus 1 Building A11, 48149, Munster, Germany
| | - Sandra Laurentino
- Centre of Reproductive Medicine and Andrology, University of Münster, Albert-Schweitzer-Campus 1 Building D11, 48149, Munster, Germany.
| |
Collapse
|
8
|
Theodoris CV, Xiao L, Chopra A, Chaffin MD, Al Sayed ZR, Hill MC, Mantineo H, Brydon EM, Zeng Z, Liu XS, Ellinor PT. Transfer learning enables predictions in network biology. Nature 2023; 618:616-624. [PMID: 37258680 PMCID: PMC10949956 DOI: 10.1038/s41586-023-06139-9] [Citation(s) in RCA: 78] [Impact Index Per Article: 78.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Accepted: 04/27/2023] [Indexed: 06/02/2023]
Abstract
Mapping gene networks requires large amounts of transcriptomic data to learn the connections between genes, which impedes discoveries in settings with limited data, including rare diseases and diseases affecting clinically inaccessible tissues. Recently, transfer learning has revolutionized fields such as natural language understanding1,2 and computer vision3 by leveraging deep learning models pretrained on large-scale general datasets that can then be fine-tuned towards a vast array of downstream tasks with limited task-specific data. Here, we developed a context-aware, attention-based deep learning model, Geneformer, pretrained on a large-scale corpus of about 30 million single-cell transcriptomes to enable context-specific predictions in settings with limited data in network biology. During pretraining, Geneformer gained a fundamental understanding of network dynamics, encoding network hierarchy in the attention weights of the model in a completely self-supervised manner. Fine-tuning towards a diverse panel of downstream tasks relevant to chromatin and network dynamics using limited task-specific data demonstrated that Geneformer consistently boosted predictive accuracy. Applied to disease modelling with limited patient data, Geneformer identified candidate therapeutic targets for cardiomyopathy. Overall, Geneformer represents a pretrained deep learning model from which fine-tuning towards a broad range of downstream applications can be pursued to accelerate discovery of key network regulators and candidate therapeutic targets.
Collapse
Affiliation(s)
- Christina V Theodoris
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA.
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Division of Genetics and Genomics, Boston Children's Hospital, Boston, MA, USA.
- Harvard Medical School Genetics Training Program, Boston, USA.
| | - Ling Xiao
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, USA
| | - Anant Chopra
- Precision Cardiology Laboratory, Bayer US LLC, Cambridge, MA, USA
| | - Mark D Chaffin
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Zeina R Al Sayed
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Matthew C Hill
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, USA
| | - Helene Mantineo
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, USA
| | | | - Zexian Zeng
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
| | - X Shirley Liu
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Patrick T Ellinor
- Cardiovascular Disease Initiative and Precision Cardiology Laboratory, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, USA.
| |
Collapse
|
9
|
Morin A, Chu ECP, Sharma A, Adrian-Hamazaki A, Pavlidis P. Characterizing the targets of transcription regulators by aggregating ChIP-seq and perturbation expression data sets. Genome Res 2023; 33:763-778. [PMID: 37308292 PMCID: PMC10317128 DOI: 10.1101/gr.277273.122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2022] [Accepted: 04/26/2023] [Indexed: 06/14/2023]
Abstract
Mapping the gene targets of chromatin-associated transcription regulators (TRs) is a major goal of genomics research. ChIP-seq of TRs and experiments that perturb a TR and measure the differential abundance of gene transcripts are a primary means by which direct relationships are tested on a genomic scale. It has been reported that there is a poor overlap in the evidence across gene regulation strategies, emphasizing the need for integrating results from multiple experiments. Although research consortia interested in gene regulation have produced a valuable trove of high-quality data, there is an even greater volume of TR-specific data throughout the literature. In this study, we show a workflow for the identification, uniform processing, and aggregation of ChIP-seq and TR perturbation experiments for the ultimate purpose of ranking human and mouse TR-target interactions. Focusing on an initial set of eight regulators (ASCL1, HES1, MECP2, MEF2C, NEUROD1, PAX6, RUNX1, and TCF4), we identified 497 experiments suitable for analysis. We used this corpus to examine data concordance, to identify systematic patterns of the two data types, and to identify putative orthologous interactions between human and mouse. We build upon commonly used strategies to forward a procedure for aggregating and combining these two genomic methodologies, assessing these rankings against independent literature-curated evidence. Beyond a framework extensible to other TRs, our work also provides empirically ranked TR-target listings, as well as transparent experiment-level gene summaries for community use.
Collapse
Affiliation(s)
- Alexander Morin
- Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
- Department of Psychiatry, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
- Graduate Program in Bioinformatics, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - Eric Ching-Pan Chu
- Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
- Department of Psychiatry, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
- Graduate Program in Bioinformatics, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - Aman Sharma
- Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - Alex Adrian-Hamazaki
- Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
- Department of Psychiatry, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
- Graduate Program in Bioinformatics, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - Paul Pavlidis
- Michael Smith Laboratories, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada;
- Department of Psychiatry, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| |
Collapse
|
10
|
Lynch AW, Theodoris CV, Long HW, Brown M, Liu XS, Meyer CA. MIRA: joint regulatory modeling of multimodal expression and chromatin accessibility in single cells. Nat Methods 2022; 19:1097-1108. [PMID: 36068320 PMCID: PMC9517733 DOI: 10.1038/s41592-022-01595-z] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 07/26/2022] [Indexed: 02/06/2023]
Abstract
Rigorously comparing gene expression and chromatin accessibility in the same single cells could illuminate the logic of how coupling or decoupling of these mechanisms regulates fate commitment. Here we present MIRA, probabilistic multimodal models for integrated regulatory analysis, a comprehensive methodology that systematically contrasts transcription and accessibility to infer the regulatory circuitry driving cells along cell state trajectories. MIRA leverages topic modeling of cell states and regulatory potential modeling of individual gene loci. MIRA thereby represents cell states in an efficient and interpretable latent space, infers high-fidelity cell state trees, determines key regulators of fate decisions at branch points and exposes the variable influence of local accessibility on transcription at distinct loci. Applied to epidermal differentiation and embryonic brain development from two different multimodal platforms, MIRA revealed that early developmental genes were tightly regulated by local chromatin landscape whereas terminal fate genes were titrated without requiring extensive chromatin remodeling.
Collapse
Affiliation(s)
- Allen W Lynch
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA.,Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, USA
| | - Christina V Theodoris
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA.,Division of Genetics and Genomics, Boston Children's Hospital, Boston, MA, USA.,Harvard Medical School Genetics Training Program, Boston, MA, USA
| | - Henry W Long
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, USA.,Department of Medical Oncology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, and Harvard Medical School, Boston, MA, USA
| | - Myles Brown
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, USA.,Department of Medical Oncology, Dana-Farber Cancer Institute, Brigham and Women's Hospital, and Harvard Medical School, Boston, MA, USA
| | - X Shirley Liu
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA. .,Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, USA. .,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA.
| | - Clifford A Meyer
- Department of Data Science, Dana-Farber Cancer Institute, Boston, MA, USA. .,Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, USA. .,Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, MA, USA.
| |
Collapse
|
11
|
Parisutham V, Chhabra S, Ali MZ, Brewster RC. Tunable transcription factor library for robust quantification of regulatory properties in Escherichia coli. Mol Syst Biol 2022; 18:e10843. [PMID: 35694815 PMCID: PMC9189660 DOI: 10.15252/msb.202110843] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2021] [Revised: 05/11/2022] [Accepted: 05/13/2022] [Indexed: 11/12/2022] Open
Abstract
Predicting the quantitative regulatory function of transcription factors (TFs) based on factors such as binding sequence, binding location, and promoter type is not possible. The interconnected nature of gene networks and the difficulty in tuning individual TF concentrations make the isolated study of TF function challenging. Here, we present a library of Escherichia coli strains designed to allow for precise control of the concentration of individual TFs enabling the study of the role of TF concentration on physiology and regulation. We demonstrate the usefulness of this resource by measuring the regulatory function of the zinc‐responsive TF, ZntR, and the paralogous TF pair, GalR/GalS. For ZntR, we find that zinc alters ZntR regulatory function in a way that enables activation of the regulated gene to be robust with respect to ZntR concentration. For GalR and GalS, we are able to demonstrate that these paralogous TFs have fundamentally distinct regulatory roles beyond differences in binding affinity.
Collapse
Affiliation(s)
- Vinuselvi Parisutham
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA, USA
| | - Shivani Chhabra
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Md Zulfikar Ali
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA, USA
| | - Robert C Brewster
- Department of Systems Biology, University of Massachusetts Chan Medical School, Worcester, MA, USA.,Department of Microbiology and Physiological Systems, University of Massachusetts Chan Medical School, Worcester, MA, USA
| |
Collapse
|
12
|
Domb K, Wang N, Hummel G, Liu C. Spatial Features and Functional Implications of Plant 3D Genome Organization. ANNUAL REVIEW OF PLANT BIOLOGY 2022; 73:173-200. [PMID: 35130445 DOI: 10.1146/annurev-arplant-102720-022810] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
The advent of high-throughput sequencing-based methods for chromatin conformation, accessibility, and immunoprecipitation assays has been a turning point in 3D genomics. Altogether, these new tools have been pushing upward the interpretation of pioneer cytogenetic evidence for a higher order in chromatin packing. Here, we review the latest development in our understanding of plant spatial genome structures and different levels of organization and discuss their functional implications. Then, we spotlight the complexity of organellar (i.e., mitochondria and plastids) genomes and discuss their 3D packing into nucleoids. Finally, we propose unaddressed research axes to investigate functional links between chromatin-like dynamics and transcriptional regulation within organellar nucleoids.
Collapse
Affiliation(s)
- Katherine Domb
- Institute of Biology, University of Hohenheim, Stuttgart, Germany;
| | - Nan Wang
- Institute of Biology, University of Hohenheim, Stuttgart, Germany;
| | - Guillaume Hummel
- Institute of Biology, University of Hohenheim, Stuttgart, Germany;
| | - Chang Liu
- Institute of Biology, University of Hohenheim, Stuttgart, Germany;
| |
Collapse
|
13
|
Human Cytomegalovirus IE2 Both Activates and Represses Initiation and Modulates Elongation in a Context-Dependent Manner. mBio 2022; 13:e0033722. [PMID: 35579393 PMCID: PMC9239164 DOI: 10.1128/mbio.00337-22] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Human cytomegalovirus (HCMV) immediate-early 2 (IE2) protein is a multifunctional transcription factor that is essential for lytic HCMV infection. IE2 functions as an activator of viral early genes, negatively regulates its own promoter, and is required for viral replication. The mechanisms by which IE2 executes these distinct functions are incompletely understood. Using PRO-Seq, which profiles nascent transcripts, and a recently developed DFF-chromatin immunoprecipitation (DFF-ChIP; employs chromatin digestion by the endonuclease DNA fragmentation factor prior to IP) approach that resolves occupancy and local chromatin environment, we show that IE2 controls viral gene transcription in three distinct capacities during late HCMV infection and reveal mechanisms that involve direct binding of IE2 to viral DNA. IE2 represses a subset of viral promoters by binding within their core promoter regions and blocking the assembly of preinitiation complexes (PICs). Remarkably, IE2 forms a repressive complex at the major immediate-early promoter region involving direct association of IE2 with nucleosomes and TBP. IE2 stimulates transcription by binding nearby, but not within, core promoter regions. In addition, IE2 functions as a direct roadblock to transcription elongation. At one locus, this function of IE2 appears to be important for the synthesis of a spliced viral RNA. Consistent with the minimal observed effects of IE2 depletion on host gene transcription, IE2 does not functionally engage the host genome. Our results reveal mechanisms of transcriptional control by IE2, uncover a previously unknown function of IE2 as a Pol II elongation modulator, and demonstrate that DFF-ChIP is a useful tool for probing transcription factor occupancy and interactions between transcription factors and nucleosomes at high resolution.
Collapse
|
14
|
Lu Y, Feng Z, Zhang S, Wang Y. Annotating regulatory elements by heterogeneous network embedding. Bioinformatics 2022; 38:2899-2911. [PMID: 35561169 PMCID: PMC9326849 DOI: 10.1093/bioinformatics/btac185] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2021] [Revised: 03/05/2022] [Accepted: 03/24/2022] [Indexed: 11/13/2022] Open
Abstract
MOTIVATION Regulatory elements (REs), such as enhancers and promoters, are known as regulatory sequences functional in a heterogeneous regulatory network to control gene expression by recruiting transcription regulators and carrying genetic variants in a context specific way. Annotating those REs relies on costly and labor-intensive next-generation sequencing and RNA-guided editing technologies in many cellular contexts. RESULTS We propose a systematic Gene Ontology Annotation method for Regulatory Elements (RE-GOA) by leveraging the powerful word embedding in natural language processing. We first assemble a heterogeneous network by integrating context specific regulations, protein-protein interactions and gene ontology (GO) terms. Then we perform network embedding and associate regulatory elements with GO terms by assessing their similarity in a low dimensional vector space. With three applications, we show that RE-GOA outperforms existing methods in annotating TFs' binding sites from ChIP-seq data, in functional enrichment analysis of differentially accessible peaks from ATAC-seq data, and in revealing genetic correlation among phenotypes from their GWAS summary statistics data. AVAILABILITY AND IMPLEMENTATION The source code and the systematic RE annotation for human and mouse are available at https://github.com/AMSSwanglab/RE-GOA. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Yurun Lu
- CEMS, NCMIS, HCMS, MADIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
- School of Mathematics, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Beijing 100049, China
| | - Zhanying Feng
- CEMS, NCMIS, HCMS, MADIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
- School of Mathematics, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Beijing 100049, China
| | - Songmao Zhang
- CEMS, NCMIS, HCMS, MADIS, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, Beijing 100190, China
| | - Yong Wang
- To whom correspondence should be addressed. or
| |
Collapse
|
15
|
Abrusán G, Ascher DB, Inouye M. Known allosteric proteins have central roles in genetic disease. PLoS Comput Biol 2022; 18:e1009806. [PMID: 35139069 PMCID: PMC10138267 DOI: 10.1371/journal.pcbi.1009806] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Revised: 04/27/2023] [Accepted: 01/05/2022] [Indexed: 12/15/2022] Open
Abstract
Allostery is a form of protein regulation, where ligands that bind sites located apart from the active site can modify the activity of the protein. The molecular mechanisms of allostery have been extensively studied, because allosteric sites are less conserved than active sites, and drugs targeting them are more specific than drugs binding the active sites. Here we quantify the importance of allostery in genetic disease. We show that 1) known allosteric proteins are central in disease networks, contribute to genetic disease and comorbidities much more than non-allosteric proteins, and there is an association between being allosteric and involvement in disease; 2) they are enriched in many major disease types like hematopoietic diseases, cardiovascular diseases, cancers, diabetes, or diseases of the central nervous system; 3) variants from cancer genome-wide association studies are enriched near allosteric proteins, indicating their importance to polygenic traits; and 4) the importance of allosteric proteins in disease is due, at least partly, to their central positions in protein-protein interaction networks, and less due to their dynamical properties.
Collapse
Affiliation(s)
- György Abrusán
- Cambridge Baker Systems Genomics Initiative, Department of Public Health and Primary Care, School of Medicine, University of Cambridge, Cambridge, United Kingdom
- * E-mail:
| | - David B. Ascher
- Department of Biochemistry, University of Cambridge, Cambridge, United Kingdom
- Structural Biology and Bioinformatics, Department of Biochemistry, Bio21 Institute, University of Melbourne, Melbourne, Australia
- Computational Biology and Clinical Informatics, Baker Heart and Diabetes Institute, Melbourne, Australia
| | - Michael Inouye
- Cambridge Baker Systems Genomics Initiative, Department of Public Health and Primary Care, School of Medicine, University of Cambridge, Cambridge, United Kingdom
- Cambridge Baker Systems Genomics Initiative, Baker Heart and Diabetes Institute, Melbourne, Australia
- British Heart Foundation Cardiovascular Epidemiology Unit, Department of Public Health and Primary Care, University of Cambridge, Cambridge, United Kingdom
- British Heart Foundation Centre of Research Excellence, University of Cambridge, Cambridge, United Kingdom
- Health Data Research UK Cambridge, Wellcome Genome Campus and University of Cambridge, Cambridge, United Kingdom
- The Alan Turing Institute, London, United Kingdom
| |
Collapse
|
16
|
Rundquist O, Nestor CE, Jenmalm MC, Hellberg S, Gustafsson M. Progesterone Inhibits the Establishment of Activation-Associated Chromatin During T H1 Differentiation. Front Immunol 2022; 13:835625. [PMID: 35185927 PMCID: PMC8848251 DOI: 10.3389/fimmu.2022.835625] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Accepted: 01/14/2022] [Indexed: 01/08/2023] Open
Abstract
TH1-mediated diseases such as multiple sclerosis (MS) and rheumatoid arthritis (RA) improve during pregnancy, coinciding with increasing levels of the pregnancy hormone progesterone (P4), highlighting P4 as a potential mediator of this immunomodulation. Here, we performed detailed characterization of how P4 affects the chromatin and transcriptomic landscape during early human TH1 differentiation, utilizing both ATAC-seq and RNA-seq. Time series analysis of the earlier events (0.5-24 hrs) during TH1 differentiation revealed that P4 counteracted many of the changes induced during normal differentiation, mainly by downregulating key regulatory genes and their upstream transcription factors (TFs) involved in the initial T-cell activation. Members of the AP-1 complex such as FOSL1, FOSL2, JUN and JUNB were particularly affected, in both in promoters and in distal regulatory elements. Moreover, the changes induced by P4 were significantly enriched for disease-associated changes related to both MS and RA, revealing several shared upstream TFs, where again JUN was highlighted to be of central importance. Our findings support an immune regulatory role for P4 during pregnancy by impeding T-cell activation, a crucial checkpoint during pregnancy and in T-cell mediated diseases, and a central event prior to T-cell lineage commitment. Indeed, P4 is emerging as a likely candidate involved in disease modulation during pregnancy and further studies evaluating P4 as a potential treatment option are needed.
Collapse
Affiliation(s)
- Olof Rundquist
- Bioinformatics, Department of Physics, Chemistry and Biology, Linköping University, Linköping, Sweden
| | - Colm E. Nestor
- Crown Princess Victoria Children’s Hospital, and Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Maria C. Jenmalm
- Division of Inflammation and Infection, Department of Biomedical and Clinical Sciences, Linköping University, Linköping, Sweden
| | - Sandra Hellberg
- Bioinformatics, Department of Physics, Chemistry and Biology, Linköping University, Linköping, Sweden
| | - Mika Gustafsson
- Bioinformatics, Department of Physics, Chemistry and Biology, Linköping University, Linköping, Sweden
| |
Collapse
|
17
|
Waymack R, Gad M, Wunderlich Z. Molecular competition can shape enhancer activity in the Drosophila embryo. iScience 2021; 24:103034. [PMID: 34568782 PMCID: PMC8449247 DOI: 10.1016/j.isci.2021.103034] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2021] [Revised: 07/27/2021] [Accepted: 08/20/2021] [Indexed: 01/12/2023] Open
Abstract
Transgenic reporters allow the measurement of regulatory DNA activity in vivo and consequently have long been useful tools for studying enhancers. Despite their utility, few studies have investigated the effects these reporters may have on the expression of other genes. Understanding these effects is required to accurately interpret reporter data and characterize gene regulatory mechanisms. By measuring the expression of Kruppel (Kr) enhancer reporters in live Drosophila embryos, we find reporters inhibit one another's expression and that of a nearby endogenous gene. Using synthetic transcription factor (TF) binding site arrays, we present evidence that competition for TFs is partially responsible for the observed transcriptional inhibition. We develop a simple thermodynamic model that predicts competition of the measured magnitude specifically when TF binding is restricted to distinct nuclear subregions. Our findings underline an unexpected role of the non-homogenous nature of the nucleus in regulating gene expression.
Collapse
Affiliation(s)
- Rachel Waymack
- Department of Developmental and Cell Biology, University of California, Irvine, CA 92697, USA
| | - Mario Gad
- Department of Developmental and Cell Biology, University of California, Irvine, CA 92697, USA
| | - Zeba Wunderlich
- Department of Developmental and Cell Biology, University of California, Irvine, CA 92697, USA
- Department of Biology, Boston University, 610 Commonwealth Ave., Boston, MA 02215, USA
- Biological Design Center, Boston University, 610 Commonwealth Avenue, Boston, MA 02215, USA
| |
Collapse
|
18
|
Wen Y, Dungan CM, Mobley CB, Valentino T, von Walden F, Murach KA. Nucleus Type-Specific DNA Methylomics Reveals Epigenetic "Memory" of Prior Adaptation in Skeletal Muscle. FUNCTION (OXFORD, ENGLAND) 2021; 2:zqab038. [PMID: 34870208 PMCID: PMC8636928 DOI: 10.1093/function/zqab038] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 07/28/2021] [Accepted: 08/03/2021] [Indexed: 02/07/2023]
Abstract
Using a mouse model of conditional and inducible in vivo fluorescent myonuclear labeling (HSA-GFP), sorting purification of nuclei, low-input reduced representation bisulfite sequencing (RRBS), and a translatable and reversible model of exercise (progressive weighted wheel running, PoWeR), we provide the first nucleus type-specific epigenetic information on skeletal muscle adaptation and detraining. Adult (>4 mo) HSA-GFP mice performed PoWeR for 8 wk then detrained for 12 wk; age-matched untrained mice were used to control for the long duration of the study. Myonuclei and interstitial nuclei from plantaris muscles were isolated for RRBS. Relative to untrained, PoWeR caused similar myonuclear CpG hypo- and hyper-methylation of promoter regions and substantial hypomethylation in interstitial nuclear promoters. Over-representation analysis of promoters revealed a larger number of hyper- versus hypo-methylated pathways in both nuclear populations after training and evidence for reciprocal regulation of methylation between nucleus types, with hypomethylation of promoter regions in Wnt signaling-related genes in myonuclei and hypermethylation in interstitial nuclei. After 12 wk of detraining, promoter CpGs in documented muscle remodeling-associated genes and pathways that were differentially methylated immediately after PoWeR were persistently differentially methylated in myonuclei, along with long-term promoter hypomethylation in interstitial nuclei. No enduring gene expression changes in muscle tissue were observed using RNA-sequencing. Upon 4 wk of retraining, mice that trained previously grew more at the whole muscle and fiber type-specific cellular level than training naïve mice, with no difference in myonuclear number. Muscle nuclei have a methylation epi-memory of prior training that may augment muscle adaptability to retraining.
Collapse
Affiliation(s)
- Yuan Wen
- Department of Physiology, University of Kentucky, Lexington, KY 40508, USA,The Center for Muscle Biology, University of Kentucky, Lexington, KY 40536, USA
| | - Cory M Dungan
- The Center for Muscle Biology, University of Kentucky, Lexington, KY 40536, USA,College of Health Sciences, University of Kentucky, Lexington, KY 40536, USA
| | - C Brooks Mobley
- Department of Physiology, University of Kentucky, Lexington, KY 40508, USA,The Center for Muscle Biology, University of Kentucky, Lexington, KY 40536, USA
| | - Taylor Valentino
- Department of Physiology, University of Kentucky, Lexington, KY 40508, USA,The Center for Muscle Biology, University of Kentucky, Lexington, KY 40536, USA
| | - Ferdinand von Walden
- Division of Pediatric Neurology, Department of Women's and Children's Health, Karolinska Institutet, Stockholm 171 77, Sweden
| | | |
Collapse
|
19
|
Cuesta-Astroz Y, Gischkow Rucatti G, Murgas L, SanMartín CD, Sanhueza M, Martin AJM. Filtering of Data-Driven Gene Regulatory Networks Using Drosophila melanogaster as a Case Study. Front Genet 2021; 12:649764. [PMID: 34394179 PMCID: PMC8355599 DOI: 10.3389/fgene.2021.649764] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Accepted: 04/30/2021] [Indexed: 01/12/2023] Open
Abstract
Gene Regulatory Networks (GRNs) allow the study of regulation of gene expression of whole genomes. Among the most relevant advantages of using networks to depict this key process, there is the visual representation of large amounts of information and the application of graph theory to generate new knowledge. Nonetheless, despite the many uses of GRNs, it is still difficult and expensive to assign Transcription Factors (TFs) to the regulation of specific genes. ChIP-Seq allows the determination of TF Binding Sites (TFBSs) over whole genomes, but it is still an expensive technique that can only be applied one TF at a time and requires replicates to reduce its noise. Once TFBSs are determined, the assignment of each TF and its binding sites to the regulation of specific genes is not trivial, and it is often performed by carrying out site-specific experiments that are unfeasible to perform in all possible binding sites. Here, we addressed these relevant issues with a two-step methodology using Drosophila melanogaster as a case study. First, our protocol starts by gathering all transcription factor binding sites (TFBSs) determined with ChIP-Seq experiments available at ENCODE and FlyBase. Then each TFBS is used to assign TFs to the regulation of likely target genes based on the TFBS proximity to the transcription start site of all genes. In the final step, to try to select the most likely regulatory TF from those previously assigned to each gene, we employ GENIE3, a random forest-based method, and more than 9,000 RNA-seq experiments from D. melanogaster. Following, we employed known TF protein-protein interactions to estimate the feasibility of regulatory events in our filtered networks. Finally, we show how known interactions between co-regulatory TFs of each gene increase after the second step of our approach, and thus, the consistency of the TF-gene assignment. Also, we employed our methodology to create a network centered on the Drosophila melanogaster gene Hr96 to demonstrate the role of this transcription factor on mitochondrial gene regulation.
Collapse
Affiliation(s)
- Yesid Cuesta-Astroz
- Colombian Institute of Tropical Medicine, CES University, Medellin, Colombia
| | | | - Leandro Murgas
- Laboratorio de Biologia de Redes, Centro de Genómica y Bioinformática, Facultad de Ciencias, Universidad Mayor, Santiago, Chile.,Programa de Doctorado en Genómica Integrativa, Vicerrectoría de Investigación, Universidad Mayor, Santiago, Chile
| | - Carol D SanMartín
- Departamento de Neurología y Neurocirugía, Hospital Clínico Universidad de Chile, Santiago, Chile.,Centro de Investigacíon Clínica Avanzada (CICA), Hospital Clínico Universidad de Chile, Santiago, Chile
| | - Mario Sanhueza
- Centro de Biología Integrativa, Facultad de Ciencias, Universidad Mayor, Santiago, Chile.,Escuela de Biotecnología, Facultad de Ciencias, Universidad Mayor, Santiago, Chile
| | - Alberto J M Martin
- Laboratorio de Biologia de Redes, Centro de Genómica y Bioinformática, Facultad de Ciencias, Universidad Mayor, Santiago, Chile.,Escuela de Biotecnología, Facultad de Ciencias, Universidad Mayor, Santiago, Chile
| |
Collapse
|
20
|
Puig RR, Boddie P, Khan A, Castro-Mondragon JA, Mathelier A. UniBind: maps of high-confidence direct TF-DNA interactions across nine species. BMC Genomics 2021; 22:482. [PMID: 34174819 PMCID: PMC8236138 DOI: 10.1186/s12864-021-07760-6] [Citation(s) in RCA: 37] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Accepted: 05/27/2021] [Indexed: 12/17/2022] Open
Abstract
BACKGROUND Transcription factors (TFs) bind specifically to TF binding sites (TFBSs) at cis-regulatory regions to control transcription. It is critical to locate these TF-DNA interactions to understand transcriptional regulation. Efforts to predict bona fide TFBSs benefit from the availability of experimental data mapping DNA binding regions of TFs (chromatin immunoprecipitation followed by sequencing - ChIP-seq). RESULTS In this study, we processed ~ 10,000 public ChIP-seq datasets from nine species to provide high-quality TFBS predictions. After quality control, it culminated with the prediction of ~ 56 million TFBSs with experimental and computational support for direct TF-DNA interactions for 644 TFs in > 1000 cell lines and tissues. These TFBSs were used to predict > 197,000 cis-regulatory modules representing clusters of binding events in the corresponding genomes. The high-quality of the TFBSs was reinforced by their evolutionary conservation, enrichment at active cis-regulatory regions, and capacity to predict combinatorial binding of TFs. Further, we confirmed that the cell type and tissue specificity of enhancer activity was correlated with the number of TFs with binding sites predicted in these regions. All the data is provided to the community through the UniBind database that can be accessed through its web-interface ( https://unibind.uio.no/ ), a dedicated RESTful API, and as genomic tracks. Finally, we provide an enrichment tool, available as a web-service and an R package, for users to find TFs with enriched TFBSs in a set of provided genomic regions. CONCLUSIONS UniBind is the first resource of its kind, providing the largest collection of high-confidence direct TF-DNA interactions in nine species.
Collapse
Affiliation(s)
- Rafael Riudavets Puig
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, 0349, Oslo, Norway
| | - Paul Boddie
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, 0349, Oslo, Norway
| | - Aziz Khan
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, 0349, Oslo, Norway
- Stanford Cancer Institute, Stanford University School of Medicine, Stanford, CA, 94305, USA
| | | | - Anthony Mathelier
- Centre for Molecular Medicine Norway (NCMM), Nordic EMBL Partnership, University of Oslo, 0349, Oslo, Norway.
- Department of Medical Genetics, Oslo University Hospital, Oslo, 0424, Norway.
| |
Collapse
|
21
|
Ong SYQ, Flyamer IM, Bickmore WA, Biddie SC. From bedside to bench: regulation of host factors in SARS-CoV-2 infection. Exp Mol Med 2021; 53:483-494. [PMID: 33828231 PMCID: PMC8024942 DOI: 10.1038/s12276-021-00595-x] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2020] [Revised: 02/01/2021] [Accepted: 02/15/2021] [Indexed: 01/08/2023] Open
Abstract
The zoonotic coronavirus SARS-CoV-2 (severe acute respiratory syndrome coronavirus-2), which causes COVID-19 (coronavirus disease-2019), has resulted in a pandemic. This has led to an urgent need to understand the molecular determinants of SARS-CoV-2 infection, factors associated with COVID-19 heterogeneity and severity, and therapeutic options for these patients. In this review, we discuss the role of host factors in SARS-CoV-2 infection and describe variations in host factor expression as mechanisms underlying the symptoms and severity of COVID-19. We focus on two host factors, angiotensin-converting enzyme 2 (ACE2) and transmembrane serine protease 2 (TMPRSS2), implicated in SARS-CoV-2 infection. We also discuss genetic variants associated with COVID-19 severity revealed in selected patients and based on genome-wide association studies (GWASs). Furthermore, we highlight important advances in cell and chromatin biology, such as single-cell RNA and chromatin sequencing and chromosomal conformation assays, as methods that may aid in the discovery of viral-host interactions in COVID-19. Understanding how regulation of host factor genes varies in physiological and pathological states might explain the heterogeneity observed in SARS-CoV-2 infection, help identify pathways for therapeutic development, and identify patients most likely to progress to severe COVID-19.
Collapse
Affiliation(s)
| | - Ilya M Flyamer
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, The University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Wendy A Bickmore
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, The University of Edinburgh, Edinburgh, EH4 2XU, UK
| | - Simon C Biddie
- Intensive Care Medicine, NHS Lothian, Edinburgh, UK.
- MRC Human Genetics Unit, Institute of Genetics and Molecular Medicine, The University of Edinburgh, Edinburgh, EH4 2XU, UK.
| |
Collapse
|
22
|
Zhou M, Li H, Wang X, Guan Y. Evidence of widespread, independent sequence signature for transcription factor cobinding. Genome Res 2021; 31:265-278. [PMID: 33303494 PMCID: PMC7849410 DOI: 10.1101/gr.267310.120] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Accepted: 12/03/2020] [Indexed: 01/03/2023]
Abstract
Transcription factors (TFs) are the vocabulary that genomes use to regulate gene expression and phenotypes. The interactions among TFs enrich this vocabulary and orchestrate diverse biological processes. Although simple models identify open chromatin and the presence of TF motifs as the two major contributors to TF binding patterns, it remains elusive what contributes to the in vivo TF cobinding landscape. In this study, we developed a machine learning algorithm to explore the contributors of the cobinding patterns. The algorithm substantially outperforms the state-of-the-field models for TF cobinding prediction. Game theory-based feature importance analysis reveals that, for most of the TF pairs we studied, independent motif sequences contribute one or more of the two TFs under investigation to their cobinding patterns. Such independent motif sequences include, but are not limited to, transcription initiation-related proteins and known TF complexes. We found the motif sequence signatures and the TFs are rarely mutual, corroborating a hierarchical and directional organization of the regulatory network and refuting the possibility of artifacts caused by shared sequence similarity with the TFs under investigation. We modeled such regulatory language with directed graphs, which reveal shared, global factors that are related to many binding and cobinding patterns.
Collapse
Affiliation(s)
- Manqi Zhou
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Hongyang Li
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Xueqing Wang
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Yuanfang Guan
- Department of Computational Medicine and Bioinformatics, University of Michigan, Ann Arbor, Michigan 48109, USA
| |
Collapse
|
23
|
Sammons MA, Nguyen TAT, McDade SS, Fischer M. Tumor suppressor p53: from engaging DNA to target gene regulation. Nucleic Acids Res 2020; 48:8848-8869. [PMID: 32797160 PMCID: PMC7498329 DOI: 10.1093/nar/gkaa666] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2020] [Revised: 07/24/2020] [Accepted: 07/30/2020] [Indexed: 12/13/2022] Open
Abstract
The p53 transcription factor confers its potent tumor suppressor functions primarily through the regulation of a large network of target genes. The recent explosion of next generation sequencing protocols has enabled the study of the p53 gene regulatory network (GRN) and underlying mechanisms at an unprecedented depth and scale, helping us to understand precisely how p53 controls gene regulation. Here, we discuss our current understanding of where and how p53 binds to DNA and chromatin, its pioneer-like role, and how this affects gene regulation. We provide an overview of the p53 GRN and the direct and indirect mechanisms through which p53 affects gene regulation. In particular, we focus on delineating the ubiquitous and cell type-specific network of regulatory elements that p53 engages; reviewing our understanding of how, where, and when p53 binds to DNA and the mechanisms through which these events regulate transcription. Finally, we discuss the evolution of the p53 GRN and how recent work has revealed remarkable differences between vertebrates, which are of particular importance to cancer researchers using mouse models.
Collapse
Affiliation(s)
- Morgan A Sammons
- Department of Biological Sciences and The RNA Institute, University at Albany, State University of New York, 1400 Washington Avenue, Albany, NY 12222, USA
| | - Thuy-Ai T Nguyen
- Genome Integrity & Structural Biology Laboratory and Immunity, Inflammation and Disease Laboratory, National Institute of Environmental Health Sciences/National Institutes of Health, 111 TW Alexander Drive, Research Triangle Park, NC 27709, USA
| | - Simon S McDade
- Patrick G Johnston Centre for Cancer Research, Queen's University Belfast, 97 Lisburn Road, Belfast BT9 7AE, UK
| | - Martin Fischer
- Computational Biology Group, Leibniz Institute on Aging – Fritz Lipmann Institute (FLI), Beutenbergstraße 11, 07745 Jena, Germany
| |
Collapse
|
24
|
Wang C, Sun D, Huang X, Wan C, Li Z, Han Y, Qin Q, Fan J, Qiu X, Xie Y, Meyer CA, Brown M, Tang M, Long H, Liu T, Liu XS. Integrative analyses of single-cell transcriptome and regulome using MAESTRO. Genome Biol 2020; 21:198. [PMID: 32767996 PMCID: PMC7412809 DOI: 10.1186/s13059-020-02116-x] [Citation(s) in RCA: 88] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2019] [Accepted: 07/23/2020] [Indexed: 12/15/2022] Open
Abstract
We present Model-based AnalysEs of Transcriptome and RegulOme (MAESTRO), a comprehensive open-source computational workflow ( http://github.com/liulab-dfci/MAESTRO ) for the integrative analyses of single-cell RNA-seq (scRNA-seq) and ATAC-seq (scATAC-seq) data from multiple platforms. MAESTRO provides functions for pre-processing, alignment, quality control, expression and chromatin accessibility quantification, clustering, differential analysis, and annotation. By modeling gene regulatory potential from chromatin accessibilities at the single-cell level, MAESTRO outperforms the existing methods for integrating the cell clusters between scRNA-seq and scATAC-seq. Furthermore, MAESTRO supports automatic cell-type annotation using predefined cell type marker genes and identifies driver regulators from differential scRNA-seq genes and scATAC-seq peaks.
Collapse
Affiliation(s)
- Chenfei Wang
- Department of Data Science, Dana-Farber Cancer Institute, Harvard T.H. Chan School of Public Health, Boston, MA, 02215, USA
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
| | - Dongqing Sun
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Science and Technology, Tongji University, Shanghai, 200433, China
| | - Xin Huang
- Beijing Institute of Radiation Medicine, Beijing, 100850, China
| | - Changxin Wan
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Science and Technology, Tongji University, Shanghai, 200433, China
| | - Ziyi Li
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Science and Technology, Tongji University, Shanghai, 200433, China
| | - Ya Han
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Science and Technology, Tongji University, Shanghai, 200433, China
| | - Qian Qin
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Science and Technology, Tongji University, Shanghai, 200433, China
| | - Jingyu Fan
- Clinical Translational Research Center, Shanghai Pulmonary Hospital, School of Life Science and Technology, Tongji University, Shanghai, 200433, China
| | - Xintao Qiu
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, 02215, USA
| | - Yingtian Xie
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, 02215, USA
| | - Clifford A Meyer
- Department of Data Science, Dana-Farber Cancer Institute, Harvard T.H. Chan School of Public Health, Boston, MA, 02215, USA
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
| | - Myles Brown
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, 02215, USA
| | - Ming Tang
- Department of Data Science, Dana-Farber Cancer Institute, Harvard T.H. Chan School of Public Health, Boston, MA, 02215, USA
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
| | - Henry Long
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Medical Oncology, Dana-Farber Cancer Institute, Harvard Medical School, Boston, MA, 02215, USA
| | - Tao Liu
- Department of Biostatistics and Bioinformatics, Roswell Park Comprehensive Cancer Center, Buffalo, NY, 14263, USA.
| | - X Shirley Liu
- Department of Data Science, Dana-Farber Cancer Institute, Harvard T.H. Chan School of Public Health, Boston, MA, 02215, USA.
- Center for Functional Cancer Epigenetics, Dana-Farber Cancer Institute, Boston, MA, 02215, USA.
| |
Collapse
|