DNase I hypersensitivity analysis of the mouse brain and retina identifies region-specific regulatory elements

Cardiac GR Mediates the Diurnal Rhythm in Ventricular Arrhythmia Susceptibility

BACKGROUND

Ventricular arrhythmias (VAs) demonstrate a prominent day-night rhythm, commonly presenting in the morning. Transcriptional rhythms in cardiac ion channels accompany this phenomenon, but their role in the morning vulnerability to VAs and the underlying mechanisms are not understood. We investigated the recruitment of transcription factors that underpins transcriptional rhythms in ion channels and assessed whether this mechanism was pertinent to the heart's intrinsic diurnal susceptibility to VA.

METHODS AND RESULTS

Assay for transposase-accessible chromatin with sequencing performed in mouse ventricular myocyte nuclei at the beginning of the animals' inactive (ZT0) and active (ZT12) periods revealed differentially accessible chromatin sites annotating to rhythmically transcribed ion channels and distinct transcription factor binding motifs in these regions. Notably, motif enrichment for the glucocorticoid receptor (GR; transcriptional effector of corticosteroid signaling) in open chromatin profiles at ZT12 was observed, in line with the well-recognized ZT12 peak in circulating corticosteroids. Molecular, electrophysiological, and in silico biophysically-detailed modeling approaches demonstrated GR-mediated transcriptional control of ion channels (including Scn5a underlying the cardiac Na+ current, Kcnh2 underlying the rapid delayed rectifier K+ current, and Gja1 responsible for electrical coupling) and their contribution to the day-night rhythm in the vulnerability to VA. Strikingly, both pharmacological block of GR and cardiomyocyte-specific genetic knockout of GR blunted or abolished ion channel expression rhythms and abolished the ZT12 susceptibility to pacing-induced VA in isolated hearts.

CONCLUSIONS

Our study registers a day-night rhythm in chromatin accessibility that accompanies diurnal cycles in ventricular myocytes. Our approaches directly implicate the cardiac GR in the myocyte excitability rhythm and mechanistically link the ZT12 surge in glucocorticoids to intrinsic VA propensity at this time.

Pathogenic variants in CRX have distinct cis-regulatory effects on enhancers and silencers in photoreceptors

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.

Rho enhancers play unexpectedly minor roles in Rhodopsin transcription and rod cell integrity

Enhancers function with a basal promoter to control the transcription of target genes. Enhancer regulatory activity is often studied using reporter-based transgene assays. However, unmatched results have been reported when selected enhancers are silenced in situ. In this study, using genomic deletion analysis in mice, we investigated the roles of two previously identified enhancers and the promoter of the Rho gene that codes for the visual pigment rhodopsin. The Rho gene is robustly expressed by rod photoreceptors of the retina, and essential for the subcellular structure and visual function of rod photoreceptors. Mutations in RHO cause severe vision loss in humans. We found that each Rho regulatory region can independently mediate local epigenomic changes, but only the promoter is absolutely required for establishing active Rho chromatin configuration and transcription and maintaining the cell integrity and function of rod photoreceptors. To our surprise, two Rho enhancers that enable strong promoter activation in reporter assays are largely dispensable for Rho expression in vivo. Only small and age-dependent impact is detectable when both enhancers are deleted. Our results demonstrate context-dependent roles of enhancers and highlight the importance of studying functions of cis-regulatory regions in the native genomic context.

iDHS-FFLG: Identifying DNase I Hypersensitive Sites by Feature Fusion and Local-Global Feature Extraction Network

The DNase I hypersensitive sites (DHSs) are active regions on chromatin that have been found to be highly sensitive to DNase I. These regions contain various cis-regulatory elements, including promoters, enhancers and silencers. Accurate identification of DHSs helps researchers better understand the transcriptional machinery of DNA and deepen the knowledge of functional DNA elements in non-coding sequences. Researchers have developed many methods based on traditional experiments and machine learning to identify DHSs. However, low prediction accuracy and robustness limit their application in genetics research. In this paper, a novel computational approach based on deep learning is proposed by feature fusion and local-global feature extraction network to identify DHSs in mouse, named iDHS-FFLG. First of all, multiple binary features of nucleotides are fused to better express sequence information. Then, a network consisting of the convolutional neural network (CNN), bidirectional long short-term memory (BiLSTM) and self-attention mechanism is designed to extract local features and global contextual associations. In the end, the prediction module is applied to distinguish between DHSs and non-DHSs. The results of several experiments demonstrate the superior performances of iDHS-FFLG compared to the latest methods.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.

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).

Disease-causing mutations in genes encoding transcription factors critical for photoreceptor development

Photoreceptor development of the vertebrate visual system is controlled by a complex transcription regulatory network. OTX2 is expressed in the mitotic retinal progenitor cells (RPCs) and controls photoreceptor genesis. CRX that is activated by OTX2 is expressed in photoreceptor precursors after cell cycle exit. NEUROD1 is also present in photoreceptor precursors that are ready to specify into rod and cone photoreceptor subtypes. NRL is required for the rod fate and regulates downstream rod-specific genes including the orphan nuclear receptor NR2E3 which further activates rod-specific genes and simultaneously represses cone-specific genes. Cone subtype specification is also regulated by the interplay of several transcription factors such as THRB and RXRG. Mutations in these key transcription factors are responsible for ocular defects at birth such as microphthalmia and inherited photoreceptor diseases such as Leber congenital amaurosis (LCA), retinitis pigmentosa (RP) and allied dystrophies. In particular, many mutations are inherited in an autosomal dominant fashion, including the majority of missense mutations in CRX and NRL. In this review, we describe the spectrum of photoreceptor defects that are associated with mutations in the above-mentioned transcription factors, and summarize the current knowledge of molecular mechanisms underlying the pathogenic mutations. At last, we deliberate the outstanding gaps in our understanding of the genotype-phenotype correlations and outline avenues for future research of the treatment strategies.

Cell fate decisions, transcription factors and signaling during early retinal development

