A single-cell transcriptome atlas of the aging human and macaque retina

Cross-species single-cell landscapes identify the pathogenic gene characteristics of inherited retinal diseases

Introduction

Inherited retinal diseases (IRDs) affect ∼4.5 million people worldwide. Elusive pathogenic variants in over 280 genes are associated with one or more clinical forms of IRDs. It is necessary to understand the complex interaction among retinal cell types and pathogenic genes by constructing a regulatory network. In this study, we attempt to establish a panoramic expression view of the cooperative work in retinal cells to understand the clinical manifestations and pathogenic bases underlying IRDs.

Methods

Single-cell RNA sequencing (scRNA-seq) data on the retinas from 35 retina samples of 3 species (human, mouse, and zebrafish) including 259,087 cells were adopted to perform a comparative analysis across species. Bioinformatic tools were used to conduct weighted gene co-expression network analysis (WGCNA), single-cell regulatory network analysis, cell-cell communication analysis, and trajectory inference analysis.

Results

The cross-species comparison revealed shared or species-specific gene expression patterns at single-cell resolution, such as the stathmin family genes, which were highly expressed specifically in zebrafish Müller glias (MGs). Thirteen gene modules were identified, of which nine were associated with retinal cell types, and Gene Ontology (GO) enrichment of module genes was consistent with cell-specific highly expressed genes. Many IRD genes were identified as hub genes and cell-specific regulons. Most IRDs, especially the retinitis pigmentosa (RP) genes, were enriched in rod-specific regulons. Integrated expression and transcription regulatory network genes, such as congenital stationary night blindness (CSNB) genes GRK1, PDE6B, and TRPM1, showed cell-specific expression and transcription characteristics in either rods or bipolar cells (BCs). IRD genes showed evolutionary conservation (GNAT2, PDE6G, and SAG) and divergence (GNAT2, MT-ND4, and PDE6A) along the trajectory of photoreceptors (PRs) among species. In particular, the Leber congenital amaurosis (LCA) gene OTX2 showed high expression at the beginning of the trajectory of both PRs and BCs.

Conclusion

We identified molecular pathways and cell types closely connected with IRDs, bridging the gap between gene expression, genetics, and pathogenesis. The IRD genes enriched in cell-specific modules and regulons suggest that these diseases share common etiological bases. Overall, mining of interspecies transcriptome data reveals conserved transcriptomic features of retinas across species and promising applications in both normal retina anatomy and retina pathology.

SPI1-mediated macrophage polarization aggravates age-related macular degeneration

Objective

This study revealed a core regulator and common upstream mechanisms for the multifaceted pathological processes of age-related macular degeneration (AMD) and provided proof-of-concept for this new therapeutic target.

Methods

Comprehensive gene expression analysis was performed using RNA sequencing of eye cup from old mice as well as laser-induced choroidal neovascularization (CNV) mouse model. Through integrative analysis and protein-protein interaction (PPI) analysis, common pathways and key transcription factor was identified simultaneously engaged in age-related retinal degeneration and CNV, the two typical pathological process of AMD. Subsequently, the expression changes of Spi1, the key regulator, as well as the alternation of the downstream mechanisms were validated in both models through qRT-PCR, Elisa, flow cytometry and immunofluorescence. Further, we assessed the impact of Spi1 knockdown in vitro and in vivo using gene intervention vectors carried by adeno-associated virus or lentivirus to test its potential as a therapeutic target.

Results

Compared to corresponding controls, we found 1,939 and 1,319 genes differentially expressed in eye cups of old and CNV mice respectively. The integrative analysis identified a total of 275 overlapping DEGs, of which 150 genes were co-upregulated. PPI analysis verified a central transcription factor, SPI1. The significant upregulation of Spi1 expression was then validated in both models, accompanied by macrophage polarization towards the M1 phenotype. Finally, SPI1 suppression significantly inhibited M1 polarization of BMDMs and attenuated neovascularization in CNV mice.

Conclusion

This study demonstrates that SPI1 exerts a pivotal role in AMD by regulation of macrophage polarization and innate immune response, offering promise as an innovative target for treating AMD.

