Search

The underrepresentation of non-Europeans in human genetic studies so far has limited the diversity of individuals in genomic datasets and led to reduced medical relevance for a large proportion of the world's population. Population-specific reference genome datasets as well as genome-wide association studies in diverse populations are needed to address this issue. Here we describe the pilot phase of the GenomeAsia 100K Project. This includes a whole-genome sequencing reference dataset from 1,739 individuals of 219 population groups and 64 countries across Asia. We catalogue genetic variation, population structure, disease associations and founder effects. We also explore the use of this dataset in imputation, to facilitate genetic studies in populations across Asia and worldwide.

ImportancePrimary 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.ObjectivesTo 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 participantsA 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.ExposuresGenetic variants associated with primary open-angle glaucoma.Main outcomes and measuresPresence of primary open-angle glaucoma. Genome-wide significance was defined as P 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

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 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

Primary open angle glaucoma (POAG) is a complex disease with a major genetic contribution. Its prevalence varies greatly among ethnic groups, and is up to five times more frequent in black African populations compared to Europeans. So far, worldwide efforts to elucidate the genetic complexity of POAG in African populations has been limited. We conducted a genome-wide association study in 1113 POAG cases and 1826 controls from Tanzanian, South African and African American study samples. Apart from confirming evidence of association at TXNRD2 (rs16984299; OR[T] 1.20; P = 0.003), we found that a genetic risk score combining the effects of the 15 previously reported POAG loci was significantly associated with POAG in our samples (OR 1.56; 95% CI 1.26-1.93; P = 4.79 × 10-5). By genome-wide association testing we identified a novel candidate locus, rs141186647, harboring EXOC4 (OR[A] 0.48; P = 3.75 × 10-8), a gene transcribing a component of the exocyst complex involved in vesicle transport. The low frequency and high degree of genetic heterogeneity at this region hampered validation of this finding in predominantly West-African replication sets. Our results suggest that established genetic risk factors play a role in African POAG, however, they do not explain the higher disease load. The high heterogeneity within Africans remains a challenge to identify the genetic commonalities for POAG in this ethnicity, and demands studies of extremely large size.

Purpose:Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk. Methods:We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]). Results:In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations. Conclusions:Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

Purpose:Sex hormones may be associated with primary open-angle glaucoma (POAG), although the mechanisms are unclear. We previously observed that gene variants involved with estrogen metabolism were collectively associated with POAG in women but not men; here we assessed gene variants related to testosterone metabolism collectively and POAG risk. Methods:We used two datasets: one from the United States (3853 cases and 33,480 controls) and another from Australia (1155 cases and 1992 controls). Both datasets contained densely called genotypes imputed to the 1000 Genomes reference panel. We used pathway- and gene-based approaches with Pathway Analysis by Randomization Incorporating Structure (PARIS) software to assess the overall association between a panel of single nucleotide polymorphisms (SNPs) in testosterone metabolism genes and POAG. In sex-stratified analyses, we evaluated POAG overall and POAG subtypes defined by maximum IOP (high-tension [HTG] or normal tension glaucoma [NTG]). Results:In the US dataset, the SNP panel was not associated with POAG (permuted P = 0.77), although there was an association in the Australian sample (permuted P = 0.018). In both datasets, the SNP panel was associated with POAG in men (permuted P ≤ 0.033) and not women (permuted P ≥ 0.42), but in gene-based analyses, there was no consistency on the main genes responsible for these findings. In both datasets, the testosterone pathway association with HTG was significant (permuted P ≤ 0.011), but again, gene-based analyses showed no consistent driver gene associations. Conclusions:Collectively, testosterone metabolism pathway SNPs were consistently associated with the high-tension subtype of POAG in two datasets.

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.

