Search

Purpose: Observational studies have suggested that individuals with pre-existing sleep apnea (SA) have up to double the risk of developing glaucoma than individuals without SA. Understanding risk factors for glaucoma is important to assist with well-structured screening, early intervention, and efficient allocation of specialist consultation. The objective of this study is therefore to use genetic data to determine whether SA is a causal risk factor for glaucoma., Methods: Two-sample Mendelian randomization (MR) analyses were performed to assess the association between genetically predicted SA and glaucoma susceptibility using genome-wide association study (GWAS) of 25,062 SA cases, 313,372 controls derived from 23andMe and summary data from a glaucoma GWAS meta-analysis (20,582 cases, 119,318 controls), including individuals of European descent, mainly from the UK Biobank., Results: Inverse variance weighted regression of genetic susceptibility for SA on risk of glaucoma revealed no strong evidence for an association between SA and glaucoma (OR = 0.95, 95% confidence intervals = 0.84-1.07), results were consistent across all MR predictors., Conclusions: We found little genetic evidence supporting a causal association between SA and glaucoma. Our results refute the possibility of a large effect (glaucoma OR > 1.5 per doubling of odds on SA) between SA and glaucoma.

Keratoconus is characterised by reduced rigidity of the cornea with distortion and focal thinning that causes blurred vision, however, the pathogenetic mechanisms are unknown. It can lead to severe visual morbidity in children and young adults and is a common indication for corneal transplantation worldwide. Here we report the first large scale genome-wide association study of keratoconus including 4,669 cases and 116,547 controls. We have identified significant association with 36 genomic loci that, for the first time, implicate both dysregulation of corneal collagen matrix integrity and cell differentiation pathways as primary disease-causing mechanisms. The results also suggest pleiotropy, with some disease mechanisms shared with other corneal diseases, such as Fuchs endothelial corneal dystrophy. The common variants associated with keratoconus explain 12.5% of the genetic variance, which shows potential for the future development of a diagnostic test to detect susceptibility to disease.

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.

Importance: Keratoconus is a condition in which the cornea progressively thins and protrudes in a conical shape, severely affecting refraction and vision. It is a major indication for corneal transplant. To discover new genetic loci associated with keratoconus and better understand the causative mechanism of this disease, we performed a genome-wide association study on patients with keratoconus., Objective: To identify genetic susceptibility regions for keratoconus in the human genome., Design, Setting, and Participants: This study was conducted with data from eye clinics in Australia, the United States, and Northern Ireland. The discovery cohort of individuals with keratoconus and control participants from Australia was genotyped using the Illumina HumanCoreExome single-nucleotide polymorphism array. After quality control and data cleaning, genotypes were imputed against the 1000 Genomes Project reference panel (phase III; version 5), and association analyses were completed using PLINK version 1.90. Single-nucleotide polymorphisms with P < 1.00 × 10-6 were assessed for replication in 3 additional cohorts. Control participants were drawn from the cohorts of the Blue Mountains Eye Study and a previous study of glaucoma. Replication cohorts were from a previous keratoconus genome-wide association study data set from the United States, a cohort of affected and control participants from Australia and Northern Ireland, and a case-control cohort from Victoria, Australia. Data were collected from January 2006 to March 2019., Main Outcomes and Measures: Associations between keratoconus and 6 252 612 genetic variants were estimated using logistic regression after adjusting for ancestry using the first 3 principal components., Results: The discovery cohort included 522 affected individuals and 655 control participants, while the replication cohorts included 818 affected individuals (222 from the United States, 331 from Australia and Northern Ireland, and 265 from Victoria, Australia) and 3858 control participants (2927 from the United States, 229 from Australia and Northern Ireland, and 702 from Victoria, Australia). Two novel loci reached genome-wide significance (defined as P < 5.00 × 10-8), with a P value of 7.46 × 10-9 at rs61876744 in patatin-like phospholipase domain-containing 2 gene (PNPLA2) on chromosome 11 and a P value of 6.35 × 10-12 at rs138380, 2.2 kb upstream of casein kinase I isoform epsilon gene (CSNK1E) on chromosome 22. One additional locus was identified with a P value less than 1.00 × 10-6 in mastermind-like transcriptional coactivator 2 (MAML2) on chromosome 11 (P = 3.91 × 10-7). The novel locus in PNPLA2 reached genome-wide significance in an analysis of all 4 cohorts (P = 2.45 × 10-8)., Conclusions and Relevance: In this relatively large keratoconus genome-wide association study, we identified a genome-wide significant locus for keratoconus in the region of PNPLA2 on chromosome 11.

Glaucoma, a disease characterized by progressive optic nerve degeneration, can be prevented through timely diagnosis and treatment. We characterize optic nerve photographs of 67,040 UK Biobank participants and use a multitrait genetic model to identify risk loci for glaucoma. A glaucoma polygenic risk score (PRS) enables effective risk stratification in unselected glaucoma cases and modifies penetrance of the MYOC variant encoding p.Gln368Ter, the most common glaucoma-associated myocilin variant. In the unselected glaucoma population, individuals in the top PRS decile reach an absolute risk for glaucoma 10 years earlier than the bottom decile and are at 15-fold increased risk of developing advanced glaucoma (top 10% versus remaining 90%, odds ratio = 4.20). The PRS predicts glaucoma progression in prospectively monitored, early manifest glaucoma cases (P = 0.004) and surgical intervention in advanced disease (P = 3.6 × 10 - 6 ). This glaucoma PRS will facilitate the development of a personalized approach for earlier treatment of high-risk individuals, with less intensive monitoring and treatment being possible for lower-risk groups.

Background & Aims: Epidemiology studies of circulating concentrations of 25 hydroxy vitamin D (25(OH)D) and risk of esophageal adenocarcinoma (EAC) have produced conflicting results. We conducted a Mendelian randomization study to determine the associations between circulating concentrations of 25(OH)D and risks of EAC and its precursor, Barrett's esophagus (BE)., Methods: We conducted a Mendelian randomization study using a 2-sample (summary data) approach. Six single-nucleotide polymorphisms (SNPs; rs3755967, rs10741657, rs12785878, rs10745742, rs8018720, and rs17216707) associated with circulating concentrations of 25(OH)D were used as instrumental variables. We collected data from 6167 patients with BE, 4112 patients with EAC, and 17,159 individuals without BE or EAC (controls) participating in the Barrett's and Esophageal Adenocarcinoma Consortium, as well as studies from Bonn, Germany, and Cambridge and Oxford, United Kingdom. Analyses were performed separately for BE and EAC., Results: Overall, we found no evidence for an association between genetically estimated 25(OH)D concentration and risk of BE or EAC. The odds ratio per 20 nmol/L increase in genetically estimated 25(OH)D concentration for BE risk estimated by combining the individual SNP association using inverse variance weighting was 1.21 (95% CI, 0.77-1.92; P = .41). The odds ratio for EAC risk, estimated by combining the individual SNP association using inverse variance weighting, was 0.68 (95% CI, 0.39-1.19; P = .18)., Conclusions: In a Mendelian randomization study, we found that low genetically estimated 25(OH)D concentrations were not associated with risk of BE or EAC., (Copyright © 2019 AGA Institute. Published by Elsevier Inc. All rights reserved.)

With the advent of very large scale genome-wide association studies (GWASs), the promise of Mendelian randomization (MR) has begun to be fulfilled. However, whilst GWASs have provided essential information on the single nucleotide polymorphisms (SNPs) associated with modifiable risk factors needed for MR, the availability of large numbers of SNP instruments raises issues of how best to use this information and how to deal with potential problems such as pleiotropy. Here we provide commentary on some of the recent advances in the MR analysis, including an overview of the different genetic architectures that are being uncovered for a variety of modifiable risk factors and how users ought to take that into consideration when designing MR studies., (© 2019 The Authors. Genetic Epidemiology Published by Wiley Periodicals, Inc.)