Age-related dysregulation of the retinal transcriptome in African turquoise killifish

Age-related vision loss caused by retinal neurodegenerative pathologies is becoming more prevalent in our ageing society. To understand the physiological and molecular impact of ageing on retinal homeostasis, we used the short-lived African turquoise killifish, a model known to naturally develop central nervous system (CNS) ageing hallmarks and vision loss. Bulk and single-cell RNA-sequencing (scRNAseq) of three age groups (6-, 12-, and 18-week-old) identified transcriptional ageing fingerprints in the killifish retina, unveiling pathways also identified in the aged brain, including oxidative stress, gliosis, and inflammageing. These findings were comparable to observations in the ageing mouse retina. Additionally, transcriptional changes in genes related to retinal diseases, such as glaucoma and age-related macular degeneration, were observed. The cellular heterogeneity in the killifish retina was characterized, confirming the presence of all typical vertebrate retinal cell types. Data integration from age-matched samples between the bulk and scRNAseq experiments revealed a loss of cellular specificity in gene expression upon ageing, suggesting potential disruption in transcriptional homeostasis. Differential expression analysis within the identified cell types highlighted the role of glial/immune cells as important stress regulators during ageing. Our work emphasizes the value of the fast-ageing killifish in elucidating molecular signatures in age-associated retinal disease and vision decline. This study contributes to the understanding of how age-related changes in molecular pathways may impact CNS health, providing insights that may inform future therapeutic strategies for age-related pathologies.

Human skeletal muscle aging atlas

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

Age-related dysregulation of the retinal transcriptome in African turquoise killifish

Age-related vision loss caused by retinal neurodegenerative pathologies is becoming more prevalent in our ageing society. To understand the physiological and molecular impact of ageing on retinal homeostasis, we used the short-lived African turquoise killifish, a model known to naturally develop central nervous system (CNS) ageing hallmarks and vision loss. Bulk and single-cell RNA-sequencing (scRNA-seq) of three age groups (6-, 12-, and 18-week-old) identified transcriptional ageing fingerprints in the killifish retina, unveiling pathways also identified in the aged brain, including oxidative stress, gliosis, and inflammageing. These findings were comparable to observations in ageing mouse retina. Additionally, transcriptional changes in genes related to retinal diseases, such as glaucoma and age-related macular degeneration, were observed. The cellular heterogeneity in the killifish retina was characterised, confirming the presence of all typical vertebrate retinal cell types. Data integration from age-matched samples between the bulk and scRNA-seq experiments revealed a loss of cellular specificity in gene expression upon ageing, suggesting potential disruption in transcriptional homeostasis. Differential expression analysis within the identified cell types highlighted the role of glial/immune cells as important stress regulators during ageing. Our work emphasises the value of the fast-ageing killifish in elucidating molecular signatures in age-associated retinal disease and vision decline. This study contributes to the understanding of how age-related changes in molecular pathways may impact CNS health, providing insights that may inform future therapeutic strategies for age-related pathologies.

Comparison of retinal parameters between rhesus and cynomolgus macaques

Nonhuman primates are important research models for basic vision research, preclinical pathogenesis, and treatment studies due to strong similarities in retinal structure and function with humans. We compared retinal parameters between 10 healthy normal rhesus macaques (Macaca mulatta) and 10 cynomolgus macaques (Macaca fascicularis) by optical coherence tomography and electroretinography. The Heidelberg Spectralis® HRA+OCT and Roland multifocal electrophysiometer were used to analyze retinal morphology, multifocal electroretinograms (mfERGs), and full-field electroretinograms (ff-ERGs). Mean retinal thickness was lowest in the central fovea of macaques and did not differ significantly between species, but the retinal thicknesses of the nerve fiber ganglion cell layer and the inner plexiform layer were significantly different. The amplitude density of the N1 wave was lower in rhesus macaques than in cynomolgus macaques in ring and quadrant areas. Dark-adapted 3.0 oscillatory potentials (reflection of amacrine cell activity) and light-adapted 30-hz flicker ERG (a sensitive cone-pathway-driven response) waveforms of the ff-ERG were similar in both species, while the times to peaks in dark-adapted 0.01 ERG (the rod-driven response of bipolar cells) and dark-adapted 3.0 ERG (combined rod and cone system responses) as well as the implicit times of the a- and b-waves in light-adapted 3.0 ERG (the single-flash cone response) were substantially different. This study provides normative retinal parameters for nonhuman primate research on basic and clinical ophthalmology, as well as a reference for researchers in the appropriate selection of rhesus or cynomolgus macaques as models for ophthalmology studies.

