GWAS meta-analysis of 16 790 patients with Barrett's oesophagus and oesophageal adenocarcinoma identifies 16 novel genetic risk loci and provides insights into disease aetiology beyond the single marker level.

OBJECTIVE: Oesophageal cancer (EC) is the sixth leading cause of cancer-related deaths. Oesophageal adenocarcinoma (EA), with Barrett's oesophagus (BE) as a precursor lesion, is the most prevalent EC subtype in the Western world. This study aims to contribute to better understand the genetic causes of BE/EA by leveraging genome wide association studies (GWAS), genetic correlation analyses and polygenic risk modelling. DESIGN: We combined data from previous GWAS with new cohorts, increasing the sample size to 16 790 BE/EA cases and 32 476 controls. We also carried out a transcriptome wide association study (TWAS) using expression data from disease-relevant tissues to identify BE/EA candidate genes. To investigate the relationship with reported BE/EA risk factors, a linkage disequilibrium score regression (LDSR) analysis was performed. BE/EA risk models were developed combining clinical/lifestyle risk factors with polygenic risk scores (PRS) derived from the GWAS meta-analysis. RESULTS: The GWAS meta-analysis identified 27 BE and/or EA risk loci, 11 of which were novel. The TWAS identified promising BE/EA candidate genes at seven GWAS loci and at five additional risk loci. The LDSR analysis led to the identification of novel genetic correlations and pointed to differences in BE and EA aetiology. Gastro-oesophageal reflux disease appeared to contribute stronger to the metaplastic BE transformation than to EA development. Finally, combining PRS with BE/EA risk factors improved the performance of the risk models. CONCLUSION: Our findings provide further insights into BE/EA aetiology and its relationship to risk factors. The results lay the foundation for future follow-up studies to identify underlying disease mechanisms and improving risk prediction.

Germline variation in the insulin-like growth factor pathway and risk of Barrett's esophagus and esophageal adenocarcinoma.

Genome-wide association studies (GWAS) of esophageal adenocarcinoma (EAC) and its precursor, Barrett's esophagus (BE), have uncovered significant genetic components of risk, but most heritability remains unexplained. Targeted assessment of genetic variation in biologically relevant pathways using novel analytical approaches may identify missed susceptibility signals. Central obesity, a key BE/EAC risk factor, is linked to systemic inflammation, altered hormonal signaling and insulin-like growth factor (IGF) axis dysfunction. Here, we assessed IGF-related genetic variation and risk of BE and EAC. Principal component analysis was employed to evaluate pathway-level and gene-level associations with BE/EAC, using genotypes for 270 single-nucleotide polymorphisms (SNPs) in or near 12 IGF-related genes, ascertained from 3295 BE cases, 2515 EAC cases and 3207 controls in the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON) GWAS. Gene-level signals were assessed using Multi-marker Analysis of GenoMic Annotation (MAGMA) and SNP summary statistics from BEACON and an expanded GWAS meta-analysis (6167 BE cases, 4112 EAC cases, 17 159 controls). Global variation in the IGF pathway was associated with risk of BE (P = 0.0015). Gene-level associations with BE were observed for GHR (growth hormone receptor; P = 0.00046, false discovery rate q = 0.0056) and IGF1R (IGF1 receptor; P = 0.0090, q = 0.0542). These gene-level signals remained significant at q < 0.1 when assessed using data from the largest available BE/EAC GWAS meta-analysis. No significant associations were observed for EAC. This study represents the most comprehensive evaluation to date of inherited genetic variation in the IGF pathway and BE/EAC risk, providing novel evidence that variation in two genes encoding cell-surface receptors, GHR and IGF1R, may influence risk of BE.

Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.

Sequencing the GRHL3 Coding Region Reveals Rare Truncating Mutations and a Common Susceptibility Variant for Nonsyndromic Cleft Palate.

