Targeted panel sequencing of pharmacogenes and oncodrivers in colorectal cancer patients reveals genes with prognostic significance.

BACKGROUND: Colorectal cancer is still the second leading cause of cancer-related deaths and thus biomarkers allowing prediction of the resistance of patients to therapy and estimating their prognosis are needed. We designed a panel of 558 genes with pharmacogenomics records related to 5-fluorouracil resistance, genes important for sensitivity to other frequently used drugs, major oncodrivers, and actionable genes. We performed a target enrichment sequencing of DNA from tumors and matched blood samples of patients, and compared the results with patient prognosis stratified by systemic adjuvant chemotherapy. RESULTS: The median number of detected variants per tumor sample was 18.5 with 4 classified as having a high predicted functional effect and 14.5 moderate effect. APC, TP53, and KRAS were the most frequent mutated genes (64%, 59%, and 42% of mutated samples, respectively) followed by FAT4 (23%), FBXW7, and PIK3CA (16% for both). Patients with advanced stage III had more frequently APC, TP53, or KRAS mutations than those in stages I or II. KRAS mutation counts followed an increasing trend with grade (G1 < G2 < G3). The response to adjuvant therapy was worse in carriers of frameshift mutations in APC or 12D variant in KRAS, but none of these oncodrivers had prognostic value. Carriage of somatic mutations in any of the genes ABCA13, ANK2, COL7A1, NAV3, or UNC80 had prognostic relevance for worse overall survival (OS) of all patients. In contrast, mutations in FLG, GLI3, or UNC80 were prognostic in the same direction for patients untreated, and mutations in COL6A3, LRP1B, NAV3, RYR1, RYR3, TCHH, or TENM4 for patients treated with adjuvant therapy. The first association was externally validated. From all germline variants with high or moderate predicted functional effects (median 326 per patient), > 5% frequency and positive Manhattan plot based on 3-year RFS, rs72753407 in NFACS, rs34621071 in ERBB4, and rs2444274 in RIF1 were significantly associated with RFS, OS or both. CONCLUSIONS: The present study identified several putative somatic and germline genetic events with prognostic potential for colorectal cancer that should undergo functional characterization.

Artificial Intelligence-Driven Prediction Revealed CFTR Associated with Therapy Outcome of Breast Cancer: A Feasibility Study.

INTRODUCTION: In silico tools capable of predicting the functional consequences of genomic differences between individuals, many of which are AI-driven, have been the most effective over the past two decades for non-synonymous single nucleotide variants (nsSNVs). When appropriately selected for the purpose of the study, a high predictive performance can be expected. In this feasibility study, we investigate the distribution of nsSNVs with an allele frequency below 5%. To classify the putative functional consequence, a tier-based filtration led by AI-driven predictors and scoring system was implemented to the overall decision-making process, resulting in a list of prioritised genes. METHODS: The study has been conducted on breast cancer patients of homogeneous ethnicity. Germline rare variants have been sequenced in genes that influence pharmacokinetic parameters of anticancer drugs or molecular signalling pathways in cancer. After AI-driven functional pathogenicity classification and data mining in pharmacogenomic (PGx) databases, variants were collapsed to the gene level and ranked according to their putative deleterious role. RESULTS: In breast cancer patients, seven of the twelve genes prioritised based on the predictions were found to be associated with response to oncotherapy, histological grade, and tumour subtype. Most importantly, we showed that the group of patients with at least one rare nsSNVs in cystic fibrosis transmembrane conductance regulator (CFTR) had significantly reduced disease-free (log rank, p = 0.002) and overall survival (log rank, p = 0.006). CONCLUSION: AI-driven in silico analysis with PGx data mining provided an effective approach navigating for functional consequences across germline genetic background, which can be easily integrated into the overall decision-making process for future studies. The study revealed a statistically significant association with numerous clinicopathological parameters, including treatment response. Our study indicates that CFTR may be involved in the processes influencing the effectiveness of oncotherapy or in the malignant progression of the disease itself.