Artificial Intelligence Models for Cell Type and Subtype Identification Based on Single-Cell RNA Sequencing Data in Vision Science

Single-cell RNA sequencing (scRNA-seq) provides a high throughput, quantitative and unbiased framework for scientists in many research fields to identify and characterize cell types within heterogeneous cell populations from various tissues. However, scRNA-seq based identification of discrete cell-types is still labor intensive and depends on prior molecular knowledge. Artificial intelligence has provided faster, more accurate, and user-friendly approaches for cell-type identification. In this review, we discuss recent advances in cell-type identification methods using artificial intelligence techniques based on single-cell and single-nucleus RNA sequencing data in vision science. The main purpose of this review paper is to assist vision scientists not only to select suitable datasets for their problems, but also to be aware of the appropriate computational tools to perform their analysis. Developing novel methods for scRNA-seq data analysis remains to be addressed in future studies.

Single-cell RNA sequencing in dissecting microenvironment of age-related macular degeneration: Challenges and perspectives

Age-related macular degeneration (AMD) is the leading cause of blindness in individuals over the age of 50 years, yet its etiology and pathogenesis largely remain uncovered. Single-cell RNA sequencing (scRNA-seq) technologies are recently developed and have a number of advantages over conventional bulk RNA sequencing techniques in uncovering the heterogeneity of complex microenvironments containing numerous cell types and cell communications during various biological processes. In this review, we summarize the latest discovered cellular components and regulatory mechanisms during AMD development revealed by scRNA-seq. In addition, we discuss the main challenges and future directions in exploring the pathophysiology of AMD equipped with single-cell technologies. Our review underscores the importance of multimodal single-cell platforms (such as single-cell spatiotemporal multi-omics and single-cell exosome omics) as new approaches for basic and clinical AMD research in identifying biomarker, characterizing cellular responses to drug treatment and environmental stimulation.

Single cell transcriptome analyses of the developing zebrafish eye- perspectives and applications

Within a relatively short period of time, single cell transcriptome analyses (SCT) have become increasingly ubiquitous with transcriptomic research, uncovering plentiful details that boost our molecular understanding of various biological processes. Stemming from SCT analyses, the ever-growing number of newly assigned genetic markers increases our understanding of general function and development, while providing opportunities for identifying genes associated with disease. SCT analyses have been carried out using tissue from numerous organisms. However, despite the great potential of zebrafish as a model organism, other models are still preferably used. In this mini review, we focus on eye research as an example of the advantages in using zebrafish, particularly its usefulness for single cell transcriptome analyses of developmental processes. As studies have already shown, the unique opportunities offered by zebrafish, including similarities to the human eye, in combination with the possibility to analyze and extract specific cells at distinct developmental time points makes the model a uniquely powerful one. Particularly the practicality of collecting large numbers of embryos and therefore isolation of sufficient numbers of developing cells is a distinct advantage compared to other model organisms. Lastly, the advent of highly efficient genetic knockouts methods offers opportunities to characterize target gene function in a more cost-efficient way. In conclusion, we argue that the use of zebrafish for SCT approaches has great potential to further deepen our molecular understanding of not only eye development, but also many other organ systems.

Spatial distribution of metabolites in the retina and its relevance to studies of metabolic retinal disorders

INTRODUCTION

The primate retina has evolved regional specialisations for specific visual functions. The macula is specialised towards high acuity vision and is an area that contains an increased density of cone photoreceptors and signal processing neurons. Different regions in the retina display unique susceptibility to pathology, with many retinal diseases primarily affecting the macula.

OBJECTIVES

