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

Human cytomegalovirus infection coopts chromatin organization to diminish TEAD1 transcription factor activity

Human cytomegalovirus (HCMV) infects up to 80% of the world's population. Here, we show that HCMV infection leads to widespread changes in human chromatin accessibility and chromatin looping, with hundreds of thousands of genomic regions affected 48 hours after infection. Integrative analyses reveal HCMV-induced perturbation of Hippo signaling through drastic reduction of TEAD1 transcription factor activity. We confirm extensive concordant loss of TEAD1 binding, active H3K27ac histone marks, and chromatin looping interactions upon infection. Our data position TEAD1 at the top of a hierarchy involving multiple altered important developmental pathways. HCMV infection reduces TEAD1 activity through four distinct mechanisms: closing of TEAD1-bound chromatin, reduction of YAP1 and phosphorylated YAP1 levels, reduction of TEAD1 transcript and protein levels, and alteration of TEAD1 exon-6 usage. Altered TEAD1-based mechanisms are highly enriched at genetic risk loci associated with eye and ear development, providing mechanistic insight into HCMV's established roles in these processes.

Rodent genetically modified models of glaucoma

Glaucoma, one of the leading causes of irreversible blindness worldwide, is a complex and heterogenous disease. While environmental factors are important, it is well-recognized that the disease has a strong heritable component. With the advent of large-cohort genome wide association studies, a myriad of genetic risk loci has been linked to different forms of glaucoma. Animal models have been an indispensable tool in characterizing these loci, especially if they lie within coding regions in the genome. Not only do these models connect genotype to phenotype, advancing our understanding of glaucoma pathogenesis in the process, they also have valuable utility as a platform for the pre-clinical testing of potential therapies. In this review, we will outline genetic models used for studying the major forms of glaucoma, including primary open angle glaucoma, normal tension glaucoma, primary angle closure glaucoma, pigmentary glaucoma, pseudoexfoliation glaucoma, and early onset glaucoma, including congenital and developmental glaucoma, and how studying these models have helped shed light on human glaucoma.

A multi-cohort genome-wide association study in African ancestry individuals reveals risk loci for primary open-angle glaucoma

Primary open-angle glaucoma (POAG), the leading cause of irreversible blindness worldwide, disproportionately affects individuals of African ancestry. We conducted a genome-wide association study (GWAS) for POAG in 11,275 individuals of African ancestry (6,003 cases; 5,272 controls). We detected 46 risk loci associated with POAG at genome-wide significance. Replication and post-GWAS analyses, including functionally informed fine-mapping, multiple trait co-localization, and in silico validation, implicated two previously undescribed variants (rs1666698 mapping to DBF4P2; rs34957764 mapping to ROCK1P1) and one previously associated variant (rs11824032 mapping to ARHGEF12) as likely causal. For individuals of African ancestry, a polygenic risk score (PRS) for POAG from our mega-analysis (African ancestry individuals) outperformed a PRS from summary statistics of a much larger GWAS derived from European ancestry individuals. This study quantifies the genetic architecture similarities and differences between African and non-African ancestry populations for this blinding disease.

Genome-Wide Meta-analysis Identifies Risk Loci and Improves Disease Prediction of Age-Related Macular Degeneration

PURPOSE

To identify age-related macular degeneration (AMD) risk loci and to establish a polygenic prediction model.

DESIGN

Genome-wide association study (GWAS) and polygenic risk score (PRS) construction.

PARTICIPANTS

We included 64 885 European patients with AMD and 568 740 control participants (with overlapped samples) in the UK Biobank, Genetic Epidemiology Research on Aging (GERA), International AMD Consortium, FinnGen, and published early AMD GWASs in meta-analyses, as well as 733 European patients with AMD and 20 487 control participants from the Canadian Longitudinal Study on Aging (CLSA) and non-Europeans from the UK Biobank and GERA for polygenic risk score validation.

METHODS

A multitrait meta-analysis of GWASs comprised 64 885 patients with AMD and 568 740 control participants; the multitrait approach accounted for sample overlap. We constructed a PRS for AMD based on both previously reported as well as unreported AMD loci. We applied the PRS to nonoverlapping data from the CLSA.

MAIN OUTCOME MEASURES

We identified several single nucleotide polymorphisms associated with AMD and established a PRS for AMD risk prediction.

RESULTS

We identified 63 AMD risk loci alongside the well-established AMD loci CFH and ARMS2, including 9 loci that were not reported in previous GWASs, some of which previously were linked to other eye diseases such as glaucoma (e.g., HIC1). We applied our PRS to nonoverlapping data from the CLSA. A new PRS was constructed using the PRS method, PRS-CS, and significantly improved the prediction accuracy of AMD risk compared with PRSs from previously published datasets. We further showed that even people who carry all the well-known AMD risk alleles at CFH and ARMS2 vary considerably in their AMD risk (ranging from close to 0 in individuals with low PRS to > 50% in individuals with high PRS). Although our PRS was derived in individuals of European ancestry, the PRS shows potential for predicting risk in people of East Asian, South Asian, and Latino ancestry.

