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.

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.

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.

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.

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.

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.

Evidence for PTGER4, PSCA, and MBOAT7 as risk genes for gastric cancer on the genome and transcriptome level.

Genetic associations between variants on chromosome 5p13 and 8q24 and gastric cancer (GC) have been previously reported in the Asian population. We aimed to replicate these findings and to characterize the associations at the genome and transcriptome level. We performed a fine-mapping association study in 1926 GC patients and 2012 controls of European descent using high dense SNP marker sets on both chromosomal regions. Next, we performed expression quantitative trait locus (eQTL) analyses using gastric transcriptome data from 143 individuals focusing on the GC associated variants. On chromosome 5p13 the strongest association was observed at rs6872282 (P = 2.53 × 10-04 ) and on chromosome 8q24 at rs2585176 (P = 1.09 × 10-09 ). On chromosome 5p13 we found cis-eQTL effects with an upregulation of PTGER4 expression in GC risk allele carrier (P = 9.27 × 10-11 ). On chromosome 8q24 we observed cis-eQTL effects with an upregulation of PSCA expression in GC risk allele carrier (P = 2.17 × 10-47 ). In addition, we found trans-eQTL effects for the same variants on 8q24 with a downregulation of MBOAT7 expression in GC risk allele carrier (P = 3.11 × 10-09 ). In summary, we confirmed and refined the previously reported GC associations at both chromosomal regions. Our data point to shared etiological factors between Asians and Europeans. Furthermore, our data imply an upregulated expression of PTGER4 and PSCA as well as a downregulated expression of MBOAT7 in gastric tissue as risk-conferring GC pathomechanisms.

Nonsyndromic cleft palate: An association study at GWAS candidate loci in a multiethnic sample.

BACKGROUND: Nonsyndromic cleft palate only (nsCPO) is a common and multifactorial form of orofacial clefting. In contrast to successes achieved for the other common form of orofacial clefting, that is, nonsyndromic cleft lip with/without cleft palate (nsCL/P), genome wide association studies (GWAS) of nsCPO have identified only one genome wide significant locus. Aim of the present study was to investigate whether common variants contribute to nsCPO and, if so, to identify novel risk loci. METHODS: We genotyped 33 SNPs at 27 candidate loci from 2 previously published nsCPO GWAS in an independent multiethnic sample. It included: (i) a family-based sample of European ancestry (n = 212); and (ii) two case/control samples of Central European (n = 94/339) and Arabian ancestry (n = 38/231), respectively. A separate association analysis was performed for each genotyped dataset, and meta-analyses were performed. RESULTS: After association analysis and meta-analyses, none of the 33 SNPs showed genome-wide significance. Two variants showed nominally significant association in the imputed GWAS dataset and exhibited a further decrease in p-value in a European and an overall meta-analysis including imputed GWAS data, respectively (rs395572: PMetaEU = 3.16 × 10-4 ; rs6809420: PMetaAll = 2.80 × 10-4 ). CONCLUSION: Our findings suggest that there is a limited contribution of common variants to nsCPO. However, the individual effect sizes might be too small for detection of further associations in the present sample sizes. Rare variants may play a more substantial role in nsCPO than in nsCL/P, for which GWAS of smaller sample sizes have identified genome-wide significant loci. Whole-exome/genome sequencing studies of nsCPO are now warranted.

Genetic estimators of DNA methylation provide insights into the molecular basis of polygenic traits.

The large biological distance between genetic risk loci and their mechanistic consequences in the tissue of interest limits the ability to establish functionality of susceptibility variants for genetically complex traits. Such a biological gap may be reduced through the systematic study of molecular mediators of genomic action, such as epigenetic modification. Here, we report the identification of robust genetic estimators of whole-blood CpG methylation, which can serve as intermediate molecular traits amenable to association testing with other genetically complex traits. We describe the relationship between these estimators and gene expression, demonstrate their genome-wide applicability to association testing even in the absence of individual genotypic data, and show that these estimators powerfully identify methylation-related genomic loci associated with polygenic traits and common diseases, such as schizophrenia. The use of genetic estimators for blood DNA methylation, which are made publically available, can serve as a valuable tool for the identification of epigenetic underpinnings of complex traits.

