Contribution of pks+ E. coli mutations to colorectal carcinogenesis.

The dominant mutational signature in colorectal cancer genomes is C > T deamination (COSMIC Signature 1) and, in a small subgroup, mismatch repair signature (COSMIC signatures 6 and 44). Mutations in common colorectal cancer driver genes are often not consistent with those signatures. Here we perform whole-genome sequencing of normal colon crypts from cancer patients, matched to a previous multi-omic tumour dataset. We analyse normal crypts that were distant vs adjacent to the cancer. In contrast to healthy individuals, normal crypts of colon cancer patients have a high incidence of pks + (polyketide synthases) E.coli (Escherichia coli) mutational and indel signatures, and this is confirmed by metagenomics. These signatures are compatible with many clonal driver mutations detected in the corresponding cancer samples, including in chromatin modifier genes, supporting their role in early tumourigenesis. These results provide evidence that pks + E.coli is a potential driver of carcinogenesis in the human gut.

The co-evolution of the genome and epigenome in colorectal cancer.

Colorectal malignancies are a leading cause of cancer-related death1 and have undergone extensive genomic study2,3. However, DNA mutations alone do not fully explain malignant transformation4-7. Here we investigate the co-evolution of the genome and epigenome of colorectal tumours at single-clone resolution using spatial multi-omic profiling of individual glands. We collected 1,370 samples from 30 primary cancers and 8 concomitant adenomas and generated 1,207 chromatin accessibility profiles, 527 whole genomes and 297 whole transcriptomes. We found positive selection for DNA mutations in chromatin modifier genes and recurrent somatic chromatin accessibility alterations, including in regulatory regions of cancer driver genes that were otherwise devoid of genetic mutations. Genome-wide alterations in accessibility for transcription factor binding involved CTCF, downregulation of interferon and increased accessibility for SOX and HOX transcription factor families, suggesting the involvement of developmental genes during tumourigenesis. Somatic chromatin accessibility alterations were heritable and distinguished adenomas from cancers. Mutational signature analysis showed that the epigenome in turn influences the accumulation of DNA mutations. This study provides a map of genetic and epigenetic tumour heterogeneity, with fundamental implications for understanding colorectal cancer biology.

Phenotypic plasticity and genetic control in colorectal cancer evolution.

Genetic and epigenetic variation, together with transcriptional plasticity, contribute to intratumour heterogeneity1. The interplay of these biological processes and their respective contributions to tumour evolution remain unknown. Here we show that intratumour genetic ancestry only infrequently affects gene expression traits and subclonal evolution in colorectal cancer (CRC). Using spatially resolved paired whole-genome and transcriptome sequencing, we find that the majority of intratumour variation in gene expression is not strongly heritable but rather 'plastic'. Somatic expression quantitative trait loci analysis identified a number of putative genetic controls of expression by cis-acting coding and non-coding mutations, the majority of which were clonal within a tumour, alongside frequent structural alterations. Consistently, computational inference on the spatial patterning of tumour phylogenies finds that a considerable proportion of CRCs did not show evidence of subclonal selection, with only a subset of putative genetic drivers associated with subclone expansions. Spatial intermixing of clones is common, with some tumours growing exponentially and others only at the periphery. Together, our data suggest that most genetic intratumour variation in CRC has no major phenotypic consequence and that transcriptional plasticity is, instead, widespread within a tumour.

Genetic susceptibility to telomere shortening through the rs2293607 polymorphism is associated with a greater risk of alcohol use disorder.

