Integration of eQTL and a Single-Cell Atlas in the Human Eye Identifies Causal Genes for Age-Related Macular Degeneration

Decreased Circulating Very Small Low-Density Lipoprotein is Likely Causal for Age-Related Macular Degeneration

Objective

Abnormal changes in metabolite levels in serum or plasma have been highlighted in several studies in age-related macular degeneration (AMD), the leading cause of irreversible vision loss. Specific changes in lipid profiles are associated with an increased risk of AMD. Metabolites could thus be used to investigate AMD disease mechanisms or incorporated into AMD risk prediction models. However, whether particular metabolites causally affect the disease has yet to be established.

Design

A 3-tiered analysis of blood metabolites in the United Kingdom (UK) Biobank cohort to identify metabolites that differ in AMD patients with evidence for a putatively causal role in AMD.

Participants

A total of 72 376 donors from the UK Biobank cohort including participants with AMD (N = 1353) and non-AMD controls (N = 71 023).

Methods

We analyzed 325 directly measured or derived blood metabolites from the UK Biobank for 72 376 donors to identify AMD-associated metabolites. Genome-wide association studies for 325 metabolites in 98 316 European participants from the UK Biobank were performed. The causal effects of these metabolites in AMD were tested using a 2-sample Mendelian randomization approach. The predictive value of these measurements together with sex and age was assessed by developing a machine learning classifier.

Main Outcome Measures

Evaluating metabolic biomarkers associated with AMD susceptibility and investigating their potential causal contribution to the development of the disease.

Results

This study noted age to be the prominent risk factor associated with AMD development. While accounting for age and sex, we identified 84 metabolic markers as significantly (false discovery rate-adjusted P value < 0.05) associated with AMD. Lipoprotein subclasses comprised the majority of the AMD-associated metabolites (39%) followed by several lipoprotein to lipid ratios. Nineteen metabolites showed a likely causative role in AMD etiology. Of these, 6 lipoproteins contain very small, very low-density lipoprotein (VLDL), and phospholipids to total lipid ratio in medium VLDL. Based on this we postulate that depletion of circulating very small VLDLs is likely causal for AMD. The risk prediction model constructed from the metabolites, age and sex, identified age as the primary predictive factor with a much smaller contribution by metabolites to AMD risk prediction.

Conclusions

This study underscores the pronounced role of lipids in AMD susceptibility and the likely causal contribution of particular subclasses of lipoproteins to AMD. Our study provides valuable insights into the metabopathological mechanisms of AMD disease development and progression.

Integrating Multi-omics to Identify Age-Related Macular Degeneration Subtypes and Biomarkers

Age-related macular degeneration (AMD) is one of the most common causes of irreversible vision loss in the elderly. Its pathogenesis is likely multifactorial, involving a complex interaction of metabolic and environmental factors, and remains poorly understood. Previous studies have shown that mitochondrial dysfunction and oxidative stress play a crucial role in the development of AMD. Oxidative damage to the retinal pigment epithelium (RPE) has been identified as one of the major mediators in the pathogenesis of age-related macular degeneration (AMD). Therefore, this article combines transcriptome sequencing (RNA-seq) and single-cell sequencing (scRNA-seq) data to explore the role of mitochondria-related genes (MRGs) in AMD. Firstly, differential expression analysis was performed on the raw RNA-seq data. The intersection of differentially expressed genes (DEGs) and MRGs was performed. This paper proposes a deep subspace nonnegative matrix factorization (DS-NMF) algorithm to perform a multi-layer nonlinear transformation on the intersection of gene expression profiles corresponding to AMD samples. The age of AMD patients is used as prior information at the network's top level to change the data distribution. The classification is based on reconstructed data with altered distribution. The types obtained significantly differ in scores of multiple immune-related pathways and immune cell infiltration abundance. Secondly, an optimal AMD diagnosis model was constructed using multiple machine learning algorithms for external and qRT-PCR verification. Finally, ten potential therapeutic drugs for AMD were identified based on cMAP analysis. The AMD subtypes identified in this article and the diagnostic model constructed can provide a reference for treating AMD and discovering new drug targets.

Aging of the eye: Lessons from cataracts and age-related macular degeneration

Aging is the greatest risk factor for chronic human diseases, including many eye diseases. Geroscience aims to understand the effects of the aging process on these diseases, including the genetic, molecular, and cellular mechanisms that underlie the increased risk of disease over the lifetime. Understanding of the aging eye increases general knowledge of the cellular physiology impacted by aging processes at various biological extremes. Two major diseases, age-related cataract and age-related macular degeneration (AMD) are caused by dysfunction of the lens and retina, respectively. Lens transparency and light refraction are mediated by lens fiber cells lacking nuclei and other organelles, which provides a unique opportunity to study a single aging hallmark, i.e., loss of proteostasis, within an environment of limited metabolism. In AMD, local dysfunction of the photoreceptors/retinal pigmented epithelium/Bruch's membrane/choriocapillaris complex in the macula leads to the loss of photoreceptors and eventually loss of central vision, and is driven by nearly all the hallmarks of aging and shares features with Alzheimer's disease, Parkinson's disease, cardiovascular disease, and diabetes. The aging eye can function as a model for studying basic mechanisms of aging and, vice versa, well-defined hallmarks of aging can be used as tools to understand age-related eye disease.

Age- and sex- divergent translatomic responses of the mouse retinal pigmented epithelium

Aging is the main risk factor for age-related macular degeneration (AMD), a retinal neurodegenerative disease that leads to irreversible blindness, particularly in people over 60 years old. Retinal pigmented epithelium (RPE) atrophy is an AMD hallmark. Genome-wide chromatin accessibility, DNA methylation, and gene expression studies of AMD and control RPE demonstrate epigenomic/transcriptomic changes occur during AMD onset and progression. However, mechanisms by which molecular alterations of normal aging impair RPE function and contribute to AMD pathogenesis are unclear. Here, we specifically interrogate the RPE translatome with advanced age and across sexes in a novel RPE reporter mouse model. We find differential age- and sex- associated transcript expression with overrepresentation of pathways related to inflammation in the RPE. Concordant with impaired RPE function, the phenotypic changes in the aged translatome suggest that aged RPE becomes immunologically active, in both males and females, with some sex-specific signatures, which supports the need for sex representation for in vivo studies.

Epigenome-metabolism nexus in the retina: implications for aging and disease

Intimate links between epigenome modifications and metabolites allude to a crucial role of cellular metabolism in transcriptional regulation. Retina, being a highly metabolic tissue, adapts by integrating inputs from genetic, epigenetic, and extracellular signals. Precise global epigenomic signatures guide development and homeostasis of the intricate retinal structure and function. Epigenomic and metabolic realignment are hallmarks of aging and highlight a link of the epigenome-metabolism nexus with aging-associated multifactorial traits affecting the retina, including age-related macular degeneration and glaucoma. Here, we focus on emerging principles of epigenomic and metabolic control of retinal gene regulation, with emphasis on their contribution to human disease. In addition, we discuss potential mitigation strategies involving lifestyle changes that target the epigenome-metabolome relationship for maintaining retinal function.

Creation of a novel zebrafish model with low DHA status to study the role of maternal nutrition during neurodevelopment

Docosahexaenoic acid (DHA), a dietary omega-3 fatty acid, is a major building block of brain cell membranes. Offspring rely on maternal DHA transfer to meet their neurodevelopmental needs, but DHA sources are lacking in the American diet. Low DHA status is linked to altered immune responses, white matter defects, impaired vision, and an increased risk of psychiatric disorders during development. However, the underlying cellular mechanisms involved are largely unknown, and advancements in the field have been limited by the existing tools and animal models. Zebrafish are an excellent model for studying neurodevelopmental mechanisms. Embryos undergo rapid external development and are optically transparent, enabling direct observation of individual cells and dynamic cell-cell interactions in a way that is not possible in rodents. Here, we create a novel DHA-deficient zebrafish model by 1) disrupting elovl2, a key gene in the DHA biosynthesis pathway, via CRISPR-Cas9 genome editing, and 2) feeding mothers a DHA-deficient diet. We show that low DHA status during development is associated with a small eye morphological phenotype and demonstrate that even the morphologically normal siblings exhibit dysregulated gene pathways related to vision and stress response. Future work using our zebrafish model could reveal the cellular and molecular mechanisms by which low DHA status leads to neurodevelopmental abnormalities and provide insight into maternal nutritional strategies that optimize infant brain health.

