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.

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.

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.

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.

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.

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.

Polygenic and multifactorial scores for pancreatic ductal adenocarcinoma risk prediction.

BACKGROUND: Most cases of pancreatic ductal adenocarcinoma (PDAC) are asymptomatic in early stages, and the disease is typically diagnosed in advanced phases, resulting in very high mortality. Tools to identify individuals at high risk of developing PDAC would be useful to improve chances of early detection. OBJECTIVE: We generated a polygenic risk score (PRS) for PDAC risk prediction, combining the effect of known risk SNPs, and carried out an exploratory analysis of a multifactorial score. METHODS: We tested the associations of the individual known risk SNPs on up to 2851 PDAC cases and 4810 controls of European origin from the PANcreatic Disease ReseArch (PANDoRA) consortium. Thirty risk SNPs were included in a PRS, which was computed on the subset of subjects that had 100% call rate, consisting of 839 cases and 2040 controls in PANDoRA and 6420 cases and 4889 controls from the previously published Pancreatic Cancer Cohort Consortium I-III and Pancreatic Cancer Case-Control Consortium genome-wide association studies. Additional exploratory multifactorial scores were constructed by complementing the genetic score with smoking and diabetes. RESULTS: The scores were associated with increased PDAC risk and reached high statistical significance (OR=2.70, 95% CI 1.99 to 3.68, p=2.54×10-10 highest vs lowest quintile of the weighted PRS, and OR=14.37, 95% CI 5.57 to 37.09, p=3.64×10-8, highest vs lowest quintile of the weighted multifactorial score). CONCLUSION: We found a highly significant association between a PRS and PDAC risk, which explains more than individual SNPs and is a step forward in the direction of the construction of a tool for risk stratification in the population.

Laser Capture Microdissection: A Gear for Pancreatic Cancer Research.

The advancement in molecular techniques has been attributed to the quality and significance of cancer research. Pancreatic cancer (PC) is one of the rare cancers with aggressive behavior and a high mortality rate. The asymptomatic nature of the disease until its advanced stage has resulted in late diagnosis as well as poor prognosis. The heterogeneous character of PC has complicated cancer development and progression studies. The analysis of bulk tissues of the disease was insufficient to understand the disease, hence, the introduction of the single-cell separating technique aided researchers to decipher more about the specific cell population of tumors. This review gives an overview of the Laser Capture Microdissection (LCM) technique, one of the single-cell separation methods used in PC research.

Associations between pancreatic expression quantitative traits and risk of pancreatic ductal adenocarcinoma.

Pancreatic ductal adenocarcinoma (PDAC) is among the most lethal cancers. Its poor prognosis is predominantly due to the fact that most patients remain asymptomatic until the disease reaches an advanced stage, alongside the lack of early markers and screening strategies. A better understanding of PDAC risk factors is essential for the identification of groups at high risk in the population. Genome-wide association studies (GWAS) have been a powerful tool for detecting genetic variants associated with complex traits, including pancreatic cancer. By exploiting functional and GWAS data, we investigated the associations between polymorphisms affecting gene function in the pancreas (expression quantitative trait loci, eQTLs) and PDAC risk. In a two-phase approach, we analysed 13 713 PDAC cases and 43 784 controls and identified a genome-wide significant association between the A allele of the rs2035875 polymorphism and increased PDAC risk (P = 7.14 × 10-10). This allele is known to be associated with increased expression in the pancreas of the keratin genes KRT8 and KRT18, whose increased levels have been reported to correlate with various tumour cell characteristics. Additionally, the A allele of the rs789744 variant was associated with decreased risk of developing PDAC (P = 3.56 × 10-6). This single nucleotide polymorphism is situated in the SRGAP1 gene and the A allele is associated with higher expression of the gene, which in turn inactivates the cyclin-dependent protein 42 (CDC42) gene expression, thus decreasing the risk of PDAC. In conclusion, we present here a functional-based novel PDAC risk locus and an additional strong candidate supported by significant associations and plausible biological mechanisms.