Nonsyndromic cleft lip with/without cleft palate (nsCL/P) and nonsyndromic cleft palate only (nsCPO) are the most frequent subphenotypes of orofacial clefts. A common syndromic form of orofacial clefting is Van der Woude syndrome (VWS) where individuals have CL/P or CPO, often but not always associated with lower lip pits. Recently, ∼5% of VWS-affected individuals were identified with mutations in the grainy head-like 3 gene (GRHL3). To investigate GRHL3 in nonsyndromic clefting, we sequenced its coding region in 576 Europeans with nsCL/P and 96 with nsCPO. Most strikingly, nsCPO-affected individuals had a higher minor allele frequency for rs41268753 (0.099) than control subjects (0.049; p = 1.24 × 10(-2)). This association was replicated in nsCPO/control cohorts from Latvia, Yemen, and the UK (pcombined = 2.63 × 10(-5); ORallelic = 2.46 [95% CI 1.6-3.7]) and reached genome-wide significance in combination with imputed data from a GWAS in nsCPO triads (p = 2.73 × 10(-9)). Notably, rs41268753 is not associated with nsCL/P (p = 0.45). rs41268753 encodes the highly conserved p.Thr454Met (c.1361C>T) (GERP = 5.3), which prediction programs denote as deleterious, has a CADD score of 29.6, and increases protein binding capacity in silico. Sequencing also revealed four novel truncating GRHL3 mutations including two that were de novo in four families, where all nine individuals harboring mutations had nsCPO. This is important for genetic counseling: given that VWS is rare compared to nsCPO, our data suggest that dominant GRHL3 mutations are more likely to cause nonsyndromic than syndromic CPO. Thus, with rare dominant mutations and a common risk variant in the coding region, we have identified an important contribution for GRHL3 in nsCPO.

No Association Between Vitamin D Status and Risk of Barrett's Esophagus or Esophageal Adenocarcinoma: A Mendelian Randomization Study.

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.

Meta-analysis Reveals Genome-Wide Significance at 15q13 for Nonsyndromic Clefting of Both the Lip and the Palate, and Functional Analyses Implicate GREM1 As a Plausible Causative Gene.

Nonsyndromic orofacial clefts are common birth defects with multifactorial etiology. The most common type is cleft lip, which occurs with or without cleft palate (nsCLP and nsCLO, respectively). Although genetic components play an important role in nsCLP, the genetic factors that predispose to palate involvement are largely unknown. In this study, we carried out a meta-analysis on genetic and clinical data from three large cohorts and identified strong association between a region on chromosome 15q13 and nsCLP (P = 8.13 × 10(-14) for rs1258763; relative risk (RR): 1.46, 95% confidence interval (CI): 1.32-1.61)) but not nsCLO (P = 0.27; RR: 1.09 (0.94-1.27)). The 5 kb region of strongest association maps downstream of Gremlin-1 (GREM1), which encodes a secreted antagonist of the BMP4 pathway. We show during mouse embryogenesis, Grem1 is expressed in the developing lip and soft palate but not in the hard palate. This is consistent with genotype-phenotype correlations between rs1258763 and a specific nsCLP subphenotype, since a more than two-fold increase in risk was observed in patients displaying clefts of both the lip and soft palate but who had an intact hard palate (RR: 3.76, CI: 1.47-9.61, Pdiff<0.05). While we did not find lip or palate defects in Grem1-deficient mice, wild type embryonic palatal shelves developed divergent shapes when cultured in the presence of ectopic Grem1 protein (P = 0.0014). The present study identified a non-coding region at 15q13 as the second, genome-wide significant locus specific for nsCLP, after 13q31. Moreover, our data suggest that the closely located GREM1 gene contributes to a rare clinical nsCLP entity. This entity specifically involves abnormalities of the lip and soft palate, which develop at different time-points and in separate anatomical regions.

Compound heterozygosity of low-frequency promoter deletions and rare loss-of-function mutations in TXNL4A causes Burn-McKeown syndrome.

