Polygenic Risk Scores Driving Clinical Change in Glaucoma.

Glaucoma is a clinically heterogeneous disease and the world's leading cause of irreversible blindness. Therapeutic intervention can prevent blindness but relies on early diagnosis, and current clinical risk factors are limited in their ability to predict who will develop sight-threatening glaucoma. The high heritability of glaucoma makes it an ideal substrate for genetic risk prediction, with the bulk of risk being polygenic in nature. Here, we summarize the foundations of glaucoma genetic risk, the development of polygenic risk prediction instruments, and emerging opportunities for genetic risk stratification. Although challenges remain, genetic risk stratification will significantly improve glaucoma screening and management.

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.

Is the disease risk and penetrance in Leber hereditary optic neuropathy actually low?

Pedigree analysis showed that a large proportion of Leber hereditary optic neuropathy (LHON) family members who carry a mitochondrial risk variant never lose vision. Mitochondrial haplotype appears to be a major factor influencing the risk of vision loss from LHON. Mitochondrial variants, including m.14484T>C and m.11778G>A, have been added to gene arrays, and thus many patients and research participants are tested for LHON mutations. Analysis of the UK Biobank and Australian cohort studies found more than 1 in 1,000 people in the general population carry either the m.14484T>C or the m.11778G>A LHON variant. None of the subset of carriers examined had visual acuity at 20/200 or worse, suggesting a very low penetrance of LHON. Haplogroup analysis of m.14484T>C carriers showed a high rate of haplogroup U subclades, previously shown to have low penetrance in pedigrees. Penetrance calculations of the general population are lower than pedigree calculations, most likely because of modifier genetic factors. This Matters Arising Response paper addresses the Watson et al. (2022) Matters Arising paper, published concurrently in The American Journal of Human Genetics.

Automated AI labeling of optic nerve head enables insights into cross-ancestry glaucoma risk and genetic discovery in >280,000 images from UKB and CLSA.

Cupping of the optic nerve head, a highly heritable trait, is a hallmark of glaucomatous optic neuropathy. Two key parameters are vertical cup-to-disc ratio (VCDR) and vertical disc diameter (VDD). However, manual assessment often suffers from poor accuracy and is time intensive. Here, we show convolutional neural network models can accurately estimate VCDR and VDD for 282,100 images from both UK Biobank and an independent study (Canadian Longitudinal Study on Aging), enabling cross-ancestry epidemiological studies and new genetic discovery for these optic nerve head parameters. Using the AI approach, we perform a systematic comparison of the distribution of VCDR and VDD and compare these with intraocular pressure and glaucoma diagnoses across various genetically determined ancestries, which provides an explanation for the high rates of normal tension glaucoma in East Asia. We then used the large number of AI gradings to conduct a more powerful genome-wide association study (GWAS) of optic nerve head parameters. Using the AI-based gradings increased estimates of heritability by ∼50% for VCDR and VDD. Our GWAS identified more than 200 loci associated with both VCDR and VDD (double the number of loci from previous studies) and uncovered dozens of biological pathways; many of the loci we discovered also confer risk for glaucoma.

Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

Optical Coherence Tomography (OCT) enables non-invasive imaging of the retina and is used to diagnose and manage ophthalmic diseases including glaucoma. We present the first large-scale genome-wide association study of inner retinal morphology using phenotypes derived from OCT images of 31,434 UK Biobank participants. We identify 46 loci associated with thickness of the retinal nerve fibre layer or ganglion cell inner plexiform layer. Only one of these loci has been associated with glaucoma, and despite its clear role as a biomarker for the disease, Mendelian randomisation does not support inner retinal thickness being on the same genetic causal pathway as glaucoma. We extracted overall retinal thickness at the fovea, representative of foveal hypoplasia, with which three of the 46 SNPs were associated. We additionally associate these three loci with visual acuity. In contrast to the Mendelian causes of severe foveal hypoplasia, our results suggest a spectrum of foveal hypoplasia, in part genetically determined, with consequences on visual function.

Correction: Genetic variation affects morphological retinal phenotypes extracted from UK Biobank optical coherence tomography images.

[This corrects the article DOI: 10.1371/journal.pgen.1009497.].

