Enhancer transcription identifies cis-regulatory elements for photoreceptor cell types

Machine Learning Prediction of Non-Coding Variant Impact in Cell-Class-Specific Human Retinal Cis-Regulatory Elements

Non-coding variants in cis-regulatory elements such as promoters and enhancers contribute to inherited retinal diseases (IRDs), however, characterizing the functional impact of most regulatory variants remains challenging. To improve identification of variants of interest, we implemented machine learning using a gapped k-mer support vector machine approach trained on single nucleus ATAC-seq data from specific cell classes of the adult and developing human retina. We developed 18 distinct ML models to predict the impact of non-coding variants on 39,437 cell-class-specific regulatory elements. These models demonstrate accuracy over 90% and a high degree of cell class specificity. Variant Impact Prediction (VIP) scores highlight specific sequences within candidate CREs, including putative transcription factor (TF) binding motifs, that are predicted to alter CRE function if mutated. Correlations to massively parallel reporter assays support the predictive value of VIP scores to model single nucleotide variants and indels in a cell-class-specific manner. These analyses demonstrate the capacity for single nucleus epigenomic data to predict the impact of non-coding sequence variants and allow for rapid prioritization of patient variants for further functional analysis.

A Genomic Link From Heart Failure to Atrial Fibrillation Risk: FOG2 Modulates a TBX5/GATA4-Dependent Atrial Gene Regulatory Network

BACKGROUND

The relationship between heart failure (HF) and atrial fibrillation (AF) is clear, with up to half of patients with HF progressing to AF. The pathophysiological basis of AF in the context of HF is presumed to result from atrial remodeling. Upregulation of the transcription factor FOG2 (friend of GATA2; encoded by ZFPM2) is observed in human ventricles during HF and causes HF in mice.

METHODS

FOG2 expression was assessed in human atria. The effect of adult-specific FOG2 overexpression in the mouse heart was evaluated by whole animal electrophysiology, in vivo organ electrophysiology, cellular electrophysiology, calcium flux, mouse genetic interactions, gene expression, and genomic function, including a novel approach for defining functional transcription factor interactions based on overlapping effects on enhancer noncoding transcription.

RESULTS

FOG2 is significantly upregulated in the human atria during HF. Adult cardiomyocyte-specific FOG2 overexpression in mice caused primary spontaneous AF before the development of HF or atrial remodeling. FOG2 overexpression generated arrhythmia substrate and trigger in cardiomyocytes, including calcium cycling defects. We found that FOG2 repressed atrial gene expression promoted by TBX5. FOG2 bound a subset of GATA4 and TBX5 co-bound genomic locations, defining a shared atrial gene regulatory network. FOG2 repressed TBX5-dependent transcription from a subset of co-bound enhancers, including a conserved enhancer at the Atp2a2 locus. Atrial rhythm abnormalities in mice caused by Tbx5 haploinsufficiency were rescued by Zfpm2 haploinsufficiency.

CONCLUSIONS

Transcriptional changes in the atria observed in human HF directly antagonize the atrial rhythm gene regulatory network, providing a genomic link between HF and AF risk independent of atrial remodeling.

Key transcription factors influence the epigenetic landscape to regulate retinal cell differentiation

How the diverse neural cell types emerge from multipotent neural progenitor cells during central nervous system development remains poorly understood. Recent scRNA-seq studies have delineated the developmental trajectories of individual neural cell types in many neural systems including the neural retina. Further understanding of the formation of neural cell diversity requires knowledge about how the epigenetic landscape shifts along individual cell lineages and how key transcription factors regulate these changes. In this study, we dissect the changes in the epigenetic landscape during early retinal cell differentiation by scATAC-seq and identify globally the enhancers, enriched motifs, and potential interacting transcription factors underlying the cell state/type specific gene expression in individual lineages. Using CUT&Tag, we further identify the enhancers bound directly by four key transcription factors, Otx2, Atoh7, Pou4f2 and Isl1, including those dependent on Atoh7, and uncover the sequential and combinatorial interactions of these factors with the epigenetic landscape to control gene expression along individual retinal cell lineages such as retinal ganglion cells (RGCs). Our results reveal a general paradigm in which transcription factors collaborate and compete to regulate the emergence of distinct retinal cell types such as RGCs from multipotent retinal progenitor cells (RPCs).

