Two novel non-coding single nucleotide variants in the DNase1 hypersensitivity site of PRDM13 causing North Carolina macular dystrophy in Korea.

Purpose: Pathogenic variants in North Carolina macular dystrophy (NCMD) have rarely been reported in the East Asian population. Herein, we reported novel variants of NCMD in 2 Korean families. Methods: The regions associated with NCMD were analyzed with genome sequencing, and variants were filtered based on the minor allele frequency (0.5%) and heterozygosity. Non-coding variants were functionally annotated using multiple computational tools. Results: We identified two rare novel variants, chr6:g.99,598,914T>C (hg38; V17) and chr6:g.99,598,926G>A (hg38; V18) upstream of PRDM13 in families A and B, respectively. In Family 1, Grade 2 NCMD and a best-corrected visual acuity of 20/25 and 20/200 in the right and left eyes, respectively, were observed. In Family B, all affected individuals had Grade 1 NCMD with characteristic confluent drusen at the fovea and a best-corrected visual acuity of 20/20 in both eyes. These two variants are 10-22 bp downstream of the reported V10 variant within the DNase1 hypersensitivity site. This site is associated with progressive bifocal chorioretinal atrophy and congenital posterior polar chorioretinal hypertrophy and lies in the putative enhancer site of PRDM13. Conclusion: We identified two novel NCMD variants in the Korean population and further validated the regulatory role of the DNase1 hypersensitivity site upstream of PRDM13.

Lipid nanoparticle-mediated delivery of mRNA into the mouse and human retina and other ocular tissues.

Purpose: Lipid nanoparticles (LNPs) show promise in their ability to introduce mRNA to drive protein expression in specific cell types of the mammalian eye. Here, we examined the ability of mRNA encapsulated in lipid nanoparticles (LNPs) with two distinct formulations to drive gene expression in mouse and human retina and other ocular tissues. Methods: We introduced mRNA carrying LNPs into two biological systems. Intravitreal injections were tested to deliver LNPs into the mouse eye. Human retinal pigment epithelium (RPE) and retinal explants were used to assess mRNA expression in human tissue. We analyzed specificity of expression using histology, immunofluorescence, and imaging. Results: In mice, mRNAs encoding GFP and ciliary neurotrophic factor (CNTF) were specifically expressed by Müller glia and retinal pigment epithelium (RPE). Acute inflammatory changes measured by microglia distribution (Iba-1) or interleukin-6 (IL-6) expression were not observed 6 hours post-injection. Human RPE also expressed high levels of GFP. Human retinal explants expressed GFP in cells with apical and basal processes consistent with Müller glia and in perivascular cells consistent with macrophages. Conclusions: We demonstrated the ability to reliably transfect subpopulations of retinal cells in mice eye tissues in vivo and in human ocular tissues. Of significance, intravitreal injections were sufficient to transfect the RPE in mice. To our knowledge we demonstrate delivery of mRNA using LNPs in human ocular tissues for the first time.

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

Purpose: Prior studies have demonstrated the significance of specific cis-regulatory variants in retinal disease; however, determining the functional impact of regulatory variants remains a major challenge. In this study, we utilized a machine learning approach, trained on epigenomic data from the adult human retina, to systematically quantify the predicted impact of cis-regulatory variants. Methods: We used human retinal DNA accessibility data (ATAC-seq) to determine a set of 18.9k high-confidence, putative cis-regulatory elements. Eighty percent of these elements were used to train a machine learning model utilizing a gapped k-mer support vector machine-based approach. In silico saturation mutagenesis and variant scoring was applied to predict the functional impact of all potential single nucleotide variants within cis-regulatory elements. Impact scores were tested in a 20% hold-out dataset and compared to allele population frequency, phylogenetic conservation, transcription factor (TF) binding motifs, and existing massively parallel reporter assay data. Results: We generated a model that distinguishes between human retinal regulatory elements and negative test sequences with 95% accuracy. Among a hold-out test set of 3.7k human retinal CREs, all possible single nucleotide variants were scored. Variants with negative impact scores correlated with higher phylogenetic conservation of the reference allele, disruption of predicted TF binding motifs, and massively parallel reporter expression. Conclusions: We demonstrated the utility of human retinal epigenomic data to train a machine learning model for the purpose of predicting the impact of non-coding regulatory sequence variants. Our model accurately scored sequences and predicted putative transcription factor binding motifs. This approach has the potential to expedite the characterization of pathogenic non-coding sequence variants in the context of unexplained retinal disease. Translational Relevance: This workflow and resulting dataset serve as a promising genomic tool to facilitate the clinical prioritization of functionally disruptive non-coding mutations in the retina.