Healthcare professionals' knowledge and attitudes towards polygenic risk testing for glaucoma.

BACKGROUND: Effective clinical implementation of polygenic risk testing for glaucoma relies on healthcare professionals' attitudes and knowledge of the test. Given the emerging applications of the test, it will likely impact a range of healthcare professionals and will require competency in polygenic risk scores concepts for all those involved in patient care. To our knowledge, this is the first study to assess healthcare professionals' views towards polygenic testing for glaucoma. METHODS: An online cross-sectional questionnaire was distributed to healthcare professionals via relevant professional organisations in Australia. The questionnaire assessed experience and confidence with genetic testing, glaucoma and genetic knowledge, recommendations for the tests, and factors affecting the decision. RESULTS: A total of 94 participants completed the questionnaire. The sample was composed of ophthalmologists (36%), optometrists (21%), orthoptists (17%), general practitioners (16%) and clinical geneticists/genetic counsellors (10%). Although familiarity with polygenic risk scores for glaucoma was low overall (11%), the majority reported a positive attitude towards recommending testing based on known risk factors such as family history (91%) and older age (57%). Over 95% indicated that ophthalmologists would be the most appropriate group to order polygenic risk testing and communicate results. The majority felt they would benefit from more training on polygenic risk scores (93%). CONCLUSIONS: Our findings indicated that multiple groups of healthcare professionals were neither familiar nor confident with the concept of glaucoma polygenic risk testing, and identified training and education needs to support the implementation of testing into clinical practice.

Mitochondrial TXNRD2 and ME3 Genetic Risk Scores Are Associated with Specific Primary Open-Angle Glaucoma Phenotypes.

PURPOSE: Genetic variants in regions that include the mitochondrial genes thioredoxin reductase 2 (TXNRD2) and malic enzyme 3 (ME3) are associated with primary open-angle glaucoma (POAG) in genome-wide association studies (GWASs). To assess their clinical impact, we investigated whether TXNRD2 and ME3 genetic risk scores (GRSs) are associated with specific glaucoma phenotypes. DESIGN: Cross-sectional study. PARTICIPANTS: A total of 2617 patients with POAG and 2634 control participants from the National Eye Institute Glaucoma Human Genetics Collaboration Hereditable Overall Operational Database (NEIGHBORHOOD) consortium. METHODS: All POAG-associated single nucleotide polymorphisms (SNPs) in the TXNRD2 and ME3 loci were identified using GWAS data (P < 0.05). Of these, 20 TXNRD2 and 24 ME3 SNPs were selected after adjusting for linkage disequilibrium. The correlation between SNP effect size and gene expression levels was investigated using the Gene-Tissue Expression database. Genetic risk scores were constructed for each individual using the unweighted sum of TXNRD2, ME3, and TXNRD2 + ME3 combined risk alleles. Age- and sex-adjusted odds ratios (ORs) for POAG diagnosis were calculated per decile for each GRS. Additionally, the clinical features of patients with POAG in the top 1%, 5%, and 10% of each GRS were compared with those in the bottom 1%, 5%, and 10%, respectively. MAIN OUTCOME MEASURES: Primary open-angle glaucoma OR per GRS decile, maximum treated intraocular pressure (IOP), and prevalence of paracentral visual field loss among patients with POAG with high versus low GRSs. RESULTS: A larger SNP effect size strongly correlated with higher TXNRD2 and lower ME3 expression levels (r = 0.95 and r = -0.97, respectively; P < 0.05 for both). Individuals in decile 10 of the TXNRD2 + ME3 GRS had the highest odds of POAG diagnosis (OR, 1.79 compared with decile 1; 95% confidence interval, 1.39-2.30; P < 0.001). Patients with POAG in the top 1% of the TXNRD2 GRS showed higher mean maximum treated IOP compared with the bottom 1% (19.9 mmHg vs. 15.6 mmHg; adjusted P = 0.03). Patients with POAG in the top 1% of the ME3 and TXNRD2 + ME3 GRS showed a higher prevalence of paracentral field loss than the bottom 1% (72.7% vs. 14.3% for ME3 GRS and 88.9% vs. 33.3% for TXNRD2+ME3 GRS; adjusted P = 0.03 for both). CONCLUSIONS: Patients with POAG with higher TXNRD2 and ME3 GRSs showed higher treated IOP and a greater prevalence of paracentral field loss. Functional studies exploring how these variants impact mitochondrial function in patients with glaucoma are warranted. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

