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.

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.

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.

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.

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.

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.

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.

MCM3AP and POMP Mutations Cause a DNA-Repair and DNA-Damage-Signaling Defect in an Immunodeficient Child.

Immunodeficiency patients with DNA repair defects exhibit radiosensitivity and proneness to leukemia/lymphoma formation. Though progress has been made in identifying the underlying mutations, in most patients the genetic basis is unknown. Two de novo mutated candidate genes, MCM3AP encoding germinal center-associated nuclear protein (GANP) and POMP encoding proteasome maturation protein (POMP), were identified by whole-exome sequencing (WES) and confirmed by Sanger sequencing in a child with complex phenotype displaying immunodeficiency, genomic instability, skin changes, and myelodysplasia. GANP was previously described to promote B-cell maturation by nuclear targeting of activation-induced cytidine deaminase (AID) and to control AID-dependent hyperrecombination. POMP is required for 20S proteasome assembly and, thus, for efficient NF-κB signaling. Patient-derived cells were characterized by impaired homologous recombination, moderate radio- and cross-linker sensitivity associated with accumulation of damage, impaired DNA damage-induced NF-κB signaling, and reduced nuclear AID levels. Complementation by wild-type (WT)-GANP normalized DNA repair and WT-POMP rescued defective NF-κB signaling. In conclusion, we identified for the first time mutations in MCM3AP and POMP in an immunodeficiency patient. These mutations lead to cooperative effects on DNA recombination and damage signaling. Digenic/polygenic mutations may constitute a novel genetic basis in immunodeficiency patients with DNA repair defects.

Supportive evidence for FOXP1, BARX1, and FOXF1 as genetic risk loci for the development of esophageal adenocarcinoma.

The Barrett's and Esophageal Adenocarcinoma Consortium (BEACON) recently performed a genome-wide association study (GWAS) on esophageal adenocarcinoma (EAC) and Barrett's esophagus. They identified genome-wide significant association for variants at three genes, namely CRTC1, FOXP1, and BARX1. Furthermore, they replicated an association at the FOXF1 gene that has been previously found in a GWAS on Barrett's esophagus. We aimed at further replicating the association at these and other loci that showed suggestive association with P < 10(-4) in the BEACON sample. In total, we tested 88 SNPs in an independent sample consisting of 1065 EAC cases and 1019 controls of German descent. We could replicate the association at FOXP1, BARX1, and FOXF1 with nominal significance and thereby confirm that genetic variants at these genes confer EAC risk. In addition, we found association of variants near the genes XRCC2 and GATA6 that were strongly (P < 10(-5) ) although not genome-wide significantly associated with the BEACON GWAS. Therefore, both variants and corresponding genes represent promising candidates for future EAC association studies on independent samples.

Characterizing the genetic basis of innate immune response in TLR4-activated human monocytes.

Toll-like receptors (TLRs) play a key role in innate immunity. Apart from their function in host defense, dysregulation in TLR signalling can confer risk to autoimmune diseases, septic shock or cancer. Here we report genetic variants and transcripts that are active only during TLR signalling and contribute to interindividual differences in immune response. Comparing unstimulated versus TLR4-stimulated monocytes reveals 1,471 expression quantitative trait loci (eQTLs) that are unique to TLR4 stimulation. Among these we find functional SNPs for the expression of NEU4, CCL14, CBX3 and IRF5 on TLR4 activation. Furthermore, we show that SNPs conferring risk to primary biliary cirrhosis (PBC), inflammatory bowel disease (IBD) and celiac disease are immune response eQTLs for PDGFB and IL18R1. Thus, PDGFB and IL18R1 represent plausible candidates for studying the pathophysiology of these disorders in the context of TLR4 activation. In summary, this study presents novel insights into the genetic basis of the innate immune response and exemplifies the value of eQTL studies in the context of exogenous cell stimulation.