Monocarboxylate Transporter 1 (MCT1) Mediates Succinate Export in the Retina.

Purpose: Succinate is exported by the retina and imported by eyecup tissue. The transporters mediating this process have not yet been identified. Recent studies showed that monocarboxylate transporter 1 (MCT1) can transport succinate across plasma membranes in cardiac and skeletal muscle. Retina and retinal pigment epithelium (RPE) both express multiple MCT isoforms including MCT1. We tested the hypothesis that MCTs facilitate retinal succinate export and RPE succinate import. Methods: We assessed retinal succinate export and eyecup succinate import in short-term ex vivo culture using gas chromatography-mass spectrometry. We tested the dependence of succinate export and import on pH, proton ionophores, conventional MCT substrates, and the MCT inhibitors AZD3965, AR-C155858, and diclofenac. Results: Succinate exits retinal tissue through MCT1 but does not enter the RPE through MCT1 or any other MCT. Intracellular succinate levels are a contributing factor that determines if an MCT1-expressing tissue will export succinate. Conclusions: MCT1 facilitates export of succinate from retinas. An unidentified, non-MCT transporter facilitates import of succinate into RPE.

Cell-specific cis-regulatory elements and mechanisms of non-coding genetic disease in human retina and retinal organoids.

Cis-regulatory elements (CREs) play a critical role in the development and disease-states of all human cell types. In the retina, CREs have been implicated in several inherited disorders. To better characterize human retinal CREs, we performed single-nucleus assay for transposase-accessible chromatin sequencing (snATAC-seq) and single-nucleus RNA sequencing (snRNA-seq) on the developing and adult human retina and on induced pluripotent stem cell (iPSC)-derived retinal organoids. These analyses identified developmentally dynamic, cell-class-specific CREs, enriched transcription-factor-binding motifs, and putative target genes. CREs in the retina and organoids are highly correlated at the single-cell level, and this supports the use of organoids as a model for studying disease-associated CREs. As a proof of concept, we disrupted a disease-associated CRE at 5q14.3, confirming its principal target gene as the miR-9-2 primary transcript and demonstrating its role in neurogenesis and gene regulation in mature glia. This study provides a resource for characterizing human retinal CREs and showcases organoids as a model to study the function of CREs that influence development and disease.

Distinct features of brain perivascular fibroblasts and mural cells revealed by in vivo two-photon imaging.

Perivascular fibroblasts (PVFs) are recognized for their pro-fibrotic role in many central nervous system disorders. Like mural cells, PVFs surround blood vessels and express Pdgfrß. However, these shared attributes hinder the ability to distinguish PVFs from mural cells. We used in vivo two-photon imaging and transgenic mice with PVF-targeting promoters (Col1a1 or Col1a2) to compare the structure and distribution of PVFs and mural cells in cerebral cortex of healthy, adult mice. We show that PVFs localize to all cortical penetrating arterioles and their offshoots (arteriole-capillary transition zone), as well as the main trunk of only larger ascending venules. However, the capillary zone is devoid of PVF coverage. PVFs display short-range mobility along the vessel wall and exhibit distinct structural features (flattened somata and thin ruffled processes) not seen with smooth muscle cells or pericytes. These findings clarify that PVFs and mural cells are distinct cell types coexisting in a similar perivascular niche.

Gene regulatory networks controlling temporal patterning, neurogenesis, and cell-fate specification in mammalian retina.

Gene regulatory networks (GRNs), consisting of transcription factors and their target sites, control neurogenesis and cell-fate specification in the developing central nervous system. In this study, we use integrated single-cell RNA and single-cell ATAC sequencing (scATAC-seq) analysis in developing mouse and human retina to identify multiple interconnected, evolutionarily conserved GRNs composed of cell-type-specific transcription factors that both activate genes within their own network and inhibit genes in other networks. These GRNs control temporal patterning in primary progenitors, regulate transition from primary to neurogenic progenitors, and drive specification of each major retinal cell type. We confirm that NFI transcription factors selectively activate expression of genes promoting late-stage temporal identity in primary retinal progenitors and identify other transcription factors that regulate rod photoreceptor specification in postnatal retina. This study inventories cis- and trans-acting factors that control retinal development and can guide cell-based therapies aimed at replacing retinal neurons lost to disease.