Primary open-angle glaucoma (POAG) is the most common chronic optic neuropathy worldwide. Epidemiological studies show a robust positive relation between intraocular pressure (IOP) and POAG and modest positive association between IOP and blood pressure (BP), while the relation between BP and POAG is controversial. The International Glaucoma Genetics Consortium (n=27 558), the International Consortium on Blood Pressure (n=69 395), and the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (n=37 333), represent genome-wide data sets for IOP, BP traits and POAG, respectively. We formed genome-wide significant variant panels for IOP and diastolic BP and found a strong relation with POAG (odds ratio and 95% confidence interval: 1.18 (1.14-1.21), P=1.8 × 10-27) for the former trait but no association for the latter (P=0.93). Next, we used linkage disequilibrium (LD) score regression, to provide genome-wide estimates of correlation between traits without the need for additional phenotyping. We also compared our genome-wide estimate of heritability between IOP and BP to an estimate based solely on direct measures of these traits in the Erasmus Rucphen Family (ERF; n=2519) study using Sequential Oligogenic Linkage Analysis Routines (SOLAR). LD score regression revealed high genetic correlation between IOP and POAG (48.5%, P=2.1 × 10-5); however, genetic correlation between IOP and diastolic BP (P=0.86) and between diastolic BP and POAG (P=0.42) were negligible. Using SOLAR in the ERF study, we confirmed the minimal heritability between IOP and diastolic BP (P=0.63). Overall, IOP shares genetic basis with POAG, whereas BP has limited shared genetic correlation with IOP or POAG.

ObjectiveSeveral attributes of female reproductive history, including age at natural menopause (ANM), have been related to primary open-angle glaucoma (POAG). We assembled 18 previously reported common genetic variants that predict ANM to determine their association with ANM or POAG.MethodsUsing data from the Nurses' Health Study (7,143 women), we validated the ANM weighted genetic risk score in relation to self-reported ANM. Subsequently, to assess the relation with POAG, we used data from 2,160 female POAG cases and 29,110 controls in the National Eye Institute Glaucoma Human Genetics Collaboration Heritable Overall Operational Database (NEIGHBORHOOD), which consists of 8 datasets with imputed genotypes to 5.6+ million markers. Associations with POAG were assessed in each dataset, and site-specific results were meta-analyzed using the inverse weighted variance method.ResultsThe genetic risk score was associated with self-reported ANM (P = 2.2 × 10) and predicted 4.8% of the variance in ANM. The ANM genetic risk score was not associated with POAG (Odds Ratio (OR) = 1.002; 95% Confidence Interval (CI): 0.998, 1.007; P = 0.28). No single genetic variant in the panel achieved nominal association with POAG (P ≥0.20). Compared to the middle 80 percent, there was also no association with the lowest 10 percentile or highest 90 percentile of genetic risk score with POAG (OR = 0.75; 95% CI: 0.47, 1.21; P = 0.23 and OR = 1.10; 95% CI: 0.72, 1.69; P = 0.65, respectively).ConclusionsA genetic risk score predicting 4.8% of ANM variation was not related to POAG; thus, genetic determinants of ANM are unlikely to explain the previously reported association between the two phenotypes.

PurposeRecent studies indicate that mitochondrial proteins may contribute to the pathogenesis of primary open-angle glaucoma (POAG). In this study, we examined the association between POAG and common variations in gene-encoding mitochondrial proteins.MethodsWe examined genetic data from 3430 POAG cases and 3108 controls derived from the combination of the GLAUGEN and NEIGHBOR studies. We constructed biological-system coherent mitochondrial nuclear-encoded protein gene-sets by intersecting the MitoCarta database with the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. We examined the mitochondrial gene-sets for association with POAG and with normal-tension glaucoma (NTG) and high-tension glaucoma (HTG) subsets using Pathway Analysis by Randomization Incorporating Structure.ResultsWe identified 22 KEGG pathways with significant mitochondrial protein-encoding gene enrichment, belonging to six general biological classes. Among the pathway classes, mitochondrial lipid metabolism was associated with POAG overall (P = 0.013) and with NTG (P = 0.0006), and mitochondrial carbohydrate metabolism was associated with NTG (P = 0.030). Examining the individual KEGG pathway mitochondrial gene-sets, fatty acid elongation and synthesis and degradation of ketone bodies, both lipid metabolism pathways, were significantly associated with POAG (P = 0.005 and P = 0.002, respectively) and NTG (P = 0.0004 and P < 0.0001, respectively). Butanoate metabolism, a carbohydrate metabolism pathway, was significantly associated with POAG (P = 0.004), NTG (P = 0.001), and HTG (P = 0.010).ConclusionsWe present an effective approach for assessing the contributions of mitochondrial genetic variation to open-angle glaucoma. Our findings support a role for mitochondria in POAG pathogenesis and specifically point to lipid and carbohydrate metabolism pathways as being important.

