Deep Sc-RNA sequencing decoding the molecular dynamic architecture of the human retina

Single-cell sequencing in diabetic retinopathy: progress and prospects

Diabetic retinopathy is a major ocular complication of diabetes, characterized by progressive retinal microvascular damage and significant visual impairment in working-age adults. Traditional bulk RNA sequencing offers overall gene expression profiles but does not account for cellular heterogeneity. Single-cell RNA sequencing overcomes this limitation by providing transcriptomic data at the individual cell level and distinguishing novel cell subtypes, developmental trajectories, and intercellular communications. Researchers can use single-cell sequencing to draw retinal cell atlases and identify the transcriptomic features of retinal cells, enhancing our understanding of the pathogenesis and pathological changes in diabetic retinopathy. Additionally, single-cell sequencing is widely employed to analyze retinal organoids and single extracellular vesicles. Single-cell multi-omics sequencing integrates omics information, whereas stereo-sequencing analyzes gene expression and spatiotemporal data simultaneously. This review discusses the protocols of single-cell sequencing for obtaining single cells from retina and accurate sequencing data. It highlights the applications and advancements of single-cell sequencing in the study of normal retinas and the pathological changes associated with diabetic retinopathy. This underscores the potential of these technologies to deepen our understanding of the pathogenesis of diabetic retinopathy that may lead to the introduction of new therapeutic strategies.

A Comprehensive Analysis of Sex-Biased Gene Expression in the Aging Human Retina Through a Combination of Single-Cell and Bulk RNA Sequencing

Purpose

Previous studies have reported divergent sexual responses to aging; however, specific variations in gene expression between aging males and females and their potential association with age-related retinal diseases remain unclear. This study collected data from public databases and developed a comprehensive comparison of retina between aging females and males.

Methods

Single-cell RNA (scRNA) and bulk RNA sequencing data of the aging retina from females and males in public databases were utilized for integrated analysis to investigate sex-biased expression in retina. Additionally, in vitro experiments were conducted on individuals with retinitis pigmentosa (RP) to validate the sex difference in degenerative retina.

Results

Bulk RNA analysis revealed sex-biased expression of specific genes in retina of aging individuals, with immune pathway-related genes exhibiting higher expression in females compared to males. The scRNA analysis demonstrated that sex-biased gene expression was cell-type specific in aging retina. Furthermore, susceptibility genes for age-related macular degeneration and RP exhibited variation across different cell types and sexes. Cell-to-cell communication unveiled an increased interaction associated with TGFB1, CCL7, and VEGFA in Müller glia, microglia, and astrocytes of female retina. Notably, we observed female-biased chemokine expression in microglia contributing to heightened susceptibility to immune inflammation in female retina. Finally, we confirmed a more pronounced inflammatory response during degeneration in female rd10 mouse retina compared to males.

Conclusions

This study provides a comprehensive comparison of retina between females and males in healthy aging human retina and highlights the significance of sex as an influential factor in retinal diseases.

A genome-wide scan of non-coding RNAs and enhancers for refractive error and myopia

Refractive error (RE) and myopia are complex polygenic conditions with the majority of genome-wide associated genetic variants in non-exonic regions. Given this, and the onset during childhood, gene-regulation is expected to play an important role in its pathogenesis. This prompted us to explore beyond traditional gene finding approaches. We performed a genetic association study between variants in non-coding RNAs and enhancers, and RE and myopia. We obtained single-nucleotide polymorphisms (SNPs) in microRNA (miRNA) genes, miRNA-binding sites, long non-coding RNAs genes (lncRNAs) and enhancers from publicly available databases: miRNASNPv2, PolymiRTS, VISTA Enhancer Browser, FANTOM5 and lncRNASNP2. We investigated whether SNPs overlapping these elements were associated with RE and myopia leveraged from a large GWAS meta-analysis (N = 160,420). With genetic risk scores (GRSs) per element, we investigated the joint effect of associated variants on RE, axial length (AL)/corneal radius (CR), and AL progression in an independent child cohort, the Generation R Study (N = 3638 children). We constructed a score for biological plausibility per SNP in highly confident miRNA-binding sites and enhancers in chromatin accessible regions. We found that SNPs in two miRNA genes, 14 enhancers and 81 lncRNA genes in chromatin accessible regions and 54 highly confident miRNA-binding sites, were in RE and myopia-associated loci. GRSs from SNPs in enhancers were significantly associated with RE, AL/CR and AL progression. GRSs from lncRNAs were significantly associated with all AL/CR and AL progression. GRSs from miRNAs were not associated with any ocular biometric measurement. GRSs from miRNA-binding sites showed suggestive but inconsistent significance. We prioritized candidate miRNA binding sites and candidate enhancers for future functional validation. Pathways of target and host genes of highly ranked variants included eye development (BMP4, MPPED2), neurogenesis (DDIT4, NTM), extracellular matrix (ANTXR2, BMP3), photoreceptor metabolism (DNAJB12), photoreceptor morphogenesis (CHDR1), neural signaling (VIPR2) and TGF-beta signaling (ANAPC16). This is the first large-scale study of non-coding RNAs and enhancers for RE and myopia. Enhancers and lncRNAs could be of large importance as they are associated with childhood myopia. We provide a confident blueprint for future functional validation by prioritizing candidate miRNA binding sites and candidate enhancers.