NMNAT2 is a druggable target to drive neuronal NAD production

Maintenance of NAD pools is critical for neuronal survival. The capacity to maintain NAD pools declines in neurodegenerative disease. We identify that low NMNAT2, the critical neuronal NAD producing enzyme, drives retinal susceptibility to neurodegenerative insults. As proof of concept, gene therapy over-expressing full length human NMNAT2 is neuroprotective. To pharmacologically target NMNAT2, we identify that epigallocatechin gallate (EGCG) can drive NAD production in neurons through an NMNAT2 and NMN dependent mechanism. We confirm this by pharmacological and genetic inhibition of the NAD-salvage pathway. EGCG is neuroprotective in rodent (mixed sex) and human models of retinal neurodegeneration. As EGCG has poor drug-like qualities, we use it as a tool compound to generate novel small molecules which drive neuronal NAD production and provide neuroprotection. This class of NMNAT2 targeted small molecules could have an important therapeutic impact for neurodegenerative disease following further drug development.

Complement Factor B Inhibition or Deletion Is Not Sufficient to Prevent Neurodegeneration in a Murine Model of Glaucoma

Purpose: Activation of the classical complement pathway is thought to contribute to the development and progression of glaucoma. The role of alternative complement or amplification pathways in glaucoma is not well understood. We evaluated complement factor B (FB) expression in postmortem human ocular tissues with or without glaucoma and the effect of FB inhibition and deletion in a mouse ocular hypertensive model of glaucoma induced by photopolymerized hyaluronic acid glycidyl methacrylate (HAGM). Methods: Human CFB mRNA in human eyes was assessed by RNAscope and TaqMan. HAGM model was performed on C57BL6/J mice. The effect of FB in HAGM model was evaluated with an oral FB inhibitor and Cfb-/- mice. Complement mRNA and proteins in mouse eyes were assessed by TaqMan and western blot, respectively. Results: CFB mRNA in human glaucomatous macular neural retina and optic nerve head was upregulated. Cfb mRNA is also upregulated in the HAGM model. Oral FB inhibitor, ED-79-GX17, dosed daily at 200 mg/kg for 3 days after intraocular pressure (IOP) induction in wild-type mice showed complement inhibition in ocular tissues and significantly inhibited systemic complement levels. Daily dosing of ED-79-GX17 for 30 days or Cfb deletion was also unable to prevent retinal ganglion cell or axon loss 30 days after IOP induction in mice. Conclusion: The alternative complement component FB may not substantially contribute to RGC loss in the HAGM mouse glaucoma model despite upregulation of Cfb expression and activation of the alternative pathway. The relevance of these findings to human glaucoma remains to be determined.

Genome-wide association study of the common retinal disorder epiretinal membrane: Significant risk loci in each of three American populations

Epiretinal membrane (ERM) is a common retinal condition characterized by the presence of fibrocellular tissue on the retinal surface, often with visual distortion and loss of visual acuity. We studied European American (EUR), African American (AFR), and Latino (admixed American, AMR) ERM participants in the Million Veteran Program (MVP) for genome-wide association analysis-a total of 38,232 case individuals and 557,988 control individuals. We completed a genome-wide association study (GWAS) in each population separately, and then results were meta-analyzed. Genome-wide significant (GWS) associations were observed in all three populations studied: 31 risk loci in EUR subjects, 3 in AFR, and 2 in AMR, with 48 in trans-ancestry meta-analysis. Many results replicated in the FinnGen sample. Several GWS variants associate to alterations in gene expression in the macula. ERM showed significant genetic correlation to multiple traits. Pathway enrichment analyses implicated collagen and collagen-adjacent mechanisms, among others. This well-powered ERM GWAS identified novel genetic associations that point to biological mechanisms for ERM.

Conditional Knockouts of Interphotoreceptor Retinoid Binding Protein Suggest Two Independent Mechanisms for Retinal Degeneration and Myopia

Purpose

Interphotoreceptor retinoid-binding protein's (IRBP) role in eye growth and its involvement in cell homeostasis remain poorly understood. One hypothesis proposes early conditional deletion of the IRBP gene could lead to a myopic response with retinal degeneration, whereas late conditional deletion (after eye size is determined) could cause retinal degeneration without myopia. Here, we sought to understand if prior myopia was required for subsequent retinal degeneration in the absence of IRBP. This study investigates if any cell type or developmental stage is more important in myopia or retinal degeneration.

Methods

IBRPfl/fl mice were bred with 5 Cre-driver lines: HRGP-Cre, Chx10-Cre, Rho-iCre75, HRGP-Cre Rho-iCre75, and Rx-Cre. Mice were analyzed for IRBP gene expression through digital droplet PCR (ddPCR). Young adult (P30) mice were tested for retinal degeneration and morphology using spectral-domain optical coherence tomography (SD-OCT) and hematoxylin and eosin (H&E) staining. Function was analyzed using electroretinograms (ERGs). Eye sizes and axial lengths were compared through external eye measurements and whole eye biometry.

Results

Across all outcome measures, when bred to IRBPfl/fl, HRGP-Cre and Chx10-Cre lines showed no differences from IRBPfl/fl alone. With the Rho-iCre75 line, small but significant reductions were seen in retinal thickness with SD-OCT imaging and postmortem H&E staining without increased axial length. Both the HRGP-Cre+Rho-iCre75 and the Rx-Cre lines showed significant decreases in retinal thickness and outer nuclear layer cell counts. Using external eye measurements and SD-OCT imaging, both lines showed an increase in eye size. Finally, function in both lines was roughly halved across scotopic, photopic, and flicker ERGs.

Conclusions

Our studies support hypotheses that for both eye size determination and retinal homeostasis, there are two critical timing windows when IRBP must be expressed in rods or cones to prevent myopia (P7-P12) and degeneration (P21 and later). The rod-specific IRBP knockout (Rho-iCre75) showed significant retinal functional losses without myopia, indicating that the two phenotypes are independent. IRBP is needed for early development of photoreceptors and eye size, whereas Rho-iCre75 IRBPfl/fl knockout results in retinal degeneration without myopia.

Vaccines and the Eye: Current Understanding of the Molecular and Immunological Effects of Vaccination on the Eye

Vaccination is a public health cornerstone that protects against numerous infectious diseases. Despite its benefits, immunization implications on ocular health warrant thorough investigation, particularly in the context of vaccine-induced ocular inflammation. This review aimed to elucidate the complex interplay between vaccination and the eye, focusing on the molecular and immunological pathways implicated in vaccine-associated ocular adverse effects. Through an in-depth analysis of recent advancements and the existing literature, we explored various mechanisms of vaccine-induced ocular inflammation, such as direct infection by live attenuated vaccines, immune complex formation, adjuvant-induced autoimmunity, molecular mimicry, hypersensitivity reactions, PEG-induced allergic reactions, Type 1 IFN activation, free extracellular RNA, and specific components. We further examined the specific ocular conditions associated with vaccination, such as uveitis, optic neuritis, and retinitis, and discussed the potential impact of novel vaccines, including those against SARS-CoV-2. This review sheds light on the intricate relationships between vaccination, the immune system, and ocular tissues, offering insights into informed discussions and future research directions aimed at optimizing vaccine safety and ophthalmological care. Our analysis underscores the importance of vigilance and further research to understand and mitigate the ocular side effects of vaccines, thereby ensuring the continued success of vaccination programs, while preserving ocular health.

