Landscape of somatic single nucleotide variants and indels in colorectal cancer and impact on survival

Colorectal cancer (CRC) is a biologically heterogeneous disease. To characterize its mutational profile, we conduct targeted sequencing of 205 genes for 2,105 CRC cases with survival data. Our data shows several findings in addition to enhancing the existing knowledge of CRC. We identify PRKCI, SPZ1, MUTYH, MAP2K4, FETUB, and TGFBR2 as additional genes significantly mutated in CRC. We find that among hypermutated tumors, an increased mutation burden is associated with improved CRC-specific survival (HR=0.42, 95% CI: 0.21-0.82). Mutations in TP53 are associated with poorer CRC-specific survival, which is most pronounced in cases carrying TP53 mutations with predicted 0% transcriptional activity (HR=1.53, 95% CI: 1.21-1.94). Furthermore, we observe differences in mutational frequency of several genes and pathways by tumor location, stage, and sex. Overall, this large study provides deep insights into somatic mutations in CRC, and their potential relationships with survival and tumor features. Large scale sequencing study is of paramount importance to unravel the heterogeneity of colorectal cancer. Here, the authors sequenced 205 cancer genes in more than 2000 tumours and identified additional mutated driver genes, determined that mutational burden and specific mutations in TP53 are associated with survival odds

A genome-wide association study of anorexia nervosa suggests a risk locus implicated in dysregulated leptin signaling

J. Kaprio, A. Palotie, A. Raevuori-Helkamaa ja S. Ripatti ovat työryhmän Eating Disorders Working Group of the Psychiatric Genomics Consortium jäseniä. Erratum in: Sci Rep. 2017 Aug 21;7(1):8379, doi: 10.1038/s41598-017-06409-3We conducted a genome-wide association study (GWAS) of anorexia nervosa (AN) using a stringently defined phenotype. Analysis of phenotypic variability led to the identification of a specific genetic risk factor that approached genome-wide significance (rs929626 in EBF1 (Early B-Cell Factor 1); P = 2.04 x 10(-7); OR = 0.7; 95% confidence interval (CI) = 0.61-0.8) with independent replication (P = 0.04), suggesting a variant-mediated dysregulation of leptin signaling may play a role in AN. Multiple SNPs in LD with the variant support the nominal association. This demonstrates that although the clinical and etiologic heterogeneity of AN is universally recognized, further careful sub-typing of cases may provide more precise genomic signals. In this study, through a refinement of the phenotype spectrum of AN, we present a replicable GWAS signal that is nominally associated with AN, highlighting a potentially important candidate locus for further investigation.Peer reviewe

Exploratory genome-wide interaction analysis of non-steroidal anti-inflammatory drugs and predicted gene expression on colorectal cancer risk

Background: Regular use of nonsteroidal anti-inflammatory drugs (NSAID) is associated with lower risk of colorectal cancer. Genome-wide interaction analysis on single variants (G × E) has identified several SNPs that may interact with NSAIDs to confer colorectal cancer risk, but variations in gene expression levels may also modify the effect of NSAID use. Therefore, we tested interactions between NSAID use and predicted gene expression levels in relation to colorectal cancer risk. Methods: Genetically predicted gene expressions were tested for interaction with NSAID use on colorectal cancer risk among 19,258 colorectal cancer cases and 18,597 controls from 21 observational studies. A Mixed Score Test for Interactions (MiSTi) approach was used to jointly assess G × E effects which are modeled via fixed interaction effects of the weighted burden within each gene set (burden) and residual G × E effects (variance). A false discovery rate (FDR) at 0.2 was applied to correct for multiple testing. Results: Among the 4,840 genes tested, genetically predicted expression levels of four genes modified the effect of any NSAID use on colorectal cancer risk, including DPP10 (PG×E = 1.96 × 10-4), KRT16 (PG×E = 2.3 × 10-4), CD14 (PG×E = 9.38 × 10-4), and CYP27A1 (PG×E = 1.44 × 10-3). There was a significant interaction between expression level of RP11-89N17 and regular use of aspirin only on colorectal cancer risk (PG×E = 3.23 × 10-5). No interactions were observed between predicted gene expression and nonaspirin NSAID use at FDR < 0.2. Conclusions: By incorporating functional information, we discovered several novel genes that interacted with NSAID use

Variation at 2q35 (PNKD and TMBIM1) influences colorectal cancer risk and identifies a pleiotropic effect with inflammatory bowel disease