Potential association between PSCA rs2976395 functional variant and pancreatic cancer risk.

Correlated regions of systemic interindividual variation (CoRSIV) represent a small proportion of the human genome showing DNA methylation patterns that are the same in all human tissues, are different among individuals, and are partially regulated by genetic variants in cis. In this study we aimed at investigating single-nucleotide polymorphisms (SNPs) within CoRSIVs and their involvement with pancreatic ductal adenocarcinoma (PDAC) risk. We analyzed 29,099 CoRSIV-SNPs and 133,615 CoRSIV-mQTLs in 14,394 cases and 247,022 controls of European and Asian descent. We observed that the A allele of the rs2976395 SNP was associated with increased PDAC risk in Europeans (p = 2.81 × 10-5). This SNP lies in the prostate stem cell antigen gene and is in perfect linkage disequilibrium with a variant (rs2294008) that has been reported to be associated with risk of many other cancer types. The A allele is associated with the DNA methylation level of the gene according to the PanCan-meQTL database and with overexpression according to QTLbase. The expression of the gene has been observed to be deregulated in many tumors of the gastrointestinal tract including pancreatic cancer; however, functional studies are needed to elucidate the function relevance of the association.

Genes divided according to the relative position of the longest intron show increased representation in different KEGG pathways.

Despite the fact that introns mean an energy and time burden for eukaryotic cells, they play an irreplaceable role in the diversification and regulation of protein production. As a common feature of eukaryotic genomes, it has been reported that in protein-coding genes, the longest intron is usually one of the first introns. The goal of our work was to find a possible difference in the biological function of genes that fulfill this common feature compared to genes that do not. Data on the lengths of all introns in genes were extracted from the genomes of six vertebrates (human, mouse, koala, chicken, zebrafish and fugu) and two other model organisms (nematode worm and arabidopsis). We showed that more than 40% of protein-coding genes have the relative position of the longest intron located in the second or third tertile of all introns. Genes divided according to the relative position of the longest intron were found to be significantly increased in different KEGG pathways. Genes with the longest intron in the first tertile predominate in a range of pathways for amino acid and lipid metabolism, various signaling, cell junctions or ABC transporters. Genes with the longest intron in the second or third tertile show increased representation in pathways associated with the formation and function of the spliceosome and ribosomes. In the two groups of genes defined in this way, we further demonstrated the difference in the length of the longest introns and the distribution of their absolute positions. We also pointed out other characteristics, namely the positive correlation between the length of the longest intron and the sum of the lengths of all other introns in the gene and the preservation of the exact same absolute and relative position of the longest intron between orthologous genes.

A Long-Term Study on the Bactericidal Effect of ZrN-Cu Nanostructured Coatings Deposited by an Industrial Physical Vapor Deposition System.

ZrN-Cu coatings containing two different amounts of Cu (~11 at.% and ~25 at.%) were deposited using an industrial physical vapor deposition (PVD) system. The as-deposited coatings exhibited 100% bactericidal efficiency against Escherichia coli CCM 3988 for an exposure time of 40 min. Subsequently, the samples were attached onto our faculty's door handles for six months to study the coatings' long-term effectiveness and durability under actual operational conditions. The samples were periodically evaluated and it was observed that the coatings with 25 at.% Cu performed better than the ones with 11 at.% Cu. For example, following 15 days of being touched, the bactericidal effectiveness of the sample containing 25 at.% Cu dropped to 65% while it fell to 42% for the sample containing 11 at.%. After 6 months, however, both samples showed bactericidal efficiency of ~16-20%. The bactericidal efficiency of the samples touched for 6 months was successfully restored by polishing them. Furthermore, a group of samples was kept untouched and was also evaluated. The untouched samples with Cu content of ~25 at.% did not show any drop in their bactericidal properties after 6 months. ZrN-Cu coatings were concluded to be promising materials for self-sanitizing application on high-touch surfaces.

