High throughput functional profiling of genes at intraocular pressure loci reveals distinct networks for glaucoma.

INTRODUCTION: Primary open angle glaucoma (POAG) is a leading cause of blindness globally. Characterized by progressive retinal ganglion cell degeneration, the precise pathogenesis remains unknown. Genome-wide association studies (GWAS) have uncovered many genetic variants associated with elevated intraocular pressure (IOP), one of the key risk factors for POAG. We aimed to identify genetic and morphological variation that can be attributed to trabecular meshwork cell (TMC) dysfunction and raised IOP in POAG. METHODS: 62 genes across 55 loci were knocked-out in a primary human TMC line. Each knockout group, including five non-targeting control groups, underwent single-cell RNA-sequencing (scRNA-seq) for differentially-expressed gene (DEG) analysis. Multiplexed fluorescence coupled with CellProfiler image analysis allowed for single-cell morphological profiling. RESULTS: Many gene knockouts invoked DEGs relating to matrix metalloproteinases and interferon-induced proteins. We have prioritized genes at four loci of interest to identify gene knockouts that may contribute to the pathogenesis of POAG, including ANGPTL2, LMX1B, CAV1, and KREMEN1. Three genetic networks of gene knockouts with similar transcriptomic profiles were identified, suggesting a synergistic function in trabecular meshwork cell physiology. TEK knockout caused significant upregulation of nuclear granularity on morphological analysis, while knockout of TRIOBP, TMCO1 and PLEKHA7 increased granularity and intensity of actin and the cell-membrane. CONCLUSION: High-throughput analysis of cellular structure and function through multiplex fluorescent single-cell analysis and scRNA-seq assays enabled the direct study of genetic perturbations at the single-cell resolution. This work provides a framework for investigating the role of genes in the pathogenesis of glaucoma and heterogenous diseases with a strong genetic basis.

A novel GSN variant outside the G2 calcium-binding domain associated with Amyloidosis of the Finnish type.

Gelsolin (GSN) variants have been implicated in amyloidosis of the Finnish type. This case series reports a novel GSN:c.1477T>C,p.(Trp493Arg) variant in a family with ocular and systemic features consistent with Finnish Amyloidosis. Exome sequencing performed on affected individuals from two families manifesting cutis laxa and polymorphic corneal stromal opacities demonstrated the classic GSN:c.654G>A,p.Asp214Asn variant in single affected individual from one family, and a previously undocumented GSN:c.1477T>C variant in three affected first-degree relatives from a separate family. Immunohistochemical studies on corneal tissue from a proband with the c.1477T>C variant identified gelsolin protein within histologically defined corneal amyloid deposits. This study reports a novel association between the predicted pathogenic GSN:c.1477T>C variant and amyloidosis of the Finnish type, and is the first to provide functional evidence of a pathological GSN variant at a locus distant to the critical G2 calcium-binding region, resulting in the phenotype of amyloidosis of the Finnish type.

Genome-wide association analysis of 95 549 individuals identifies novel loci and genes influencing optic disc morphology.

Optic nerve head morphology is affected by several retinal diseases. We measured the vertical optic disc diameter (DD) of the UK Biobank (UKBB) cohort (N = 67 040) and performed the largest genome-wide association study (GWAS) of DD to date. We identified 81 loci (66 novel) for vertical DD. We then replicated the novel loci in International Glaucoma Genetic Consortium (IGGC, N = 22 504) and European Prospective Investigation into Cancer-Norfolk (N = 6005); in general the concordance in effect sizes was very high (correlation in effect size estimates 0.90): 44 of the 66 novel loci were significant at P < 0.05, with 19 remaining significant after Bonferroni correction. We identified another 26 novel loci in the meta-analysis of UKBB and IGGC data. Gene-based analyses identified an additional 57 genes. Human ocular tissue gene expression analysis showed that most of the identified genes are enriched in optic nerve head tissue. Some of the identified loci exhibited pleiotropic effects with vertical cup-to-disc ratio, intraocular pressure, glaucoma and myopia. These results can enhance our understanding of the genetics of optic disc morphology and shed light on the genetic findings for other ophthalmic disorders such as glaucoma and other optic nerve diseases.

Biallelic CPAMD8 Variants Are a Frequent Cause of Childhood and Juvenile Open-Angle Glaucoma.