CONCLUSIONS

Our findings improve the knowledge of the genetic architecture of AMD and help achieve better accuracy in AMD prediction.

FINANCIAL DISCLOSURE(S)

Proprietary or commercial disclosure may be found in the Footnotes and Disclosures at the end of this article.

Psychological stress induces moderate pathology in the ganglion cell layer in mice

Purpose

Primary open-angle glaucoma (POAG) is a condition with unclear pathogenesis. Researchers have observed an increased incidence of young Chinese POAG patients who manifest significant psychological stress while their intraocular pressure (IOP) is normal or close to normal; we hypothesize that psychological stress may play a causal role in initiating POAG.

Methods

Twenty-four male C57BL/6 mice were included and divided randomly into two groups. A chronic unpredictable mild stress (CUMS) mouse model was established to evaluate the effect of psychological stress on glaucoma-related retinal pathologies. Body weight and IOP were recorded weekly. At 5 weeks after the CUMS procedure, a behavior test, serum corticosterone level, retinal nerve fiber layer (RNFL) thickness, retinal ganglion cell (RGC) number and neurotrophic factor expression were evaluated and compared between the CUMS group and the control group.

Results

CUMS exposure induced depression-like behaviors, lighter body weight, and increased serum corticosterone levels in mice. RNFL thinning and neural cell loss in the ganglion cell layer (GCL) were observed in CUMS mice without significant IOP elevation. Decreased mRNA expression and protein levels of neurotropic factors in retinas of CUMS mice were observed, especially brain-derived neurotrophic factor (BDNF).

Conclusions

The CUMS mouse model demonstrated that psychological stress induced glaucoma-like changes in the retinas of CUMS mice. The mechanism by which psychological stress induces retina defects may be due to a reduced expression of retinal neurotropic factors. Thus, we conclude that psychological stress is causally associated with POAG.

Genetic variants associated with glaucomatous visual field loss in primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is characterized by a progressive optic neuropathy with visual field loss. To investigate the genetic variants associated with visual field loss in POAG, Japanese POAG patients (n = 426) and control subjects (n = 246) were genotyped for 22 genetic variants predisposing to POAG that can be classified into those associated with intraocular pressure (IOP) elevation (IOP-related genetic variants) and optic nerve vulnerability independent of IOP (optic nerve-related genetic variants). The genetic risk score (GRS) of the 17 IOP-related and five optic nerve-related genetic variants was calculated, and the associations between the GRS and the mean deviation (MD) of automated static perimetry as an indicator of the severity of visual field loss and pattern standard deviation (PSD) as an indicator of the focal disturbance were evaluated. There was a significant association (Beta = - 0.51, P = 0.0012) between the IOP-related GRS and MD. The severity of visual field loss may depend on the magnitude of IOP elevation induced by additive effects of IOP-related genetic variants. A significant association (n = 135, Beta = 0.65, P = 0.0097) was found between the optic nerve-related, but not IOP-related, GRS and PSD. The optic nerve-related (optic nerve vulnerability) and IOP-related (IOP elevation) genetic variants may play an important role in the focal and diffuse visual field loss respectively. To our knowledge, this is the first report to show an association between additive effects of genetic variants predisposing to POAG and glaucomatous visual field loss, including severity and focal/diffuse disturbance of visual field loss, in POAG.

Development and Evaluation of a Rules-based Algorithm for Primary Open-Angle Glaucoma in the VA Million Veteran Program

The availability of electronic health record (EHR)-linked biobank data for research presents opportunities to better understand complex ocular diseases. Developing accurate computable phenotypes for ocular diseases for which gold standard diagnosis includes imaging remains inaccessible in most biobank-linked EHRs. The objective of this study was to develop and validate a computable phenotype to identify primary open-angle glaucoma (POAG) through accessing the Department of Veterans Affairs (VA) Computerized Patient Record System (CPRS) and Million Veteran Program (MVP) biobank. Accessing CPRS clinical ophthalmology data from VA Medical Center Eye Clinic (VAMCEC) patients, we developed and iteratively refined POAG case and control algorithms based on clinical, prescription, and structured diagnosis data (ICD-CM codes). Refinement was performed via detailed chart review, initially at a single VAMCEC (n = 200) and validated at two additional VAMCECs (n = 100 each). Positive and negative predictive values (PPV, NPV) were computed as the proportion of CPRS patients correctly classified with POAG or without POAG, respectively, by the algorithms, validated by ophthalmologists and optometrists with access to gold-standard clinical diagnosis data. The final algorithms performed better than previously reported approaches in assuring the accuracy and reproducibility of POAG classification (PPV >83% and NPV >97%) with consistent performance in Black or African American and in White Veterans. Applied to the MVP to identify cases and controls, genetic analysis of a known POAG-associated locus further validated the algorithms. We conclude that ours is a viable approach to use combined EHR-genetic data to study patients with complex diseases that require imaging confirmation.