High Polygenic Risk Is Associated with Earlier Initiation and Escalation of Treatment in Early Primary Open-Angle Glaucoma.

PURPOSE: To assess the association between a glaucoma polygenic risk score (PRS) and treatment outcomes in primary open-angle glaucoma. DESIGN: Prospective, observational cohort study. PARTICIPANTS: Participants from the Progression Risk of Glaucoma: Relevant SNPs with Significant Association Study were divided into a cohort with suspect glaucoma who were treatment naive at enrollment and one with early manifest and suspect glaucoma receiving treatment at enrollment. METHODS: A per-allele weighted glaucoma PRS was calculated for 1107 participants. Multivariable mixed-effects Cox proportional regression analysis assessed the association between PRS and time to commencement of intraocular pressure (IOP)-lowering therapy in 416 patients with suspect glaucoma who were treatment naive at study enrollment. Secondary analysis evaluated the association between PRS and escalation of IOP-lowering therapy among 691 patients with suspect and early manifest glaucoma who were receiving IOP-lowering therapy at enrollment. MAIN OUTCOME MEASURES: Commencement or escalation of IOP-lowering therapy. RESULTS: A higher PRS was associated with a greater risk of commencing IOP-lowering therapy within 5 years (hazard ratio [HR], 1.45 per 1 standard deviation [/SD]; 95% confidence interval [CI], 1.27-1.62; P < 0.001). Participants in the upper population-based quintile showed a 3.3 times greater risk of commencing therapy by 5 years than those in the lowest quintile (HR, 3.30; 95% CI, 1.63-6,70; P < 0.001) and a 5.4 times greater risk of commencing IOP-lowering therapy by 2 years than the those in the lowest quintile (HR, 5.45; 95% CI, 2.08-14.25; P < 0.001). A higher PRS was associated with a greater risk of treatment escalation among patients receiving treatment at enrollment (HR, 1.19/SD; 95% CI, 1.09-1.31; P < 0.001). In combined analysis of all participants, participants in the top population-based quintile were at 2.3 times greater risk of requiring initiation or escalation of IOP-lowering therapy than those in the lowest quintile (HR, 2.33; 95% CI, 1.75-3.01; P < 0.001). CONCLUSIONS: This study demonstrated novel associations between glaucoma polygenic risk and risk of commencement or escalation of IOP-lowering therapy, building on previous work highlighting the potential clinical usefulness of genetic risk stratification in glaucoma. FINANCIAL DISCLOSURE(S): Proprietary or commercial disclosure may be found after the references.

Genetic Risk Assessment of Degenerative Eye Disease (GRADE): study protocol of a prospective assessment of polygenic risk scores to predict diagnosis of glaucoma and age-related macular degeneration.

BACKGROUND: Glaucoma and age-related macular degeneration (AMD) account for a substantial portion of global blindness. Both conditions are highly heritable, with recognised monogenic and polygenic inheritance patterns. Current screening guidelines lack decisive recommendations. Polygenic risk scores (PRS) allow for cost-effective broad population risk stratification for these conditions. The predictive potential of PRS could facilitate earlier diagnosis and treatment, and prevent unnecessary vision loss. METHODS: The Genetic Risk Assessment of Degenerative Eye disease (GRADE) study is a prospective study designed to generate high-quality evidence about the feasibility of PRS to stratify individuals from the general population, enabling identification of those at highest risk of developing glaucoma or AMD. The targeted recruitment is 1000 individuals aged over 50 years, from which blood or saliva samples will be used for genotyping and an individual PRS for glaucoma and AMD will be derived. Individuals with PRS values in the bottom decile (n = 100), top decile (n = 100) and middle 80% (n = 100) for both glaucoma and AMD will undergo a detailed eye examination for glaucoma and/or AMD. DISCUSSION: The primary objective will be to compare the prevalence of glaucoma and AMD cases between low, intermediate, and high PRS risk groups. We expect to find a higher prevalence of both diseases in the high PRS risk group, as compared to the middle and low risk groups. This prospective study will assess the clinical validity of a PRS for glaucoma and AMD in the general Australian population. Positive findings will support the implementation of PRS into clinical practice.