PURPOSE: Developmental abnormalities of the ocular anterior segment in some cases can lead to ocular hypertension and glaucoma. CPAMD8 is a gene of unknown function recently associated with ocular anterior segment dysgenesis, myopia, and ectopia lentis. We sought to assess the contribution of biallelic CPAMD8 variants to childhood and juvenile open-angle glaucoma. DESIGN: Retrospective, multicenter case series. PARTICIPANTS: A total of 268 probands and their relatives with a diagnosis of childhood or juvenile open-angle glaucoma. PURPOSE: Developmental abnormalities of the ocular anterior segment in some cases can lead to ocular hypertension and glaucoma. CPAMD8 is a gene of unknown function recently associated with ocular anterior segment dysgenesis, myopia, and ectopia lentis. We sought to assess the contribution of biallelic CPAMD8 variants to childhood and juvenile open-angle glaucoma. METHODS: Patients underwent a comprehensive ophthalmic assessment, with DNA from patients and their relatives subjected to genome, exome, or capillary sequencing. CPAMD8 RNA expression analysis was performed on tissues dissected from cadaveric human eyes. MAIN OUTCOME MEASURES: Diagnostic yield within a cohort of childhood and juvenile open-angle glaucoma, prevalence and risk of ophthalmic phenotypes, and relative expression of CPAMD8 in the human eye. RESULTS: We identified rare (allele frequency < 4×10-5) biallelic CPAMD8 variants in 5.7% (5/88) of probands with childhood glaucoma and 2.1% (2/96) of probands with juvenile open-angle glaucoma. When including family members, we identified 11 individuals with biallelic variants in CPAMD8 from 7 unrelated families. Nine of these individuals were diagnosed with glaucoma (9/11, 81.8%), with a mean age at diagnosis of 9.22±14.89 years, and all individuals with glaucoma required 1 or more incisional procedures to control high intraocular pressure. Iris abnormalities were observed in 9 of 11 individuals, cataract was observed in 8 of 11 individuals (72.7%), and retinal detachment was observed in 3 of 11 individuals (27.3%). CPAMD8 expression was highest in neural crest-derived tissues of the adult anterior segment, suggesting that CPAMD8 variation may cause malformation or obstruction of key drainage structures. CONCLUSIONS: Biallelic CPAMD8 variation was associated with a highly heterogeneous phenotype and in our cohorts was the second most common inherited cause of childhood glaucoma after CYP1B1 and juvenile open-angle glaucoma after MYOC. CPAMD8 sequencing should be considered in the investigation of both childhood and juvenile open-angle glaucoma, particularly when associated with iris abnormalities, cataract, or retinal detachment.

New insights into the genetics of primary open-angle glaucoma based on meta-analyses of intraocular pressure and optic disc characteristics.

Primary open-angle glaucoma (POAG), the most common optic neuropathy, is a heritable disease. Siblings of POAG cases have a ten-fold increased risk of developing the disease. Intraocular pressure (IOP) and optic nerve head characteristics are used clinically to predict POAG risk. We conducted a genome-wide association meta-analysis of IOP and optic disc parameters and validated our findings in multiple sets of POAG cases and controls. Using imputation to the 1000 genomes (1000G) reference set, we identified 9 new genomic regions associated with vertical cup-disc ratio (VCDR) and 1 new region associated with IOP. Additionally, we found 5 novel loci for optic nerve cup area and 6 for disc area. Previously it was assumed that genetic variation influenced POAG either through IOP or via changes to the optic nerve head; here we present evidence that some genomic regions affect both IOP and the disc parameters. We characterized the effect of the novel loci through pathway analysis and found that pathways involved are not entirely distinct as assumed so far. Further, we identified a novel association between CDKN1A and POAG. Using a zebrafish model we show that six6b (associated with POAG and optic nerve head variation) alters the expression of cdkn1a. In summary, we have identified several novel genes influencing the major clinical risk predictors of POAG and showed that genetic variation in CDKN1A is important in POAG risk.

Autosomal dominant nanophthalmos and high hyperopia associated with a C-terminal frameshift variant in MYRF.