The genetics of glaucoma: Disease associations, personalised risk assessment and therapeutic opportunities-A review

Glaucoma refers to a heterogenous group of disorders characterised by progressive loss of retinal ganglion cells and associated visual field loss. Both early-onset and adult-onset forms of the disease have a strong genetic component. Here, we summarise the known genetic associations for various forms of glaucoma and the possible functional roles for these genes in disease pathogenesis. We also discuss efforts to translate genetic knowledge into clinical practice, including gene-based tests for disease diagnosis and risk-stratification as well as gene-based therapies.

Association analysis of variants rs35934224 in TXNRD2 and rs6478746 in LMX1B in primary angle-closure and pseudoexfoliation glaucoma

OBJECTIVE

Previous genome-wide studies have demonstrated significant pathogenic association between variants rs35934224 within TXNRD2 and rs6478746 near LMX1B in primary open-angle glaucoma. We investigated the association between these variants in primary angle-closure glaucoma (PACG) and pseudoexfoliation glaucoma (PXG) patients of Saudi origin.

METHODS

In a case-control study, DNA samples from 249 controls (135 men and 114 women), 100 PACG cases (44 men and 56 women), and 95 PXG cases (61 men and 34 women) were genotyped by TaqMan^® based real-time PCR. Statistical tests were performed to evaluate genetic association with glaucoma types and related clinical indices.

RESULTS

The allele frequencies of rs35934224 and rs6478746 did not show significant variation in PACG and PXG than controls, except that the rs35934224[T] allele was found to be significantly low among PXG women (0.10) as compared to controls (0.21) (odds ratio = 0.38, 95% confidence interval = 0.16-0.94, p = 0.024). Rs35934224 genotypes showed a nominal-to-borderline protective association with PACG and PXG among women in different genetic models. However, except for the over-dominant model in PACG (p = 0.0095), none of the effects survived Bonferroni's correction (p < 0.01). Rs6478746 showed no significant genotype or allelic association with PACG and PXG. Regression analysis showed no influence on disease outcome, and neither showed any correlation with intraocular pressure and cup/disk ratio in both PACG and PXG.

CONCLUSIONS

Variants rs35934224 in TXNRD2 and rs6478746 near LMX1B are not associated with PACG and PXG in the Saudi cohort, but rs35934224 may confer modest protection among women. Further population-based studies are needed to validate these results.

Lack of Association Between Polymorphisms in TXNRD2 and LMX1B and Primary Open-Angle Glaucoma in a Saudi Cohort

Objective: Recent studies have demonstrated an association of single nucleotide polymorphisms (SNPs) rs35934224 in TXNRD2 and rs6478746 near LMX1B genes in primary open-angle glaucoma (POAG) among Europeans. We performed a retrospective, case-control study to investigate the association between the rs35934224 (TXNRD2) and rs6478746 (LMX1B) and POAG in a middle-eastern population from Saudi Arabia. Methods: DNA from 399 participants consisting of 150 POAG cases (83 males and 67 females) and 249 controls (135 males and 114 females) were genotyped using TaqMan® real-time PCR. Statistical tests were performed to evaluate genetic association with POAG and related clinical indices. Results: The minor allele frequency (MAF) of rs35934224[T] was 0.19 and 0.20 in POAG and controls, respectively. The difference was non-significant (odds ratio [OR] = 1.08, 95% confidence interval [CI] = 0.75-1.55, p = 0.663). Likewise, rs6478746[G] MAF was 0.12 in both cases and controls with no statistical significance (OR = 1.02, 95% CI = 0.67-1.56, p = 0.910). Genotype analysis showed no association with POAG for both the SNPs in combined and gender-stratified groups. Regression analysis showed no significant effect of risk factors such as age, sex, rs35934224, and rs6478746 genotypes on POAG outcome. Furthermore, both the SNPs showed no significant genotype effect on clinical indices such as intraocular pressure (IOP) and cup/disc ratio in POAG patients. Conclusions: Rs35934224 in TXNRD2 and rs6478746 near LMX1B genes are not associated with POAG or related clinical indices such as IOP and cup/disc ratio in a Saudi cohort. Since the study is limited by sample size further investigations are needed to confirm these results in a larger cohort.

LMX1B Locus Associated with Low-Risk Baseline Glaucomatous Features in the POAAGG Study