To identify new risk loci for colorectal cancer (CRC), we conducted a meta-analysis of seven genome-wide association studies (GWAS) with independent replication, totalling 13 656 CRC cases and 21 667 controls of European ancestry. The combined analysis identified a new risk association for CRC at 2q35 marked by rs992157 (P = 3.15 x 10(-8), odds ratio = 1.10, 95% confidence interval = 1.06-1.13), which is intronic to PNKD (paroxysmal non-kinesigenic dyskinesia) and TMBIM1 (transmembrane BAX inhibitor motif containing 1). Intriguingly this susceptibility single-nucleotide polymorphism (SNP) is in strong linkage disequilibrium (r(2) = 0.90, D' = 0.96) with the previously discovered GWAS SNP rs2382817 for inflammatory bowel disease (IBD). Following on from this observation we examined for pleiotropy, or shared genetic susceptibility, between CRC and the 200 established IBD risk loci, identifying an additional 11 significant associations (false discovery rate [FDR]) <0.05). Our findings provide further insight into the biological basis of inherited genetic susceptibility to CRC, and identify risk factors that may influence the development of both CRC and IBD.Peer reviewe

CCAT2, a novel noncoding RNA mapping to 8q24, underlies metastatic progression and chromosomal instability in colon cancer

The functional roles of SNPs within the 8q24 gene desert in the cancer phenotype are not yet well understood. Here, we report that CCAT2, a novel long noncoding RNA transcript (lncRNA) encompassing the rs6983267 SNP, is highly overexpressed in microsatellite-stable colorectal cancer and promotes tumor growth, metastasis, and chromosomal instability. We demonstrate that MYC, miR-17-5p, and miR-20a are up-regulated by CCAT2 through TCF7L2-mediated transcriptional regulation. We further identify the physical interaction between CCAT2 and TCF7L2 resulting in an enhancement of WNT signaling activity. We show that CCAT2 is itself a WNT downstream target, which suggests the existence of a feedback loop. Finally, we demonstrate that the SNP status affects CCAT2 expression and the risk allele G produces more CCAT2 transcript. Our results support a new mechanism of MYC and WNT regulation by the novel lncRNA CCAT2 in colorectal cancer pathogenesis, and provide an alternative explanation of the SNP-conferred cancer risk

Mendelian randomisation analysis strongly implicates adiposity with risk of developing colorectal cancer

Background: Observational studies have associated adiposity with an increased risk of colorectal cancer (CRC). However, such studies do not establish a causal relationship. To minimise bias from confounding we performed a Mendelian randomisation (MR) analysis to examine the relationship between adiposity and CRC. Methods: We used SNPs associated with adult body mass index (BMI), waist-hip ratio (WHR), childhood obesity and birth weight as instrumental variables in a MR analysis of 9254 CRC cases and 18 386 controls. Results: In the MR analysis, the odds ratios (ORs) of CRC risk per unit increase in BMI, WHR and childhood obesity were 1.23 (95% CI: 1.02-1.49, P = 0.033), 1.59 (95% CI: 1.08-2.34, P = 0.019) and 1.07 (95% CI: 1.03-1.13, P = 0.018), respectively. There was no evidence for association between birth weight and CRC (OR = 1.22, 95% CI: 0.89-1.67, P = 0.22). Combining these data with a concurrent MR-based analysis for BMI and WHR with CRC risk (totalling to 18 190 cases, 27 617 controls) provided increased support, ORs for BMI and WHR were 1.26 (95% CI: 1.10-1.44, P = 7.7 x 10(-4)) and 1.40 (95% CI: 1.14-1.72, P = 1.2 x 10(-3)), respectively. Conclusions: These data provide further evidence for a strong causal relationship between adiposity and the risk of developing CRC highlighting the urgent need for prevention and treatment of adiposity.Peer reviewe

Retrospective evaluation of whole exome and genome mutation calls in 746 cancer samples

Funder: NCI U24CA211006Abstract: The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) curated consensus somatic mutation calls using whole exome sequencing (WES) and whole genome sequencing (WGS), respectively. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2,658 cancers across 38 tumour types, we compare WES and WGS side-by-side from 746 TCGA samples, finding that ~80% of mutations overlap in covered exonic regions. We estimate that low variant allele fraction (VAF < 15%) and clonal heterogeneity contribute up to 68% of private WGS mutations and 71% of private WES mutations. We observe that ~30% of private WGS mutations trace to mutations identified by a single variant caller in WES consensus efforts. WGS captures both ~50% more variation in exonic regions and un-observed mutations in loci with variable GC-content. Together, our analysis highlights technological divergences between two reproducible somatic variant detection efforts

Comparison of the clinical prediction model PREMM1,2,6 and molecular testing for the systematic identification of Lynch syndrome in colorectal cancer