Circulating levels of adiponectin, an adipocyte-secreted protein associated with cardiovascular and metabolic risk, are highly heritable. To gain insights into the biology that regulates adiponectin levels, we performed an exome array meta-analysis of 265,780 genetic variants in 67,739 individuals of European, Hispanic, African American, and East Asian ancestry. We identified 20 loci associated with adiponectin, including 11 that had been reported previously (p < 2 × 10 -7 ). Comparison of exome array variants to regional linkage disequilibrium (LD) patterns and prior genome-wide association study (GWAS) results detected candidate variants (r 2 > .60) spanning as much as 900 kb. To identify potential genes and mechanisms through which the previously unreported association signals act to affect adiponectin levels, we assessed cross-trait associations, expression quantitative trait loci in subcutaneous adipose, and biological pathways of nearby genes. Eight of the nine loci were also associated (p < 1 × 10 -4 ) with at least one obesity or lipid trait. Candidate genes include PRKAR2A, PTH1R, and HDAC9, which have been suggested to play roles in adipocyte differentiation or bone marrow adipose tissue. Taken together, these findings provide further insights into the processes that influence circulating adiponectin levels., (Copyright © 2019 American Society of Human Genetics. All rights reserved.)

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.

Purpose: To identify genetic variants conferring susceptibility to esotropia. Esotropia is the most common form of comitant strabismus, has its highest incidence in European ancestry populations, and is believed to be inherited as a complex trait., Methods: White European American discovery cohorts with nonaccommodative (826 cases and 2991 controls) or accommodative (224 cases and 749 controls) esotropia were investigated. White European Australian and United Kingdom cohorts with nonaccommodative (689 cases and 1448 controls) or accommodative (66 cases and 264 controls) esotropia were tested for replication. We performed a genome-wide case-control association study using a mixed linear additive model. Meta-analyses of discovery and replication cohorts were then conducted., Results: A significant association with nonaccommodative esotropia was discovered (odds ratio [OR] = 1.41, P = 2.84 × 10-09) and replicated (OR = 1.23, P = 0.01) at rs2244352 [T] located within intron 1 of the WRB (tryptophan rich basic protein) gene on chromosome 21 (meta-analysis OR = 1.33, P = 9.58 × 10-11). This single nucleotide polymorphism (SNP) is differentially methylated, and there is a statistically significant skew toward paternal inheritance in the discovery cohort. Meta-analysis of the accommodative discovery and replication cohorts identified an association with rs912759 [T] (OR = 0.59, P = 1.89 × 10-08), an intergenic SNP on chromosome 1p31.1., Conclusions: This is the first genome-wide association study (GWAS) to identify significant associations in esotropia and suggests a parent-of-origin effect. Additional cohorts will permit replication and extension of these findings. Future studies of rs2244352 and WRB should provide insight into pathophysiological mechanisms underlying comitant strabismus.

Central corneal thickness (CCT) is a highly heritable trait associated with complex eye diseases such as keratoconus and glaucoma. We perform a genome-wide association meta-analysis of CCT and identify 19 novel regions. In addition to adding support for known connective tissue-related pathways, pathway analyses uncover previously unreported gene sets. Remarkably, >20% of the CCT-loci are near or within Mendelian disorder genes. These included FBN1, ADAMTS2 and TGFB2 which associate with connective tissue disorders (Marfan, Ehlers-Danlos and Loeys-Dietz syndromes), and the LUM-DCN-KERA gene complex involved in myopia, corneal dystrophies and cornea plana. Using index CCT-increasing variants, we find a significant inverse correlation in effect sizes between CCT and keratoconus (r = -0.62, P = 5.30 × 10 -5 ) but not between CCT and primary open-angle glaucoma (r = -0.17, P = 0.2). Our findings provide evidence for shared genetic influences between CCT and keratoconus, and implicate candidate genes acting in collagen and extracellular matrix regulation.

Purpose: To identify genetic risk factors contributing to central corneal thickness (CCT) in individuals from South India, a population with a high prevalence of ocular disorders., Methods: One hundred ninety-five individuals from 15 large South Indian pedigrees were genotyped using the Omni2.5 bead array. Family-based association for CCT was conducted using the score test in MERLIN., Results: Genome-wide association study (GWAS) identified strongest association for single nucleotide polymorphisms (SNPs) in the first intron of WNT7B and CCT (top SNP rs9330813; β = -0.57, 95% confidence interval [CI]: -0.78 to -0.36; P = 1.7 × 10-7). We further investigated rs9330813 in a Latino cohort and four independent European cohorts. A meta-analysis of these data sets demonstrated statistically significant association between rs9330813 and CCT (β = -3.94, 95% CI: -5.23 to -2.66; P = 1.7 × 10-9). WNT7B SNPs located in the same genomic region that includes rs9330813 have previously been associated with CCT in Latinos but with other ocular quantitative traits related to myopia (corneal curvature and axial length) in a Japanese population (rs10453441 and rs200329677). To evaluate the specificity of the observed WNT7B association with CCT in the South Indian families, we completed an ocular phenome-wide association study (PheWAS) for the top WNT7B SNPs using 45 ocular traits measured in these same families including corneal curvature and axial length. The ocular PheWAS results indicate that in the South Indian families WNT7B SNPs are primarily associated with CCT., Conclusions: The results indicate robust evidence for association between WNT7B SNPs and CCT in South Indian pedigrees, and suggest that WNT7B SNPs can have population-specific effects on ocular quantitative traits.

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.

Cohort studies investigating aging and dementia require APOE genotyping. We compared directly measured APOE genotypes to 'hard-call' genotypes derived from imputing genome-wide genotyping data from a range of platforms using several imputation panels. Older GWAS arrays imputed to 1000 Genomes Project (1KGP) phases and the Haplotype Reference Consortium (HRC) reference panels were able to achieve concordance rates of over 98% with stringent quality control (hard-call-threshold 0.8). However, this resulted in high levels of missingness (>12% with 1KGP and 5% with HRC). With recent GWAS arrays, concordance of 99% could be obtained with relatively lenient QC, resulting in no missingness.

Previous genome-wide association studies have identified a melanoma-associated locus at 1q42.1 that encompasses a ∼100-kb region spanning the PARP1 gene. Expression quantitative trait locus (eQTL) analysis in multiple cell types of the melanocytic lineage consistently demonstrated that the 1q42.1 melanoma risk allele (rs3219090[G]) is correlated with higher PARP1 levels. In silico fine-mapping and functional validation identified a common intronic indel, rs144361550 (-/GGGCCC; r 2 = 0.947 with rs3219090), as displaying allele-specific transcriptional activity. A proteomic screen identified RECQL as binding to rs144361550 in an allele-preferential manner. In human primary melanocytes, PARP1 promoted cell proliferation and rescued BRAF V600E -induced senescence phenotypes in a PARylation-independent manner. PARP1 also transformed TERT-immortalized melanocytes expressing BRAF V600E . PARP1-mediated senescence rescue was accompanied by transcriptional activation of the melanocyte-lineage survival oncogene MITF, highlighting a new role for PARP1 in melanomagenesis.

