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.

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.

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.

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.

Genetic variations in microRNA-binding sites of solute carrier transporter genes as predictors of clinical outcome in colorectal cancer.

One of the principal mechanisms of chemotherapy resistance in highly frequent solid tumors, such as colorectal cancer (CRC), is the decreased activity of drug transport into tumor cells due to low expression of important membrane proteins, such as solute carrier (SLC) transporters. Sequence complementarity is a major determinant for target gene recognition by microRNAs (miRNAs). Single-nucleotide polymorphisms (SNPs) in target sequences transcribed into messenger RNA may therefore alter miRNA binding to these regions by either creating a new site or destroying an existing one. miRSNPs may explain the modulation of expression levels in association with increased/decreased susceptibility to common diseases as well as in chemoresistance and the consequent inter-individual variability in drug response. In the present study, we investigated whether miRSNPs in SLC transporter genes may modulate CRC susceptibility and patient's survival. Using an in silico approach for functional predictions, we analyzed 26 miRSNPs in 9 SLC genes in a cohort of 1368 CRC cases and 698 controls from the Czech Republic. After correcting for multiple tests, we found several miRSNPs significantly associated with patient's survival. SNPs in SLCO3A1, SLC22A2 and SLC22A3 genes were defined as prognostic factors in the classification and regression tree analysis. In contrast, we did not observe any significant association between miRSNPs and CRC risk. To the best of our knowledge, this is the first study investigating miRSNPs potentially affecting miRNA binding to SLC transporter genes and their impact on CRC susceptibility or patient's prognosis.

Novel morphological multi-scale evaluation system for quality assessment of decellularized liver scaffolds.

Decellularized scaffolds can serve as an excellent three-dimensional environment for cell repopulation. They maintain tissue-specific microarchitecture of extracellular matrix proteins with important spatial cues for cell adhesion, migration, growth, and differentiation. However, criteria for quality assessment of the three-dimensional structure of decellularized scaffolds are rather fragmented, usually study-specific, and mostly semi-quantitative. Thus, we aimed to develop a robust structural assessment system for decellularized porcine liver scaffolds. Five scaffolds of different quality were used to establish the new evaluation system. We combined conventional semi-quantitative scoring criteria with a quantitative scaffold evaluation based on automated image analysis. For the quantitation, we developed a specific open source software tool (ScaffAn) applying algorithms designed for texture analysis, segmentation, and skeletonization. ScaffAn calculates selected parameters characterizing structural features of porcine liver scaffolds such as the sinusoidal network. After evaluating individual scaffolds, the total scores predicted scaffold interaction with cells in terms of cell adhesion. Higher scores corresponded to higher numbers of cells attached to the scaffolds. Moreover, our analysis revealed that the conventional system could not identify fine differences between good quality scaffolds while the additional use of ScaffAn allowed discrimination. This led us to the conclusion that only using the combined score resulted in the best discrimination between different quality scaffolds. Overall, our newly defined evaluation system has the potential to select the liver scaffolds most suitable for recellularization, and can represent a step toward better success in liver tissue engineering.

Expression quantitative trait loci in ABC transporters are associated with survival in 5-FU treated colorectal cancer patients.

The chemotherapeutic efficacy in colorectal cancer (CRC) is limited due to the inter-individual variability in drug response and the development of tumour resistance. ATP-binding cassette (ABC) transporters are crucial in the development of resistance by the efflux of anticancer agents from cancer cells. In this study, we identified 14 single nucleotide polymorphisms (SNPs) in 11 ABC transporter genes acting as an expression of quantitative trait loci (eQTLs), i.e. whose variation influence the expression of many downstream genes. These SNPs were genotyped in a case-control study comprising 1098 cases and 1442 healthy controls and analysed in relation to CRC development risk and patient survival. Considering a strict correction for multiple tests, we did not observe any significant association between SNPs and CRC risk. The rs3819720 polymorphism in the ABCB3/TAP2 gene was statistically significantly associated with shorter overall survival (OS) in the codominant, and dominant models [GA vs. GG, hazard ratio (HR) = 1.48; P = 0.002; AA vs. GG, HR = 1.70; P = 0.004 and GA + AA vs. GG, HR = 1.52; P = 0.0006]. Additionally, GA carriers of the same SNP displayed worse OS after receiving 5-FU based chemotherapy. The variant allele of rs3819720 polymorphism statistically significantly affected the expression of 36 downstream genes. Screening for eQTL polymorphisms in relevant genes such as ABC transporters that can regulate the expression of several other genes may help to identify the genetic background involved in the individual response to the treatment of CRC patients.