The development of the vertebrate eyes is a complex process starting from anterior-posterior and dorso-ventral patterning of the anterior neural tube, resulting in the formation of the eye field. Symmetrical separation of the eye field at the anterior neural plate is followed by two symmetrical evaginations to generate a pair of optic vesicles. Next, reciprocal invagination of the optic vesicles with surface ectoderm-derived lens placodes generates double-layered optic cups. The inner and outer layers of the optic cups develop into the neural retina and retinal pigment epithelium (RPE), respectively. In vitro produced retinal tissues, called retinal organoids, are formed from human pluripotent stem cells, mimicking major steps of retinal differentiation in vivo. This review article summarizes recent progress in our understanding of early eye development, focusing on the formation the eye field, optic vesicles, and early optic cups. Recent single-cell transcriptomic studies are integrated with classical in vivo genetic and functional studies to uncover a range of cellular mechanisms underlying early eye development. The functions of signal transduction pathways and lineage-specific DNA-binding transcription factors are dissected to explain cell-specific regulatory mechanisms underlying cell fate determination during early eye development. The functions of homeodomain (HD) transcription factors Otx2, Pax6, Lhx2, Six3 and Six6, which are required for early eye development, are discussed in detail. Comprehensive understanding of the mechanisms of early eye development provides insight into the molecular and cellular basis of developmental ocular anomalies, such as optic cup coloboma. Lastly, modeling human development and inherited retinal diseases using stem cell-derived retinal organoids generates opportunities to discover novel therapies for retinal diseases.

An enhancer located in a Pde6c intron drives transient expression in the cone photoreceptors of developing mouse and human retinas

How cone photoreceptors are formed during retinal development is only partially known. This is in part because we do not fully understand the gene regulatory network responsible for cone genesis. We reasoned that cis-regulatory elements (enhancers) active in nascent cones would be regulated by the same upstream network that controls cone formation. To dissect this network, we searched for enhancers active in developing cones. By electroporating enhancer-driven fluorescent reporter plasmids, we observed that a sequence within an intron of the cone-specific Pde6c gene acted as an enhancer in developing mouse cones. Similar fluorescent reporter plasmids were used to generate stable transgenic human induced pluripotent stem cells that were then grown into three-dimensional human retinal organoids. These organoids contained fluorescently labeled cones, demonstrating that the Pde6c enhancer was also active in human cones. We observed that enhancer activity was transient and labeled a minor population of developing rod photoreceptors in both mouse and human systems. This cone-enriched pattern argues that the Pde6c enhancer is activated in cells poised between rod and cone fates. Additionally, it suggests that the Pde6c enhancer is activated by the same regulatory network that selects or stabilizes cone fate choice. To further understand this regulatory network, we identified essential enhancer sequence regions through a series of mutagenesis experiments. This suggested that the Pde6c enhancer was regulated by transcription factor binding at five or more locations. Binding site predictions implicated transcription factor families known to control photoreceptor formation and families not previously associated with cone development. These results provide a framework for deciphering the gene regulatory network that controls cone genesis in both human and mouse systems. Our new transgenic human stem cell lines provide a tool for determining which cone developmental mechanisms are shared and distinct between mice and humans.

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.

Construction and characterization of EGFP reporter plasmid harboring putative human RAX promoter for in vitro monitoring of retinal progenitor cells identity

BACKGROUND

In retinal degenerative disease, progressive and debilitating conditions result in deterioration of retinal cells and visual loss. In human, retina lacks the inherent capacity for regeneration. Therefore, regeneration of retinal layer from human retinal progenitor cells (hRPCs) is a challenging task and restricted in vitro maintenance of hRPCs remains as the main hurdle. Retina and anterior neural fold homeobox gene (RAX) play critical roles in developing retina and maintenance of hRPCs. In this study, for the first time regulatory regions of human RAX gene with potential promoter activity were experimentally investigated.

RESULTS

For this purpose, after in silico analysis of regulatory regions of human RAX gene, the expression of EGFP reporter derived by putative promoter sequences was first evaluated in 293 T cells and then in hRPCS derived from human embryonic stem cells. The candidate region (RAX-3258 bp) showed the highest EGFP expression in hRPCs. This reporter construct can be used for in vitro monitoring of hRPC identity and verification of an efficient culture medium for maintenance of these cells.

CONCLUSIONS

Furthermore, our findings provide a platform for better insight into regulatory regions of human RAX gene and molecular mechanisms underlying its vital functions in retina development.

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.

The Atoh7 remote enhancer provides transcriptional robustness during retinal ganglion cell development

The retinal ganglion cell (RGC) competence factor ATOH7 is dynamically expressed during retinal histogenesis. ATOH7 transcription is controlled by a promoter-adjacent primary enhancer and a remote shadow enhancer (SE). Deletion of the ATOH7 human SE causes nonsyndromic congenital retinal nonattachment (NCRNA) disease, characterized by optic nerve aplasia and total blindness. We used genome editing to model NCRNA in mice. Deletion of the murine SE reduces Atoh7 messenger RNA (mRNA) fivefold but does not recapitulate optic nerve loss; however, SEdel/knockout (KO) trans heterozygotes have thin optic nerves. By analyzing Atoh7 mRNA and protein levels, RGC development and survival, and chromatin landscape effects, we show that the SE ensures robust Atoh7 transcriptional output. Combining SE deletion and KO and wild-type alleles in a genotypic series, we determined the amount of Atoh7 needed to produce a normal complement of adult RGCs, and the secondary consequences of graded reductions in Atoh7 dosage. Together, these data reveal the workings of an evolutionary fail-safe, a duplicate enhancer mechanism that is hard-wired in the machinery of vertebrate retinal ganglion cell genesis.

Developmental changes in the accessible chromatin, transcriptome and Ascl1-binding correlate with the loss in Müller Glial regenerative potential

