Abstract
Primary open-angle glaucoma (POAG), characterized by retinal ganglion cell death, is a leading cause of irreversible blindness worldwide; however, the 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 GWAS loci and of overlapping eQTLs and sQTLs in 49 GTEx tissues and retina prioritized causal genes for 60% of loci. These genes were enriched in pathways implicated in extracellular matrix organization, cell adhesion, and vascular development. Analysis of single-nucleus RNA-seq of glaucoma-relevant eye tissues revealed that the colocalizing genes and genome-wide POAG and IOP associations were enriched in specific cell types in the aqueous outflow pathways, retina, optic nerve head, peripapillary sclera, and choroid. This study nominated IOP-dependent and independent regulatory mechanisms, genes, and cell types that may contribute to POAG pathogenesis.
Competing Interest Statement
Joshua R. Sanes is a consultant for Biogen, but none of this work benefits or benefited from that relationship. Tavé van Zyl is an employee of Regeneron. Anthony P. Khawaja has acted as a paid consultant or lecturer to Abbvie, Aerie, Allergan, Google Health, Heidelberg Engineering, Novartis, Reichert, Santen and Thea.
Funding Statement
This work was funded by NIH/NEI R01 EY031424-01 (AVS, ARH, JR, PAM), NIH/NEI P30 EY014104 (JLW, AVS), NIH/NEI EY032559-01 (JLW, AVS), and the Chan Zuckerberg Initiative (CZI) Seed Network for the Human Cell Atlas awards CZF2019-002459 (JRS, AVS, WY, AM) and CZF2019-002425 (QL, RC). XJ is supported by the University of Edinburgh and University of Helsinki joint PhD studentship program in Human Genomics. VV is supported by an MRC University Unit Programme grant (MC_UU_00007/10) (QTL in Health and Disease). APK is supported by a UK Research and Innovation Future Leaders Fellowship, an Alcon Research Institute Young Investigator Award and a Lister Institute for Preventive Medicine Award. This research was supported by the NIHR Biomedical Research Centre at Moorfields Eye Hospital and the UCL Institute of Ophthalmology.
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
This work used publicly available summary statistics data of genome-wide association studies of complex human traits. The eye tissues samples subjected to single nucleus RNA-sequencing that was analyzed in this study were collected from de-identified, deceased donors under HIPAA Privacy Rules. The institutional review board (IRB) approvals for post-mortem tissue acquisition at Mass General Brigham (Rapid Autopsy Program) and University of Utah (partnership with Lions International) gave ethical approval for research with the tissue samples. No IRB approval is required for transcriptomic analysis of post-mortem tissues.
I confirm that all necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived, and that any patient/participant/sample identifiers included were not known to anyone (e.g., hospital staff, patients or participants themselves) outside the research group so cannot be used to identify individuals.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Footnotes
In the revised manuscript, we have performed multiple new analyses, providing additional support for the regulatory mechanisms, genes and cell types that may contribute to POAG risk in an IOP-dependent or independent manner, and to IOP regulation, thus strengthening the conclusions of our paper. These added analyses include (i) two-sample Mendelian Randomization analysis of all significantly colocalizing eQTLs and sQTLs with primary open angle glaucoma (POAG) and intraocular pressure (IOP) GWAS loci, (ii) integration of retinal genome topology (Hi-C) and epigenetic data with colocalizing e/sQTLs to provide additional support for regulatory variants and genes that may affect POAG in the retina, (iii) cell type enrichment analysis of POAG and IOP colocalizing e/sGenes using an updated single nucleus expression atlas of the optic nerve head, optic nerve and surrounding posterior tissues, including the peripapillary sclera, peripheral sclera and choroid, generated by the Sanes lab (Monavarfeshani*, Yan* et al., bioRxiv 2023, in press, PNAS 2023), (iv) cell type enrichment analysis of POAG-only loci, IOP-only loci or shared loci in anterior and posterior eye tissue single nucleus expression data, using ECLISPER, (v) cell type enrichment of POAG and IOP GWAS associations, using two additional methods, S-LDSC and MAGMA, to increase confidence in the ECLIPSER POAG and IOP cell type enrichment results, and (vi) conditional analysis, implemented in MAGMA, used to test whether the different cell type enrichment signals for POAG and IOP in the eye tissue single cell data were independent of each other in each ocular tissue. Our revised manuscript provides new insights into the regulatory effects, genes and ocular cell types that may contribute to POAG, a leading cause of blindness worldwide, and its major risk factor, IOP, which might help guide development of new IOP reducing drugs and neuroprotective treatments.
Data availability
All GTEx protected data are available through the database of Genotypes and Phenotypes (dbGaP) (accession no. phs000424.v8). The GTEx eQTL and sQTL and EyeGEx retina eQTL summary statistics are available on the GTEx portal (https://gtexportal.org/home/datasets). The snRNA-seq data for the anterior segment and macula are available in Gene Expression Omnibus (GEO) accession number GSE199013, for the optic nerve head and posterior tissues in GSE236566, and for the retina in GSE226108. The processed data of the anterior and posterior segments can be visualized in the Broad Institute’s Single Cell Portal at https://singlecell.broadinstitute.org/single_cell/study/SCP1841 and https://singlecell.broadinstitute.org/single_cell/study/SCP2298. The retina Hi-C data is accessible in GEO accession number GSE202471. The GWAS summary statistics for the POAG cross-ancestry GWAS meta-analysis and European subset meta-analysis are accessible in GEO under accession numbers GCST90011770 and GCST90011766, respectively, and for IOP are available from the corresponding publication (Khawaja et al., Nature Genetics 2018). The GWAS loci for complex traits analyzed in this study were downloaded from Open Targets Genetics (https://genetics.opentargets.org/). The gene sets taken from MSigDB were downloaded from: http://www.gsea-msigdb.org/gsea/msigdb/collections.jsp, and the mouse phenotype ontology gene sets from the Mouse Genome Informatics (MGI) website (http://www.informatics.jax.org/).
All data produced in the current study are either contained in the manuscript, deposited in online databases as stated in the manuscript, or can be obtained upon request to the authors. All GTEx protected data are available through the database of Genotypes and Phenotypes (dbGaP) (accession no. phs000424.v8). The single nucleus RNA-sequencing data of the anterior and posterior eye tissues are available in Gene Expression Omnibus (GEO) accession numbers GSE199013, GSE236566, and GSE226108. The retina Hi-C data is accessible in GEO accession number GSE202471. The GWAS summary statistics for the POAG cross-ancestry GWAS meta-analysis and European subset meta-analysis are accessible in GEO accession numbers GCST90011770 and GCST90011766, respectively, and for IOP GWAS are available from the corresponding publication (Khawaja et al., Nature Genetics 2018).