Primary open-angle glaucoma (POAG) is the leading cause of irreversible blindness worldwide and has been associated with multiple genetic risk factors. The LMX1B gene is a genetic susceptibility factor for POAG, and several single-nucleotide polymorphisms (SNPs) were shown to be associated with POAG in our own prior Primary Open-Angle African American Glaucoma Genetics (POAAGG) study genome-wide association study (GWAS). This study evaluated the association of the LMX1B locus with baseline optic disc and clinical phenotypic characteristics of glaucoma patients from our African American cohort. Compared to the GG genotype in SNP rs187699205, the GC genotype in this SNP was found to be significantly associated with a smaller cup-to-disc ratio (CDR) and increased (better) visual field mean deviation (MD) in glaucoma cases. None of the glaucoma cases with the GC genotype had disc hemorrhages, disc notching, or beanpot disc appearance. In conclusion, glaucoma phenotypes differed significantly by LMX1B variant in African American patients with POAG, and a SNP variant was associated with certain disease features considered lower risk.

MicroRNAs in the aqueous humor of patients with different types of glaucoma

Purpose

The aim of the study was to compare the frequency and the level of expression of selected miRNAs in the aqueous humor of patients with various types of glaucoma.

Methods

The studied group consisted of 42 patients with glaucoma: 19 with primary open-angle glaucoma (POAG), 14 with pseudoexfoliation glaucoma (PEXG), 9 with primary angle closure glaucoma (PACG), and the control group of 36 patients with senile cataract without glaucoma. The real-time polymerase chain reaction method was used to analyze the expression of miRNAs.

Results

There were no significant differences in the frequency and the level of miRNA expression between various types of glaucoma. There was a tendency for hsa-miR-6722-3p and hsa-miR-184 to be expressed more frequently in PEXG and hsa-miR-1260b in POAG. The expression levels of hsa-miR-1260b and hsa-miR-6515-3p were correlated with age in POAG. Target annotation and functional analyses showed that genes targeted by the most frequently expressed miRNAs (hsa-miR-1202, -1260b, -184, -187-5p, -6515-3p, -6722-3p, and hsa-mir-4634) are involved mainly in response to hypoxia, cardiovascular system development, and apoptosis.

Conclusion

Hsa-miR-1260b was the most abundantly expressed among studied miRNAs and may be a potential biomarker of clinical status in PEXG and PACG.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00417-021-05214-z.

Genome-wide meta-analysis identifies 127 open-angle glaucoma loci with consistent effect across ancestries

Primary open-angle glaucoma (POAG), is a heritable common cause of blindness world-wide. To identify risk loci, we conduct a large multi-ethnic meta-analysis of genome-wide association studies on a total of 34,179 cases and 349,321 controls, identifying 44 previously unreported risk loci and confirming 83 loci that were previously known. The majority of loci have broadly consistent effects across European, Asian and African ancestries. Cross-ancestry data improve fine-mapping of causal variants for several loci. Integration of multiple lines of genetic evidence support the functional relevance of the identified POAG risk loci and highlight potential contributions of several genes to POAG pathogenesis, including SVEP1, RERE, VCAM1, ZNF638, CLIC5, SLC2A12, YAP1, MXRA5, and SMAD6. Several drug compounds targeting POAG risk genes may be potential glaucoma therapeutic candidates.

Genetic background modifies vulnerability to glaucoma-related phenotypes in Lmx1b mutant mice

Variants in the LIM homeobox transcription factor 1-beta (LMX1B) gene predispose individuals to elevated intraocular pressure (IOP), a key risk factor for glaucoma. However, the effect of LMX1B mutations varies widely between individuals. To better understand the mechanisms underlying LMX1B-related phenotypes and individual differences, we backcrossed the Lmx1b^V265D (also known as Lmx1b^Icst ) allele onto the C57BL/6J (B6), 129/Sj (129), C3A/BLiA-Pde6b⁺ /J (C3H) and DBA/2J-Gpnmb⁺ (D2-G) mouse strain backgrounds. Strain background had a significant effect on the onset and severity of ocular phenotypes in Lmx1b^V265D/+ mutant mice. Mice of the B6 background were the most susceptible to developing abnormal IOP distribution, severe anterior segment developmental anomalies (including malformed eccentric pupils, iridocorneal strands and corneal abnormalities) and glaucomatous nerve damage. By contrast, Lmx1b^V265D mice of the 129 background were the most resistant to developing anterior segment abnormalities, had less severe IOP elevation than B6 mutants at young ages and showed no detectable nerve damage. To identify genetic modifiers of susceptibility to Lmx1b^V265D -induced glaucoma-associated phenotypes, we performed a mapping cross between mice of the B6 (susceptible) and 129 (resistant) backgrounds. We identified a modifier locus on Chromosome 18, with the 129 allele(s) substantially lessening severity of ocular phenotypes, as confirmed by congenic analysis. By demonstrating a clear effect of genetic background in modulating Lmx1b-induced phenotypes, providing a panel of strains with different phenotypic severities and identifying a modifier locus, this study lays a foundation for better understanding the roles of LMX1B in glaucoma with the goal of developing new treatments.

MLIP genotype as a predictor of pharmacological response in primary open-angle glaucoma and ocular hypertension

