The Dynamic Epigenetic Landscape of the Retina During Development, Reprogramming, and Tumorigenesis

Bioinformatics analysis of multi-omics data identifying molecular biomarker candidates and epigenetically regulatory targets associated with retinoblastoma

Retinoblastoma (RB) is the commonest malignant tumor of the infant retina. Besides genetic changes, epigenetic events are also considered to implicate the occurrence of RB. This study aimed to identify significantly altered protein-coding genes, DNA methylation, microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and their molecular functions and pathways associated with RB, and investigate the epigenetically regulatory mechanism of DNA methylation modification and non-coding RNAs on key genes of RB via bioinformatics method.We obtained multi-omics data on protein-coding genes, DNA methylation, miRNAs, and lncRNAs from the Gene Expression Omnibus database. We identified differentially expressed genes (DEGs) using the Limma package in R, discerned their biological functions and pathways using enrichment analysis, and conducted the modular analysis based on protein-protein interaction network to identify hub genes of RB. Survival analyses based on The Cancer Genome Atlas clinical database were performed to analyze prognostic values of key genes of RB. Subsequently, we identified the differentially methylated genes, differentially expressed miRNAs (DEMs) and lncRNAs (DELs), and intersected them with key genes to analyze possible targets of the underlying epigenetic regulatory mechanisms. Finally, the ceRNA network of lncRNAs-miRNAs-mRNAs was constructed using Cytoscape.A total of 193 DEGs, 74 differentially methylated-DEGs (DM-DEGs), 45 DEMs, 5 DELs were identified. The molecular pathways of DEGs were enriched in cell cycle, p53 signaling pathway, and DNA replication. A total of 10 key genes were identified and found significantly associated with poor survival outcome based on survival analyses, including CDK1, BUB1, CCNB2, TOP2A, CCNB1, RRM2, KIF11, KIF20A, NDC80, and TTK. We further found that hub genes MCM6 and KIF14 were differentially methylated, key gene RRM2 was targeted by DEMs, and key genes TTK, RRM2, and CDK1 were indirectly regulated by DELs. Additionally, the ceRNA network with 222 regulatory associations was constructed to visualize the correlations between lncRNAs-miRNAs-mRNAs.This study presents an integrated bioinformatics analysis of genetic and epigenetic changes that may be associated with the development of RB. Findings may yield many new insights into the molecular biomarker candidates and epigenetically regulatory targets of RB.

Organoids and organ chips in ophthalmology

Recent advances have driven the development of stem cell-derived, self-organizing, three-dimensional miniature organs, termed organoids, which mimic different eye tissues including the retina, cornea, and lens. Organoids and engineered microfluidic organ-on-chips (organ chips) are transformative technologies that show promise in simulating the architectural and functional complexity of native organs. Accordingly, they enable exploration of facets of human disease and development not accurately recapitulated by animal models. Together, these technologies will increase our understanding of the basic physiology of different eye structures, enable us to interrogate unknown aspects of ophthalmic disease pathogenesis, and serve as clinically-relevant surrogates for the evaluation of ocular therapeutics. Both the burden and prevalence of monogenic and multifactorial ophthalmic diseases, which can cause visual impairment or blindness, in the human population warrants a paradigm shift towards organoids and organ chips that can provide sensitive, quantitative, and scalable phenotypic assays. In this article, we review the current situation of organoids and organ chips in ophthalmology and discuss how they can be leveraged for translational applications.

On the Generation and Regeneration of Retinal Ganglion Cells

Retinal development follows a conserved neurogenic program in vertebrates to orchestrate the generation of specific cell types from multipotent progenitors in sequential but overlapping waves. In this program, retinal ganglion cells (RGCs) are the first cell type generated. RGCs are the final output neurons of the retina and are essential for vision and circadian rhythm. Key molecular steps have been defined in multiple vertebrate species to regulate competence, specification, and terminal differentiation of this cell type. This involves neuronal-specific transcription factor networks, regulators of chromatin dynamics and miRNAs. In mammals, RGCs and their optic nerve axons undergo neurodegeneration and loss in glaucoma and other optic neuropathies, resulting in irreversible vision loss. The incapacity of RGCs and axons to regenerate reinforces the need for the design of efficient RGC replacement strategies. Here we describe the essential molecular pathways for the differentiation of RGCs in vertebrates, as well as experimental manipulations that extend the competence window for generation of this early cell type from late progenitors. We discuss recent advances in regeneration of retinal neurons in vivo in both mouse and zebrafish and discuss possible strategies and barriers to achieving RGC regeneration as a therapeutic approach for vision restoration in blinding diseases such as glaucoma.

Core transcriptional regulatory circuitries in cancer

Transcription factors (TFs) coordinate the on-and-off states of gene expression typically in a combinatorial fashion. Studies from embryonic stem cells and other cell types have revealed that a clique of self-regulated core TFs control cell identity and cell state. These core TFs form interconnected feed-forward transcriptional loops to establish and reinforce the cell-type-specific gene-expression program; the ensemble of core TFs and their regulatory loops constitutes core transcriptional regulatory circuitry (CRC). Here, we summarize recent progress in computational reconstitution and biologic exploration of CRCs across various human malignancies, and consolidate the strategy and methodology for CRC discovery. We also discuss the genetic basis and therapeutic vulnerability of CRC, and highlight new frontiers and future efforts for the study of CRC in cancer. Knowledge of CRC in cancer is fundamental to understanding cancer-specific transcriptional addiction, and should provide important insight to both pathobiology and therapeutics.

Hyperacetylated chromatin domains mark cell type-specific genes and suggest distinct modes of enhancer function

Stratification of enhancers by signal strength in ChIP-seq assays has resulted in the establishment of super-enhancers as a widespread and useful tool for identifying cell type-specific, highly expressed genes and associated pathways. We examine a distinct method of stratification that focuses on peak breadth, termed hyperacetylated chromatin domains (HCDs), which classifies broad regions exhibiting histone modifications associated with gene activation. We find that this analysis serves to identify genes that are both more highly expressed and more closely aligned to cell identity than super-enhancer analysis does using multiple data sets. Moreover, genetic manipulations of selected gene loci suggest that some enhancers located within HCDs work at least in part via a distinct mechanism involving the modulation of histone modifications across domains and that this activity can be imported into a heterologous gene locus. In addition, such genetic dissection reveals that the super-enhancer concept can obscure important functions of constituent elements.

Super-enhancer are usually defined by high levels of chromatin modification and associate with cell-specific gene expression. Here, the authors define hyperacetylated chromatin domains (HCDs) by using histone hyperacetylation peak breadth information and show that HCDs associated more closely with cell identity than super-enhancers.

DNA demethylation is a driver for chick retina regeneration