Purpose: Nanophthalmos is a rare subtype of microphthalmia associated with high hyperopia and an increased risk of angle-closure glaucoma. We investigated the genetic cause of nanophthalmos and high hyperopia in an autosomal dominant kindred. Methods: A proband with short axial length, high hyperopia, and dextrocardia was subjected to exome sequencing. Human and rodent gene expression data sets were used to investigate the expression of relevant genes. Results: We identified a segregating heterozygous frameshift variant at the 3' end of the penultimate exon of MYRF. Using Myc-MYRF chromatin immunoprecipitation data from rat oligodendrocytes, MYRF was found to bind immediately upstream of the transcriptional start site of Tmem98, a gene that itself has been implicated in autosomal dominant nanophthalmos. MYRF and TMEM98 were found to be expressed in the human retina, with a similar pattern of expression across several dissected human eye tissues. Conclusions: C-terminal variants in MYRF, which are expected to escape nonsense-mediated decay, represent a rare cause of autosomal dominant nanophthalmos with or without dextrocardia or congenital diaphragmatic hernia.

Reduced expression of apolipoprotein E and immunoglobulin heavy constant gamma 1 proteins in Fuchs endothelial corneal dystrophy.

BACKGROUND: Fuchs endothelial corneal dystrophy (FECD) is a progressive and potentially a sight threatening disease, and a common indication for corneal grafting in the elderly. Aberrant thickening of Descemet's membrane, formation of microscopic excrescences (guttae) and gradual loss of corneal endothelial cells are the hallmarks of the disease. The aim of this study was to identify differentially abundant proteins between FECD-affected and unaffected Descemet's membrane. METHODS: Label-free quantitative proteomics using nanoscale ultra-performance liquid chromatography-mass spectrometry (nUPLC-MSE ) was employed on affected and unaffected Descemet's membrane extracts, and interesting findings were further investigated using quantitative reverse transcription-polymerase chain reaction and immunohistochemical techniques. RESULTS: Quantitative proteomics revealed significantly lower abundance of apolipoprotein E (APOE) and immunoglobulin heavy constant gamma 1 protein (IGHG1) in affected Descemet's membrane. The difference in the distribution of APOE between affected and unaffected Descemet's membrane and of IGHG1 detected by immunohistochemistry support their down-regulation in the disease. Comparative gene expression analysis showed significantly lower APOE mRNA levels in FECD-affected than unaffected corneal endothelium. IGHG1 gene is expressed at extremely low levels in the corneal endothelium, precluding relative expression analysis. CONCLUSIONS: This is the first study to report comparative proteomics of Descemet's membrane tissue, and implicates dysregulation of APOE and IGHG1 proteins in the pathogenesis of Fuchs endothelial corneal dystrophy.

Novel protein constituents of pathological ocular pseudoexfoliation syndrome deposits identified with mass spectrometry.

Purpose: Pseudoexfoliation (PEX) syndrome is an age-related progressive disease of the extracellular matrix with ocular manifestations. PEX is clinically diagnosed by the presence of extracellular exfoliative deposits on the anterior surface of the ocular lens. PEX syndrome is a major risk factor for developing glaucoma, the leading cause of irreversible blindness in the world, and is often associated with the development of cataract. PEX reportedly coexists with Alzheimer disease and increases the risk of heart disease and stroke. PEX material deposited on the anterior surface of the ocular lens is highly proteinaceous, complex, and insoluble, making deciphering the protein composition of the material challenging. Thus, to date, only a small proportion of the protein composition of PEX material is known. The aim of this study was to decipher the protein composition of pathological PEX material deposited on the ocular lens in patients and advance the understanding of pathophysiology of PEX syndrome. Methods: Liquid-chromatography and tandem mass spectrometry (LC-MS/MS) was employed to discover novel proteins in extracts of neat PEX material surgically isolated from patients (n = 4) with PEX syndrome undergoing cataract surgery. A sub-set of the identified proteins was validated with immunohistochemistry using lens capsule specimens from independent patients (n=3); lens capsules from patients with cataract but without PEX syndrome were used as controls (n=4). Expression of transcripts of the validated proteins in the human lens epithelium was analyzed with reverse transcription PCR (RT-PCR). Functional relationships among the proteins identified in this study and genes and proteins previously implicated in the disease were bioinformatically determined using InnateDB. Results: Peptides corresponding to 66 proteins, including ten proteins previously known to be present in PEX material, were identified. Thirteen newly identified proteins were chosen for validation. Of those proteins, 12 were found to be genuine components of the material. The novel protein constituents include apolipoproteins (APOA1 and APOA4), stress response proteins (CRYAA and PRDX2), and blood-related proteins (fibrinogen and hemoglobin subunits), including iron-free hemoglobin. The gene expression data suggest that the identified stress-response proteins and hemoglobin are contributed by the lens epithelium and apolipoproteins and fibrinogen by the aqueous humor to the PEX material. Pathway analysis of the identified novel protein constituents and genes or proteins previously implicated in the disease reiterated the involvement of extracellular matrix organization and degradation, elastic fiber formation, and complement cascade in PEX syndrome. Network analysis suggested a central role of fibronectin in the pathophysiology of the disease. The identified novel protein constituents of PEX material also shed light on the molecular basis of the association of PEX syndrome with heart disease, stroke, and Alzheimer disease. Conclusions: This study expands the understanding of the protein composition of pathological PEX material deposited on the ocular lens in patients with PEX syndrome and provides useful insights into the pathophysiology of this disease. This study together with the previous study by our group (Sharma et al. Experimental Eye Research 2009;89(4):479-85) demonstrate that using neat PEX material, devoid of the underlying lens capsule, for proteomics analysis is an effective approach for deciphering the protein composition of complex and highly insoluble extracellular pathological ocular deposits present in patients with PEX syndrome.

