Integrative genomic analyses of European intrahepatic cholangiocarcinoma: Novel ROS1 fusion gene and PBX1 as prognostic marker.

BACKGROUND: Cholangiocarcinoma (CCA) is a fatal cancer of the bile duct with a poor prognosis owing to limited therapeutic options. The incidence of intrahepatic CCA (iCCA) is increasing worldwide, and its molecular basis is emerging. Environmental factors may contribute to regional differences in the mutation spectrum of European patients with iCCA, which are underrepresented in systematic genomic and transcriptomic studies of the disease. METHODS: We describe an integrated whole-exome sequencing and transcriptomic study of 37 iCCAs patients in Germany. RESULTS: We observed as most frequently mutated genes ARID1A (14%), IDH1, BAP1, TP53, KRAS, and ATM in 8% of patients. We identified FGFR2::BICC1 fusions in two tumours, and FGFR2::KCTD1 and TMEM106B::ROS1 as novel fusions with potential therapeutic implications in iCCA and confirmed oncogenic properties of TMEM106B::ROS1 in vitro. Using a data integration framework, we identified PBX1 as a novel central regulatory gene in iCCA. We performed extended screening by targeted sequencing of an additional 40 CCAs. In the joint analysis, IDH1 (13%), BAP1 (10%), TP53 (9%), KRAS (7%), ARID1A (7%), NF1 (5%), and ATM (5%) were the most frequently mutated genes, and we found PBX1 to show copy gain in 20% of the tumours. According to other studies, amplifications of PBX1 tend to occur in European iCCAs in contrast to liver fluke-associated Asian iCCAs. CONCLUSIONS: By analyzing an additional European cohort of iCCA patients, we found that PBX1 protein expression was a marker of poor prognosis. Overall, our findings provide insight into key molecular alterations in iCCA, reveal new targetable fusion genes, and suggest that PBX1 is a novel modulator of this disease.

The genetic regulation of the gastric transcriptome is associated with metabolic and obesity-related traits and diseases.

Tissue-specific gene expression and gene regulation lead to a better understanding of tissue-specific physiology and pathophysiology. We analyzed the transcriptome and genetic regulatory profiles of two distinct gastric sites, corpus and antrum, to identify tissue-specific gene expression and its regulation. Gastric corpus and antrum mucosa biopsies were collected during routine gastroscopies from up to 431 healthy individuals. We obtained genotype and transcriptome data and performed transcriptome profiling and expression quantitative trait locus (eQTL) studies. We further used data from genome-wide association studies (GWAS) of various diseases and traits to partition their heritability and to perform transcriptome-wide association studies (TWAS). The transcriptome data from corpus and antral mucosa highlights the heterogeneity of gene expression in the stomach. We identified enriched pathways revealing distinct and common physiological processes in gastric corpus and antrum. Furthermore, we found an enrichment of the single nucleotide polymorphism (SNP)-based heritability of metabolic, obesity-related, and cardiovascular traits and diseases by considering corpus- and antrum-specifically expressed genes. Particularly, we could prioritize gastric-specific candidate genes for multiple metabolic traits, like NQO1 which is involved in glucose metabolism, MUC1 which contributes to purine and protein metabolism or RAB27B being a regulator of weight and body composition. Our findings show that gastric corpus and antrum vary in their transcriptome and genetic regulatory profiles indicating physiological differences which are mostly related to digestion and epithelial protection. Moreover, our findings demonstrate that the genetic regulation of the gastric transcriptome is linked to biological mechanisms associated with metabolic, obesity-related, and cardiovascular traits and diseases. NEW & NOTEWORTHY We analyzed the transcriptomes and genetic regulatory profiles of gastric corpus and for the first time also of antrum mucosa in 431 healthy individuals. Through tissue-specific gene expression and eQTL analyses, we uncovered unique and common physiological processes across both primary gastric sites. Notably, our findings reveal that stomach-specific eQTLs are enriched in loci associated with metabolic traits and diseases, highlighting the pivotal role of gene expression regulation in gastric physiology and potential pathophysiology.

Pathogen-specific innate immune response patterns are distinctly affected by genetic diversity.

Innate immune responses vary by pathogen and host genetics. We analyze quantitative trait loci (eQTLs) and transcriptomes of monocytes from 215 individuals stimulated by fungal, Gram-negative or Gram-positive bacterial pathogens. We identify conserved monocyte responses to bacterial pathogens and a distinct antifungal response. These include 745 response eQTLs (reQTLs) and corresponding genes with pathogen-specific effects, which we find first in samples of male donors and subsequently confirm for selected reQTLs in females. reQTLs affect predominantly upregulated genes that regulate immune response via e.g., NOD-like, C-type lectin, Toll-like and complement receptor-signaling pathways. Hence, reQTLs provide a functional explanation for individual differences in innate response patterns. Our identified reQTLs are also associated with cancer, autoimmunity, inflammatory and infectious diseases as shown by external genome-wide association studies. Thus, reQTLs help to explain interindividual variation in immune response to infection and provide candidate genes for variants associated with a range of diseases.