Identification of the Major Effector StSROs in Potato: A Potential StWRKY-SRO6 Regulatory Pathway Enhances Plant Tolerance to Cadmium Stress

SIMILAR TO RCD-ONE (SRO) family members and transcription factors (TFs) often improve plant antioxidant capacity through interaction and co-regulation and participate in plant resistance to drought and high-salt stress. However, whether SROs are involved in the response to heavy metal stress, especially SRO genes with a specific response and tolerance characteristics to cadmium (Cd) stress, remains unclear. We first identified six SRO genes in the potato genome by PARP and RST domains. Special and conserved StSROs were found, and the spatio temporal tissue-specific expression patterns and co-expression network diagrams of StSROs under the stress of 5 heavy metals were constructed. Second, we identified StSRO6 as a major effector gene (StSRO6-MEG) and StSRO5 as a secondary effector gene (StSRO5-SEG) through a comprehensive analysis. Interestingly, they may hold true for various physiological or stress responses in plants. In addition, using systematic genomics and comparative omics techniques, the key gene StSRO6 that affects the difference in Cd accumulation was discovered, cloned in the low-Cd accumulation "Yunshu 505", and transformed into the yeast mutant ycf1 for overexpression. The results proved that StSRO6 could confer Cd tolerance. Finally, through transient expression and in vitro culture tests, we hypothesized that StSROs 5/6 are regulated by the transcription factor StWRKY6 and mediates the reactive oxygen species (ROS) system to confer Cd tolerance. These findings offer a new perspective for understanding the mechanisms underlying Cd tolerance in plants, and simultaneously provide clues for the development of biological agents for preventing and controlling Cd migration and transformation.

Machine Learning Prediction of Non-Coding Variant Impact in Human Retinal cis-Regulatory Elements

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.

Deciphering the Retinal Epigenome during Development, Disease and Reprogramming: Advancements, Challenges and Perspectives

Retinal neurogenesis is driven by concerted actions of transcription factors, some of which are expressed in a continuum and across several cell subtypes throughout development. While seemingly redundant, many factors diversify their regulatory outcome on gene expression, by coordinating variations in chromatin landscapes to drive divergent retinal specification programs. Recent studies have furthered the understanding of the epigenetic contribution to the progression of age-related macular degeneration, a leading cause of blindness in the elderly. The knowledge of the epigenomic mechanisms that control the acquisition and stabilization of retinal cell fates and are evoked upon damage, holds the potential for the treatment of retinal degeneration. Herein, this review presents the state-of-the-art approaches to investigate the retinal epigenome during development, disease, and reprogramming. A pipeline is then reviewed to functionally interrogate the epigenetic and transcriptional networks underlying cell fate specification, relying on a truly unbiased screening of open chromatin states. The related work proposes an inferential model to identify gene regulatory networks, features the first footprinting analysis and the first tentative, systematic query of candidate pioneer factors in the retina ever conducted in any model organism, leading to the identification of previously uncharacterized master regulators of retinal cell identity, such as the nuclear factor I, NFI. This pipeline is virtually applicable to the study of genetic programs and candidate pioneer factors in any developmental context. Finally, challenges and limitations intrinsic to the current next-generation sequencing techniques are discussed, as well as recent advances in super-resolution imaging, enabling spatio-temporal resolution of the genome.

Building a Mammalian Retina: An Eye on Chromatin Structure

Regulation of gene expression by chromatin structure has been under intensive investigation, establishing nuclear organization and genome architecture as a potent and effective means of regulating developmental processes. The substantial growth in our knowledge of the molecular mechanisms underlying retinogenesis has been powered by several genome-wide based tools that mapped chromatin organization at multiple cellular and biochemical levels. Studies profiling the retinal epigenome and transcriptome have allowed the systematic annotation of putative cis-regulatory elements associated with transcriptional programs that drive retinal neural differentiation, laying the groundwork to understand spatiotemporal retinal gene regulation at a mechanistic level. In this review, we outline recent advances in our understanding of the chromatin architecture in the mammalian retina during development and disease. We focus on the emerging roles of non-coding regulatory elements in controlling retinal cell-type specific transcriptional programs, and discuss potential implications in untangling the etiology of eye-related disorders.

Identification of direct transcriptional targets of NFATC2 that promote β cell proliferation