scDAC: deep adaptive clustering of single-cell transcriptomic data with coupled autoencoder and Dirichlet process mixture model

MOTIVATION

Clustering analysis for single-cell RNA sequencing (scRNA-seq) data is an important step in revealing cellular heterogeneity. Many clustering methods have been proposed to discover heterogenous cell types from scRNA-seq data. However, adaptive clustering with accurate cluster number reflecting intrinsic biology nature from large-scale scRNA-seq data remains quite challenging.

RESULTS

Here, we propose a single-cell Deep Adaptive Clustering (scDAC) model by coupling the Autoencoder (AE) and the Dirichlet Process Mixture Model (DPMM). By jointly optimizing the model parameters of AE and DPMM, scDAC achieves adaptive clustering with accurate cluster numbers on scRNA-seq data. We verify the performance of scDAC on five subsampled datasets with different numbers of cell types and compare it with 15 widely used clustering methods across nine scRNA-seq datasets. Our results demonstrate that scDAC can adaptively find accurate numbers of cell types or subtypes and outperforms other methods. Moreover, the performance of scDAC is robust to hyperparameter changes.

AVAILABILITY AND IMPLEMENTATION

The scDAC is implemented in Python. The source code is available at https://github.com/labomics/scDAC.

Identification of key lncRNAs in age-related macular degeneration through integrated bioinformatics and experimental validation

This study aimed to identify key long noncoding RNAs (lncRNAs) in age-related macular degeneration (AMD) patients and to identify relevant pathological mechanisms of AMD development. We identified 407 differentially expressed mRNAs and 429 differentially expressed lncRNAs in retinal pigment epithelium (RPE) and retina in the macular region of AMD patients versus controls (P < 0.05 and |log2FC| > 0.585) from GSE135092. A total of 14 key differentially expressed mRNAs were obtained through external data validation from GSE115828. A miRNA-mRNA and miRNA-lncRNA network containing 52 lncRNA nodes, 49 miRNA nodes, 14 mRNA nodes and 351 edges was constructed via integrated analysis of these components. Finally, the LINC00276-miR-619-5p-IFIT3 axis was identified via protein-protein network analysis. In the t-BH-induced ARPE-19 senescent cell model, LINC00276 and IFIT3 were downregulated. Overexpression of LINC00276 could accelerate cell migration in combination with IFIT3 upregulation. This compelling finding suggests that LINC00276 plays an influential role in the progression of AMD, potentially through modulating senescence processes, thereby setting a foundation for future investigative efforts to verify this relationship.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation provides a crucial epigenetic mark linking genetic variations to environmental influence. We have analyzed array-based DNA methylation profiles of 160 human retinas with co-measured RNA-seq and >8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 methylation quantitative trait loci and 12,505 expression quantitative trait loci) and 13,747 DNA methylation loci affecting gene expression, with over one-third specific to the retina. Methylation and expression quantitative trait loci show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration. Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of macular degeneration pathology by genotype-environment interaction in retina.

Whole genome sequencing of 4,787 individuals identifies gene-based rare variants in age-related macular degeneration

Genome-wide association studies have contributed extensively to the discovery of disease-associated common variants. However, the genetic contribution to complex traits is still largely difficult to interpret. We report a genome-wide association study of 2394 cases and 2393 controls for age-related macular degeneration (AMD) via whole-genome sequencing, with 46.9 million genetic variants. Our study reveals significant single-variant association signals at four loci and independent gene-based signals in CFH, C2, C3, and NRTN. Using data from the Exome Aggregation Consortium (ExAC) for a gene-based test, we demonstrate an enrichment of predicted rare loss-of-function variants in CFH, CFI, and an as-yet unreported gene in AMD, ORMDL2. Our method of using a large variant list without individual-level genotypes as an external reference provides a flexible and convenient approach to leverage the publicly available variant datasets to augment the search for rare variant associations, which can explain additional disease risk in AMD.

5-Aza-2'-Deoxycytidine Ameliorates Choroidal Neovascularization by Inhibiting the Wnt/β-Catenin Signaling Pathway

Purpose

Choroidal neovascularization (CNV) can constitute the final pathology of many ocular diseases and result in severe vision loss. Studies have demonstrated that DNA methylation is critical in retinal development, aging, and disorders. The current work investigated the effects and underlying mechanism of 5-Aza-2'-deoxycytidine (5-aza-dC), a suppressor of DNA methylation, in the pathological progression of CNV.

Methods

The DNA methylation profiles of retinal pigment epithelial (RPE)/choroidal complexes in normal and laser-induced CNV mice were assessed by Arraystar Mouse RefSeq Promoter Arrays. The CNV area and blood flow density and intensity were observed by optical coherence tomography angiography, and fluorescence leakage was examined by fundus fluorescein angiography in CNV mice with systemic administration of 5-aza-dC. The effects of 5-aza-dC on the biological functions of bEnd.3 cells were estimated by related assays. Notum gene promoter methylation was measured using bisulfite sequencing PCR. Methyltransferases and Wnt signaling-related genes were detected in animal and cell culture experiments by real-time PCR and immunoblot.

Results

Methyltransferases were upregulated, but Notum (a secretion inhibitor of Wnt signaling) was downregulated in the RPE/choroidal complexes of mice with experimental CNV. Intraperitoneal injection of 5-aza-dC inactivated the Wnt pathway and ameliorated the lesion area and the intensity and density of blood flow, as well as the degree of leakage in CNV. In vitro, vascular endothelial growth factor A (VEGFA) stimulation promoted methyltransferases expression and suppressed Notum expression, consequently activating Wnt signaling, whereas exogenous 5-aza-dC reversed VEGFA-induced hyperpermeability, proliferation, migration, and tube formation in bEnd.3 cells via demethylation of Notum promoter.

Conclusions

We observed that 5-aza-dC attenuates the growth of CNV by inhibiting the Wnt signaling pathway via promoter demethylation of the Wnt antagonist Notum. These findings provide a theoretical basis for methylation-based treatment with the Notum gene as a potential target for CNV treatment.

Roles of transmembrane protein 135 in mitochondrial and peroxisomal functions - implications for age-related retinal disease

Aging is the most significant risk factor for age-related diseases in general, which is true for age-related diseases in the eye including age-related macular degeneration (AMD). Therefore, in order to identify potential therapeutic targets for these diseases, it is crucial to understand the normal aging process and how its mis-regulation could cause age-related diseases at the molecular level. Recently, abnormal lipid metabolism has emerged as one major aspect of age-related symptoms in the retina. Animal models provide excellent means to identify and study factors that regulate lipid metabolism in relation to age-related symptoms. Central to this review is the role of transmembrane protein 135 (TMEM135) in the retina. TMEM135 was identified through the characterization of a mutant mouse strain exhibiting accelerated retinal aging and positional cloning of the responsible mutation within the gene, indicating the crucial role of TMEM135 in regulating the normal aging process in the retina. Over the past decade, the molecular functions of TMEM135 have been explored in various models and tissues, providing insights into the regulation of metabolism, particularly lipid metabolism, through its action in multiple organelles. Studies indicated that TMEM135 is a significant regulator of peroxisomes, mitochondria, and their interaction. Here, we provide an overview of the molecular functions of TMEM135 which is crucial for regulating mitochondria, peroxisomes, and lipids. The review also discusses the age-dependent phenotypes in mice with TMEM135 perturbations, emphasizing the importance of a balanced TMEM135 function for the health of the retina and other tissues including the heart, liver, and adipose tissue. Finally, we explore the potential roles of TMEM135 in human age-related retinal diseases, connecting its functions to the pathobiology of AMD.