Mutations in components of the major spliceosome have been described in disorders with craniofacial anomalies, e.g., Nager syndrome and mandibulofacial dysostosis type Guion-Almeida. The U5 spliceosomal complex of eight highly conserved proteins is critical for pre-mRNA splicing. We identified biallelic mutations in TXNL4A, a member of this complex, in individuals with Burn-McKeown syndrome (BMKS). This rare condition is characterized by bilateral choanal atresia, hearing loss, cleft lip and/or palate, and other craniofacial dysmorphisms. Mutations were found in 9 of 11 affected families. In 8 families, affected individuals carried a rare loss-of-function mutation (nonsense, frameshift, or microdeletion) on one allele and a low-frequency 34 bp deletion (allele frequency 0.76%) in the core promoter region on the other allele. In a single highly consanguineous family, formerly diagnosed as oculo-oto-facial dysplasia, the four affected individuals were homozygous for a 34 bp promoter deletion, which differed from the promoter deletion in the other families. Reporter gene and in vivo assays showed that the promoter deletions led to reduced expression of TXNL4A. Depletion of TXNL4A (Dib1) in yeast demonstrated reduced assembly of the tri-snRNP complex. Our results indicate that BMKS is an autosomal-recessive condition, which is frequently caused by compound heterozygosity of low-frequency promoter deletions in combination with very rare loss-of-function mutations.

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

BACKGROUND: Oesophageal adenocarcinoma represents one of the fastest rising cancers in high-income countries. Barrett's oesophagus is the premalignant precursor of oesophageal adenocarcinoma. However, only a few patients with Barrett's oesophagus develop adenocarcinoma, which complicates clinical management in the absence of valid predictors. Within an international consortium investigating the genetics of Barrett's oesophagus and oesophageal adenocarcinoma, we aimed to identify novel genetic risk variants for the development of Barrett's oesophagus and oesophageal adenocarcinoma. METHODS: We did a meta-analysis of all genome-wide association studies of Barrett's 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 Barrett's 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 Barrett's 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 Barrett's 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 Barrett's 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 Barrett's oesophagus to oesophageal adenocarcinoma. Fine-mapping and functional studies of new risk loci could lead to identification of key molecules in the development of Barrett's 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.

Replication analysis of 15 susceptibility loci for nonsyndromic cleft lip with or without cleft palate in an italian population.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is one of the most common congenital malformations in humans. Its average global incidence is 1.7 per 1000 live births, with wide variation according to geographical location and ethnicity. Its etiology involves both genetic and environmental factors. The aim of the present study was to confirm genetic association of a selection of 15 candidate nsCL/P loci using an independent sample of the Italian population. METHODS: At least one single-nucleotide polymorphism (SNP) for each locus was genotyped in 380 nuclear trios. RESULTS: Transmission disequilibrium analysis revealed significant associations for three variants at two loci (8q24 and 1p22). Two SNPs at 8q24 showed the strongest level of association, the rs987525 (PTDT = 6.81 × 10(-6) ; homozygous relative risk = 3.60 [95% confidence interval, 2.12-6.13]), and the rs17241253 (PTDT = 1.03 × 10(-5) ; homozygous relative risk = 3.75 [95% confidence interval, 2.10-6.67]). Four additional loci (at 1q32, 3q12, 8q21, and 10q25) achieved nominally significant p-values. Two SNPs at 1p36 achieved p-values of < 0.1. The present data suggest that 6 of the 15 analyzed nsCL/P risk loci contribute significantly to nsCL/P risk in the Italian population. These include the 1p22 locus, which previous research has implicated predominantly in Asian populations. CONCLUSION: Different loci, including 8q24 and 1p22 have been found associated with nsCL/P in multiple populations. Further efforts are needed to identify causative variants and transfer knowledge to clinical application, such as personal genetic risk assessment.

Further evidence for deletions in 7p14.1 contributing to nonsyndromic cleft lip with or without cleft palate.