The transcription factor NFATC2 induces β cell proliferation in mouse and human islets. However, the genomic targets that mediate these effects have not been identified. We expressed active forms of Nfatc2 and Nfatc1 in human islets. By integrating changes in gene expression with genomic binding sites for NFATC2, we identified approximately 2200 transcriptional targets of NFATC2. Genes induced by NFATC2 were enriched for transcripts that regulate the cell cycle and for DNA motifs associated with the transcription factor FOXP. Islets from an endocrine-specific Foxp1, Foxp2, and Foxp4 triple-knockout mouse were less responsive to NFATC2-induced β cell proliferation, suggesting the FOXP family works to regulate β cell proliferation in concert with NFATC2. NFATC2 induced β cell proliferation in both mouse and human islets, whereas NFATC1 did so only in human islets. Exploiting this species difference, we identified approximately 250 direct transcriptional targets of NFAT in human islets. This gene set enriches for cell cycle-associated transcripts and includes Nr4a1. Deletion of Nr4a1 reduced the capacity of NFATC2 to induce β cell proliferation, suggesting that much of the effect of NFATC2 occurs through its induction of Nr4a1. Integration of noncoding RNA expression, chromatin accessibility, and NFATC2 binding sites enabled us to identify NFATC2-dependent enhancer loci that mediate β cell proliferation.

Initiation of Otx2 expression in the developing mouse retina requires a unique enhancer and either Ascl1 or Neurog2 activity

During retinal development, a large subset of progenitors upregulates the transcription factor Otx2, which is required for photoreceptor and bipolar cell formation. How these retinal progenitor cells initially activate Otx2 expression is unclear. To address this, we investigated the cis-regulatory network that controls Otx2 expression in mice. We identified a minimal enhancer element, DHS-4D, that drove expression in newly formed OTX2+ cells. CRISPR/Cas9-mediated deletion of DHS-4D reduced OTX2 expression, but this effect was diminished in postnatal development. Systematic mutagenesis of the enhancer revealed that three basic helix-loop-helix (bHLH) transcription factor-binding sites were required for its activity. Single cell RNA-sequencing of nascent Otx2+ cells identified the bHLH factors Ascl1 and Neurog2 as candidate regulators. CRISPR/Cas9 targeting of these factors showed that only the simultaneous loss of Ascl1 and Neurog2 prevented OTX2 expression. Our findings suggest that Ascl1 and Neurog2 act either redundantly or in a compensatory fashion to activate the DHS-4D enhancer and Otx2 expression. We observed redundancy or compensation at both the transcriptional and enhancer utilization levels, suggesting that the mechanisms governing Otx2 regulation in the retina are flexible and robust.

Cis-regulatory dissection of cone development reveals a broad role for Otx2 and Oc transcription factors

The vertebrate retina is generated by retinal progenitor cells (RPCs), which produce >100 cell types. Although some RPCs produce many cell types, other RPCs produce restricted types of daughter cells, such as a cone photoreceptor and a horizontal cell (HC). We used genome-wide assays of chromatin structure to compare the profiles of a restricted cone/HC RPC and those of other RPCs in chicks. These data nominated regions of regulatory activity, which were tested in tissue, leading to the identification of many cis-regulatory modules (CRMs) active in cone/HC RPCs and developing cones. Two transcription factors, Otx2 and Oc1, were found to bind to many of these CRMs, including those near genes important for cone development and function, and their binding sites were required for activity. We also found that Otx2 has a predicted autoregulatory CRM. These results suggest that Otx2, Oc1 and possibly other Onecut proteins have a broad role in coordinating cone development and function. The many newly discovered CRMs for cones are potentially useful reagents for gene therapy of cone diseases.

The Impact of Modern Technologies on Molecular Diagnostic Success Rates, with a Focus on Inherited Retinal Dystrophy and Hearing Loss

The identification of pathogenic variants in monogenic diseases has been of interest to researchers and clinicians for several decades. However, for inherited diseases with extremely high genetic heterogeneity, such as hearing loss and retinal dystrophies, establishing a molecular diagnosis requires an enormous effort. In this review, we use these two genetic conditions as examples to describe the initial molecular genetic identification approaches, as performed since the early 90s, and subsequent improvements and refinements introduced over the years. Next, the history of DNA sequencing from conventional Sanger sequencing to high-throughput massive parallel sequencing, a.k.a. next-generation sequencing, is outlined, including their advantages and limitations and their impact on identifying the remaining genetic defects. Moreover, the development of recent technologies, also coined "third-generation" sequencing, is reviewed, which holds the promise to overcome these limitations. Furthermore, we outline the importance and complexity of variant interpretation in clinical diagnostic settings concerning the massive number of different variants identified by these methods. Finally, we briefly mention the development of novel approaches such as optical mapping and multiomics, which can help to further identify genetic defects in the near future.

