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 MBD4 deficiency causes a multi-tumor predisposition syndrome.

We report an autosomal recessive, multi-organ tumor predisposition syndrome, caused by bi-allelic loss-of-function germline variants in the base excision repair (BER) gene MBD4. We identified five individuals with bi-allelic MBD4 variants within four families and these individuals had a personal and/or family history of adenomatous colorectal polyposis, acute myeloid leukemia, and uveal melanoma. MBD4 encodes a glycosylase involved in repair of G:T mismatches resulting from deamination of 5'-methylcytosine. The colorectal adenomas from MBD4-deficient individuals showed a mutator phenotype attributable to mutational signature SBS1, consistent with the function of MBD4. MBD4-deficient polyps harbored somatic mutations in similar driver genes to sporadic colorectal tumors, although AMER1 mutations were more common and KRAS mutations less frequent. Our findings expand the role of BER deficiencies in tumor predisposition. Inclusion of MBD4 in genetic testing for polyposis and multi-tumor phenotypes is warranted to improve disease management.

The clinical features of polymerase proof-reading associated polyposis (PPAP) and recommendations for patient management.

Pathogenic germline exonuclease domain (ED) variants of POLE and POLD1 cause the Mendelian dominant condition polymerase proof-reading associated polyposis (PPAP). We aimed to describe the clinical features of all PPAP patients with probably pathogenic variants. We identified patients with a variants mapping to the EDs of POLE or POLD1 from cancer genetics clinics, a colorectal cancer (CRC) clinical trial, and systematic review of the literature. We used multiple evidence sources to separate ED variants into those with strong evidence of pathogenicity and those of uncertain importance. We performed quantitative analysis of the risk of CRC, colorectal adenomas, endometrial cancer or any cancer in the former group. 132 individuals carried a probably pathogenic ED variant (105 POLE, 27 POLD1). The earliest malignancy was colorectal cancer at 14. The most common tumour types were colorectal, followed by endometrial in POLD1 heterozygotes and duodenal in POLE heterozygotes. POLD1-mutant cases were at a significantly higher risk of endometrial cancer than POLE heterozygotes. Five individuals with a POLE pathogenic variant, but none with a POLD1 pathogenic variant, developed ovarian cancer. Nine patients with POLE pathogenic variants and one with a POLD1 pathogenic variant developed brain tumours. Our data provide important evidence for PPAP management. Colonoscopic surveillance is recommended from age 14 and upper-gastrointestinal surveillance from age 25. The management of other tumour risks remains uncertain, but surveillance should be considered. In the absence of strong genotype-phenotype associations, these recommendations should apply to all PPAP patients.

An Evaluation of the Diagnostic Accuracy of a Panel of Variants in DPYD and a Single Variant in ENOSF1 for Predicting Common Capecitabine Related Toxicities.

Efficacy of 5-Fluorouracil (5-FU)-based chemotherapy is limited by significant toxicity. Tests based upon variants in the Clinical Pharmacogenetics Implementation Consortium (CPIC) guidelines with high level evidence of a link to dihydropyrimidine dehydrogenase (DPD) phenotype and 5-FU toxicity are available to identify patients at high risk of severe adverse events (AEs). We previously reported associations between rs1213215, rs2612091, and NM_000110.3:c.1906-14763G>A (rs12022243) and capecitabine induced toxicity in clinical trial QUASAR 2. We also identified patients with DPD deficiency alleles NM_000110.3: c.1905+1G>A, NM_000110.3: c.2846C>T, NM_000110.3:c.1679T>G and NM_000110.3:c.1651G>A. We have now assessed the frequency of thirteen additional DPYD deficiency variants in 888 patients from the QUASAR 2 clinical trial. We also compared the area under the curve (AUC)-a measure of diagnostic accuracy-of the high-level evidence variants from the CPIC guidelines plus and minus additional DPYD deficiency variants and or common variants associated with 5-FU toxicity. Including additional DPYD deficiency variants retained good diagnostic accuracy for serious adverse events (AEs) and improved sensitivity for predicting grade 4 haematological toxicities (sensitivity 0.75, specificity 0.94) but the improvement in AUC for this toxicity was not significant. Larger datasets will be required to determine the benefit of including additional DPYD deficiency variants not observed here. Genotyping two common alleles statistically significantly improves AUC for prediction of risk of HFS and may be clinically useful (AUC difference 0.177, sensitivity 0.84, specificity 0.31).

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.

The polymorphic variant rs1800734 influences methylation acquisition and allele-specific TFAP4 binding in the MLH1 promoter leading to differential mRNA expression.

