Combining Gene-Disease Associations with Single-Cell Gene Expression Data Provides Anatomy-Specific Subnetworks in...

Background: Age-related macular degeneration (AMD) is the most common cause of visual impairment in the developed world. Despite some treatment options for late AMD, there is no intervention that blocks early AMD proceeding to the late and blinding forms. This is partly due to the lack of precise drug targets, despite great advances in genetics, epidemiology, and protein-protein interaction (PPI) networks proposed to be driving the disease pathology. A systems approach to narrow down PPI networks to specific protein drug targets would provide new therapeutic options. Materials and Methods: In this study we analyzed single cell RNAseq (RNA sequencing) datasets of 17 cell types present in choroidal, retinal pigment epithelium (RPE), and neural retina (NR) tissues to explore if a more granular analysis incorporating different cell types exposes more specific pathways and relationships. Furthermore, we developed a novel and systematic gene ontology database (SysGO) to explore if a subcellular classification of processes will further enhance the understanding of the pathogenesis of this complex disorder and its comorbidities with other age-related diseases. Results: We found that 57% of the AMD (risk) genes are among the top 25% expressed genes in ∼1 of the 17 choroidal/RPE/NR cell types, and 9% were among the top 1% of expressed genes. Using SysGO, we identified an enrichment of AMD genes in cell membrane and extracellular anatomical locations, and we found both functional enrichments (e.g., cell adhesion) and cell types (e.g., fibroblasts, microglia) not previously associated with AMD pathogenesis. We reconstructed PPI networks among the top expressed AMD genes for all 17 choroidal/RPE/NR cell types, which provides molecular and anatomical definitions of AMD phenotypes that can guide therapeutic approaches to target this complex disease. Conclusion: We provide mechanism-based AMD endophenotypes that can be exploited in vitro, using computational models and for drug discovery/repurposing.