Puya Gharahkhani

Genome-wide association analysis identifies TXNRD2, ATXN2 and FOXC1 as susceptibility loci for primary open-angle glaucoma

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.

Common variants near ABCA1 , AFAP1 and GMDS confer risk of primary open-angle glaucoma

Primary open-angle glaucoma (POAG) is a major cause of irreversible blindness worldwide. We performed a genome-wide association study in an Australian discovery cohort comprising 1,155 cases with advanced POAG and 1,992 controls. We investigated the association of the top SNPs from the discovery stage in two Australian replication cohorts (932 cases and 6,862 controls total) and two US replication cohorts (2,616 cases and 2,634 controls total). Meta-analysis of all cohorts identified three loci newly associated with development of POAG. These loci are located upstream of ABCA1 (rs2472493[G], odds ratio (OR) = 1.31, P = 2.1 × 10−19), within AFAP1 (rs4619890[G], OR = 1.20, P = 7.0 × 10−10) and within GMDS (rs11969985[G], OR = 1.31, P = 7.7 × 10−10). Using RT-PCR and immunolabeling, we show that these genes are expressed within human retina, optic nerve and trabecular meshwork and that ABCA1 and AFAP1 are also expressed in retinal ganglion cells.

ARHGEF12 influences the risk of glaucoma by increasing intraocular pressure

Primary open-angle glaucoma (POAG) is a blinding disease. Two important risk factors for this disease are a positive family history and elevated intraocular pressure (IOP), which is also highly heritable. Genes found to date associated with IOP and POAG are ABCA1, CAV1/CAV2, GAS7 and TMCO1. However, these genes explain only a small part of the heritability of IOP and POAG. We performed a genome-wide association study of IOP in the population-based Rotterdam Study I and Rotterdam Study II using single nucleotide polymorphisms (SNPs) imputed to 1000 Genomes. In this discovery cohort (n = 8105), we identified a new locus associated with IOP. The most significantly associated SNP was rs58073046 (β = 0.44, P-value = 1.87 × 10(-8), minor allele frequency = 0.12), within the gene ARHGEF12. Independent replication in five population-based studies (n = 7471) resulted in an effect size in the same direction that was significantly associated (β = 0.16, P-value = 0.04). The SNP was also significantly associated with POAG in two independent case-control studies [n = 1225 cases and n = 4117 controls; odds ratio (OR) = 1.53, P-value = 1.99 × 10(-8)], especially with high-tension glaucoma (OR = 1.66, P-value = 2.81 × 10(-9); for normal-tension glaucoma OR = 1.29, P-value = 4.23 × 10(-2)). ARHGEF12 plays an important role in the RhoA/RhoA kinase pathway, which has been implicated in IOP regulation. Furthermore, it binds to ABCA1 and links the ABCA1, CAV1/CAV2 and GAS7 pathway to Mendelian POAG genes (MYOC, OPTN, WDR36). In conclusion, this study identified a novel association between IOP and ARHGEF12.

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

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

Genome-wide association studies in oesophageal adenocarcinoma and Barrett's oesophagus: a large-scale meta-analysis

