To identify specific mutations causing North Carolina macular dystrophy (NCMD).
Whole-genome sequencing coupled with reverse transcription polymerase chain reaction (RT-PCR) analysis of gene ...expression in human retinal cells.
A total of 141 members of 12 families with NCMD and 261 unrelated control individuals.
Genome sequencing was performed on 8 affected individuals from 3 families affected with chromosome 6–linked NCMD (MCDR1) and 2 individuals affected with chromosome 5–linked NCMD (MCDR3). Variants observed in the MCDR1 locus with frequencies <1% in published databases were confirmed using Sanger sequencing. Confirmed variants absent from all published databases were sought in 8 additional MCDR1 families and 261 controls. The RT-PCR analysis of selected genes was performed in stem cell–derived human retinal cells.
Co-segregation of rare genetic variants with disease phenotype.
Five sequenced individuals with MCDR1-linked NCMD shared a haplotype of 14 rare variants spanning 1 Mb of the disease-causing allele. One of these variants (V1) was absent from all published databases and all 261 controls, but was found in 5 additional NCMD kindreds. This variant lies in a DNase 1 hypersensitivity site (DHS) upstream of both the PRDM13 and CCNC genes. Sanger sequencing of 1 kb centered on V1 was performed in the remaining 4 NCMD probands, and 2 additional novel single nucleotide variants (V2 in 3 families and V3 in 1 family) were identified in the DHS within 134 bp of the location of V1. A complete duplication of the PRDM13 gene was also discovered in a single family (V4). The RT-PCR analysis of PRDM13 expression in developing retinal cells revealed marked developmental regulation. Next-generation sequencing of 2 individuals with MCDR3-linked NCMD revealed a 900-kb duplication that included the entire IRX1 gene (V5). The 5 mutations V1 to V5 segregated perfectly in the 102 affected and 39 unaffected members of the 12 NCMD families.
We identified 5 rare mutations, each capable of arresting human macular development. Four of these strongly implicate the involvement of PRDM13 in macular development, whereas the pathophysiologic mechanism of the fifth remains unknown but may involve the developmental dysregulation of IRX1.
Gene expression provides valuable insight into cell function. As such, vision researchers have frequently employed gene expression studies to better understand retinal physiology and disease. With ...the advent of single-cell RNA sequencing, expression experiments provide an unparalleled resolution of information. Instead of studying aggregated gene expression across all cells in a heterogenous tissue, single-cell technology maps RNA to an individual cell, which facilitates grouping of retinal and choroidal cell types for further study. Single-cell RNA sequencing has been quickly adopted by both basic and translational vision researchers, and single-cell level gene expression has been studied in the visual systems of animal models, retinal organoids, and primary human retina, RPE, and choroid. These experiments have generated detailed atlases of gene expression and identified new retinal cell types. Likewise, single-cell RNA sequencing investigations have characterized how gene expression changes in the setting of many retinal diseases, including how choroidal endothelial cells are altered in age-related macular degeneration. In addition, this technology has allowed vision researchers to discover drivers of retinal development and model rare retinal diseases with induced pluripotent stem cells. In this review, we will overview the growing number of single-cell RNA sequencing studies in the field of vision research. We will summarize experimental considerations for designing single-cell RNA sequencing experiments and highlight important advancements in retinal, RPE, choroidal, and retinal organoid biology driven by this technology. Finally, we generalize these findings to genes involved in retinal degeneration and outline the future of single-cell expression experiments in studying retinal disease.
•Single-cell RNA sequencing allows for high resolution transcriptome profiling of heterogenous ocular tissues.•Human cone photoreceptor cells and Muller cells have different expression profiles in the fovea versus the periphery.•Distinct populations of choriocapillaris, arterial, and venous endothelial cells can be identified with scRNA-seq.•scRNA-seq of patient-derived induced pluripotent stem cell models promotes discovery of new disease pathogenesis models.
The identification of mutations in genes that cause human diseases has largely been accomplished through the use of positional cloning, which relies on linkage mapping. In studies of rare diseases, ...the resolution of linkage mapping is limited by the number of available meioses and informative marker density. One recent advance is the development of high-density SNP microarrays for genotyping. The SNP arrays overcome low marker informativity by using a large number of markers to achieve greater coverage at finer resolution. We used SNP microarray genotyping for homozygosity mapping in a small consanguineous Israeli Bedouin family with autosomal recessive Bardet-Biedl syndrome (BBS; obesity, pigmentary retinopathy, polydactyly, hypogonadism, renal and cardiac abnormalities, and cognitive impairment) in which previous linkage studies using short tandem repeat polymorphisms failed to identify a disease locus. SNP genotyping revealed a homozygous candidate region. Mutation analysis in the region of homozygosity identified a conserved homozygous missense mutation in the TRIM32 gene, a gene coding for an E3 ubiquitin ligase. Functional analysis of this gene in zebrafish and expression correlation analyses among other BBS genes in an expression quantitative trait loci data set demonstrate that TRIM32 is a BBS gene. This study shows the value of high-density SNP genotyping for homozygosity mapping and the use of expression correlation data for evaluation of candidate genes and identifies the proteasome degradation pathway as a pathway involved in BBS.
Mutations in myocilin (MYOC) are the leading known genetic cause of primary open-angle glaucoma, responsible for about 4% of all cases. Mutations in MYOC cause a gain-of-function phenotype in which ...mutant myocilin accumulates in the endoplasmic reticulum (ER) leading to ER stress and trabecular meshwork (TM) cell death. Therefore, knocking out myocilin at the genome level is an ideal strategy to permanently cure the disease. We have previously utilized CRISPR/Cas9 genome editing successfully to target MYOC using adenovirus 5 (Ad5). However, Ad5 is not a suitable vector for clinical use. Here, we sought to determine the efficacy of adeno-associated viruses (AAVs) and lentiviruses (LVs) to target the TM. First, we examined the TM tropism of single-stranded (ss) and self-complimentary (sc) AAV serotypes as well as LV expressing GFP via intravitreal (IVT) and intracameral (IC) injections. We observed that LV_GFP expression was more specific to the TM injected via the IVT route. IC injections of Trp-mutant scAAV2 showed a prominent expression of GFP in the TM. However, robust GFP expression was also observed in the ciliary body and retina. We next constructed lentiviral particles expressing Cas9 and guide RNA (gRNA) targeting MYOC (crMYOC) and transduction of TM cells stably expressing mutant myocilin with LV_crMYOC significantly reduced myocilin accumulation and its associated chronic ER stress. A single IVT injection of LV_crMYOC in Tg-MYOC
mice decreased myocilin accumulation in TM and reduced elevated IOP significantly. Together, our data indicates, LV_crMYOC targets MYOC gene editing in TM and rescues a mouse model of myocilin-associated glaucoma.
Proper spatial differentiation of retinal cell types is necessary for normal human vision. Many retinal diseases, such as Best disease and male germ cell associated kinase (MAK)-associated retinitis ...pigmentosa, preferentially affect distinct topographic regions of the retina. While much is known about the distribution of cell types in the retina, the distribution of molecular components across the posterior pole of the eye has not been well-studied. To investigate regional difference in molecular composition of ocular tissues, we assessed differential gene expression across the temporal, macular, and nasal retina and retinal pigment epithelium (RPE)/choroid of human eyes using RNA-Seq. RNA from temporal, macular, and nasal retina and RPE/choroid from four human donor eyes was extracted, poly-A selected, fragmented, and sequenced as 100 bp read pairs. Digital read files were mapped to the human genome and analyzed for differential expression using the Tuxedo software suite. Retina and RPE/choroid samples were clearly distinguishable at the transcriptome level. Numerous transcription factors were differentially expressed between regions of the retina and RPE/choroid. Photoreceptor-specific genes were enriched in the peripheral samples, while ganglion cell and amacrine cell genes were enriched in the macula. Within the RPE/choroid, RPE-specific genes were upregulated at the periphery while endothelium associated genes were upregulated in the macula. Consistent with previous studies, BEST1 expression was lower in macular than extramacular regions. The MAK gene was expressed at lower levels in macula than in extramacular regions, but did not exhibit a significant difference between nasal and temporal retina. The regional molecular distinction is greatest between macula and periphery and decreases between different peripheral regions within a tissue. Datasets such as these can be used to prioritize candidate genes for possible involvement in retinal diseases with regional phenotypes.
•We examined gene expression in temporal, macular, and nasal regions of human retina and RPE/choroid using RNA-Seq.•Expression differences between macula and periphery in both tissues reflect the distribution of cell types.•Nasal and temporal regions of neural retina are indistinguishable in our analysis.
Single-cell RNA sequencing has revolutionized ocular gene expression studies. This technology has enabled researchers to identify expression signatures for rare cell types and characterize how gene ...expression changes across biological conditions, such as topographic region or disease status. However, sharing single-cell RNA sequencing results remains a major obstacle, particular for individuals without a computational background. To address these limitations, we developed Spectacle, an interactive web-based resource for exploring previously published single-cell RNA sequencing data from ocular studies. Spectacle is powered by a locally developed R package, cellcuratoR, which utilizes the Shiny framework in R to generate interactive visualizations for single-cell expression data. Spectacle contains five pre-processed ocular single-cell RNA sequencing data sets and is accessible via the web at OcularGeneExpression.org/singlecell. With Spectacle, users can interactively identify which cell types express a gene of interest, detect transcriptomic subpopulations within a cell type, and perform highly flexible differential expression analyses. The freely-available Spectacle system reduces the bioinformatic barrier for interacting with rich single-cell RNA sequencing studies from ocular tissues, making it easy to quickly identify cell types that express a gene of interest.
•Spectacle is a resource for exploring ocular single-cell RNA sequencing datasets.•Spectacle includes five datasets from the retina, RPE, and choroid.•Spectacle facilitates interactive expression analyses on pre-classified cell types.•Spectacle is freely available at OcularGeneExpression.org/singlecell.
The human choroid is a heterogeneous, highly vascular connective tissue that dysfunctions in age-related macular degeneration (AMD). In this study, we performed single-cell RNA sequencing on 21 human ...choroids, 11 of which were derived from donors with early atrophic or neovascular AMD. Using this large donor cohort, we identified new gene expression signatures and immunohistochemically characterized discrete populations of resident macrophages, monocytes/inflammatory macrophages and dendritic cells. These three immune populations demonstrated unique expression patterns for AMD genetic risk factors, with dendritic cells possessing the highest expression of the neovascular AMD-associated MMP9 gene. Additionally, we performed trajectory analysis to model transcriptomic changes across the choroidal vasculature, and we identified expression signatures for endothelial cells from choroidal arterioles and venules. Finally, we performed differential expression analysis between control, early atrophic AMD, and neovascular AMD samples, and we observed that early atrophic AMD samples had high expression of SPARCL1, a gene that has been shown to increase in response to endothelial damage. Likewise, neovascular endothelial cells harbored gene expression changes consistent with endothelial cell damage and demonstrated increased expression of the sialomucins CD34 and ENCM, which were also observed at the protein level within neovascular membranes. Overall, this study characterizes the molecular features of new populations of choroidal endothelial cells and mononuclear phagocytes in a large cohort of AMD and control human donors.
We used expression quantitative trait locus mapping in the laboratory rat (Rattus norvegicus) to gain a broad perspective of gene regulation in the mammalian eye and to identify genetic variation ...relevant to human eye disease. Of >31,000 gene probes represented on an Affymetrix expression microarray, 18,976 exhibited sufficient signal for reliable analysis and at least 2-fold variation in expression among 120 F₂ rats generated from an SR/JrHsd x SHRSP intercross. Genome-wide linkage analysis with 399 genetic markers revealed significant linkage with at least one marker for 1,300 probes (α = 0.001; estimated empirical false discovery rate = 2%). Both contiguous and noncontiguous loci were found to be important in regulating mammalian eye gene expression. We investigated one locus of each type in greater detail and identified putative transcription-altering variations in both cases. We found an inserted cREL binding sequence in the 5' flanking sequence of the Abca4 gene associated with an increased expression level of that gene, and we found a mutation of the gene encoding thyroid hormone receptor β2 associated with a decreased expression level of the gene encoding short-wavelength sensitive opsin (Opn1sw). In addition to these positional studies, we performed a pairwise analysis of gene expression to identify genes that are regulated in a coordinated manner and used this approach to validate two previously undescribed genes involved in the human disease Bardet-Biedl syndrome. These data and analytical approaches can be used to facilitate the discovery of additional genes and regulatory elements involved in human eye disease.
The human neural retina is a light sensitive tissue with remarkable spatial and cellular organization. Compared with the periphery, the central retina contains more densely packed cone photoreceptor ...cells with unique morphologies and synaptic wiring. Some regions of the central retina exhibit selective degeneration or preservation in response to retinal disease and the basis for this variation is unknown. In this study, we used both bulk and single-cell RNA sequencing to compare gene expression within concentric regions of the central retina. We identified unique gene expression patterns of foveal cone photoreceptor cells, including many foveal-enriched transcription factors. In addition, we found that the genes RORB1, PPFIA1 and KCNAB2 are differentially spliced in the foveal, parafoveal and macular regions. These results provide a highly detailed spatial characterization of the retinal transcriptome and highlight unique molecular features of different retinal regions.