Specifications of the ACMG/AMP variant curation guidelines for myocilin: Recommendations from the clingen glaucoma expert panel.

The standardization of variant curation criteria is essential for accurate interpretation of genetic results and clinical care of patients. The variant curation guidelines developed by the American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) in 2015 are widely used but are not gene specific. To address this issue, the Clinical Genome Resource (ClinGen) Variant Curation Expert Panels (VCEP) have been tasked with developing gene-specific variant curation guidelines. The Glaucoma VCEP was created to develop rule specifications for genes associated with primary glaucoma, including myocilin (MYOC), the most common cause of Mendelian glaucoma. Of the 28 ACMG/AMP criteria, the Glaucoma VCEP adapted 15 rules to MYOC and determined 13 rules not applicable. Key specifications included determining minor allele frequency thresholds, developing an approach to counting probands and segregations, and reviewing functional assays. The rules were piloted on 81 variants and led to a change in classification in 40% of those that were classified in ClinVar, with functional evidence influencing the classification of 18 variants. The standardized variant curation guidelines for MYOC provide a framework for the consistent application of the rules between laboratories, to improve MYOC genetic testing in the management of glaucoma.

Clinical Outcome and Underlying Genetic Cause of Functional Terminal Complement Pathway Deficiencies in a Multicenter UK Cohort.

BACKGROUND: Terminal complement pathway deficiencies often present with severe and recurrent infections. There is a lack of good-quality data on these rare conditions. This study investigated the clinical outcome and genetic variation in a large UK multi-center cohort with primary and secondary terminal complement deficiencies. METHODS: Clinicians from seven UK centers provided anonymised demographic, clinical, and laboratory data on patients with terminal complement deficiencies, which were collated and analysed. RESULTS: Forty patients, median age 19 (range 3-62) years, were identified with terminal complement deficiencies. Ten (62%) of 16 patients with low serum C5 concentrations had underlying pathogenic CFH or CFI gene variants. Two-thirds were from consanguineous Asian families, and 80% had an affected family member. The median age of the first infection was 9 years. Forty-three percent suffered meningococcal serotype B and 43% serotype Y infections. Nine (22%) were treated in intensive care for meningococcal septicaemia. Two patients had died, one from intercurrent COVID-19. Twenty-one (52%) were asymptomatic and diagnosed based on family history. All but one patient had received booster meningococcal vaccines and 70% were taking prophylactic antibiotics. DISCUSSION: The genetic etiology and clinical course of patients with primary and secondary terminal complement deficiency are variable. Patients with low antigenic C5 concentrations require genetic testing, as the low level may reflect consumption secondary to regulatory defects in the pathway. Screening of siblings is important. Only half of the patients develop septicaemia, but all should have a clear management plan.

The phenotypic spectrum of ADAMTSL4-associated ectopia lentis: Additional cases, complications, and review of literature.

Purpose: ADAMTSL4-associated ectopia lentis is a rare autosomal recessive condition that is primarily associated with crystalline lens displacement. However, the prevalence of other ocular and systemic manifestations of this condition is poorly understood. In this study, we summarize the ocular and systemic phenotypic spectrum of this condition. Methods: A cross-sectional case study series of four individuals with biallelic pathogenic or likely pathogenic ADAMTSL4 variants was performed alongside a literature review of individuals with ADAMTSL4-associated ectopia lentis on September 29, 2021. Ocular and systemic findings, complications, and genetic findings of all four individuals were collected and summarized. Results: The phenotypic spectrum across 91 individuals sourced from literature and four individuals from this case study series was highly variable. The main ocular phenotypes included ectopia lentis (95/95, 100%), ectopia lentis et pupillae (18/95, 19%), iris transillumination (13/95, 14%), iridodonesis (12/95, 13%), persistent pupillary membrane (12/95, 13%), and early-onset cataract or lens opacities (12/95, 13%). Anterior segment features other than ectopia lentis appeared to be exclusively associated with biallelic loss of function variants (p<0.001). Pupillary block glaucoma had a prevalence of 1%. Post-lensectomy complications included retinal detachment (6/41, 15%), elevated intraocular pressure (4/41, 10%), and aphakic glaucoma (1/41, 2%). Most individuals were not reported to have had systemic features (69/95, 73%). Conclusions: The clinical phenotype of ADAMTSL4-associated ectopia lentis was summarized and expanded. Clinicians should be aware of the varied ocular phenotype and the risks of retinal detachment, ocular hypertension, and glaucoma in the diagnosis and management of this condition.