Summary Background Oesophageal adenocarcinoma represents one of the fastest rising cancers in high-income countries. Barretts oesophagus is the premalignant precursor of oesophageal adenocarcinoma. However, only a few patients with Barretts oesophagus develop adenocarcinoma, which complicates clinical management in the absence of valid predictors. Within an international consortium investigating the genetics of Barretts oesophagus and oesophageal adenocarcinoma, we aimed to identify novel genetic risk variants for the development of Barretts oesophagus and oesophageal adenocarcinoma. Methods We did a meta-analysis of all genome-wide association studies of Barretts oesophagus and oesophageal adenocarcinoma available in PubMed up to Feb 29, 2016; all patients were of European ancestry and disease was confirmed histopathologically. All participants were from four separate studies within Europe, North America, and Australia and were genotyped on high-density single nucleotide polymorphism (SNP) arrays. Meta-analysis was done with a fixed-effects inverse variance-weighting approach and with a standard genome-wide significance threshold (p<5 × 10−8). We also did an association analysis after reweighting of loci with an approach that investigates annotation enrichment among genome-wide significant loci. Furthermore, the entire dataset was analysed with bioinformatics approaches—including functional annotation databases and gene-based and pathway-based methods—to identify pathophysiologically relevant cellular mechanisms. Findings Our sample comprised 6167 patients with Barretts oesophagus and 4112 individuals with oesophageal adenocarcinoma, in addition to 17 159 representative controls from four genome-wide association studies in Europe, North America, and Australia. We identified eight new risk loci associated with either Barretts oesophagus or oesophageal adenocarcinoma, within or near the genes CFTR (rs17451754; p=4·8 × 10−10), MSRA (rs17749155; p=5·2 × 10−10), LINC00208 and BLK (rs10108511; p=2·1 × 10−9), KHDRBS2 (rs62423175; p=3·0 × 10−9), TPPP and CEP72 (rs9918259; p=3·2 × 10−9), TMOD1 (rs7852462; p=1·5 × 10−8), SATB2 (rs139606545; p=2·0 × 10−8), and HTR3C and ABCC5 (rs9823696; p=1·6 × 10−8). The locus identified near HTR3C and ABCC5 (rs9823696) was associated specifically with oesophageal adenocarcinoma (p=1·6 × 10−8) and was independent of Barretts oesophagus development (p=0·45). A ninth novel risk locus was identified within the gene LPA (rs12207195; posterior probability 0·925) after reweighting with significantly enriched annotations. The strongest disease pathways identified (p<10−6) belonged to muscle cell differentiation and to mesenchyme development and differentiation. Interpretation Our meta-analysis of genome-wide association studies doubled the number of known risk loci for Barretts oesophagus and oesophageal adenocarcinoma and revealed new insights into causes of these diseases. Furthermore, the specific association between oesophageal adenocarcinoma and the locus near HTR3C and ABCC5 might constitute a novel genetic marker for prediction of the transition from Barretts oesophagus to oesophageal adenocarcinoma. Fine-mapping and functional studies of new risk loci could lead to identification of key molecules in the development of Barretts oesophagus and oesophageal adenocarcinoma, which might encourage development of advanced prevention and intervention strategies. Funding US National Cancer Institute, US National Institutes of Health, National Health and Medical Research Council of Australia, Swedish Cancer Society, Medical Research Council UK, Cambridge NIHR Biomedical Research Centre, Cambridge Experimental Cancer Medicine Centre, Else Kröner Fresenius Stiftung, Wellcome Trust, Cancer Research UK, AstraZeneca UK, University Hospitals of Leicester, University of Oxford, Australian Research Council.

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

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

Chronic gastroesophageal reflux disease shares genetic background with esophageal adenocarcinoma and Barrett's esophagus

Esophageal adenocarcinoma (EA) is a rapidly fatal cancer with rising incidence in the developed world. Most EAs arise in a metaplastic epithelium, Barretts esophagus (BE), which is associated with greatly increased risk of EA. One of the key risk factors for both BE and EA is chronic gastroesophageal reflux disease (GERD). This study used the linkage disequilibrium (LD) score regression and genomic profile risk scoring approaches to investigate the contribution of multiple common single-nucleotide polymorphisms (SNPs) to the risk of GERD, and the extent of genetic overlap between GERD and BE or EA. Using LD score regression, we estimated an overall phenotypic variance of 7% (95% CI 3–11%) for GERD explained by all the genotyped SNPs. A genetic correlation of 77% (s.e. = 24%, P = 0.0012) between GERD and BE and 88% between GERD and EA (s.e. = 25%, P = 0.0004) was estimated using the LD score regression approach. Results from the genomic profile risk scoring approach, as a robustness check, were broadly similar to those from the LD score regression. This study provides the first evidence for a polygenic basis for GERD and supports for a polygenic overlap between GERD and BE, and GERD and EA.

Genome-wide association study identifies WNT7B as a novel locus for central corneal thickness in Latinos