BACKGROUND: Nonsyndromic cleft with or without cleft palate (nsCL/P) is a common birth defect. Although genome-wide association studies (GWAS) have identified numerous risk variants, a considerable fraction of the genetic heritability remains unknown. The aim of the present study was to replicate a previous finding that de novo deletions in a 62 kb region of chromosome 7p14 are a risk factor for nsCL/P, using an independent cohort. METHODS: Data from a published case-control GWAS cohort of 399 patients and 1318 controls were used. Copy number variant (CNV) detection in the 62 kb candidate region of 7p14 was performed using QuantiSNP. Putative CNVs in probands were verified and validated by quantitative polymerase chain reaction. Segregation analyses were performed in family members for whom DNA was available. RESULTS: Within the 62 kb candidate region, a deletion of 7.4 kb showed association with nsCL/P (13/387 cases, 20/1300 controls, plowest = 0.024, odds ratio = 2.22). In all families with a sporadic case (n = 3), the deletion occurred de novo. In multiplex families, both incomplete segregation and incomplete penetrance were observed. CONCLUSION: The present data support the hypothesis that deletions at 7p14 are a common risk factor for nsCL/P. Genome-wide CNV analyses in nsCL/P cohorts are warranted to explore the functional relevance of these deletions and their contribution to nsCL/P, and to determine exact breakpoints. Birth Defects Research (Part A) 106:767-772, 2016. © 2016 Wiley Periodicals, Inc.

Investigation of the involvement of MIR185 and its target genes in the development of schizophrenia.

BACKGROUND: Schizophrenia is a complex neuropsychiatric disorder of unclear etiology. The strongest known genetic risk factor is the 22q11.2 microdeletion. Research has yet to confirm which genes within the deletion region are implicated in schizophrenia. The minimal 1.5 megabase deletion contains MIR185, which encodes microRNA 185. METHODS: We determined miR-185 expression in embryonic and adult mouse brains. Common and rare variants at this locus were then investigated using a human genetics approach. First, we performed gene-based analyses for MIR185 common variants and target genes using Psychiatric Genomics Consortium genome-wide association data. Second, MIR185 was resequenced in German patients (n = 1000) and controls (n = 500). We followed up promising variants by genotyping an additional European sample (patients, n = 3598; controls, n = 4082). RESULTS: In situ hybridization in mice revealed miR-185 expression in brain regions implicated in schizophrenia. Gene-based tests revealed association between common variants in 3 MIR185 target genes (ATAT1, SH3PXD2A, NTRK3) and schizophrenia. Further analyses in mice revealed overlapping expression patterns for these target genes and miR-185. Resequencing identified 2 rare patient-specific novel variants flanking MIR185. However, follow-up genotyping provided no further evidence of their involvement in schizophrenia. LIMITATIONS: Power to detect rare variant associations was limited. CONCLUSION: Human genetic analyses generated no evidence of the involvement of MIR185 in schizophrenia. However, the expression patterns of miR-185 and its target genes in mice, and the genetic association results for the 3 target genes, suggest that further research into the involvement of miR-185 and its downstream pathways in schizophrenia is warranted.

Nonsyndromic cleft lip with or without cleft palate in arab populations: genetic analysis of 15 risk loci in a novel case-control sample recruited in Yemen.

