Abstract: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.

Genetic analyses of human fetal retinal pigment epithelium gene expression suggest ocular disease mechanisms

Dissecting the Transcriptome Landscape of the Human Fetal Neural Retina and Retinal Pigment Epithelium by Single-Cell RNA-seq Analysis

Integrative Single‐Cell Transcriptomics and Epigenomics Mapping of the Fetal Retina Developmental Dynamics

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

RNASeq profiling of retinal pigment epithelial cells derived from induced pluripotent stem cells revealed 3 genes involved in lipid homeostasis in age-related macular degeneration

Single-cell transcriptomics of the human retinal pigment epithelium and choroid in health and macular degeneration

A Single-Cell Transcriptome Atlas of the Human Retinal Pigment Epithelium

Single-Cell Analysis of Human Retina Identifies Evolutionarily Conserved and Species-Specific Mechanisms Controlling Development

( ) is a critical early regulator of retinal pigment epithelial development

Insights into the genetic basis of retinal detachment

Multi-omic Analysis of Developing Human Retina and Organoids Reveals Cell-Specific Cis-Regulatory Elements and Mechanisms of Non-Coding Genetic Disease Risk

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

Genetic Analysis Of Families With Retinal Dystrophies

Identification of Gene Mutations in Atypical Retinopathy of Prematurity Cases.

Whole Genome Expression Profiling of Normal Human Fetal and Adult Ocular Tissues.

RPE65 Mutation Frequency and Phenotypic Variation According to Exome Sequencing in a Tertiary Centre for Genetic Eye Diseases in China

Genetic Risk of Reticular Pseudodrusen in Age-Related Macular Degeneration: HTRA1/lncRNA BX842242.1 dominates, with no evidence for Complement Cascade involvement

Pigmentation and retinal pigment epithelium thickness: a study of the phenotypic and genotypic relationships between ocular and extraocular pigmented tissues.

Single-cell RNA-seq analysis maps the development of human fetal retina

High-resolution genome topology of human retina uncovers super enhancer-promoter interactions at tissue-specific and multifactorial disease loci

Identification of Human Retinal Organoid Cell Differentiation-Related Genes via Single-Cell Sequencing Data Analysis