Dissecting the genetic heterogeneity of gastric cancer.

BACKGROUND: Gastric cancer (GC) is clinically heterogenous according to location (cardia/non-cardia) and histopathology (diffuse/intestinal). We aimed to characterize the genetic risk architecture of GC according to its subtypes. Another aim was to examine whether cardia GC and oesophageal adenocarcinoma (OAC) and its precursor lesion Barrett's oesophagus (BO), which are all located at the gastro-oesophageal junction (GOJ), share polygenic risk architecture. METHODS: We did a meta-analysis of ten European genome-wide association studies (GWAS) of GC and its subtypes. All patients had a histopathologically confirmed diagnosis of gastric adenocarcinoma. For the identification of risk genes among GWAS loci we did a transcriptome-wide association study (TWAS) and expression quantitative trait locus (eQTL) study from gastric corpus and antrum mucosa. To test whether cardia GC and OAC/BO share genetic aetiology we also used a European GWAS sample with OAC/BO. FINDINGS: Our GWAS consisting of 5816 patients and 10,999 controls highlights the genetic heterogeneity of GC according to its subtypes. We newly identified two and replicated five GC risk loci, all of them with subtype-specific association. The gastric transcriptome data consisting of 361 corpus and 342 antrum mucosa samples revealed that an upregulated expression of MUC1, ANKRD50, PTGER4, and PSCA are plausible GC-pathomechanisms at four GWAS loci. At another risk locus, we found that the blood-group 0 exerts protective effects for non-cardia and diffuse GC, while blood-group A increases risk for both GC subtypes. Furthermore, our GWAS on cardia GC and OAC/BO (10,279 patients, 16,527 controls) showed that both cancer entities share genetic aetiology at the polygenic level and identified two new risk loci on the single-marker level. INTERPRETATION: Our findings show that the pathophysiology of GC is genetically heterogenous according to location and histopathology. Moreover, our findings point to common molecular mechanisms underlying cardia GC and OAC/BO. FUNDING: German Research Foundation (DFG).

Characteristics of familial pancreatic cancer families with additional colorectal carcinoma.

Familial pancreatic cancer (FPC) is a rare hereditary tumor entity with broad phenotypic heterogeneity, including colorectal carcinoma (CRC) in some families. The underlying factors for this co-occurrence are still not well evaluated. FPC families in the National Case Collection of Familial Pancreatic Cancer with an additional occurrence of CRC were analyzed regarding the phenotype, genotype and recommendation for a clinical screening program. The total cohort of 272 FPC families included 30 (11%) families with at least one CRC case. The proportion of affected family members with PDAC was 16.1% (73/451) compared to 9.3% of family members with CRC (42/451, p < 0.01). Females were affected with PDAC in 49% (36/73) and CRC in 38% (16/42). The median age of PDAC was 63 compared to 66 years in CRC, whereas 8 (26.6%) of families had an early onset of PDAC and 2 (6.7%) of CRC. Seventeen families had 2 or more affected generations with PDAC and 6 families with CRC. Eleven (9.6%) of affected patients had both PDAC and CRC. Potentially causative germline mutations (2 ATM, 1 CDKN2a, 1 MLH1, 1 PALB2) were detected in 5 of 18 (27.7%) analyzed cases. These findings provide a step forward to include the phenotypic and genotypic characteristics of FPC-CRC families for the genetic counseling and management of these families. Nevertheless, results need to be verified in a larger patient cohort beforehand.

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.

First genome-wide association study of esophageal atresia identifies three genetic risk loci at CTNNA3, FOXF1/FOXC2/FOXL1, and HNF1B.

Esophageal atresia with or without tracheoesophageal fistula (EA/TEF) is the most common congenital malformation of the upper digestive tract. This study represents the first genome-wide association study (GWAS) to identify risk loci for EA/TEF. We used a European case-control sample comprising 764 EA/TEF patients and 5,778 controls and observed genome-wide significant associations at three loci. On chromosome 10q21 within the gene CTNNA3 (p = 2.11 × 10-8; odds ratio [OR] = 3.94; 95% confidence interval [CI], 3.10-5.00), on chromosome 16q24 next to the FOX gene cluster (p = 2.25 × 10-10; OR = 1.47; 95% CI, 1.38-1.55) and on chromosome 17q12 next to the gene HNF1B (p = 3.35 × 10-16; OR = 1.75; 95% CI, 1.64-1.87). We next carried out an esophageal/tracheal transcriptome profiling in rat embryos at four selected embryonic time points. Based on these data and on already published data, the implicated genes at all three GWAS loci are promising candidates for EA/TEF development. We also analyzed the genetic EA/TEF architecture beyond the single marker level, which revealed an estimated single-nucleotide polymorphism (SNP)-based heritability of around 37% ± 14% standard deviation. In addition, we examined the polygenicity of EA/TEF and found that EA/TEF is less polygenic than other complex genetic diseases. In conclusion, the results of our study contribute to a better understanding on the underlying genetic architecture of ET/TEF with the identification of three risk loci and candidate genes.