Cellular reprogramming resets the epigenetic landscape to drive shifts in transcriptional programmes and cell identity. The embryonic chick can regenerate a complete neural retina, after retinectomy, via retinal pigment epithelium (RPE) reprogramming in the presence of FGF2. In this study, we systematically analysed the reprogramming competent chick RPE prior to injury, and during different stages of reprogramming. In addition to changes in the expression of genes associated with epigenetic modifications during RPE reprogramming, we observed dynamic changes in histone marks associated with bivalent chromatin (H3K27me3/H3K4me3) and intermediates of the process of DNA demethylation including 5hmC and 5caC. Comprehensive analysis of the methylome by whole-genome bisulphite sequencing (WGBS) confirmed extensive rearrangements of DNA methylation patterns including differentially methylated regions (DMRs) found at promoters of genes associated with chromatin organization and fibroblast growth factor production. We also identified Tet methylcytosine dioxygenase 3 (TET3) as an important factor for DNA demethylation and retina regeneration, capable of reprogramming RPE in the absence of exogenous FGF2. In conclusion, we demonstrate that injury early in RPE reprogramming triggers genome-wide dynamic changes in chromatin, including bivalent chromatin and DNA methylation. In the presence of FGF2, these dynamic modifications are further sustained in the commitment to form a new retina. Our findings reveal active DNA demethylation as an important process that may be applied to remove the epigenetic barriers in order to regenerate retina in mammals.

ABBREVIATIONS

bp: Base pair; DMR: Differentially methylated region; DMC: Differentially methylated cytosines; GFP: Green fluorescent protein; PCR: Polymerase chain reaction. TET: Ten-eleven translocation; RPE: retinal pigment epithelium.

Jmjd3 Plays Pivotal Roles in the Proper Development of Early-Born Retinal Lineages: Amacrine, Horizontal, and Retinal Ganglion Cells

Purpose

Trimethylation of histone H3 at lysine 27 (H3K27me3) is a critical mediator of transcriptional gene repression, and Jmjd3 and Utx are the demethylases specific to H3K27me3. Using an in vitro retinal explant culture system, we previously revealed the role of Jmjd3 in the development of rod bipolar cells; however, the roles of Jmjd3 in the development of early-born retinal cells are unknown due to limitations concerning the use of retinal explant culture systems. In this study, we investigated the roles of Jmjd3 in the development of early-born retinal cells.

Methods

We examined retina-specific conditional Jmjd3 knockout (Jmjd3-cKO) mice using immunohistochemistry and quantitative reverse transcription PCR and JMJD3 binding to a target locus by chromatin immunoprecipitation analysis.

Results

We observed reductions in amacrine cells (ACs) and horizontal cells (HCs), as well as lowered expression levels of several transcription factors involved in the development of ACs and HCs in the Jmjd3-cKO mouse retina. JMJD3 bound the promoter regions of these transcription factors. Notably, an elevated number of retinal ganglion cells (RGCs) was observed at embryonic stages, whereas RGCs were moderately reduced at later postnatal stages in the Jmjd3-cKO retina. We also observed reduced expression of Eomes, which is required for the maintenance of RGCs, as well as lower H3K27me3 level and lower JMJD3 binding in the promoter region of Eomes in RGC-enriched cells.

Conclusions

The results indicated that Jmjd3 has critical roles in the development of early-born retinal subtypes, and suggested biphasic roles of Jmjd3 in RGC production and maintenance.

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.

Dystrophin Is Required for the Proper Timing in Retinal Histogenesis: A Thorough Investigation on the mdx Mouse Model of Duchenne Muscular Dystrophy

Duchenne muscular dystrophy (DMD) is a lethal X-linked muscular disease caused by defective expression of the cytoskeletal protein dystrophin (Dp427). Selected autonomic and central neurons, including retinal neurons, express Dp427 and/or dystrophin shorter isoforms. Because of this, DMD patients may also experience different forms of cognitive impairment, neurological and autonomic disorders, and specific visual defects. DMD-related damages to the nervous system are established during development, suggesting a role for all dystrophin isoforms in neural circuit development and differentiation; however, to date, their function in retinogenesis has never been investigated. In this large-scale study, we analyzed whether the lack of Dp427 affects late retinogenesis in the mdx mouse, the most well studied animal model of DMD. Retinal gene expression and layer maturation, as well as neural cell proliferation, apoptosis, and differentiation, were evaluated in E18 and/or P0, P5, P10, and adult mice. In mdx mice, expression of Capn3, Id3 (E18-P5), and Dtnb (P5) genes, encoding proteins involved in different aspects of retina development and synaptogenesis (e.g., Calpain 3, DNA-binding protein inhibitor-3, and β-dystrobrevin, respectively), was transiently reduced compared to age-matched wild type mice. Concomitantly, a difference in the time required for the retinal ganglion cell layer to reach appropriate thickness was observed (P0–P5). Immunolabeling for specific cell markers also evidenced a significant dysregulation in the number of GABAergic amacrine cells (P5–P10), a transient decrease in the area immunopositive for the Vesicular Glutamate Transporter 1 (VGluT1) during ribbon synapse maturation (P10) and a reduction in the number of calretinin⁺ retinal ganglion cells (RGCs) (adults). Finally, the number of proliferating retinal progenitor cells (P5–P10) and apoptotic cells (P10) was reduced. These results support the hypothesis of a role for Dp427 during late retinogenesis different from those proposed in consolidated neural circuits. In particular, Dp427 may be involved in shaping specific steps of retina differentiation. Notably, although most of the above described quantitative alterations recover over time, the number of calretinin⁺ RGCs is reduced only in the mature retina. This suggests that alterations subtler than the timing of retinal maturation may occur, a hypothesis that demands further in-depth functional studies.

Chromatin regulators in retinoblastoma: Biological roles and therapeutic applications

Retinoblastoma (RB) is a pediatric ocular tumor mostly occurring due to the biallelic loss of RB1 gene in the developing retina. Early studies of genomic aberrations in RB have provided a valuable insight into how RB can progress following the tumor-initiating RB1 mutations and have established a notion that inactivation of RB1 gene is critical to initiate RB but this causative genetic lesion alone is not sufficient for malignant progression. With the advent of high-throughput sequencing technologies, we now have access to the comprehensive genomic and epigenetic landscape of RB and have come to appreciate that RB tumorigenesis requires both genetic and epigenetic alterations that might be directly or indirectly driven by RB1 loss. This integrative perspective on RB tumorigenesis has inspired research efforts to better understand the types and functions of epigenetic mechanisms contributing to RB development, leading to the identification of multiple epigenetic regulators misregulated in RB in recent years. A complete understanding of the intricate network of genetic and epigenetic factors in modulation of gene expression during RB tumorigenesis remains a major challenge but would be crucial to translate these findings into therapeutic interventions. In this review, we will provide an overview of chromatin regulators identified to be misregulated in human RB among the numerous epigenetic factors implicated in RB development. For a subset of these chromatin regulators, recent findings on their functions in RB development and potential therapeutic applications are discussed.