Multi-omics in exploring the pathophysiology of diabetic retinopathy

Diabetic retinopathy (DR) is a leading global cause of vision impairment, with its prevalence increasing alongside the rising rates of diabetes mellitus (DM). Despite the retina's complex structure, the underlying pathology of DR remains incompletely understood. Single-cell RNA sequencing (scRNA-seq) and recent advancements in multi-omics analyses have revolutionized molecular profiling, enabling high-throughput analysis and comprehensive characterization of complex biological systems. This review highlights the significant contributions of scRNA-seq, in conjunction with other multi-omics technologies, to DR research. Integrated scRNA-seq and transcriptomic analyses have revealed novel insights into DR pathogenesis, including alternative transcription start site events, fluctuations in cell populations, altered gene expression profiles, and critical signaling pathways within retinal cells. Furthermore, by integrating scRNA-seq with genetic association studies and multi-omics analyses, researchers have identified novel biomarkers, susceptibility genes, and potential therapeutic targets for DR, emphasizing the importance of specific retinal cell types in disease progression. The integration of scRNA-seq with metabolomics has also been instrumental in identifying specific metabolites and dysregulated pathways associated with DR. It is highly conceivable that the continued synergy between scRNA-seq and other multi-omics approaches will accelerate the discovery of underlying mechanisms and the development of novel therapeutic interventions for DR.

Identification of LRRC46 as a novel candidate gene for high myopia

High myopia (HM) is the primary cause of blindness, with the microstructural organization and composition of collagenous fibers in the cornea and sclera playing a crucial role in the biomechanical behavior of these tissues. In a previously reported myopic linkage region, MYP5 (17q21-22), a potential candidate gene, LRRC46 (c.C235T, p.Q79X), was identified in a large Han Chinese pedigree. LRRC46 is expressed in various eye tissues in humans and mice, including the retina, cornea, and sclera. In subsequent cell experiments, the mutation (c.C235T) decreased the expression of LRRC46 protein in human corneal epithelial cells (HCE-T). Further investigation revealed that Lrrc46-/- mice (KO) exhibited a classical myopia phenotype. The thickness of the cornea and sclera in KO mice became thinner and more pronounced with age, the activity of limbal stem cells decreased, and microstructural changes were observed in the fibroblasts of the sclera and cornea. We performed RNA-seq on scleral and corneal tissues of KO and normal control wild-type (WT) mice, which indicated a significant downregulation of the collagen synthesis-related pathway (extracellular matrix, ECM) in KO mice. Subsequent in vitro studies further indicated that LRRC46, a member of the important LRR protein family, primarily affected the formation of collagens. This study suggested that LRRC46 is a novel candidate gene for HM, influencing collagen protein VIII (Col8a1) formation in the eye and gradually altering the biomechanical structure of the cornea and sclera, thereby promoting the occurrence and development of HM.

Senolytic and senomorphic agent procyanidin C1 alleviates structural and functional decline in the aged retina

Increased cellular senescence burden contributes in part to age-related organ dysfunction and pathologies. In our study, using mouse models of natural aging, we observed structural and functional decline in the aged retina, which was accompanied by the accumulation of senescent cells and senescence-associated secretory phenotype factors. We further validated the senolytic and senomorphic properties of procyanidin C1 (PCC1) both in vitro and in vivo, the long-term treatment of which ameliorated age-related retinal impairment. Through high-throughput single-cell RNA sequencing (scRNA-seq), we comprehensively characterized the retinal landscape after PCC1 administration and deciphered the molecular basis underlying the senescence burden increment and elimination. By exploring the scRNA-seq database of age-related retinal disorders, we revealed the role of cellular senescence and the therapeutic potential of PCC1 in these pathologies. Overall, these results indicate the therapeutic effects of PCC1 on the aged retina and its potential use for treating age-related retinal disorders.

RNA fusion in human retinal development