Catalytic isoforms of AMP-activated protein kinase differentially regulate IMPDH activity and photoreceptor neuron function

AMP-activated protein kinase (AMPK) plays a crucial role in maintaining ATP homeostasis in photoreceptor neurons. AMPK is a heterotrimeric protein consisting of α, β, and γ subunits. The independent functions of the 2 isoforms of the catalytic α subunit, PRKAA1 and PRKAA2, are uncharacterized in specialized neurons, such as photoreceptors. Here, we demonstrate in mice that rod photoreceptors lacking PRKAA2, but not PRKAA1, showed altered levels of cGMP, GTP, and ATP, suggesting isoform-specific regulation of photoreceptor metabolism. Furthermore, PRKAA2-deficient mice displayed visual functional deficits on electroretinography and photoreceptor outer segment structural abnormalities on transmission electron microscopy consistent with neuronal dysfunction, but not neurodegeneration. Phosphoproteomics identified inosine monophosphate dehydrogenase (IMPDH) as a molecular driver of PRKAA2-specific photoreceptor dysfunction, and inhibition of IMPDH improved visual function in Prkaa2 rod photoreceptor-knockout mice. These findings highlight a therapeutically targetable PRKAA2 isoform-specific function of AMPK in regulating photoreceptor metabolism and function through a potentially previously uncharacterized mechanism affecting IMPDH activity.

Integrating genetic regulation and single-cell expression with GWAS prioritizes causal genes and cell types for glaucoma

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.

Systemic Lipopolysaccharide Exposure Exacerbates Choroidal Neovascularization in Mice

This study aims to investigate the effect of a systemic lipopolysaccharide (LPS) stimulus in the course of laser-induced choroidal neovascularization (CNV) in C57BL/6 J mice. A group of CNV-subjected mice received 1 mg/kg LPS via the tail vein immediately after CNV induction. Mouse eyes were monitored in vivo with fluorescein angiography for 2 weeks. In situ hybridization and flow cytometry were performed in the retina at different time points. LPS led to increased fluorescein leakage 3 days after CNV, correlated with a large influx of monocyte-derived macrophages and increase of pro-inflammatory microglia/macrophages in the retina. Additionally, LPS enhanced Vegfα mRNA expression by Glul-expressing cells but not Aif1 positive microglia/macrophages in the laser lesion. These findings suggest that systemic LPS exposure has transient detrimental effects in the course of CNV through activation of microglia/macrophages to a pro-inflammatory phenotype and supports the important role of these cells in the CNV course.

Geographic atrophy: pathophysiology and current therapeutic strategies

Geographic atrophy (GA) is an advanced stage of age-related macular degeneration (AMD) that leads to gradual and permanent vision loss. GA is characterized by the loss of photoreceptor cells and retinal pigment epithelium (RPE), leading to distinct atrophic patches in the macula, which tends to increase with time. Patients with geographic atrophy often experience a gradual and painless loss of central vision, resulting in difficulty reading, recognizing faces, or performing activities that require detailed vision. The primary risk factor for the development of geographic atrophy is advanced age; however, other risk factors, such as family history, smoking, and certain genetic variations, are also associated with AMD. Diagnosis is usually based on a comprehensive eye examination, including imaging tests such as fundus photography, optical coherence tomography (OCT), and fluorescein angiography. Numerous clinical trials are underway, targeting identified molecular pathways associated with GA that are promising. Recent approvals of Syfovre and Izervay by the FDA for the treatment of GA provide hope to affected patients. Administration of these drugs resulted in slowing the rate of progression of the disease. Though these products provide treatment benefits to the patients, they do not offer a cure for geographic atrophy and are limited in efficacy. Considering these safety concerns and limited treatment benefits, there is still a significant need for therapeutics with improved efficacy, safety profiles, and better patient compliance. This comprehensive review discusses pathophysiology, currently approved products, their limitations, and potential future treatment strategies for GA.

Patterns of Gene Expression, Splicing, and Allele-Specific Expression Vary among Macular Tissues and Clinical Stages of Age-Related Macular Degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness, and elucidating its underlying disease mechanisms is vital to the development of appropriate therapeutics. We identified differentially expressed genes (DEGs) and differentially spliced genes (DSGs) across the clinical stages of AMD in disease-affected tissue, the macular retina pigment epithelium (RPE)/choroid and the macular neural retina within the same eye. We utilized 27 deeply phenotyped donor eyes (recovered within a 6 h postmortem interval time) from Caucasian donors (60-94 years) using a standardized published protocol. Significant findings were then validated in an independent set of well-characterized donor eyes (n = 85). There was limited overlap between DEGs and DSGs, suggesting distinct mechanisms at play in AMD pathophysiology. A greater number of previously reported AMD loci overlapped with DSGs compared to DEGs between disease states, and no DEG overlap with previously reported loci was found in the macular retina between disease states. Additionally, we explored allele-specific expression (ASE) in coding regions of previously reported AMD risk loci, uncovering a significant imbalance in C3 rs2230199 and CFH rs1061170 in the macular RPE/choroid for normal eyes and intermediate AMD (iAMD), and for CFH rs1061147 in the macular RPE/choroid for normal eyes and iAMD, and separately neovascular AMD (NEO). Only significant DEGs/DSGs from the macular RPE/choroid were found to overlap between disease states. STAT1, validated between the iAMD vs. normal comparison, and AGTPBP1, BBS5, CERKL, FGFBP2, KIFC3, RORα, and ZNF292, validated between the NEO vs. normal comparison, revealed an intricate regulatory network with transcription factors and miRNAs identifying potential upstream and downstream regulators. Findings regarding the complement genes C3 and CFH suggest that coding variants at these loci may influence AMD development via an imbalance of gene expression in a tissue-specific manner. Our study provides crucial insights into the multifaceted genomic underpinnings of AMD (i.e., tissue-specific gene expression changes, potential splice variation, and allelic imbalance), which may open new avenues for AMD diagnostics and therapies specific to iAMD and NEO.

PolySialic acid-nanoparticles inhibit macrophage mediated inflammation through Siglec agonism: a potential treatment for age related macular degeneration

Age-related macular degeneration (AMD) is a chronic, progressive retinal disease characterized by an inflammatory response mediated by activated macrophages and microglia infiltrating the inner layer of the retina. In this study, we demonstrate that inhibition of macrophages through Siglec binding in the AMD eye can generate therapeutically useful effects. We show that Siglecs-7, -9 and -11 are upregulated in AMD associated M0 and M1 macrophages, and that these can be selectively targeted using polysialic acid (PolySia)-nanoparticles (NPs) to control dampen AMD-associated inflammation. In vitro studies showed that PolySia-NPs bind to macrophages through human Siglecs-7, -9, -11 as well as murine ortholog Siglec-E. Following treatment with PolySia-NPs, we observed that the PolySia-NPs bound and agonized the macrophage Siglecs resulting in a significant decrease in the secretion of IL-6, IL-1β, TNF-α and VEGF, and an increased secretion of IL-10. In vivo intravitreal (IVT) injection of PolySia-NPs was found to be well-tolerated and safe making it effective in preventing thinning of the retinal outer nuclear layer (ONL), inhibiting macrophage infiltration, and restoring electrophysiological retinal function in a model of bright light-induced retinal degeneration. In a clinically validated, laser-induced choroidal neovascularization (CNV) model of exudative AMD, PolySia-NPs reduced the size of neovascular lesions with associated reduction in macrophages. The PolySia-NPs described herein are therefore a promising therapeutic strategy for repolarizing pro-inflammatory macrophages to a more anti-inflammatory, non-angiogenic phenotype, which play a key role in the pathophysiology of non-exudative AMD.

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.

Exploring the contribution of ARMS2 and HTRA1 genetic risk factors in age-related macular degeneration