PurposeNoncoding microRNAs (miRNAs) have been implicated in the pathogenesis of glaucoma. We aimed to identify common variants in miRNA coding genes (MIR) associated with primary open-angle glaucoma (POAG).MethodsUsing the NEIGHBORHOOD data set (3853 cases/33,480 controls with European ancestry), we first assessed the relation between 85 variants in 76 MIR genes and overall POAG. Subtype-specific analyses were performed in high-tension glaucoma (HTG) and normal-tension glaucoma subsets. Second, we examined the expression of miR-182, which was associated with POAG, in postmortem human ocular tissues (ciliary body, cornea, retina, and trabecular meshwork [TM]), using miRNA sequencing (miRNA-Seq) and droplet digital PCR (ddPCR). Third, miR-182 expression was also examined in human aqueous humor (AH) by using miRNA-Seq. Fourth, exosomes secreted from primary human TM cells were examined for miR-182 expression by using miRNA-Seq. Fifth, using ddPCR we compared miR-182 expression in AH between five HTG cases and five controls.ResultsOnly rs76481776 in MIR182 gene was associated with POAG after adjustment for multiple comparisons (odds ratio [OR] = 1.23, 95% confidence interval [CI]: 1.11-1.42, P = 0.0002). Subtype analysis indicated that the association was primarily in the HTG subset (OR = 1.26, 95% CI: 1.08-1.47, P = 0.004). The risk allele T has been associated with elevated miR-182 expression in vitro. Data from ddPCR and miRNA-Seq confirmed miR-182 expression in all examined ocular tissues and TM-derived exosomes. Interestingly, miR-182 expression in AH was 2-fold higher in HTG patients than nonglaucoma controls (P = 0.03) without controlling for medication treatment.ConclusionsOur integrative study is the first to associate rs76481776 with POAG via elevated miR-182 expression.

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. To identify new susceptibility loci, we performed meta-analysis on genome-wide association study (GWAS) results from eight independent studies from the United States (3,853 cases and 33,480 controls) and investigated the most significantly associated SNPs in two Australian studies (1,252 cases and 2,592 controls), three European studies (875 cases and 4,107 controls) and a Singaporean Chinese study (1,037 cases and 2,543 controls). A meta-analysis of the top SNPs identified three new associated loci: rs35934224[T] in TXNRD2 (odds ratio (OR) = 0.78, P = 4.05 × 10(-11)) encoding a mitochondrial protein required for redox homeostasis; rs7137828[T] in ATXN2 (OR = 1.17, P = 8.73 × 10(-10)); and rs2745572[A] upstream of FOXC1 (OR = 1.17, P = 1.76 × 10(-10)). Using RT-PCR and immunohistochemistry, we show TXNRD2 and ATXN2 expression in retinal ganglion cells and the optic nerve head. These results identify new pathways underlying POAG susceptibility and suggest new targets for preventative therapies.

Exfoliation syndrome (XFS) is the most common recognizable cause of open-angle glaucoma worldwide. To better understand the etiology of XFS, we conducted a genome-wide association study (GWAS) of 1,484 cases and 1,188 controls from Japan and followed up the most significant findings in a further 6,901 cases and 20,727 controls from 17 countries across 6 continents. We discovered a genome-wide significant association between a new locus (CACNA1A rs4926244) and increased susceptibility to XFS (odds ratio (OR) = 1.16, P = 3.36 × 10(-11)). Although we also confirmed overwhelming association at the LOXL1 locus, the key SNP marker (LOXL1 rs4886776) demonstrated allelic reversal depending on the ancestry group (Japanese: OR(A allele) = 9.87, P = 2.13 × 10(-217); non-Japanese: OR(A allele) = 0.49, P = 2.35 × 10(-31)). Our findings represent the first genetic locus outside of LOXL1 surpassing genome-wide significance for XFS and provide insight into the biology and pathogenesis of the disease.

