Single-cell RNA-seq is a powerful tool in decoding the heterogeneity in complex tissues by generating transcriptomic profiles of the individual cell. Here, we report a single-nuclei RNA-seq (snRNA-seq) transcriptomic study on human retinal tissue, which is composed of multiple cell types with distinct functions. Six samples from three healthy donors are profiled and high-quality RNA-seq data is obtained for 5873 single nuclei. All major retinal cell types are observed and marker genes for each cell type are identified. The gene expression of the macular and peripheral retina is compared to each other at cell-type level. Furthermore, our dataset shows an improved power for prioritizing genes associated with human retinal diseases compared to both mouse single-cell RNA-seq and human bulk RNA-seq results. In conclusion, we demonstrate that obtaining single cell transcriptomes from human frozen tissues can provide insight missed by either human bulk RNA-seq or animal models.

Inherited retinal degenerations are a common cause of untreatable blindness worldwide, with retinitis pigmentosa and cone dystrophy affecting approximately 1 in 3500 and 1 in 10,000 individuals, respectively. A major limitation to the development of effective therapies is the lack of availability of animal models that fully replicate the human condition. Particularly for cone disorders, rodent, canine, and feline models with no true macula have substantive limitations. By contrast, the cone-rich macula of a nonhuman primate (NHP) closely mirrors that of the human retina. Consequently, well-defined NHP models of heritable retinal diseases, particularly cone disorders that are predictive of human conditions, are necessary to more efficiently advance new therapies for patients. We have identified 4 related NHPs at the California National Primate Research Center with visual impairment and findings from clinical ophthalmic examination, advanced retinal imaging, and electrophysiology consistent with achromatopsia. Genetic sequencing confirmed a homozygous R565Q missense mutation in the catalytic domain of PDE6C, a cone-specific phototransduction enzyme associated with achromatopsia in humans. Biochemical studies demonstrate that the mutant mRNA is translated into a stable protein that displays normal cellular localization but is unable to hydrolyze cyclic GMP (cGMP). This NHP model of a cone disorder will not only serve as a therapeutic testing ground for achromatopsia gene replacement, but also for optimization of gene editing in the macula and of cone cell replacement in general.

Cell classes in the human retina are highly heterogeneous with their abundance varying by several orders of magnitude. Here, we generated and integrated a multi-omics single-cell atlas of the adult human retina, including more than 250,000 nuclei for single-nuclei RNA-seq and 137,000 nuclei for single-nuclei ATAC-seq. Cross-species comparison of the retina atlas among human, monkey, mice, and chicken revealed relatively conserved and non-conserved types. Interestingly, the overall cell heterogeneity in primate retina decreases compared with that of rodent and chicken retina. Through integrative analysis, we identified 35,000 distal cis-element-gene pairs, constructed transcription factor (TF)-target regulons for more than 200 TFs, and partitioned the TFs into distinct co-active modules. We also revealed the heterogeneity of the cis-element-gene relationships in different cell types, even from the same class. Taken together, we present a comprehensive single-cell multi-omics atlas of the human retina as a resource that enables systematic molecular characterization at individual cell-type resolution.

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide. However, its molecular and cellular causes are not well understood. Elevated intraocular pressure (IOP) is a major risk factor, but many patients have normal IOP. Colocalization and Mendelian randomization analysis of >240 POAG and IOP genome-wide association study (GWAS) loci and overlapping expression and splicing quantitative trait loci (e/sQTLs) in 49 GTEx tissues and retina prioritizes causal genes for 60% of loci. These genes are enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues reveals that the POAG and IOP colocalizing genes and genome-wide associations are enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominates IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

Integration of single-cell multiomics profiles generated by different single-cell technologies from the same biological sample is still challenging. Previous approaches based on shared features have only provided approximate solutions. Here, we present a novel mathematical solution named bi-order canonical correlation analysis (bi-CCA), which extends the widely used CCA approach to iteratively align the rows and the columns between data matrices. Bi-CCA is generally applicable to combinations of any two single-cell modalities. Validations using co-assayed ground truth data and application to a CAR-NK study and a fetal muscle atlas demonstrate its capability in generating accurate multimodal co-embeddings and discovering cellular identity.

EGCG protects transduced PC12 cells against α-Syn-induced cytotoxicity by inhibiting the overexpression and fibrillation of α-Syn in the cells.