Clinical Relevance of Gastroesophageal Cancer Associated SNPs for Oncologic Outcome After Curative Surgery.

BACKGROUND: Gastric and esophageal cancers are malignant diseases with rising importance in Western countries. To improve oncologic outcome after surgery, it is essential to understand the relevance of germline mutations. The aim of the study was to identify and distinguish clinically relevant single-nucleotide polymorphisms (SNPs). PATIENTS AND METHODS: In total, 190 patients with curative oncological resections of gastric and distal esophageal adenocarcinomas at Heidelberg University Hospital were eligible for this study. Outcome differences were determined for each SNP by analysis of clinical variables, survival, and mRNA expression levels. RESULTS: Significant survival differences were found on univariate analysis for usual prognostic variables (such as pTNM) and for six SNPs. On multivariate survival analysis, the SNPs rs12268840 (intron variant of MGMT, p = 0.045) and rs9972882 (intron variant of STARD3 and eQTL of PGAP3, p = 0.030) were independent and significant survival predictors along with R status and pT/pN category. Group TT of rs12268840 had the highest rate of second primary carcinoma (30.4%, p = 0.0003), lowest expression of MGMT based on cis-eQTL analysis in normal gastroesophageal tissue (p = 1.99 × 10-17), and worst oncologic outcome. Group AA of rs9972882 had the highest rate of distant metastases pM1 (42.9%, p = 0.0117), highest expression of PGAP3 (p = 1.29 × 10-15), and worst oncologic outcome. CONCLUSIONS: Two intron variant SNPs of MGMT and STARD3 were identified that were significant survival predictors and may influence tumor biology. The data indicate that DNA methylation (MGMT) and malfunction of GPI anchoring (PGAP3) are distinct mechanisms that are relevant for tumor progression and relapse.

Genome-wide association study of panic disorder reveals genetic overlap with neuroticism and depression.

Panic disorder (PD) has a lifetime prevalence of 2-4% and heritability estimates of 40%. The contributory genetic variants remain largely unknown, with few and inconsistent loci having been reported. The present report describes the largest genome-wide association study (GWAS) of PD to date comprising genome-wide genotype data of 2248 clinically well-characterized PD patients and 7992 ethnically matched controls. The samples originated from four European countries (Denmark, Estonia, Germany, and Sweden). Standard GWAS quality control procedures were conducted on each individual dataset, and imputation was performed using the 1000 Genomes Project reference panel. A meta-analysis was then performed using the Ricopili pipeline. No genome-wide significant locus was identified. Leave-one-out analyses generated highly significant polygenic risk scores (PRS) (explained variance of up to 2.6%). Linkage disequilibrium (LD) score regression analysis of the GWAS data showed that the estimated heritability for PD was 28.0-34.2%. After correction for multiple testing, a significant genetic correlation was found between PD and major depressive disorder, depressive symptoms, and neuroticism. A total of 255 single-nucleotide polymorphisms (SNPs) with p < 1 × 10-4 were followed up in an independent sample of 2408 PD patients and 228,470 controls from Denmark, Iceland and the Netherlands. In the combined analysis, SNP rs144783209 showed the strongest association with PD (pcomb = 3.10 × 10-7). Sign tests revealed a significant enrichment of SNPs with a discovery p-value of <0.0001 in the combined follow up cohort (p = 0.048). The present integrative analysis represents a major step towards the elucidation of the genetic susceptibility to PD.

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.

Toll-like Receptor Signaling Rewires Macrophage Metabolism and Promotes Histone Acetylation via ATP-Citrate Lyase.

Toll-like receptor (TLR) activation induces inflammatory responses in macrophages by activating temporally defined transcriptional cascades. Whether concurrent changes in the cellular metabolism that occur upon TLR activation influence the quality of the transcriptional responses remains unknown. Here, we investigated how macrophages adopt their metabolism early after activation to regulate TLR-inducible gene induction. Shortly after TLR4 activation, macrophages increased glycolysis and tricarboxylic acid (TCA) cycle volume. Metabolic tracing studies revealed that TLR signaling redirected metabolic fluxes to generate acetyl-Coenzyme A (CoA) from glucose resulting in augmented histone acetylation. Signaling through the adaptor proteins MyD88 and TRIF resulted in activation of ATP-citrate lyase, which in turn facilitated the induction of distinct LPS-inducible gene sets. We postulate that metabolic licensing of histone acetylation provides another layer of control that serves to fine-tune transcriptional responses downstream of TLR activation. Our work highlights the potential of targeting the metabolic-epigenetic axis in inflammatory settings.