Diseases and damage to the retina lead to losses in retinal neurons and eventual visual impairment. Although the mammalian retina has no inherent regenerative capabilities, fish have robust regeneration from Müller glia (MG). Recently, we have shown that driving expression of Ascl1 in adult mouse MG stimulates neural regeneration. The regeneration observed in the mouse is limited in the variety of neurons that can be derived from MG; Ascl1-expressing MG primarily generate bipolar cells. To better understand the limits of MG-based regeneration in mouse retinas, we used ATAC- and RNA-seq to compare newborn progenitors, immature MG (P8-P12), and mature MG. Our analysis demonstrated developmental differences in gene expression and accessible chromatin between progenitors and MG, primarily in neurogenic genes. Overexpression of Ascl1 is more effective in reprogramming immature MG, than mature MG, consistent with a more progenitor-like epigenetic landscape in the former. We also used ASCL1 ChIPseq to compare the differences in ASCL1 binding in progenitors and reprogrammed MG. We find that bipolar-specific accessible regions are more frequently linked to bHLH motifs and ASCL1 binding. Overall, our analysis indicates a loss of neurogenic gene expression and motif accessibility during glial maturation that may prevent efficient reprogramming.

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.

Prdm1 overexpression causes a photoreceptor fate-shift in nascent, but not mature, bipolar cells

The transcription factors Prdm1 (Blimp1) and Vsx2 (Chx10) work downstream of Otx2 to regulate photoreceptor and bipolar cell fates in the developing retina. Mice that lack Vsx2 fail to form bipolar cells while Prdm1 mutants form excess bipolars at the direct expense of photoreceptors. Excess bipolars in Prdm1 mutants appear to derive from rods, suggesting that photoreceptor fate remains mutable for some time after cells become specified. Here we tested whether bipolar cell fate is also plastic during development. To do this, we created a system to conditionally misexpress Prdm1 at different stages of bipolar cell development. We found that Prdm1 blocks bipolar cell formation if expressed before the fate choice decision occurred. When we misexpressed Prdm1 just after the decision to become a bipolar cell was made, some cells were reprogrammed into photoreceptors. In contrast, Prdm1 misexpression in mature bipolar cells did not affect cell fate. We also provide evidence that sustained misexpression of Prdm1 was selectively toxic to photoreceptors. Our data show that bipolar fate is malleable, but only for a short temporal window following fate specification. Prdm1 and Vsx2 act by stabilizing photoreceptor and bipolar fates in developing OTX2+ cells of the retina.

Mapping the cis-regulatory architecture of the human retina reveals noncoding genetic variation in disease

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.

Targeted deletion of an NRL- and CRX-regulated alternative promoter specifically silences FERM and PDZ domain containing 1 (Frmpd1) in rod photoreceptors

Regulation of cell type-specific gene expression is critical for generating neuronal diversity. Transcriptome analyses have unraveled extensive heterogeneity of transcribed sequences in retinal photoreceptors because of alternate splicing and/or promoter usage. Here we show that Frmpd1 (FERM and PDZ domain containing 1) is transcribed from an alternative promoter specifically in the retina. Electroporation of Frmpd1 promoter region, -505 to +382 bp, activated reporter gene expression in mouse retina in vivo. A proximal promoter sequence (-8 to +33 bp) of Frmpd1 binds to neural retina leucine zipper (NRL) and cone-rod homeobox protein (CRX), two rod-specific differentiation factors, and is necessary for activating reporter gene expression in vitro and in vivo. Clustered regularly interspaced short palindromic repeats/Cas9-mediated deletion of the genomic region, including NRL and CRX binding sites, in vivo completely eliminated Frmpd1 expression in rods and dramatically reduced expression in rod bipolar cells, thereby overcoming embryonic lethality caused by germline Frmpd1 deletion. Our studies demonstrate that a cell type-specific regulatory control region is a credible target for creating loss-of-function alleles of widely expressed genes.

STAT Signaling Modifies Ascl1 Chromatin Binding and Limits Neural Regeneration from Muller Glia in Adult Mouse Retina

Müller glia (MG) serve as sources for retinal regeneration in non-mammalian vertebrates. We find that this process can be induced in mouse MG, after injury, by transgenic expression of the proneural transcription factor Ascl1 and the HDAC inhibitor TSA. However, new neurons are generated only from a subset of MG. Identifying factors that limit Ascl1-mediated MG reprogramming could make this process more efficient. In this study, we test whether injury-induced STAT activation hampers the ability of Ascl1 to reprogram MG into retinal neurons. Single-cell RNA-seq shows that progenitor-like cells derived from Ascl1-expressing MG have a higher level of STAT signaling than do those cells that become neurons. Ascl1-ChIPseq and ATAC-seq show that STAT potentially directs Ascl1 to developmentally inappropriate targets. Using a STAT inhibitor, in combination with our previously described reprogramming paradigm, we found a large increase in the ability of MG to generate neurons.

Enhancer transcription identifies cis-regulatory elements for photoreceptor cell types

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.

Sleep-wake-driven and circadian contributions to daily rhythms in gene expression and chromatin accessibility in the murine cortex

The timing and duration of sleep results from the interaction between a homeostatic sleep-wake-driven process and a periodic circadian process, and involves changes in gene regulation and expression. Unraveling the contributions of both processes and their interaction to transcriptional and epigenomic regulatory dynamics requires sampling over time under conditions of unperturbed and perturbed sleep. We profiled mRNA expression and chromatin accessibility in the cerebral cortex of mice over a 3-d period, including a 6-h sleep deprivation (SD) on day 2. We used mathematical modeling to integrate time series of mRNA expression data with sleep-wake history, which established that a large proportion of rhythmic genes are governed by the homeostatic process with varying degrees of interaction with the circadian process, sometimes working in opposition. Remarkably, SD caused long-term effects on gene-expression dynamics, outlasting phenotypic recovery, most strikingly illustrated by a damped oscillation of most core clock genes, including Arntl/Bmal1, suggesting that enforced wakefulness directly impacts the molecular clock machinery. Chromatin accessibility proved highly plastic and dynamically affected by SD. Dynamics in distal regions, rather than promoters, correlated with mRNA expression, implying that changes in expression result from constitutively accessible promoters under the influence of enhancers or repressors. Serum response factor (SRF) was predicted as a transcriptional regulator driving immediate response, suggesting that SRF activity mirrors the build-up and release of sleep pressure. Our results demonstrate that a single, short SD has long-term aftereffects at the genomic regulatory level and highlights the importance of the sleep-wake distribution to diurnal rhythmicity and circadian processes.

Open chromatin dynamics in prosensory cells of the embryonic mouse cochlea