Age-related macular degeneration (AMD) is the leading cause of severe irreversible central vision loss in individuals over 65 years old. Genome-wide association studies (GWASs) have shown that the region at chromosome 10q26, where the age-related maculopathy susceptibility (ARMS2/LOC387715) and HtrA serine peptidase 1 (HTRA1) genes are located, represents one of the strongest associated loci for AMD. However, the underlying biological mechanism of this genetic association has remained elusive. In this article, we extensively review the literature by us and others regarding the ARMS2/HTRA1 risk alleles and their functional significance. We also review the literature regarding the presumed function of the ARMS2 protein and the molecular processes of the HTRA1 protein in AMD pathogenesis in vitro and in vivo, including those of transgenic mice overexpressing HtrA1/HTRA1 which developed Bruch's membrane (BM) damage, choroidal neovascularization (CNV), and polypoidal choroidal vasculopathy (PCV), similar to human AMD patients. The elucidation of the molecular mechanisms of the ARMS2 and HTRA1 susceptibility loci has begun to untangle the complex biological pathways underlying AMD pathophysiology, pointing to new testable paradigms for treatment.

Limitations of direct-to-consumer genetic testing for age-related macular degeneration

The availability of direct-to-consumer genetic testing (DTCGT) for age-related macular degeneration (AMD) provides the public with access to disease risk estimations that may be used to guide lifestyle adjustments. However, AMD development risk is more complex than can be captured by gene mutations alone. The methodologies employed by current DTCGTs to estimate AMD risk vary and are limited in several ways. Genotyping-based DTCGT is biased toward European ancestry and only considers a limited number of genes. Whole genome sequencing based DTCGTs uncovers several genetic variations with unknown relevance, making risk interpretation challenging. In this perspective, we describe the limitations of the DTCGT for AMD.

Human retinal ganglion cell neurons generated by synchronous BMP inhibition and transcription factor mediated reprogramming

In optic neuropathies, including glaucoma, retinal ganglion cells (RGCs) die. Cell transplantation and endogenous regeneration offer strategies for retinal repair, however, developmental programs required for this to succeed are incompletely understood. To address this, we explored cellular reprogramming with transcription factor (TF) regulators of RGC development which were integrated into human pluripotent stem cells (PSCs) as inducible gene cassettes. When the pioneer factor NEUROG2 was combined with RGC-expressed TFs (ATOH7, ISL1, and POU4F2) some conversion was observed and when pre-patterned by BMP inhibition, RGC-like induced neurons (RGC-iNs) were generated with high efficiency in just under a week. These exhibited transcriptional profiles that were reminiscent of RGCs and exhibited electrophysiological properties, including AMPA-mediated synaptic transmission. Additionally, we demonstrated that small molecule inhibitors of DLK/LZK and GCK-IV can block neuronal death in two pharmacological axon injury models. Combining developmental patterning with RGC-specific TFs thus provided valuable insight into strategies for cell replacement and neuroprotection.

CellAnn: a comprehensive, super-fast, and user-friendly single-cell annotation web server

MOTIVATION

Single-cell sequencing technology has become a routine in studying many biological problems. A core step of analyzing single-cell data is the assignment of cell clusters to specific cell types. Reference-based methods are proposed for predicting cell types for single-cell clusters. However, the scalability and lack of preprocessed reference datasets prevent them from being practical and easy to use.

RESULTS

Here, we introduce a reference-based cell annotation web server, CellAnn, which is super-fast and easy to use. CellAnn contains a comprehensive reference database with 204 human and 191 mouse single-cell datasets. These reference datasets cover 32 organs. Furthermore, we developed a cluster-to-cluster alignment method to transfer cell labels from the reference to the query datasets, which is superior to the existing methods with higher accuracy and higher scalability. Finally, CellAnn is an online tool that integrates all the procedures in cell annotation, including reference searching, transferring cell labels, visualizing results, and harmonizing cell annotation labels. Through the user-friendly interface, users can identify the best annotation by cross-validating with multiple reference datasets. We believe that CellAnn can greatly facilitate single-cell sequencing data analysis.

AVAILABILITY AND IMPLEMENTATION

The web server is available at www.cellann.io, and the source code is available at https://github.com/Pinlyu3/CellAnn_shinyapp.

TNFα induced by DNA-sensing in macrophage compromises retinal pigment epithelial (RPE) barrier function

Increasing evidence suggests that chronic inflammation plays an important role in the pathogenesis of age-related macular degeneration (AMD); however, the precise pathogenic stressors and sensors, and their impact on disease progression remain unclear. Several studies have demonstrated that type I interferon (IFN) response is activated in the retinal pigment epithelium (RPE) of AMD patients. Previously, we demonstrated that human RPE cells can initiate RNA-mediated type I IFN responses through RIG-I, yet are unable to directly sense and respond to DNA. In this study, we utilized a co-culture system combining primary human macrophage and iPS-derived RPE to study how each cell type responds to nucleic acids challenges and their effect on RPE barrier function in a homotypic and heterotypic manner. We find that DNA-induced macrophage activation induces an IFN response in the RPE, and compromises RPE barrier function via tight-junction remodeling. Investigation of the secreted cytokines responsible for RPE dysfunction following DNA-induced macrophages activation indicates that neutralization of macrophage-secreted TNFα, but not IFNβ, is sufficient to rescue RPE morphology and barrier function. Our data reveals a novel mechanism of intercellular communication by which DNA induces RPE dysfunction via macrophage-secreted TNFa, highlighting the complexity and potential pathological relevance of RPE and macrophage interactions.

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.

Complement factor I: Regulatory nexus, driver of immunopathology, and therapeutic

Complement factor I (FI) is the nexus for classical, lectin and alternative pathway complement regulation. FI is an 88 kDa plasma protein that circulates in an inactive configuration until it forms a trimolecular complex with its cofactor and substrate whereupon a structural reorganization allows the catalytic triad to cleave its substrates, C3b and C4b. In keeping with its role as the master complement regulatory enzyme, deficiency has been linked to immunopathology. In the setting of complete FI deficiency, a consumptive C3 deficiency results in recurrent infections with encapsulated microorganisms. Aseptic cerebral inflammation and vasculitic presentations are also less commonly observed. Heterozygous mutations in the factor I gene (CFI) have been demonstrated to be enriched in atypical haemolytic uraemic syndrome, albeit with a very low penetrance. Haploinsufficiency of CFI has also been associated with decreased retinal thickness and is a strong risk factor for the development of age-related macular degeneration. Supplementation of FI using plasma purified or recombinant protein has long been postulated, however, technical difficulties prevented progression into clinical trials. It is only using gene therapy that CFI supplementation has reached the clinic with GT005 in phase I/II clinical trials for geographic atrophy.

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.

Integration of human stem cell-derived in vitro systems and mouse preclinical models identifies complex pathophysiologic mechanisms in retinal dystrophy

Rare DRAM2 coding variants cause retinal dystrophy with early macular involvement via unknown mechanisms. We found that DRAM2 is ubiquitously expressed in the human eye and expression changes were observed in eyes with more common maculopathy such as Age-related Macular Degeneration (AMD). To gain insights into pathogenicity of DRAM2-related retinopathy, we used a combination of in vitro and in vivo models. We found that DRAM2 loss in human pluripotent stem cell (hPSC)-derived retinal organoids caused the presence of additional mesenchymal cells. Interestingly, Dram2 loss in mice also caused increased proliferation of cells from the choroid in vitro and exacerbated choroidal neovascular lesions in vivo. Furthermore, we observed that DRAM2 loss in human retinal pigment epithelial (RPE) cells resulted in increased susceptibility to stress-induced cell death in vitro and that Dram2 loss in mice caused age-related photoreceptor degeneration. This highlights the complexity of DRAM2 function, as its loss in choroidal cells provided a proliferative advantage, whereas its loss in post-mitotic cells, such as photoreceptor and RPE cells, increased degeneration susceptibility. Different models such as human pluripotent stem cell-derived systems and mice can be leveraged to study and model human retinal dystrophies; however, cell type and species-specific expression must be taken into account when selecting relevant systems.