Targeted DNA sequencing of high-grade serous ovarian carcinoma reveals association of TP53 mutations with platinum resistance when combined with gene expression.

High-grade serous ovarian carcinoma (HGSC) is the most common subtype of ovarian cancer and is among the most fatal gynecological malignancies worldwide, due to late diagnosis at advanced stages and frequent therapy resistance. In 47 HGSC patients, we assessed somatic and germline genetic variability of a custom panel of 144 known or suspected HGSC-related genes by high-coverage targeted DNA sequencing to identify the genetic determinants associated with resistance to platinum-based therapy. In the germline, the most mutated genes were DNAH14 (17%), RAD51B (17%), CFTR (13%), BRCA1 (11%), and RAD51 (11%). Somatically, the most mutated gene was TP53 (98%), followed by CSMD1/2/3 (19/19/36%), and CFTR (23%). Results were compared with those from whole exome sequencing of a similar set of 35 HGSC patients. Somatic variants in TP53 were also validated using GENIE data of 1287 HGSC samples. Our approach showed increased prevalence of high impact somatic and germline mutations, especially those affecting splice sites of TP53, compared to validation datasets. Furthermore, nonsense TP53 somatic mutations were negatively associated with patient survival. Elevated TP53 transcript levels were associated with platinum resistance and presence of TP53 missense mutations, while decreased TP53 levels were found in tumors carrying mutations with predicted high impact, which was confirmed in The Cancer Genome Atlas data (n = 260). Targeted DNA sequencing of TP53 combined with transcript quantification may contribute to the concept of precision oncology of HGSC. Future studies should explore targeting the p53 pathway based on specific mutation types and co-analyze the expression and mutational profiles of other key cancer genes.

Third-generation taxanes SB-T-121605 and SB-T-121606 are effective in pancreatic ductal adenocarcinoma.

Pancreatic cancer is a severe malignancy with increasing incidence and high mortality due to late diagnosis and low sensitivity to treatments. Search for the most appropriate drugs and therapeutic regimens is the most promising way to improve the treatment outcomes of the patients. This study aimed to compare (1) in vitro efficacy and (2) in vivo antitumor effects of conventional paclitaxel and the newly synthesized second (SB-T-1216) and third (SB-T-121605 and SB-T-121606) generation taxanes in KRAS wild type BxPC-3 and more aggressive KRAS G12V mutated Paca-44 pancreatic cancer cell line models. In vitro, paclitaxel efficacy was 27.6 ± 1.7 nM, while SB-Ts showed 1.7-7.4 times higher efficacy. Incorporation of SB-T-121605 and SB-T-121606 into in vivo therapeutic regimens containing paclitaxel was effective in suppressing tumor growth in Paca-44 tumor-bearing mice at small doses (≤3 mg/kg). SB-T-121605 and SB-T-121606 in combination with paclitaxel are promising candidates for the next phase of preclinical testing.

Polymorphisms in transcription factor binding sites and enhancer regions and pancreatic ductal adenocarcinoma risk.

Genome-wide association studies (GWAS) are a powerful tool for detecting variants associated with complex traits and can help risk stratification and prevention strategies against pancreatic ductal adenocarcinoma (PDAC). However, the strict significance threshold commonly used makes it likely that many true risk loci are missed. Functional annotation of GWAS polymorphisms is a proven strategy to identify additional risk loci. We aimed to investigate single-nucleotide polymorphisms (SNP) in regulatory regions [transcription factor binding sites (TFBSs) and enhancers] that could change the expression profile of multiple genes they act upon and thereby modify PDAC risk. We analyzed a total of 12,636 PDAC cases and 43,443 controls from PanScan/PanC4 and the East Asian GWAS (discovery populations), and the PANDoRA consortium (replication population). We identified four associations that reached study-wide statistical significance in the overall meta-analysis: rs2472632(A) (enhancer variant, OR 1.10, 95%CI 1.06,1.13, p = 5.5 × 10-8), rs17358295(G) (enhancer variant, OR 1.16, 95%CI 1.10,1.22, p = 6.1 × 10-7), rs2232079(T) (TFBS variant, OR 0.88, 95%CI 0.83,0.93, p = 6.4 × 10-6) and rs10025845(A) (TFBS variant, OR 1.88, 95%CI 1.50,1.12, p = 1.32 × 10-5). The SNP with the most significant association, rs2472632, is located in an enhancer predicted to target the coiled-coil domain containing 34 oncogene. Our results provide new insights into genetic risk factors for PDAC by a focused analysis of polymorphisms in regulatory regions and demonstrating the usefulness of functional prioritization to identify loci associated with PDAC risk.