The Use of Induced Pluripotent Stem Cells as a Model for Developmental Eye Disorders

Approximately one-third of childhood blindness is attributed to developmental eye disorders, of which 80% have a genetic cause. Eye morphogenesis is tightly regulated by a highly conserved network of transcription factors when disrupted by genetic mutations can result in severe ocular malformation. Human-induced pluripotent stem cells (hiPSCs) are an attractive tool to study early eye development as they are more physiologically relevant than animal models, can be patient-specific and their use does not elicit the ethical concerns associated with human embryonic stem cells. The generation of self-organizing hiPSC-derived optic cups is a major advancement to understanding mechanisms of ocular development and disease. Their development in vitro has been found to mirror that of the human eye and these early organoids have been used to effectively model microphthalmia caused by a VSX2 variant. hiPSC-derived optic cups, retina, and cornea organoids are powerful tools for future modeling of disease phenotypes and will enable a greater understanding of the pathophysiology of many other developmental eye disorders. These models will also provide an effective platform for identifying molecular therapeutic targets and for future clinical applications.

Ybx1 fine-tunes PRC2 activities to control embryonic brain development

Chromatin modifiers affect spatiotemporal gene expression programs that underlie organismal development. The Polycomb repressive complex 2 (PRC2) is a crucial chromatin modifier in executing neurodevelopmental programs. Here, we find that PRC2 interacts with the nucleic acid-binding protein Ybx1. In the mouse embryo in vivo, Ybx1 is required for forebrain specification and restricting mid-hindbrain growth. In neural progenitor cells (NPCs), Ybx1 controls self-renewal and neuronal differentiation. Mechanistically, Ybx1 highly overlaps PRC2 binding genome-wide, controls PRC2 distribution, and inhibits H3K27me3 levels. These functions are consistent with Ybx1-mediated promotion of genes involved in forebrain specification, cell proliferation, or neuronal differentiation. In Ybx1-knockout NPCs, H3K27me3 reduction by PRC2 enzymatic inhibitor or genetic depletion partially rescues gene expression and NPC functions. Our findings suggest that Ybx1 fine-tunes PRC2 activities to regulate spatiotemporal gene expression in embryonic neural development and uncover a crucial epigenetic mechanism balancing forebrain-hindbrain lineages and self-renewal-differentiation choices in NPCs.

Core-clock genes Period 1 and 2 regulate visual cascade and cell cycle components during mouse eye development

The retinas from Period 1 (Per1) and Period 2 (Per2) double-mutant mice (Per1^-/-Per2^Brdm1) display abnormal blue-cone distribution associated with a reduction in cone opsin mRNA and protein levels, up to 1 year of age. To reveal the molecular mechanisms by which Per1 and Per2 control retina development, we analyzed genome-wide gene expression differences between wild-type (WT) and Per1^-/-Per2^Brdm1 mice across ocular developmental stages (E15, E18 and P3). All clock genes displayed changes in transcript levels along with normal eye development. RNA-Seq data show major gene expression changes between WT and mutant eyes, with the number of differentially expressed genes (DEG) increasing with developmental age. Functional annotation of the genes showed that the most significant changes in expression levels in mutant mice involve molecular pathways relating to circadian rhythm signaling at E15 and E18. At P3, the visual cascade and the cell cycle were respectively higher and lower expressed compared to WT eyes. Overall, our study provides new insights into signaling pathways -phototransduction and cell cycle- controlled by the circadian clock in the eye during development.

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.

Development of Stem Cell Therapies for Retinal Degeneration

Degenerative retinal disease is the major cause of sight loss in the developed world and currently there is a lack of effective treatments. As the loss of vision is directly the result of the loss of retinal cells, effective cell replacement through stem-cell-based therapies may have the potential to treat a great number of retinal diseases whatever their underlying etiology. The eye is an ideal organ to develop cell therapies as it is immune privileged, and modern surgical techniques enable precise delivery of cells to the retina. Furthermore, a range of noninvasive diagnostic tests and high-resolution imaging techniques facilitate the evaluation of any therapeutic intervention. In this review, we evaluate the progress to date of current cell therapy strategies for retinal repair, focusing on transplantation of pluripotent stem-cell-derived retinal pigment epithelium (RPE) and photoreceptor cells.

UTX Regulates Human Neural Differentiation and Dendritic Morphology by Resolving Bivalent Promoters

UTX, a H3K27me3 demethylase, plays an important role in mouse brain development. However, so little is known about the function of UTX in human neural differentiation and dendritic morphology. In this study, we generated UTX-null human embryonic stem cells using CRISPR/Cas9, and differentiated them into neural progenitor cells and neurons to investigate the effects of UTX loss of function on human neural development. The results showed that the number of differentiated neurons significantly reduced after loss of UTX, and that the dendritic morphology of UTX KO neurons tended to be simplified. The electrophysiological recordings showed that most of the UTX KO neurons were immature. Finally, RNA sequencing identified dozens of differentially expressed genes involved in neural differentiation and synaptic function in UTX KO neurons and our results demonstrated that UTX regulated these critical genes by resolving bivalent promoters. In summary, we establish a reference for the important role of UTX in human neural differentiation and dendritic morphology.

•

Loss of UTX in hESCs reduces their neural differentiation potential

•

The dendritic morphology of UTX KO neurons tends to be simplified

•

UTX regulates human neural development depending on its demethylation

•

UTX regulates the expression of genes by resolving bivalent promoters

dnmt1 function is required to maintain retinal stem cells within the ciliary marginal zone of the zebrafish eye

The ciliary marginal zone (CMZ) of the zebrafish retina contains a population of actively proliferating resident stem cells, which generate retinal neurons throughout life. The maintenance methyltransferase, dnmt1, is expressed within the CMZ. Loss of dnmt1 function results in gene misregulation and cell death in a variety of developmental contexts, however, its role in retinal stem cell (RSC) maintenance is currently unknown. Here, we demonstrate that zebrafish dnmt1^s872 mutants possess severe defects in RSC maintenance within the CMZ. Using a combination of immunohistochemistry, in situ hybridization, and a transgenic reporter assay, our results demonstrate a requirement for dnmt1 activity in the regulation of RSC proliferation, gene expression and in the repression of endogenous retroelements (REs). Ultimately, cell death is elevated in the dnmt1^−/− CMZ, but in a p53-independent manner. Using a transgenic reporter for RE transposition activity, we demonstrate increased transposition in the dnmt1^−/− CMZ. Taken together our data identify a critical role for dnmt1 function in RSC maintenance in the vertebrate eye.

Comprehensive identification of mRNA isoforms reveals the diversity of neural cell-surface molecules with roles in retinal development and disease