Meta-analysis of Genome-Wide Association Studies Identifies Novel Loci Associated With Optic Disc Morphology.

Primary open-angle glaucoma is the most common optic neuropathy and an important cause of irreversible blindness worldwide. The optic nerve head or optic disc is divided in two parts: a central cup (without nerve fibers) surrounded by the neuroretinal rim (containing axons of the retinal ganglion cells). The International Glaucoma Genetics Consortium conducted a meta-analysis of genome-wide association studies consisting of 17,248 individuals of European ancestry and 6,841 individuals of Asian ancestry. The outcomes of the genome-wide association studies were disc area and cup area. These specific measurements describe optic nerve morphology in another way than the vertical cup-disc ratio, which is a clinically used measurement, and may shed light on new glaucoma mechanisms. We identified 10 new loci associated with disc area (CDC42BPA, F5, DIRC3, RARB, ABI3BP, DCAF4L2, ELP4, TMTC2, NR2F2, and HORMAD2) and another 10 new loci associated with cup area (DHRS3, TRIB2, EFEMP1, FLNB, FAM101, DDHD1, ASB7, KPNB1, BCAS3, and TRIOBP). The new genes participate in a number of pathways and future work is likely to identify more functions related to the pathogenesis of glaucoma.

CYP1B1 copy number variation is not a major contributor to primary congenital glaucoma.

PURPOSE: To evaluate the prevalence and the diagnostic utility of testing for CYP1B1 copy number variation (CNV) in primary congenital glaucoma (PCG) cases unexplained by CYP1B1 point mutations in The Australian and New Zealand Registry of Advanced Glaucoma. METHODS: In total, 50 PCG cases either heterozygous for disease-causing variants or with no CYP1B1 sequence variants were included in the study. CYP1B1 CNV was analyzed by Multiplex Ligation-dependent Probe Amplification (MLPA). RESULTS: No deletions or duplications were found in any of the cases. CONCLUSION: This is the first study to report on CYP1B1 CNV in PCG cases. Our findings show that this mechanism is not a major contributor to the phenotype and is of limited diagnostic utility.

Predictive genetic testing experience for myocilin primary open-angle glaucoma using the Australian and New Zealand Registry of Advanced Glaucoma.

PURPOSE: Predictive genetic testing of relatives of known myocilin (MYOC) gene mutation carriers is an appropriate strategy to identify individuals at risk for glaucoma. It is likely to prevent irreversible blindness in this high-risk group because this treatable condition might otherwise be diagnosed late. The Australian and New Zealand Registry of Advanced Glaucoma has established genetic testing protocols for known glaucoma genes, including MYOC. METHODS: Through the Australian and New Zealand Registry of Advanced Glaucoma, we investigated the experience of 40 unaffected individuals who had undergone predictive genetic testing for MYOC mutations through questionnaires. RESULTS: The main motivations for being tested were (i) to make appropriate interventions and (ii) to reduce uncertainty. All our respondents perceived strong benefits, either medical or emotional, in being tested. However, different concerns were raised by the respondents that need to be addressed during counseling. Greater family awareness was reported by the majority of the respondents, and the ability to provide information to children was a strong motivation for being tested. CONCLUSION: This study provides valuable information on the personal and familial impacts of having predictive genetic testing for glaucoma, which will help health professionals to better address the issues faced by patients and provide them adequate support.