Telomere shortening is usually considered a biomarker of ageing. Harmful alcohol use promotes accelerated biological ageing and alcohol use disorders (AUDs) are associated with short telomere length (TL). This study was conducted to examine the relationship of TL to AUD and determine whether single nucleotide polymorphisms (SNPs) in TERC and TERT modulate this association. For this purpose, we genotyped TERC SNPs rs2293607, rs12696304, and rs16847897 and TERT SNPs rs2735940, rs2736100, and rs2736098 in 308 male patients with AUD and 255 sex-matched healthy controls and measured TL in a subset of 99 patients and 99 controls paired by age and smoking status. Our results showed that the mean TL was shorter in patients with AUD than in controls. The area under the ROC curve was 0.70 (P < 0.001). The GG genotype of TERC rs2293607 was more common among patients with AUD than among controls (9.8% vs. 5.1%; P = 0.038). No difference was found for the other SNPs. Carriers of the GG genotype of rs2293607 had shorter telomeres than did allele A carriers. In conclusion, patients with AUD had shorter telomeres. Genetic susceptibility to telomere shortening through the rs2293607 SNP is associated with a greater risk of AUD.

EGFR amplification and outcome in a randomised phase III trial of chemotherapy alone or chemotherapy plus panitumumab for advanced gastro-oesophageal cancers.

OBJECTIVE: Epidermal growth factor receptor (EGFR) inhibition may be effective in biomarker-selected populations of advanced gastro-oesophageal adenocarcinoma (aGEA) patients. Here, we tested the association between outcome and EGFR copy number (CN) in pretreatment tissue and plasma cell-free DNA (cfDNA) of patients enrolled in a randomised first-line phase III clinical trial of chemotherapy or chemotherapy plus the anti-EGFR monoclonal antibody panitumumab in aGEA (NCT00824785). DESIGN: EGFR CN by either fluorescence in situ hybridisation (n=114) or digital-droplet PCR in tissues (n=250) and plasma cfDNAs (n=354) was available for 474 (86%) patients in the intention-to-treat (ITT) population. Tissue and plasma low-pass whole-genome sequencing was used to screen for coamplifications in receptor tyrosine kinases. Interaction between chemotherapy and EGFR inhibitors was modelled in patient-derived organoids (PDOs) from aGEA patients. RESULTS: EGFR amplification in cfDNA correlated with poor survival in the ITT population and similar trends were observed when the analysis was conducted in tissue and plasma by treatment arm. EGFR inhibition in combination with chemotherapy did not correlate with improved survival, even in patients with significant EGFR CN gains. Addition of anti-EGFR inhibitors to the chemotherapy agent epirubicin in PDOs, resulted in a paradoxical increase in viability and accelerated progression through the cell cycle, associated with p21 and cyclin B1 downregulation and cyclin E1 upregulation, selectively in organoids from EGFR-amplified aGEA. CONCLUSION: EGFR CN can be accurately measured in tissue and liquid biopsies and may be used for the selection of aGEA patients. EGFR inhibitors may antagonise the antitumour effect of anthracyclines with important implications for the design of future combinatorial trials.

Oncogenic driver mutations predict outcome in a cohort of head and neck squamous cell carcinoma (HNSCC) patients within a clinical trial.

234 diagnostic formalin-fixed paraffin-embedded (FFPE) blocks from homogeneously treated patients with locally advanced head and neck squamous cell carcinoma (HNSCC) within a multicentre phase III clinical trial were characterised. The mutational spectrum was examined by next generation sequencing in the 26 most frequent oncogenic drivers in cancer and correlated with treatment response and survival. Human papillomavirus (HPV) status was measured by p16INK4a immunohistochemistry in oropharyngeal tumours. Clinicopathological features and response to treatment were measured and compared with the sequencing results. The results indicated TP53 as the most mutated gene in locally advanced HNSCC. HPV-positive oropharyngeal tumours were less mutated than HPV-negative tumours in TP53 (p < 0.01). Mutational and HPV status influences patient survival, being mutated or HPV-negative tumours associated with poor overall survival (p < 0.05). No association was found between mutations and clinicopathological features. This study confirmed and expanded previously published genomic characterization data in HNSCC. Survival analysis showed that non-mutated HNSCC tumours associated with better prognosis and lack of mutations can be identified as an important biomarker in HNSCC. Frequent alterations in PI3K pathway in HPV-positive HNSCC could define a promising pathway for pharmacological intervention in this group of tumours.

