ABSTRACT
Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.
Subject(s)
Enhancer Elements, Genetic/genetics , High-Throughput Screening Assays/methods , Polydactyly/genetics , Animals , Enhancer Elements, Genetic/physiology , Gene Expression Regulation, Developmental/genetics , Gene Knock-In Techniques/methods , Hedgehog Proteins/genetics , Hedgehog Proteins/metabolism , Humans , Mice , Mutation , Phenotype , Polydactyly/metabolism , RNA, Untranslated/geneticsABSTRACT
The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP.
Subject(s)
Biological Evolution , Enhancer Elements, Genetic , Extremities/growth & development , Hedgehog Proteins/genetics , Snakes/genetics , Animals , Base Sequence , Evolution, Molecular , Gene Knock-In Techniques , Mice , Mice, Transgenic , Mutation , Phylogeny , Snakes/classificationABSTRACT
Uncovering the cis-regulatory code that governs when and how much each gene is transcribed in a given genome and cellular state remains a central goal of biology. Here, we discuss major layers of regulation that influence how transcriptional outputs are encoded by DNA sequence and cellular context. We first discuss how transcription factors bind specific DNA sequences in a dosage-dependent and cooperative manner and then proceed to the cofactors that facilitate transcription factor function and mediate the activity of modular cis-regulatory elements such as enhancers, silencers, and promoters. We then consider the complex and poorly understood interplay of these diverse elements within regulatory landscapes and its relationships with chromatin states and nuclear organization. We propose that a mechanistically informed, quantitative model of transcriptional regulation that integrates these multiple regulatory layers will be the key to ultimately cracking the cis-regulatory code.
Subject(s)
Enhancer Elements, Genetic , Transcription Factors , Transcription Factors/genetics , Transcription Factors/metabolism , Promoter Regions, Genetic , Gene Expression Regulation , Base Sequence , Chromatin/geneticsABSTRACT
Gene expression is in part controlled by cis-regulatory elements (CREs) such as enhancers and repressive elements. Anecdotal evidence has indicated that a CRE and a promoter need to be biochemically compatible for promoter regulation to occur, but this compatibility has remained poorly characterized in mammalian cells. We used high-throughput combinatorial reporter assays to test thousands of CRE-promoter pairs from three Mb-sized genomic regions in mouse cells. This revealed that CREs vary substantially in their promoter compatibility, ranging from striking specificity to broad promiscuity. More than half of the tested CREs exhibit significant promoter selectivity. Housekeeping promoters tend to have similar CRE preferences, but other promoters exhibit a wide diversity of compatibilities. Higher-order transcription factors (TF) motif combinations may account for compatibility. CRE-promoter selectivity does not correlate with looping interactions in the native genomic context, suggesting that chromatin folding and compatibility are two orthogonal mechanisms that confer specificity to gene regulation.
Subject(s)
Enhancer Elements, Genetic , Genome , Promoter Regions, Genetic , Transcription Factors , Animals , Enhancer Elements, Genetic/genetics , Gene Expression Regulation , Genome/genetics , Genomics , Mammals/metabolism , Mice , Promoter Regions, Genetic/genetics , Regulatory Sequences, Nucleic Acid/genetics , Transcription Factors/genetics , Transcription Factors/metabolismABSTRACT
Development and function of conventional dendritic cell (cDC) subsets, cDC1 and cDC2, depend on transcription factors (TFs) IRF8 and IRF4, respectively. Since IRF8 and IRF4 can each interact with TF BATF3 at AP1-IRF composite elements (AICEs) and with TF PU.1 at Ets-IRF composite elements (EICEs), it is unclear how these factors exert divergent actions. Here, we determined the basis for distinct effects of IRF8 and IRF4 in cDC development. Genes expressed commonly by cDC1 and cDC2 used EICE-dependent enhancers that were redundantly activated by low amounts of either IRF4 or IRF8. By contrast, cDC1-specific genes relied on AICE-dependent enhancers, which required high IRF concentrations, but were activated by either IRF4 or IRF8. IRF8 was specifically required only by a minority of cDC1-specific genes, such as Xcr1, which could distinguish between IRF8 and IRF4 DNA-binding domains. Thus, these results explain how BATF3-dependent Irf8 autoactivation underlies emergence of the cDC1-specific transcriptional program.
Subject(s)
Dendritic Cells/metabolism , Enhancer Elements, Genetic/genetics , Interferon Regulatory Factors/genetics , Animals , Gene Expression Regulation/genetics , Mice , Mice, Inbred C57BL , Receptors, Chemokine/genetics , Transcription, Genetic/geneticsABSTRACT
To withstand a hostile cellular environment and replicate, viruses must sense, interpret, and respond to many internal and external cues. Retroviruses and DNA viruses can intercept these cues impinging on host transcription factors via cis-regulatory elements (CREs) in viral genomes, allowing them to sense and coordinate context-specific responses to varied signals. Here, we explore the characteristics of viral CREs, the classes of signals and host transcription factors that regulate them, and how this informs outcomes of viral replication, immune evasion, and latency. We propose that viral CREs constitute central hubs for signal integration from multiple pathways and that sequence variation between viral isolates can rapidly rewire sensing mechanisms, contributing to the variability observed in patient outcomes.
ABSTRACT
Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.
Subject(s)
Alzheimer Disease , Animals , Humans , Alzheimer Disease/genetics , Evolution, Molecular , Primates/genetics , BrainABSTRACT
Cellular lineage determination is controlled by combinations of lineage-selective transcription factors (TFs) and associated coregulators that bind to cis-regulatory elements in DNA and regulate gene expression. The ability of these factors to regulate transcription is determined not only by their cooperativity, but also by biochemical and structural properties of the chromatin, sculpting higher-order genome organization. Here, we review recent advances in the understanding of the interplay between chromatin topology and transcription. Studies from many different fields, including adipocyte lineage determination, indicate that lineage determination and differentiation are dependent on elaborate crosstalk between cis-regulatory elements, leading to the formation of transcriptional hubs. Chromatin topology appears to provide a dynamic and supportive, rather than a deterministic, scaffold for this crosstalk.
Subject(s)
Chromatin , Enhancer Elements, Genetic , Cell Lineage/genetics , Chromatin/genetics , DNA , Transcription Factors/genetics , Transcription Factors/metabolismABSTRACT
BACKGROUND: The endothelial-to-hematopoietic transition (EHT) process during definitive hematopoiesis is highly conserved in vertebrates. Stage-specific expression of transposable elements (TEs) has been detected during zebrafish EHT and may promote hematopoietic stem cell (HSC) formation by activating inflammatory signaling. However, little is known about how TEs contribute to the EHT process in human and mouse. RESULTS: We reconstructed the single-cell EHT trajectories of human and mouse and resolved the dynamic expression patterns of TEs during EHT. Most TEs presented a transient co-upregulation pattern along the conserved EHT trajectories, coinciding with the temporal relaxation of epigenetic silencing systems. TE products can be sensed by multiple pattern recognition receptors, triggering inflammatory signaling to facilitate HSC emergence. Interestingly, we observed that hypoxia-related signals were enriched in cells with higher TE expression. Furthermore, we constructed the hematopoietic cis-regulatory network of accessible TEs and identified potential TE-derived enhancers that may boost the expression of specific EHT marker genes. CONCLUSIONS: Our study provides a systematic vision of how TEs are dynamically controlled to promote the hematopoietic fate decisions through transcriptional and cis-regulatory networks, and pre-train the immunity of nascent HSCs.
Subject(s)
DNA Transposable Elements , Hematopoiesis , Hematopoietic Stem Cells , Single-Cell Analysis , Animals , DNA Transposable Elements/genetics , Single-Cell Analysis/methods , Mice , Hematopoiesis/genetics , Humans , Hematopoietic Stem Cells/metabolism , Endothelial Cells/metabolismABSTRACT
Brain development and function are governed by precisely regulated protein expressions in different regions. To date, multiregional brain proteomes have been systematically analyzed only for adult human and mouse brains. To understand the underpinnings of brain development and function, we generated proteomes from six regions of the postnatal brain at three developmental stages of domestic dogs (Canis familiaris), which are special among animals in terms of their remarkable human-like social cognitive abilities. Quantitative analysis of the spatiotemporal proteomes identified region-enriched synapse types at different developmental stages and differential myelination progression in different brain regions. Through integrative analysis of inter-regional expression patterns of orthologous proteins and genome-wide cis-regulatory element frequencies, we found that proteins related with myelination and hippocampus were highly correlated between dog and human but not between mouse and human, although mouse is phylogenetically closer to human. Moreover, the global expression patterns of neurodegenerative disease and autism spectrum disorder-associated proteins in dog brain more resemble human brain than in mouse brain. The high similarity of myelination and hippocampus-related pathways in dog and human at both proteomic and genetic levels may contribute to their shared social cognitive abilities. The inter-regional expression patterns of disease-associated proteins in the brain of different species provide important information to guide mechanistic and translational study using appropriate animal models.
Subject(s)
Autism Spectrum Disorder , Neurodegenerative Diseases , Adult , Animals , Brain , Dogs , Humans , Mice , Proteome , ProteomicsABSTRACT
Genome-wide association studies have identified the chromosome 10q26 (Chr10) locus, which contains the age-related maculopathy susceptibility 2 (ARMS2) and high temperature requirement A serine peptidase 1 (HTRA1) genes, as the strongest genetic risk factor for age-related macular degeneration (AMD) [L.G. Fritsche et al., Annu. Rev. Genomics Hum. Genet. 15, 151-171, (2014)]. To date, it has been difficult to assign causality to any specific single nucleotide polymorphism (SNP), haplotype, or gene within this region because of high linkage disequilibrium among the disease-associated variants [J. Jakobsdottir et al. Am. J. Hum. Genet. 77, 389-407 (2005); A. Rivera et al. Hum. Mol. Genet. 14, 3227-3236 (2005)]. Here, we show that HTRA1 messenger RNA (mRNA) is reduced in retinal pigment epithelium (RPE) but not in neural retina or choroid tissues derived from human donors with homozygous risk at the 10q26 locus. This tissue-specific decrease is mediated by the presence of a noncoding, cis-regulatory element overlapping the ARMS2 intron, which contains a potential Lhx2 transcription factor binding site that is disrupted by risk variant rs36212733. HtrA1 protein increases with age in the RPE-Bruch's membrane (BM) interface in Chr10 nonrisk donors but fails to increase in donors with homozygous risk at the 10q26 locus. We propose that HtrA1, an extracellular chaperone and serine protease, functions to maintain the optimal integrity of the RPE-BM interface during the aging process and that reduced expression of HTRA1 mRNA and protein in Chr10 risk donors impairs this protective function, leading to increased risk of AMD pathogenesis. HtrA1 augmentation, not inhibition, in high-risk patients should be considered as a potential therapy for AMD.
Subject(s)
Genetic Predisposition to Disease , High-Temperature Requirement A Serine Peptidase 1/metabolism , Macular Degeneration/genetics , Retinal Pigment Epithelium/metabolism , Choroid/metabolism , Genetic Variation , High-Temperature Requirement A Serine Peptidase 1/genetics , Humans , Linkage Disequilibrium , RNA, Messenger/genetics , RNA, Messenger/metabolism , Retina/metabolismABSTRACT
3'-untranslated regions (UTRs) specify post-transcriptional fates of mammalian messenger RNAs (mRNAs), yet knowledge of the underlying sequences and mechanisms is largely incomplete. Here, we identify two related novel 3' UTR motifs in mammals that specify transcript degradation. These motifs are interchangeable and active only within 3' UTRs, where they are often preferentially conserved; furthermore, they are found in hundreds of transcripts, many encoding regulatory proteins. We found that degradation occurs via mRNA deadenylation, mediated by the CCR4-NOT complex. We purified trans factors that recognize the motifs and identified heterogeneous nuclear ribonucleoproteins (hnRNPs) A1 and A2/B1, which are required for transcript degradation, acting in a previously unknown manner. We used RNA sequencing (RNA-seq) to confirm hnRNP A1 and A2/B1 motif-dependent roles genome-wide, profiling cells depleted of these factors singly and in combination. Interestingly, the motifs are most active within the distal portion of 3' UTRs, suggesting that their role in gene regulation can be modulated by alternative processing, resulting in shorter 3' UTRs.
Subject(s)
Gene Expression Regulation/genetics , Heterogeneous-Nuclear Ribonucleoprotein Group A-B/metabolism , RNA Stability/genetics , 3' Untranslated Regions/genetics , A549 Cells , Amino Acid Motifs/genetics , Animals , COS Cells , Cell Line , Chlorocebus aethiops , HEK293 Cells , Heterogeneous Nuclear Ribonucleoprotein A1 , Humans , MCF-7 Cells , Mice , Regulatory Elements, Transcriptional/genetics , TranscriptomeABSTRACT
BACKGROUND: The auxin indole-3-acetic acid (IAA) is a vital phytohormone that influences plant growth and development. Our previous work showed that IAA content decreased during flower development in the medicinally important orchid Dendrobium officinale, while Aux/IAA genes were downregulated. However, little information about auxin-responsive genes and their roles in D. officinale flower development exists. RESULTS: This study validated 14 DoIAA and 26 DoARF early auxin-responsive genes in the D. officinale genome. A phylogenetic analysis classified the DoIAA genes into two subgroups. An analysis of cis-regulatory elements indicated that they were related by phytohormones and abiotic stresses. Gene expression profiles were tissue-specific. Most DoIAA genes (except for DoIAA7) were sensitive to IAA (10 µmol/L) and were downregulated during flower development. Four DoIAA proteins (DoIAA1, DoIAA6, DoIAA10 and DoIAA13) were mainly localized in the nucleus. A yeast two-hybrid assay showed that these four DoIAA proteins interacted with three DoARF proteins (DoARF2, DoARF17, DoARF23). CONCLUSIONS: The structure and molecular functions of early auxin-responsive genes in D. officinale were investigated. The DoIAA-DoARF interaction may play an important role in flower development via the auxin signaling pathway.
Subject(s)
Dendrobium , Dendrobium/genetics , Dendrobium/metabolism , Phylogeny , Plant Proteins/metabolism , Indoleacetic Acids/metabolism , Plant Growth Regulators/metabolism , Flowers/genetics , Flowers/metabolism , Gene Expression Regulation, PlantABSTRACT
BACKGROUND: Triple-negative breast cancer (TNBC) is an aggressive subtype that exhibits a high incidence of distant metastases and lacks targeted therapeutic options. Here we explored how the epigenome contributes to matrix metalloprotease (MMP) dysregulation impacting tumor invasion, which is the first step of the metastatic process. METHODS: We combined RNA expression and chromatin interaction data to identify insulator elements potentially associated with MMP gene expression and invasion. We employed CRISPR/Cas9 to disrupt the CCCTC-Binding Factor (CTCF) binding site on an insulator element downstream of the MMP8 gene (IE8) in two TNBC cellular models. We characterized these models by combining Hi-C, ATAC-seq, and RNA-seq with functional experiments to determine invasive ability. The potential of our findings to predict the progression of ductal carcinoma in situ (DCIS), was tested in data from clinical specimens. RESULTS: We explored the clinical relevance of an insulator element located within the Chr11q22.2 locus, downstream of the MMP8 gene (IE8). This regulatory element resulted in a topologically associating domain (TAD) boundary that isolated nine MMP genes into two anti-correlated expression clusters. This expression pattern was associated with worse relapse-free (HR = 1.57 [1.06 - 2.33]; p = 0.023) and overall (HR = 2.65 [1.31 - 5.37], p = 0.005) survival of TNBC patients. After CRISPR/Cas9-mediated disruption of IE8, cancer cells showed a switch in the MMP expression signature, specifically downregulating the pro-invasive MMP1 gene and upregulating the antitumorigenic MMP8 gene, resulting in reduced invasive ability and collagen degradation. We observed that the MMP expression pattern predicts DCIS that eventually progresses into invasive ductal carcinomas (AUC = 0.77, p < 0.01). CONCLUSION: Our study demonstrates how the activation of an IE near the MMP8 gene determines the regional transcriptional regulation of MMP genes with opposing functional activity, ultimately influencing the invasive properties of aggressive forms of breast cancer.
Subject(s)
Breast Neoplasms , Carcinoma, Intraductal, Noninfiltrating , Triple Negative Breast Neoplasms , Humans , Female , Breast Neoplasms/genetics , Breast Neoplasms/pathology , Carcinoma, Intraductal, Noninfiltrating/genetics , Carcinoma, Intraductal, Noninfiltrating/pathology , Chromatin , Matrix Metalloproteinase 8/genetics , Triple Negative Breast Neoplasms/genetics , Neoplasm Recurrence, Local/genetics , Multigene FamilyABSTRACT
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.
Subject(s)
Basic-Leucine Zipper Transcription Factors/metabolism , Eye Proteins/metabolism , Regulatory Sequences, Nucleic Acid/genetics , Retinal Cone Photoreceptor Cells/metabolism , Retinal Rod Photoreceptor Cells/metabolism , Transcription, Genetic/genetics , Animals , Basic-Leucine Zipper Transcription Factors/genetics , Eye Proteins/genetics , Female , Gene Expression Regulation, Developmental , Homeodomain Proteins/genetics , Homeodomain Proteins/metabolism , Male , Mice , Mice, Knockout , RNA, Untranslated/genetics , RNA, Untranslated/metabolism , Trans-Activators/genetics , Trans-Activators/metabolism , TranscriptomeABSTRACT
Changes in gene expression are a prominent feature of morphological evolution. These changes occur to hierarchical gene regulatory networks (GRNs) of transcription factor genes that regulate the expression of trait-building differentiation genes. While changes in the expression of differentiation genes are essential to phenotypic evolution, they can be caused by mutations within cis-regulatory elements (CREs) that drive their expression (cis-evolution) or within genes for CRE-interacting transcription factors (trans-evolution). Locating these mutations remains a challenge, especially when experiments are limited to one species that possesses the ancestral or derived phenotype. We investigated CREs that control the expression of the differentiation genes tan and yellow, the expression of which evolved during the gain, modification, and loss of dimorphic pigmentation among Sophophora fruit flies. We show these CREs to be necessary components of a pigmentation GRN, as deletion from Drosophila melanogaster (derived dimorphic phenotype) resulted in lost expression and lost male-specific pigmentation. We evaluated the ability of orthologous CRE sequences to drive reporter gene expression in species with modified (Drosophila auraria), secondarily lost (Drosophila ananassae), and ancestrally absent (Drosophila willistoni) pigmentation. We show that the transgene host frequently determines CRE activity, implicating trans-evolution as a significant factor for this trait's diversity. We validated the gain of dimorphic Bab transcription factor expression as a trans-change contributing to the dimorphic trait. Our findings suggest an amenability to change for the landscape of trans-regulators and begs for an explanation as to why this is so common compared to the evolution of differentiation gene CREs.
Subject(s)
Drosophila Proteins , Drosophila melanogaster , Male , Animals , Drosophila melanogaster/genetics , Drosophila Proteins/genetics , Drosophila/genetics , Drosophila/metabolism , Transcription Factors/genetics , Pigmentation/genetics , Phenotype , Evolution, MolecularABSTRACT
Organismal phenotypes result largely from inherited developmental programs, usually executed during embryonic and juvenile life stages. These programs are not blank slates onto which natural selection can draw arbitrary forms. Rather, the mechanisms of development play an integral role in shaping phenotypic diversity and help determine the evolutionary trajectories of species. Modern evolutionary biology must, therefore, account for these mechanisms in both theory and in practice. The gene regulatory network (GRN) concept represents a potent tool for achieving this goal whose utility has grown in tandem with advances in "omic" technologies and experimental techniques. However, while the GRN concept is widely utilized, it is often less clear what practical implications it has for conducting research in evolutionary developmental biology. In this Perspective, we attempt to provide clarity by discussing how experiments and projects can be designed in light of the GRN concept. We first map familiar biological notions onto the more abstract components of GRN models. We then review how diverse functional genomic approaches can be directed toward the goal of constructing such models and discuss current methods for functionally testing evolutionary hypotheses that arise from them. Finally, we show how the major steps of GRN model construction and experimental validation suggest generalizable workflows that can serve as a scaffold for project design. Taken together, the practical implications that we draw from the GRN concept provide a set of guideposts for studies aiming at unraveling the molecular basis of phenotypic diversity.
Subject(s)
Biological Evolution , Gene Regulatory Networks , Animals , Phenotype , Genomics , Developmental BiologyABSTRACT
The CRISPR-Cas-based genome editing field in plants is expanding rapidly. Editing plant promoters to obtain cis-regulatory alleles with altered expression levels or patterns of target genes is a highly promising topic. However, primarily used CRISPR-Cas9 has significant limitations when editing noncoding sequences like promoters, which have unique structures and regulatory mechanisms, including A-T richness, repetitive redundancy, difficulty in identifying key regulatory regions, and a higher frequency of DNA structure, epigenetic modification, and protein binding accessibility issues. Researchers urgently require efficient and feasible editing tools and strategies to address these obstacles, enhance promoter editing efficiency, increase diversity in promoter polymorphism, and, most importantly, enable 'non-silent' editing events that achieve precise target gene expression regulation. This article provides insights into the key challenges and references for implementing promoter editing-based research in plants.
Subject(s)
CRISPR-Cas Systems , Gene Editing , CRISPR-Cas Systems/genetics , Plants/genetics , Promoter Regions, Genetic/genetics , Regulatory Sequences, Nucleic Acid , Genome, PlantABSTRACT
The interplay of transcription factors and cis-regulatory elements (CREs) orchestrates the dynamic and diverse genetic programs that assemble the human central nervous system (CNS) during development and maintain its function throughout life. Genetic variation within CREs plays a central role in phenotypic variation in complex traits including the risk of developing disease. We took advantage of the retina, a well-characterized region of the CNS known to be affected by pathogenic variants in CREs, to establish a roadmap for characterizing regulatory variation in the human CNS. This comprehensive analysis of tissue-specific regulatory elements, transcription factor binding, and gene expression programs in three regions of the human visual system (retina, macula, and retinal pigment epithelium/choroid) reveals features of regulatory element evolution that shape tissue-specific gene expression programs and defines regulatory elements with the potential to contribute to Mendelian and complex disorders of human vision.
Subject(s)
Evolution, Molecular , Gene Expression Regulation, Developmental , Regulatory Sequences, Nucleic Acid/genetics , Retina/pathology , Retinal Diseases/genetics , Adult , Animals , DNA Mutational Analysis , Epigenomics , Female , Genetic Variation , Humans , Male , Mice , Middle Aged , Mutation , RNA-Seq , Retina/growth & development , Retinal Diseases/pathology , Species SpecificityABSTRACT
Co-option of transposable elements (TEs) to become part of existing or new enhancers is an important mechanism for evolution of gene regulation. However, contributions of lineage-specific TE insertions to recent regulatory adaptations remain poorly understood. Gibbons present a suitable model to study these contributions as they have evolved a lineage-specific TE called LAVA (LINE-AluSz-VNTR-AluLIKE), which is still active in the gibbon genome. The LAVA retrotransposon is thought to have played a role in the emergence of the highly rearranged structure of the gibbon genome by disrupting transcription of cell cycle genes. In this study, we investigated whether LAVA may have also contributed to the evolution of gene regulation by adopting enhancer function. We characterized fixed and polymorphic LAVA insertions across multiple gibbons and found 96 LAVA elements overlapping enhancer chromatin states. Moreover, LAVA was enriched in multiple transcription factor binding motifs, was bound by an important transcription factor (PU.1), and was associated with higher levels of gene expression in cis We found gibbon-specific signatures of purifying/positive selection at 27 LAVA insertions. Two of these insertions were fixed in the gibbon lineage and overlapped with enhancer chromatin states, representing putative co-opted LAVA enhancers. These putative enhancers were located within genes encoding SETD2 and RAD9A, two proteins that facilitate accurate repair of DNA double-strand breaks and prevent chromosomal rearrangement mutations. Co-option of LAVA in these genes may have influenced regulation of processes that preserve genome integrity. Our findings highlight the importance of considering lineage-specific TEs in studying evolution of gene regulatory elements.