To better understand the properties of different retinal areas we studied the differential distribution of metabolites across the retina.

METHODS

We conducted an untargeted metabolomics analysis on full-thickness punches from three different regions (macula, temporal peri-macula and periphery) of healthy primate retina.

RESULTS

Nearly half of all metabolites identified showed differential abundance in at least one comparison between the three regions. Furthermore, mapping metabolomics results from macula-specific eye diseases onto our region-specific metabolite distributions revealed differential abundance defining systemic metabolic dysregulations that were region specific.

CONCLUSIONS

The unique metabolic phenotype of different retinal regions is likely due to the differential distribution of different cell types in these regions reflecting the specific metabolic requirements of each cell type. Our results may help to better understand the pathobiology of retinal diseases with region specificity.

HBEGF-TNF induce a complex outer retinal pathology with photoreceptor cell extrusion in human organoids

Human organoids could facilitate research of complex and currently incurable neuropathologies, such as age-related macular degeneration (AMD) which causes blindness. Here, we establish a human retinal organoid system reproducing several parameters of the human retina, including some within the macula, to model a complex combination of photoreceptor and glial pathologies. We show that combined application of TNF and HBEGF, factors associated with neuropathologies, is sufficient to induce photoreceptor degeneration, glial pathologies, dyslamination, and scar formation: These develop simultaneously and progressively as one complex phenotype. Histologic, transcriptome, live-imaging, and mechanistic studies reveal a previously unknown pathomechanism: Photoreceptor neurodegeneration via cell extrusion. This could be relevant for aging, AMD, and some inherited diseases. Pharmacological inhibitors of the mechanosensor PIEZO1, MAPK, and actomyosin each avert pathogenesis; a PIEZO1 activator induces photoreceptor extrusion. Our model offers mechanistic insights, hypotheses for neuropathologies, and it could be used to develop therapies to prevent vision loss or to regenerate the retina in patients suffering from AMD and other diseases.

No evidence for age-related alterations in the marmoset retina

The physiological aging process of the retina is accompanied by various and sometimes extensive changes: Macular degeneration, retinopathies and glaucoma are the most common findings in the elderly and can potentially lead to irreversible visual disablements up to blindness. To study the aging process and to identify possible therapeutic targets to counteract these diseases, the use of appropriate animal models is mandatory. Besides the most commonly used rodent species, a non-human primate, the common marmoset (Callithrix jacchus) emerged as a promising animal model of human aging over the last years. However, the visual aging process in this species is only partially characterized, especially with regard to retinal aberrations. Therefore, we assessed here for the first time potential changes in retinal morphology of the common marmoset of different age groups. By cell type specific immunolabeling, we analyzed different cell types and distributions, potential photoreceptor and ganglion cell loss, and structural reorganization. We detected no signs of age-related differences in staining patterns or densities of various cell populations. For example, there were no signs of photoreceptor degeneration, and there was only minimal sprouting of rod bipolar cells in aged retinas. Altogether, we describe here the maintenance of a stable neuronal architecture, distribution and number of different cell populations with only mild aberrations during the aging process in the common marmoset retina. These findings are in stark contrast to previously reported findings in rodent species and humans and deserve further investigations to identify the underlying mechanisms and possible therapeutic targets.

A single-nucleus transcriptomic atlas of the dog hippocampus reveals the potential relationship between specific cell types and domestication

The process of domestication has led to dramatic differences in behavioral traits between domestic dogs and gray wolves. Whole-genome research found that a class of putative positively selected genes were related to various aspects of learning and memory, such as long-term potentiation and long-term depression. In this study, we constructed a single-nucleus transcriptomic atlas of the dog hippocampus to illustrate its cell types, cell lineage and molecular features. Using the transcriptomes of 105 057 nuclei from the hippocampus of a Beagle dog, we identified 26 cell clusters and a putative trajectory of oligodendrocyte development. Comparative analysis revealed a significant convergence between dog differentially expressed genes (DEGs) and putative positively selected genes (PSGs). Forty putative PSGs were DEGs in glutamatergic neurons, especially in Cluster 14, which is related to the regulation of nervous system development. In summary, this study provides a blueprint to understand the cellular mechanism of dog domestication.