Previous studies have identified many genetic loci for refractive error and myopia. We aimed to investigate the effect of these loci on ocular biometry as a function of age in children, adolescents, and adults. The study population consisted of three age groups identified from the international CREAM consortium: 5,490 individuals aged <10 years; 5,000 aged 10-25 years; and 16,274 aged >25 years. All participants had undergone standard ophthalmic examination including measurements of axial length (AL) and corneal radius (CR). We examined the lead SNP at all 39 currently known genetic loci for refractive error identified from genome-wide association studies (GWAS), as well as a combined genetic risk score (GRS). The beta coefficient for association between SNP genotype or GRS versus AL/CR was compared across the three age groups, adjusting for age, sex, and principal components. Analyses were Bonferroni-corrected. In the age group <10 years, three loci (GJD2, CHRNG, ZIC2) were associated with AL/CR. In the age group 10-25 years, four loci (BMP2, KCNQ5, A2BP1, CACNA1D) were associated; and in adults 20 loci were associated. Association with GRS increased with age; β = 0.0016 per risk allele (P = 2 × 10 -8 ) in <10 years, 0.0033 (P = 5 × 10 -15 ) in 10- to 25-year-olds, and 0.0048 (P = 1 × 10 -72 ) in adults. Genes with strongest effects (LAMA2, GJD2) had an early effect that increased with age. Our results provide insights on the age span during which myopia genes exert their effect. These insights form the basis for understanding the mechanisms underlying high and pathological myopia., (© 2016 WILEY PERIODICALS, INC.)

Background: In vitro and observational epidemiological studies suggest that vitamin D may play a role in cancer prevention. However, the relationship between vitamin D and ovarian cancer is uncertain, with observational studies generating conflicting findings. A potential limitation of observational studies is inadequate control of confounding. To overcome this problem, we used Mendelian randomization (MR) to evaluate the association between single nucleotide polymorphisms (SNPs) associated with circulating 25-hydroxyvitamin D [25(OH)D] concentration and risk of ovarian cancer., Methods: We employed SNPs with well-established associations with 25(OH)D concentration as instrumental variables for MR: rs7944926 (DHCR7), rs12794714 (CYP2R1) and rs2282679 (GC). We included 31 719 women of European ancestry (10 065 cases, 21 654 controls) from the Ovarian Cancer Association Consortium, who were genotyped using customized Illumina Infinium iSelect (iCOGS) arrays. A two-sample (summary data) MR approach was used and analyses were performed separately for all ovarian cancer (10 065 cases) and for high-grade serous ovarian cancer (4121 cases)., Results: The odds ratio for epithelial ovarian cancer risk (10 065 cases) estimated by combining the individual SNP associations using inverse variance weighting was 1.27 (95% confidence interval: 1.06 to 1.51) per 20 nmol/L decrease in 25(OH)D concentration. The estimated odds ratio for high-grade serous epithelial ovarian cancer (4121 cases) was 1.54 (1.19, 2.01)., Conclusions: Genetically lowered 25-hydroxyvitamin D concentrations were associated with higher ovarian cancer susceptibility in Europeans. These findings suggest that increasing plasma vitamin D levels may reduce risk of ovarian cancer., (© The Author 2016; all rights reserved. Published by Oxford University Press on behalf of the International Epidemiological Association.)

Purpose: Many genome-wide association studies have identified common single nucleotide polymorphisms (SNPs) at the 9p21 glaucoma locus (CDKN2B/CDKN2B-AS1) to be significantly associated with primary open-angle glaucoma (POAG), with association being stronger in normal tension glaucoma (NTG) and advanced glaucoma. We aimed to determine whether any observed differences in genetic association at the 9p21 locus are influenced by sex., Methods: Sex was assessed as a risk factor for POAG for 2241 glaucoma participants from the Australian and New Zealand Registry of Advanced Glaucoma, the Glaucoma Inheritance Study in Tasmania, and the Flinders Medical Centre. A total of 3176 controls were drawn from the Blue Mountains Eye Study and South Australia: 1523 advanced POAG and 718 nonadvanced POAG cases were genotyped along with 3176 controls. We selected 13 SNPs at the 9p21 locus, and association results were subanalyszd by sex for high-tension glaucoma (HTG) and NTG. Odds ratios (ORs) between sexes were compared., Results: A sex bias was present within advanced NTG cases (57.1% female versus 42.9% male, P = 0.0026). In all POAG cases, the strongest associated SNP at 9p21 was rs1063192 (OR, 1.43; P = 4 × 10-18). This association was stronger in females (OR, 1.5; P = 5 × 10-13) than in males (OR, 1.35; P = 7 × 10-7), with a statistically significant difference in female to male OR comparison (P = 1.0 × 10-2). An NTG to HTG subanalysis yielded statistically significant results only in females (OR, 1.63; P = 1.5 × 10-4) but not in males (OR, 1.15; P = 2.8 × 10-1), with a statistically significant difference in female to male OR comparison (P = 1.4 × 10-4)., Conclusions: This study demonstrated that female sex is a risk factor for developing advanced NTG. The stronger genetic signals at the 9p21 locus among females may contribute at least in part to the observed sex bias for NTG.

Purpose: Thyroid-associated orbitopathy (TO) is an autoimmune-mediated orbital inflammation that can lead to disfigurement and blindness. Multiple genetic loci have been associated with Graves' disease, but the genetic basis for TO is largely unknown. This study aimed to identify loci associated with TO in individuals with Graves' disease, using a genome-wide association scan (GWAS) for the first time to our knowledge in TO., Methods: Genome-wide association scan was performed on pooled DNA from an Australian Caucasian discovery cohort of 265 participants with Graves' disease and TO (cases) and 147 patients with Graves' disease without TO (controls). Top-ranked single nucleotide polymorphisms (SNPs) then were genotyped in individual DNA samples from the discovery cohort, and two replication cohorts totaling 584 cases and 367 controls., Results: In the GWAS of pooled DNA samples, several SNPs showed suggestive association with TO at genome-wide P ≤ 10-6; rs953128 located on chr10q21.1, rs2867161 on chr7q11.22, rs13360861 on chr5q12.3, rs7636326 on chr3q26.2, rs10266576 on chr 7q11.22, rs60457622 on chr3q23, and rs6110809 on chr20p12.1. However, the only SNP consistently associated with TO on individual genotyping in the discovery and replication cohorts was rs6110809, located within MACROD2 on chromosome 20p12.1. On combined analysis of discovery and replication cohorts, the minor A allele of rs6110809 was more frequent in TO than in Graves' disease controls without TO (P = 4.35 × 10-5; odds ratio [OR] = 1.77; 95% confidence interval [CI], 1.35-2.32) after adjusting for age, sex, duration of Graves' disease, and smoking., Conclusions: In patients with Graves' disease, a common genetic variant in MACROD2 may increase susceptibility for thyroid-associated orbitopathy. This association now requires confirmation in additional independent cohorts.

The strong male predominance in Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) remains inadequately explained, but sex hormones might be involved. We hypothesized that single nucleotide polymorphisms (SNPs) in the androgen pathway influence risk of developing BE and EAC. This genetic-epidemiological analysis included 14 studies from Australia, Europe and North America. Polymorphisms in 16 genes coding for the androgen pathway were analyzed using a gene-based approach: versatile gene-based test association study. This method evaluates associations between a trait and all SNPs within a specific gene rather than each SNP marker individually as in a conventional GWAS. The data were stratified for sex, body-mass index, waist-to-hip ratio, tobacco smoking and gastroesophageal reflux status. Included were data from 1,508 EAC patients, 2,383 BE patients and 2,170 control participants. SNPs within the gene CYP17A1 were associated with risk of BE in the sexes combined (p = 0.002) and in males (p = 0.003), but not in females separately (p = 0.3). This association was found in tobacco smokers (p = 0.003) and in BE patients without reflux (p = 0.004), but not in nonsmokers (p = 0.2) or those with reflux (p = 0.036). SNPs within JMJD1C were associated with risk of EAC in females (p = 0.001). However, none of these associations replicated in a subsequent sample. Fourteen other genes studied did not reach statistically significant levels of association with BE, EAC or the combination of BE and EAC, after correcting for the number of genes included in the analysis. In conclusion, genetic variants in the androgen-related genes CYP17A1 and JMJD1C might be associated with risk of BE and EAC, respectively, but replication data with larger sample sizes are needed., (© 2015 UICC.)