Investigation of dominant and recessive inheritance models in genome-wide association studies data of nonsyndromic cleft lip with or without cleft palate.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is one of the most common congenital malformation worldwide, and its etiology involves both genetic and environmental factors. Recent genome-wide and targeted genetic studies of nsCL/P have identified numerous genetic risk loci, under the hypothesis of a multiplicative mode of inheritance. The present study investigated whether novel nsCL/P risk loci could be identified by analyzing dominant/recessive genetic effects in single nucleotide polymorphism (SNP) data from genome-wide association studies. For this purpose, a genome-wide investigation of dominant/recessive common SNP effects was performed in our previously published meta-analysis data set. Twenty-four loci were identified as candidate regions. In a subsequent association analysis in an independent study cohort of 224 nsCL/P patients and 986 controls of European descent, none of the loci could be replicated. Therefore, our strategy of identifying novel loci by applying different genetic models did not yield any novel findings, suggesting that recessive/dominant common variation only make a limited contribution to nsCL/P in Europeans. However, we cannot rule out that such effects are present at some of the loci that have previously been identified, or are present in different populations.

Further evidence for genetic variation at the serotonin transporter gene SLC6A4 contributing toward anxiety.

OBJECTIVES: Social anxiety disorder (SAD) is a common and heritable psychiatric disorder. However, genetic studies in SAD are rare and only a few candidate genes have been implicated so far. In the present study, we investigated whether single-nucleotide polymorphisms (SNPs) associated with other psychiatric disorders also contribute toward the development of SAD and followed up variants associated with SAD on the phenotypic level. PATIENTS AND METHODS: We genotyped a total of 24 SNPs in a German sample of 321 SAD patients and 804 controls. We carried out single-marker analyses as well as quantitative association analyses of SAD severity and harm avoidance. RESULTS: None of the variants investigated showed an association with SAD in our case-control sample after Bonferroni correction. Two SNPs reached nominal significance (rs818702, P=0.032; rs140701, P=0.048). Of these, only rs140701 showed an association in the same allelic direction as reported previously. This SNP is located within the serotonin transporter gene SLC6A4, which is the primary target of selective-serotonin reuptake inhibitors used for the treatment of depressive and anxiety disorders. The quantitative association analysis of all cases with available data on symptom severity showed four SNPs with a nominal significant association. Among these SNPs, rs10994359 showed the strongest association (P=0.001) and was located near the ANK3 gene. In addition, rs10994359 was nominally associated with harm avoidance scores (P=0.001). CONCLUSION: Our results provide further evidence for an involvement of the serotonin transporter gene SLC6A4 in the etiology of anxiety-related traits. Furthermore, our study implicates that genetic variation at the genome-wide associated bipolar disorder locus ANK3 might influence anxiety-related personality traits.

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.

Nonsyndromic cleft lip with or without cleft palate and cancer: Evaluation of a possible common genetic background through the analysis of GWAS data.

Previous research suggests a genetic overlap between nonsyndromic cleft lip with or without cleft palate (NSCL/P) and cancer. The aim of the present study was to identify common genetic risk loci for NSCL/P and cancer entities that have been reported to co-occur with orofacial clefting. This was achieved through the investigation of large genome-wide association study datasets. Investigations of 12 NSCL/P single nucleotide polymorphisms (SNPs) in 32 cancer datasets, and 204 cancer SNPs in two NSCL/P datasets, were performed. The SNPs rs13041247 (20q12) and rs6457327 (6p21.33) showed suggestive evidence for an association with both NSCL/P and a specific cancer entity. These loci harbor genes of biological relevance to oncogenesis (MAFB and OCT4, respectively). This study is the first to characterize possible pleiotropic risk loci for NSCL/P and cancer in a systematic manner. The data represent a starting point for future research by identifying a genetic link between NSCL/P and cancer.

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.

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.

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.

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.

The Barrett-associated variants at GDF7 and TBX5 also increase esophageal adenocarcinoma risk.

Barrett's esophagus (BE) and esophageal adenocarcinoma (EAC) represent two stages within the esophagitis-metaplasia-dysplasia-adenocarcinoma sequence. Previously genetic risk factors have been identified that confer risk to BE and EAC development. However, to which extent the genetic variants confer risk to different stages of the BE/EAC sequence remains mainly unknown. In this study we analyzed three most recently identified BE variants at the genes GDF7 (rs3072), TBX5 (rs2701108), and ALDH1A2 (rs3784262) separately in BE and EAC samples in order to determine their risk effects during BE/EAC sequence. Our data show that rs3072 at GDF7 and rs2701108 at TBX5 are also associated with EAC and conclude that both loci confer disease risk also at later stages of the BE/EAC sequence. In contrast, rs3784262 at ALDH1A2 was highly significantly associated with BE, but showed no association with EAC. Our data do not provide evidence that the ALDH1A2 locus confers equal risk in early and late stages of BE/EAC sequence.