Predicting the therapeutic response to ocular hypotensive drugs is crucial for the clinical treatment and management of glaucoma. Our aim was to identify a possible genetic contribution to the response to current pharmacological treatments of choice in a white Mediterranean population with primary open-angle glaucoma (POAG) or ocular hypertension (OH). We conducted a prospective, controlled, randomized, partial crossover study that included 151 patients of both genders, aged 18 years and older, diagnosed with and requiring pharmacological treatment for POAG or OH in one or both eyes. We sought to identify copy number variants (CNVs) associated with differences in pharmacological response, using a DNA pooling strategy of carefully phenotyped treatment responders and non-responders, treated for a minimum of 6 weeks with a beta-blocker (timolol maleate) and/or prostaglandin analog (latanoprost). Diurnal intraocular pressure reduction and comparative genome wide CNVs were analyzed. Our finding that copy number alleles of an intronic portion of the MLIP gene is a predictor of pharmacological response to beta blockers and prostaglandin analogs could be used as a biomarker to guide first-tier POAG and OH treatment. Our finding improves understanding of the genetic factors modulating pharmacological response in POAG and OH, and represents an important contribution to the establishment of a personalized approach to the treatment of glaucoma.

Association of Glaucoma Risk Genes with Retinal Nerve Fiber Layer in a Multi-ethnic Asian Population: The Singapore Epidemiology of Eye Diseases Study

Purpose

Genome-wide association studies have identified several genes associated with glaucoma. However, their roles in the pathogenesis of glaucoma remain unclear, particularly their effects on retinal nerve fiber layer (RNFL) thickness. The aim of this study was to investigate the associations between the identified glaucoma risk genes and RNFL thickness.

Methods

A total of 3843 participants (7,020 healthy eyes) were enrolled from the Singapore Epidemiology of Eye Diseases (SEED) study, a population-based study composing of three major ethnic groups—Malay, Indian, and Chinese—in Singapore. Ocular examinations were performed, and spectral-domain optical coherence tomography (SD-OCT) was used to measure circumpapillary RNFL thickness. We selected 35 independent glaucoma-associated genetic loci for analysis. An linear regression model was conducted to determine the association of these variants with circumpapillary RNFL, assuming an additive genetic model. We conducted association analysis in each of the three ethnic groups, followed by a meta-analysis of them.

Results

The mean age of the included participants was 59.4 ± 8.9 years, and the mean RFNL thickesss is 92.3 ± 11.2 µm. In the meta-analyses, of the 35 glacuoma loci, we found that only SIX6 was significantly associated with reduction in global RNFL thickness (rs33912345; β = −1.116 um per risk allele, P = 1.64E-05), and the effect size was larger in the inferior RNFL quadrant (β = −2.015 µm, P = 2.9E-6), and superior RNFL quadrant (β = −1.646 µm, P = 6.54E-5). The SIX6 association were consistently observed across all three ethnic groups. Other than RNFL, we also found several genetic varaints associated with vertical cuo-to-disc ratio (ATOH7, CDKN2B-AS1, and TGFBR3-CDC7), rim area (SIX6 and CDKN2B-AS1), and disc area (SIX6, ATOH7, and TGFBR3-CDC7). The association of SIX6 rs33912345 with NRFL thickness remained similar after further adjusting for disc area and 3 other disc parameter associated SNPs (ATOH7, CDKN2B-AS1, and TGFBR3-CDC7).

Conclusions

Of the 35 glaucoma identified risk loci, only SIX6 is significantly and independently associated with thinner RNFL. Our study further supports the involvement of SIX6 with RNFL thickness and pathogensis of glaucoma.

Insights into the regulatory molecules involved in glaucoma pathogenesis

Glaucoma is an important cause of irreversible blindness, characterized by optic nerve anomalies. Increased intraocular pressure (IOP) and aging are major risk factors. Retinal ganglion cells and trabecular meshwork cells are certainly involved in the etiology of glaucoma. Glaucoma is usually a complex disease, and various genes and functions may contribute to its etiology. Among these may be genes that encode regulatory molecules. In this review, regulatory molecules including 18 transcription factors (TFs), 195 microRNAs (miRNAs), 106 long noncoding RNAs (lncRNAs), and two circular RNAs (circRNAs) that are reasonable candidates for having roles in glaucoma pathogenesis are described. The targets of the regulators are reported. Glaucoma-related features including apoptosis, stress responses, immune functions, ECM properties, IOP, and eye development are affected by the targeted genes. The targeted genes that are frequently targeted by multiple regulators most often affect apoptosis and the related features of cell death and cell survival. BCL2, CDKN1A, and TP53 are among the frequent targets of three types of glaucoma-relevant regulators, TFs, miRNAs, and lncRNAs. TP53 was itself identified as a glaucoma-relevant TF. Several of the glaucoma-relevant TFs are themselves among frequent targets of regulatory molecules, which is consistent with existence of a complex network involved in glaucoma pathogenesis.