Multi-ancestry GWAS analysis identifies two novel loci associated with diabetic eye disease and highlights APOL1 as a high risk locus in patients with diabetic macular edema

Diabetic retinopathy (DR) is a common complication of diabetes. Approximately 20% of DR patients have diabetic macular edema (DME) characterized by fluid leakage into the retina. There is a genetic component to DR and DME risk, but few replicable loci. Because not all DR cases have DME, we focused on DME to increase power, and conducted a multi-ancestry GWAS to assess DME risk in a total of 1,502 DME patients and 5,603 non-DME controls in discovery and replication datasets. Two loci reached GWAS significance (p<5x10-8). The strongest association was rs2239785, (K150E) in APOL1. The second finding was rs10402468, which co-localized to PLVAP and ANKLE1 in vascular / endothelium tissues. We conducted multiple sensitivity analyses to establish that the associations were specific to DME status and did not reflect diabetes status or other diabetic complications. Here we report two novel loci for risk of DME which replicated in multiple clinical trial and biobank derived datasets. One of these loci, containing the gene APOL1, is a risk factor in African American DME and DKD patients, indicating that this locus plays a broader role in diabetic complications for multiple ancestries. Trial Registration: NCT00473330, NCT00473382, NCT03622580, NCT03622593, NCT04108156.

Cross-species scRNA-seq reveals the cellular landscape of retina and early alterations in type 2 diabetes mice

Single-cell RNA sequencing (scRNA-seq) analysis have provided an unprecedented resolution for the studies on diabetic retinopathy (DR). However, the early changes in the retina in diabetes remain unclear. A total of 8 human and mouse scRNA-seq datasets, containing 276,402 cells were analyzed individually to comprehensively delineate the retinal cell atlas. The neural retinas were isolated from the type 2 diabetes (T2D) and control mice, and scRNA-seq analysis was conducted to evaluate the early effects of diabetes on the retina. Bipolar cell (BC) heterogeneity were identified. We found some stable BCs across multiple datasets, and explored their biological functions. A new RBC subtype (Car8_RBC) in the mouse retina was validated using the multi-color immunohistochemistry. AC149090.1 was significantly upregulated in the rod cells, ON cone BCs (CBCs), OFF CBCs, and RBCs in T2D mice. Additionally, the interneurons, especially BCs, were the most vulnerable cells to diabetes by integrating scRNA-seq and genome-wide association studies (GWAS) analyses. In conclusion, this study delineated a cross-species retinal cell atlas and uncovered the early pathological alterations in the retina of T2D mice.

QTL mapping of human retina DNA methylation identifies 87 gene-epigenome interactions in age-related macular degeneration

DNA methylation (DNAm) provides a crucial epigenetic mark linking genetic variations to environmental influence. We analyzed array-based DNAm profiles of 160 human retinas with co-measured RNA-seq and > 8 million genetic variants, uncovering sites of genetic regulation in cis (37,453 mQTLs and 12,505 eQTLs) and 13,747 eQTMs (DNAm loci affecting gene expression), with over one-third specific to the retina. mQTLs and eQTMs show non-random distribution and enrichment of biological processes related to synapse, mitochondria, and catabolism. Summary data-based Mendelian randomization and colocalization analyses identify 87 target genes where methylation and gene-expression changes likely mediate the genotype effect on age-related macular degeneration (AMD). Integrated pathway analysis reveals epigenetic regulation of immune response and metabolism including the glutathione pathway and glycolysis. Our study thus defines key roles of genetic variations driving methylation changes, prioritizes epigenetic control of gene expression, and suggests frameworks for regulation of AMD pathology by genotype-environment interaction in retina.

A multi-omics atlas of the human retina at single-cell resolution

Cell classes in the human retina are highly heterogeneous with their abundance varying by several orders of magnitude. Here, we generated and integrated a multi-omics single-cell atlas of the adult human retina, including more than 250,000 nuclei for single-nuclei RNA-seq and 137,000 nuclei for single-nuclei ATAC-seq. Cross-species comparison of the retina atlas among human, monkey, mice, and chicken revealed relatively conserved and non-conserved types. Interestingly, the overall cell heterogeneity in primate retina decreases compared with that of rodent and chicken retina. Through integrative analysis, we identified 35,000 distal cis-element-gene pairs, constructed transcription factor (TF)-target regulons for more than 200 TFs, and partitioned the TFs into distinct co-active modules. We also revealed the heterogeneity of the cis-element-gene relationships in different cell types, even from the same class. Taken together, we present a comprehensive single-cell multi-omics atlas of the human retina as a resource that enables systematic molecular characterization at individual cell-type resolution.

A systems biology approach uncovers novel disease mechanisms in age-related macular degeneration

Age-related macular degeneration (AMD) is a leading cause of blindness, affecting 200 million people worldwide. To identify genes that could be targeted for treatment, we created a molecular atlas at different stages of AMD. Our resource is comprised of RNA sequencing (RNA-seq) and DNA methylation microarrays from bulk macular retinal pigment epithelium (RPE)/choroid of clinically phenotyped normal and AMD donor eyes (n = 85), single-nucleus RNA-seq (164,399 cells), and single-nucleus assay for transposase-accessible chromatin (ATAC)-seq (125,822 cells) from the retina, RPE, and choroid of 6 AMD and 7 control donors. We identified 23 genome-wide significant loci differentially methylated in AMD, over 1,000 differentially expressed genes across different disease stages, and an AMD Müller state distinct from normal or gliosis. Chromatin accessibility peaks in genome-wide association study (GWAS) loci revealed putative causal genes for AMD, including HTRA1 and C6orf223. Our systems biology approach uncovered molecular mechanisms underlying AMD, including regulators of WNT signaling, FRZB and TLE2, as mechanistic players in disease.

Reliable Identification and Interpretation of Single-Cell Molecular Heterogeneity and Transcriptional Regulation using Dynamic Ensemble Pruning

Unsupervised clustering is an essential step in identifying cell types from single-cell RNA sequencing (scRNA-seq) data. However, a common issue with unsupervised clustering models is that the optimization direction of the objective function and the final generated clustering labels in the absence of supervised information may be inconsistent or even arbitrary. To address this challenge, a dynamic ensemble pruning framework (DEPF) is proposed to identify and interpret single-cell molecular heterogeneity. In particular, a silhouette coefficient-based indicator is developed to determine the optimization direction of the bi-objective function. In addition, a hierarchical autoencoder is employed to project the high-dimensional data onto multiple low-dimensional latent space sets, and then a clustering ensemble is produced in the latent space by the basic clustering algorithm. Following that, a bi-objective fruit fly optimization algorithm is designed to prune dynamically the low-quality basic clustering in the ensemble. Multiple experiments are conducted on 28 real scRNA-seq datasets and one large real scRNA-seq dataset from diverse platforms and species to validate the effectiveness of the DEPF. In addition, biological interpretability and transcriptional and post-transcriptional regulatory are conducted to explore biological patterns from the cell types identified, which could provide novel insights into characterizing the mechanisms.

Transfer learning enables predictions in network biology

Mapping gene networks requires large amounts of transcriptomic data to learn the connections between genes, which impedes discoveries in settings with limited data, including rare diseases and diseases affecting clinically inaccessible tissues. Recently, transfer learning has revolutionized fields such as natural language understanding1,2 and computer vision3 by leveraging deep learning models pretrained on large-scale general datasets that can then be fine-tuned towards a vast array of downstream tasks with limited task-specific data. Here, we developed a context-aware, attention-based deep learning model, Geneformer, pretrained on a large-scale corpus of about 30 million single-cell transcriptomes to enable context-specific predictions in settings with limited data in network biology. During pretraining, Geneformer gained a fundamental understanding of network dynamics, encoding network hierarchy in the attention weights of the model in a completely self-supervised manner. Fine-tuning towards a diverse panel of downstream tasks relevant to chromatin and network dynamics using limited task-specific data demonstrated that Geneformer consistently boosted predictive accuracy. Applied to disease modelling with limited patient data, Geneformer identified candidate therapeutic targets for cardiomyopathy. Overall, Geneformer represents a pretrained deep learning model from which fine-tuning towards a broad range of downstream applications can be pursued to accelerate discovery of key network regulators and candidate therapeutic targets.