Genes encoding cell-surface proteins control nervous system development and are implicated in neurological disorders. These genes produce alternative mRNA isoforms which remain poorly characterized, impeding understanding of how disease-associated mutations cause pathology. Here we introduce a strategy to define complete portfolios of full-length isoforms encoded by individual genes. Applying this approach to neural cell-surface molecules, we identify thousands of unannotated isoforms expressed in retina and brain. By mass spectrometry we confirm expression of newly-discovered proteins on the cell surface in vivo. Remarkably, we discover that the major isoform of a retinal degeneration gene, CRB1, was previously overlooked. This CRB1 isoform is the only one expressed by photoreceptors, the affected cells in CRB1 disease. Using mouse mutants, we identify a function for this isoform at photoreceptor-glial junctions and demonstrate that loss of this isoform accelerates photoreceptor death. Therefore, our isoform identification strategy enables discovery of new gene functions relevant to disease.

Transcriptomic Analysis of the Developmental Similarities and Differences Between the Native Retina and Retinal Organoids

Purpose

We performed a bioinformatic transcriptome analysis to determine the alteration of gene expression between the native retina and retinal organoids in both mice and humans.

Methods

The datasets of mouse native retina (GSE101986), mouse retinal organoids (GSE102794), human native retina (GSE104827), and human retinal organoids (GSE119320) were obtained from Gene Expression Omnibus. After normalization, a principal component analysis was performed to categorize the samples. The genes were clustered to classify them. A functional analysis was performed using the bioinformatics tool Gene ontology enrichment to analyze the biological processes of selected genes and cellular components.

Results

The development of retinal organoids is slower than that in the native retina. In the early stage, cell proliferation predominates. Subsequently, neural differentiation is dominant. In the later stage, the dominant differentiated cells are photoreceptors. Additionally, the fatty acid metabolic process and mitochondria-related genes are upregulated over time, and the glycogen catabolic process and activin receptors are gradually downregulated in human retinal organoids. However, these trends are opposite in mouse retinal organoids. There are two peaks in mitochondria-related genes, one in the early development period and another during the photoreceptor development period. It takes about five times longer for human retinal development to achieve similar levels of mouse retinal development.

Conclusions

Our study reveals the similarities and differences in the developmental features of retinal organoids as well as the corresponding relationship between mouse and human retinal development.

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.

Genome-wide Profiling Identifies DNA Methylation Signatures of Aging in Rod Photoreceptors Associated with Alterations in Energy Metabolism

Aging-associated functional decline is accompanied by alterations in the epigenome. To explore DNA modifications that could influence visual function with age, we perform whole-genome bisulfite sequencing of purified mouse rod photoreceptors at four ages and identify 2,054 differentially methylated regions (DMRs). We detect many DMRs during early stages of aging and in rod regulatory regions, and some of these cluster at chromosomal hotspots, especially on chromosome 10, which includes a longevity interactome. Integration of methylome to age-related transcriptome changes, chromatin signatures, and first-order protein-protein interactions uncover an enrichment of DMRs in altered pathways that are associated with rod function, aging, and energy metabolism. In concordance, we detect reduced basal mitochondrial respiration and increased fatty acid dependency with retinal age in ex vivo assays. Our study reveals age-dependent genomic and chromatin features susceptible to DNA methylation changes in rod photoreceptors and identifies a link between DNA methylation and energy metabolism in aging.

Clinical and genetical features of probands and affected family members with familial exudative vitreoretinopathy in a large Chinese cohort

Aims

To explore the clinical and genetical features of families with strictly confirmed familial exudative vitreoretinopathy (FEVR) in a large Chinese cohort.

Methods

A retrospective chart review study was conducted on the FEVR families diagnosed by both angiography and targeted next-generation sequencing in six FEVR known genes (FZD4, LRP5, TSPAN12, NDP, KIF11, ZNF408) in the probands and at least one first-degree family member. Variation in expressivity and severity was evaluated in different gene groups.

Results

105 FEVR families (223 FEVR affected subjects with 434 eyes) met the inclusion criteria. There were 105 probands with mean age of 3.8 years old and 118 affected family members of 32.7 years old averagely. Mutations in FZD4 were most prevalent (33.33%), followed by LRP5 (29.52%), TSPAN12 (22.86%), NDP (5.71%), KIF11 (1.9%) and ZNF408 (0.95%). 81% of the probands were classified as stage 4 or worse which most prevalently contributed to FZD4 mutations. All of the three affected family members with stage 4 or worse carried FZD4 variants. More than half (51.43%) of the probands in FZD4 group showed asymmetry. Unilateral FEVR was detected in 11 (10.5%) families consisting of six probands and six affected relatives, and FZD4 mutations accounted for 63.64% of all the cases with variant (c.1282_1285del, p. D428fs) identified in three families.

Conclusions

Genotype-phenotype correlation in FEVR was complex with family dependent. Mutations in FZD4 might initiate the most diverse and asymmetric phenotypes.

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.

Chromatin remodeling protein HELLS is critical for retinoblastoma tumor initiation and progression

Retinoblastoma is an aggressive childhood cancer of the developing retina that initiates by biallelic RB1 gene inactivation. Tumor progression in retinoblastoma is driven by epigenetics, as retinoblastoma genomes are stable, but the mechanism(s) that drive these epigenetic changes remain unknown. Lymphoid-specific helicase (HELLS) protein is an epigenetic modifier directly regulated by the RB/E2F pathway. In this study, we used novel genetically engineered mouse models to investigate the role of HELLS during retinal development and tumorigenesis. Our results indicate that Hells-null retinal progenitor cells divide, undergo cell-fate specification, and give rise to fully laminated retinae with minor bipolar cells defects, but normal retinal function. Despite the apparent nonessential role of HELLS in retinal development, failure to transcriptionally repress Hells during retinal terminal differentiation due to retinoblastoma (RB) family loss significantly contributes to retinal tumorigenesis. Loss of HELLS drastically reduced ectopic division of differentiating cells in Rb1/p107-null retinae, significantly decreased the incidence of retinoblastoma, delayed tumor progression, and increased overall survival. Despite its role in heterochromatin formation, we found no evidence that Hells loss directly affected chromatin accessibility in the retina but functioned as transcriptional co-activator of E2F3, decreasing expression of cell cycle genes. We propose that HELLS is a critical downstream mediator of E2F-dependent ectopic proliferation in RB-null retinae. Together with the nontoxic effect of HELLS loss in the developing retina, our results suggest that HELLS and its downstream pathways could serve as potential therapeutic targets for retinoblastoma.

MYCN amplification and ATRX mutations are incompatible in neuroblastoma