[The role of metalloproteinases in the development of primary open-angle glaucoma]

The article reviews literature on the role of matrix metalloproteinases (MMP) in the development of eye pathologies, specifically primary open-angle glaucoma (POAG). Regulation of the extracellular matrix is carried out by proteolytic enzymes - metalloproteinases in particular - as well as specific inhibitors of their activity (tissue metalloproteinases inhibitors). The review also reveals its involvement in the pathogenesis of various types of glaucoma, decrease in stability of the intraocular liquid outflow, proves the role of MMR in the ganglionic apoptosis, remodeling of the optic disk and change of lamina cribrosa in primary open-angle glaucoma. There is research devoted to studying the interrelations of MMP genes with POAG development. Associations of single nucleotide polymorphic loci of MMP-1, MMP-2, MMP-9 genes with POAG are shown.

Genetic Variants Associated With the Onset and Progression of Primary Open-Angle Glaucoma

PURPOSE

We sought to investigate the genetic variants associated with the onset and progression of primary open-angle glaucoma (POAG).

DESIGN

Case-control genetic association study.

METHODS

Japanese POAG patients (n = 505) and control subjects (n = 246) were genotyped for 22 genetic variants predisposing to POAG that can be classified into those associated with intraocular pressure (IOP) elevation (IOP-related genetic variants) and optic nerve vulnerability independent of IOP (non-IOP-related genetic variants). The total number of risk alleles of the 17 IOP-related and 5 non-IOP-related genetic variants were calculated as the genetic risk score (GRS), and the associations between the GRS and family history of glaucoma as an indicator of POAG onset and age at the diagnosis of glaucoma as an indicator of POAG progression were evaluated.

RESULTS

There was a significant association (P = .014; odds ratio 1.26 per GRS) between the non-IOP-related GRS, but not IOP-related GRS, and a family history of glaucoma in POAG. As the non-IOP-related GRS increased, the risk of a family history of glaucoma increased. In contrast, a significant association (P = .0014; β = -0.14) was found between the IOP-related GRS, but not non-IOP-related GRS, and age at the diagnosis of glaucoma. As the IOP-related GRS increased, age at the diagnosis of glaucoma decreased.

CONCLUSION

The results indicate that non-IOP-related (optic nerve vulnerability) rather than IOP-related (IOP elevation) genetic variants may play an important role in the onset of POAG (family history of glaucoma) and that IOP-related rather than non-IOP-related genetic variants may play an important role in its progression (age at the diagnosis of glaucoma).

Transcriptomic stratification of late-onset Alzheimer's cases reveals novel genetic modifiers of disease pathology

Late-Onset Alzheimer's disease (LOAD) is a common, complex genetic disorder well-known for its heterogeneous pathology. The genetic heterogeneity underlying common, complex diseases poses a major challenge for targeted therapies and the identification of novel disease-associated variants. Case-control approaches are often limited to examining a specific outcome in a group of heterogenous patients with different clinical characteristics. Here, we developed a novel approach to define relevant transcriptomic endophenotypes and stratify decedents based on molecular profiles in three independent human LOAD cohorts. By integrating post-mortem brain gene co-expression data from 2114 human samples with LOAD, we developed a novel quantitative, composite phenotype that can better account for the heterogeneity in genetic architecture underlying the disease. We used iterative weighted gene co-expression network analysis (WGCNA) to reduce data dimensionality and to isolate gene sets that are highly co-expressed within disease subtypes and represent specific molecular pathways. We then performed single variant association testing using whole genome-sequencing data for the novel composite phenotype in order to identify genetic loci that contribute to disease heterogeneity. Distinct LOAD subtypes were identified for all three study cohorts (two in ROSMAP, three in Mayo Clinic, and two in Mount Sinai Brain Bank). Single variant association analysis identified a genome-wide significant variant in TMEM106B (p-value < 5×10-8, rs1990620G) in the ROSMAP cohort that confers protection from the inflammatory LOAD subtype. Taken together, our novel approach can be used to stratify LOAD into distinct molecular subtypes based on affected disease pathways.

Genome-wide association meta-analysis of corneal curvature identifies novel loci and shared genetic influences across axial length and refractive error

Corneal curvature, a highly heritable trait, is a key clinical endophenotype for myopia - a major cause of visual impairment and blindness in the world. Here we present a trans-ethnic meta-analysis of corneal curvature GWAS in 44,042 individuals of Caucasian and Asian with replication in 88,218 UK Biobank data. We identified 47 loci (of which 26 are novel), with population-specific signals as well as shared signals across ethnicities. Some identified variants showed precise scaling in corneal curvature and eye elongation (i.e. axial length) to maintain eyes in emmetropia (i.e. HDAC11/FBLN2 rs2630445, RBP3 rs11204213); others exhibited association with myopia with little pleiotropic effects on eye elongation. Implicated genes are involved in extracellular matrix organization, developmental process for body and eye, connective tissue cartilage and glycosylation protein activities. Our study provides insights into population-specific novel genes for corneal curvature, and their pleiotropic effect in regulating eye size or conferring susceptibility to myopia.