The cornea is the outermost layer of the eye and is a vital component of focusing incoming light on the retina. Central corneal thickness (CCT) is now recognized to have a significant role in ocular health and is a risk factor for various ocular diseases, such as keratoconus and primary open angle glaucoma. Most previous genetic studies utilized European and Asian subjects to identify genetic loci associated with CCT. Minority populations, such as Latinos, may aid in identifying additional loci and improve our understanding of the genetic architecture of CCT. In this study, we conducted a genome-wide association study (GWAS) in Latinos, a traditionally understudied population in genetic research, to further identify loci contributing to CCT. Study participants were genotyped using either the Illumina OmniExpress BeadChip (∼730K markers) or the Illumina Hispanic/SOL BeadChip (∼2.5 million markers). All study participants were 40 years of age and older. We assessed the association between individual single nucleotide polymorphisms (SNPs) and CCT using linear regression, adjusting for age, gender and principal components of genetic ancestry. To expand genomic coverage and to interrogate additional SNPs, we imputed SNPs from the 1000 Genomes Project reference panels. We identified a novel SNP, rs10453441 (P = 6.01E-09), in an intron of WNT7B that is associated with CCT. Furthermore, WNT7B is expressed in the human cornea. We also replicated 11 previously reported loci, including IBTK, RXRA-COL5A1, COL5A1, FOXO1, LRRK1 and ZNF469 (P < 1.25E-3). These findings provide further insight into the genetic architecture of CCT and illustrate that the use of minority groups in GWAS will help identify additional loci.

Germline variation in inflammation-related pathways and risk of Barrett's oesophagus and oesophageal adenocarcinoma

Objective Oesophageal adenocarcinoma (OA) incidence has risen sharply in Western countries over recent decades. Local and systemic inflammation is considered an important contributor to OA pathogenesis. Established risk factors for OA and its precursor, Barretts oesophagus (BE), include symptomatic reflux, obesity and smoking. The role of inherited genetic susceptibility remains an area of active investigation. Here, we explore whether germline variation related to inflammatory processes influences susceptibility to BE/OA. Design We used data from a genomewide association study of 2515 OA cases, 3295 BE cases and 3207 controls. Our analysis included 7863 single-nucleotide polymorphisms (SNPs) in 449 genes assigned to five pathways: cyclooxygenase (COX), cytokine signalling, oxidative stress, human leucocyte antigen and nuclear factor-κB. A principal components-based analytic framework was employed to evaluate pathway-level and gene-level associations with disease risk. Results We identified a significant signal for the COX pathway in relation to BE risk (p=0.0059, false discovery rate q=0.03), and in gene-level analyses found an association with microsomal glutathione-S-transferase 1 (MGST1); (p=0.0005, q=0.005). Assessment of 36 MGST1 SNPs identified 14 variants associated with elevated BE risk (q<0.05). Four of these were subsequently confirmed (p<5.5×10−5) in a meta-analysis encompassing an independent set of 1851 BE cases and 3496 controls, and are known strong expression quantitative trait loci for MGST1. Three such variants were associated with similar elevations in OA risk. Conclusions This study provides the most comprehensive evaluation of inflammation-related germline variation in relation to risk of BE/OA and suggests that variants in MGST1 influence disease susceptibility.

Accurate Imputation-Based Screening of Gln368Ter Myocilin Variant in Primary Open-Angle Glaucoma

PURPOSE Myocilin (MYOC) is a well-established primary open-angle glaucoma (POAG) risk gene, with rare variants known to have high penetrance. The most common clinically relevant risk variant, Gln368Ter, has an allele frequency of 0.1% to 0.3% in populations of European ancestry. Detection of rare MYOC variants has traditionally been conducted using Sanger sequencing. Here we report the use of genotyping arrays and imputation to assess whether rare variants including Gln368Ter can be reliably detected. METHODS A total of 1155 cases with advanced POAG and 1992 unscreened controls genotyped on common variant arrays participated in this study. Accuracy of imputation of Gln368Ter variants was compared with direct sequencing. A genome-wide association study was performed using additive model adjusted for sex and the first six principal components. RESULTS We found that although the arrays we used were designed to tag common variants, we could reliably impute the Gln368Ter variant (rs74315329). When tested in 1155 POAG cases and 1992 controls, rs74315329 was strongly associated with risk (odds ratio = 15.53, P = 1.07 × 10-9). All POAG samples underwent full sequencing of the MYOC gene, and we found a sensitivity of 100%, specificity of 99.91%, positive predictive value of 95.65%, and negative predictive value of 100% between imputation and sequencing. Gln368Ter was also accurately imputed in a further set of 1801 individuals without POAG. Among the total set of 3793 (1992 + 1801) individuals without POAG, six were predicted (probability > 95%) to carry the risk variant. CONCLUSIONS We demonstrate that some clinically important rare variants can be reliably detected using arrays and imputation. These results have important implications for the detection of clinically relevant incidental findings in ongoing and future studies using arrays.