Genome-wide scan of long noncoding RNA single nucleotide polymorphisms and pancreatic cancer susceptibility.

Pancreatic ductal adenocarcinoma (PDAC) is projected to become the second cancer-related cause of death by 2030. Identifying novel risk factors, including genetic risk loci, could be instrumental in risk stratification and implementation of prevention strategies. Long noncoding RNAs (lncRNAs) are involved in regulation of key biological processes, and the possible role of their genetic variability has been unexplored so far. Combining genome wide association studies and functional data, we investigated the genetic variability in all lncRNAs. We analyzed 9893 PDAC cases and 9969 controls and identified a genome-wide significant association between the rs7046076 SNP and risk of developing PDAC (P = 9.73 × 10-9 ). This SNP is located in the NONHSAG053086.2 (lnc-SMC2-1) gene and the risk allele is predicted to disrupt the binding of the lncRNA with the micro-RNA (miRNA) hsa-mir-1256 that regulates several genes involved in cell cycle, such as CDKN2B. The CDKN2B region is pleiotropic and its genetic variants have been associated with several human diseases, possibly though an imperfect interaction between lncRNA and miRNA. We present a novel PDAC risk locus, supported by a genome-wide statistical significance and a plausible biological mechanism.

Role of Genetic Variation in Cytochromes P450 in Breast Cancer Prognosis and Therapy Response.

Breast cancer is the most frequent cancer in the female population worldwide. The role of germline genetic variability in cytochromes P450 (CYP) in breast cancer prognosis and individualized therapy awaits detailed elucidation. In the present study, we used the next-generation sequencing to assess associations of germline variants in the coding and regulatory sequences of all human CYP genes with response of the patients to the neoadjuvant cytotoxic chemotherapy and disease-free survival (n = 105). A total of 22 prioritized variants associating with a response or survival in the above evaluation phase were then analyzed by allelic discrimination in the large confirmation set (n = 802). Associations of variants in CYP1B1, CYP4F12, CYP4X1, and TBXAS1 with the response to the neoadjuvant cytotoxic chemotherapy were replicated by the confirmation phase. However, just association of variant rs17102977 in CYP4X1 passed the correction for multiple testing and can be considered clinically and statistically validated. Replicated associations for variants in CYP4X1, CYP24A1, and CYP26B1 with disease-free survival of all patients or patients stratified to subgroups according to therapy type have not passed a false discovery rate test. Although statistically not confirmed by the present study, the role of CYP genes in breast cancer prognosis should not be ruled out. In conclusion, the present study brings replicated association of variant rs17102977 in CYP4X1 with the response of patients to the neoadjuvant cytotoxic chemotherapy and warrants further research of genetic variation CYPs in breast cancer.

Genome-wide association study identifies an early onset pancreatic cancer risk locus.

Early onset pancreatic cancer (EOPC) is a rare disease with a very high mortality rate. Almost nothing is known on the genetic susceptibility of EOPC, therefore, we performed a genome-wide association study (GWAS) to identify novel genetic variants specific for patients diagnosed with pancreatic ductal adenocarcinoma (PDAC) at younger ages. In the first phase, conducted on 821 cases with age of onset ≤60 years, of whom 198 with age of onset ≤50, and 3227 controls from PanScan I-II, we observed four SNPs (rs7155613, rs2328991, rs4891017 and rs12610094) showing an association with EOPC risk (P < 1 × 10-4 ). We replicated these SNPs in the PANcreatic Disease ReseArch (PANDoRA) consortium and used additional in silico data from PanScan III and PanC4. Among these four variants rs2328991 was significant in an independent set of 855 cases with age of onset ≤60 years, of whom 265 with age of onset ≤50, and 4142 controls from the PANDoRA consortium while in the in silico data, we observed no statistically significant association. However, the resulting meta-analysis supported the association (P = 1.15 × 10-4 ). In conclusion, we propose a novel variant rs2328991 to be involved in EOPC risk. Even though it was not possible to find a mechanistic link between the variant and the function, the association is supported by a solid statistical significance obtained in the largest study on EOPC genetics present so far in the literature.