First genotype-phenotype study reveals HLA-DQß1 insertion heterogeneity in high-resolution manometry achalasia subtypes.

Background: Achalasia is a primary oesophageal motility disorder. Although aetiology remains mainly unknown, a genetic risk variant, rs28688207 in HLA-DQB1, showed strong achalasia association suggesting involvement of immune-mediated processes in the pathogenesis. High-resolution manometry recognises three types of achalasia. The aim of our study was to perform the first genotype-phenotype analysis investigating the frequency of rs28688207 across the high-resolution manometry subtypes. Methods: This was a cross-sectional retrospective study. Achalasia patients from tertiary centres in the Czech Republic (n = 163), Germany (n = 114), Greece (n = 70) and controls were enrolled. All subjects were genotyped for the rs28688207 insertion. The Kruskal-Wallis test was used for the genotype-phenotype analysis. Results: A total of 347 achalasia patients (type I - 89, II - 210, III - 48) were included. The overall frequency of the rs28688207 was 10.3%. The distribution of the insertion was significantly different across the high-resolution manometry subtypes (p = 0.038), being most prevalent in type I (14.6%), followed by type II (9.5%) and III (6.3%). Conclusion: The frequency of the HLA-DQB1 insertion differs among high-resolution manometry achalasia subtypes. The insertion is most prevalent in type I, suggesting that immune-mediated mechanisms triggered by the insertion may play a more prominent role in the pathogenesis of this subtype.

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.

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.

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.

Comprehensive epidemiological and genotype-phenotype analyses in a large European sample with idiopathic achalasia.

BACKGROUND AND AIM: Although an eight-residue insertion in HLA-DQß1 has been recently identified as a genetic risk factor for idiopathic achalasia, other risk factors are still unknown. In the present study, we carried out an epidemiological survey and a genotype-phenotype (G×P) analysis to gain further insights into the etiology of achalasia. METHODS: We obtained medical data from 696 achalasia patients and 410 controls, as well as their first-degree relatives (2543 of patients and 1497 of controls). For the G×P analysis, we stratified the patients into HLA-DQß1 insertion carriers and noncarriers. RESULTS: Our data show that patients are more often affected by viral infections before achalasia onset (P<0.0001, most significantly for varicella zoster virus infections). In addition, allergic (P=0.0005) and autoimmune disorders (P=0.0007, most significantly for psoriasis and Sjögren's syndrome) represent comorbid disease conditions. First-degree relatives of patients also show higher prevalence rates of allergic disorders (P=0.0007) and psoriasis (P=0.016) compared with control relatives. Moreover, the G×P analysis reveals that achalasia is triggered by pregnancies in female HLA-DQß1 insertion carriers (P=0.031). CONCLUSION: Our data point to a role of viral infections in the development of achalasia. In addition, they provide evidence for a relationship between achalasia and allergic, as well as autoimmune, disorders. Furthermore, pregnancy seems to be a disease-triggering factor in female HLA-DQß1 insertion carriers, which points to hormonal and/or immunosuppressive factors influencing disease development.

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.

The HLA-DQß1 insertion is a strong achalasia risk factor and displays a geospatial north-south gradient among Europeans.

Idiopathic achalasia is a severe motility disorder of the esophagus and is characterized by a failure of the lower esophageal sphincter to relax due to a loss of neurons in the myenteric plexus. Most recently, we identified an eight-amino-acid insertion in the cytoplasmic tail of HLA-DQß1 as strong achalasia risk factor in a sample set from Central Europe, Italy and Spain. Here, we tested whether the HLA-DQß1 insertion also confers achalasia risk in the Polish and Swedish population. We could replicate the initial findings and the insertion shows strong achalasia association in both samples (Poland P=1.84 × 10(-04), Sweden P=7.44 × 10(-05)). Combining all five European data sets - Central Europe, Italy, Spain, Poland and Sweden - the insertion is achalasia associated with Pcombined=1.67 × 10(-35). In addition, we observe that the frequency of the insertion shows a geospatial north-south gradient. The insertion is less common in northern (around 6-7% in patients and 2% in controls from Sweden and Poland) compared with southern Europeans (~16% in patients and 8% in controls from Italy) and shows a stronger attributable risk in the southern European population. Our study provides evidence that the prevalence of achalasia may differ between populations.