Exploiting evolutionary steering to induce collateral drug sensitivity in cancer.

Drug resistance mediated by clonal evolution is arguably the biggest problem in cancer therapy today. However, evolving resistance to one drug may come at a cost of decreased fecundity or increased sensitivity to another drug. These evolutionary trade-offs can be exploited using 'evolutionary steering' to control the tumour population and delay resistance. However, recapitulating cancer evolutionary dynamics experimentally remains challenging. Here, we present an approach for evolutionary steering based on a combination of single-cell barcoding, large populations of 108-109 cells grown without re-plating, longitudinal non-destructive monitoring of cancer clones, and mathematical modelling of tumour evolution. We demonstrate evolutionary steering in a lung cancer model, showing that it shifts the clonal composition of the tumour in our favour, leading to collateral sensitivity and proliferative costs. Genomic profiling revealed some of the mechanisms that drive evolved sensitivity. This approach allows modelling evolutionary steering strategies that can potentially control treatment resistance.

Measuring single cell divisions in human tissues from multi-region sequencing data.

Both normal tissue development and cancer growth are driven by a branching process of cell division and mutation accumulation that leads to intra-tissue genetic heterogeneity. However, quantifying somatic evolution in humans remains challenging. Here, we show that multi-sample genomic data from a single time point of normal and cancer tissues contains information on single-cell divisions. We present a new theoretical framework that, applied to whole-genome sequencing data of healthy tissue and cancer, allows inferring the mutation rate and the cell survival/death rate per division. On average, we found that cells accumulate 1.14 mutations per cell division in healthy haematopoiesis and 1.37 mutations per division in brain development. In both tissues, cell survival was maximal during early development. Analysis of 131 biopsies from 16 tumours showed 4 to 100 times increased mutation rates compared to healthy development and substantial inter-patient variation of cell survival/death rates.

Genetic Susceptibility in Head and Neck Squamous Cell Carcinoma in a Spanish Population.

Despite classical environmental risk factors like tobacco, alcohol or viral infection, not all individuals develop head and neck cancer. Therefore, identification of the genetic susceptibility produced by single nucleotide polymorphisms (SNPs) is an important task. A total of 296 human papillomavirus negative head and neck cancer (HNC) patients (126 laryngeal, 100 pharyngeal and 70 oral cavity) were included in the study, involving 29 candidate SNPs in genes within important carcinogenic pathways (oncogenesis and tumour suppression, DNA repair, inflammation, oxidation and apoptosis). Genotyping was performed using TaqMan probes or restriction fragment length assays in peripheral blood DNA. In addition, 259 paired controls were also evaluated with the same risk factors for each specific location. Nine SNPs in DNA repair (ERCC1 rs11615, ERCC2 rs13181), inflammatory (IL2 rs2069762, IL6 rs1800795), oxidative (NFE2L2 rs13035806 and rs2706110) and apoptotic genes (TP53 rs1042522, MDM2 rs2279744, BCL2 rs2279115) were differently associated with HNSCC susceptibility by location. Some of these SNPs were not described before in this tumour type. In conclusion, we describe several SNPs associated with HNC in a Spanish population.

Corrigendum: Association of VAV2 and VAV3 polymorphisms with cardiovascular risk factors.

This corrects the article DOI: 10.1038/srep41875.

Patient-derived organoids model treatment response of metastatic gastrointestinal cancers.