Evaluation of human papillomavirus DNA in colorectal cancer and adjacent mucosal tissue samples.

BACKGROUND: Although the role of viral agents, such as human papillomavirus (e.g. HPV16, HPV18) in colorectal cancer (CRC) has been previously investigated, results remain inconclusive. METHODS: To further evaluate the involvement of oncogenic HPV types in CRC, 40 frozen neoplastic and 40 adjacent colonic tissues collected from Italian patients were analyzed by Luminex-based assays that detect a broad spectrum of HPV types, i.e. Alpha (n = 21), Beta (n = 46) and Gamma HPVs (n = 52). In addition, 125 frozen CRC samples and 70 surrounding mucosal tissues were collected from Czech patients and analyzed by broad spectrum PCR protocols: (i) FAP59/64, (ii) FAPM1 and (iii) CUT combined with Next Generation Sequencing (NGS). RESULTS: Using Luminex-basedassays, DNA from HPV16 was detected in 5% (2/40) CRC tissues from Italian patients. One HPV16 DNA-positive CRC case was subsequently confirmed positive for E6*I mRNA. Cutaneous beta HPV types were detected in 10% (4/40) adjacent tissues only, namely HPV111 (n = 3) and HPV120 (n = 1), while gamma HPV168 (n = 1) and HPV199 (n = 1) types were detected in adjacent and in tumor tissues, respectively. The NGS analysis of the CRC Czech samples identified HPV sequences from mucosal alpha-3 (HPV89), alpha-7 (HPV18, 39, 68 and 70) and alpha-10 species (HPV11), as well as cutaneous beta-1 (HPV20, 24, 93, 98, 105,124) beta-2 (HPV23), beta-3 (HPV49) and gamma-1 species (HPV205). CONCLUSIONS: Our findings indicate that HPV types belonging to the mucosal alpha, and the 'cutaneous' beta and gamma genera can be detected in the colonic mucosal samples with a low prevalence rate and a low number of HPV reads by Luminex and NGS, respectively. However, additional studies are required to corroborate these findings.

Prognostic potential of whole exome sequencing in the clinical management of metachronous colorectal cancer liver metastases.

BACKGROUND: Colorectal cancer is a highly prevalent and deadly. The most common metastatic site is the liver. We performed a whole exome sequencing analysis of a series of metachronous colorectal cancer liver metastases (mCLM) and matched non-malignant liver tissues to investigate the genomic profile of mCLM and explore associations with the patients' prognosis and therapeutic modalities. METHODS: DNA samples from mCLM and non-malignant liver tissue pairs (n = 41) were sequenced using whole exome target enrichment and their germline and somatic genetic variability, copy number variations, and mutational signatures were assessed for associations with relapse-free (RFS) and overall survival (OS). RESULTS: Our genetic analysis could stratify all patients into existing targeted therapeutic regimens. The most commonly mutated genes in mCLM were TP53, APC, and KRAS together with PIK3CA and several passenger genes like ABCA13, FAT4, PCLO, and UNC80. Patients with somatic alterations in genes from homologous recombination repair, Notch, and Hedgehog pathways had significantly prolonged RFS, while those with altered MYC pathway genes had poor RFS. Additionally, alterations in the JAK-STAT pathway were prognostic of longer OS. Patients bearing somatic variants in VIPR2 had significantly shorter OS and those with alterations in MUC16 prolonged OS. Carriage of the KRAS-12D variant was associated with shortened survival in our and external datasets. On the other hand, tumor mutation burden, mismatch repair deficiency, microsatellite instability, mutational signatures, or copy number variation in mCLM had no prognostic value. CONCLUSIONS: The results encourage further molecular profiling for personalized treatment of colorectal cancer liver metastases discerning metachronous from synchronous scenarios.

