Person:

Wiggs, Janey

Background: Loss of vision in glaucoma is due to apoptotic retinal ganglion cell loss. While p53 modulates apoptosis, gene association studies between p53 variants and glaucoma have been inconsistent. In this study we evaluate the association between a p53 variant functionally known to influence apoptosis (codon 72 Pro/Arg) and the subset of primary open angle glaucoma (POAG) patients with early loss of central visual field. Methods: Genotypes for the p53 codon 72 polymorphism (Pro/Arg) were obtained for 264 POAG patients and 400 controls from the U.S. and in replication studies for 308 POAG patients and 178 controls from Australia (GIST). The glaucoma patients were divided into two groups according to location of initial visual field defect (either paracentral or peripheral). All cases and controls were Caucasian with European ancestry. Results: The p53-PRO/PRO genotype was more frequent in the U.S. POAG patients with early visual field defects in the paracentral regions compared with those in the peripheral regions or control group (p = 2.7×(10^{−5})). We replicated this finding in the GIST cohort (p = 7.3×(10^{−3}), and in the pooled sample (p = 6.6×(10^{−7})) and in a meta-analysis of both the US and GIST datasets (1.3×(10^{−6}), OR 2.17 (1.58–2.98 for the PRO allele). Conclusions: These results suggest that the p53 codon 72 PRO/PRO genotype is potentially associated with early paracentral visual field defects in primary open-angle glaucoma patients.

Primary open angle glaucoma (POAG) is a leading cause of blindness worldwide. The molecular signaling involved in the pathogenesis of POAG remains unknown. Here, we report that mice lacking the (α_1) subunit of the nitric oxide receptor soluble guanylate cyclase represent a novel and translatable animal model of POAG, characterized by thinning of the retinal nerve fiber layer and loss of optic nerve axons in the context of an open iridocorneal angle. The optic neuropathy associated with soluble guanylate cyclase (α_1)–deficiency was accompanied by modestly increased intraocular pressure and retinal vascular dysfunction. Moreover, data from a candidate gene association study suggests that a variant in the locus containing the genes encoding for the (α_1) and (β_1) subunits of soluble guanylate cyclase is associated with POAG in patients presenting with initial paracentral vision loss, a disease subtype thought to be associated with vascular dysregulation. These findings provide new insights into the pathogenesis and genetics of POAG and suggest new therapeutic strategies for POAG.

Glaucoma is a leading cause of blindness worldwide. Primary open-angle glaucoma (POAG) is the most common subtype and is a complex trait with multigenic inheritance. Genome-wide association studies have previously identified a significant association between POAG and the SIX6 locus (rs10483727, odds ratio (OR) = 1.32, p = 3.87×10−11). SIX6 plays a role in ocular development and has been associated with the morphology of the optic nerve. We sequenced the SIX6 coding and regulatory regions in 262 POAG cases and 256 controls and identified six nonsynonymous coding variants, including five rare and one common variant, Asn141His (rs33912345), which was associated significantly with POAG (OR = 1.27, p = 4.2×10−10) in the NEIGHBOR/GLAUGEN datasets. These variants were tested in an in vivo Danio rerio (zebrafish) complementation assay to evaluate ocular metrics such as eye size and optic nerve structure. Five variants, found primarily in POAG cases, were hypomorphic or null, while the sixth variant, found only in controls, was benign. One variant in the SIX6 enhancer increased expression of SIX6 and disrupted its regulation. Finally, to our knowledge for the first time, we have identified a clinical feature in POAG patients that appears to be dependent upon SIX6 genotype: patients who are homozygous for the SIX6 risk allele (His141) have a statistically thinner retinal nerve fiber layer than patients homozygous for the SIX6 non-risk allele (Asn141). Our results, in combination with previous SIX6 work, lead us to hypothesize that SIX6 risk variants disrupt the development of the neural retina, leading to a reduced number of retinal ganglion cells, thereby increasing the risk of glaucoma-associated vision loss.

Primary open angle glaucoma (POAG) is a complex disease and is one of the major leading causes of blindness worldwide. Genome-wide association studies have successfully identified several common variants associated with glaucoma; however, most of these variants only explain a small proportion of the genetic risk. Apart from the standard approach to identify main effects of variants across the genome, it is believed that gene-gene interactions can help elucidate part of the missing heritability by allowing for the test of interactions between genetic variants to mimic the complex nature of biology. To explain the etiology of glaucoma, we first performed a genome-wide association study (GWAS) on glaucoma case-control samples obtained from electronic medical records (EMR) to establish the utility of EMR data in detecting non-spurious and relevant associations; this analysis was aimed at confirming already known associations with glaucoma and validating the EMR derived glaucoma phenotype. Our findings from GWAS suggest consistent evidence of several known associations in POAG. We then performed an interaction analysis for variants found to be marginally associated with glaucoma (SNPs with main effect p-value <0.01) and observed interesting findings in the electronic MEdical Records and GEnomics Network (eMERGE) network dataset. Genes from the top epistatic interactions from eMERGE data (Likelihood Ratio Test i.e. LRT p-value <1e-05) were then tested for replication in the NEIGHBOR consortium dataset. To replicate our findings, we performed a gene-based SNP-SNP interaction analysis in NEIGHBOR and observed significant gene-gene interactions (p-value <0.001) among the top 17 gene-gene models identified in the discovery phase. Variants from gene-gene interaction analysis that we found to be associated with POAG explain 3.5% of additional genetic variance in eMERGE dataset above what is explained by the SNPs in genes that are replicated from previous GWAS studies (which was only 2.1% variance explained in eMERGE dataset); in the NEIGHBOR dataset, adding replicated SNPs from gene-gene interaction analysis explain 3.4% of total variance whereas GWAS SNPs alone explain only 2.8% of variance. Exploring gene-gene interactions may provide additional insights into many complex traits when explored in properly designed and powered association studies.

The Nurses’ Health Study (NHS), Nurses’ Health Study II (NHSII), Health Professionals Follow Up Study (HPFS) and the Physicians Health Study (PHS) have collected detailed longitudinal data on multiple exposures and traits for approximately 310,000 study participants over the last 35 years. Over 160,000 study participants across the cohorts have donated a DNA sample and to date, 20,691 subjects have been genotyped as part of genome-wide association studies (GWAS) of twelve primary outcomes. However, these studies utilized six different GWAS arrays making it difficult to conduct analyses of secondary phenotypes or share controls across studies. To allow for secondary analyses of these data, we have created three new datasets merged by platform family and performed imputation using a common reference panel, the 1,000 Genomes Phase I release. Here, we describe the methodology behind the data merging and imputation and present imputation quality statistics and association results from two GWAS of secondary phenotypes (body mass index (BMI) and venous thromboembolism (VTE)). We observed the strongest BMI association for the FTO SNP rs55872725 (β = 0.45, p = 3.48x10-22), and using a significance level of p = 0.05, we replicated 19 out of 32 known BMI SNPs. For VTE, we observed the strongest association for the rs2040445 SNP (OR = 2.17, 95% CI: 1.79–2.63, p = 2.70x10-15), located downstream of F5 and also observed significant associations for the known ABO and F11 regions. This pooled resource can be used to maximize power in GWAS of phenotypes collected across the cohorts and for studying gene-environment interactions as well as rare phenotypes and genotypes.

Shorter axial length observed in patients with primary angle closure glaucoma (PACG) might be due to altered matrix metalloproteinase-9 (MMP9) activity resulting in ECM remodeling during eye growth and development. This study aimed to evaluate common variants in MMP9 for association with PACG. Six tag SNPs of MMP9 were genotyped in a Chinese sample of 1,030 cases, including 572 PACG and 458 primary angle closure (PAC), and 499 controls. None of 6 SNPs were significantly associated with overall PAC/PACG (P > 0.07) or with PAC/PACG subgroups (Pc > 0.18). Meta-analysis of two non-Chinese studies revealed significant association between rs17576 and PACG (ORs = 0.56, P < 0.0001); however, meta-analysis of our dataset with 4 Chinese datasets did not replicate this association (ORs = 1.23, P = 0.29). Prior significant association for rs3918249 in one Caucasian study (OR = 0.63, P = 0.006) was not replicated in meta-analysis of 3 Chinese studies including this study (ORs = 0.91, P = 0.13). Significant heterogeneity between non-Chinese and Chinese datasets precluded overall meta-analysis for rs17576 and rs3918249 (Q = 0.001 and 0.04 respectively). rs17577 was nominally associated with PACG in one Caucasian study (OR = 1.71, P = 0.02), but not in 3 Chinese studies including our study (ORs = 1.20, P = 0.07). Overall meta-analysis revealed nominal association for rs17577 and PAC/PACG (ORs = 1.26, Pc = 0.05). Meta-analysis did not show significant association between the other SNPs and PAC/PACG (P > 0.47). The largest association study to date did not find significant association between MMP9 and PAC/PACG in Chinese; meta-analysis with other Chinese datasets did not produce significant association. In most instances combination with non-Chinese datasets was not possible except for one variant showing nominally significant association. More work is needed to define the role of MMP9 variants in PACG.

Central corneal thickness (CCT) is one of the most heritable ocular traits and it is also a phenotypic risk factor for primary open angle glaucoma (POAG). The present study uses the BXD Recombinant Inbred (RI) strains to identify novel quantitative trait loci (QTLs) modulating CCT in the mouse with the potential of identifying a molecular link between CCT and risk of developing POAG. The BXD RI strain set was used to define mammalian genomic loci modulating CCT, with a total of 818 corneas measured from 61 BXD RI strains (between 60–100 days of age). The mice were anesthetized and the eyes were positioned in front of the lens of the Phoenix Micron IV Image-Guided OCT system or the Bioptigen OCT system. CCT data for each strain was averaged and used to QTLs modulating this phenotype using the bioinformatics tools on GeneNetwork (www.genenetwork.org). The candidate genes and genomic loci identified in the mouse were then directly compared with the summary data from a human POAG genome wide association study (NEIGHBORHOOD) to determine if any genomic elements modulating mouse CCT are also risk factors for POAG.This analysis revealed one significant QTL on Chr 13 and a suggestive QTL on Chr 7. The significant locus on Chr 13 (13 to 19 Mb) was examined further to define candidate genes modulating this eye phenotype. For the Chr 13 QTL in the mouse, only one gene in the region (Pou6f2) contained nonsynonymous SNPs. Of these five nonsynonymous SNPs in Pou6f2, two resulted in changes in the amino acid proline which could result in altered secondary structure affecting protein function. The 7 Mb region under the mouse Chr 13 peak distributes over 2 chromosomes in the human: Chr 1 and Chr 7. These genomic loci were examined in the NEIGHBORHOOD database to determine if they are potential risk factors for human glaucoma identified using meta-data from human GWAS. The top 50 hits all resided within one gene (POU6F2), with the highest significance level of p = 10−6 for SNP rs76319873. POU6F2 is found in retinal ganglion cells and in corneal limbal stem cells. To test the effect of POU6F2 on CCT we examined the corneas of a Pou6f2-null mice and the corneas were thinner than those of wild-type littermates. In addition, these POU6F2 RGCs die early in the DBA/2J model of glaucoma than most RGCs. Using a mouse genetic reference panel, we identified a transcription factor, Pou6f2, that modulates CCT in the mouse. POU6F2 is also found in a subset of retinal ganglion cells and these RGCs are sensitive to injury.