Patient-derived organoids (PDOs) have recently emerged as robust preclinical models; however, their potential to predict clinical outcomes in patients has remained unclear. We report on a living biobank of PDOs from metastatic, heavily pretreated colorectal and gastroesophageal cancer patients recruited in phase 1/2 clinical trials. Phenotypic and genotypic profiling of PDOs showed a high degree of similarity to the original patient tumors. Molecular profiling of tumor organoids was matched to drug-screening results, suggesting that PDOs could complement existing approaches in defining cancer vulnerabilities and improving treatment responses. We compared responses to anticancer agents ex vivo in organoids and PDO-based orthotopic mouse tumor xenograft models with the responses of the patients in clinical trials. Our data suggest that PDOs can recapitulate patient responses in the clinic and could be implemented in personalized medicine programs.

Biomarkers and polymorphisms in pancreatic neuroendocrine tumors treated with sunitinib.

Several circulating biomarkers and single nucleotide polymorphisms (SNPs) have been correlated with efficacy and tolerability to antiangiogenic agents. These associations remain unexplored in well-differentiated, metastatic pancreatic neuroendocrine tumors treated with the multitargeted tyrosine kinase inhibitor sunitinib. We have assessed the effect on tumor response at 6 months, overall survival, progression-free survival and safety of 14 SNPs, and 6 soluble proteins. Forty-three patients were recruited. Two SNPs in the vascular endothelial growth factor receptor 3 (VEGFR-3) gene predicted lower overall survival: rs307826 with hazard ratio (HR) 3.67 (confidence interval [CI] 95%, 1.35-10.00) and rs307821 with HR 3.84 (CI 95%, 1.47-10.0). Interleukin-6 was associated with increased mortality: HR 1.06 (CI 95%, 1.01-1.12), and osteopontin was associated with shorter PFS: HR 1.087 (1.01-1.16), independently of Ki-67. Furthermore, levels of osteopontin remained higher at the end of the study in patients considered non-responders: 38.5 ng/mL vs. responders: 18.7 ng/mL, p-value=0.039. Dynamic upward variations were also observed with respect to IL-8 levels in sunitinib-refractory individuals: 28.5 pg/mL at baseline vs. 38.3 pg/mL at 3 months, p-value=0.024. In conclusion, two VEGFR-3 SNPs as well as various serum biomarkers were associated with diverse clinical outcomes in patients with well-differentiated pancreatic neuroendocrine tumors treated with sunitinib.

Analysis of autophagy gene polymorphisms in Spanish patients with head and neck squamous cell carcinoma.

Head and neck squamous cell carcinoma (HNSCC) is the sixth cancer on incidence worldwide. Tobacco and alcohol consumption are the most classical risk factors associated with its development. Autophagy process has a dual effect both in tumourigenesis and tumour suppressing activity. To investigate the importance of this pathway in HNSCC susceptibility, a risk factor matched case-control association study was performed with four candidate polymorphisms in autophagy genes (ATG2B, ATG5, ATG10, ATG16L1). We found an association between the variant in ATG10 rs1864183 and a higher susceptibility to develop laryngeal cancer, ATG2B rs3759601 and pharyngeal cancer and ATG16L1 rs2241880 and oral carcinoma. ATG5 rs2245214 SNP was not associated with any location. Overall, our results indicate the importance of the autophagy pathway in the susceptibility of head and neck squamous cell carcinoma and demonstrate the heterogeneity between its locations encompassed under a single terminology.

Association of VAV2 and VAV3 polymorphisms with cardiovascular risk factors.

Hypertension, diabetes and obesity are cardiovascular risk factors closely associated to the development of renal and cardiovascular target organ damage. VAV2 and VAV3, members of the VAV family proto-oncogenes, are guanosine nucleotide exchange factors for the Rho and Rac GTPase family, which is related with cardiovascular homeostasis. We have analyzed the relationship between the presence of VAV2 rs602990 and VAV3 rs7528153 polymorphisms with cardiovascular risk factors and target organ damage (heart, vessels and kidney) in 411 subjects. Our results show that being carrier of the T allele in VAV2 rs602990 polymorphism is associated with an increased risk of obesity, reduced levels of ankle-brachial index and diastolic blood pressure and reduced retinal artery caliber. In addition, being carrier of T allele is associated with increased risk of target organ damage in males. On the other hand, being carrier of the T allele in VAV3 rs7528153 polymorphism is associated with a decreased susceptibility of developing a pathologic state composed by the presence of hypertension, diabetes, obesity or cardiovascular damage, and with an increased risk of developing altered basal glycaemia. This is the first report showing an association between VAV2 and VAV3 polymorphisms with cardiovascular risk factors and target organ damage.