Lipid nanoparticles for gene therapy in ocular diseases

OBJECTIVES

Lipid nanoparticles, as a nucleic acid delivery system, have been used as an alternative to treat ocular diseases, since they can cross the ocular barrier and efficiently transfecting nucleic acids to various cells of the eye. The size influences the transfection of genes, biological distribution, diffusion, and cellular uptake. It is therefore important to establish a relationship between size, formulation, and encapsulation percentage.

EVIDENCE ACQUISITION

In this review, we used a search strategy to compare studies of nanomedicine systems aimed at eye diseases where the size of the nanoparticles and the efficiency of encapsulation of genetic material are reported based on the criteria of Preferred Reporting Items for Systematic Reviews (PRISMA ScR 2020 guidelines).

RESULTS

Out of the initial 5932, 169 studies met the inclusion criteria and were included to form the basis of the analysis. Nanoparticles reported are composed mainly of PEG-modified lipids, cholesterol, and cationic lipids, that in combination with messenger or interference RNA, allow the formulation of a nanoparticle with an encapsulation efficiency greater than 95%. The diseases treated mainly focus on conditions related to the retina and cornea. Certain characteristics of nanoparticles increase encapsulation efficiency, such as the size of the nanoparticle and the charge of the outer layer of the nanoparticle.

CONCLUSION

It is still unknown what characteristics lipid nanoparticles should have to successfully treat human eye illnesses. The in vitro and in vivo investigations covered in this review, however, present encouraging results. To improve encapsulation effectiveness and disease gene silencing, nanoparticle formulation is essential. The most stable nanoparticles are those made mostly of cationic lipids, PEG lipids, and cholesterol, which also effectively encapsulate RNA. The encapsulation efficiency is not only influenced by size, but also by other factors such as methods of preparation.

Circular RNA expression and the competitive endogenous RNA network in pathological, age-related macular degeneration events: A cross-platform normalization study

Age-related macular degeneration (AMD) causes irreversible blindness in people aged over 50 worldwide. The dysfunction of the retinal pigment epithelium is the primary cause of atrophic AMD. In the current study, we used the ComBat and Training Distribution Matching method to integrate data obtained from the Gene Expression Omnibus database. We analyzed the integrated sequencing data by the Gene Set Enrichment Analysis. Peroxisome and tumor necrosis factor-α (TNF-α) signaling and nuclear factor kappa B (NF-κB) were among the top 10 pathways, and thus we selected them to construct AMD cell models to identify differentially expressed circular RNAs (circRNAs). We then constructed a competing endogenous RNA network, which is related to differentially expressed circRNAs. This network included seven circRNAs, 15 microRNAs, and 82 mRNAs. The Kyoto Encyclopedia of Genes and Genomes analysis of mRNAs in this network showed that the hypoxia-inducible factor-1 (HIF-1) signaling pathway was a common downstream event. The results of the current study may provide insights into the pathological processes of atrophic AMD.

Transcriptomic Analysis of the Ocular Posterior Segment Completes a Cell Atlas of the Human Eye

Although the visual system extends through the brain, most vision loss originates from defects in the eye. Its central element is the neural retina, which senses light, processes visual signals, and transmits them to the rest of the brain through the optic nerve (ON). Surrounding the retina are numerous other structures, conventionally divided into anterior and posterior segments. Here we used high-throughput single nucleus RNA sequencing (snRNA-seq) to classify and characterize cells in the extraretinal components of the posterior segment: ON, optic nerve head (ONH), peripheral sclera, peripapillary sclera (PPS), choroid, and retinal pigment epithelium (RPE). Defects in each of these tissues are associated with blinding diseases - for example, glaucoma (ONH and PPS), optic neuritis (ON), retinitis pigmentosa (RPE), and age-related macular degeneration (RPE and choroid). From ∼151,000 single nuclei, we identified 37 transcriptomically distinct cell types, including multiple types of astrocytes, oligodendrocytes, fibroblasts, and vascular endothelial cells. Our analyses revealed a differential distribution of many cell types among distinct structures. Together with our previous analyses of the anterior segment and retina, the new data complete a "Version 1" cell atlas of the human eye. We used this atlas to map the expression of >180 genes associated with the risk of developing glaucoma, which is known to involve ocular tissues in both anterior and posterior segments as well as neural retina. Similar methods can be used to investigate numerous additional ocular diseases, many of which are currently untreatable.

Transcriptomic Changes Predict Metabolic Alterations in LC3 Associated Phagocytosis in Aged Mice

LC3b (Map1lc3b) plays an essential role in canonical autophagy and is one of several components of the autophagy machinery that mediates non-canonical autophagic functions. Phagosomes are often associated with lipidated LC3b to promote phagosome maturation in a process called LC3-associated phagocytosis (LAP). Specialized phagocytes, such as mammary epithelial cells, retinal pigment epithelial (RPE) cells, and sertoli cells, utilize LAP for optimal degradation of phagocytosed material, including debris. In the visual system, LAP is critical to maintain retinal function, lipid homeostasis, and neuroprotection. In a mouse model of retinal lipid steatosis-mice lacking LC3b (LC3b-/-), we observed increased lipid deposition, metabolic dysregulation, and enhanced inflammation. Herein, we present a non-biased approach to determine if loss of LAP mediated processes modulate the expression of various genes related to metabolic homeostasis, lipid handling, and inflammation. A comparison of the RPE transcriptome of WT and LC3b-/- mice revealed 1533 DEGs, with ~73% upregulated and 27% downregulated. Enriched gene ontology (GO) terms included inflammatory response (upregulated DEGs), fatty acid metabolism, and vascular transport (downregulated DEGs). Gene set enrichment analysis (GSEA) identified 34 pathways; 28 were upregulated (dominated by inflammation/related pathways) and 6 were downregulated (dominated by metabolic pathways). Analysis of additional gene families identified significant differences for genes in the solute carrier family, RPE signature genes, and genes with a potential role in age-related macular degeneration. These data indicate that loss of LC3b induces robust changes in the RPE transcriptome contributing to lipid dysregulation and metabolic imbalance, RPE atrophy, inflammation, and disease pathophysiology.

Genome-Wide Association Study of Age-Related Macular Degeneration Reveals 2 New Loci Implying Shared Genetic Components with Central Serous Chorioretinopathy

PURPOSE

To investigate the genetic architecture of age-related macular degeneration (AMD) in a Japanese population.

DESIGN

Genome-wide association study (GWAS).

PARTICIPANTS

Three thousand seven hundred seventy-two patients with AMD and 16 770 control participants from the Japanese population were enrolled in the association analyses.

METHODS

We conducted a meta-analysis of 2 independent GWASs that included a total of 2663 patients with AMD and 9471 control participants using the imputation reference panel for genotype imputation specified for the Japanese population (n = 3541). A replication study was performed using an independent set of 1109 patients with AMD and 7299 control participants.

MAIN OUTCOME MEASURES

Associations of genetic variants with AMD.

RESULTS