PurposeWe examined the role of DNA copy number variants (CNVs) of known glaucoma genes in relation to primary open angle glaucoma (POAG).MethodsOur study included DNA samples from two studies (NEIGHBOR and GLAUGEN). All the samples were genotyped with the Illumina Human660W_Quad_v1 BeadChip. After removing non-blood-derived and amplified DNA samples, we applied quality control steps based on the mean Log R Ratio and the mean B allele frequency. Subsequently, data from 3057 DNA samples (1599 cases and 1458 controls) were analyzed with PennCNV software. We defined CNVs as those ≥5 kilobases (kb) in size and interrogated by ≥5 consecutive probes. We further limited our investigation to CNVs in known POAG-related genes, including CDKN2B-AS1, TMCO1, SIX1/SIX6, CAV1/CAV2, the LRP12-ZFPM2 region, GAS7, ATOH7, FNDC3B, CYP1B1, MYOC, OPTN, WDR36, SRBD1, TBK1, and GALC.ResultsGenomic duplications of CDKN2B-AS1 and TMCO1 were each found in a single case. Two cases carried duplications in the GAS7 region. Genomic deletions of SIX6 and ATOH7 were each identified in one case. One case carried a TBK1 deletion and another case carried a TBK1 duplication. No controls had duplications or deletions in these six genes. A single control had a duplication in the MYOC region. Deletions of GALC were observed in five cases and two controls.ConclusionsThe CNV analysis of a large set of cases and controls revealed the presence of rare CNVs in known POAG susceptibility genes. Our data suggest that these rare CNVs may contribute to POAG pathogenesis and merit functional evaluation.

Primary open-angle glaucoma (POAG) is a leading cause of blindness worldwide. Using genome-wide association single-nucleotide polymorphism data from the Glaucoma Genes and Environment study and National Eye Institute Glaucoma Human Genetics Collaboration comprising 3,108 cases and 3,430 controls, we assessed biologic pathways as annotated in the KEGG database for association with risk of POAG. After correction for genic overlap among pathways, we found 4 pathways, butanoate metabolism (hsa00650), hematopoietic cell lineage (hsa04640), lysine degradation (hsa00310) and basal transcription factors (hsa03022) related to POAG with permuted p

PurposeThe CAV1/CAV2 (caveolin 1 and caveolin 2) genomic region previously was associated with primary open-angle glaucoma (POAG), although replication among independent studies has been variable. The aim of this study was to assess the association between CAV1/CAV2 single nucleotide polymorphisms (SNPs) and POAG in a large case-control dataset and to explore associations by gender and pattern of visual field (VF) loss further.DesignCase-control study.ParticipantsWe analyzed 2 large POAG data sets: the Glaucoma Genes and Environment (GLAUGEN) study (976 cases, 1140 controls) and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium (2132 cases, 2290 controls).MethodsWe studied the association between 70 SNPs located within the CAV1/CAV2 genomic region in the GLAUGEN and NEIGHBOR studies, both genotyped on the Illumina Human 660WQuadv1C BeadChip array and imputed with the Markov Chain Haplotyping algorithm using the HapMap 3 reference panel. We used logistic regression models of POAG in the overall population and separated by gender, as well as by POAG subtypes defined by type of VF defect (peripheral or paracentral). Results from GLAUGEN and NEIGHBOR were meta-analyzed, and a Bonferroni-corrected significance level of 7.7 × 10(-4) was used to account for multiple comparisons.Main outcome measuresOverall POAG, overall POAG by gender, and POAG subtypes defined by pattern of early VF loss.ResultsWe found significant associations between 10 CAV1/CAV2 SNPs and POAG (top SNP, rs4236601; pooled P = 2.61 × 10(-7)). Of these, 9 were significant only in women (top SNP, rs4236601; pooled P = 1.59 × 10(-5)). Five of the 10 CAV1/CAV2 SNPs were associated with POAG with early paracentral VF (top SNP, rs17588172; pooled P = 1.07 × 10(-4)), and none of the 10 were associated with POAG with peripheral VF loss only or POAG among men.ConclusionsCAV1/CAV2 SNPs were associated significantly with POAG overall, particularly among women. Furthermore, we found an association between CAV1/CAV2 SNPs and POAG with paracentral VF defects. These data support a role for caveolin 1, caveolin 2, or both in POAG and suggest that the caveolins particularly may affect POAG pathogenesis in women and in patients with early paracentral VF defects.

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.

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.