Protein expression and localization of ABC transporters in pancreatic adenocarcinoma: Prognostic role of ABCC8.

BACKGROUND: ATP-binding cassette (ABC) transporters translocate various substances across cellular membranes. Their deregulation may cause cancer drug resistance or perturbations in the supply of building blocks for cancer cells and modify patients' prognosis. This study investigated protein expression and cellular localization of the previously suggested putative prognostic biomarkers - ABCB2/TAP1, ABCC7/CFTR, ABCC8/SUR1, and ABCD4 in patients with pancreatic ductal adenocarcinoma (PDAC). METHODS: Protein expression and localization were assessed by immunohistochemistry in formalin-fixed paraffin-embedded primary tumor tissue blocks of 61 PDAC patients and associated with clinical data and the survival of patients. RESULTS: No CFTR protein expression was observed in PDAC, while TAP1 and ABCC8 were expressed predominantly in the cytoplasm of tumor cells. Most samples (81 %) had detectable both membranous and cytoplasmic ABCD4 staining and 42 % had ABCD4 expressed in the apical orientation. Negative membranous ABCD4 staining was significantly more frequent in advanced stage III or IV tumors (p = 0.022). Small or medium counts of individual ABCC8-positive cells in the stroma surrounding tumor tubules were also more often found in stage III or IV (p = 0.044). Patients with moderate or strong ABCC8 cytoplasmic staining intensity in tumor cells had a 3.5-fold higher risk of disease progression than those with weak staining (p = 0.002). CONCLUSIONS: The study shows for the first time that the cytoplasmic ABCC8 protein expression has prognostic value in PDAC.

The potential of exome sequencing of paired colorectal tumors and synchronous liver metastases to improve treatment.

Tweetable abstract Sequencing exomes of synchronous and metachronous liver metastases of colorectal cancer has potential to enhance treatment, since they can have molecular profiles distinct from primary tumors.

Association of lncRNA and transcriptome intersections with response to targeted therapy in metastatic renal cell carcinoma.

Long non-coding RNAs (lncRNAs) serve an important role in cancer progression and may be used as efficient molecular biomarkers. The present study aimed to identify lncRNAs associated with the response to the receptor tyrosine kinase inhibitor sunitinib and transcriptome profile and clinical features of metastatic renal cell carcinoma (mRCC). The gene expression of 84 cancer-associated lncRNAs in tumor and non-malignant tissue samples of 38 patients with mRCC was evaluated using quantitative PCR. In addition, the coding transcriptome was estimated using RNA sequencing in a subgroup of 20 patients and mRNA-lncRNA intersections were identified. In total, 37 and 13 lncRNAs were down- and upregulated, respectively, in tumor compared with non-malignant adjacent tissue samples. A total of 10 and 4 lncRNAs were up- and downregulated, respectively, in good responders to sunitinib compared with poor responders. High expression of HNF1A-AS1 and IPW lncRNAs was associated with prolonged progression-free survival of patients and a high expression of the TUSC7 lncRNA was associated with poor response and worse survival. Significant associations of dysregulated MEG3 and SNHG16 lncRNAs with expression of protein-coding genes representing various pathways, were identified. Furthermore, a significantly higher expression of CLIP4 gene was observed in good responders. The present study revealed promising candidates for predictive and prognostic biomarkers with further therapeutic potential.

Exploring the Neandertal legacy of pancreatic ductal adenocarcinoma risk in Eurasians.