Correction: Study of Association between Pre-Senile Cataracts and the Polymorphisms rs2228000 in XPC and rs1042522 in p53 in Spanish Population.

[This corrects the article DOI: 10.1371/journal.pone.0156317.].

Study of association between pre-senile cataracts and rs11615 of ERCC1, rs13181 of ERCC2, and rs25487 of XRCC1 polymorphisms in a Spanish population.

BACKGROUND: The purpose of this study was to attempt to determine if the presence of certain polymorphisms in the DNA repair genes (ERCC1, ERCC2, and XRCC1) is associated with pre-senile cataract development. MATERIALS AND METHODS: We performed a retrospective study over three groups of patients. The first group with pre-senile cataract was formed by 72 patients younger than 55 years with cataract surgery. The second group with senile cataract was formed by 101 patients older than 55 years with cataract surgery. And the third group, without cataract, was formed by 42 subjects older than 55 years without lens opacities. We analyzed the presence of SNP rs11615 from ERCC1, rs13181 from ERCC2, and rs25487 from XRCC1 and the relationship between risk factors such as smoking, alcohol intake, hypertension, and diabetes. RESULTS: The comparison of the genotype distribution in ERCC1 and ERCC2 did not show any statistically significant association in any of our analyses (p > 0.05). The comparison of the genotype distribution in XRCC1 within the different groups did not show any statistically significant associations (p > 0.05), except for the comparison between the pre-senile cataract group and the group without cataract, where an increased risk of developing pre-senile cataract for the genotype Gln/Gln (p = 0.029; OR = 1.02-40.67) in recessive inheritance models was observed when adjusting for risk factors. CONCLUSIONS: Allelic variants in ERCC1 and ERCC2 are not associated with an increased risk of developing pre-senile cataract. The presence of Gln/Gln in XRCC1 in the pre-senile cataract group with regard to the group without cataract is associated with a major risk of developing pre-senile cataract.

Study of Association between Pre-Senile Cataracts and the Polymorphisms rs2228000 in XPC and rs1042522 in p53 in Spanish Population.

PURPOSE: To determine if the presence of certain polymorphisms in the DNA repair gene XPC and the apoptosis inductor gene p53 is associated with pre-senile cataract development. METHODS: We have performed a retrospective study over three groups of patients. The group with pre-senile cataract formed by 72 patients younger than 55 with cataract surgery. The group with senile cataract formed by 101 patients older than 55 with cataract surgery. The group without cataract was formed by 42 subjects older than 55 without lens opacities. We analyzed the presence of SNP rs2228000 from XPC and rs1042522 from p53; and the relationship between risk factors such as smoking, alcohol intake, hypertension or diabetes. RESULTS: The comparison of the genotype distribution in XPC, within the different groups, did not show any statistically significant association in any of our analysis (p>0,05). The comparison of the genotype distribution in p53 within the different groups did not show any statistically significant association (p>0,05); except for the comparison between the pre-senile cataract group and the group with senile cataract where the genotype Pro/Pro (C/C) in the recessive inheritance model showed a higher risk for developing pre-senile cataract (p = 0,031; OR = 1.04-15.97). This association decreased when we performed the analysis adjusting by the studied risk factors (p = 0.056). CONCLUSIONS: Allelic variants in the gene XPC are not associated with an increased risk for developing pre-senile cataract. The presence of the genotype Pro/Pro in p53 might be associated with a major risk for developing pre-senile cataract.