PurposeCirculating estrogen levels are relevant in glaucoma phenotypic traits. We assessed the association between an estrogen metabolism single nucleotide polymorphism (SNP) panel in relation to primary open angle glaucoma (POAG), accounting for gender.MethodsWe included 3,108 POAG cases and 3,430 controls of both genders from the Glaucoma Genes and Environment (GLAUGEN) study and the National Eye Institute Glaucoma Human Genetics Collaboration (NEIGHBOR) consortium genotyped on the Illumina 660W-Quad platform. We assessed the relation between the SNP panels representative of estrogen metabolism and POAG using pathway- and gene-based approaches with the Pathway Analysis by Randomization Incorporating Structure (PARIS) software. PARIS executes a permutation algorithm to assess statistical significance relative to the pathways and genes of comparable genetic architecture. These analyses were performed using the meta-analyzed results from the GLAUGEN and NEIGHBOR data sets. We evaluated POAG overall as well as two subtypes of POAG defined as intraocular pressure (IOP) ≥22 mmHg (high-pressure glaucoma [HPG]) or IOP 0.99). Among women, gene-based analyses revealed that the catechol-O-methyltransferase gene showed strong associations with HTG (permuted gene p≤0.001) and NPG (permuted gene p=0.01).ConclusionsThe estrogen SNP pathway was associated with POAG among women.

Optic nerve degeneration caused by glaucoma is a leading cause of blindness worldwide. Patients affected by the normal-pressure form of glaucoma are more likely to harbor risk alleles for glaucoma-related optic nerve disease. We have performed a meta-analysis of two independent genome-wide association studies for primary open angle glaucoma (POAG) followed by a normal-pressure glaucoma (NPG, defined by intraocular pressure (IOP) less than 22 mmHg) subgroup analysis. The single-nucleotide polymorphisms that showed the most significant associations were tested for association with a second form of glaucoma, exfoliation-syndrome glaucoma. The overall meta-analysis of the GLAUGEN and NEIGHBOR dataset results (3,146 cases and 3,487 controls) identified significant associations between two loci and POAG: the CDKN2BAS region on 9p21 (rs2157719 [G], OR = 0.69 [95%CI 0.63-0.75], p = 1.86×10⁻¹⁸), and the SIX1/SIX6 region on chromosome 14q23 (rs10483727 [A], OR = 1.32 [95%CI 1.21-1.43], p = 3.87×10⁻¹¹). In sub-group analysis two loci were significantly associated with NPG: 9p21 containing the CDKN2BAS gene (rs2157719 [G], OR = 0.58 [95% CI 0.50-0.67], p = 1.17×10⁻¹²) and a probable regulatory region on 8q22 (rs284489 [G], OR = 0.62 [95% CI 0.53-0.72], p = 8.88×10⁻¹⁰). Both NPG loci were also nominally associated with a second type of glaucoma, exfoliation syndrome glaucoma (rs2157719 [G], OR = 0.59 [95% CI 0.41-0.87], p = 0.004 and rs284489 [G], OR = 0.76 [95% CI 0.54-1.06], p = 0.021), suggesting that these loci might contribute more generally to optic nerve degeneration in glaucoma. Because both loci influence transforming growth factor beta (TGF-beta) signaling, we performed a genomic pathway analysis that showed an association between the TGF-beta pathway and NPG (permuted p = 0.009). These results suggest that neuro-protective therapies targeting TGF-beta signaling could be effective for multiple forms of glaucoma.

Although glaucoma is a disease modulated by eye pressure, the mechanisms of pressure sensing in the eye are not well understood. Here, we investigated associations between mechanosensitive ion channel gene variants and primary open-angle glaucoma (POAG). Common (minor allele frequency > 5%) single nucleotide polymorphisms located within the genomic regions of 20 mechanosensitive ion channel genes in the K2P, TMEM63, PIEZO and TRP channel families were assessed using genotype data from the NEIGHBORHOOD consortium of 3853 cases and 33,480 controls. Rare (minor allele frequency < 1%) coding variants were assessed using exome array genotyping data for 2606 cases and 2606 controls. Association with POAG was analyzed using logistic regression adjusting for age and sex. Two rare PIEZO1 coding variants with protective effects were identified in the NEIGHBOR dataset: R1527H, (OR 0.17, P = 0.0018) and a variant that alters a canonical splice donor site, g.16-88737727-C-G Hg38 (OR 0.38, P = 0.02). Both variants showed similar effects in the UK Biobank and the R1527H also in the FinnGen database. Several common variants also reached study-specific thresholds for association in the NEIGHBORHOOD dataset. These results identify novel variants in several mechanosensitive channel genes that show associations with POAG, suggesting that these channels may be potential therapeutic targets. [ABSTRACT FROM AUTHOR]