Hearing loss is often due to the absence or the degeneration of hair cells in the cochlea. Understanding the mechanisms regulating the generation of hair cells may therefore lead to better treatments for hearing disorders. To elucidate the transcriptional control mechanisms specifying the progenitor cells (i.e. prosensory cells) that generate the hair cells and support cells critical for hearing function, we compared chromatin accessibility using ATAC-seq in sorted prosensory cells (Sox2-EGFP⁺) and surrounding cells (Sox2-EGFP⁻) from E12, E14.5 and E16 cochlear ducts. In Sox2-EGFP⁺, we find greater accessibility in and near genes restricted in expression to the prosensory region of the cochlear duct including Sox2, Isl1, Eya1 and Pou4f3. Furthermore, we find significant enrichment for the consensus binding sites of Sox2, Six1 and Gata3—transcription factors required for prosensory development—in the open chromatin regions. Over 2,200 regions displayed differential accessibility with developmental time in Sox2-EGFP⁺ cells, with most changes in the E12-14.5 window. Open chromatin regions detected in Sox2-EGFP⁺ cells map to over 48,000 orthologous regions in the human genome that include regions in genes linked to deafness. Our results reveal a dynamic landscape of open chromatin in prosensory cells with potential implications for cochlear development and disease.

Epigenetics in neuronal regeneration

Damage to neuronal tissues in mammals leads to permanent loss of tissue function that can have major health consequences. While mammals have no inherent regenerative capacity to functionally repair neuronal tissue, other species such as amphibians and teleost fish readily replace damaged tissue. The exploration of development and native regeneration can thus inform the process of inducing regeneration in non-regenerative systems, which can be used to develop new therapeutics. Increasing evidence points to an epigenetic component in the regulation of the changes in cellular gene expression necessary for regeneration. In this review, we compare evidence of epigenetic roles in development and regeneration of neuronal tissue. We have focused on three key systems of important clinical significance: the neural retina, the inner ear, and the spinal cord in regenerative and non-regenerative species. While evidence for epigenetic regulation of regeneration is still limited, changes in DNA accessibility, histone acetylation and DNA methylation have all emerged as key elements in this process. To date, most studies have used broadly acting experimental manipulations to establish a role for epigenetics in regeneration, but the advent of more targeted approaches to modify the epigenome will be critical to dissecting the relative contributions of these regulatory factors in this process and the development of methods to stimulate the regeneration in those organisms like ourselves where only limited regeneration occurs in these neural systems.

HOX Loci Focused CRISPR/sgRNA Library Screening Identifying Critical CTCF Boundaries

CCCTC-binding factor (CTCF)-mediated stable topologically associating domains (TADs) play a critical role in constraining interactions of DNA elements that are located in neighboring TADs. CTCF plays an important role in regulating the spatial and temporal expression of HOX genes that control embryonic development, body patterning, hematopoiesis, and leukemogenesis. However, it remains largely unknown whether and how HOX loci associated CTCF boundaries regulate chromatin organization and HOX gene expression. In the current protocol, a specific sgRNA pooled library targeting all CTCF binding sites in the HOXA/B/C/D loci has been generated to examine the effects of disrupting CTCF-associated chromatin boundaries on TAD formation and HOX gene expression. Through CRISPR-Cas9 genetic screening, the CTCF binding site located between HOXA7/HOXA9 genes (CBS7/9) has been identified as a critical regulator of oncogenic chromatin domain, as well as being important for maintaining ectopic HOX gene expression patterns in MLL-rearranged acute myeloid leukemia (AML). Thus, this sgRNA library screening approach provides novel insights into CTCF mediated genome organization in specific gene loci and also provides a basis for the functional characterization of the annotated genetic regulatory elements, both coding and noncoding, during normal biological processes in the post-human genome project era.

Lmx1a drives Cux2 expression in the cortical hem through activation of a conserved intronic enhancer

During neocortical development, neurons are produced by a diverse pool of neural progenitors. A subset of progenitors express the Cux2 gene and are fate restricted to produce certain neuronal subtypes; however, the upstream pathways that specify these progenitor fates remain unknown. To uncover the transcriptional networks that regulate Cux2 expression in the forebrain, we characterized a conserved Cux2 enhancer that recapitulates Cux2 expression specifically in the cortical hem. Using a bioinformatic approach, we identified putative transcription factor (TF)-binding sites for cortical hem-patterning TFs. We found that the homeobox TF Lmx1a can activate the Cux2 enhancer in vitro Furthermore, we showed that Lmx1a-binding sites were required for enhancer activity in the cortical hem in vivo Mis-expression of Lmx1a in hippocampal progenitors caused an increase in Cux2 enhancer activity outside the cortical hem. Finally, we compared several human enhancers with cortical hem-restricted activity and found that recurrent Lmx1a-binding sites are a top shared feature. Uncovering the network of TFs involved in regulating Cux2 expression will increase our understanding of the mechanisms pivotal in establishing Cux2 lineage fates in the developing forebrain.

A massively parallel reporter assay reveals context-dependent activity of homeodomain binding sites in vivo

Cone-rod homeobox (CRX) is a paired-like homeodomain transcription factor (TF) and a master regulator of photoreceptor development in vertebrates. The in vitro DNA binding preferences of CRX have been described in detail, but the degree to which in vitro binding affinity is correlated with in vivo enhancer activity is not known. In addition, paired-class homeodomain TFs can bind DNA cooperatively as both homodimers and heterodimers at inverted TAAT half-sites separated by 2 or 3 nucleotides. This dimeric configuration is thought to mediate target specificity, but whether monomeric and dimeric sites encode distinct levels of activity is not known. Here, we used a massively parallel reporter assay to determine how local sequence context shapes the regulatory activity of CRX binding sites in mouse photoreceptors. We assayed inactivating mutations in more than 1700 TF binding sites and found that dimeric CRX binding sites act as stronger enhancers than monomeric CRX binding sites. Furthermore, the activity of dimeric half-sites is cooperative, dependent on a strict 3-bp spacing, and tuned by the identity of the spacer nucleotides. Saturating single-nucleotide mutagenesis of 195 CRX binding sites showed that, on average, changes in TF binding site affinity are correlated with changes in regulatory activity, but this relationship is obscured when considering mutations across multiple cis-regulatory elements (CREs). Taken together, these results demonstrate that the activity of CRX binding sites is highly dependent on sequence context, providing insight into photoreceptor gene regulation and illustrating functional principles of homeodomain binding sites that may be conserved in other cell types.