A meta-analysis of the 2 GWASs identified 6 loci significantly associated with AMD (P < 5.0 × 10^-8). Of these loci, 4 were known to be associated with AMD (CFH, C2/FB, TNFRSF10A, and ARMS2), and 2 were novel (rs4147157 near WBP1L and rs76228488 near GATA5). The newly identified associations were confirmed in a replication study (P < 0.01). After the meta-analysis of all datasets, we observed strong associations in these loci (P = 1.88 × 10^-12 and P = 1.35 × 10^-9 for meta-analysis for rs4147157 and rs76228488, respectively). When we looked up the associations in the reported central serous chorioretinopathy (CSC) GWAS conducted in the Japanese population, both loci were associated significantly with CSC (P = 4.86 × 10^-3 and P = 4.28 × 10^-3 for rs4147157 and rs76228488, respectively). We performed a genetic colocalization analysis for these loci and estimated that the posterior probabilities of shared causal variants between AMD and CSC were 0.39 and 0.60 for WBP1L and GATA5, respectively. Genetic correlation analysis focusing on the epidemiologically suggested clinical risk factors implicated shared polygenic architecture between AMD and smoking cessation (r_g [the measure of genetic correlation] = -0.33; P = 0.01; false discovery rate, 0.099).

CONCLUSIONS

Our findings imply shared genetic components conferring the risk of both AMD and CSC.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found after the references.

Transcriptomic changes predict metabolic alterations in LC3 associated phagocytosis in aged mice

LC3b ( Map1lc3b ) plays an essential role in canonical autophagy and is one of several components of the autophagy machinery that mediates non-canonical autophagic functions. Phagosomes are often associated with lipidated LC3b, to pro-mote phagosome maturation in a process called LC3-associated phagocytosis (LAP). Specialized phagocytes such as mammary epithelial cells, retinal pigment epithelial (RPE) cells, and sertoli cells utilize LAP for optimal degradation of phagocytosed material, including debris. In the visual system, LAP is critical to maintain retinal function, lipid homeostasis and neuroprotection. In a mouse model of retinal lipid steatosis - mice lacking LC3b ( LC3b ^-/- ), we observed increased lipid deposition, metabolic dysregulation and enhanced inflammation. Herein we present a non-biased approach to determine if loss of LAP mediated processes modulate the expression of various genes related to metabolic homeostasis, lipid handling, and inflammation. A comparison of the RPE transcriptome of WT and LC3b ^-/- mice revealed 1533 DEGs, with ~73% upregulated and 27% down-regulated. Enriched gene ontology (GO) terms included inflammatory response (upregulated DEGs), fatty acid metabolism and vascular transport (downregulated DEGs). Gene set enrichment analysis (GSEA) identified 34 pathways; 28 were upregulated (dominated by inflammation/related pathways) and 6 were downregulated (dominated by metabolic pathways). Analysis of additional gene families identified significant differences for genes in the solute carrier family, RPE signature genes, and genes with potential role in age-related macular degeneration. These data indicate that loss of LC3b induces robust changes in the RPE transcriptome contributing to lipid dysregulation and metabolic imbalance, RPE atrophy, inflammation, and disease pathophysiology.

Genome-wide analysis of genetic pleiotropy and causal genes across three age-related ocular disorders

Age-related macular degeneration (AMD), cataract, and glaucoma are leading causes of blindness worldwide. Previous genome-wide association studies (GWASs) have revealed a variety of susceptible loci associated with age-related ocular disorders, yet the genetic pleiotropy and causal genes across these diseases remain poorly understood. By leveraging large-scale genetic and observational data from ocular disease GWASs and UK Biobank (UKBB), we found significant pairwise genetic correlations and consistent epidemiological associations among these ocular disorders. Cross-disease meta-analysis uncovered seven pleiotropic loci, three of which were replicated in an additional cohort. Integration of variants in pleiotropic loci and multiple single-cell omics data identified that Müller cells and astrocytes were likely trait-related cell types underlying ocular comorbidity. In addition, we comprehensively integrated eye-specific gene expression quantitative loci (eQTLs), epigenomic profiling, and 3D genome data to prioritize causal pleiotropic genes. We found that pleiotropic genes were essential in nerve development and eye pigmentation, and targetable by aflibercept and pilocarpine for the treatment of AMD and glaucoma. These findings will not only facilitate the mechanistic research of ocular comorbidities but also benefit the therapeutic optimization of age-related ocular diseases.

Ambra1 haploinsufficiency in CD1 mice results in metabolic alterations and exacerbates age-associated retinal degeneration

Macroautophagy/autophagy is a key process in the maintenance of cellular homeostasis. The age-dependent decline in retinal autophagy has been associated with photoreceptor degeneration. Retinal dysfunction can also result from damage to the retinal pigment epithelium (RPE), as the RPE-retina constitutes an important metabolic ecosystem that must be finely tuned to preserve visual function. While studies of mice lacking essential autophagy genes have revealed a predisposition to retinal degeneration, the consequences of a moderate reduction in autophagy, similar to that which occurs during physiological aging, remain unclear. Here, we described a retinal phenotype consistent with accelerated aging in mice carrying a haploinsufficiency for Ambra1, a pro-autophagic gene. These mice showed protein aggregation in the retina and RPE, metabolic underperformance, and premature vision loss. Moreover, Ambra1^+/gt mice were more prone to retinal degeneration after RPE stress. These findings indicate that autophagy provides crucial support to RPE-retinal metabolism and protects the retina against stress and physiological aging.Abbreviations : 4-HNE: 4-hydroxynonenal; AMBRA1: autophagy and beclin 1 regulator 1, AMD: age-related macular degeneration;; GCL: ganglion cell layer; GFAP: glial fibrillary acidic protein; GLUL: glutamine synthetase/glutamate-ammonia ligase; HCL: hierarchical clustering; INL: inner nuclear layer; IPL: inner plexiform layer; LC/GC-MS: liquid chromatography/gas chromatography-mass spectrometry; MA: middle-aged; MTDR: MitoTracker Deep Red; MFI: mean fluorescence intensity; NL: NH₄Cl and leupeptin; Nqo: NAD(P)H quinone dehydrogenase; ONL: outer nuclear layer; OPL: outer plexiform layer; OP: oscillatory potentials; OXPHOS: oxidative phosphorylation; PCR: polymerase chain reaction; PRKC/PKCα: protein kinase C; POS: photoreceptor outer segment; RGC: retinal ganglion cells; RPE: retinal pigment epithelium; SI: sodium iodate; TCA: tricarboxylic acid.

Zinc Supplementation Induced Transcriptional Changes in Primary Human Retinal Pigment Epithelium: A Single-Cell RNA Sequencing Study to Understand Age-Related Macular Degeneration

Zinc supplementation has been shown to be beneficial to slow the progression of age-related macular degeneration (AMD). However, the molecular mechanism underpinning this benefit is not well understood. This study used single-cell RNA sequencing to identify transcriptomic changes induced by zinc supplementation. Human primary retinal pigment epithelial (RPE) cells could mature for up to 19 weeks. After 1 or 18 weeks in culture, we supplemented the culture medium with 125 µM added zinc for one week. RPE cells developed high transepithelial electrical resistance, extensive, but variable pigmentation, and deposited sub-RPE material similar to the hallmark lesions of AMD. Unsupervised cluster analysis of the combined transcriptome of the cells isolated after 2, 9, and 19 weeks in culture showed considerable heterogeneity. Clustering based on 234 pre-selected RPE-specific genes divided the cells into two distinct clusters, we defined as more and less differentiated cells. The proportion of more differentiated cells increased with time in culture, but appreciable numbers of cells remained less differentiated even at 19 weeks. Pseudotemporal ordering identified 537 genes that could be implicated in the dynamics of RPE cell differentiation (FDR < 0.05). Zinc treatment resulted in the differential expression of 281 of these genes (FDR < 0.05). These genes were associated with several biological pathways with modulation of ID1/ID3 transcriptional regulation. Overall, zinc had a multitude of effects on the RPE transcriptome, including several genes involved in pigmentation, complement regulation, mineralization, and cholesterol metabolism processes associated with AMD.