1
|
Shepherdson JL, Friedman RZ, Zheng Y, Sun C, Oh IY, Granas DM, Cohen BA, Chen S, White MA. Pathogenic variants in CRX have distinct cis-regulatory effects on enhancers and silencers in photoreceptors. Genome Res 2024; 34:243-255. [PMID: 38355306 PMCID: PMC10984388 DOI: 10.1101/gr.278133.123] [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: 05/27/2023] [Accepted: 02/01/2024] [Indexed: 02/16/2024]
Abstract
Dozens of variants in the gene for the homeodomain transcription factor (TF) cone-rod homeobox (CRX) are linked with human blinding diseases that vary in their severity and age of onset. How different variants in this single TF alter its function in ways that lead to a range of phenotypes is unclear. We characterized the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in mouse retina explants carrying knock-ins of two variants, one in the DNA-binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation in these mutant Crx retinas corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, and p.E168d2 has distinct effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are derepressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci are partially predictive of episomal MPRA activity, and distal elements whose accessibility increases later in retinal development are enriched for CREs with silencer activity. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers while having a qualitatively different impact on silencers.
Collapse
Affiliation(s)
- James L Shepherdson
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
- Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Ryan Z Friedman
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
- Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Yiqiao Zheng
- Department of Ophthalmology and Visual Sciences, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Chi Sun
- Department of Ophthalmology and Visual Sciences, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Inez Y Oh
- Department of Ophthalmology and Visual Sciences, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - David M Granas
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
- Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Barak A Cohen
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
- Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Shiming Chen
- Department of Ophthalmology and Visual Sciences, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA;
- Department of Developmental Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| | - Michael A White
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA;
- Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri 63110, USA
| |
Collapse
|
2
|
Zheng Y, Stormo GD, Chen S. Aberrant homeodomain-DNA cooperative dimerization underlies distinct developmental defects in two dominant CRX retinopathy models. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.03.12.584677. [PMID: 38559186 PMCID: PMC10979960 DOI: 10.1101/2024.03.12.584677] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 04/04/2024]
Abstract
Paired-class homeodomain transcription factors (HD TFs) play essential roles in vertebrate development, and their mutations are linked to human diseases. One unique feature of paired-class HD is cooperative dimerization on specific palindrome DNA sequences. Yet, the functional significance of HD cooperative dimerization in animal development and its dysregulation in diseases remain elusive. Using the retinal TF Cone-rod Homeobox (CRX) as a model, we have studied how blindness-causing mutations in the paired HD, p.E80A and p.K88N, alter CRX's cooperative dimerization, lead to gene misexpression and photoreceptor developmental deficits in dominant manners. CRXE80A maintains binding at monomeric WT CRX motifs but is deficient in cooperative binding at dimeric motifs. CRXE80A's cooperativity defect impacts the exponential increase of photoreceptor gene expression in terminal differentiation and produces immature, non-functional photoreceptors in the CrxE80A retinas. CRXK88N is highly cooperative and localizes to ectopic genomic sites with strong enrichment of dimeric HD motifs. CRXK88N's altered biochemical properties disrupt CRX's ability to direct dynamic chromatin remodeling during development to activate photoreceptor differentiation programs and silence progenitor programs. Our study here provides in vitro and in vivo molecular evidence that paired-class HD cooperative dimerization regulates neuronal development and dysregulation of cooperative binding contributes to severe dominant blinding retinopathies.
Collapse
Affiliation(s)
- Yiqiao Zheng
- Molecular Genetics and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Washington University in St Louis, Saint Louis, Missouri, 63110, USA
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, 63110, USA
| | - Gary D. Stormo
- Department of Genetics, Washington University in St Louis, Saint Louis, Missouri, 63110, USA
| | - Shiming Chen
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, 63110, USA
- Department of Developmental Biology, Washington University in St Louis, Saint Louis, Missouri, 63110, USA
| |
Collapse
|
3
|
Zheng Y, Chen S. Transcriptional precision in photoreceptor development and diseases - Lessons from 25 years of CRX research. Front Cell Neurosci 2024; 18:1347436. [PMID: 38414750 PMCID: PMC10896975 DOI: 10.3389/fncel.2024.1347436] [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: 11/30/2023] [Accepted: 01/19/2024] [Indexed: 02/29/2024] Open
Abstract
The vertebrate retina is made up of six specialized neuronal cell types and one glia that are generated from a common retinal progenitor. The development of these distinct cell types is programmed by transcription factors that regulate the expression of specific genes essential for cell fate specification and differentiation. Because of the complex nature of transcriptional regulation, understanding transcription factor functions in development and disease is challenging. Research on the Cone-rod homeobox transcription factor CRX provides an excellent model to address these challenges. In this review, we reflect on 25 years of mammalian CRX research and discuss recent progress in elucidating the distinct pathogenic mechanisms of four CRX coding variant classes. We highlight how in vitro biochemical studies of CRX protein functions facilitate understanding CRX regulatory principles in animal models. We conclude with a brief discussion of the emerging systems biology approaches that could accelerate precision medicine for CRX-linked diseases and beyond.
Collapse
Affiliation(s)
- Yiqiao Zheng
- Molecular Genetics and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Saint Louis, MO, United States
- Department of Ophthalmology and Visual Sciences, Saint Louis, MO, United States
| | - Shiming Chen
- Molecular Genetics and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Saint Louis, MO, United States
- Department of Ophthalmology and Visual Sciences, Saint Louis, MO, United States
- Department of Developmental Biology, Washington University in St. Louis, Saint Louis, MO, United States
| |
Collapse
|
4
|
Loell KJ, Friedman RZ, Myers CA, Corbo JC, Cohen BA, White MA. Transcription factor interactions explain the context-dependent activity of CRX binding sites. PLoS Comput Biol 2024; 20:e1011802. [PMID: 38227575 PMCID: PMC10817189 DOI: 10.1371/journal.pcbi.1011802] [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: 03/06/2023] [Revised: 01/26/2024] [Accepted: 01/06/2024] [Indexed: 01/18/2024] Open
Abstract
The effects of transcription factor binding sites (TFBSs) on the activity of a cis-regulatory element (CRE) depend on the local sequence context. In rod photoreceptors, binding sites for the transcription factor (TF) Cone-rod homeobox (CRX) occur in both enhancers and silencers, but the sequence context that determines whether CRX binding sites contribute to activation or repression of transcription is not understood. To investigate the context-dependent activity of CRX sites, we fit neural network-based models to the activities of synthetic CREs composed of photoreceptor TFBSs. The models revealed that CRX binding sites consistently make positive, independent contributions to CRE activity, while negative homotypic interactions between sites cause CREs composed of multiple CRX sites to function as silencers. The effects of negative homotypic interactions can be overcome by the presence of other TFBSs that either interact cooperatively with CRX sites or make independent positive contributions to activity. The context-dependent activity of CRX sites is thus determined by the balance between positive heterotypic interactions, independent contributions of TFBSs, and negative homotypic interactions. Our findings explain observed patterns of activity among genomic CRX-bound enhancers and silencers, and suggest that enhancers may require diverse TFBSs to overcome negative homotypic interactions between TFBSs.
Collapse
Affiliation(s)
- Kaiser J. Loell
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
| | - Ryan Z. Friedman
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
| | - Connie A. Myers
- Department of Pathology and Immunology, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
| | - Joseph C. Corbo
- Department of Pathology and Immunology, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
| | - Barak A. Cohen
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
| | - Michael A. White
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine in St. Louis, St. Louis, Missouri, United States of America
| |
Collapse
|
5
|
Pan D, Zhang X, Jin K, Jin ZB. CRX haploinsufficiency compromises photoreceptor precursor translocation and differentiation in human retinal organoids. Stem Cell Res Ther 2023; 14:346. [PMID: 38049871 PMCID: PMC10696917 DOI: 10.1186/s13287-023-03590-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2023] [Accepted: 11/28/2023] [Indexed: 12/06/2023] Open
Abstract
BACKGROUND The CRX-associated autosomal dominant retinopathies suggest a possible pathogenic mechanism of gene haploinsufficiency. However, based on reported human patient cases and studies with mouse models, it is hard to confirm the specific weight of haploinsufficiency in pathogenesis due to the interspecies gaps between gene expression and function. METHODS We created monoallelic CRX by replacing one allele with tdTomato in human embryonic stem cells (hESCs) and subsequently dissect pathogenesis in hESCs-derived retinal organoids. We used transcriptome and immunofluorescence analyses to dissect phenotypic differences between CRX-monoallelic knockout and control wildtype organoids. For location analysis of CRX+ cells, a CRX-expression-tracing system was constructed in control hESCs. We implemented long-term live-cell imaging to describe the translocation of CRX+ cells between two groups in early organoid differentiation. The expression pattern of these dynamic differences was validated using RNA-seq and immunofluorescence assays. RESULTS We identified delayed differentiation of outer nuclear layer (ONL) stratification along with thinner ONL, serious loss of photoreceptor outer segments, as well as downregulated expression of gene for phototransduction and inner/outer segment formation. By live-cell imaging and immunostaining, we observed the overtension of actomyosin network and the arrested translocation of monoallelic CRX+ cells in the early stage of retinal differentiation. CONCLUSIONS We confirmed that gene haploinsufficiency is the mechanism for the dominant pathogenicity of CRX and discovered that CRX regulated postmitotic photoreceptor precursor translocation in addition to its specification of photoreceptor cell fates during human retinal development. These findings revealed a new underlying mechanism of CRX dominant pathogenesis and provided a new clue for the treatment of CRX-associated human retinopathies.
Collapse
Affiliation(s)
- Deng Pan
- Beijing Institute of Ophthalmology, Beijing Tongren Hospital, Capital Medical University, Beijing, 100730, China.
| | - Xiao Zhang
- Beijing Institute of Ophthalmology, Beijing Tongren Hospital, Capital Medical University, Beijing, 100730, China
| | - Kangxin Jin
- Beijing Institute of Ophthalmology, Beijing Tongren Hospital, Capital Medical University, Beijing, 100730, China
| | - Zi-Bing Jin
- Beijing Institute of Ophthalmology, Beijing Tongren Hospital, Capital Medical University, Beijing, 100730, China.
| |
Collapse
|
6
|
Shepherdson JL, Friedman RZ, Zheng Y, Sun C, Oh IY, Granas DM, Cohen BA, Chen S, White MA. Pathogenic variants in Crx have distinct cis-regulatory effects on enhancers and silencers in photoreceptors. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.05.27.542576. [PMID: 37292699 PMCID: PMC10245955 DOI: 10.1101/2023.05.27.542576] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Dozens of variants in the photoreceptor-specific transcription factor (TF) CRX are linked with human blinding diseases that vary in their severity and age of onset. It is unclear how different variants in this single TF alter its function in ways that lead to a range of phenotypes. We examined the effects of human disease-causing variants on CRX cis-regulatory function by deploying massively parallel reporter assays (MPRAs) in live mouse retinas carrying knock-ins of two variants, one in the DNA binding domain (p.R90W) and the other in the transcriptional effector domain (p.E168d2). The degree of reporter gene dysregulation caused by the variants corresponds with their phenotypic severity. The two variants affect similar sets of enhancers, while p.E168d2 has stronger effects on silencers. Cis-regulatory elements (CREs) near cone photoreceptor genes are enriched for silencers that are de-repressed in the presence of p.E168d2. Chromatin environments of CRX-bound loci were partially predictive of episomal MPRA activity, and silencers were notably enriched among distal elements whose accessibility increases later in retinal development. We identified a set of potentially pleiotropic regulatory elements that convert from silencers to enhancers in retinas that lack a functional CRX effector domain. Our findings show that phenotypically distinct variants in different domains of CRX have partially overlapping effects on its cis-regulatory function, leading to misregulation of similar sets of enhancers, while having a qualitatively different impact on silencers.
Collapse
Affiliation(s)
- James L. Shepherdson
- Department of Genetics
- Edison Family Center for Genome Sciences & Systems Biology
| | - Ryan Z. Friedman
- Department of Genetics
- Edison Family Center for Genome Sciences & Systems Biology
| | | | - Chi Sun
- Department of Ophthalmology and Visual Sciences
| | - Inez Y. Oh
- Department of Ophthalmology and Visual Sciences
| | - David M. Granas
- Department of Genetics
- Edison Family Center for Genome Sciences & Systems Biology
| | - Barak A. Cohen
- Department of Genetics
- Edison Family Center for Genome Sciences & Systems Biology
| | - Shiming Chen
- Department of Ophthalmology and Visual Sciences
- Department of Developmental Biology, Washington University School of Medicine in St. Louis, St. Louis, MO 63110
| | - Michael A. White
- Department of Genetics
- Edison Family Center for Genome Sciences & Systems Biology
| |
Collapse
|
7
|
Li J, Wang J, Ibarra IL, Cheng X, Luecken MD, Lu J, Monavarfeshani A, Yan W, Zheng Y, Zuo Z, Colborn SLZ, Cortez BS, Owen LA, Tran NM, Shekhar K, Sanes JR, Stout JT, Chen S, Li Y, DeAngelis MM, Theis FJ, Chen R. Integrated multi-omics single cell atlas of the human retina. RESEARCH SQUARE 2023:rs.3.rs-3471275. [PMID: 38014002 PMCID: PMC10680922 DOI: 10.21203/rs.3.rs-3471275/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Single-cell sequencing has revolutionized the scale and resolution of molecular profiling of tissues and organs. Here, we present an integrated multimodal reference atlas of the most accessible portion of the mammalian central nervous system, the retina. We compiled around 2.4 million cells from 55 donors, including 1.4 million unpublished data points, to create a comprehensive human retina cell atlas (HRCA) of transcriptome and chromatin accessibility, unveiling over 110 types. Engaging the retina community, we annotated each cluster, refined the Cell Ontology for the retina, identified distinct marker genes, and characterized cis-regulatory elements and gene regulatory networks (GRNs) for these cell types. Our analysis uncovered intriguing differences in transcriptome, chromatin, and GRNs across cell types. In addition, we modeled changes in gene expression and chromatin openness across gender and age. This integrated atlas also enabled the fine-mapping of GWAS and eQTL variants. Accessible through interactive browsers, this multimodal cross-donor and cross-lab HRCA, can facilitate a better understanding of retinal function and pathology.
Collapse
Affiliation(s)
- Jin Li
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas, United States
| | - Jun Wang
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas, United States
| | - Ignacio L Ibarra
- Institute of Computational Biology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
| | - Xuesen Cheng
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas, United States
| | - Malte D Luecken
- Institute of Computational Biology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
- Institute of Lung Health & Immunity, Helmholtz Munich; Member of the German Center for Lung Research (DZL), Munich, Germany
| | - Jiaxiong Lu
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas, United States
- Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, Texas, United States
| | - Aboozar Monavarfeshani
- Center for Brain Science and Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
| | - Wenjun Yan
- Center for Brain Science and Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
| | - Yiqiao Zheng
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, United States
| | - Zhen Zuo
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
| | | | | | - Leah A Owen
- John A. Moran Eye Center, Department of Ophthalmology and Visual Sciences, University of Utah, Salt Lake City, Utah, United States
| | - Nicholas M Tran
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
| | - Karthik Shekhar
- Department of Chemical and Biomolecular Engineering; Helen Wills Neuroscience Institute; Center for Computational Biology; California Institute for Quantitative Biosciences, QB3, University of California, Berkeley, Berkeley, California, United States
| | - Joshua R Sanes
- Center for Brain Science and Department of Molecular and Cellular Biology, Harvard University, Cambridge, United States
| | - J Timothy Stout
- Department of Ophthalmology, Cullen Eye Institute, Baylor College of Medicine, Houston, Texas, United States
| | - Shiming Chen
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, United States
- Department of Developmental Biology, Washington University in St Louis, Saint Louis, Missouri, United States
| | - Yumei Li
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas, United States
| | - Margaret M DeAngelis
- Department of Ophthalmology, Ross Eye Institute, Jacobs School of Medicine and Biomedical Sciences, State University of New York at Buffalo, Buffalo, New York, United States
| | - Fabian J Theis
- Institute of Computational Biology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
| | - Rui Chen
- Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, Texas, United States
- Human Genome Sequencing Center, Baylor College of Medicine, Houston, Texas, United States
- Department of Biochemistry and Molecular Biology, Baylor College of Medicine, Houston, Texas, United States
| |
Collapse
|
8
|
Zheng Y, Sun C, Zhang X, Ruzycki PA, Chen S. Missense mutations in CRX homeodomain cause dominant retinopathies through two distinct mechanisms. eLife 2023; 12:RP87147. [PMID: 37963072 PMCID: PMC10645426 DOI: 10.7554/elife.87147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2023] Open
Abstract
Homeodomain transcription factors (HD TFs) are instrumental to vertebrate development. Mutations in HD TFs have been linked to human diseases, but their pathogenic mechanisms remain elusive. Here, we use Cone-Rod Homeobox (CRX) as a model to decipher the disease-causing mechanisms of two HD mutations, p.E80A and p.K88N, that produce severe dominant retinopathies. Through integrated analysis of molecular and functional evidence in vitro and in knock-in mouse models, we uncover two novel gain-of-function mechanisms: p.E80A increases CRX-mediated transactivation of canonical CRX target genes in developing photoreceptors; p.K88N alters CRX DNA-binding specificity resulting in binding at ectopic sites and severe perturbation of CRX target gene expression. Both mechanisms produce novel retinal morphological defects and hinder photoreceptor maturation distinct from loss-of-function models. This study reveals the distinct roles of E80 and K88 residues in CRX HD regulatory functions and emphasizes the importance of transcriptional precision in normal development.
Collapse
Affiliation(s)
- Yiqiao Zheng
- Molecular Genetic and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Washington University in St LouisSaint LouisUnited States
- Department of Ophthalmology and Visual Sciences, Washington University in St LouisSaint LouisUnited States
| | - Chi Sun
- Molecular Genetic and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Washington University in St LouisSaint LouisUnited States
- Department of Ophthalmology and Visual Sciences, Washington University in St LouisSaint LouisUnited States
| | - Xiaodong Zhang
- Department of Ophthalmology and Visual Sciences, Washington University in St LouisSaint LouisUnited States
| | - Philip A Ruzycki
- Department of Ophthalmology and Visual Sciences, Washington University in St LouisSaint LouisUnited States
- Department of Genetics, Washington University in St LouisSaint LouisUnited States
| | - Shiming Chen
- Molecular Genetic and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Washington University in St LouisSaint LouisUnited States
- Department of Ophthalmology and Visual Sciences, Washington University in St LouisSaint LouisUnited States
- Department of Developmental Biology, Washington University in St LouisSaint LouisUnited States
| |
Collapse
|
9
|
Friedman RZ, Ramu A, Lichtarge S, Myers CA, Granas DM, Gause M, Corbo JC, Cohen BA, White MA. Active learning of enhancer and silencer regulatory grammar in photoreceptors. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.21.554146. [PMID: 37662358 PMCID: PMC10473580 DOI: 10.1101/2023.08.21.554146] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/05/2023]
Abstract
Cis-regulatory elements (CREs) direct gene expression in health and disease, and models that can accurately predict their activities from DNA sequences are crucial for biomedicine. Deep learning represents one emerging strategy to model the regulatory grammar that relates CRE sequence to function. However, these models require training data on a scale that exceeds the number of CREs in the genome. We address this problem using active machine learning to iteratively train models on multiple rounds of synthetic DNA sequences assayed in live mammalian retinas. During each round of training the model actively selects sequence perturbations to assay, thereby efficiently generating informative training data. We iteratively trained a model that predicts the activities of sequences containing binding motifs for the photoreceptor transcription factor Cone-rod homeobox (CRX) using an order of magnitude less training data than current approaches. The model's internal confidence estimates of its predictions are reliable guides for designing sequences with high activity. The model correctly identified critical sequence differences between active and inactive sequences with nearly identical transcription factor binding sites, and revealed order and spacing preferences for combinations of motifs. Our results establish active learning as an effective method to train accurate deep learning models of cis-regulatory function after exhausting naturally occurring training examples in the genome.
Collapse
Affiliation(s)
- Ryan Z. Friedman
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine, Saint Louis, MO, 63110
- Department of Genetics, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Avinash Ramu
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine, Saint Louis, MO, 63110
- Department of Genetics, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Sara Lichtarge
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine, Saint Louis, MO, 63110
- Department of Genetics, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Connie A. Myers
- Department of Pathology and Immunology, Washington University School of Medicine, Saint Louis, MO, 63110
| | - David M. Granas
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine, Saint Louis, MO, 63110
- Department of Genetics, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Maria Gause
- Department of Pathology and Immunology, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Joseph C. Corbo
- Department of Pathology and Immunology, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Barak A. Cohen
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine, Saint Louis, MO, 63110
- Department of Genetics, Washington University School of Medicine, Saint Louis, MO, 63110
| | - Michael A. White
- The Edison Family Center for Genome Sciences & Systems Biology, Washington University School of Medicine, Saint Louis, MO, 63110
- Department of Genetics, Washington University School of Medicine, Saint Louis, MO, 63110
| |
Collapse
|
10
|
Cain B, Webb J, Yuan Z, Cheung D, Lim HW, Kovall R, Weirauch MT, Gebelein B. Prediction of cooperative homeodomain DNA binding sites from high-throughput-SELEX data. Nucleic Acids Res 2023; 51:6055-6072. [PMID: 37114997 PMCID: PMC10325903 DOI: 10.1093/nar/gkad318] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Revised: 04/12/2023] [Accepted: 04/25/2023] [Indexed: 04/29/2023] Open
Abstract
Homeodomain proteins constitute one of the largest families of metazoan transcription factors. Genetic studies have demonstrated that homeodomain proteins regulate many developmental processes. Yet, biochemical data reveal that most bind highly similar DNA sequences. Defining how homeodomain proteins achieve DNA binding specificity has therefore been a long-standing goal. Here, we developed a novel computational approach to predict cooperative dimeric binding of homeodomain proteins using High-Throughput (HT) SELEX data. Importantly, we found that 15 of 88 homeodomain factors form cooperative homodimer complexes on DNA sites with precise spacing requirements. Approximately one third of the paired-like homeodomain proteins cooperatively bind palindromic sequences spaced 3 bp apart, whereas other homeodomain proteins cooperatively bind sites with distinct orientation and spacing requirements. Combining structural models of a paired-like factor with our cooperativity predictions identified key amino acid differences that help differentiate between cooperative and non-cooperative factors. Finally, we confirmed predicted cooperative dimer sites in vivo using available genomic data for a subset of factors. These findings demonstrate how HT-SELEX data can be computationally mined to predict cooperativity. In addition, the binding site spacing requirements of select homeodomain proteins provide a mechanism by which seemingly similar AT-rich DNA sequences can preferentially recruit specific homeodomain factors.
Collapse
Affiliation(s)
- Brittany Cain
- Department of Biomedical Engineering, University of Cincinnati, Cincinnati, OH 45221, USA
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, 3333 Burnet Ave, MLC 7007, Cincinnati, OH 45229, USA
| | - Jordan Webb
- Department of Molecular Genetics, Biochemistry and Microbiology, University of Cincinnati College of Medicine, Cincinnati, OH 45267, USA
| | - Zhenyu Yuan
- Department of Molecular Genetics, Biochemistry and Microbiology, University of Cincinnati College of Medicine, Cincinnati, OH 45267, USA
| | - David Cheung
- Graduate Program in Molecular and Developmental Biology, Cincinnati Children's Hospital Research Foundation, Cincinnati, OH 45229, USA
| | - Hee-Woong Lim
- Division of Biomedical Informatics, Cincinnati Children's Hospital Medical Center, Cincinnati, OH 45229, USA
- Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, OH 45229, USA
| | - Rhett A Kovall
- Department of Molecular Genetics, Biochemistry and Microbiology, University of Cincinnati College of Medicine, Cincinnati, OH 45267, USA
| | - Matthew T Weirauch
- Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, OH 45229, USA
- Divisions of Human Genetics, Biomedical Informatics and Developmental Biology, Center for Autoimmune Genomics and Etiology (CAGE), Cincinnati Children's Hospital Medical Center, Cincinnati, OH 45229, USA
| | - Brian Gebelein
- Division of Developmental Biology, Cincinnati Children's Hospital Medical Center, 3333 Burnet Ave, MLC 7007, Cincinnati, OH 45229, USA
- Department of Pediatrics, University of Cincinnati College of Medicine, Cincinnati, OH 45229, USA
| |
Collapse
|
11
|
Zheng Y, Sun C, Zhang X, Ruzycki PA, Chen S. Missense mutations in CRX homeodomain cause dominant retinopathies through two distinct mechanisms. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.01.526652. [PMID: 36778408 PMCID: PMC9915647 DOI: 10.1101/2023.02.01.526652] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Homeodomain transcription factors (HD TFs) are instrumental to vertebrate development. Mutations in HD TFs have been linked to human diseases, but their pathogenic mechanisms remain elusive. Here we use Cone-Rod Homeobox (CRX) as a model to decipher the disease-causing mechanisms of two HD mutations, p.E80A and p.K88N, that produce severe dominant retinopathies. Through integrated analysis of molecular and functional evidence in vitro and in knock-in mouse models, we uncover two novel gain-of-function mechanisms: p.E80A increases CRX-mediated transactivation of canonical CRX target genes in developing photoreceptors; p.K88N alters CRX DNA-binding specificity resulting in binding at ectopic sites and severe perturbation of CRX target gene expression. Both mechanisms produce novel retinal morphological defects and hinder photoreceptor maturation distinct from loss-of-function models. This study reveals the distinct roles of E80 and K88 residues in CRX HD regulatory functions and emphasizes the importance of transcriptional precision in normal development.
Collapse
Affiliation(s)
- Yiqiao Zheng
- Molecular Genetic and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
| | - Chi Sun
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
| | - Xiaodong Zhang
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
| | - Philip A. Ruzycki
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
- Department of Genetics, Washington University in St Louis, Saint Louis, Missouri, USA
| | - Shiming Chen
- Molecular Genetic and Genomics Graduate Program, Division of Biological and Biomedical Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
- Department of Ophthalmology and Visual Sciences, Washington University in St Louis, Saint Louis, Missouri, USA
- Department of Developmental Biology, Washington University in St Louis, Saint Louis, Missouri, USA
| |
Collapse
|
12
|
Zheng Y, VanDusen NJ. Massively Parallel Reporter Assays for High-Throughput In Vivo Analysis of Cis-Regulatory Elements. J Cardiovasc Dev Dis 2023; 10:jcdd10040144. [PMID: 37103023 PMCID: PMC10146671 DOI: 10.3390/jcdd10040144] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 03/24/2023] [Accepted: 03/27/2023] [Indexed: 03/31/2023] Open
Abstract
The rapid improvement of descriptive genomic technologies has fueled a dramatic increase in hypothesized connections between cardiovascular gene expression and phenotypes. However, in vivo testing of these hypotheses has predominantly been relegated to slow, expensive, and linear generation of genetically modified mice. In the study of genomic cis-regulatory elements, generation of mice featuring transgenic reporters or cis-regulatory element knockout remains the standard approach. While the data obtained is of high quality, the approach is insufficient to keep pace with candidate identification and therefore results in biases introduced during the selection of candidates for validation. However, recent advances across a range of disciplines are converging to enable functional genomic assays that can be conducted in a high-throughput manner. Here, we review one such method, massively parallel reporter assays (MPRAs), in which the activities of thousands of candidate genomic regulatory elements are simultaneously assessed via the next-generation sequencing of a barcoded reporter transcript. We discuss best practices for MPRA design and use, with a focus on practical considerations, and review how this emerging technology has been successfully deployed in vivo. Finally, we discuss how MPRAs are likely to evolve and be used in future cardiovascular research.
Collapse
|
13
|
A single-cell massively parallel reporter assay detects cell-type-specific gene regulation. Nat Genet 2023; 55:346-354. [PMID: 36635387 PMCID: PMC9931678 DOI: 10.1038/s41588-022-01278-7] [Citation(s) in RCA: 24] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 12/05/2022] [Indexed: 01/14/2023]
Abstract
Massively parallel reporter gene assays are key tools in regulatory genomics but cannot be used to identify cell-type-specific regulatory elements without performing assays serially across different cell types. To address this problem, we developed a single-cell massively parallel reporter assay (scMPRA) to measure the activity of libraries of cis-regulatory sequences (CRSs) across multiple cell types simultaneously. We assayed a library of core promoters in a mixture of HEK293 and K562 cells and showed that scMPRA is a reproducible, highly parallel, single-cell reporter gene assay that detects cell-type-specific cis-regulatory activity. We then measured a library of promoter variants across multiple cell types in live mouse retinas and showed that subtle genetic variants can produce cell-type-specific effects on cis-regulatory activity. We anticipate that scMPRA will be widely applicable for studying the role of CRSs across diverse cell types.
Collapse
|
14
|
Bian F, Daghsni M, Lu F, Liu S, Gross JM, Aldiri I. Functional analysis of the Vsx2 super-enhancer uncovers distinct cis-regulatory circuits controlling Vsx2 expression during retinogenesis. Development 2022; 149:dev200642. [PMID: 35831950 PMCID: PMC9440754 DOI: 10.1242/dev.200642] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 06/24/2022] [Indexed: 11/20/2022]
Abstract
Vsx2 is a transcription factor essential for retinal proliferation and bipolar cell differentiation, but the molecular mechanisms underlying its developmental roles are unclear. Here, we have profiled VSX2 genomic occupancy during mouse retinogenesis, revealing extensive retinal genetic programs associated with VSX2 during development. VSX2 binds and transactivates its enhancer in association with the transcription factor PAX6. Mice harboring deletions in the Vsx2 regulatory landscape exhibit specific abnormalities in retinal proliferation and in bipolar cell differentiation. In one of those deletions, a complete loss of bipolar cells is associated with a bias towards photoreceptor production. VSX2 occupies cis-regulatory elements nearby genes associated with photoreceptor differentiation and homeostasis in the adult mouse and human retina, including a conserved region nearby Prdm1, a factor implicated in the specification of rod photoreceptors and suppression of bipolar cell fate. VSX2 interacts with the transcription factor OTX2 and can act to suppress OTX2-dependent enhancer transactivation of the Prdm1 enhancer. Taken together, our analyses indicate that Vsx2 expression can be temporally and spatially uncoupled at the enhancer level, and they illuminate important mechanistic insights into how VSX2 is engaged with gene regulatory networks that are essential for retinal proliferation and cell fate acquisition.
Collapse
Affiliation(s)
- Fuyun Bian
- Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| | - Marwa Daghsni
- Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| | - Fangfang Lu
- Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
- Department of Ophthalmology, The Second Xiangya Hospital, Central South University, Changsha, Hunan 410011, China
| | - Silvia Liu
- Department of Pathology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| | - Jeffrey M Gross
- Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
- Department of Developmental Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
- Louis J. Fox Center for Vision Restoration, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| | - Issam Aldiri
- Department of Ophthalmology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
- Department of Developmental Biology, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
- Louis J. Fox Center for Vision Restoration, University of Pittsburgh School of Medicine, Pittsburgh, PA 15213, USA
| |
Collapse
|
15
|
VandenBosch LS, Luu K, Timms AE, Challam S, Wu Y, Lee AY, Cherry TJ. Machine Learning Prediction of Non-Coding Variant Impact in Human Retinal cis-Regulatory Elements. Transl Vis Sci Technol 2022; 11:16. [PMID: 35435921 PMCID: PMC9034719 DOI: 10.1167/tvst.11.4.16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2021] [Accepted: 03/25/2022] [Indexed: 11/24/2022] Open
Abstract
Purpose Prior studies have demonstrated the significance of specific cis-regulatory variants in retinal disease; however, determining the functional impact of regulatory variants remains a major challenge. In this study, we utilized a machine learning approach, trained on epigenomic data from the adult human retina, to systematically quantify the predicted impact of cis-regulatory variants. Methods We used human retinal DNA accessibility data (ATAC-seq) to determine a set of 18.9k high-confidence, putative cis-regulatory elements. Eighty percent of these elements were used to train a machine learning model utilizing a gapped k-mer support vector machine-based approach. In silico saturation mutagenesis and variant scoring was applied to predict the functional impact of all potential single nucleotide variants within cis-regulatory elements. Impact scores were tested in a 20% hold-out dataset and compared to allele population frequency, phylogenetic conservation, transcription factor (TF) binding motifs, and existing massively parallel reporter assay data. Results We generated a model that distinguishes between human retinal regulatory elements and negative test sequences with 95% accuracy. Among a hold-out test set of 3.7k human retinal CREs, all possible single nucleotide variants were scored. Variants with negative impact scores correlated with higher phylogenetic conservation of the reference allele, disruption of predicted TF binding motifs, and massively parallel reporter expression. Conclusions We demonstrated the utility of human retinal epigenomic data to train a machine learning model for the purpose of predicting the impact of non-coding regulatory sequence variants. Our model accurately scored sequences and predicted putative transcription factor binding motifs. This approach has the potential to expedite the characterization of pathogenic non-coding sequence variants in the context of unexplained retinal disease. Translational Relevance This workflow and resulting dataset serve as a promising genomic tool to facilitate the clinical prioritization of functionally disruptive non-coding mutations in the retina.
Collapse
Affiliation(s)
- Leah S. VandenBosch
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
| | - Kelsey Luu
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
| | - Andrew E. Timms
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
| | - Shriya Challam
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
| | - Yue Wu
- University of Washington Department of Ophthalmology, Seattle, WA, USA
| | - Aaron Y. Lee
- University of Washington Department of Ophthalmology, Seattle, WA, USA
- Brotman Baty Institute for Precision Medicine, Seattle, WA, USA
| | - Timothy J. Cherry
- Center for Developmental Biology and Regenerative Medicine, Seattle Children's Research Institute, Seattle, WA, USA
- Brotman Baty Institute for Precision Medicine, Seattle, WA, USA
- University of Washington Department of Pediatrics, Seattle, WA, USA
| |
Collapse
|
16
|
Shimai K, Veeman M. Quantitative Dissection of the Proximal Ciona brachyury Enhancer. Front Cell Dev Biol 2022; 9:804032. [PMID: 35127721 PMCID: PMC8814421 DOI: 10.3389/fcell.2021.804032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/28/2021] [Accepted: 12/27/2021] [Indexed: 11/13/2022] Open
Abstract
A major goal in biology is to understand the rules by which cis-regulatory sequences control spatially and temporally precise expression patterns. Here we present a systematic dissection of the proximal enhancer for the notochord-specific transcription factor brachyury in the ascidian chordate Ciona. The study uses a quantitative image-based reporter assay that incorporates a dual-reporter strategy to control for variable electroporation efficiency. We identified and mutated multiple predicted transcription factor binding sites of interest based on statistical matches to the JASPAR binding motif database. Most sites (Zic, Ets, FoxA, RBPJ) were selected based on prior knowledge of cell fate specification in both the primary and secondary notochord. We also mutated predicted Brachyury sites to investigate potential autoregulation as well as Fos/Jun (AP1) sites that had very strong matches to JASPAR. Our goal was to quantitatively define the relative importance of these different sites, to explore the importance of predicted high-affinity versus low-affinity motifs, and to attempt to design mutant enhancers that were specifically expressed in only the primary or secondary notochord lineages. We found that the mutation of all predicted high-affinity sites for Zic, FoxA or Ets led to quantifiably distinct effects. The FoxA construct caused a severe loss of reporter expression whereas the Ets construct had little effect. A strong Ets phenotype was only seen when much lower-scoring binding sites were also mutated. This supports the enhancer suboptimization hypothesis proposed by Farley and Levine but suggests that it may only apply to some but not all transcription factor families. We quantified reporter expression separately in the two notochord lineages with the expectation that Ets mutations and RBPJ mutations would have distinct effects given that primary notochord is induced by Ets-mediated FGF signaling whereas secondary notochord is induced by RBPJ/Su(H)-mediated Notch/Delta signaling. We found, however, that ETS mutations affected primary and secondary notochord expression relatively equally and that RBPJ mutations were only moderately more severe in their effect on secondary versus primary notochord. Our results point to the promise of quantitative reporter assays for understanding cis-regulatory logic but also highlight the challenge of arbitrary statistical thresholds for predicting potentially important sites.
Collapse
|
17
|
Waters CT, Gisselbrecht SS, Sytnikova YA, Cafarelli TM, Hill DE, Bulyk ML. Quantitative-enhancer-FACS-seq (QeFS) reveals epistatic interactions among motifs within transcriptional enhancers in developing Drosophila tissue. Genome Biol 2021; 22:348. [PMID: 34930411 PMCID: PMC8686523 DOI: 10.1186/s13059-021-02574-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2020] [Accepted: 12/10/2021] [Indexed: 11/16/2022] Open
Abstract
Understanding the contributions of transcription factor DNA binding sites to transcriptional enhancers is a significant challenge. We developed Quantitative enhancer-FACS-Seq for highly parallel quantification of enhancer activities from a genomically integrated reporter in Drosophila melanogaster embryos. We investigate the contributions of the DNA binding motifs of four poorly characterized TFs to the activities of twelve embryonic mesodermal enhancers. We measure quantitative changes in enhancer activity and discover a range of epistatic interactions among the motifs, both synergistic and alleviating. We find that understanding the regulatory consequences of TF binding motifs requires that they be investigated in combination across enhancer contexts.
Collapse
Affiliation(s)
- Colin T Waters
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, 02115, USA
- Program in Biological and Biomedical Sciences, Harvard University, Cambridge, MA, 02138, USA
| | - Stephen S Gisselbrecht
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, 02115, USA
| | - Yuliya A Sytnikova
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, 02115, USA
| | - Tiziana M Cafarelli
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Genetics, Blavatnik Institute, Harvard Medical School, Boston, MA, 02115, USA
| | - David E Hill
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Cancer Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA
- Department of Genetics, Blavatnik Institute, Harvard Medical School, Boston, MA, 02115, USA
| | - Martha L Bulyk
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, 02115, USA.
- Program in Biological and Biomedical Sciences, Harvard University, Cambridge, MA, 02138, USA.
- Center for Cancer Systems Biology, Dana-Farber Cancer Institute, Boston, MA, 02215, USA.
- Department of Pathology, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, 02115, USA.
| |
Collapse
|
18
|
Friedman RZ, Granas DM, Myers CA, Corbo JC, Cohen BA, White MA. Information content differentiates enhancers from silencers in mouse photoreceptors. eLife 2021; 10:67403. [PMID: 34486522 PMCID: PMC8492058 DOI: 10.7554/elife.67403] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2021] [Accepted: 09/03/2021] [Indexed: 12/12/2022] Open
Abstract
Enhancers and silencers often depend on the same transcription factors (TFs) and are conflated in genomic assays of TF binding or chromatin state. To identify sequence features that distinguish enhancers and silencers, we assayed massively parallel reporter libraries of genomic sequences targeted by the photoreceptor TF cone-rod homeobox (CRX) in mouse retinas. Both enhancers and silencers contain more TF motifs than inactive sequences, but relative to silencers, enhancers contain motifs from a more diverse collection of TFs. We developed a measure of information content that describes the number and diversity of motifs in a sequence and found that, while both enhancers and silencers depend on CRX motifs, enhancers have higher information content. The ability of information content to distinguish enhancers and silencers targeted by the same TF illustrates how motif context determines the activity of cis-regulatory sequences. Different cell types are established by activating and repressing the activity of specific sets of genes, a process controlled by proteins called transcription factors. Transcription factors work by recognizing and binding short stretches of DNA in parts of the genome called cis-regulatory sequences. A cis-regulatory sequence that increases the activity of a gene when bound by transcription factors is called an enhancer, while a sequence that causes a decrease in gene activity is called a silencer. To establish a cell type, a particular transcription factor will act on both enhancers and silencers that control the activity of different genes. For example, the transcription factor cone-rod homeobox (CRX) is critical for specifying different types of cells in the retina, and it acts on both enhancers and silencers. In rod photoreceptors, CRX activates rod genes by binding their enhancers, while repressing cone photoreceptor genes by binding their silencers. However, CRX always recognizes and binds to the same DNA sequence, known as its binding site, making it unclear why some cis-regulatory sequences bound to CRX act as silencers, while others act as enhancers. Friedman et al. sought to understand how enhancers and silencers, both bound by CRX, can have different effects on the genes they control. Since both enhancers and silencers contain CRX binding sites, the difference between the two must lie in the sequence of the DNA surrounding these binding sites. Using retinas that have been explanted from mice and kept alive in the laboratory, Friedman et al. tested the activity of thousands of CRX-binding sequences from the mouse genome. This showed that both enhancers and silencers have more copies of CRX-binding sites than sequences of the genome that are inactive. Additionally, the results revealed that enhancers have a diverse collection of binding sites for other transcription factors, while silencers do not. Friedman et al. developed a new metric they called information content, which captures the diverse combinations of different transcription binding sites that cis-regulatory sequences can have. Using this metric, Friedman et al. showed that it is possible to distinguish enhancers from silencers based on their information content. It is critical to understand how the DNA sequences of cis-regulatory regions determine their activity, because mutations in these regions of the genome can cause disease. However, since every person has thousands of benign mutations in cis-regulatory sequences, it is a challenge to identify specific disease-causing mutations, which are relatively rare. One long-term goal of models of enhancers and silencers, such as Friedman et al.’s information content model, is to understand how mutations can affect cis-regulatory sequences, and, in some cases, lead to disease.
Collapse
Affiliation(s)
- Ryan Z Friedman
- Edison Family Center for Genome Sciences and Systems Biology, Washington University School of Medicine, St. Louis, United States.,Department of Genetics, Washington University School of Medicine, St. Louis, United States
| | - David M Granas
- Edison Family Center for Genome Sciences and Systems Biology, Washington University School of Medicine, St. Louis, United States.,Department of Genetics, Washington University School of Medicine, St. Louis, United States
| | - Connie A Myers
- Department of Pathology and Immunology, Washington University School of Medicine, St Louis, United States
| | - Joseph C Corbo
- Department of Pathology and Immunology, Washington University School of Medicine, St Louis, United States
| | - Barak A Cohen
- Edison Family Center for Genome Sciences and Systems Biology, Washington University School of Medicine, St. Louis, United States.,Department of Genetics, Washington University School of Medicine, St. Louis, United States
| | - Michael A White
- Edison Family Center for Genome Sciences and Systems Biology, Washington University School of Medicine, St. Louis, United States.,Department of Genetics, Washington University School of Medicine, St. Louis, United States
| |
Collapse
|
19
|
Guerrero-Santoro J, Khor JM, Açıkbaş AH, Jaynes JB, Ettensohn CA. Analysis of the DNA-binding properties of Alx1, an evolutionarily conserved regulator of skeletogenesis in echinoderms. J Biol Chem 2021; 297:100901. [PMID: 34157281 PMCID: PMC8319359 DOI: 10.1016/j.jbc.2021.100901] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Revised: 06/09/2021] [Accepted: 06/18/2021] [Indexed: 11/24/2022] Open
Abstract
Alx1, a homeodomain-containing transcription factor, is a highly conserved regulator of skeletogenesis in echinoderms. In sea urchins, Alx1 plays a central role in the differentiation of embryonic primary mesenchyme cells (PMCs) and positively regulates the transcription of most biomineralization genes expressed by these cells. The alx1 gene arose via duplication and acquired a skeletogenic function distinct from its paralog (alx4) through the exonization of a 41–amino acid motif (the D2 domain). Alx1 and Alx4 contain glutamine-50 paired-type homeodomains, which interact preferentially with palindromic binding sites in vitro. Chromatin immunoprecipitation sequencing (ChIP-seq) studies have shown, however, that Alx1 binds both to palindromic and half sites in vivo. To address this apparent discrepancy and explore the function of the D2 domain, we used an endogenous cis-regulatory module associated with Sp-mtmmpb, a gene that encodes a PMC-specific metalloprotease, to analyze the DNA-binding properties of Alx1. We find that Alx1 forms dimeric complexes on TAAT-containing half sites by a mechanism distinct from the well-known mechanism of dimerization on palindromic sites. We used transgenic reporter assays to analyze the functional roles of half sites in vivo and demonstrate that two sites with partially redundant functions are essential for the PMC-specific activity of the Sp-mtmmpb cis-regulatory module. Finally, we show that the D2 domain influences the DNA-binding properties of Alx1 in vitro, suggesting that the exonization of this motif may have facilitated the acquisition of new transcriptional targets and consequently a novel developmental function.
Collapse
Affiliation(s)
| | - Jian Ming Khor
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA
| | - Ayşe Haruka Açıkbaş
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA
| | - James B Jaynes
- Department of Biochemistry and Molecular Biology, Thomas Jefferson University, Philadelphia, Pennsylvania, USA
| | - Charles A Ettensohn
- Department of Biological Sciences, Carnegie Mellon University, Pittsburgh, Pennsylvania, USA.
| |
Collapse
|
20
|
Mulvey B, Lagunas T, Dougherty JD. Massively Parallel Reporter Assays: Defining Functional Psychiatric Genetic Variants Across Biological Contexts. Biol Psychiatry 2021; 89:76-89. [PMID: 32843144 PMCID: PMC7938388 DOI: 10.1016/j.biopsych.2020.06.011] [Citation(s) in RCA: 26] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/11/2020] [Revised: 06/09/2020] [Accepted: 06/10/2020] [Indexed: 12/18/2022]
Abstract
Neuropsychiatric phenotypes have long been known to be influenced by heritable risk factors, directly confirmed by the past decade of genetic studies that have revealed specific genetic variants enriched in disease cohorts. However, the initial hope that a small set of genes would be responsible for a given disorder proved false. The more complex reality is that a given disorder may be influenced by myriad small-effect noncoding variants and/or by rare but severe coding variants, many de novo. Noncoding genomic sequences-for which molecular functions cannot usually be inferred-harbor a large portion of these variants, creating a substantial barrier to understanding higher-order molecular and biological systems of disease. Fortunately, novel genetic technologies-scalable oligonucleotide synthesis, RNA sequencing, and CRISPR (clustered regularly interspaced short palindromic repeats)-have opened novel avenues to experimentally identify biologically significant variants en masse. Massively parallel reporter assays (MPRAs) are an especially versatile technique resulting from such innovations. MPRAs are powerful molecular genetics tools that can be used to screen thousands of untranscribed or untranslated sequences and their variants for functional effects in a single experiment. This approach, though underutilized in psychiatric genetics, has several useful features for the field. We review methods for assaying putatively functional genetic variants and regions, emphasizing MPRAs and the opportunities they hold for dissection of psychiatric polygenicity. We discuss literature applying functional assays in neurogenetics, highlighting strengths, caveats, and design considerations-especially regarding disease-relevant variables (cell type, neurodevelopment, and sex), and we ultimately propose applications of MPRA to both computational and experimental neurogenetics of polygenic disease risk.
Collapse
Affiliation(s)
- Bernard Mulvey
- Division of Biology and Biomedical Sciences, Washington University School of Medicine in St. Louis, St. Louis, Missouri; Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri; Department of Psychiatry, Washington University School of Medicine in St. Louis, St. Louis, Missouri
| | - Tomás Lagunas
- Division of Biology and Biomedical Sciences, Washington University School of Medicine in St. Louis, St. Louis, Missouri; Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri; Department of Psychiatry, Washington University School of Medicine in St. Louis, St. Louis, Missouri
| | - Joseph D Dougherty
- Department of Genetics, Washington University School of Medicine in St. Louis, St. Louis, Missouri; Department of Psychiatry, Washington University School of Medicine in St. Louis, St. Louis, Missouri.
| |
Collapse
|
21
|
High-Throughput Analysis of Retinal Cis-Regulatory Networks by Massively Parallel Reporter Assays. ADVANCES IN EXPERIMENTAL MEDICINE AND BIOLOGY 2020; 1185:359-364. [PMID: 31884638 DOI: 10.1007/978-3-030-27378-1_59] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
Abstract
Inherited retinal degenerations are diverse and debilitating blinding diseases. Genetic tests and exome sequencing have identified mutations in many protein-coding genes associated with such diseases, but causal sequence variants remain to be found in many retinopathy cases. Since 99% of our genome does not code for protein but contains cis-regulatory elements (CREs) that regulate the expression of essential genes, CRE variants might hold the answer for some of these cases. However, identifying functional CREs within the noncoding genome and predicting the pathogenicity of CRE variants pose a significant challenge. Here, we review the development of massively parallel reporter assays in the mouse retina, its use in dissecting retinal cis-regulatory networks, and its potential application for developing therapies.
Collapse
|
22
|
Perez-Cervantes C, Smith LA, Nadadur RD, Hughes AEO, Wang S, Corbo JC, Cepko C, Lonfat N, Moskowitz IP. Enhancer transcription identifies cis-regulatory elements for photoreceptor cell types. Development 2020; 147:dev184432. [PMID: 31915147 PMCID: PMC7033740 DOI: 10.1242/dev.184432] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2019] [Accepted: 12/13/2019] [Indexed: 12/30/2022]
Abstract
Identification of cell type-specific cis-regulatory elements (CREs) is crucial for understanding development and disease, although identification of functional regulatory elements remains challenging. We hypothesized that context-specific CREs could be identified by context-specific non-coding RNA (ncRNA) profiling, based on the observation that active CREs produce ncRNAs. We applied ncRNA profiling to identify rod and cone photoreceptor CREs from wild-type and mutant mouse retinas, defined by presence or absence, respectively, of the rod-specific transcription factor (TF) NrlNrl-dependent ncRNA expression strongly correlated with epigenetic profiles of rod and cone photoreceptors, identified thousands of candidate rod- and cone-specific CREs, and identified motifs for rod- and cone-specific TFs. Colocalization of NRL and the retinal TF CRX correlated with rod-specific ncRNA expression, whereas CRX alone favored cone-specific ncRNA expression, providing quantitative evidence that heterotypic TF interactions distinguish cell type-specific CRE activity. We validated the activity of novel Nrl-dependent ncRNA-defined CREs in developing cones. This work supports differential ncRNA profiling as a platform for the identification of cell type-specific CREs and the discovery of molecular mechanisms underlying TF-dependent CRE activity.
Collapse
Affiliation(s)
- Carlos Perez-Cervantes
- Departments of Pediatrics, Pathology, and Human Genetics, University of Chicago, Chicago, IL 60637, USA
| | - Linsin A Smith
- Departments of Pediatrics, Pathology, and Human Genetics, University of Chicago, Chicago, IL 60637, USA
| | - Rangarajan D Nadadur
- Departments of Pediatrics, Pathology, and Human Genetics, University of Chicago, Chicago, IL 60637, USA
| | - Andrew E O Hughes
- Department of Pathology and Immunology, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Sui Wang
- Departments of Genetics and Ophthalmology, Howard Hughes Medical Institute, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Joseph C Corbo
- Department of Pathology and Immunology, Washington University School of Medicine, St. Louis, MO 63110, USA
| | - Constance Cepko
- Departments of Genetics and Ophthalmology, Howard Hughes Medical Institute, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Nicolas Lonfat
- Departments of Genetics and Ophthalmology, Howard Hughes Medical Institute, Blavatnik Institute, Harvard Medical School, Boston, MA 02115, USA
| | - Ivan P Moskowitz
- Departments of Pediatrics, Pathology, and Human Genetics, University of Chicago, Chicago, IL 60637, USA
| |
Collapse
|
23
|
Zhou W, Dorrity MW, Bubb KL, Queitsch C, Fields S. Binding and Regulation of Transcription by Yeast Ste12 Variants To Drive Mating and Invasion Phenotypes. Genetics 2020; 214:397-407. [PMID: 31810988 PMCID: PMC7017024 DOI: 10.1534/genetics.119.302929] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/27/2019] [Accepted: 11/25/2019] [Indexed: 12/31/2022] Open
Abstract
Amino acid substitutions are commonly found in human transcription factors, yet the functional consequences of much of this variation remain unknown, even in well-characterized DNA-binding domains. Here, we examine how six single-amino acid variants in the DNA-binding domain of Ste12-a yeast transcription factor regulating mating and invasion-alter Ste12 genome binding, motif recognition, and gene expression to yield markedly different phenotypes. Using a combination of the "calling-card" method, RNA sequencing, and HT-SELEX (high throughput systematic evolution of ligands by exponential enrichment), we find that variants with dissimilar binding and expression profiles can converge onto similar cellular behaviors. Mating-defective variants led to decreased expression of distinct subsets of genes necessary for mating. Hyper-invasive variants also decreased expression of subsets of genes involved in mating, but increased the expression of other subsets of genes associated with the cellular response to osmotic stress. While single-amino acid changes in the coding region of this transcription factor result in complex regulatory reconfiguration, the major phenotypic consequences for the cell appear to depend on changes in the expression of a small number of genes with related functions.
Collapse
Affiliation(s)
- Wei Zhou
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195
- Molecular and Cellular Biology Program, University of Washington, Seattle, Washington 98195
| | - Michael W Dorrity
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195
| | - Kerry L Bubb
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195
| | - Christine Queitsch
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195
| | - Stanley Fields
- Department of Genome Sciences, University of Washington, Seattle, Washington 98195
- Department of Medicine, University of Washington, Seattle, Washington 98195
| |
Collapse
|
24
|
Murphy DP, Hughes AEO, Lawrence KA, Myers CA, Corbo JC. Cis-regulatory basis of sister cell type divergence in the vertebrate retina. eLife 2019; 8:e48216. [PMID: 31633482 PMCID: PMC6802965 DOI: 10.7554/elife.48216] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2019] [Accepted: 09/19/2019] [Indexed: 12/27/2022] Open
Abstract
Multicellular organisms evolved via repeated functional divergence of transcriptionally related sister cell types, but the mechanisms underlying sister cell type divergence are not well understood. Here, we study a canonical pair of sister cell types, retinal photoreceptors and bipolar cells, to identify the key cis-regulatory features that distinguish them. By comparing open chromatin maps and transcriptomic profiles, we found that while photoreceptor and bipolar cells have divergent transcriptomes, they share remarkably similar cis-regulatory grammars, marked by enrichment of K50 homeodomain binding sites. However, cell class-specific enhancers are distinguished by enrichment of E-box motifs in bipolar cells, and Q50 homeodomain motifs in photoreceptors. We show that converting K50 motifs to Q50 motifs represses reporter expression in bipolar cells, while photoreceptor expression is maintained. These findings suggest that partitioning of Q50 motifs within cell type-specific cis-regulatory elements was a critical step in the evolutionary divergence of the bipolar transcriptome from that of photoreceptors.
Collapse
Affiliation(s)
- Daniel P Murphy
- Department of Pathology and ImmunologyWashington University School of MedicineSt. LouisUnited States
| | - Andrew EO Hughes
- Department of Pathology and ImmunologyWashington University School of MedicineSt. LouisUnited States
| | - Karen A Lawrence
- Department of Pathology and ImmunologyWashington University School of MedicineSt. LouisUnited States
| | - Connie A Myers
- Department of Pathology and ImmunologyWashington University School of MedicineSt. LouisUnited States
| | - Joseph C Corbo
- Department of Pathology and ImmunologyWashington University School of MedicineSt. LouisUnited States
| |
Collapse
|