CRX directs photoreceptor differentiation by accelerating chromatin remodeling at specific target sites

Background

Recent technological advances have delivered the genome-wide targets of many important transcription factors (TFs). However, increasing evidence suggests that not all target sites mediate regulatory function, raising the questions of how to determine which sites are active, what are the epigenetic consequences of TF binding at these sites, and how the specificity is coded. To address these questions, we focused on CRX, a disease-associated homeodomain TF required for photoreceptor gene expression and development. Since CRX binds more than 6000 sites across the genome in the retina, we profiled chromatin landscape changes at each binding site during normal development and in the absence of CRX and interpreted the results by thorough investigation of other epigenomic datasets and sequence features.

Results

CRX is required for chromatin remodeling at only a subset of its binding sites, which undergo retina or neuronal specific activation during photoreceptor differentiation. Genes near these “CRX Dependent” sites code for proteins important for photoreceptor physiology and function, and their transcription is significantly reduced in Crx deficient retinas. In addition, the nucleotide and motif content distinguish these CRX Dependent sites from other CRX-bound sites.

Conclusions

Together, our results suggest that CRX acts only at select, uniquely-coded binding sites to accelerate chromatin remodeling during photoreceptor differentiation. This study emphasizes the importance of connecting TF binding with its functional consequences and provides a framework for making such a connection using comparative analyses of available genomic datasets. Finally, this study prioritizes sets of non-coding DNA sites for future functional interrogation and identification of mutations associated with retinal disease.

Electronic supplementary material

The online version of this article (10.1186/s13072-018-0212-2) contains supplementary material, which is available to authorized users.

Open chromatin dynamics reveals stage-specific transcriptional networks in hiPSC-based neurodevelopmental model

Chromatin accessibility to transcription factors (TFs) strongly influences gene transcription and cell differentiation. However, a mechanistic understanding of the transcriptional control during the neuronal differentiation of human induced pluripotent stem cells (hiPSCs), a promising cellular model for mental disorders, remains elusive. Here, we carried out additional analyses on our recently published open chromatin regions (OCRs) profiling at different stages of hiPSC neuronal differentiation. We found that the dynamic changes of OCR during neuronal differentiation highlighted cell stage-specific gene networks, and the chromatin accessibility at the core promoter region of a gene correlates with the corresponding transcript abundance. Within the cell stage-specific OCRs, we identified the binding of cell stage-specific TFs and observed a lag of a neuronal TF binding behind the mRNA expression of the corresponding TF. Interestingly, binding footprints of NEUROD1 and NEUROG2, both of which induce high efficient conversion of hiPSCs to glutamatergic neurons, were among those most enriched in the relatively mature neurons. Furthermore, TF network analysis showed that both NEUROD1 and NEUROG2 were present in the same core TF network specific to more mature neurons, suggesting a pivotal mechanism of epigenetic control of neuronal differentiation and maturation. Our study provides novel insights into the epigenetic control of glutamatergic neurogenesis in the context of TF networks, which may be instrumental to improving hiPSC modeling of neuropsychiatric disorders.

Global analysis of primary mesenchyme cell cis-regulatory modules by chromatin accessibility profiling

Background

The developmental gene regulatory network (GRN) that underlies skeletogenesis in sea urchins and other echinoderms is a paradigm of GRN structure, function, and evolution. This transcriptional network is deployed selectively in skeleton-forming primary mesenchyme cells (PMCs) of the early embryo. To advance our understanding of this model developmental GRN, we used genome-wide chromatin accessibility profiling to identify and characterize PMC cis-regulatory modules (CRMs).

Results

ATAC-seq (Assay for Transposase-Accessible Chromatin using sequencing) analysis of purified PMCs provided a global picture of chromatin accessibility in these cells. We used both ATAC-seq and DNase-seq (DNase I hypersensitive site sequencing) to identify > 3000 sites that exhibited increased accessibility in PMCs relative to other embryonic cell lineages, and provide both computational and experimental evidence that a large fraction of these sites represent bona fide skeletogenic CRMs. Putative PMC CRMs were preferentially located near genes differentially expressed by PMCs and consensus binding sites for two key transcription factors in the PMC GRN, Alx1 and Ets1, were enriched in these CRMs. Moreover, a high proportion of candidate CRMs drove reporter gene expression specifically in PMCs in transgenic embryos. Surprisingly, we found that PMC CRMs were partially open in other embryonic lineages and exhibited hyperaccessibility as early as the 128-cell stage.

Conclusions

Our work provides a comprehensive picture of chromatin accessibility in an early embryonic cell lineage. By identifying thousands of candidate PMC CRMs, we significantly enhance the utility of the sea urchin skeletogenic network as a general model of GRN architecture and evolution. Our work also shows that differential chromatin accessibility, which has been used for the high-throughput identification of enhancers in differentiated cell types, is a powerful approach for the identification of CRMs in early embryonic cells. Lastly, we conclude that in the sea urchin embryo, CRMs that control the cell type-specific expression of effector genes are hyperaccessible several hours in advance of gene activation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4542-z) contains supplementary material, which is available to authorized users.

Prdm13 is required for Ebf3+ amacrine cell formation in the retina

Amacrine interneurons play a critical role in the processing of visual signals within the retina. They are highly diverse, representing 30 or more distinct subtypes. Little is known about how amacrine subtypes acquire their unique gene expression and morphological features. We characterized the gene expression pattern of the zinc-finger transcription factor Prdm13 in the mouse. Consistent with a developmental role, Prdm13 was expressed by Ptf1a+ amacrine and horizontal precursors. Over time, Prdm13 expression diverged from the transiently expressed Ptf1a and marked just a subset of amacrine cells in the adult retina. While heterogeneous, we show that most of these Prdm13+ amacrine cells express the transcription factor Ebf3 and the calcium binding protein calretinin. Loss of Prdm13 did not affect the number of amacrine cells formed during development. However, we observed a modest loss of amacrine cells and increased apoptosis that correlated with the onset timing of Ebf3 expression. Adult Prdm13 loss-of-function mice had 25% fewer amacrine cells, altered calretinin expression, and a lack of Ebf3+ amacrines. Forcing Prdm13 expression in retinal progenitor cells did not significantly increase amacrine cell formation, Ebf3 or calretinin expression, and appeared detrimental to the survival of photoreceptors. Our data show that Prdm13 is not required for amacrine fate as a class, but is essential for the formation of Ebf3+ amacrine cell subtypes. Rather than driving subtype identity, Prdm13 may act by restricting competing fate programs to maintain identity and survival.