Revival of light signalling in the postmortem mouse and human retina

Death is defined as the irreversible cessation of circulatory, respiratory or brain activity. Many peripheral human organs can be transplanted from deceased donors using protocols to optimize viability. However, tissues from the central nervous system rapidly lose viability after circulation ceases^1,2, impeding their potential for transplantation. The time course and mechanisms causing neuronal death and the potential for revival remain poorly defined. Here, using the retina as a model of the central nervous system, we systemically examine the kinetics of death and neuronal revival. We demonstrate the swift decline of neuronal signalling and identify conditions for reviving synchronous in vivo-like trans-synaptic transmission in postmortem mouse and human retina. We measure light-evoked responses in human macular photoreceptors in eyes removed up to 5 h after death and identify modifiable factors that drive reversible and irreversible loss of light signalling after death. Finally, we quantify the rate-limiting deactivation reaction of phototransduction, a model G protein signalling cascade, in peripheral and macular human and macaque retina. Our approach will have broad applications and impact by enabling transformative studies in the human central nervous system, raising questions about the irreversibility of neuronal cell death, and providing new avenues for visual rehabilitation.

Molecular and cytological profiling of biological aging of mouse cochlear inner and outer hair cells

Age-related hearing loss (ARHL) negatively impacts quality of life in the elderly population. The prevalent cause of ARHL is loss of mechanosensitive cochlear hair cells (HCs). The molecular and cellular mechanisms of HC degeneration remain poorly understood. Using RNA-seq transcriptomic analyses of inner and outer HCs isolated from young and aged mice, we show that HC aging is associated with changes in key molecular processes, including transcription, DNA damage, autophagy, and oxidative stress, as well as genes related to HC specialization. At the cellular level, HC aging is characterized by loss of stereocilia, shrinkage of HC soma, and reduction in outer HC mechanical properties, suggesting that functional decline in mechanotransduction and cochlear amplification precedes HC loss and contributes to ARHL. Our study reveals molecular and cytological profiles of aging HCs and identifies genes such as Sod1, Sirt6, Jund, and Cbx3 as biomarkers and potential therapeutic targets for ameliorating ARHL.

Using RNA-seq, advanced imaging, and electrophysiology, Liu et al. reveal molecular and cytological profiles of aging cochlear hair cells. Their study also suggests that a functional decline in mechanotransduction and cochlear amplification precedes hair cell loss and contributes to age-related hearing loss.

Microglia: Key Players in Retinal Ageing and Neurodegeneration

Microglia are the resident immune cells of the central nervous system (CNS) and play a key role in maintaining the normal function of the retina and brain. During early development, microglia migrate into the retina, transform into a highly ramified phenotype, and scan their environment constantly. Microglia can be activated by any homeostatic disturbance that may endanger neurons and threaten tissue integrity. Once activated, the young microglia exhibit a high diversity in their phenotypes as well as their functions, which relate to either beneficial or harmful consequences. Microglial activation is associated with the release of cytokines, chemokines, and growth factors that can determine pathological outcomes. As the professional phagocytes in the retina, microglia are responsible for the clearance of pathogens, dead cells, and protein aggregates. However, their phenotypic diversity and phagocytic capacity is compromised with ageing. This may result in the accumulation of protein aggregates and myelin debris leading to retinal neuroinflammation and neurodegeneration. In this review, we describe microglial phenotypes and functions in the context of the young and ageing retina, and the mechanisms underlying changes in ageing. Additionally, we review microglia-mediated retinal neuroinflammation and discuss the mechanisms of microglial involvement in retinal neurodegenerative diseases.

Synaptic or Non-synaptic? Different Intercellular Interactions with Retinal Ganglion Cells in Optic Nerve Regeneration