Expression of the mismatch repair gene MutL homolog 1 (MLH1) is silenced in a clinically important subgroup of sporadic colorectal cancers. These cancers exhibit hypermutability with microsatellite instability (MSI) and differ from microsatellite-stable (MSS) colorectal cancers in both prognosis and response to therapies. Loss of MLH1 is usually due to epigenetic silencing with associated promoter methylation; coding somatic mutations rarely occur. Here we use the presence of a colorectal cancer (CRC) risk variant (rs1800734) within the MLH1 promoter to investigate the poorly understood mechanisms of MLH1 promoter methylation and loss of expression. We confirm the association of rs1800734 with MSI+ but not MSS cancer risk in our own data and by meta-analysis. Using sensitive allele-specific detection methods, we demonstrate that MLH1 is the target gene for rs1800734 mediated cancer risk. In normal colon tissue, small allele-specific differences exist only in MLH1 promoter methylation, but not gene expression. In contrast, allele-specific differences in both MLH1 methylation and expression are present in MSI+ cancers. We show that MLH1 transcriptional repression is dependent on DNA methylation and can be reversed by a methylation inhibitor. The rs1800734 allele influences the rate of methylation loss and amount of re-expression. The transcription factor TFAP4 binds to the rs1800734 region but with much weaker binding to the risk than the protective allele. TFAP4 binding is absent on both alleles when promoter methylation is present. Thus we propose that TFAP4 binding shields the protective rs1800734 allele of the MLH1 promoter from BRAF induced DNA methylation more effectively than the risk allele.

The evolutionary landscape of colorectal tumorigenesis.

The evolutionary events that cause colorectal adenomas (benign) to progress to carcinomas (malignant) remain largely undetermined. Using multi-region genome and exome sequencing of 24 benign and malignant colorectal tumours, we investigate the evolutionary fitness landscape occupied by these neoplasms. Unlike carcinomas, advanced adenomas frequently harbour sub-clonal driver mutations-considered to be functionally important in the carcinogenic process-that have not swept to fixation, and have relatively high genetic heterogeneity. Carcinomas are distinguished from adenomas by widespread aneusomies that are usually clonal and often accrue in a 'punctuated' fashion. We conclude that adenomas evolve across an undulating fitness landscape, whereas carcinomas occupy a sharper fitness peak, probably owing to stabilizing selection.

Somatic POLE exonuclease domain mutations are early events in sporadic endometrial and colorectal carcinogenesis, determining driver mutational landscape, clonal neoantigen burden and immune response.

Genomic instability, which is a hallmark of cancer, is generally thought to occur in the middle to late stages of tumourigenesis, following the acquisition of permissive molecular aberrations such as TP53 mutation or whole genome doubling. Tumours with somatic POLE exonuclease domain mutations are notable for their extreme genomic instability (their mutation burden is among the highest in human cancer), distinct mutational signature, lymphocytic infiltrate, and excellent prognosis. To what extent these characteristics are determined by the timing of POLE mutations in oncogenesis is unknown. Here, we have shown that pathogenic POLE mutations are detectable in non-malignant precursors of endometrial and colorectal cancer. Using genome and exome sequencing, we found that multiple driver mutations in POLE-mutant cancers show the characteristic POLE mutational signature, including those in genes conventionally regarded as initiators of tumourigenesis. In POLE-mutant cancers, the proportion of monoclonal predicted neoantigens was similar to that in other cancers, but the absolute number was much greater. We also found that the prominent CD8+ T-cell infiltrate present in POLE-mutant cancers was evident in their precursor lesions. Collectively, these data indicate that somatic POLE mutations are early, quite possibly initiating, events in the endometrial and colorectal cancers in which they occur. The resulting early onset of genomic instability may account for the striking immune response and excellent prognosis of these tumours, as well as their early presentation. © 2018 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of Pathological Society of Great Britain and Ireland.

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

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

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.

Polymorphisms near TBX5 and GDF7 are associated with increased risk for Barrett's esophagus.

BACKGROUND & AIMS: Barrett's esophagus (BE) increases the risk of esophageal adenocarcinoma (EAC). We found the risk to be BE has been associated with single nucleotide polymorphisms (SNPs) on chromosome 6p21 (within the HLA region) and on 16q23, where the closest protein-coding gene is FOXF1. Subsequently, the Barrett's and Esophageal Adenocarcinoma Consortium (BEACON) identified risk loci for BE and esophageal adenocarcinoma near CRTC1 and BARX1, and within 100 kb of FOXP1. We aimed to identify further SNPs that increased BE risk and to validate previously reported associations. METHODS: We performed a genome-wide association study (GWAS) to identify variants associated with BE and further analyzed promising variants identified by BEACON by genotyping 10,158 patients with BE and 21,062 controls. RESULTS: We identified 2 SNPs not previously associated with BE: rs3072 (2p24.1; odds ratio [OR] = 1.14; 95% CI: 1.09-1.18; P = 1.8 × 10(-11)) and rs2701108 (12q24.21; OR = 0.90; 95% CI: 0.86-0.93; P = 7.5 × 10(-9)). The closest protein-coding genes were respectively GDF7 (rs3072), which encodes a ligand in the bone morphogenetic protein pathway, and TBX5 (rs2701108), which encodes a transcription factor that regulates esophageal and cardiac development. Our data also supported in BE cases 3 risk SNPs identified by BEACON (rs2687201, rs11789015, and rs10423674). Meta-analysis of all data identified another SNP associated with BE and esophageal adenocarcinoma: rs3784262, within ALDH1A2 (OR = 0.90; 95% CI: 0.87-0.93; P = 3.72 × 10(-9)). CONCLUSIONS: We identified 2 loci associated with risk of BE and provided data to support a further locus. The genes we found to be associated with risk for BE encode transcription factors involved in thoracic, diaphragmatic, and esophageal development or proteins involved in the inflammatory response.