Cis-regulatory landscapes of four cell types of the retina

The retina is composed of ∼50 cell-types with specific functions for the process of vision. Identification of the cis-regulatory elements active in retinal cell-types is key to elucidate the networks controlling this diversity. Here, we combined transcriptome and epigenome profiling to map the regulatory landscape of four cell-types isolated from mouse retinas including rod and cone photoreceptors as well as rare inter-neuron populations such as horizontal and starburst amacrine cells. Integration of this information reveals sequence determinants and candidate transcription factors for controlling cellular specialization. Additionally, we refined parallel reporter assays to enable studying the transcriptional activity of large collection of sequences in individual cell-types isolated from a tissue. We provide proof of concept for this approach and its scalability by characterizing the transcriptional capacity of several hundred putative regulatory sequences within individual retinal cell-types. This generates a catalogue of cis-regulatory regions active in retinal cell types and we further demonstrate their utility as potential resource for cellular tagging and manipulation.

Combinatorial regulation of a Blimp1 (Prdm1) enhancer in the mouse retina

The mouse retina comprises seven major cell types that exist in differing proportions. They are generated from multipotent progenitors in a stochastic manner, such that the relative frequency of any given type generated changes over time. The mechanisms determining the proportions of each cell type are only partially understood. Photoreceptors and bipolar interneurons are derived from cells that express Otx2. Within this population, Blimp1 (Prdm1) helps set the balance between photoreceptors and bipolar cells by suppressing bipolar identity in most of the cells. How only a subset of these Otx2+ cells decides to upregulate Blimp1 and adopt photoreceptor fate is unknown. To understand this, we investigated how Blimp1 transcription is regulated. We identified several potential Blimp1 retinal enhancer elements using DNase hypersensitivity sequencing. Only one of the elements recapitulated Blimp1 spatial and temporal expression in cultured explant assays and within the retinas of transgenic mice. Mutagenesis of this retinal Blimp1 enhancer element revealed four discrete sequences that were each required for its activity. These included highly conserved Otx2 and ROR (retinoic acid receptor related orphan receptor) binding sites. The other required sequences do not appear to be controlled by Otx2 or ROR factors, increasing the complexity of the Blimp1 gene regulatory network. Our results show that the intersection of three or more transcription factors is required to correctly regulate the spatial and temporal features of Blimp1 enhancer expression. This explains how Blimp1 expression can diverge from Otx2 and set the balance between photoreceptor and bipolar fates.

Limb-Enhancer Genie: An accessible resource of accurate enhancer predictions in the developing limb

Epigenomic mapping of enhancer-associated chromatin modifications facilitates the genome-wide discovery of tissue-specific enhancers in vivo. However, reliance on single chromatin marks leads to high rates of false-positive predictions. More sophisticated, integrative methods have been described, but commonly suffer from limited accessibility to the resulting predictions and reduced biological interpretability. Here we present the Limb-Enhancer Genie (LEG), a collection of highly accurate, genome-wide predictions of enhancers in the developing limb, available through a user-friendly online interface. We predict limb enhancers using a combination of >50 published limb-specific datasets and clusters of evolutionarily conserved transcription factor binding sites, taking advantage of the patterns observed at previously in vivo validated elements. By combining different statistical models, our approach outperforms current state-of-the-art methods and provides interpretable measures of feature importance. Our results indicate that including a previously unappreciated score that quantifies tissue-specific nuclease accessibility significantly improves prediction performance. We demonstrate the utility of our approach through in vivo validation of newly predicted elements. Moreover, we describe general features that can guide the type of datasets to include when predicting tissue-specific enhancers genome-wide, while providing an accessible resource to the general biological community and facilitating the functional interpretation of genetic studies of limb malformations.

Identification and prediction of alternative transcription start sites that generate rod photoreceptor-specific transcripts from ubiquitously expressed genes

Transcriptome complexity is substantially increased by the use of multiple transcription start sites for a given gene. By utilizing a rod photoreceptor-specific chromatin signature, and the RefSeq database of established transcription start sites, we have identified essentially all known rod photoreceptor genes as well as a group of novel genes that have a high probability of being expressed in rod photoreceptors. Approximately half of these novel rod genes are transcribed into multiple mRNA and/or protein isoforms through alternative transcriptional start sites (ATSS), only one of which has a rod-specific epigenetic signature and gives rise to a rod transcript. This suggests that, during retina development, some genes use ATSS to regulate cell type and temporal specificity, effectively generating a rod transcript from otherwise ubiquitously expressed genes. Biological confirmation of the relationship between epigenetic signatures and gene expression, as well as comparison of our genome-wide chromatin signature maps with available data sets for retina, namely a ChIP-on-Chip study of Polymerase-II (Pol-II) binding sites, ChIP-Seq studies for NRL- and CRX- binding sites and DHS (University of Washington data, available on UCSC mouse Genome Browser as a part of ENCODE project) fully support our hypothesis and together accurately identify and predict an array of new rod transcripts. The same approach was used to identify a number of TSS that are not currently in RefSeq. Biological conformation of the use of some of these TSS suggests that this method will be valuable for exploring the range of transcriptional complexity in many tissues. Comparison of mouse and human genome-wide data indicates that most of these alternate TSS appear to be present in both species, indicating that our approach can be useful for identification of regulatory regions that might play a role in human retinal disease.

Cell Type-Specific Epigenomic Analysis Reveals a Uniquely Closed Chromatin Architecture in Mouse Rod Photoreceptors