Background Lynch syndrome is caused by germline mismatch repair (MMR) gene mutations. The PREMM1,2,6 model predicts the likelihood of a MMR gene mutation based on personal and family cancer history. Objective To compare strategies using PREMM1,2,6 and tumour testing (microsatellite instability (MSI) and/or immunohistochemistry (IHC) staining) to identify mutation carriers. Design Data from population-based or clinic-based patients with colorectal cancers enrolled through the Colon Cancer Family Registry were analysed. Evaluation included MSI, IHC and germline mutation analysis for MLH1, MSH2, MSH6 and PMS2. Personal and family cancer histories were used to calculate PREMM1,2,6 predictions. Discriminative ability to identify carriers from non-carriers using the area under the receiver operating characteristic curve (AUC) was assessed. Predictions were based on logistic regression models for (1) cancer assessment using PREMM1,2,6, (2) MSI, (3) IHC for loss of any MMR protein expression, (4) MSI+IHC, (5) PREMM1,2,6+MSI, (6) PREMM1,2,6+IHC, (7) PREMM1,2,6+IHC+MSI. Results Among 1651 subjects, 239 (14%) had mutations (90 MLH1, 125 MSH2, 24 MSH6). PREMM1,2,6 discriminated well with AUC 0.90 (95% CI 0.88 to 0.92). MSI alone, IHC alone, or MSI+IHC each had lower AUCs: 0.77, 0.82 and 0.82, respectively. The added value of IHC+PREMM1,2,6 was slightly greater than PREMM1,2,6+MSI (AUC 0.94 vs 0.93). Adding MSI to PREMM1,2,6+IHC did not improve discrimination. Conclusion PREMM1,2,6 and IHC showed excellent performance in distinguishing mutation carriers from noncarriers and performed best when combined. MSI may have a greater role in distinguishing Lynch syndrome from other familial colorectal cancer subtypes among cases with high PREMM1,2,6 scores where genetic evaluation does not disclose a MMR mutation

Long-term outcomes of laparoscopic liver resection for hepatocellular carcinoma: A propensity score matched analysis of a high-volume North American center

Background: Laparoscopic liver resections for malignancy are increasing worldwide, and yet data from North America are lacking. We aimed to assess the long-term outcomes of patients undergoing laparoscopic liver resection and open liver resection as a treatment for hepatocellular carcinoma. Methods: Patients undergoing liver resection for hepatocellular carcinoma between January 2008 and December 2019 were retrospectively studied. A propensity score matching was performed using patient demographics, laboratory parameters, etiology of liver disease, liver function, and tumor characteristics. Primary outcomes included overall survival and cumulative incidence of recurrence. Kaplan-Meier and competing risk cumulative incidence were used for survival analyses. Multivariable Cox regression and Fine-Gray proportional hazard regression were performed to determine hazard for death and recurrence, respectively. Results: Three hundred and ninety-one patients were identified (laparoscopic liver resection: 110; open liver resection: 281). After propensity score matching, 149 patients remained (laparoscopic liver resection: 57; open liver resection: 92). There were no significant differences between groups with regard to extent of hepatectomy performed and tumor characteristics. The laparoscopic liver resection group experienced a lower proportion of ≥Clavien-Dindo grade III complications (14% vs 29%; P = .01). In the matched cohort, the 1-, 3-, and 5-year overall survival rate in the laparoscopic liver resection versus open liver resection group was 90.9%, 79.3%, 70.5% vs 91.3%, 88.5%, 83.1% (P = .26), and the cumulative incidence of recurrence 31.1%, 59.7%, 62.9% vs 18.9%, 40.6%, 49.2% (P = .06), respectively. Conclusion: This study represents the largest single institutional study from North America comparing long-term oncologic outcomes of laparoscopic liver resection and open liver resection as a treatment for primary hepatocellular carcinoma. The combination of reduced short-term complications and equivalent long-term oncologic outcomes favor the laparoscopic approach when feasible

Genome-wide search for determinants of survival in 1,926 patients with advanced colorectal cancer with follow-up in over 22,000 patients

Background While genome-wide association studies (GWAS) have identified germline variants influencing the risk of developing colorectal cancer (CRC), there has been limited examination of the possible role of inherited variation as a determinant of patient outcome. Patients and methods We performed a GWAS for overall survival (OS) in 1926 patients with advanced CRC from the COIN and COIN-B clinical trials. For single nucleotide polymorphisms (SNPs) showing an association with OS (P < 1.0 × 10−5), we conducted sensitivity analyses based on the time from diagnosis to death and sought independent replications in 5675 patients from the Study of Colorectal Cancer in Scotland (SOCCS) and 16,964 patients from the International Survival Analysis in Colorectal cancer Consortium (ISACC). We analysed the Human Protein Atlas to determine if ERBB4 expression was associated with survival in 438 patients with colon adenocarcinomas. Results The most significant SNP associated with OS was rs79612564 in ERBB4 (hazard ratio [HR] = 1.24, 95% confidence interval [CI] = 1.16–1.32, P = 1.9 × 10−7). SNPs at 17 loci had suggestive associations for OS and all had similar effects on the time from diagnosis to death. No lead SNPs were independently replicated in the meta-analysis of all patients from SOCCS and ISACC. However, rs79612564 was significant in stage-IV patients from SOCCS (P = 2.1 × 10−2) but not ISACC (P = 0.89) and SOCCS combined with COIN and COIN-B attained genome-wide significance (P = 1.7 × 10−8). Patients with high ERBB4 expression in their colon adenocarcinomas had worse survival (HR = 1.50, 95% CI = 1.1–1.9, P = 4.6 × 10−2). Conclusions Genetic and expression data support a potential role for rs79612564 in the receptor tyrosine kinase ERBB4 as a predictive biomarker of survival