Women with epithelial ovarian cancer (EOC) are usually treated with platinum/taxane therapy after cytoreductive surgery but there is considerable inter-individual variation in response. To identify germline single-nucleotide polymorphisms (SNPs) that contribute to variations in individual responses to chemotherapy, we carried out a multi-phase genome-wide association study (GWAS) in 1,244 women diagnosed with serous EOC who were treated with the same first-line chemotherapy, carboplatin and paclitaxel. We identified two SNPs (rs7874043 and rs72700653) in TTC39B (best P=7x10-5, HR=1.90, for rs7874043) associated with progression-free survival (PFS). Functional analyses show that both SNPs lie in a putative regulatory element (PRE) that physically interacts with the promoters of PSIP1, CCDC171 and an alternative promoter of TTC39B. The C allele of rs7874043 is associated with poor PFS and showed increased binding of the Sp1 transcription factor, which is critical for chromatin interactions with PSIP1. Silencing of PSIP1 significantly impaired DNA damage-induced Rad51 nuclear foci and reduced cell viability in ovarian cancer lines. PSIP1 (PC4 and SFRS1 Interacting Protein 1) is known to protect cells from stress-induced apoptosis, and high expression is associated with poor PFS in EOC patients. We therefore suggest that the minor allele of rs7874043 confers poor PFS by increasing PSIP1 expression.