Rod photoreceptors are specialized neurons that mediate vision in dim light and are the predominant photoreceptor type in nocturnal mammals. The rods of nocturnal mammals are unique among vertebrate cell types in having an ‘inverted’ nuclear architecture, with a dense mass of heterochromatin in the center of the nucleus rather than dispersed clumps at the periphery. To test if this unique nuclear architecture is correlated with a unique epigenomic landscape, we performed ATAC-seq on mouse rods and their most closely related cell type, cone photoreceptors. We find that thousands of loci are selectively closed in rods relative to cones as well as >60 additional cell types. Furthermore, we find that the open chromatin profile of photoreceptors lacking the rod master regulator Nrl is nearly indistinguishable from that of native cones, indicating that Nrl is required for selective chromatin closure in rods. Finally, we identified distinct enrichments of transcription factor binding sites in rods and cones, revealing key differences in the cis-regulatory grammar of these cell types. Taken together, these data provide insight into the development and maintenance of photoreceptor identity, and highlight rods as an attractive system for studying the relationship between nuclear organization and local changes in gene regulation.

Retinoid-Related Orphan Receptor β and Transcriptional Control of Neuronal Differentiation

The ability to generate neuronal diversity is central to the function of the nervous system. Here we discuss the key neurodevelopmental roles of retinoid-related orphan receptor β (RORβ) encoded by the Rorb (Nr1f2) gene. Recent studies have reported loss of function of the human RORB gene in cases of familial epilepsy and intellectual disability. Principal sites of expression of the Rorb gene in model species include sensory organs, the spinal cord, and brain regions that process sensory and circadian information. Genetic analyses in mice have indicated functions in circadian behavior, vision, and, at the cellular level, the differentiation of specific neuronal cell types. Studies in the retina and sensory areas of the cerebral cortex suggest that this orphan nuclear receptor acts at decisive steps in transcriptional hierarchies that determine neuronal diversity.

Histone Deacetylase 1 Is Essential for Rod Photoreceptor Differentiation by Regulating Acetylation at Histone H3 Lysine 9 and Histone H4 Lysine 12 in the Mouse Retina

Histone acetylation has a regulatory role in gene expression and is necessary for proper tissue development. To investigate the specific roles of histone deacetylases (HDACs) in rod differentiation in neonatal mouse retinas, we used a pharmacological approach that showed that inhibition of class I but not class IIa HDACs caused the same phenotypic changes seen with broad spectrum HDAC inhibitors, most notably a block in the differentiation of rod photoreceptors. Inhibition of HDAC1 resulted in increase of acetylation of lysine 9 of histone 3 (H3K9) and lysine 12 of histone 4 (H4K12) but not lysine 27 of histone 3 (H3K27) and led to maintained expression of progenitor-specific genes such as Vsx2 and Hes1 with concomitant block of expression of rod-specific genes. ChiP experiments confirmed these changes in the promoters of a group of progenitor genes. Based on our results, we suggest that HDAC1-specific inhibition prevents progenitor cells of the retina from exiting the cell cycle and differentiating. HDAC1 may be an essential epigenetic regulator of the transition from progenitor cells to terminally differentiated photoreceptors.

Identification of Regulatory DNA Elements Using Genome-wide Mapping of DNase I Hypersensitive Sites during Tomato Fruit Development

Development and ripening of tomato fruit are precisely controlled by transcriptional regulation, which depends on the orchestrated accessibility of regulatory proteins to promoters and other cis-regulatory DNA elements. This accessibility and its effect on gene expression play a major role in defining the developmental process. To understand the regulatory mechanism and functional elements modulating morphological and anatomical changes during fruit development, we generated genome-wide high-resolution maps of DNase I hypersensitive sites (DHSs) from the fruit tissues of the tomato cultivar "Moneymaker" at 20 days post anthesis as well as break stage. By exploring variation of DHSs across fruit development stages, we pinpointed the most likely hypersensitive sites related to development-specific genes. By detecting binding motifs on DHSs of these development-specific genes or genes in the ascorbic acid biosynthetic pathway, we revealed the common regulatory elements contributing to coordinating gene transcription of plant ripening and specialized metabolic pathways. Our results contribute to a better understanding of the regulatory dynamics of genes involved in tomato fruit development and ripening.

Towards genome-wide prediction and characterization of enhancers in plants

Enhancers are important cis-regulatory DNA elements that regulate transcription programs by recruiting transcription factors and directing them to the promoters of target genes in a cell-type/tissue-specific manner. The expression of a gene can be regulated by one or multiple enhancers at different developmental stages and/or in different tissues. Enhancers are difficult to identify because of their unpredictable positions relative to their cognate promoters. Remarkably, only a handful of enhancers have been identified in plant species largely due to the lack of general approaches for enhancer identification. Extensive genomic and epigenomic research in mammalian species has revealed that the genomic locations of enhancers can be predicted based on the binding sites of transcriptional co-factors and several distinct features associated with open chromatin. Here we review the methodologies used in enhancer prediction in mammalian species. We also review the recent applications of these methodologies in Arabidopsis thaliana and discuss the future directions of enhancer identification in plants. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Next generation sequencing technology and genomewide data analysis: Perspectives for retinal research

The advent of high throughput next generation sequencing (NGS) has accelerated the pace of discovery of disease-associated genetic variants and genomewide profiling of expressed sequences and epigenetic marks, thereby permitting systems-based analyses of ocular development and disease. Rapid evolution of NGS and associated methodologies presents significant challenges in acquisition, management, and analysis of large data sets and for extracting biologically or clinically relevant information. Here we illustrate the basic design of commonly used NGS-based methods, specifically whole exome sequencing, transcriptome, and epigenome profiling, and provide recommendations for data analyses. We briefly discuss systems biology approaches for integrating multiple data sets to elucidate gene regulatory or disease networks. While we provide examples from the retina, the NGS guidelines reviewed here are applicable to other tissues/cell types as well.

Temporal establishment of neural cell identity in vivo and in vitro