Deep Learning-Based Identification of Intraocular Pressure-Associated Genes Influencing Trabecular Meshwork Cell Morphology.

Purpose: Genome-wide association studies have recently uncovered many loci associated with variation in intraocular pressure (IOP). Artificial intelligence (AI) can be used to interrogate the effect of specific genetic knockouts on the morphology of trabecular meshwork cells (TMCs) and thus, IOP regulation. Design: Experimental study. Subjects: Primary TMCs collected from human donors. Methods: Sixty-two genes at 55 loci associated with IOP variation were knocked out in primary TMC lines. All cells underwent high-throughput microscopy imaging after being stained with a 5-channel fluorescent cell staining protocol. A convolutional neural network was trained to distinguish between gene knockout and normal control cell images. The area under the receiver operator curve (AUC) metric was used to quantify morphological variation in gene knockouts to identify potential pathological perturbations. Main Outcome Measures: Degree of morphological variation as measured by deep learning algorithm accuracy of differentiation from normal controls. Results: Cells where LTBP2 or BCAS3 had been perturbed demonstrated the greatest morphological variation from normal TMCs (AUC 0.851, standard deviation [SD] 0.030; and AUC 0.845, SD 0.020, respectively). Of 7 multigene loci, 5 had statistically significant differences in AUC (P < 0.05) between genes, allowing for pathological gene prioritization. The mitochondrial channel most frequently showed the greatest degree of morphological variation (33.9% of cell lines). Conclusions: We demonstrate a robust method for functionally interrogating genome-wide association signals using high-throughput microscopy and AI. Genetic variations inducing marked morphological variation can be readily identified, allowing for the gene-based dissection of loci associated with complex traits. Financial Disclosures: Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

RNA Sequencing of Lens Capsular Epithelium Implicates Novel Pathways in Pseudoexfoliation Syndrome.

Purpose: Pseudoexfoliation syndrome (PEX) is a common systemic disease that results in severe and often irreversible vision loss. Despite considerable research effort, PEX remains incompletely understood. This study sought to perform the first RNAseq study in elucidate the pathophysiology of PEX, and contribute a publicly available transcriptomic data resource for future research. Methods: Human ocular lens capsular epithelium samples were collected from 25 patients with PEX and 39 non-PEX controls undergoing cataract surgery. RNA extracted from these specimens was subjected to polyadenylated (mRNA) selection and deep bulk RNA sequencing. Differential expression analysis investigated protein-coding gene transcripts. Exploratory analyses used pathway analysis tools, and curated class- and disease-specific gene sets. Results: Differential expression analysis demonstrated that 2882 genes were differentially expressed according to PEX status. Genes associated with viral gene expression pathways were among the most upregulated, alongside genes encoding ribosomal and mitochondrial respiratory transport chain proteins. Cell adhesion protein transcripts including type 4 collagen subunits were downregulated. Conclusions: This comparative transcriptomic dataset highlights novel and previously recognized pathogenic pathways in PEX and provides the first comprehensive transcriptomic resource, adding an additional layer to build further understanding of PEX pathophysiology.

Multitrait analysis of glaucoma identifies new risk loci and enables polygenic prediction of disease susceptibility and progression.

Glaucoma, a disease characterized by progressive optic nerve degeneration, can be prevented through timely diagnosis and treatment. We characterize optic nerve photographs of 67,040 UK Biobank participants and use a multitrait genetic model to identify risk loci for glaucoma. A glaucoma polygenic risk score (PRS) enables effective risk stratification in unselected glaucoma cases and modifies penetrance of the MYOC variant encoding p.Gln368Ter, the most common glaucoma-associated myocilin variant. In the unselected glaucoma population, individuals in the top PRS decile reach an absolute risk for glaucoma 10 years earlier than the bottom decile and are at 15-fold increased risk of developing advanced glaucoma (top 10% versus remaining 90%, odds ratio = 4.20). The PRS predicts glaucoma progression in prospectively monitored, early manifest glaucoma cases (P = 0.004) and surgical intervention in advanced disease (P = 3.6 × 10-6). This glaucoma PRS will facilitate the development of a personalized approach for earlier treatment of high-risk individuals, with less intensive monitoring and treatment being possible for lower-risk groups.