Aggressive cancers often have activating mutations in growth-controlling oncogenes and inactivating mutations in tumor-suppressor genes. In neuroblastoma, amplification of the MYCN oncogene and inactivation of the ATRX tumor-suppressor gene correlate with high-risk disease and poor prognosis. Here we show that ATRX mutations and MYCN amplification are mutually exclusive across all ages and stages in neuroblastoma. Using human cell lines and mouse models, we found that elevated MYCN expression and ATRX mutations are incompatible. Elevated MYCN levels promote metabolic reprogramming, mitochondrial dysfunction, reactive-oxygen species generation, and DNA-replicative stress. The combination of replicative stress caused by defects in the ATRX-histone chaperone complex, and that induced by MYCN-mediated metabolic reprogramming, leads to synthetic lethality. Therefore, ATRX and MYCN represent an unusual example, where inactivation of a tumor-suppressor gene and activation of an oncogene are incompatible. This synthetic lethality may eventually be exploited to improve outcomes for patients with high-risk neuroblastoma.

Chromatin accessibility analysis reveals regulatory dynamics of developing human retina and hiPSC-derived retinal organoids

Retinal organoids (ROs) derived from human induced pluripotent stem cells (hiPSCs) provide potential opportunities for studying human retinal development and disorders; however, to what extent ROs recapitulate the epigenetic features of human retinal development is unknown. In this study, we systematically profiled chromatin accessibility and transcriptional dynamics over long-term human retinal and RO development. Our results showed that ROs recapitulated the human retinogenesis to a great extent, but divergent chromatin features were also discovered. We further reconstructed the transcriptional regulatory network governing human and RO retinogenesis in vivo. Notably, NFIB and THRA were identified as regulators in human retinal development. The chromatin modifications between developing human and mouse retina were also cross-analyzed. Notably, we revealed an enriched bivalent modification of H3K4me3 and H3K27me3 in human but not in murine retinogenesis, suggesting a more dedicated epigenetic regulation on human genome.

Retina Development in Vertebrates: Systems Biology Approaches to Understanding Genetic Programs: On the Contribution of Next-Generation Sequencing Methods to the Characterization of the Regulatory Networks Controlling Vertebrate Eye Development

The ontogeny of the vertebrate retina has been a topic of interest to developmental biologists and human geneticists for many decades. Understanding the unfolding of the genetic program that transforms a field of progenitors cells into a functionally complex and multi-layered sensory organ is a formidable challenge. Although classical genetic studies succeeded in identifying the key regulators of retina specification, understanding the architecture of their gene network and predicting their behavior are still a distant hope. The emergence of next-generation sequencing platforms revolutionized the field unlocking the access to genome-wide datasets. Emerging techniques such as RNA-seq, ChIP-seq, ATAC-seq, or single cell RNA-seq are used to characterize eye developmental programs. These studies provide valuable information on the transcriptional and cis-regulatory profiles of precursors and differentiated cells, outlining the trajectories that connect each intermediate state. Here, recent systems biology efforts are reviewed to understand the genetic programs shaping the vertebrate retina.

Pou2f1 and Pou2f2 cooperate to control the timing of cone photoreceptor production in the developing mouse retina

Multipotent retinal progenitor cells (RPCs) generate various cell types in a precise chronological order, but how exactly cone photoreceptor production is restricted to early stages remains unclear. Here, we show that the POU-homeodomain factors Pou2f1/Pou2f2, the homologs of Drosophila temporal identity factors nub/pdm2, regulate the timely production of cones in mice. Forcing sustained expression of Pou2f1 or Pou2f2 in RPCs expands the period of cone production, whereas misexpression in late-stage RPCs triggers ectopic cone production at the expense of late-born fates. Mechanistically, we report that Pou2f1 induces Pou2f2 expression, which binds to a POU motif in the promoter of the rod-inducing factor Nrl to repress its expression. Conversely, conditional inactivation of Pou2f2 in RPCs increases Nrl expression and reduces cone production. Finally, we provide evidence that Pou2f1 is part of a cross-regulatory cascade with the other temporal identity factors Ikzf1 and Casz1. These results uncover Pou2f1/2 as regulators of the temporal window for cone genesis and, given their widespread expression in the nervous system, raise the possibility of a general role in temporal patterning.

Muscle-derived Dpp regulates feeding initiation via endocrine modulation of brain dopamine biosynthesis

In animals, the brain regulates feeding behavior in response to local energy demands of peripheral tissues, which secrete orexigenic and anorexigenic hormones. Although skeletal muscle is a key peripheral tissue, it remains unknown whether muscle-secreted hormones regulate feeding. In Drosophila, we found that decapentaplegic (dpp), the homolog of human bone morphogenetic proteins BMP2 and BMP4, is a muscle-secreted factor (a myokine) that is induced by nutrient sensing and that circulates and signals to the brain. Muscle-restricted dpp RNAi promotes foraging and feeding initiation, whereas dpp overexpression reduces it. This regulation of feeding by muscle-derived Dpp stems from modulation of brain tyrosine hydroxylase (TH) expression and dopamine biosynthesis. Consistently, Dpp receptor signaling in dopaminergic neurons regulates TH expression and feeding initiation via the downstream transcriptional repressor Schnurri. Moreover, pharmacologic modulation of TH activity rescues the changes in feeding initiation due to modulation of dpp expression in muscle. These findings indicate that muscle-to-brain endocrine signaling mediated by the myokine Dpp regulates feeding behavior.

Identification of Genes With Enriched Expression in Early Developing Mouse Cone Photoreceptors

Purpose

The early transcriptional events that occur in newly generated cone photoreceptors are not well described. Knowledge of these events is critical to provide benchmarks for in vitro-derived cone photoreceptors and to understand the process of cone and rod photoreceptor diversification. We sought to identify genes with differential gene expression in embryonic mouse cone photoreceptors.

Methods

The specificity of expression of the LHX4 transcription factor in developing cone photoreceptors was examined using immunofluorescence visualization in both mouse and chicken retinas. A LHX4 transgenic reporter line with high specificity for developing mouse cone photoreceptors was identified and used to purify early-stage cone photoreceptors for profiling by single-cell RNA sequencing. Comparisons were made to previous datasets targeting photoreceptors.

Results

The LHX4 transcription factor and a transgenic reporter were determined to be highly specific to early developing cone photoreceptors in the mouse. Single-cell transcriptional profiling identified new genes with enriched expression in cone photoreceptors relative to concurrent cell populations. Comparison to previous profiling datasets allowed for further characterization of these genes across developmental time, species, photoreceptor type, and gene regulatory network.

Conclusions

The LHX4 gene is highly enriched in developing cone photoreceptors as are several new genes identified through transcriptional profiling, some of which are expressed in subclusters of cones. Many of these cone-enriched genes do not show obvious de-repression in profiling of retinas mutant for the rod-specific transcription factor NRL, highlighting differences between endogenous cones and those induced in NRL mutants.

Eye in a Disk: eyeIntegration Human Pan-Eye and Body Transcriptome Database Version 1.0

Purpose

