Wenjun Yan

Single-cell sequencing has revolutionized the scale and resolution of molecular profiling of tissues and organs. Here we present an integrated dual-modal reference atlas of the most accessible portion of the mammalian central nervous system, the retina. We compiled around 3.9 million cells from 125 donors of diverse ancestral backgrounds, including 8 published studies and 2.7 million unpublished data points, to create a comprehensive human retina cell atlas (HRCA) with more than 130 cell types identified. We annotated each cluster, identified marker genes and characterized cis-regulatory elements and gene regulatory networks. Our analysis uncovered differences in transcriptome, chromatin and gene regulatory networks across cell types. We modeled changes in gene expression and chromatin accessibility across age, ancestry and tissue region. This integrated atlas enhanced the fine-mapping of genome-wide association study and expression quantitative trait loci variants. Accessible through interactive browsers, this multimodal multidonor and multilab HRCA can facilitate a better understanding of retinal function and pathology.

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.

The basic plan of the retina is conserved across vertebrates, yet species differ profoundly in their visual needs¹. Retinal cell types may have evolved to accommodate these varied needs, but this has not been systematically studied. Here we generated and integrated single-cell transcriptomic atlases of the retina from 17 species: humans, two non-human primates, four rodents, three ungulates, opossum, ferret, tree shrew, a bird, a reptile, a teleost fish and a lamprey. We found high molecular conservation of the six retinal cell classes (photoreceptors, horizontal cells, bipolar cells, amacrine cells, retinal ganglion cells (RGCs) and Müller glia), with transcriptomic variation across species related to evolutionary distance. Major subclasses were also conserved, whereas variation among cell types within classes or subclasses was more pronounced. However, an integrative analysis revealed that numerous cell types are shared across species, based on conserved gene expression programmes that are likely to trace back to an early ancestral vertebrate. The degree of variation among cell types increased from the outer retina (photoreceptors) to the inner retina (RGCs), suggesting that evolution acts preferentially to shape the retinal output. Finally, we identified rodent orthologues of midget RGCs, which comprise more than 80% of RGCs in the human retina, subserve high-acuity vision, and were previously believed to be restricted to primates². By contrast, the mouse orthologues have large receptive fields and comprise around 2% of mouse RGCs. Projections of both primate and mouse orthologous types are overrepresented in the thalamus, which supplies the primary visual cortex. We suggest that midget RGCs are not primate innovations, but are descendants of evolutionarily ancient types that decreased in size and increased in number as primates evolved, thereby facilitating high visual acuity and increased cortical processing of visual information.

The intrinsic mechanisms that regulate neurotoxic versus neuroprotective astrocyte phenotypes and their effects on central nervous system (CNS) degeneration and repair remain poorly understood. Here, we show injured white matter astrocytes differentiate into two distinct C3-positive and C3-negative reactive populations, previously simplified as neurotoxic (A1) and neuroprotective (A2)^1,2, which can be further subdivided into unique subpopulations defined by proliferation and differential gene expression signatures. We find the balance of neurotoxic versus neuroprotective astrocytes is regulated by discrete pools of compartmented cAMP derived from soluble adenylyl cyclase (sAC) and show proliferating neuroprotective astrocytes inhibit microglial activation and downstream neurotoxic astrocyte differentiation to promote retinal ganglion cell (RGC) survival. Finally, we report a new, therapeutically tractable viral vector to specifically target optic nerve head astrocytes and show elevating nuclear or depleting cytoplasmic cAMP in reactive astrocytes inhibits deleterious microglial/macrophage cell activation and promotes RGC survival after optic nerve injury. Thus, soluble adenylyl cyclase and compartmented, nuclear- and cytoplasmic-localized cAMP in reactive astrocytes act as a molecular switch for neuroprotective astrocyte reactivity that can be targeted to inhibit microglial activation and neurotoxic astrocyte differentiation to therapeutic effect. These data expand upon and define new reactive astrocyte subtypes and represent a novel step toward the development of gliotherapeutics for the treatment of glaucoma and other optic neuropathies.

Single-cell sequencing has revolutionized the scale and resolution of molecular profiling of tissues and organs. Here, we present an integrated multimodal reference atlas of the most accessible portion of the mammalian central nervous system, the retina. We compiled around 2.4 million cells from 55 donors, including 1.4 million unpublished data points, to create a comprehensive human retina cell atlas (HRCA) of transcriptome and chromatin accessibility, unveiling over 110 types. Engaging the retina community, we annotated each cluster, refined the Cell Ontology for the retina, identified distinct marker genes, and characterized cis-regulatory elements and gene regulatory networks (GRNs) for these cell types. Our analysis uncovered intriguing differences in transcriptome, chromatin, and GRNs across cell types. In addition, we modeled changes in gene expression and chromatin openness across gender and age. This integrated atlas also enabled the fine-mapping of GWAS and eQTL variants. Accessible through interactive browsers, this multimodal cross-donor and cross-lab HRCA, can facilitate a better understanding of retinal function and pathology.

Non-neuronal cells are key to the complex cellular interplay that follows central nervous system insult. To understand this interplay, we generated a single-cell atlas of immune, glial and retinal pigment epithelial cells from adult mouse retina before and at multiple time points after axonal transection. We identified rare subsets in naive retina, including interferon (IFN)-response glia and border-associated macrophages, and delineated injury-induced changes in cell composition, expression programs and interactions. Computational analysis charted a three-phase multicellular inflammatory cascade after injury. In the early phase, retinal macroglia and microglia were reactivated, providing chemotactic signals concurrent with infiltration of CCR2⁺ monocytes from the circulation. These cells differentiated into macrophages in the intermediate phase, while an IFN-response program, likely driven by microglia-derived type I IFN, was activated across resident glia. The late phase indicated inflammatory resolution. Our findings provide a framework to decipher cellular circuitry, spatial relationships and molecular interactions following tissue injury.

Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide; however, the 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 GWAS loci and of overlapping eQTLs and sQTLs in 49 GTEx tissues and retina prioritized causal genes for 60% of loci. These genes were enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues revealed that the colocalizing genes and genome-wide POAG and IOP associations were enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominated IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.

There is a gap in understanding the effect of the essential ω-3 and ω-6 long-chain polyunsaturated fatty acids (LCPUFA) on Phase I retinopathy of prematurity (ROP), which precipitates proliferative ROP. Postnatal hyperglycemia contributes to Phase I ROP by delaying retinal vascularization. In mouse neonates with hyperglycemia-associated Phase I retinopathy, dietary ω-3 (vs. ω-6 LCPUFA) supplementation promoted retinal vessel development. However, ω-6 (vs. ω-3 LCPUFA) was also developmentally essential, promoting neuronal growth and metabolism as suggested by a strong metabolic shift in almost all types of retinal neuronal and glial cells identified with single-cell transcriptomics. Loss of adiponectin (APN) in mice (mimicking the low APN levels in Phase I ROP) decreased LCPUFA levels (including ω-3 and ω-6) in retinas under normoglycemic and hyperglycemic conditions. ω-3 (vs. ω-6) LCPUFA activated the APN pathway by increasing the circulating APN levels and inducing expression of the retinal APN receptor. Our findings suggested that both ω-3 and ω-6 LCPUFA are crucial in protecting against retinal neurovascular dysfunction in a Phase I ROP model; adequate ω-6 LCPUFA levels must be maintained in addition to ω-3 supplementation to prevent retinopathy. Activation of the APN pathway may further enhance the ω-3 and ω-6 LCPUFA’s protection against ROP.