Rare, potentially pathogenic variants in 21 keratoconus candidate genes are not enriched in cases in a large Australian cohort of European descent.

Many genes have been suggested as candidate genes for keratoconus based on their function, their proximity to associated polymorphisms or due to the identification of putative causative variants within the gene. However, very few of these genes have been assessed for rare variation in keratoconus more broadly. In contrast, VSX1 and SOD1 have been widely assessed, however, the vast majority of studies have been small and the findings conflicting. In a cohort of Australians of European descent, consisting of 385 keratoconus cases and 396 controls, we screened 21 keratoconus candidate genes: BANP, CAST, COL4A3, COL4A4, COL5A1, FOXO1, FNDC3B, HGF, IL1A, IL1B, ILRN, IMMP2L, MPDZ, NFIB, RAB3GAP1, RAD51, RXRA, SLC4A11, SOD1, TF and VSX1. The candidate genes were sequenced in these individuals by either whole exome sequencing or targeted gene sequencing. Variants were filtered to identify rare (minor allele frequency <1%), potentially pathogenic variants. A total of 164 such variants were identified across the two groups with no variants fulfilling these criteria in cases in IL1RN, BANP, IL1B, RAD51 or SOD1. The frequency of variants was compared between cases and controls using chi-square or Fishers' Exact tests for each gene with at least one rare potentially pathogenic variant identified in the case cohort. The number of rare potentially pathogenic variants per gene ranged from three (RXRA) to 102 (MPDZ), however for all genes, there was no difference in the frequency between the cases and controls. We conclude that rare potentially pathogenic variation in the 21 candidate genes assessed do not play a major role in keratoconus susceptibility and pathogenesis.

Prevalence and type of artefact with spectral domain optical coherence tomography macular ganglion cell imaging in glaucoma surveillance.

PURPOSE: The ganglion cell analysis (GCA) of the CIRRUSTM HD-OCT (Carl Zeiss, Meditec; Dublin, CA) provides measurement of the macular ganglion cell-inner plexiform layer (GCIPL) thickness. This study determined the frequency of scan artefacts and errors in GCIPL imaging in individuals undergoing HD-OCT surveillance for glaucoma. METHOD: A total of 1439 eyes from 721 subjects enrolled in a prospective study assessing predictors of glaucoma progression underwent macular GCIPL imaging with the CIRRUS HD-OCT at recruitment. The prevalence of acquisition errors, segmentation errors, and co-morbid macular pathology was determined. RESULTS: A total of 87 (6.0%) of the 1439 scans had either acquisition errors, segmentation artefacts, or other macular pathology. The most common co-morbid macular pathology was epiretinal membrane in 2.2% of eyes. CONCLUSION: The macular GCIPL scan was artefact free in 94% of eyes. However, epiretinal membrane and high myopia can cause scan artefact and should be considered when interpreting the results.

Analysis combining correlated glaucoma traits identifies five new risk loci for open-angle glaucoma.

Open-angle glaucoma (OAG) is a major cause of blindness worldwide. To identify new risk loci for OAG, we performed a genome-wide association study in 3,071 OAG cases and 6,750 unscreened controls, and meta-analysed the results with GWAS data for intraocular pressure (IOP) and optic disc parameters (the overall meta-analysis sample size varying between 32,000 to 48,000 participants), which are glaucoma-related traits. We identified and independently validated four novel genome-wide significant associations within or near MYOF and CYP26A1, LINC02052 and CRYGS, LMX1B, and LMO7 using single variant tests, one additional locus (C9) using gene-based tests, and two genetic pathways - "response to fluid shear stress" and "abnormal retina morphology" - in pathway-based tests. Interestingly, some of the new risk loci contribute to risk of other genetically-correlated eye diseases including myopia and age-related macular degeneration. To our knowledge, this study is the first integrative study to combine genetic data from OAG and its correlated traits to identify new risk variants and genetic pathways, highlighting the future potential of combining genetic data from genetically-correlated eye traits for the purpose of gene discovery and mapping.

Genome-wide association study of intraocular pressure uncovers new pathways to glaucoma.