Axons of adult neurons in the mammalian central nervous system generally fail to regenerate by themselves, and few if any therapeutic options exist to reverse this situation. Due to a weak intrinsic potential for axon growth and the presence of strong extrinsic inhibitors, retinal ganglion cells (RGCs) cannot regenerate their axons spontaneously after optic nerve injury and eventually undergo apoptosis, resulting in permanent visual dysfunction. Regarding the extracellular environment, research to date has generally focused on glial cells and inflammatory cells, while few studies have discussed the potentially significant role of interneurons that make direct connections with RGCs as part of the complex retinal circuitry. In this study, we provide a novel angle to summarize these extracellular influences following optic nerve injury as "intercellular interactions" with RGCs and classify these interactions as synaptic and non-synaptic. By discussing current knowledge of non-synaptic (glial cells and inflammatory cells) and synaptic (mostly amacrine cells and bipolar cells) interactions, we hope to accentuate the previously neglected but significant effects of pre-synaptic interneurons and bring unique insights into future pursuit of optic nerve regeneration and visual function recovery.

Variability in Retinal Neuron Populations and Associated Variations in Mass Transport Systems of the Retina in Health and Aging

Aging is associated with a broad range of visual impairments that can have dramatic consequences on the quality of life of those impacted. These changes are driven by a complex series of alterations affecting interactions between multiple cellular and extracellular elements. The resilience of many of these interactions may be key to minimal loss of visual function in aging; yet many of them remain poorly understood. In this review, we focus on the relation between retinal neurons and their respective mass transport systems. These metabolite delivery systems include the retinal vasculature, which lies within the inner portion of the retina, and the choroidal vasculature located externally to the retinal tissue. A framework for investigation is proposed and applied to identify the structures and processes determining retinal mass transport at the cellular and tissue levels. Spatial variability in the structure of the retina and changes observed in aging are then harnessed to explore the relation between variations in neuron populations and those seen among retinal metabolite delivery systems. Existing data demonstrate that the relation between inner retinal neurons and their mass transport systems is different in nature from that observed between the outer retina and choroid. The most prominent structural changes observed across the eye and in aging are seen in Bruch’s membrane, which forms a selective barrier to mass transfers at the interface between the choroidal vasculature and the outer retina.

Generating and Using Transcriptomically Based Retinal Cell Atlases

It has been known for over a century that the basic organization of the retina is conserved across vertebrates. It has been equally clear that retinal cells can be classified into numerous types, but only recently have methods been devised to explore this diversity in unbiased, scalable, and comprehensive ways. Advances in high-throughput single-cell RNA-sequencing (scRNA-seq) have played a pivotal role in this effort. In this article, we outline the experimental and computational components of scRNA-seq and review studies that have used them to generate retinal atlases of cell types in several vertebrate species. These atlases have enabled studies of retinal development, responses of retinal cells to injury, expression patterns of genes implicated in retinal disease, and the evolution of cell types. Recently, the inquiry has expanded to include the entire eye and visual centers in the brain. These studies have enhanced our understanding of retinal function and dysfunction and provided tools and insights for exploring neural diversity throughout the brain. Expected final online publication date for the Annual Review of Vision Science, Volume 7 is September 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Antioxidant Defenses in the Human Eye: A Focus on Metallothioneins

The human eye, the highly specialized organ of vision, is greatly influenced by oxidants of endogenous and exogenous origin. Oxidative stress affects all structures of the human eye with special emphasis on the ocular surface, the lens, the retina and its retinal pigment epithelium, which are considered natural barriers of antioxidant protection, contributing to the onset and/or progression of eye diseases. These ocular structures contain a complex antioxidant defense system slightly different along the eye depending on cell tissue. In addition to widely studied enzymatic antioxidants, including superoxide dismutase, glutathione peroxidase, catalase, peroxiredoxins and selenoproteins, inter alia, metallothioneins (MTs) are considered antioxidant proteins of growing interest with further cell-mediated functions. This family of cysteine rich and low molecular mass proteins captures and neutralizes free radicals in a redox-dependent mechanism involving zinc binding and release. The state of the art of MTs, including the isoforms classification, the main functions described to date, the Zn-MT redox cycle as antioxidant defense system, and the antioxidant activity of Zn-MTs in the ocular surface, lens, retina and its retinal pigment epithelium, dependent on the number of occupied zinc-binding sites, will be comprehensively reviewed.