Transient expression of a GABA receptor subunit during early development is critical for inhibitory synapse maturation and function.

Developing neural circuits, including GABAergic circuits, switch receptor types. But the role of early GABA receptor expression for establishment of functional inhibitory circuits remains unclear. Tracking the development of GABAergic synapses across axon terminals of retinal bipolar cells (BCs), we uncovered a crucial role of early GABAA receptor expression for the formation and function of presynaptic inhibitory synapses. Specifically, early α3-subunit-containing GABAA (GABAAα3) receptors are a key developmental organizer. Before eye opening, GABAAα3 gives way to GABAAα1 at individual BC presynaptic inhibitory synapses. The developmental downregulation of GABAAα3 is independent of GABAAα1 expression. Importantly, lack of early GABAAα3 impairs clustering of GABAAα1 and formation of functional GABAA synapses across mature BC terminals. This impacts the sensitivity of visual responses transmitted through the circuit. Lack of early GABAAα3 also perturbs aggregation of LRRTM4, the organizing protein at GABAergic synapses of rod BC terminals, and their arrangement of output ribbon synapses.

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.

Chromatin Environment and Cellular Context Specify Compensatory Activity of Paralogous MEF2 Transcription Factors.

Compensation among paralogous transcription factors (TFs) confers genetic robustness of cellular processes, but how TFs dynamically respond to paralog depletion on a genome-wide scale in vivo remains incompletely understood. Using single and double conditional knockout of myocyte enhancer factor 2 (MEF2) family TFs in granule neurons of the mouse cerebellum, we find that MEF2A and MEF2D play functionally redundant roles in cerebellar-dependent motor learning. Although both TFs are highly expressed in granule neurons, transcriptomic analyses show MEF2D is the predominant genomic regulator of gene expression in vivo. Strikingly, genome-wide occupancy analyses reveal upon depletion of MEF2D, MEF2A occupancy robustly increases at a subset of sites normally bound to MEF2D. Importantly, sites experiencing compensatory MEF2A occupancy are concentrated within open chromatin and undergo functional compensation for genomic activation and gene expression. Finally, motor activity induces a switch from non-compensatory to compensatory MEF2-dependent gene regulation. These studies uncover genome-wide functional interdependency between paralogous TFs in the brain.

Correction: Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

The original version of this Article contained an incorrect version of Fig. 3, which included two variants initially shown in black text in Fig. 3a that the authors removed from the final manuscript. The correct version of Fig. 3 without the two variants now appears in the PDF and HTML versions of the Article.

ABCA4-associated disease as a model for missing heritability in autosomal recessive disorders: novel noncoding splice, cis-regulatory, structural, and recurrent hypomorphic variants.

PURPOSE: ABCA4-associated disease, a recessive retinal dystrophy, is hallmarked by a large proportion of patients with only one pathogenic ABCA4 variant, suggestive for missing heritability. METHODS: By locus-specific analysis of ABCA4, combined with extensive functional studies, we aimed to unravel the missing alleles in a cohort of 67 patients (p), with one (p = 64) or no (p = 3) identified coding pathogenic variants of ABCA4. RESULTS: We identified eight pathogenic (deep-)intronic ABCA4 splice variants, of which five are novel and six structural variants, four of which are novel, including two duplications. Together, these variants account for the missing alleles in 40.3% of patients. Furthermore, two novel variants with a putative cis-regulatory effect were identified. The common hypomorphic variant c.5603A>T p.(Asn1868Ile) was found as a candidate second allele in 43.3% of patients. Overall, we have elucidated the missing heritability in 83.6% of our cohort. In addition, we successfully rescued three deep-intronic variants using antisense oligonucleotide (AON)-mediated treatment in HEK 293-T cells and in patient-derived fibroblast cells. CONCLUSION: Noncoding pathogenic variants, novel structural variants, and a common hypomorphic allele of the ABCA4 gene explain the majority of unsolved cases with ABCA4-associated disease, rendering this retinopathy a model for missing heritability in autosomal recessive disorders.