We develop an accessible and reliable RNA sequencing (RNA-seq) transcriptome database of healthy human eye tissues and a matching reactive web application to query gene expression in eye and body tissues.

Methods

We downloaded the raw sequence data for 1375 RNA-seq samples across 54 tissues in the Genotype-Tissue Expression (GTEx) project as a noneye reference set. We then queried several public repositories to find all healthy, nonperturbed, human eye-related tissue RNA-seq samples. The 916 eye and 1375 GTEx samples were sent into a Snakemake-based reproducible pipeline we wrote to quantify all known transcripts and genes, removes samples with poor sequence quality and mislabels, normalizes expression values across each tissue, perform 882 differential expression tests, calculate GO term enrichment, and output all as a single SQLite database file: the Eye in a Disk (EiaD) dataset. Furthermore, we rewrote the web application eyeIntegration (available in the public domain at https://eyeIntegration.nei.nih.gov) to display EiaD.

Results

The new eyeIntegration portal provides quick visualization of human eye-related transcriptomes published to date by database version, gene/transcript, 19 eye tissues, and 54 body tissues. As a test of the value of this unified pan-eye dataset, we showed that fetal and organoid retina are highly similar at a pan-transcriptome level, but display distinct differences in certain pathways and gene families, such as protocadherin and HOXB family members.

Conclusions

The eyeIntegration v1.0 web app serves the pan-human eye and body transcriptome dataset, EiaD. This offers the eye community a powerful and quick means to test hypotheses on human gene and transcript expression across 54 body and 19 eye tissues.

Genetic and epigenetic control of retinal development in zebrafish

The vertebrate retina is a complex structure composed of seven cell types (six neuron and one glia), and all of which originate from a seemingly homogeneous population of proliferative multipotent retinal progenitor cells (RPCs) that exit the cell cycle and differentiate in a spatio-temporally regulated and stereotyped fashion. This neurogenesis process requires intricate genetic regulation involving a combination of cell intrinsic transcription factors and extrinsic signaling molecules, and many critical factors have been identified that influence the timing and composition of the developing retina. Adding complexity to the process, over the past decade, a variety of epigenetic regulatory mechanisms have been shown to influence neurogenesis, and these include changes in histone modifications and the chromatin landscape and changes in DNA methylation and hydroxymethylation patterns. This review summarizes recent findings in the genetic and epigenetic regulation of retinal development, with an emphasis on the zebrafish model system, and it outlines future areas of investigation that will continue to push the field forward into the epigenomics era.

Retinogenesis of the Human Fetal Retina: An Apical Polarity Perspective

The Crumbs complex has prominent roles in the control of apical cell polarity, in the coupling of cell density sensing to downstream cell signaling pathways, and in regulating junctional structures and cell adhesion. The Crumbs complex acts as a conductor orchestrating multiple downstream signaling pathways in epithelial and neuronal tissue development. These pathways lead to the regulation of cell size, cell fate, cell self-renewal, proliferation, differentiation, migration, mitosis, and apoptosis. In retinogenesis, these are all pivotal processes with important roles for the Crumbs complex to maintain proper spatiotemporal cell processes. Loss of Crumbs function in the retina results in loss of the stratified appearance resulting in retinal degeneration and loss of visual function. In this review, we begin by discussing the physiology of vision. We continue by outlining the processes of retinogenesis and how well this is recapitulated between the human fetal retina and human embryonic stem cell (ESC) or induced pluripotent stem cell (iPSC)-derived retinal organoids. Additionally, we discuss the functionality of in utero and preterm human fetal retina and the current level of functionality as detected in human stem cell-derived organoids. We discuss the roles of apical-basal cell polarity in retinogenesis with a focus on Leber congenital amaurosis which leads to blindness shortly after birth. Finally, we discuss Crumbs homolog (CRB)-based gene augmentation.

Improved Retinal Organoid Differentiation by Modulating Signaling Pathways Revealed by Comparative Transcriptome Analyses with Development In Vivo

Stem cell-derived retinal organoids recapitulate many landmarks of in vivo differentiation but lack functional maturation of distinct cell types, especially photoreceptors. Using comprehensive temporal transcriptome analyses, we show that transcriptome shift from postnatal day 6 (P6) to P10, associated with morphogenesis and synapse formation during mouse retina development, was not evident in organoids, and co-expression clusters with similar patterns included different sets of genes. Furthermore, network analysis identified divergent regulatory dynamics between developing retina in vivo and in organoids, with temporal dysregulation of specific signaling pathways and delayed or reduced expression of genes involved in photoreceptor function(s) and survival. Accordingly, addition of docosahexaenoic acid and fibroblast growth factor 1 to organoid cultures specifically promoted the maturation of photoreceptors, including cones. Our study thus identifies regulatory signals deficient in developing retinal organoids and provides experimental validation by producing a more mature retina in vitro, thereby facilitating investigations in disease modeling and therapies.

ATRX In-Frame Fusion Neuroblastoma Is Sensitive to EZH2 Inhibition via Modulation of Neuronal Gene Signatures

ATRX alterations occur at high frequency in neuroblastoma of adolescents and young adults. Particularly intriguing are the large N-terminal deletions of ATRX (Alpha Thalassemia/Mental Retardation, X-linked) that generate in-frame fusion (IFF) proteins devoid of key chromatin interaction domains, while retaining the SWI/SNF-like helicase region. We demonstrate that ATRX IFF proteins are redistributed from H3K9me3-enriched chromatin to promoters of active genes and identify REST as an ATRX IFF target whose activation promotes silencing of neuronal differentiation genes. We further show that ATRX IFF cells display sensitivity to EZH2 inhibitors, due to derepression of neurogenesis genes, including a subset of REST targets. Taken together, we demonstrate that ATRX structural alterations are not loss-of-function and put forward EZH2 inhibitors as a potential therapy for ATRX IFF neuroblastoma.

Nucleome Dynamics during Retinal Development

More than 8,000 genes are turned on or off as progenitor cells produce the 7 classes of retinal cell types during development. Thousands of enhancers are also active in the developing retinae, many having features of cell- and developmental stage-specific activity. We studied dynamic changes in the 3D chromatin landscape important for precisely orchestrated changes in gene expression during retinal development by ultra-deep in situ Hi-C analysis on murine retinae. We identified developmental-stage-specific changes in chromatin compartments and enhancer-promoter interactions. We developed a machine learning-based algorithm to map euchromatin and heterochromatin domains genome-wide and overlaid it with chromatin compartments identified by Hi-C. Single-cell ATAC-seq and RNA-seq were integrated with our Hi-C and previous ChIP-seq data to identify cell- and developmental-stage-specific super-enhancers (SEs). We identified a bipolar neuron-specific core regulatory circuit SE upstream of Vsx2, whose deletion in mice led to the loss of bipolar neurons.

Characterization of LSD1 Expression Within the Murine Eye

Purpose

The purpose of this study was to extend the current understanding of endogenous lysine-specific demethylase 1 (LSD1) expression spatially and temporally in the retina. Toward that end, we determined the localization and levels of LSD1 and its substrates H3K4me1 and H3K4me2 (H3K4me1/2) within the murine eye.

Methods

Immunofluorescent microscopy for LSD1, H3K4me1, and H3K4me2 was conducted on murine formalin-fixed paraffin-embedded eye sections across development in addition to Western immunoblotting to assess localization and protein levels.

Results

Retinal LSD1 protein levels were highest at postnatal day 7 (P7), whereas its substrates H3K4me1 and H3K4me2 had equally high levels at P2 and P14. Concentrations of all three proteins gradually decreased over developmental time until reaching a basement level of ∼60% of maximum at P36. LSD1 and H3K4me1/2 were expressed uniformly in all retinal progenitor cells. By P36, there was variability in LSD1 expression in the ganglion cell layer, uniform expression in the inner nuclear layer, and dichotomous expression between photoreceptors in the outer nuclear layer. This contrasted with H3K4me1/2 expression, which remained uniform. Additionally, LSD1 was widely expressed in the lens, cornea, and retinal pigment epithelium.

Conclusions

Consistent with its known role in neuronal differentiation, LSD1 is highly and uniformly expressed throughout all retinal progenitor cells. Variability in LSD1 expression, particularly in photoreceptors, may be indicative of their unique transcriptomes and epigenetic patterns of rods and cones. Murine rod nuclei exhibit LSD1 expression in a ring or shell, rather than throughout the nucleus, consistent with their unique inverted chromatin organization. LSD1 has substantial expression throughout adulthood, especially in cone nuclei. By providing insight into endogenous LSD1 expression, our current findings could directly inform future studies to determine the exact role of Lsd1 in the development and maintenance of specific structures and cell types within the eye.

Long noncoding RNAs sustain high expression levels of exogenous octamer-binding protein 4 by sponging regulatory microRNAs during cellular reprogramming

Long noncoding RNAs (lncRNAs) modulate gene expression as competing endogenous RNAs (ceRNAs) that sponge regulatory microRNAs (miRNAs). During cellular reprogramming, genes associated with pluripotency establishment need to be up-regulated, and developmental genes need to be silenced. However, how ceRNAs control cellular reprogramming still awaits full elucidation. Here, we used doxycycline-inducible expression of the four transcription factors octamer-binding protein 4 (OCT4), SRY-box 2 (SOX2), Krüppel-like factor 4 (KLF4), and proto-oncogene c-Myc (c-Myc) to generate induced pluripotent stem cells (iPSCs) from mouse embryonic fibroblasts (MEFs). Using RNA-Seq and bioinformatics approaches, we found that the expression levels of miRNAs from MEFs remain high from day 0 to 6 after the doxycycline induction. Many genes targeted by these miRNAs were up-regulated, and long intergenic noncoding RNAs (lincRNAs) and circular RNAs (circRNAs), which have complementary binding sites to these miRNAs, were highly expressed, indicating lincRNAs and circRNAs may function as ceRNAs. Intriguingly, knockdown of the linc/circRNAs that sponge the miRNAs, which target OCT4 down-regulated exogenous OCT4, decreased reprogramming efficiency, and resulted in low-grade iPSCs. Our results suggest that the ceRNA network plays an important role in cellular reprogramming.

Intrinsic mechanisms of neuronal axon regeneration

Permanent disabilities following CNS injuries result from the failure of injured axons to regenerate and rebuild functional connections with their original targets. By contrast, injury to peripheral nerves is followed by robust regeneration, which can lead to recovery of sensory and motor functions. This regenerative response requires the induction of widespread transcriptional and epigenetic changes in injured neurons. Considerable progress has been made in recent years in understanding how peripheral axon injury elicits these widespread changes through the coordinated actions of transcription factors, epigenetic modifiers and, to a lesser extent, microRNAs. Although many questions remain about the interplay between these mechanisms, these new findings provide important insights into the pivotal role of coordinated gene expression and chromatin remodelling in the neuronal response to injury.

dbCoRC: a database of core transcriptional regulatory circuitries modeled by H3K27ac ChIP-seq signals

Core transcription regulatory circuitry (CRC) is comprised of a small group of self-regulated transcription factors (TFs) and their interconnected regulatory loops. Studies from embryonic stem cells and other cellular models have revealed the elementary roles of CRCs in transcriptional control of cell identity and cellular fate. Systematic identification and subsequent archiving of CRCs across diverse cell types and tissues are needed to explore both cell/tissue type-specific and disease-associated transcriptional networks. Here, we present a comprehensive and interactive database (dbCoRC, http://dbcorc.cam-su.org) of CRC models which are computationally inferred from mapping of super-enhancer and prediction of TF binding sites. The current version of dbCoRC contains CRC models for 188 human and 50 murine cell lines/tissue samples. In companion with CRC models, this database also provides: (i) super enhancer, typical enhancer, and H3K27ac landscape for individual samples, (ii) putative binding sites of each core TF across the super-enhancer regions within CRC and (iii) expression of each core TF in normal or cancer cells/tissues. The dbCoRC will serve as a valuable resource for the scientific community to explore transcriptional control and regulatory circuitries in biological processes related to, but not limited to lineage specification, tissue homeostasis and tumorigenesis.

DNA Methylation Dynamics During the Differentiation of Retinal Progenitor Cells Into Retinal Neurons Reveal a Role for the DNA Demethylation Pathway

To evaluate the contribution of the DNA methylation and DNA demethylation pathways in retinal development, we studied DNA methylation in retinal progenitor cells (RPCs) and retinal neurons using a combination of whole genome bisulfite sequencing (WGBS) data obtained in our study and WGBS data collected from previous studies. The data was analyzed using Hidden Markov Model- and change point-based methods to identify methylome states in different segments of the studied genomes following genome annotation. We found that promoters of rod and cone phototransduction genes and rod photoreceptor genes, but not genes required for the development and function of other retinal phenotypes, were highly methylated in DNA isolated from human and murine fetal retinas (which mostly contain RPCs) and postnatal murine RPCs. While these highly methylated genomic regions were inherited by non-photoreceptor phenotypes during RPC differentiation, the methylation of these promoters was significantly reduced during RPC differentiation into photoreceptors and accompanied by increased expression of these genes. Our analysis of DNA methylation during embryogenesis revealed low methylation levels in genomic regions containing photoreceptor genes at the inner cell mass stage, but a sharp increase in methylation at the epiblast stage, which remained the same later on (except for DNA demethylation in photoreceptors). Thus, our data suggest that the DNA demethylation pathway is required for photoreceptor phenotypes in the developing retina. Meanwhile, the role of the DNA methylation and DNA demethylation pathways during RPC differentiation into non-photoreceptor retinal phenotypes might be less important.

Dissecting the transcriptome landscape of the human fetal neural retina and retinal pigment epithelium by single-cell RNA-seq analysis

The developmental pathway of the neural retina (NR) and retinal pigment epithelium (RPE) has been revealed by extensive research in mice. However, the molecular mechanisms underlying the development of the human NR and RPE, as well as the interactions between these two tissues, have not been well defined. Here, we analyzed 2,421 individual cells from human fetal NR and RPE using single-cell RNA sequencing (RNA-seq) technique and revealed the tightly regulated spatiotemporal gene expression network of human retinal cells. We identified major cell classes of human fetal retina and potential crucial transcription factors for each cell class. We dissected the dynamic expression patterns of visual cycle- and ligand-receptor interaction-related genes in the RPE and NR. Moreover, we provided a map of disease-related genes for human fetal retinal cells and highlighted the importance of retinal progenitor cells as potential targets of inherited retinal diseases. Our findings captured the key in vivo features of the development of the human NR and RPE and offered insightful clues for further functional studies.

Single-Cell RNA-Seq Analysis of Retinal Development Identifies NFI Factors as Regulating Mitotic Exit and Late-Born Cell Specification

Precise temporal control of gene expression in neuronal progenitors is necessary for correct regulation of neurogenesis and cell fate specification. However, the cellular heterogeneity of the developing CNS has posed a major obstacle to identifying the gene regulatory networks that control these processes. To address this, we used single-cell RNA sequencing to profile ten developmental stages encompassing the full course of retinal neurogenesis. This allowed us to comprehensively characterize changes in gene expression that occur during initiation of neurogenesis, changes in developmental competence, and specification and differentiation of each major retinal cell type. We identify the NFI transcription factors (Nfia, Nfib, and Nfix) as selectively expressed in late retinal progenitor cells and show that they control bipolar interneuron and Müller glia cell fate specification and promote proliferative quiescence.

From Retina to Stem Cell and Back Again: Memories of a Chromatin Journey

The use of retinal organoids requires efficient differentiation from induced pluripotent stem cells (iPSCs). In this issue of Cell Reports, Wang et al. (2018) examine how the chromatin landscape after iPSC programming predicts their ability to differentiate into retinal tissue.

Retinal Cell Type DNA Methylation and Histone Modifications Predict Reprogramming Efficiency and Retinogenesis in 3D Organoid Cultures

Diverse cell types can be reprogrammed into pluripotent stem cells by ectopic expression of Oct4 (Pou5f1), Klf4, Sox3, and Myc. Many of these induced pluripotent stem cells (iPSCs) retain memory, in terms of DNA methylation and histone modifications (epigenetic memory), of their cellular origins, and this may bias subsequent differentiation. Neurons are difficult to reprogram, and there has not been a systematic side-by-side characterization of reprogramming efficiency or epigenetic memory across different neuronal subtypes. Here, we compare reprogramming efficiency of five different retinal cell types at two different stages of development. Retinal differentiation from each iPSC line was measured using a quantitative standardized scoring system called STEM-RET and compared to the epigenetic memory. Neurons with the lowest reprogramming efficiency produced iPSC lines with the best retinal differentiation and were more likely to retain epigenetic memory of their cellular origins. In addition, we identified biomarkers of iPSCs that are predictive of retinal differentiation.

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.

Generation, transcriptome profiling, and functional validation of cone-rich human retinal organoids

Rod and cone photoreceptors are light-sensing cells in the human retina. Rods are dominant in the peripheral retina, whereas cones are enriched in the macula, which is responsible for central vision and visual acuity. Macular degenerations affect vision the most and are currently incurable. Here we report the generation, transcriptome profiling, and functional validation of cone-rich human retinal organoids differentiated from hESCs using an improved retinal differentiation system. Induced by extracellular matrix, aggregates of hESCs formed single-lumen cysts composed of epithelial cells with anterior neuroectodermal/ectodermal fates, including retinal cell fate. Then, the cysts were en bloc-passaged, attached to culture surface, and grew, forming colonies in which retinal progenitor cell patches were found. Following gentle cell detachment, retinal progenitor cells self-assembled into retinal epithelium-retinal organoid-that differentiated into stratified cone-rich retinal tissue in agitated cultures. Electron microscopy revealed differentiating outer segments of photoreceptor cells. Bulk RNA-sequencing profiling of time-course retinal organoids demonstrated that retinal differentiation in vitro recapitulated in vivo retinogenesis in temporal expression of cell differentiation markers and retinal disease genes, as well as in mRNA alternative splicing. Single-cell RNA-sequencing profiling of 8-mo retinal organoids identified cone and rod cell clusters and confirmed the cone enrichment initially revealed by quantitative microscopy. Notably, cones from retinal organoids and human macula had similar single-cell transcriptomes, and so did rods. Cones in retinal organoids exhibited electrophysiological functions. Collectively, we have established cone-rich retinal organoids and a reference of transcriptomes that are valuable resources for retinal studies.

Identification of factors associated with duplicate rate in ChIP-seq data

Chromatin immunoprecipitation and sequencing (ChIP-seq) has been widely used to map DNA-binding proteins, histone proteins and their modifications. ChIP-seq data contains redundant reads termed duplicates, referring to those mapping to the same genomic location and strand. There are two main sources of duplicates: polymerase chain reaction (PCR) duplicates and natural duplicates. Unlike natural duplicates that represent true signals from sequencing of independent DNA templates, PCR duplicates are artifacts originating from sequencing of identical copies amplified from the same DNA template. In analysis, duplicates are removed from peak calling and signal quantification. Nevertheless, a significant portion of the duplicates is believed to represent true signals. Obviously, removing all duplicates will underestimate the signal level in peaks and impact the identification of signal changes across samples. Therefore, an in-depth evaluation of the impact from duplicate removal is needed. Using eight public ChIP-seq datasets from three narrow-peak and two broad-peak marks, we tried to understand the distribution of duplicates in the genome, the extent by which duplicate removal impacts peak calling and signal estimation, and the factors associated with duplicate level in peaks. The three PCR-free histone H3 lysine 4 trimethylation (H3K4me3) ChIP-seq data had about 40% duplicates and 97% of them were within peaks. For the other datasets generated with PCR amplification of ChIP DNA, as expected, the narrow-peak marks have a much higher proportion of duplicates than the broad-peak marks. We found that duplicates are enriched in peaks and largely represent true signals, more conspicuous in those with high confidence. Furthermore, duplicate level in peaks is strongly correlated with the target enrichment level estimated using nonredundant reads, which provides the basis to properly allocate duplicates between noise and signal. Our analysis supports the feasibility of retaining the portion of signal duplicates into downstream analysis, thus alleviating the limitation of complete deduplication.