BACKGROUND: The genomes of present-day non-Africans are composed of 1-3% of Neandertal-derived DNA as a consequence of admixture events between Neandertals and anatomically modern humans about 50-60 thousand years ago. Neandertal-introgressed single nucleotide polymorphisms (aSNPs) have been associated with modern human disease-related traits, which are risk factors for pancreatic ductal adenocarcinoma (PDAC), such as obesity, type 2 diabetes, and inflammation. In this study, we aimed at investigating the role of aSNPs in PDAC in three Eurasian populations. RESULTS: The high-coverage Vindija Neandertal genome was used to select aSNPs in non-African populations from 1000 Genomes project phase 3 data. Then, the association between aSNPs and PDAC risk was tested independently in Europeans and East Asians, using existing GWAS data on more than 200 000 individuals. We did not find any significant associations between aSNPs and PDAC in samples of European descent, whereas, in East Asians, we observed that the Chr10p12.1-rs117585753-T allele (MAF = 10%) increased the risk to develop PDAC (OR = 1.35, 95%CI 1.19-1.54, P = 3.59 × 10-6), with a P-value close to a threshold that takes into account multiple testing. CONCLUSIONS: Our results show only a minimal contribution of Neandertal SNPs to PDAC risk.

Integrative analysis of mRNA and miRNA expression profiles and somatic variants in oxysterol signaling in early-stage luminal breast cancer.

Oxysterols, oxidized derivatives of cholesterol, act in breast cancer (BC) as selective estrogen receptor modulators and affect cholesterol homeostasis, drug transport, nuclear and cell receptors, and other signaling proteins. Using data from three highly overlapping sets of patients (N = 162 in total) with early-stage estrogen-receptor-positive luminal BC-high-coverage targeted DNA sequencing (113 genes), mRNA sequencing, and full micro-RNA (miRNA) transcriptome microarrays-we describe complex oxysterol-related interaction (correlation) networks, with validation in public datasets (n = 538) and 11 databases. The ESR1-CH25H-INSIG1-ABCA9 axis was the most prominent, interconnected through miR-125b-5p, miR-99a-5p, miR-100-5p, miR-143-3p, miR-199b-5p, miR-376a-3p, and miR-376c-3p. Mutations in SC5D, CYP46A1, and its functionally linked gene set were associated with multiple differentially expressed oxysterol-related genes. STARD5 was upregulated in patients with positive lymph node status. High expression of hsa-miR-19b-3p was weakly associated with poor survival. This is the first study of oxysterol-related genes in BC that combines DNA, mRNA, and miRNA multiomics with detailed clinical data. Future studies should provide links between intratumoral oxysterol signaling depicted here, circulating oxysterol levels, and therapy outcomes, enabling eventual clinical exploitation of present findings.

Effects of 7-ketocholesterol on tamoxifen efficacy in breast carcinoma cell line models in vitro.

Oxysterols play significant roles in many physiological and pathological processes including cancer. They modulate some of the cancer hallmarks pathways, influence the efficacy of anti-cancer drugs, and associate with patient survival. In this study, we aimed to analyze the role of 7-ketocholesterol (7-KC) in breast carcinoma cells and its potential modulation of the tamoxifen effect. 7-KC effects were studied in two estrogen receptor (ER)-positive (MCF-7 and T47D) and one ER-negative (BT-20) breast cancer cell lines. First, we tested the viability of cells in the presence of 7-KC. Next, we co-incubated cells with tamoxifen and sublethal concentrations of 7-KC. We also tested changes in caspase 3/7 activity, deregulation of the cell cycle, and changes in expression of selected genes/proteins in the presence of tamoxifen, 7-KC, or their combination. Finally, we analyzed the effect of 7-KC on cellular migration and invasion. We found that the presence of 7-KC slightly decreases the efficacy of tamoxifen in MCF-7 cells, while an increased effect of tamoxifen and higher caspase 3/7 activity was observed in the BT-20 cell line. In the T47D cell line, we did not find any modulation of tamoxifen efficacy by the presence of 7-KC. Expression analysis showed the deregulation in CYP1A1 and CYP1B1 with the opposite trend in MCF-7 and BT-20 cells. Moreover, 7-KC increased cellular migration and invasion potential regardless of the ER status. This study shows that 7-KC can modulate tamoxifen efficacy as well as cellular migration and invasion, making 7-KC a promising candidate for future studies.