BACKGROUND: Nonsyndromic orofacial clefting (nsOFC) is among the most common of all congenital disorders and has a genetically complex etiology. Based on embryological and epidemiological data, the phenotype can be differentiated into nonsyndromic cleft lip with or without cleft palate (nsCL/P) and nonsyndromic cleft palate only, with nsCL/P being the most frequent form. Recent genetic research, predominantly performed in populations from Europe and Asia, has identified numerous genetic susceptibility loci for nsCL/P. As only few data are available concerning genetic susceptibility to nsCL/P in Arab populations, we investigated a newly recruited nsOFC sample from Yemen. METHODS: For each of the 15 currently known nsCL/P risk loci, the top single-nucleotide polymorphism (plus nine back-up variants) were genotyped in 242 nsCL/P cases and 420 healthy controls. RESULTS: Single-marker association analysis revealed significant associations for four loci (8q24, 9q22, 10q25, 13q31). The strongest association was for the European high risk locus at 8q24 (Pcorrected = 5.09 × 10(-4) ; heterozygous odds ratio = 1.74 (1.22-2.47), homozygous odds ratio = 2.47 (1.55-3.93). Five additional loci (1q32.2, 3q12, 8q21, 17q22, 20q12) showed nominal significance that did not withstand correction for multiple testing. Although the six remaining loci (1p22, 1p36, 2p21, 3p11, 15q22, 17p13) failed to reach nominal significance, the risk alleles were in the same direction as in the discovery studies. CONCLUSION: The results suggest that four of the 15 analyzed nsCL/P risk loci which were identified in European and Asian ethnicities significantly confer risk for nsCL/P in Arab populations.

Evaluating eight newly identified susceptibility loci for nonsyndromic cleft lip with or without cleft palate in a Mesoamerican population.

BACKGROUND: Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is among the most frequently occurring congenital malformations worldwide. The number of genetic loci identified as being involved in NSCL/P etiology was recently increased by a large genome-wide meta-analysis of European and Asian samples. This meta-analysis confirmed all six previously recognized genetic susceptibility loci and identified six novel ones. METHODS: To investigate which of these 12 loci contribute to NSCL/P risk in an independent sample of distinct ethnicity, we performed a case-control association analysis in a sample of the Mesoamerican population. A total of 153 individuals with NSCL/P (cases) and 337 unaffected controls were included. Top single-nucleotide polymorphisms (SNPs) at 8 of the 12 loci (1p22.1, 1p36, 2p21, 3p11.1, 8q21.3, 13q31.1, 15q22, and 20q12) were analyzed using mass spectroscopy and restriction-length-fragment polymorphism analyses. In a previous study, we had analyzed the remaining four NSCL/P susceptibility regions (IRF6, 8q24, 10q25, and 17q22) in the same sample. RESULTS: Single-marker association analyses applying allelic, dominant, and recessive models revealed nominal significant associations for four of the eight loci, with two additional loci showing at least a trend of association in the hypothesized direction. CONCLUSION: In combination with results from our previous study using the same sample, our data suggest that the majority of the known NSCL/P susceptibility regions identified to date also confer risk for this malformation in the Mesoamerican population. Birth Defects Research (Part A) 100:43-47, 2014. © 2013 Wiley Periodicals, Inc.

Nonsyndromic cleft lip with or without cleft palate: Increased burden of rare variants within Gremlin-1, a component of the bone morphogenetic protein 4 pathway.

BACKGROUND: The genes Gremlin-1 (GREM1) and Noggin (NOG) are components of the bone morphogenetic protein 4 pathway, which has been implicated in craniofacial development. Both genes map to recently identified susceptibility loci (chromosomal region 15q13, 17q22) for nonsyndromic cleft lip with or without cleft palate (nsCL/P). The aim of the present study was to determine whether rare variants in either gene are implicated in nsCL/P etiology. METHODS: The complete coding regions, untranslated regions, and splice sites of GREM1 and NOG were sequenced in 96 nsCL/P patients and 96 controls of Central European ethnicity. Three burden and four nonburden tests were performed. Statistically significant results were followed up in a second case-control sample (n = 96, respectively). For rare variants observed in cases, segregation analyses were performed. RESULTS: In NOG, four rare sequence variants (minor allele frequency < 1%) were identified. Here, burden and nonburden analyses generated nonsignificant results. In GREM1, 33 variants were identified, 15 of which were rare. Of these, five were novel. Significant p-values were generated in three nonburden analyses. Segregation analyses revealed incomplete penetrance for all variants investigated. CONCLUSION: Our study did not provide support for NOG being the causal gene at 17q22. However, the observation of a significant excess of rare variants in GREM1 supports the hypothesis that this is the causal gene at chr. 15q13. Because no single causal variant was identified, future sequencing analyses of GREM1 should involve larger samples and the investigation of regulatory elements.

Analysis of susceptibility loci for nonsyndromic orofacial clefting in a European trio sample.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P), the most common type of orofacial clefting, is one of the most frequent congenital defects. Based on epidemiological data, NSCL/P can be distinguished from nonsyndromic cleft palate only (NSCPO). Both phenotypes have a complex etiology and environmental and genetic factors are involved in their development. To date, genome-wide association studies have identified 12 genetic factors that increase the risk for NSCL/P in Europeans. Six of them have been independently replicated in samples derived from the same population. The aim of the present study was to replicate the remaining six NSCL/P risk loci in chromosomal regions 1p22.1, 1p36, 3p11.1, 8q21.3, 15q22.2, and 20q12 in a family-based sample of European descent. Each of the top-associated SNPs (single nucleotide polymorphisms) was genotyped in 343 NSCL/P and 266 NSCPO nuclear trios. Single-marker association analysis in the NSCL/P sample showed a significant association with SNP rs742071 (1p36, Pcorrected = 3.74 × 10(-3) ), which is located in the intronic region of PAX7, a gene known to be functionally implicated in craniofacial development. Two additional loci, 1p22.1 and 20q12, were nominally significant, but did not withstand correction for multiple testing. There was no evidence that the NSCL/P risk alleles contribute to the etiology of NSCPO, further supporting that these two subtypes of orofacial clefting are primarily etiologically distinct.