Contribution of noncoding pathogenic variants to RPGRIP1-mediated inherited retinal degeneration.

PURPOSE: With the advent of gene therapies for inherited retinal degenerations (IRDs), genetic diagnostics will have an increasing role in clinical decision-making. Yet the genetic cause of disease cannot be identified using exon-based sequencing for a significant portion of patients. We hypothesized that noncoding pathogenic variants contribute significantly to the genetic causality of IRDs and evaluated patients with single coding pathogenic variants in RPGRIP1 to test this hypothesis. METHODS: IRD families underwent targeted panel sequencing. Unsolved cases were explored by exome and genome sequencing looking for additional pathogenic variants. Candidate pathogenic variants were then validated by Sanger sequencing, quantitative polymerase chain reaction, and in vitro splicing assays in two cell lines analyzed through amplicon sequencing. RESULTS: Among 1722 families, 3 had biallelic loss-of-function pathogenic variants in RPGRIP1 while 7 had a single disruptive coding pathogenic variants. Exome and genome sequencing revealed potential noncoding pathogenic variants in these 7 families. In 6, the noncoding pathogenic variants were shown to lead to loss of function in vitro. CONCLUSION: Noncoding pathogenic variants were identified in 6 of 7 families with single coding pathogenic variants in RPGRIP1. The results suggest that noncoding pathogenic variants contribute significantly to the genetic causality of IRDs and RPGRIP1-mediated IRDs are more common than previously thought.

Correction: Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

The original version of this Article contained an error in the spelling of the author Anja K. Mayer, which was incorrectly given as Anja Kathrin Mayer. This has now been corrected in both the PDF and HTML versions of the Article.

Biallelic sequence and structural variants in RAX2 are a novel cause for autosomal recessive inherited retinal disease.

PURPOSE: RAX2 encodes a homeobox-containing transcription factor, in which four monoallelic pathogenic variants have been described in autosomal dominant cone-dominated retinal disease. METHODS: Exome sequencing in a European cohort with inherited retinal disease (IRD) (n = 2086) was combined with protein structure modeling of RAX2 missense variants, bioinformatics analysis of deletion breakpoints, haplotyping of RAX2 variant c.335dup, and clinical assessment of biallelic RAX2-positive cases and carrier family members. RESULTS: Biallelic RAX2 sequence and structural variants were found in five unrelated European index cases, displaying nonsyndromic autosomal recessive retinitis pigmentosa (ARRP) with an age of onset ranging from childhood to the mid-40s (average mid-30s). Protein structure modeling points to loss of function of the novel recessive missense variants and to a dominant-negative effect of the reported dominant RAX2 alleles. Structural variants were fine-mapped to disentangle their underlying mechanisms. Haplotyping of c.335dup in two cases suggests a common ancestry. CONCLUSION: This study supports a role for RAX2 as a novel disease gene for recessive IRD, broadening the mutation spectrum from sequence to structural variants and revealing a founder effect. The identification of biallelic RAX2 pathogenic variants in five unrelated families shows that RAX2 loss of function may be a nonnegligible cause of IRD in unsolved ARRP cases.

Mapping the genomic landscape of inherited retinal disease genes prioritizes genes prone to coding and noncoding copy-number variations.

PurposePart of the hidden genetic variation in heterogeneous genetic conditions such as inherited retinal diseases (IRDs) can be explained by copy-number variations (CNVs). Here, we explored the genomic landscape of IRD genes listed in RetNet to identify and prioritize those genes susceptible to CNV formation.MethodsRetNet genes underwent an assessment of genomic features and of CNV occurrence in the Database of Genomic Variants and literature. CNVs identified in an IRD cohort were characterized using targeted locus amplification (TLA) on extracted genomic DNA.ResultsExhaustive literature mining revealed 1,345 reported CNVs in 81 different IRD genes. Correlation analysis between rankings of genomic features and CNV occurrence demonstrated the strongest correlation between gene size and CNV occurrence of IRD genes. Moreover, we identified and delineated 30 new CNVs in IRD cases, 13 of which are novel and three of which affect noncoding, putative cis-regulatory regions. Finally, the breakpoints of six complex CNVs were determined using TLA in a hypothesis-neutral manner.ConclusionWe propose a ranking of CNV-prone IRD genes and demonstrate the efficacy of TLA for the characterization of CNVs on extracted DNA. Finally, this IRD-oriented CNV study can serve as a paradigm for other genetically heterogeneous Mendelian diseases with hidden genetic variation.