Chimeric RNAs have been found in both cancerous and healthy human cells. They have regulatory effects on human stem/progenitor cell differentiation, stemness maintenance, and central nervous system development. However, whether they are present in human retinal cells and their physiological functions in the retinal development remain unknown. Based on the human embryonic stem cell-derived retinal organoids (ROs) spanning from days 0 to 120, we present the expression atlas of chimeric RNAs throughout the developing ROs. We confirmed the existence of some common chimeric RNAs and also discovered many novel chimeric RNAs during retinal development. We focused on CTNNBIP1-CLSTN1 (CTCL) whose downregulation caused precocious neuronal differentiation and a marked reduction of neural progenitors in human cerebral organoids. CTCL is universally present in human retinas, ROs, and retinal cell lines, and its loss-of-function biases the progenitor cells toward retinal pigment epithelial cell fate at the expense of retinal cells. Together, this work provides a landscape of chimeric RNAs and reveals evidence for their critical role in human retinal development.

Identification and validation of key biomarkers and potential therapeutic targets for primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is a prevalent cause of blindness worldwide, resulting in degeneration of retinal ganglion cells and permanent damage to the optic nerve. However, the underlying pathogenetic mechanisms of POAG are currently indistinct, and there has been no effective nonsurgical treatment regimen. The objective of this study is to identify novel biomarkers and potential therapeutic targets for POAG. The mRNA expression microarray datasets GSE27276 and GSE138125, as well as the single-cell high-throughput RNA sequencing (scRNA-seq) dataset GSE148371 were utilized to screen POAG-related differentially expressed genes (DEGs). Functional enrichment analyses, protein-protein interaction (PPI) analysis, and weighted gene co-expression network analysis (WGCNA) of the DEGs were performed. Subsequently, the hub genes were validated at a single-cell level, where trabecular cells were annotated, and the mRNA expression levels of target genes in different cell clusters were analyzed. Immunofluorescence and quantitative real-time PCR (qPCR) were performed for further validation. DEGs analysis identified 43 downregulated and 32 upregulated genes in POAG, which were mainly enriched in immune-related pathways, oxidative stress, and endoplasmic reticulum (ER) stress. PPI networks showed that FN1 and DUSP1 were the central hub nodes, while GPX3 and VAV3 were screened out as hub genes through WGCNA and subsequently validated by qPCR. Finally, FN1, GPX3, and VAV3 were determined to be pivotal core genes via single-cell validation. The relevant biomarkers involved in the pathogenesis of POAG, may serve as potential therapeutic targets. Further studies are necessary to unveil the mechanisms underlying the expression variations of these genes in POAG.

Dynamic human retinal pigment epithelium (RPE) and choroid architecture based on single-cell transcriptomic landscape analysis

The retinal pigment epithelium (RPE) and choroid are located behind the human retina and have multiple functions in the human visual system. Knowledge of the RPE and choroid cells and their gene expression profiles are fundamental for understanding retinal disease mechanisms and therapeutic strategies. Here, we sequenced the RNA of about 0.3 million single cells from human RPE and choroids across two regions and seven ages, revealing regional and age differences within the human RPE and choroid. Cell-cell interactions highlight the broad connectivity networks between the RPE and different choroid cell types. Moreover, the transcription factors and their target genes change during aging. The coding of somatic variations increases during aging in the human RPE and choroid at the single-cell level. Moreover, we identified ELN as a candidate for improving RPE degeneration and choroidal structure during aging. The mapping of the molecular architecture of the human RPE and choroid improves our understanding of the human vision support system and offers potential insights into the intervention targets for retinal diseases.

Single-cell transcriptome reveals diversity of Müller cells with different metabolic-mitochondrial signatures in normal and degenerated macula

Müller cell is the most abundant glial cell in mammalian retina, supporting the functions of photoreceptors and other retinal neurons via maintaining environmental homeostasis. In response to injury and/or neuronal degeneration, Müller cells undergo morphological and functional alternations, known as reactive gliosis documented in multiple retinal diseases, including age-related macular degeneration (AMD), retinitis pigmentosa, diabetic retinopathy, and traumatic retinal detachment. But the functional consequences of Müller glia cell reactivation or even the regulatory networks of the retinal gliosis are still controversial. In this study, we reveal different subpopulations of Müller cells with distinct metabolic-mitochondrial signatures by integrating single cell transcriptomic data from Early AMD patients and healthy donors. Our results show that a portion of Müller cells exhibits low mitochondrial DNA (mtDNA) expressions, reduced protein synthesis, impaired homeostatic regulation, decreased proliferative ability but enhanced proangiogenic function. Interestingly, the major alternation of Müller cells in Early AMD retina is the change of subpopulation abundance, rather than generation of new subcluster. Transcription factor enrichment analysis further highlights the key regulators of metabolic-mitochondrial states of Müller glias in Early AMD patients especially. Our study demonstrates new characteristics of retinal gliosis associated with Early AMD and suggests the possibility to prevent degeneration by intervening mitochondrial functions of Müller cells.