eQTL Set-Based Association Analysis Identifies Novel Susceptibility Loci for Barrett Esophagus and Esophageal Adenocarcinoma.

BACKGROUND: Over 20 susceptibility single-nucleotide polymorphisms (SNP) have been identified for esophageal adenocarcinoma (EAC) and its precursor, Barrett esophagus (BE), explaining a small portion of heritability. METHODS: Using genetic data from 4,323 BE and 4,116 EAC patients aggregated by international consortia including the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON), we conducted a comprehensive transcriptome-wide association study (TWAS) for BE/EAC, leveraging Genotype Tissue Expression (GTEx) gene-expression data from six tissue types of plausible relevance to EAC etiology: mucosa and muscularis from the esophagus, gastroesophageal (GE) junction, stomach, whole blood, and visceral adipose. Two analytical approaches were taken: standard TWAS using the predicted gene expression from local expression quantitative trait loci (eQTL), and set-based SKAT association using selected eQTLs that predict the gene expression. RESULTS: Although the standard approach did not identify significant signals, the eQTL set-based approach identified eight novel associations, three of which were validated in independent external data (eQTL SNP sets for EXOC3, ZNF641, and HSP90AA1). CONCLUSIONS: This study identified novel genetic susceptibility loci for EAC and BE using an eQTL set-based genetic association approach. IMPACT: This study expanded the pool of genetic susceptibility loci for EAC and BE, suggesting the potential of the eQTL set-based genetic association approach as an alternative method for TWAS analysis.

Shared Genetic Etiology of Obesity-Related Traits and Barrett's Esophagus/Adenocarcinoma: Insights from Genome-Wide Association Studies.