Understanding cell fate specification is particularly useful because it enables biologists to generate specific neural cell types for treating currently untreatable neurological diseases. Traditionally, lineage-specific progenitors are generated in vitro from pluripotent cells, after which they may be channeled into more mature cell types in a stage-specific manner, which is similar to the way cells behave during development. However, the emergence of induced pluripotent stem cells means that specific cell types can be generated directly from fibroblasts or other somatic cell types, thus bypassing all of the necessary steps that happen in vivo. Based on this information, the present review first explores the regulatory circuitry that drives cell fate specification over time in vivo. In particular, it describes how the appearance of specific neuronal and glial cell types is governed by an intrinsic biological clock, followed by a discussion of how this can be achieved through the temporal expression of intracellular regulators in relation to cell-specific Dnase I hypersensitivity sites, promoters and enhancers. Cell fate acquisition in vitro was then examined in an attempt to evaluate whether the temporal regulation neural cell fate in vivo is still relevant to the generation of reprogrammed neural stem cells and neurons. Copyright © 2016 John Wiley & Sons, Ltd.

Epigenomic landscapes of retinal rods and cones

Rod and cone photoreceptors are highly similar in many respects but they have important functional and molecular differences. Here, we investigate genome-wide patterns of DNA methylation and chromatin accessibility in mouse rods and cones and correlate differences in these features with gene expression, histone marks, transcription factor binding, and DNA sequence motifs. Loss of NR2E3 in rods shifts their epigenomes to a more cone-like state. The data further reveal wide differences in DNA methylation between retinal photoreceptors and brain neurons. Surprisingly, we also find a substantial fraction of DNA hypo-methylated regions in adult rods that are not in active chromatin. Many of these regions exhibit hallmarks of regulatory regions that were active earlier in neuronal development, suggesting that these regions could remain undermethylated due to the highly compact chromatin in mature rods. This work defines the epigenomic landscapes of rods and cones, revealing features relevant to photoreceptor development and function.

DOI:http://dx.doi.org/10.7554/eLife.11613.001

Vision in humans is made possible by a light-sensing sheet of cells at the back of the eye called the retina. The surface of the retina is populated by specialized sensory cells, known as rods and cones. The rod cells detect very dim light, while the cones are less sensitive to light but are used to detect color. Together, the rods and cones gather the information needed to create a picture that is then transmitted to the brain.

Rods and cones have been studied for decades, and genetic analyses have revealed the patterns of gene expression that lead a cell to develop into either a rod or a cone. Researchers have also identified several key regulatory genes that control these patterns, but less is known about the role of other factors that control the expression of genes.

Chemical modifications to DNA or modifications to the proteins associated with DNA – which are collectively called epigenetic modifications – can either promote or inhibit the activation of nearby genes. Now, Mo et al. have shown that rods and cones from mice have very different patterns of epigenetic modifications. The experiments also revealed that many sections of DNA that are marked to promote gene activation contain known rod-specific or cone-specific genes; and that rod cells need a known regulatory gene to develop their specific pattern of epigenetic modifications. Finally, Mo et al. showed that epigenetic regulation differed between brain cells and rods and cones.

These insights into epigenetic regulation of rod and cone genes may help explain why some people with eye diseases caused by the same genetic mutation may develop symptoms at different ages or lose vision at different rates. The new information about gene regulation may also help scientists to reprogram stem cells to become healthy rods or cones that could be transplanted into people with eye disease to restore their vision.

DOI:http://dx.doi.org/10.7554/eLife.11613.002

Massively parallel cis-regulatory analysis in the mammalian central nervous system

Cis-regulatory elements (CREs, e.g., promoters and enhancers) regulate gene expression, and variants within CREs can modulate disease risk. Next-generation sequencing has enabled the rapid generation of genomic data that predict the locations of CREs, but a bottleneck lies in functionally interpreting these data. To address this issue, massively parallel reporter assays (MPRAs) have emerged, in which barcoded reporter libraries are introduced into cells, and the resulting barcoded transcripts are quantified by next-generation sequencing. Thus far, MPRAs have been largely restricted to assaying short CREs in a limited repertoire of cultured cell types. Here, we present two advances that extend the biological relevance and applicability of MPRAs. First, we adapt exome capture technology to instead capture candidate CREs, thereby tiling across the targeted regions and markedly increasing the length of CREs that can be readily assayed. Second, we package the library into adeno-associated virus (AAV), thereby allowing delivery to target organs in vivo. As a proof of concept, we introduce a capture library of about 46,000 constructs, corresponding to roughly 3500 DNase I hypersensitive (DHS) sites, into the mouse retina by ex vivo plasmid electroporation and into the mouse cerebral cortex by in vivo AAV injection. We demonstrate tissue-specific cis-regulatory activity of DHSs and provide examples of high-resolution truncation mutation analysis for multiplex parsing of CREs. Our approach should enable massively parallel functional analysis of a wide range of CREs in any organ or species that can be infected by AAV, such as nonhuman primates and human stem cell–derived organoids.

Transgenic expression of the proneural transcription factor Ascl1 in Müller glia stimulates retinal regeneration in young mice

Müller glial cells are the source of retinal regeneration in fish and birds; although this process is efficient in fish, it is less so in birds and very limited in mammals. It has been proposed that factors necessary for providing neurogenic competence to Müller glia in fish and birds after retinal injury are not expressed in mammals. One such factor, the proneural transcription factor Ascl1, is necessary for retinal regeneration in fish but is not expressed after retinal damage in mice. We previously reported that forced expression of Ascl1 in vitro reprograms Müller glia to a neurogenic state. We now test whether forced expression of Ascl1 in mouse Müller glia in vivo stimulates their capacity for retinal regeneration. We find that transgenic expression of Ascl1 in adult Müller glia in undamaged retina does not overtly affect their phenotype; however, when the retina is damaged, the Ascl1-expressing glia initiate a response that resembles the early stages of retinal regeneration in zebrafish. The reaction to injury is even more pronounced in Müller glia in young mice, where the Ascl1-expressing Müller glia give rise to amacrine and bipolar cells and photoreceptors. DNaseI-seq analysis of the retina and Müller glia shows progressive reduction in accessibility of progenitor gene cis-regulatory regions consistent with the reduction in their reprogramming. These results show that at least one of the differences between mammal and fish Müller glia that bears on their difference in regenerative potential is the proneural transcription factor Ascl1.