The effect of insulin on response to intravitreal anti-VEGF injection in diabetic macular edema in type 2 diabetes mellitus.

OBJECTIVES: To assess whether insulin therapy impacts the effectiveness of anti-vascular endothelial growth factor (anti-VEGF) injection for the treatment of diabetic macular edema (DME) in type 2 diabetes mellitus. METHODS: This was a retrospective multi-center analysis. The best-corrected visual acuity (BCVA) at 12 months, BCVA change, central macular thickness (CMT), CMT change, and cumulative injection number were compared between the insulin and the oral hypoglycemic agent (OHA) groups. RESULTS: The mean final BCVA and CMT improved in both the insulin (N = 137; p < 0.001; p < 0.001, respectively) and the OHA group (N = 61; p = 0.199; p < 0.001, respectively). The two treatment groups were comparable for final BCVA (p = 0.263), BCVA change (p = 0.184), final CMT (p = 0.741), CMT change (p = 0.458), and the cumulative injections received (p = 0.594). The results were comparable between the two groups when stratified by baseline vision (p > 0.05) and baseline HbA1c (p > 0.05). CONCLUSION: Insulin therapy does not alter treatment outcomes for anti-VEGF therapy in DME.

Quantification of localised vascular wedge-shaped defects in glaucoma.

BACKGROUND: Vascular dysfunction plays a considerable role in glaucoma pathogenesis. Previous glaucoma case studies described localised wedge-shaped vascular defects, similar to retinal nerve fibre layer (RNFL) wedge defects. This study investigates the prevalence and quantification of this vessel loss, in relation to primary open angle glaucoma (POAG) parameters. METHODS: This study included 608 eyes (351 participants): 192 PROGRESSA study participants (342 eyes) with suspect, preperimetric or early manifest POAG, observed for vascular wedge defect presence (cohort one); an additional 114 individuals (cohort two-208 eyes) with POAG at various stages of progression for wedge characterisation; and 38 controls (56 eyes). Vascular wedge defects were observed using optical coherence tomography angiography (OCTA). Wedge parameters and vessel densities were quantified using ImageJ software. RNFL and ganglion cell layer inner plexiform layer (GCLIPL) from OCT scans, and mean deviation (Humphrey visual field 24-2) were also assessed. RESULTS: Vascular wedge defects were found in 45/342 eyes (13.2%) in cohort one, in 41/208 eyes (19.7%) in cohort two and were not found in controls. Wedge defects were mostly inferotemporal (80%), and present at all disease stages. They were associated with visual field loss in the opposite hemisphere, thinner RNFL (p < 0.001), thinner GCLIPL (p = 0.003), and focal RNFL loss corresponding with the vascular defect region. CONCLUSION: Vascular wedge defects are present at all POAG stages even before functional change and are strongly concordant with focal RNFL loss. Further research is needed to explore these defects in particular their temporal relationship with clinical measures of POAG.

Identifying Genetic Biomarkers Predicting Response to Anti-Vascular Endothelial Growth Factor Injections in Diabetic Macular Edema.