Multi-trait genome-wide association study identifies new loci associated with optic disc parameters

A new avenue of mining published genome-wide association studies includes the joint analysis of related traits. The power of this approach depends on the genetic correlation of traits, which reflects the number of pleiotropic loci, i.e. genetic loci influencing multiple traits. Here, we applied new meta-analyses of optic nerve head (ONH) related traits implicated in primary open-angle glaucoma (POAG); intraocular pressure and central corneal thickness using Haplotype reference consortium imputations. We performed a multi-trait analysis of ONH parameters cup area, disc area and vertical cup-disc ratio. We uncover new variants; rs11158547 in PPP1R36-PLEKHG3 and rs1028727 near SERPINE3 at genome-wide significance that replicate in independent Asian cohorts imputed to 1000 Genomes. At this point, validation of these variants in POAG cohorts is hampered by the high degree of heterogeneity. Our results show that multi-trait analysis is a valid approach to identify novel pleiotropic variants for ONH.

Myocilin Gene Gln368Ter Variant Penetrance and Association With Glaucoma in Population-Based and Registry-Based Studies

Importance

The p.Gln368Ter (rs74315329) risk allele in the myocilin gene (MYOC) was initially reported to have high penetrance in glaucoma registry-based studies, but much lower estimates were recently obtained from population-based studies. We investigated this disparity using data from Australia and the United Kingdom.

Objectives

To examine the penetrance and effect size of the MYOC p.Gln368Ter variant with glaucoma and ocular hypertension (OHT).

Design, Setting, and Participants

This cross-sectional study within the UK Biobank (UKBB) included participants of white British ancestry. Glaucoma cases were defined by International Classification of Diseases, Ninth Revision (ICD-9) and Tenth Revision (ICD-10) diagnoses and self-reported questionnaires. Carriers of the MYOC p.Gln368Ter variant were identified using genotype imputation from arrays. In contrast, 2 Australian registry-based studies, the Australian and New Zealand Registry of Advanced Glaucoma and the Glaucoma Inheritance Study in Tasmania, ascertained glaucoma cases referred by eye care clinicians, with historic control participants recruited from other Australian studies. Samples were either directly sequenced or had genotypes determined by imputation (for the Australian registry and historic control participants). Recruitment to the UKBB occurred between 2006 and 2010, and data analysis occurred from September 2017 to July 2018.

Main Outcomes and Measures

The penetrance and odds ratio (OR) were estimated for the MYOC p.Gln368Ter variants in participants with glaucoma and OHT.

Results

A total of 411 337 UKBB participants of white British ancestry (mean [SD] age, 56.6 [8.0] years) were included, plus 3071 Australian registry and 6750 historic control participants. In the UKBB, the minor allele frequency of the MYOC p.Gln368Ter variant was 1 in 786 individuals (0.13%). The odds ratio of p.Gln368Ter in patients with primary open-angle glaucoma (POAG) was 6.76 (95% CI, 4.05-11.29); glaucoma (POAG, self-reported glaucoma, and unspecified glaucoma), 4.40 (95% CI, 3.38-5.71); OHT, 3.56 (95% CI, 2.53-4.92); and OHT and glaucoma combined, 4.18 (95% CI, 3.05-5.67). The penetrance of the MYOC p.Gln368Ter variant was 7.6% in patients with glaucoma, 24.3% in patients with OHT, and 30.8% in patients with OHT and glaucoma combined. In the Australian registry studies, the odds of MYOC p.Gln368Ter variant were 12.16 (95% CI, 6.34-24.97) in patients with advanced glaucoma and 3.97 (95% CI, 1.55-9.75) in those with nonadvanced glaucoma; the penetrance of glaucoma was 56.1%, and penetrance in those considered to have glaucoma or be glaucoma suspects was 69.5%.

Conclusions and Relevance

The MYOC p.Gln368Ter variant confers a very high-risk effect size for advanced glaucoma; the risk is lower in nonadvanced glaucoma and OHT. In the general population sample, approximately 50% of MYOC p.Gln368Ter carriers 65 years and older had glaucoma or OHT, with higher prevalence in the Australian registry studies.

Association of Genetic Variants With Primary Open-Angle Glaucoma Among Individuals With African Ancestry

Importance

Primary open-angle glaucoma presents with increased prevalence and a higher degree of clinical severity in populations of African ancestry compared with European or Asian ancestry. Despite this, individuals of African ancestry remain understudied in genomic research for blinding disorders.

Objectives

To perform a genome-wide association study (GWAS) of African ancestry populations and evaluate potential mechanisms of pathogenesis for loci associated with primary open-angle glaucoma.

Design, Settings, and Participants