5' Untranslated Region Elements Show High Abundance and Great Variability in Homologous ABCA Subfamily Genes.

The 12 members of the ABCA subfamily in humans are known for their ability to transport cholesterol and its derivatives, vitamins, and xenobiotics across biomembranes. Several ABCA genes are causatively linked to inborn diseases, and the role in cancer progression and metastasis is studied intensively. The regulation of translation initiation is implicated as the major mechanism in the processes of post-transcriptional modifications determining final protein levels. In the current bioinformatics study, we mapped the features of the 5' untranslated regions (5'UTR) known to have the potential to regulate translation, such as the length of 5'UTRs, upstream ATG codons, upstream open-reading frames, introns, RNA G-quadruplex-forming sequences, stem loops, and Kozak consensus motifs, in the DNA sequences of all members of the subfamily. Subsequently, the conservation of the features, correlations among them, ribosome profiling data as well as protein levels in normal human tissues were examined. The 5'UTRs of ABCA genes contain above-average numbers of upstream ATGs, open-reading frames and introns, as well as conserved ones, and these elements probably play important biological roles in this subfamily, unlike RG4s. Although we found significant correlations among the features, we did not find any correlation between the numbers of 5'UTR features and protein tissue distribution and expression scores. We showed the existence of single nucleotide variants in relation to the 5'UTR features experimentally in a cohort of 105 breast cancer patients. 5'UTR features presumably prepare a complex playground, in which the other elements such as RNA binding proteins and non-coding RNAs play the major role in the fine-tuning of protein expression.

Role of Genetic Variation in ABC Transporters in Breast Cancer Prognosis and Therapy Response.

Breast cancer is the most common cancer in women in the world. The role of germline genetic variability in ATP-binding cassette (ABC) transporters in cancer chemoresistance and prognosis still needs to be elucidated. We used next-generation sequencing to assess associations of germline variants in coding and regulatory sequences of all human ABC genes with response of the patients to the neoadjuvant cytotoxic chemotherapy and disease-free survival (n = 105). A total of 43 prioritized variants associating with response or survival in the above testing phase were then analyzed by allelic discrimination in the large validation set (n = 802). Variants in ABCA4, ABCA9, ABCA12, ABCB5, ABCC5, ABCC8, ABCC11, and ABCD4 associated with response and variants in ABCA7, ABCA13, ABCC4, and ABCG8 with survival of the patients. No association passed a false discovery rate test, however, the rs17822931 (Gly180Arg) in ABCC11, associating with response, and the synonymous rs17548783 in ABCA13 (survival) have a strong support in the literature and are, thus, interesting for further research. Although replicated associations have not reached robust statistical significance, the role of ABC transporters in breast cancer should not be ruled out. Future research and careful validation of findings will be essential for assessment of genetic variation which was not in the focus of this study, e.g., non-coding sequences, copy numbers, and structural variations together with somatic mutations.

ABC Transporters and Their Role in the Neoadjuvant Treatment of Esophageal Cancer.

The prognosis of esophageal cancer (EC) is poor, despite considerable effort of both experimental scientists and clinicians. The tri-modality treatment consisting of neoadjuvant chemoradiation followed by surgery has remained the gold standard over decades, unfortunately, without significant progress in recent years. Suitable prognostic factors indicating which patients will benefit from this tri-modality treatment are missing. Some patients rapidly progress on the neoadjuvant chemoradiotherapy, which is thus useless and sometimes even harmful. At the same time, other patients achieve complete remission on neoadjuvant chemoradiotherapy and subsequent surgery may increase their risk of morbidity and mortality. The prognosis of patients ranges from excellent to extremely poor. Considering these differences, the role of drug metabolizing enzymes and transporters, among other factors, in the EC response to chemotherapy may be more important compared, for example, with pancreatic cancer where all patients progress on chemotherapy regardless of the treatment or disease stage. This review surveys published literature describing the potential role of ATP-binding cassette transporters, the genetic polymorphisms, epigenetic regulations, and phenotypic changes in the prognosis and therapy of EC. The review provides knowledge base for further research of potential predictive biomarkers that will allow the stratification of patients into defined groups for optimal therapeutic outcome.