Intraocular anti-vascular endothelial growth factor (VEGF) therapies are the front-line treatment for diabetic macular edema (DME); however, treatment response varies widely. This study aimed to identify genetic determinants associated with anti-VEGF treatment response in DME. We performed a genome-wide association study on 220 Australian patients with DME treated with anti-VEGF therapy, genotyped on the Illumina Global Screening Array, and imputed to the Haplotype Reference Consortium panel. The primary outcome measures were changes in central macular thickness (CMT in microns) and best-corrected visual acuity (BCVA in ETDRS letters) after 12 months. Association between single nucleotide polymorphism (SNP) genotypes and DME outcomes were evaluated by linear regression, adjusting for the first three principal components, age, baseline CMT/BCVA, duration of diabetic retinopathy, and HbA1c. Two loci reached genome-wide significance (p < 5 × 10−8) for association with increased CMT: a single SNP on chromosome 6 near CASC15 (rs78466540, p = 1.16 × 10−9) and a locus on chromosome 12 near RP11-116D17.1 (top SNP rs11614480, p = 2.69 × 10−8). Four loci were significantly associated with reduction in BCVA: two loci on chromosome 11, downstream of NTM (top SNP rs148980760, p = 5.30 × 10−9) and intronic in RP11-744N12.3 (top SNP rs57801753, p = 1.71 × 10−8); one near PGAM1P1 on chromosome 5 (rs187876551, p = 1.52 × 10−8); and one near TBC1D32 on chromosome 6 (rs118074968, p = 4.94 × 10−8). In silico investigations of each locus identified multiple expression quantitative trait loci and potentially relevant candidate genes warranting further analysis. Thus, we identified multiple genetic loci predicting treatment outcomes for anti-VEGF therapies in DME. This work may potentially lead to managing DME using personalized treatment approaches.

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.

Exome-based investigation of the genetic basis of human pigmentary glaucoma.

BACKGROUND: Glaucoma is a leading cause of visual disability and blindness. Release of iris pigment within the eye, pigment dispersion syndrome (PDS), can lead to one type of glaucoma known as pigmentary glaucoma. PDS has a genetic component, however, the genes involved with this condition are largely unknown. We sought to discover genes that cause PDS by testing cohorts of patients and controls for mutations using a tiered analysis of exome data. RESULTS: Our primary analysis evaluated melanosome-related genes that cause dispersion of iris pigment in mice (TYRP1, GPNMB, LYST, DCT, and MITF). We identified rare mutations, but they were not statistically enriched in PDS patients. Our secondary analyses examined PMEL (previously linked with PDS), MRAP, and 19 other genes. Four MRAP mutations were identified in PDS cases but not in controls (p = 0.016). Immunohistochemical analysis of human donor eyes revealed abundant MRAP protein in the iris, the source of pigment in PDS. However, analysis of MRAP in additional cohorts (415 cases and 1645 controls) did not support an association with PDS. We also did not confirm a link between PMEL and PDS in our cohorts due to lack of reported mutations and similar frequency of the variants in PDS patients as in control subjects. CONCLUSIONS: We did not detect a statistical enrichment of mutations in melanosome-related genes in human PDS patients and we found conflicting data about the likely pathogenicity of MRAP mutations. PDS may have a complex genetic basis that is not easily unraveled with exome analyses.

Non-Synonymous variants in premelanosome protein (PMEL) cause ocular pigment dispersion and pigmentary glaucoma.

Pigmentary glaucoma (PG) is a common glaucoma subtype that results from release of pigment from the iris, called pigment dispersion syndrome (PDS), and its deposition throughout the anterior chamber of the eye. Although PG has a substantial heritable component, no causative genes have yet been identified. We used whole exome sequencing of two independent pedigrees to identify two premelanosome protein (PMEL) variants associated with heritable PDS/PG. PMEL encodes a key component of the melanosome, the organelle essential for melanin synthesis, storage and transport. Targeted screening of PMEL in three independent cohorts (n = 394) identified seven additional PDS/PG-associated non-synonymous variants. Five of the nine variants exhibited defective processing of the PMEL protein. In addition, analysis of PDS/PG-associated PMEL variants expressed in HeLa cells revealed structural changes to pseudomelanosomes indicating altered amyloid fibril formation in five of the nine variants. Introduction of 11-base pair deletions to the homologous pmela in zebrafish by the clustered regularly interspaced short palindromic repeats (CRISPR)-Cas9 method caused profound pigmentation defects and enlarged anterior segments in the eye, further supporting PMEL's role in ocular pigmentation and function. Taken together, these data support a model in which missense PMEL variants represent dominant negative mutations that impair the ability of PMEL to form functional amyloid fibrils. While PMEL mutations have previously been shown to cause pigmentation and ocular defects in animals, this research is the first report of mutations in PMEL causing human disease.

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.