BACKGROUND: Obesity is a major risk factor for esophageal adenocarcinoma (EA) and its precursor Barrett's esophagus (BE). Research suggests that individuals with high genetic risk to obesity have a higher BE/EA risk. To facilitate understanding of biological factors that lead to progression from BE to EA, the present study investigated the shared genetic background of BE/EA and obesity-related traits. METHODS: Cross-trait linkage disequilibrium score regression was applied to summary statistics from genome-wide association meta-analyses on BE/EA and on obesity traits. Body mass index (BMI) was used as a proxy for general obesity, and waist-to-hip ratio (WHR) for abdominal obesity. For single marker analyses, all genome-wide significant risk alleles for BMI and WHR were compared with summary statistics of the BE/EA meta-analyses. RESULTS: Sex-combined analyses revealed a significant genetic correlation between BMI and BE/EA (rg = 0.13, P = 2 × 10-04) and a rg of 0.12 between WHR and BE/EA (P = 1 × 10-02). Sex-specific analyses revealed a pronounced genetic correlation between BMI and EA in females (rg = 0.17, P = 1.2 × 10-03), and WHR and EA in males (rg = 0.18, P = 1.51 × 10-02). On the single marker level, significant enrichment of concordant effects was observed for BMI and BE/EA risk variants (P = 8.45 × 10-03) and WHR and BE/EA risk variants (P = 2 × 10-02). CONCLUSIONS: Our study provides evidence for sex-specific genetic correlations that might reflect specific biological mecha-nisms. The data demonstrate that shared genetic factors are particularly relevant in progression from BE to EA. IMPACT: Our study quantifies the genetic correlation between BE/EA and obesity. Further research is now warranted to elucidate these effects and to understand the shared pathophysiology.

Author Correction: Gastroesophageal reflux GWAS identifies risk loci that also associate with subsequent severe esophageal diseases.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Gastroesophageal reflux GWAS identifies risk loci that also associate with subsequent severe esophageal diseases.

Gastroesophageal reflux disease (GERD) is caused by gastric acid entering the esophagus. GERD has high prevalence and is the major risk factor for Barrett's esophagus (BE) and esophageal adenocarcinoma (EA). We conduct a large GERD GWAS meta-analysis (80,265 cases, 305,011 controls), identifying 25 independent genome-wide significant loci for GERD. Several of the implicated genes are existing or putative drug targets. Loci discovery is greatest with a broad GERD definition (including cases defined by self-report or medication data). Further, 91% of the GERD risk-increasing alleles also increase BE and/or EA risk, greatly expanding gene discovery for these traits. Our results map genes for GERD and related traits and uncover potential new drug targets for these conditions.

Identification of loci of functional relevance to Barrett's esophagus and esophageal adenocarcinoma: Cross-referencing of expression quantitative trait loci data from disease-relevant tissues with genetic association data.

Esophageal adenocarcinoma (EA) and its precancerous condition Barrett's esophagus (BE) are multifactorial diseases with rising prevalence rates in Western populations. A recent meta-analysis of genome-wide association studies (GWAS) data identified 14 BE/EA risk loci located in non-coding genomic regions. Knowledge about the impact of non-coding variation on disease pathology is incomplete and needs further investigation. The aim of the present study was (i) to identify candidate genes of functional relevance to BE/EA at known risk loci and (ii) to find novel risk loci among the suggestively associated variants through the integration of expression quantitative trait loci (eQTL) and genetic association data. eQTL data from two BE/EA-relevant tissues (esophageal mucosa and gastroesophageal junction) generated within the context of the GTEx project were cross-referenced with the GWAS meta-analysis data. Variants representing an eQTL in at least one of the two tissues were categorized into genome-wide significant loci (P < 5×10-8) and novel candidate loci (5×10-8 ≤ P ≤ 5×10-5). To follow up these novel candidate loci, a genetic association study was performed in a replication cohort comprising 1,993 cases and 967 controls followed by a combined analysis with the GWAS meta-analysis data. The cross-referencing of eQTL and genetic data yielded 2,180 variants that represented 25 loci. Among the previously reported genome-wide significant loci, 22 eQTLs were identified in esophageal mucosa and/or gastroesophageal junction tissue. The regulated genes, most of which have not been linked to BE/EA etiology so far, included C2orf43/LDAH, ZFP57, and SLC9A3. Among the novel candidate loci, replication was achieved for two variants (rs7754014, Pcombined = 3.16×10-7 and rs1540, Pcombined = 4.16×10-6) which represent eQTLs for CFDP1 and SLC22A3, respectively. In summary, the present approach identified candidate genes whose expression was regulated by risk variants in disease-relevant tissues. These findings may facilitate the elucidation of BE/EA pathophysiology.