MEF2D drives photoreceptor development through a genome-wide competition for tissue-specific enhancers.

Organismal development requires the precise coordination of genetic programs to regulate cell fate and function. MEF2 transcription factors (TFs) play essential roles in this process but how these broadly expressed factors contribute to the generation of specific cell types during development is poorly understood. Here we show that despite being expressed in virtually all mammalian tissues, in the retina MEF2D binds to retina-specific enhancers and controls photoreceptor cell development. MEF2D achieves specificity by cooperating with a retina-specific factor CRX, which recruits MEF2D away from canonical MEF2 binding sites and redirects it to retina-specific enhancers that lack the consensus MEF2-binding sequence. Once bound to retina-specific enhancers, MEF2D and CRX co-activate the expression of photoreceptor-specific genes that are critical for retinal function. These findings demonstrate that broadly expressed TFs acquire specific functions through competitive recruitment to enhancers by tissue-specific TFs and through selective activation of these enhancers to regulate tissue-specific genes.

NeuroD factors regulate cell fate and neurite stratification in the developing retina.

Members of the basic helix-loop-helix (bHLH) family of transcription factors have been shown to control critical aspects of development in many tissues. To identify bHLH genes that might regulate specific aspects of retinal cell development, we investigated the expression of bHLH genes in single, developing mouse retinal cells, with particular emphasis on the NeuroD family. Two of these factors, NeuroD2 and NeuroD6/NEX, had not been previously reported as expressed in the retina. A series of loss- and gain-of-function experiments was performed, which suggested that NeuroD genes have both similarities and differences in their activities. Notably, misexpression of NeuroD genes can direct amacrine cell processes to two to three specific sublaminae in the inner plexiform layer. This effect is specific to cell type and NeuroD gene, as the AII amacrine cell type is refractory to the effects of NeuroD1 and NeuroD6, but uniquely sensitive to the effect of NeuroD2 on neurite targeting. Additionally, NeuroD2 is endogenously expressed in AII amacrine cells, among others, and loss of NeuroD2 function results in a partial loss of AII amacrine cells. The effects of misexpressing NeuroD genes on retinal cell fate determination also suggested shared and divergent functions. Remarkably, NeuroD2 misexpression induced ganglion cell production even after the normal developmental window of ganglion cell genesis. Together, these data suggest that members of the NeuroD family are important for neuronal cell type identity and may be involved in several cell type-specific aspects of retinal development, including fate determination, differentiation, morphological development, and circuit formation.

Development and diversification of retinal amacrine interneurons at single cell resolution.

The vertebrate retina uses diverse neuronal cell types arrayed into complex neural circuits to extract, process, and relay information from the visual scene to the higher order processing centers of the brain. Amacrine cells, a class of interneurons, are thought to mediate much of the processing of the visual signal that occurs within the retina. Although amacrine cells display extensive morphological diversity, the molecular nature of this diversity is largely unknown. Furthermore, it is not known how this diversity arises during development. Here, we have combined in vivo genetic labeling, single cell genome-wide expression profiling, and classical birthdating to (i) identify specific molecular types of amacrine cells, (ii) demonstrate the molecular diversity of the amacrine cell class, and (iii) show that amacrine cell diversity arises at least in part through temporal patterning.

Early asymmetry of gene transcription in embryonic human left and right cerebral cortex.

The human left and right cerebral hemispheres are anatomically and functionally asymmetric. To test whether human cortical asymmetry has a molecular basis, we studied gene expression levels between the left and right embryonic hemispheres using serial analysis of gene expression (SAGE). We identified and verified 27 differentially expressed genes, which suggests that human cortical asymmetry is accompanied by early, marked transcriptional asymmetries. LMO4 is consistently more highly expressed in the right perisylvian human cerebral cortex than in the left and is essential for cortical development in mice, suggesting that human left-right specialization reflects asymmetric cortical development at early stages.