DNA Repair Pathway in Ovarian Cancer Patients Treated with HIPEC.

DNA repair pathways are essential for maintaining genome stability, and understanding the regulation of these mechanisms may help in the design of new strategies for treatments, the prevention of platinum-based chemoresistance, and the prolongation of overall patient survival not only with respect to ovarian cancer. The role of hyperthermic intraperitoneal chemotherapy (HIPEC) together with cytoreductive surgery (CRS) and adjuvant systemic chemotherapy is receiving more interest in ovarian cancer (OC) treatment because of the typical peritoneal spread of the disease. The aim of our study was to compare the expression level of 84 genes involved in the DNA repair pathway in tumors and the paired peritoneal metastasis tissue of patients treated with CRS/platinum-based HIPEC with respect to overall patient survival, presence of peritoneal carcinomatosis, treatment response, and alterations in the BRCA1 and BRCA2 genes. Tumors and metastatic tissue from 28 ovarian cancer patients collected during cytoreductive surgery before HIPEC with cisplatin were used for RNA isolation and subsequent cDNA synthesis. Quantitative real-time PCR followed. The most interesting findings of our study are undoubtedly the gene interactions among the genes CCNH, XPA, SLK, RAD51C, XPA, NEIL1, and ATR for primary tumor tissue and ATM, ATR, BRCA2, CDK7, MSH2, MUTYH, POLB, and XRCC4 for metastases. Another interesting finding is the correlation between gene expression and overall survival (OS), where a low expression correlates with a worse OS.

Pathobionts in the tumour microbiota predict survival following resection for colorectal cancer.

BACKGROUND AND AIMS: The gut microbiota is implicated in the pathogenesis of colorectal cancer (CRC). We aimed to map the CRC mucosal microbiota and metabolome and define the influence of the tumoral microbiota on oncological outcomes. METHODS: A multicentre, prospective observational study was conducted of CRC patients undergoing primary surgical resection in the UK (n = 74) and Czech Republic (n = 61). Analysis was performed using metataxonomics, ultra-performance liquid chromatography-mass spectrometry (UPLC-MS), targeted bacterial qPCR and tumour exome sequencing. Hierarchical clustering accounting for clinical and oncological covariates was performed to identify clusters of bacteria and metabolites linked to CRC. Cox proportional hazards regression was used to ascertain clusters associated with disease-free survival over median follow-up of 50 months. RESULTS: Thirteen mucosal microbiota clusters were identified, of which five were significantly different between tumour and paired normal mucosa. Cluster 7, containing the pathobionts Fusobacterium nucleatum and Granulicatella adiacens, was strongly associated with CRC (PFDR = 0.0002). Additionally, tumoral dominance of cluster 7 independently predicted favourable disease-free survival (adjusted p = 0.031). Cluster 1, containing Faecalibacterium prausnitzii and Ruminococcus gnavus, was negatively associated with cancer (PFDR = 0.0009), and abundance was independently predictive of worse disease-free survival (adjusted p = 0.0009). UPLC-MS analysis revealed two major metabolic (Met) clusters. Met 1, composed of medium chain (MCFA), long-chain (LCFA) and very long-chain (VLCFA) fatty acid species, ceramides and lysophospholipids, was negatively associated with CRC (PFDR = 2.61 × 10-11); Met 2, composed of phosphatidylcholine species, nucleosides and amino acids, was strongly associated with CRC (PFDR = 1.30 × 10-12), but metabolite clusters were not associated with disease-free survival (p = 0.358). An association was identified between Met 1 and DNA mismatch-repair deficiency (p = 0.005). FBXW7 mutations were only found in cancers predominant in microbiota cluster 7. CONCLUSIONS: Networks of pathobionts in the tumour mucosal niche are associated with tumour mutation and metabolic subtypes and predict favourable outcome following CRC resection. Video Abstract.