Enhancer grammar in development, evolution, and disease: dependencies and interplay

Each language has standard books describing that language's grammatical rules. Biologists have searched for similar, albeit more complex, principles relating enhancer sequence to gene expression. Here, we review the literature on enhancer grammar. We introduce dependency grammar, a model where enhancers encode information based on dependencies between enhancer features shaped by mechanistic, evolutionary, and biological constraints. Classifying enhancers based on the types of dependencies may identify unifying principles relating enhancer sequence to gene expression. Such rules would allow us to read the instructions for development within genomes and pinpoint causal enhancer variants underlying disease and evolutionary changes.

Cell type- and stage-specific expression of Otx2 is regulated by multiple transcription factors and cis-regulatory modules in the retina

Transcription factors (TFs) are often used repeatedly during development and homeostasis to control distinct processes in the same and/or different cellular contexts. Considering the limited number of TFs in the genome and the tremendous number of events that need to be regulated, re-use of TFs is necessary. We analyzed how the expression of the homeobox TF, orthodenticle homeobox 2 (Otx2), is regulated in a cell type- and stage-specific manner during development in the mouse retina. We identified seven Otx2 cis-regulatory modules (CRMs), among which the O5, O7 and O9 CRMs mark three distinct cellular contexts of Otx2 expression. We discovered that Otx2, Crx and Sox2, which are well-known TFs regulating retinal development, bind to and activate the O5, O7 or O9 CRMs, respectively. The chromatin status of these three CRMs was found to be distinct in vivo in different retinal cell types and at different stages. We conclude that retinal cells use a cohort of TFs with different expression patterns and multiple CRMs with different chromatin configurations to regulate the expression of Otx2 precisely.

Cis-regulatory analysis of Onecut1 expression in fate-restricted retinal progenitor cells

Background

The vertebrate retina consists of six major classes of neuronal cells. During development, these cells are generated from a pool of multipotent retinal progenitor cells (RPCs) that express the gene Vsx2. Fate-restricted RPCs have recently been identified, with limited mitotic potential and cell fate possibilities compared to multipotent RPCs. One population of fate-restricted RPCs, marked by activity of the regulatory element ThrbCRM1, gives rise to both cone photoreceptors and horizontal cells. These cells do not express Vsx2, but co-express the transcription factors (TFs) Onecut1 and Otx2, which bind to ThrbCRM1. The components of the gene regulatory networks that control the transition from multipotent to fate-restricted gene expression are not known. This work aims to identify and evaluate cis-regulatory elements proximal to Onecut1 to identify the gene regulatory networks involved in RPC fate-restriction.

Method

We identified regulatory elements through ATAC-seq and conservation, followed by reporter assays to screen for activity based on temporal and spatial criteria. The regulatory elements of interest were subject to deletion and mutation analysis to identify functional sequences and evaluated by quantitative flow cytometry assays. Finally, we combined the enhancer::reporter assays with candidate TF overexpression to evaluate the relationship between the TFs, the enhancers, and early vertebrate retinal development. Statistical tests included ANOVA, Kruskal-Wallis, or unpaired t-tests.

Results

Two regulatory elements, ECR9 and ECR65, were identified to be active in ThrbCRM1(+) restricted RPCs. Candidate bHLH binding sites were identified as critical sequences in both elements. Overexpression of candidate bHLH TFs revealed specific enhancer-bHLH interactions. Nhlh1 overexpression expanded ECR65 activity into the Vsx2(+) RPC population, and overexpression of NeuroD1/NeuroG2/NeuroD4 had a similar effect on ECR9. Furthermore, bHLHs that were able to activate ectopic ECR9 reporter were able to induce endogenous Otx2 expression.

Conclusions

This work reports a large-scale screen to identify spatiotemporally specific regulatory elements near the Onecut1 locus. These elements were used to identify distinct populations in the developing retina. In addition, fate-restricted regulatory elements responded differentially to bHLH factors, and suggest a role for retinal bHLHs upstream of the Otx2 and Onecut1 genes during the formation of restricted RPCs from multipotent RPCs.