A 2-stage GWAS with a discovery data set of 2320 individuals with primary open-angle glaucoma and 2121 control individuals without primary open-angle glaucoma. The validation stage included an additional 6937 affected individuals and 14 917 unaffected individuals using multicenter clinic- and population-based participant recruitment approaches. Study participants were recruited from Ghana, Nigeria, South Africa, the United States, Tanzania, Britain, Cameroon, Saudi Arabia, Brazil, the Democratic Republic of the Congo, Morocco, Peru, and Mali from 2003 to 2018. Individuals with primary open-angle glaucoma had open iridocorneal angles and displayed glaucomatous optic neuropathy with visual field defects. Elevated intraocular pressure was not included in the case definition. Control individuals had no elevated intraocular pressure and no signs of glaucoma.

Exposures

Genetic variants associated with primary open-angle glaucoma.

Main Outcomes and Measures

Presence of primary open-angle glaucoma. Genome-wide significance was defined as P < 5 × 10-8 in the discovery stage and in the meta-analysis of combined discovery and validation data.

Results

A total of 2320 individuals with primary open-angle glaucoma (mean [interquartile range] age, 64.6 [56-74] years; 1055 [45.5%] women) and 2121 individuals without primary open-angle glaucoma (mean [interquartile range] age, 63.4 [55-71] years; 1025 [48.3%] women) were included in the discovery GWAS. The GWAS discovery meta-analysis demonstrated association of variants at amyloid-β A4 precursor protein-binding family B member 2 (APBB2; chromosome 4, rs59892895T>C) with primary open-angle glaucoma (odds ratio [OR], 1.32 [95% CI, 1.20-1.46]; P = 2 × 10-8). The association was validated in an analysis of an additional 6937 affected individuals and 14 917 unaffected individuals (OR, 1.15 [95% CI, 1.09-1.21]; P < .001). Each copy of the rs59892895*C risk allele was associated with increased risk of primary open-angle glaucoma when all data were included in a meta-analysis (OR, 1.19 [95% CI, 1.14-1.25]; P = 4 × 10-13). The rs59892895*C risk allele was present at appreciable frequency only in African ancestry populations. In contrast, the rs59892895*C risk allele had a frequency of less than 0.1% in individuals of European or Asian ancestry.

Conclusions and Relevance

In this genome-wide association study, variants at the APBB2 locus demonstrated differential association with primary open-angle glaucoma by ancestry. If validated in additional populations this finding may have implications for risk assessment and therapeutic strategies.

Clinical implications of recent advances in primary open-angle glaucoma genetics

Over the last decade, genetic studies, including genome-wide association studies (GWAS), have accelerated the discovery of genes and genomic regions contributing to primary open-angle glaucoma (POAG), a leading cause of irreversible vision loss. Here, we review the findings of genetic studies of POAG published in English prior to September 2019. In total, 74 genomic regions have been associated at a genome-wide level of significance with POAG susceptibility. Recent POAG GWAS provide not only insight into global and ethnic-specific genetic risk factors for POAG susceptibility across populations of diverse ancestry, but also important functional insights underlying biological mechanisms of glaucoma pathogenesis. In this review, we also summarize the genetic overlap between POAG, glaucoma endophenotypes, such as intraocular pressure and vertical cup-disc ratio (VCDR), and other eye disorders. We also discuss approaches recently developed to increase power for POAG locus discovery and to predict POAG risk. Finally, we discuss the recent development of POAG gene-based therapies and future strategies to treat glaucoma effectively. Understanding the genetic architecture of POAG is essential for an earlier diagnosis of this common eye disorder, predictive testing of at-risk patients, and design of gene-based targeted medical therapies none of which are currently available.

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 &lt; 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.

Mendelian genes in primary open angle glaucoma

Mutations in each of three genes, myocilin (MYOC), optineurin (OPTN), and TANK binding kinase 1 (TBK1), may cause primary open-angle glaucoma (POAG) that is inherited as a Mendelian trait. MYOC mutations cause 3-4% of POAG cases with IOP >21 mmHg, while mutations in OPTN, TBK1, and MYOC each cause ∼1% of POAG with IOP ≤21 mmHg, i.e. normal tension glaucoma. Identification of these disease-causing genes has provided insights into glaucoma pathogenesis. Mutations in MYOC cause a cascade of abnormalities in the trabecular meshwork including intracellular retention of MYOC protein, decreased aqueous outflow, higher intraocular pressure, and glaucoma. Investigation of MYOC mutations demonstrated that abnormal retention of intracellular MYOC and stimulation of endoplasmic reticular (ER) stress may be important steps in the development of MYOC-associated glaucoma. Mutations in OPTN and TBK1 cause a dysregulation of autophagy which may directly cause retinal ganglion cell damage and normal tension glaucoma. Discovery of these Mendelian causes of glaucoma has also provided a new set of potential therapeutic targets that may ultimately lead to novel, gene-directed glaucoma treatments.