Intraocular pressure (IOP) is currently the sole modifiable risk factor for primary open-angle glaucoma (POAG), one of the leading causes of blindness worldwide1. Both IOP and POAG are highly heritable2. We report a combined analysis of participants from the UK Biobank (n = 103,914) and previously published data from the International Glaucoma Genetic Consortium (n = 29,578)3,4 that identified 101 statistically independent genome-wide-significant SNPs for IOP, 85 of which have not been previously reported4-12. We examined these SNPs in 11,018 glaucoma cases and 126,069 controls, and 53 SNPs showed evidence of association. Gene-based tests implicated an additional 22 independent genes associated with IOP. We derived an allele score based on the IOP loci and loci influencing optic nerve head morphology. In 1,734 people with advanced glaucoma and 2,938 controls, participants in the top decile of the allele score were at increased risk (odds ratio (OR) = 5.6; 95% confidence interval (CI): 4.1-7.6) of glaucoma relative to the bottom decile.

Rare, Potentially Pathogenic Variants in ZNF469 Are Not Enriched in Keratoconus in a Large Australian Cohort of European Descent.

Purpose: The Zinc Finger Protein 469 (ZNF469) gene has been proposed as a candidate gene for keratoconus due to the association of an upstream polymorphism (rs9938149) with the disease in two independent studies, and the role of the gene in the autosomal recessive disease Brittle Cornea Syndrome. Coding variants in ZNF469 have been assessed for association with keratoconus in several small studies, with conflicting results. We assessed rare, potentially pathogenic variants in ZNF469 for enrichment in keratoconus patients in a cohort larger than all previous studies combined. Methods: ZNF469 was sequenced in 385 Australian keratoconus patients of European descent, 346 population controls, and 230 ethnically matched screened controls by either whole exome sequencing or targeted gene sequencing. The frequency of rare and very rare potentially pathogenic variants was compared between cases and controls using χ2 or Fisher's exact tests and further explored using a gene based test (Sequence Kernel Association Test [SKAT]), weighting on the rarity of variants. Results: A total of 49 rare, including 33 very rare, potentially pathogenic variants were identified across all groups. No enrichment of rare or very rare potentially pathogenic variants in ZNF469 was observed in our cases compared to the control groups following analysis using χ2 or Fisher's exact tests. This finding was further supported by the SKAT results, which found no significant difference in the frequency of variants predicted to be damaging between cases and either control group (P = 0.06). Conclusions: Rare variants in ZNF469 do not contribute to keratoconus susceptibility and do not account for the association at rs9938149.

Contribution of Mutations in Known Mendelian Glaucoma Genes to Advanced Early-Onset Primary Open-Angle Glaucoma.

Purpose: Primary open-angle glaucoma (POAG) and primary congenital glaucoma (PCG) with Mendelian inheritance are caused by mutations in at least nine genes. Utilizing whole-exome sequencing, we examined the disease burden accounted for by these known Mendelian glaucoma genes in a cohort of individuals with advanced early-onset POAG. Methods: The cases exhibited advanced POAG with young age of diagnosis. Cases and examined local controls were subjected to whole-exome sequencing. Nine hundred ninety-three previously sequenced exomes of Australian controls were called jointly with our dataset. Qualifying variants were selected based on predicted pathogenicity and rarity in public domain gene variant databases. Case-control mutational burdens were calculated for glaucoma-linked genes. Results: Two hundred eighteen unrelated POAG participants and 103 nonglaucomatous controls were included in addition to 993 unexamined controls. Fifty-eight participants (26.6%) harbored rare potentially pathogenic variants in known glaucoma genes. Enrichment of qualifying variants toward glaucoma was present in all genes except WDR36, in which controls harbored more variants, and TBK1, in which no qualifying variants were detected in cases or controls. After multiple testing correction, only MYOC showed statistically significant enrichment of qualifying variants (odds ratio [OR] = 16.62, P = 6.31×10-16). Conclusions: Rare, potentially disease-causing variants in Mendelian POAG genes that showed enrichment in our dataset were found in 22.9% of advanced early-onset POAG cases. MYOC variants represented the largest monogenic cause in POAG. The association between WDR36 and POAG was not supported, and the majority of POAG cases did not harbor a potentially disease-causing variant in the remaining Mendelian genes.