Transcript expression and genetic variability analysis of caspases in breast carcinomas suggests CASP9 as the most interesting target.

BACKGROUND: Apoptosis plays a critical role in cancer cell survival and tumor development. We provide a hypothesis-generating screen for further research by exploring the expression profile and genetic variability of caspases (2, 3, 7, 8, 9, and 10) in breast carcinoma patients. This study addressed isoform-specific caspase transcript expression and genetic variability in regulatory sequences of caspases 2 and 9. METHODS: Gene expression profiling was performed by quantitative real-time PCR in tumor and paired non-malignant tissues of two independent groups of patients. Genetic variability was determined by high resolution melting, allelic discrimination, and sequencing analysis in tumor and peripheral blood lymphocyte DNA of the patients. RESULTS: CASP3 A+B and S isoforms were over-expressed in tumors of both patient groups. The CASP9 transcript was down-regulated in tumors of both groups of patients and significantly associated with expression of hormonal receptors and with the presence of rs4645978-rs2020903-rs4646034 haplotype in the CASP9 gene. Patients with a low intratumoral CASP9A/B isoform expression ratio (predicted to shift equilibrium towards anti-apoptotic isoform) subsequently treated with adjuvant chemotherapy had a significantly shorter disease-free survival than those with the high ratio (p=0.04). Inheritance of CC genotype of rs2020903 in CASP9 was associated with progesterone receptor expression in tumors (p=0.003). CONCLUSIONS: Genetic variability in CASP9 and expression of its splicing variants present targets for further study.

Role of family D ATP-binding cassette transporters (ABCD) in cancer.

ATP-binding cassette (ABC) transporters, belonging to the family D, are expressed in peroxisomes, endoplasmic reticulum or lysosomes. ABCD transporters play a role in transport of lipids, bile acids and vitamin B12 and associate with peroxisomal disorders. ABCD1 performs transport of coenzyme A esters of very-long-chain fatty acids (VLCFA) in peroxisomes and a number of mutations in ABCD1 gene were linked to an X-linked adrenoleucodystrophy (X-ALD). The role of ABCD transporters in tumour growth has not been studied in detail, but there is some evidence that ABCDs levels differ between undifferentiated stem or tumour cells and differentiated cells suggesting a possible link to tumorigenesis. In this mini-review, we discuss the available information about the role of ABCD transporters in cancer.

Importance of transcript levels of caspase-2 isoforms S and L for breast carcinoma progression.

AIM: A role of caspase-2 in chemotherapy-induced apoptosis has been suggested. Our study aimed to evaluate the prognostic and predictive importance of caspase-2 isoforms in breast cancer patients. MATERIALS & METHODS: Caspase-2L and -2S transcript levels were determined in paired tumor and non-malignant control tissues from 64 patients after neoadjuvant chemotherapy and 100 pretreatment patients (general set) by real-time PCR with absolute quantification. RESULTS: Low but statistically significant upregulation of caspase-2L in tumor versus control tissues was observed in both sets. Significant associations of the levels of caspase-2L, -2S or S/L ratio with clinical prognostic factors were observed. However, none of these associations were confirmed in both sets. Levels of caspase-2 isoforms or the S/L ratio did not significantly associate with progression-free survival in the general set or with chemotherapy response in the neoadjuvant set. CONCLUSION: Our results suggest that the role of caspase-2 isoforms in the progression of breast cancer may considerably differ between pre- and post-chemotherapy patients.

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.

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.