Chromosomal aberrations in tire plant workers and interaction with polymorphisms of biotransformation and DNA repair genes.

We evaluated chromosomal aberrations in lymphocytes of 177 workers exposed to xenobiotics in a tire plant and in 172 controls, in relation to their genetic background. Nine polymorphisms in genes encoding biotransformation enzymes and nine polymorphisms in genes involved in main DNA repair pathways were investigated for possible modulation of chromosomal damage. Chromosomal aberration frequencies were the highest among exposed smokers and the lowest in non-smoking unexposed individuals (2.5+/-1.8% vs. 1.7+/-1.2%, respectively). The differences between groups (ANOVA) were borderline significant (F=2.6, P=0.055). Chromosomal aberrations were higher in subjects with GSTT1-null (2.4+/-1.7%) than in those with GSTT1-plus genotype (1.8+/-1.4%; F=7.2, P=0.008). Considering individual groups, this association was significant in smoking exposed workers (F=4.4, P=0.040). Individuals with low activity EPHX1 genotype exhibited significantly higher chromosomal aberrations (2.3+/-1.6%) in comparison with those bearing medium (1.7+/-1.2%) and high activity genotype (1.5+/-1.2%; F=4.7, P=0.010). Both chromatid- and chromosome-type aberration frequencies were mainly affected by exposure and smoking status. Binary logistic regression analysis revealed that frequencies of chromatid-type aberrations were modulated by NBS1 Glu185Gln (OR 4.26, 95%CI 1.38-13.14, P=0.012), and to a moderate extent, by XPD Lys751Gln (OR 0.16, 95%CI 0.02-1.25, P=0.081) polymorphisms. Chromosome-type aberrations were lowest in individuals bearing the EPHX1 genotype conferring the high activity (OR 0.38, 95%CI 0.15-0.98, P=0.045). Present results show that exposed individuals in the tire production, who smoke, exhibit higher chromosomal aberrations frequencies, and the extent of chromosomal damage may additionally be modified by relevant polymorphisms.

Single nucleotide polymorphism in 5'-flanking region of BCL6 is not associated with increased risk of non-Hodgkin's lymphoma.

Mutations in the 5'-regulatory region of BCL6 were suggested to play a role in non-Hodgkin's lymphoma (NHL) progression and in the transformation of follicular lymphoma to more aggressive diffuse large B-cell type. The aim of this study was to explore association between polymorphism G397C in the 5'-region of BCL6 and both incidence and progression of NHL in 154 NHL cases and 362 controls. Neither frequencies of the rare BCL6 allele 397C nor particular genotypes differed significantly between NHL cases and controls. There was no significant association of histological type of NHL and clinical characteristics with this polymorphism.

Breast cancer: role of polymorphisms in biotransformation enzymes.

We aimed at determining whether any association exists between genetic polymorphisms in epoxide hydrolase (EPHX1), NADPH-quinone oxidoreductase (NQO1), glutathione S-transferases (GSTM1/P1/T1) and individual susceptibility to breast cancer. Polymerase chain reaction-restriction fragment length polymorphism-based genotyping assays were used to determine the frequency of polymorphisms in EPHX1 (exons 3 and 4), NQO1 (exon 6), GSTM1 (deletion), GSTP1 (exon 5), and GSTT1 (deletion) in a case-control study comprised of 238 patients with breast cancer and 313 healthy individuals. The distribution of genotypes in exon 6 of NQO1 was significantly different between the control group and breast cancer cases. Age-adjusted odds ratio (OR) for variant genotype NQO1*2/*2 was 3.68 (confidence interval (CI) = 1.41-9.62, P = 0.008). Association of GSTP1*2/*2 genotype as well as that of low EPHX1 activity deduced by combinations of genotypes in exons 3 and 4 with breast cancer was suggestive, but nonsignificant. Individuals simultaneously lacking GSTM1 and carrying at least one GSTP1 variant allele were at significantly higher risk of breast cancer (OR = 2.03, CI = 1.18-3.50, P = 0.010). Combinations of either GSTM1null or GSTP1*2 with low activity of EPHX1 presented significant risk of breast cancer (OR = 1.88, CI = 1.00-3.52, P = 0.049 and OR = 2.40, CI = 1.15-5.00, P = 0.019, respectively) as well. In conclusion, the results suggest that genetic polymorphisms in biotransformation enzymes may play a significant role in the development of breast cancer.