Polymorphic variants involved in methylation regulation: a strategy to discover risk loci for pancreatic ductal adenocarcinoma.

INTRODUCTION: Only a small number of risk factors for pancreatic ductal adenocarcinoma (PDAC) has been established. Several studies identified a role of epigenetics and of deregulation of DNA methylation. DNA methylation is variable across a lifetime and in different tissues; nevertheless, its levels can be regulated by genetic variants like methylation quantitative trait loci (mQTLs), which can be used as a surrogate. MATERIALS AND METHODS: We scanned the whole genome for mQTLs and performed an association study in 14 705 PDAC cases and 246 921 controls. The methylation data were obtained from whole blood and pancreatic cancer tissue through online databases. We used the Pancreatic Cancer Cohort Consortium and the Pancreatic Cancer Case-Control Consortium genome-wide association study (GWAS) data as discovery phase and the Pancreatic Disease Research consortium, the FinnGen project and the Japan Pancreatic Cancer Research consortium GWAS as replication phase. RESULTS: The C allele of 15q26.1-rs12905855 showed an association with a decreased risk of PDAC (OR=0.90, 95% CI 0.87 to 0.94, p=4.93×10-8 in the overall meta-analysis), reaching genome-level statistical significance. 15q26.1-rs12905855 decreases the methylation of a 'C-phosphate-G' (CpG) site located in the promoter region of the RCCD1 antisense (RCCD1-AS1) gene which, when expressed, decreases the expression of the RCC1 domain-containing (RCCD1) gene (part of a histone demethylase complex). Thus, it is possible that the rs12905855 C-allele has a protective role in PDAC development through an increase of RCCD1 gene expression, made possible by the inactivity of RCCD1-AS1. CONCLUSION: We identified a novel PDAC risk locus which modulates cancer risk by controlling gene expression through DNA methylation.

Exome Sequencing of Paired Colorectal Carcinomas and Synchronous Liver Metastases for Prognosis and Therapy Prediction.

PURPOSE: Analysis of somatic variant profiles in retrospectively collected pairs of primary tumors and synchronous liver metastases from surgically treated patients with colorectal carcinomas. Mutational profiles were compared between groups of patients stratified by response to chemotherapy and survival. PATIENTS AND METHODS: The study used whole-exome sequencing of tumor sample pairs from 20 patients diagnosed and treated at a single center. The Cancer Genome Atlas COAD-READ data set (n = 380) was used for validation in silico, where possible. RESULTS: The most frequently altered oncodrivers were APC (55% in primaries and 60% in metastases), TP53 (50/45), TRIP11 (30/5), FAT4 (20/25), and KRAS (15/25). Harboring variants with a high or moderate predicted functional effect in KRAS in primary tumors was significantly associated with poor relapse-free survival in both our sample set and the validation data set. We found a number of additional prognostic associations, including mutational load, alterations in individual genes, oncodriver pathways, and single base substitution (SBS) signatures in primary tissues, which were not confirmed by validation. Altered ATM, DNAH11, and MUC5AC, or a higher share of SBS24 signature in metastases seemed to represent poor prognostic factors, but because of a lack of suitable validation data sets, these results must be treated with extreme caution. No gene or profile was significantly associated with response to chemotherapy. CONCLUSION: Taken together, we report subtle differences in exome mutational profiles between paired primary tumors and synchronous liver metastases and a distinct prognostic relevance of KRAS in primary tumors. Although the general scarcity of primary tumor-synchronous metastasis sample pairs with high-quality clinical data makes robust validation difficult, this study provides potentially valuable data for utilization in precision oncology and could serve as a springboard for larger studies.