Epidemiological studies have demonstrated associations between endometriosis and certain histotypes of ovarian cancer, including clear cell, low-grade serous and endometrioid carcinomas. We aimed to determine whether the observed associations might be due to shared genetic aetiology. To address this, we used two endometriosis datasets genotyped on common arrays with full-genome coverage (3194 cases and 7060 controls) and a large ovarian cancer dataset genotyped on the customized Illumina Infinium iSelect (iCOGS) arrays (10 065 cases and 21 663 controls). Previous work has suggested that a large number of genetic variants contribute to endometriosis and ovarian cancer (all histotypes combined) susceptibility. Here, using the iCOGS data, we confirmed polygenic architecture for most histotypes of ovarian cancer. This led us to evaluate if the polygenic effects are shared across diseases. We found evidence for shared genetic risks between endometriosis and all histotypes of ovarian cancer, except for the intestinal mucinous type. Clear cell carcinoma showed the strongest genetic correlation with endometriosis (0.51, 95% CI = 0.18-0.84). Endometrioid and low-grade serous carcinomas had similar correlation coefficients (0.48, 95% CI = 0.07-0.89 and 0.40, 95% CI = 0.05-0.75, respectively). High-grade serous carcinoma, which often arises from the fallopian tubes, showed a weaker genetic correlation with endometriosis (0.25, 95% CI = 0.11-0.39), despite the absence of a known epidemiological association. These results suggest that the epidemiological association between endometriosis and ovarian adenocarcinoma may be attributable to shared genetic susceptibility loci., (© The Author 2015. Published by Oxford University Press. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Background: The strong male predominance in oesophageal adenocarcinoma (OAC) and Barrett's oesophagus (BO) continues to puzzle. Hormonal influence, e.g. oestrogen or oxytocin, might contribute., Methods: This genetic-epidemiological study pooled 14 studies from three continents, Australia, Europe, and North America. Polymorphisms in 3 key genes coding for the oestrogen pathway (receptor alpha (ESR1), receptor beta (ESR2), and aromatase (CYP19A1)), and 3 key genes of the oxytocin pathway (the oxytocin receptor (OXTR), oxytocin protein (OXT), and cyclic ADP ribose hydrolase glycoprotein (CD38)), were analysed using a gene-based approach, versatile gene-based test association study (VEGAS)., Results: Among 1508 OAC patients, 2383 BO patients, and 2170 controls, genetic variants within ESR1 were associated with BO in males (p = 0.0058) and an increased risk of OAC and BO combined in males (p = 0.0023). Genetic variants within OXTR were associated with an increased risk of BO in both sexes combined (p = 0.0035) and in males (p = 0.0012). We followed up these suggestive findings in a further smaller data set, but found no replication. There were no significant associations between the other 4 genes studied and risk of OAC, BO, separately on in combination, in males and females combined or in males only., Conclusion: Genetic variants in the oestrogen receptor alpha and the oxytocin receptor may be associated with an increased risk of BO or OAC, but replication in other large samples are needed.

Primary open angle glaucoma (POAG), a major cause of blindness worldwide, is a complex disease with a significant genetic contribution. We performed Exome Array (Illumina) analysis on 3504 POAG cases and 9746 controls with replication of the most significant findings in 9173 POAG cases and 26 780 controls across 18 collections of Asian, African and European descent. Apart from confirming strong evidence of association at CDKN2B-AS1 (rs2157719 [G], odds ratio [OR] = 0.71, P = 2.81 × 10(-33)), we observed one SNP showing significant association to POAG (CDC7-TGFBR3 rs1192415, ORG-allele = 1.13, Pmeta = 1.60 × 10(-8)). This particular SNP has previously been shown to be strongly associated with optic disc area and vertical cup-to-disc ratio, which are regarded as glaucoma-related quantitative traits. Our study now extends this by directly implicating it in POAG disease pathogenesis., (© The Author 2015. Published by Oxford University Press.)

Importance: Conjunctival UV autofluorescence (CUVAF) photography was developed to detect and characterize preclinical sunlight-induced ocular damage. Ocular sun exposure has been related to cases of pterygia and was recently negatively correlated with myopia. Hence, CUVAF has excellent potential as an objective biomarker of sun exposure. However, much variation in CUVAF has been observed, and the relative contributions of genes and environment to this variation have not yet been identified., Objective: To investigate sources of variation in CUVAF in relation to its potential clinical relevance., Design, Setting, and Participants: We performed a cross-sectional analysis of 3 population-based cohort studies in the general community, including the Twins Eye Study in Tasmania, the Brisbane Adolescent Twin Study, and the Western Australian Pregnancy Cohort (Raine) Study. The twin studies were conducted between 2001 and 2009, and the 20-year follow-up of the Raine Study was completed between March 2010 and February 2012. We included genotypic and phenotypic data from 295 Australian families in the Tasmanian and Brisbane twin studies and from 661 participants in the 20-year follow-up of the Raine Study. We compared CUVAF levels in the 3 cohorts and performed a classic twin study to partition variation in CUVAF. We also conducted a genome-wide association analysis to identify specific genetic variants associated with CUVAF., Main Outcomes and Measures: The total area of CUVAF, heritability of CUVAF, and single-nucleotide polymorphisms (SNPs) associated with CUVAF from the genome-wide association study., Results: Within twin cohorts, individuals living closer to the equator (latitude, 27.47° S) had higher levels of CUVAF compared with individuals from southern regions (latitude, 42.88° S) (median [interquartile range], 45.4 [26.8-68.5] vs 28.7 [15.0-42.3] mm2; P < .001). The variation in CUVAF explained by the additive genetic component was 0.37 (95% CI, 0.22-0.56), whereas the variation due to the common environment was 0.50 (95% CI; 0.29-0.71). The SNP rs1060043, located approximately 800 base pairs away from the SLC1A5 gene, a member of the solute carrier family 1, had a genome-wide significant association with a P value of 3.2 × 10-8. Gene-based analysis did not improve our power to detect association with other genes., Conclusions and Relevance: Our findings confirm that, although a large environmental component to CUVAF (equivalent of sun exposure) exists, genes also play a significant role. We identified a SNP (rs1060043) as being significantly associated with CUVAF; replication of this finding in future studies is warranted.

Gene-based tests such as versatile gene-based association study (VEGAS) are commonly used following per-single nucleotide polymorphism (SNP) GWAS (genome-wide association studies) analysis. Two limitations of VEGAS were that the HapMap2 reference set was used to model the correlation between SNPs and only autosomal genes were considered. HapMap2 has now been superseded by the 1,000 Genomes reference set, and whereas early GWASs frequently ignored the X chromosome, it is now commonly included. Here we have developed VEGAS2, an extension that uses 1,000 Genomes data to model SNP correlations across the autosomes and chromosome X. VEGAS2 allows greater flexibility when defining gene boundaries. VEGAS2 offers both a user-friendly, web-based front end and a command line Linux version. The online version of VEGAS2 can be accessed through https://vegas2.qimrberghofer.edu.au/. The command line version can be downloaded from https://vegas2.qimrberghofer.edu.au/zVEGAS2offline.tgz. The command line version is developed in Perl, R and shell scripting languages; source code is available for further development.

Study Question: Are single-nucleotide polymorphisms (SNPs) at the interleukin 1A (IL1A) gene locus associated with endometriosis risk?, Summary Answer: We found evidence for strong association between IL1A SNPs and endometriosis risk., What Is Known Already: Genetic factors contribute substantially to the complex aetiology of endometriosis and the disease has an estimated heritability of ∼51%. We, and others, have conducted genome-wide association (GWA) studies for endometriosis, which identified a total of nine independent risk loci. Recently, two small Japanese studies reported eight SNPs (rs6542095, rs11677416, rs3783550, rs3783525, rs3783553, rs2856836, rs1304037 and rs17561) at the IL1A gene locus as suggestively associated with endometriosis risk. There is also evidence of a link between inflammation and endometriosis., Study Design, Size, Duration: We sought to further investigate the eight IL1A SNPs for association with endometriosis using an independent sample of 3908 endometriosis cases and 8568 controls of European and Japanese ancestry. The study was conducted between October 2013 and July 2014., Participants/materials, Setting, Methods: By leveraging GWA data from our previous multi-ethnic GWA meta-analysis for endometriosis, we imputed variants in the IL1A region, using a recent 1000 Genomes reference panel. After combining summary statistics for the eight SNPs from our European and Japanese imputed data with the published results, a fixed-effect meta-analysis was performed. An additional meta-analysis restricted to endometriosis cases with moderate-to-severe (revised American Fertility Society stage 3 or 4) disease versus controls was also performed., Main Results and the Role of Chance: All eight IL1A SNPs successfully replicated at P < 0.014 in the European imputed data with concordant direction and similar size to the effects reported in the original Japanese studies. Of these, three SNPs (rs6542095, rs3783550 and rs3783525) also showed association with endometriosis at a nominal P < 0.05 in our independent Japanese sample. Fixed-effect meta-analysis of the eight SNPs for moderate-to-severe endometriosis produced a genome-wide significant association for rs6542095 (odds ratio = 1.21; 95% confidence interval = 1.13-1.29; P = 3.43 × 10(-8))., Limitations, Reasons for Caution: The meta-analysis for moderate-to-severe endometriosis included results of moderate-to-severe endometriosis cases from our European data sets and all endometriosis cases from the Japanese data sets, as disease stage information was not available for endometriosis cases in the Japanese data sets., Wider Implications of the Findings: SNP rs6542095 is located ∼2.3 kb downstream of the IL1A gene and ∼6.9 kb upstream of cytoskeleton-associated protein 2-like (CKAP2L) gene. The IL1A gene encodes the IL1a protein, a member of the interleukin 1 cytokine family which is involved in various immune responses and inflammatory processes. These results provide important replication in an independent Japanese sample and, for the first time, association of the IL1A locus in endometriosis patients of European ancestry. SNPs within the IL1A locus may regulate other genes, but if IL1A is the target, our results provide supporting evidence for a link between inflammatory responses and the pathogenesis of endometriosis., Study Funding/competing Interests: The research was funded by grants from the Australian National Health and Medical Research Council and Wellcome Trust. None of the authors has competing interests for the study., (© The Author 2014. Published by Oxford University Press on behalf of the European Society of Human Reproduction and Embryology. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

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.

Background: Data from observational studies suggest that body mass index (BMI) is causally related to esophageal adenocarcinoma (EAC) and its precursor, Barrett's esophagus (BE). However, the relationships may be affected by bias and confounding., Methods: We used data from the Barrett's and Esophageal Adenocarcinoma Genetic Susceptibility Study: 999 patients with EAC, 2061 patients with BE, and 2169 population controls. We applied the two-stage control function instrumental variable method of the Mendelian randomization approach to estimate the unbiased, unconfounded effect of BMI on risk of EAC and BE. This was performed using a genetic risk score, derived from 29 genetic variants shown to be associated with BMI, as an instrument for lifetime BMI. A higher score indicates propensity to obesity. All tests were two-sided., Results: The genetic risk score was not associated with potential confounders, including gastroesophageal reflux symptoms and smoking. In the instrumental variable analyses (IV), EAC risk increased by 16% (IV-odds ratio [OR] = 1.16, 95% confidence interval [CI] = 1.01 to 1.33) and BE risk increased by 12% (IV-OR = 1.12, 95% CI = 1.00 to 1.25) per 1kg/m(2) increase in BMI. BMI was statistically significantly associated with EAC and BE in conventional epidemiologic analyses., Conclusions: People with a high genetic propensity to obesity have higher risks of esophageal metaplasia and neoplasia than people with low genetic propensity. These analyses provide the strongest evidence to date that obesity is independently associated with BE and EAC, and is not due to confounding or bias inherent in conventional epidemiologic analyses., (© The Author 2014. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oup.com.)

Telomere length has been associated with risk of many cancers, but results are inconsistent. Seven single nucleotide polymorphisms (SNPs) previously associated with mean leukocyte telomere length were either genotyped or well-imputed in 11108 case patients and 13933 control patients from Europe, Israel, the United States and Australia, four of the seven SNPs reached a P value under .05 (two-sided). A genetic score that predicts telomere length, derived from these seven SNPs, is strongly associated (P = 8.92x10(-9), two-sided) with melanoma risk. This demonstrates that the previously observed association between longer telomere length and increased melanoma risk is not attributable to confounding via shared environmental effects (such as ultraviolet exposure) or reverse causality. We provide the first proof that multiple germline genetic determinants of telomere length influence cancer risk., (© The Author 2014. Published by Oxford University Press.)

Background: ATP-binding cassette (ABC) transporters play various roles in cancer biology and drug resistance, but their association with outcomes in serous epithelial ovarian cancer (EOC) is unknown., Methods: The relationship between clinical outcomes and ABC transporter gene expression in two independent cohorts of high-grade serous EOC tumors was assessed with real-time quantitative polymerase chain reaction, analysis of expression microarray data, and immunohistochemistry. Associations between clinical outcomes and ABCA transporter gene single nucleotide polymorphisms were tested in a genome-wide association study. Impact of short interfering RNA-mediated gene suppression was determined by colony forming and migration assays. Association with survival was assessed with Kaplan-Meier analysis and log-rank tests. All statistical tests were two-sided., Results: Associations with outcome were observed with ABC transporters of the "A" subfamily, but not with multidrug transporters. High-level expression of ABCA1, ABCA6, ABCA8, and ABCA9 in primary tumors was statistically significantly associated with reduced survival in serous ovarian cancer patients. Low levels of ABCA5 and the C-allele of rs536009 were associated with shorter overall survival (hazard ratio for death = 1.50; 95% confidence interval [CI] =1.26 to 1.79; P = 6.5e-6). The combined expression pattern of ABCA1, ABCA5, and either ABCA8 or ABCA9 was associated with particularly poor outcome (mean overall survival in group with adverse ABCA1, ABCA5 and ABCA9 gene expression = 33.2 months, 95% CI = 26.4 to 40.1; vs 55.3 months in the group with favorable ABCA gene expression, 95% CI = 49.8 to 60.8; P = .001), independently of tumor stage or surgical debulking status. Suppression of cholesterol transporter ABCA1 inhibited ovarian cancer cell growth and migration in vitro, and statin treatment reduced ovarian cancer cell migration., Conclusions: Expression of ABCA transporters was associated with poor outcome in serous ovarian cancer, implicating lipid trafficking as a potentially important process in EOC., (© The Author 2014. Published by Oxford University Press. All rights reserved.)

Purpose: A recent genome-wide association study (GWAS) identified six loci associated with central corneal thickness that also conferred associated risk of keratoconus (KC). We aimed to assess whether genetic associations existed for these loci with KC or corneal curvature in an independent cohort of European ancestry., Methods: In total, 157 patients with KC were recruited from public and private clinics in Melbourne, Australia, and 673 individuals without KC were identified through the Genes in Myopia study from Australia. The following six single-nucleotide polymorphisms (SNPs) that showed a statistically significant association with KC in a recent GWAS study were selected for genotyping in our cohort: rs4894535 (FNDC3B), rs1324183 (MPDZ-NF1B), rs1536482 (RXRA-COL5A1), rs7044529 (COL5A), rs2721051 (FOXO1), and rs9938149 (BANP-ZNF469). The SNPs were assessed for their association with KC or corneal curvature using logistic or linear regression methods, with age and sex included as covariates. Bonferroni corrections were applied to account for multiple testing., Results: Genotyping data were available for five of the SNPs. Statistically significant associations with KC were found for the SNPs rs1324183 (P = 0.001; odds ratio [OR], 1.68) and rs9938149 (P = 0.010; OR, 1.47). Meta-analysis of previous studies yielded genome-wide significant evidence of an association for rs1324183, firmly establishing it as a KC risk variant. None of the SNPs were significantly associated with corneal curvature., Conclusions: The SNPs rs1324183 in the MPDZ-NF1B gene and rs9938149 (between BANP and ZNF4659) were associated with KC in this independent cohort, but their association was via a non-corneal curvature route.

Background: Esophageal adenocarcinoma (EA) is an increasingly common cancer with poor survival. Barrett's esophagus (BE) is the main precursor to EA, and every year 0.12% to 0.5% of BE patients progress to EA. BE typically arises on a background of chronic gastroesophageal reflux (GERD), one of the risk factors for EA., Methods: We used genome-wide association data to investigate the genetic architecture underlying GERD, BE, and EA. We applied a method to estimate the variance explained (array heritability, h(2)g) and the genetic correlation (rg) between GERD, BE, and EA by considering all single nucleotide polymorphisms (SNPs) simultaneously. We also estimated the polygenic overlap between GERD, BE, and EA using a prediction approach. All tests were two-sided, except in the case of variance-explained estimation where one-sided tests were used., Results: We estimated a statistically significant genetic variance explained for BE (h(2)g = 35%; standard error [SE] = 6%; one-sided P = 1 × 10(-9)) and for EA (h(2)g = 25 %; SE = 5%; one-sided P = 2 × 10(-7)). The genetic correlation between BE and EA was found to be high (rg = 1.0; SE = 0.37). We also estimated a statistically significant polygenic overlap between BE and EA (one-sided P = 1 × 10(-6)), which suggests, together with the high genetic correlation, that shared genes underlie the development of BE and EA. Conversely, no statistically significant results were obtained for GERD., Conclusions: We have demonstrated that risk to BE and EA is influenced by many germline genetic variants of small effect and that shared polygenic effects contribute to risk of these two diseases.

To mine possibly hidden causal single-nucleotide polymorphisms (SNPs) of melanoma, we investigated the association of SNPs in 76 M/G1 transition genes with melanoma risk using our published genome-wide association study (GWAS) data set with 1804 melanoma cases and 1026 cancer-free controls. We found multiple SNPs with P < 0.01 and performed validation studies for 18 putative functional SNPs in PSMB9 in two other GWAS data sets. Two SNPs (rs1351383 and rs2127675) were associated with melanoma risk in the GenoMEL data set (P = 0.013 and 0.004, respectively), but failed in validation using the Australian data set. Genotype-phenotype analysis revealed these two SNPs were significantly correlated with mRNA expression level of PSMB9. Further experiments revealed that SNP rs2071480, which is in high LD with rs1351383 and rs2127675, may have a weak effect on the promoter activity of PSMB9. Taken together, our data suggested that functional variants in PSMB9 may contribute to melanoma susceptibility., (© 2013 John Wiley & Sons A/S.)

We report the results of an association study of melanoma that is based on the genome-wide imputation of the genotypes of 1,353 cases and 3,566 controls of European origin conducted by the GenoMEL consortium. This revealed an association between several SNPs in intron 8 of the FTO gene, including rs16953002, which replicated using 12,313 cases and 55,667 controls of European ancestry from Europe, the USA and Australia (combined P = 3.6 × 10(-12), per-allele odds ratio for allele A = 1.16). In addition to identifying a new melanoma-susceptibility locus, this is to our knowledge the first study to identify and replicate an association with SNPs in FTO not related to body mass index (BMI). These SNPs are not in intron 1 (the BMI-related region) and exhibit no association with BMI. This suggests FTO's function may be broader than the existing paradigm that FTO variants influence multiple traits only through their associations with BMI and obesity.

Genome-wide association studies (GWASs) have mainly focused on top significant single nucleotide polymorphisms (SNPs), most of which did not have clear biological functions but were just surrogates for unknown causal variants. Studying SNPs with modest association and putative functions in biologically plausible pathways has become one complementary approach to GWASs. To unravel the key roles of mitogen-activated protein kinase (MAPK) pathways in cutaneous melanoma (CM) risk, we re-evaluated the associations between 47 818 SNPs in 280 MAPK genes and CM risk using our published GWAS dataset with 1804 CM cases and 1026 controls. We initially found 105 SNPs with P ≤ 0.001, more than expected by chance, 26 of which were predicted to be putatively functional SNPs. The risk associations with 16 SNPs around DUSP14 (rs1051849) and a previous reported melanoma locus MAFF/PLA2G6 (proxy SNP rs4608623) were replicated in the GenoMEL dataset (P < 0.01) but failed in the Australian dataset. Meta-analysis showed that rs1051849 in the 3' untranslated regions of DUSP14 was associated with a reduced risk of melanoma (odds ratio = 0.89, 95% confidence interval: 0.82-0.96, P = 0.003, false discovery rate = 0.056). Further genotype-phenotype correlation analysis using the 90 HapMap lymphoblastoid cell lines from Caucasians showed significant correlations between two SNPs (rs1051849 and rs4608623) and messenger RNA expression levels of DUSP14 and MAFF (P = 0.025 and P = 0.010, respectively). Gene-based tests also revealed significant SNPs were over-represented in MAFF, PLA2G6, DUSP14 and other 16 genes. Our results suggest that functional SNPs in MAPK pathways may contribute to CM risk. Further studies are warranted to validate our findings.

Purpose: Astigmatism is a common refractive error that reduces vision, where the curvature and refractive power of the cornea in one meridian are less than those of the perpendicular axis. It is a complex trait likely to be influenced by both genetic and environmental factors. Twin studies of astigmatism have found approximately 60% of phenotypic variance is explained by genetic factors. This study aimed to identify susceptibility loci for astigmatism., Methods: We performed a meta-analysis of seven genome-wide association studies that included 22,100 individuals of European descent, where astigmatism was defined as the number of diopters of cylinder prescription, using fixed effect inverse variance-weighted methods., Results: A susceptibility locus was identified with lead single nucleotide polymorphism rs3771395 on chromosome 2p13.3 (meta-analysis, P = 1.97 × 10(-7)) in the VAX2 gene. VAX2 plays an important role in the development of the dorsoventral axis of the eye. Animal studies have shown a gradient in astigmatism along the vertical plane, with corresponding changes in refraction, particularly in the ventral field., Conclusions: This finding advances the understanding of refractive error, and provides new potential pathways to be evaluated with regard to the development of astigmatism.

Barrett's esophagus is an increasingly common disease that is strongly associated with reflux of stomach acid and usually a hiatus hernia, and it strongly predisposes to esophageal adenocarcinoma (EAC), a tumor with a very poor prognosis. We report the first genome-wide association study on Barrett's esophagus, comprising 1,852 UK cases and 5,172 UK controls in the discovery stage and 5,986 cases and 12,825 controls in the replication stage. Variants at two loci were associated with disease risk: chromosome 6p21, rs9257809 (Pcombined=4.09×10(-9); odds ratio (OR)=1.21, 95% confidence interval (CI)=1.13-1.28), within the major histocompatibility complex locus, and chromosome 16q24, rs9936833 (Pcombined=2.74×10(-10); OR=1.14, 95% CI=1.10-1.19), for which the closest protein-coding gene is FOXF1, which is implicated in esophageal development and structure. We found evidence that many common variants of small effect contribute to genetic susceptibility to Barrett's esophagus and that SNP alleles predisposing to obesity also increase risk for Barrett's esophagus.

Recent Genome-Wide Association Studies (GWAS) have identified four low-penetrance ovarian cancer susceptibility loci. We hypothesized that further moderate- or low-penetrance variants exist among the subset of single-nucleotide polymorphisms (SNPs) not well tagged by the genotyping arrays used in the previous studies, which would account for some of the remaining risk. We therefore conducted a time- and cost-effective stage 1 GWAS on 342 invasive serous cases and 643 controls genotyped on pooled DNA using the high-density Illumina 1M-Duo array. We followed up 20 of the most significantly associated SNPs, which are not well tagged by the lower density arrays used by the published GWAS, and genotyping them on individual DNA. Most of the top 20 SNPs were clearly validated by individually genotyping the samples used in the pools. However, none of the 20 SNPs replicated when tested for association in a much larger stage 2 set of 4,651 cases and 6,966 controls from the Ovarian Cancer Association Consortium. Given that most of the top 20 SNPs from pooling were validated in the same samples by individual genotyping, the lack of replication is likely to be due to the relatively small sample size in our stage 1 GWAS rather than due to problems with the pooling approach. We conclude that there are unlikely to be any moderate or large effects on ovarian cancer risk untagged by less dense arrays. However, our study lacked power to make clear statements on the existence of hitherto untagged small-effect variants.

Purpose: Keratoconus is a progressive disorder of the cornea that can lead to severe visual impairment or blindness. Although several genomic regions have been linked to rare familial forms of keratoconus, no genes have yet been definitively identified for common forms of the disease., Methods: Two genome-wide association scans were undertaken in parallel. The first used pooled DNA from an Australian cohort, followed by typing of top-ranked single-nucleotide polymorphisms (SNPs) in individual DNA samples. The second was conducted in individually genotyped patients, and controls from the USA. Tag SNPs around the hepatocyte growth factor (HGF) gene were typed in three additional replication cohorts. Serum levels of HGF protein in normal individuals were assessed with ELISA and correlated with genotype., Results: The only SNP observed to be associated in both the pooled discovery and primary replication cohort was rs1014091, located upstream of the HGF gene. The nearby SNP rs3735520 was found to be associated in the individually typed discovery cohort (P = 6.1 × 10(-7)). Genotyping of tag SNPs around HGF revealed association at rs3735520 and rs17501108/rs1014091 in four of the five cohorts. Meta-analysis of all five datasets together yielded suggestive P values for rs3735520 (P = 9.9 × 10(-7)) and rs17501108 (P = 9.9 × 10(-5)). In addition, SNP rs3735520 was found to be associated with serum HGF level in normal individuals (P = 0.036)., Conclusions: Taken together, these results implicate genetic variation at the HGF locus with keratoconus susceptibility.

Purpose: Cell-based approaches were used to identify genetic markers predictive of patients' risk for poor response prior to chemotherapy., Experimental Design: We conducted genome-wide association studies (GWAS) to identify single-nucleotide polymorphisms (SNP) associated with cellular sensitivity to carboplatin through their effects on mRNA expression using International HapMap lymphoblastoid cell lines (LCL) and replicated them in additional LCLs. SNPs passing both stages of the cell-based study were tested for association with progression-free survival (PFS) in patients. Phase 1 validation was based on 377 ovarian cancer patients receiving at least four cycles of carboplatin and paclitaxel from the Australian Ovarian Cancer Study (AOCS). Positive associations were then assessed in phase 2 validation analysis of 1,326 patients from the Ovarian Cancer Association Consortium and The Cancer Genome Atlas., Results: In the initial GWAS, 342 SNPs were associated with carboplatin-induced cytotoxicity, of which 18 unique SNPs were retained after assessing their association with gene expression. One SNP (rs1649942) was replicated in an independent LCL set (Bonferroni adjusted P < 0.05). It was found to be significantly associated with decreased PFS in phase 1 AOCS patients (P(per-allele) = 2 × 10(-2)), with a stronger effect in the subset of women with optimally debulked tumors (P(per-allele) = 4 × 10(-3)). rs1649942 was also associated with poorer overall survival in women with optimally debulked tumors (P(per-allele) = 9 × 10(-3)). However, this SNP was not significant in phase 2 validation analysis with patients from numerous cohorts., Conclusion: This study shows the potential of cell-based, genome-wide approaches to identify germline predictors of treatment outcome and highlights the need for extensive validation in patients to assess their clinical effect., (©2011 AACR.)

Human iris patterns are highly variable. The origins of this variation are of interest in the study of iris-related eye diseases and forensics, as well as from an embryological developmental perspective, with regard to their possible relationship to fundamental processes of neurodevelopment. We have performed genome-wide association scans on four iris characteristics (crypt frequency, furrow contractions, presence of peripupillary pigmented ring, and number of nevi) in three Australian samples of European descent. Both the discovery (n = 2121) and replication (n = 499 and 73) samples showed evidence for association between (1) crypt frequency and variants in the axonal guidance gene SEMA3A (p = 6.6 × 10(-11)), (2) furrow contractions and variants within the cytoskeleton gene TRAF3IP1 (p = 2.3 × 10(-12)), and (3) the pigmented ring and variants in the well-known pigmentation gene SLC24A4 (p = 7.6 × 10(-21)). These replicated findings individually accounted for around 1.5%-3% of the variance for these iris characteristics. Because both SEMA3A and TRAFIP1 are implicated in pathways that control neurogenesis, neural migration, and synaptogenesis, we also examined the evidence of enhancement among such genes, finding enrichment for crypts and furrows. These findings suggest that genes involved in normal neuronal pattern development may also influence tissue structures in the human iris., (Copyright © 2011 The American Society of Human Genetics. Published by Elsevier Inc. All rights reserved.)

Genetic variation at the TERT-CLPTM1L locus at 5p15.33 is associated with susceptibility to several cancers, including epithelial ovarian cancer (EOC). We have carried out fine-mapping of this region in EOC which implicates an association with a single nucleotide polymorphism (SNP) within the TERT promoter. We demonstrate that the minor alleles at rs2736109, and at an additional TERT promoter SNP, rs2736108, are associated with decreased breast cancer risk, and that the combination of both SNPs substantially reduces TERT promoter activity.

There is increasing evidence that the microcirculation plays an important role in the pathogenesis of cardiovascular diseases. Changes in retinal vascular caliber reflect early microvascular disease and predict incident cardiovascular events. We performed a genome-wide association study to identify genetic variants associated with retinal vascular caliber. We analyzed data from four population-based discovery cohorts with 15,358 unrelated Caucasian individuals, who are members of the Cohort for Heart and Aging Research in Genomic Epidemiology (CHARGE) consortium, and replicated findings in four independent Caucasian cohorts (n  =  6,652). All participants had retinal photography and retinal arteriolar and venular caliber measured from computer software. In the discovery cohorts, 179 single nucleotide polymorphisms (SNP) spread across five loci were significantly associated (p<5.0×10(-8)) with retinal venular caliber, but none showed association with arteriolar caliber. Collectively, these five loci explain 1.0%-3.2% of the variation in retinal venular caliber. Four out of these five loci were confirmed in independent replication samples. In the combined analyses, the top SNPs at each locus were: rs2287921 (19q13; p  =  1.61×10(-25), within the RASIP1 locus), rs225717 (6q24; p = 1.25×10(-16), adjacent to the VTA1 and NMBR loci), rs10774625 (12q24; p  =  2.15×10(-13), in the region of ATXN2,SH2B3 and PTPN11 loci), and rs17421627 (5q14; p = 7.32×10(-16), adjacent to the MEF2C locus). In two independent samples, locus 12q24 was also associated with coronary heart disease and hypertension. Our population-based genome-wide association study demonstrates four novel loci associated with retinal venular caliber, an endophenotype of the microcirculation associated with clinical cardiovascular disease. These data provide further insights into the contribution and biological mechanisms of microcirculatory changes that underlie cardiovascular disease., Competing Interests: The authors have declared that no competing interests exist.

We conducted a genome-wide association study for primary open-angle glaucoma (POAG) in 1,263 affected individuals (cases) and 34,877 controls from Iceland. We identified a common sequence variant at 7q31 (rs4236601[A], odds ratio (OR) = 1.36, P = 5.0 × 10⁻¹⁰). We then replicated the association in sample sets of 2,175 POAG cases and 2,064 controls from Sweden, the UK and Australia (combined OR = 1.18, P = 0.0015) and in 299 POAG cases and 580 unaffected controls from Hong Kong and Shantou, China (combined OR = 5.42, P = 0.0021). The risk variant identified here is located close to CAV1 and CAV2, both of which are expressed in the trabecular meshwork and retinal ganglion cells that are involved in the pathogenesis of POAG.

We hypothesized that variants in genes expressed as a consequence of interactions between ovarian cancer cells and the host micro-environment could contribute to cancer susceptibility. We therefore used a two-stage approach to evaluate common single nucleotide polymorphisms (SNPs) in 173 genes involved in stromal epithelial interactions in the Ovarian Cancer Association Consortium (OCAC). In the discovery stage, cases with epithelial ovarian cancer (n=675) and controls (n=1,162) were genotyped at 1,536 SNPs using an Illumina GoldenGate assay. Based on Positive Predictive Value estimates, three SNPs-PODXL rs1013368, ITGA6 rs13027811, and MMP3 rs522616-were selected for replication using TaqMan genotyping in up to 3,059 serous invasive cases and 8,905 controls from 16 OCAC case-control studies. An additional 18 SNPs with Pper-allele<0.05 in the discovery stage were selected for replication in a subset of five OCAC studies (n=1,233 serous invasive cases; n=3,364 controls). The discovery stage associations in PODXL, ITGA6, and MMP3 were attenuated in the larger replication set (adj. Pper-allele>or=0.5). However genotypes at TERT rs7726159 were associated with ovarian cancer risk in the smaller, five-study replication study (Pper-allele=0.03). Combined analysis of the discovery and replication sets for this TERT SNP showed an increased risk of serous ovarian cancer among non-Hispanic whites [adj. ORper-allele 1.14 (1.04-1.24) p=0.003]. Our study adds to the growing evidence that, like the 8q24 locus, the telomerase reverse transcriptase locus at 5p15.33, is a general cancer susceptibility locus., Competing Interests: The authors have declared that no competing interests exist.

Background: Genome-wide association studies have found type 2 diabetes-associated variants in the HNF1B gene to exhibit reciprocal associations with prostate cancer risk. We aimed to identify whether these variants may have an effect on cancer risk in general versus a specific effect on prostate cancer only., Methodology/principal Findings: In a collaborative analysis, we collected data from GWAS of cancer phenotypes for the frequently reported variants of HNF1B, rs4430796 and rs7501939, which are in linkage disequilibrium (r(2) = 0.76, HapMap CEU). Overall, the analysis included 16 datasets on rs4430796 with 19,640 cancer cases and 21,929 controls; and 21 datasets on rs7501939 with 26,923 cases and 49,085 controls. Malignancies other than prostate cancer included colorectal, breast, lung and pancreatic cancers, and melanoma. Meta-analysis showed large between-dataset heterogeneity that was driven by different effects in prostate cancer and other cancers. The per-T2D-risk-allele odds ratios (95% confidence intervals) for rs4430796 were 0.79 (0.76, 0.83)] per G allele for prostate cancer (p<10(-15) for both); and 1.03 (0.99, 1.07) for all other cancers. Similarly for rs7501939 the per-T2D-risk-allele odds ratios (95% confidence intervals) were 0.80 (0.77, 0.83) per T allele for prostate cancer (p<10(-15) for both); and 1.00 (0.97, 1.04) for all other cancers. No malignancy other than prostate cancer had a nominally statistically significant association., Conclusions/significance: The examined HNF1B variants have a highly specific effect on prostate cancer risk with no apparent association with any of the other studied cancer types.

Central corneal thickness (CCT), one of the most highly heritable human traits (h(2) typically>0.9), is important for the diagnosis of glaucoma and a potential risk factor for glaucoma susceptibility. We conducted genome-wide association studies in five cohorts from Australia and the United Kingdom (total N = 5058). Three cohorts were based on individually genotyped twin collections, with the remaining two cohorts genotyped on pooled samples from singletons with extreme trait values. The pooled sample findings were validated by individual genotyping the pooled samples together with additional samples also within extreme quantiles. We describe methods for efficient combined analysis of the results from these different study designs. We have identified and replicated quantitative trait loci on chromosomes 13 and 16 for association with CCT. The locus on chromosome 13 (nearest gene FOXO1) had an overall meta-analysis p-value for all the individually genotyped samples of 4.6x10(-10). The locus on chromosome 16 was associated with CCT with p = 8.95x10(-11). The nearest gene to the associated chromosome 16 SNPs was ZNF469, a locus recently implicated in Brittle Cornea Syndrome (BCS), a very rare disorder characterized by abnormal thin corneas. Our findings suggest that in addition to rare variants in ZNF469 underlying CCT variation in BCS patients, more common variants near this gene may contribute to CCT variation in the general population., Competing Interests: The authors have declared that no competing interests exist.

Genome-wide association studies (GWAS) have now successfully identified important genetic variants associated with many human traits and diseases. The high cost of genotyping arrays in large data sets remains the major barrier to wider utilization of GWAS. We have developed a novel method in which whole blood from cases and controls, respectively, is pooled prior to DNA extraction for genotyping. We demonstrate proof of principle by clearly identifying the associated variants for eye color, age-related macular degeneration, and pseudoexfoliation syndrome in cohorts not previously studied. Blood pooling has the potential to reduce GWAS cost by several orders of magnitude and dramatically shorten gene discovery time. This method has profound implications for translation of modern genetic approaches to a multitude of diseases and traits yet to be analyzed by GWAS, and will enable developing nations to participate in GWAS.