Evolutionary mode and timing of dissemination of high-grade serous carcinomas.

Dissemination within the peritoneal cavity is a main determinant of poor patient outcomes from high-grade serous carcinomas (HGSCs). The dissemination process is poorly understood from a cancer evolutionary perspective. We reconstructed the evolutionary trajectories across a median of 5 tumor sites and regions from each of 23 patients based on deep whole-exome sequencing. Polyclonal cancer origin was detected in 1 patient. Ovarian tumors had more complex subclonal architectures than other intraperitoneal tumors in each patient, which indicated that tumors developed earlier in the ovaries. Three common modes of dissemination were identified, including monoclonal or polyclonal dissemination of monophyletic (linear) or polyphyletic (branched) subclones. Mutation profiles of initial or disseminated clones varied greatly among cancers, but recurrent mutations were found in 7 cancer-critical genes, including TP53, BRCA1, BRCA2, and DNMT3A, and in the PI3K/AKT1 pathway. Disseminated clones developed late in the evolutionary trajectory models of most cancers, in particular in cancers with DNA damage repair deficiency. Polyclonal dissemination was predicted to occur predominantly as a single and rapid wave, but chemotherapy exposure was associated with higher genomic diversity of disseminated clones. In conclusion, we described three common evolutionary dissemination modes across HGSCs and proposed factors associated with dissemination diversity.

Deviating Alternative Splicing as a Molecular Subtype of Microsatellite Stable Colorectal Cancer.

PURPOSE: Colorectal cancer (CRC) is the third most commonly diagnosed cancer worldwide. Biomarkers to aid in prognostication and treatment decisions are in high demand, and to facilitate their development, a better understanding of the underlying biology of the highly heterogeneous disease is needed. METHODS: A genome-scale alternative splicing (AS) analysis using RNA-sequencing data from primary microsatellite stable (MSS) CRCs from 127 patients was performed. Splice variant-specific expression levels of individual cancer samples were compared with the total set of samples, and a metric for a tumor sample's global amount of deviating AS was developed. This metric varied considerably across the cohort and ranged from 6 to 282 deviating AS events per tumor sample. A threshold of 45 or more deviating events was set to distinguish cancers with high (n = 44) and low (n = 83) levels of deviating AS. RESULTS: Patients with high amounts of AS deviations had significantly shorter time to relapse compared with patients with fewer deviations (P = .04). Furthermore, differential gene expression analysis revealed nine known cancer-critical genes that are significantly upregulated in samples with high amounts of deviating AS. Validation of the results in an independent cohort of MSS CRCs showed the same tendency toward shorter progression-free survival among the high-deviation group. In both cohorts, enrichment for growth factors was identified among upregulated genes associated with this phenotype. CONCLUSION: There is a large variation in the amount of deviating AS among MSS CRCs, and we provide evidence that those with high amounts of deviations represent different cancer biology.

Somatic mutations reveal complex metastatic seeding from multifocal primary prostate cancer.

Primary prostate cancer shows a striking intraorgan molecular heterogeneity, with multiple spatially separated malignant foci in the majority of patients. Metastatic prostate cancer, however, typically reveals more homogenous molecular profiles, suggesting a monoclonal origin of the metastatic lesions. Longitudinal mutational spectra, comparing multiple primary lesions with metastases from the same patients remain poorly defined. We have here analyzed somatic mutations in multisampled, spatio-temporal biobanked lesions (38 samples from primary foci and 1 sample from each of 8 metastases from seven prostate cancer patients) applying a custom-designed panel targeting 68 prostate cancer relevant genes. The metastatic samples were taken at time of primary surgery and up to 7 years later, and sampling included circulating tumor DNA in plasma or solid metastatic tissue samples. A total of 282 somatic mutations were detected, with a range of 0 to 25 mutations per sample. Although seven samples had solely private mutations, the remaining 39 samples had both private and shared mutations. Seventy-four percent of mutations in metastases were not found in any primary samples, and vice versa, 96% of mutations in primary cancers were not found in any metastatic samples. However, for three patients, shared mutations were found suggesting the focus of origin, including mutations in AKT1, FOXA1, HOXB13, RB1 and TP53. In conclusion, the spatio-temporal heterogeneous nature of multifocal disease is emphasized in our study, and underlines the importance of testing a recent sample in genomics-based precision medicine for metastatic prostate cancer.

NRF2 drives an oxidative stress response predictive of breast cancer.

Many breast cancer patients are diagnosed with small, well-differentiated, hormone receptor-positive tumors. Risk of relapse is not easily identified in these patients, resulting in overtreatment. To identify metastasis-related gene expression patterns, we compared the transcriptomes of the non-metastatic 67NR and metastatic 66cl4 cell lines from the murine 4T1 mammary tumor model. The transcription factor nuclear factor, erythroid 2-like 2 (NRF2, encoded by NFE2L2) was constitutively activated in the metastatic cells and tumors, and correspondingly a subset of established NRF2-regulated genes was also upregulated. Depletion of NRF2 increased basal levels of reactive oxygen species (ROS) and severely reduced ability to form primary tumors and lung metastases. Consistently, a set of NRF2-controlled genes was elevated in breast cancer biopsies. Sixteen of these were combined into a gene expression signature that significantly improves the PAM50 ROR score, and is an independent, strong predictor of prognosis, even in hormone receptor-positive tumors.

Expressed prognostic biomarkers for primary prostate cancer independent of multifocality and transcriptome heterogeneity.

The majority of prostate cancer patients are diagnosed with multiple primary malignant foci. The distinct foci are exceptionally heterogeneous with regard to DNA mutations, but whether this is recapitulated at the transcriptome level remains unknown. In this study, inter- and intrafocal heterogeneity has been assessed by whole-transcriptome sequencing of 87 tissue samples from 23 patients with localized prostate cancer treated with radical prostatectomy. From each patient, multiple samples were taken from one or more malignant foci, in addition to one sample from benign prostate tissue. Transcriptomic profiles of different malignant foci from the same patient showed a similar level of heterogeneity as tumors from different patients. This applies to expression of genes, fusion genes, and somatic mutations. Within-patient pair-wise analyses identified expression patterns linked to ETS status and extraprostatic extension. A set of 62 genes were found with low intrapatient heterogeneity and high interpatient heterogeneity, retaining stable expression profiles across foci within the same patient. Among these, 16 genes are associated with biochemical recurrence in a separately published study and are therefore nominated as biomarkers with prognostic value regardless of which malignant focus is sampled. In conclusion, an extensive heterogeneity in multifocal prostate cancer is confirmed at the gene expression level. Diagnostic biomarkers were identified for ETS positive samples and samples from extraprostatic extensions. Finally, prognostic biomarkers independent of multifocal heterogeneity were found.

Collision tumors revealed by prospectively assessing subtype-defining molecular alterations in 904 individual prostate cancer foci.

BACKGROUNDProstate cancer is multifocal with distinct molecular subtypes. The utility of genomic subtyping has been challenged due to inter- and intrafocal heterogeneity. We sought to characterize the subtype-defining molecular alterations of primary prostate cancer across all tumor foci within radical prostatectomy (RP) specimens and determine the prevalence of collision tumors.METHODSFrom the Early Detection Research Network cohort, we identified 333 prospectively collected RPs from 2010 to 2014 and assessed ETS-related gene (ERG), serine peptidase inhibitor Kazal type 1 (SPINK1), phosphatase and tensin homolog (PTEN), and speckle type BTB/POZ protein (SPOP) molecular status. We utilized dual ERG/SPINK1 immunohistochemistry and fluorescence in situ hybridization to confirm ERG rearrangements and characterize PTEN deletion, as well as high-resolution melting curve analysis and Sanger sequencing to determine SPOP mutation status.RESULTSBased on index focus alone, ERG, SPINK1, PTEN, and SPOP alterations were identified in 47.5%, 10.8%, 14.3%, and 5.1% of RP specimens, respectively. In 233 multifocal RPs with ERG/SPINK1 status in all foci, 139 (59.7%) had discordant molecular alterations between foci. Collision tumors, as defined by discrepant ERG/SPINK1 status within a single focus, were identified in 29 (9.4%) RP specimens.CONCLUSIONInterfocal molecular heterogeneity was identified in about 60% of multifocal RP specimens, and collision tumors were present in about 10%. We present this phenomenon as a model for the intrafocal heterogeneity observed in previous studies and propose that future genomic studies screen for collision tumors to better characterize molecular heterogeneity.FUNDINGEarly Detection Research Network US National Cancer Institute (NCI) 5U01 CA111275-09, Center for Translational Pathology at Weill Cornell Medicine (WCM) Department of Pathology and Laboratory Medicine, US NCI (WCM SPORE in Prostate Cancer, P50CA211024-01), R37CA215040, Damon Runyon Cancer Research Foundation, US MetLife Foundation Family Clinical Investigator Award, Norwegian Cancer Society (grant 208197), and South-Eastern Norway Regional Health Authority (grant 2019016 and 2020063).

The expressed mutational landscape of microsatellite stable colorectal cancers.

BACKGROUND: Colorectal cancer is the 2nd leading cause of cancer-related deaths with few patients benefiting from biomarker-guided therapy. Mutation expression is essential for accurate interpretation of mutations as biomarkers, but surprisingly, little has been done to analyze somatic cancer mutations on the expression level. We report a large-scale analysis of allele-specific mutation expression. METHODS: Whole-exome and total RNA sequencing was performed on 137 samples from 121 microsatellite stable colorectal cancers, including multiregional samples of primary and metastatic tumors from 4 patients. Data were integrated with allele-specific resolution. Results were validated in an independent set of 241 colon cancers. Therapeutic associations were explored by pharmacogenomic profiling of 15 cell lines or patient-derived organoids. RESULTS: The median proportion of expressed mutations per tumor was 34%. Cancer-critical mutations had the highest expression frequency (gene-wise mean of 58%), independent of frequent allelic imbalance. Systematic deviation from the general pattern of expression levels according to allelic frequencies was detected, including preferential expression of mutated alleles dependent on the mutation type and target gene. Translational relevance was suggested by correlations of KRAS/NRAS or TP53 mutation expression levels with downstream oncogenic signatures (p < 0.03), overall survival among patients with stage II and III cancer (KRAS/NRAS: hazard ratio 6.1, p = 0.0070), and targeted drug sensitivity. The latter was demonstrated for EGFR and MDM2 inhibition in pre-clinical models. CONCLUSIONS: Only a subset of mutations in microsatellite stable colorectal cancers were expressed, and the "expressed mutation dose" may provide an opportunity for more fine-tuned biomarker interpretations.

Molecular correlates of sensitivity to PARP inhibition beyond homologous recombination deficiency in pre-clinical models of colorectal cancer point to wild-type TP53 activity.

BACKGROUND: PARP inhibitors are active in various tumour types beyond BRCA-mutant cancers, but their activity and molecular correlates in colorectal cancer (CRC) are not well studied. METHODS: Mutations and genome-wide mutational patterns associated with homologous recombination deficiency (HRD) were investigated in 255 primary CRCs with whole-exome sequencing and/or DNA copy number data. Efficacy of five PARP inhibitors and their molecular correlates were evaluated in 93 CRC cell lines partly annotated with mutational-, DNA copy number-, and/or gene expression profiles. Post-treatment gene expression profiling and specific protein expression analyses were performed in two pairs of PARP inhibitor sensitive and resistant cell lines. FINDINGS: A subset of microsatellite stable (MSS) CRCs had truncating mutations in homologous recombination-related genes, but these were not associated with genomic signatures of HRD. Eight CRC cell lines (9%) were sensitive to PARP inhibition, but sensitivity was not predicted by HRD-related genomic and transcriptomic signatures. In contrast, drug sensitivity in MSS cell lines was strongly associated with TP53 wild-type status (odds ratio 15.7, p = 0.023) and TP53-related expression signatures. Increased downstream TP53 activity was among the primary response mechanisms, and TP53 inhibition antagonized the effect of PARP inhibitors. Wild-type TP53-mediated suppression of RAD51 was identified as a possible mechanism of action for sensitivity to PARP inhibition. INTERPRETATION: PARP inhibitors are active in a subset of CRC cell lines and preserved TP53 function may increase the likelihood of response.

Frequent copy number gains of SLC2A3 and ETV1 in testicular embryonal carcinomas.

Testicular germ cell tumours (TGCTs) appear as different histological subtypes or mixtures of these. They show similar, multiple DNA copy number changes, where gain of 12p is pathognomonic. However, few high-resolution analyses have been performed and focal DNA copy number changes with corresponding candidate target genes remain poorly described for individual subtypes. We present the first high-resolution DNA copy number aberration (CNA) analysis on the subtype embryonal carcinomas (ECs), including 13 primary ECs and 5 EC cell lines. We identified recurrent gains and losses and allele-specific CNAs. Within these regions, we nominate 30 genes that may be of interest to the EC subtype. By in silico analysis of data from 150 TGCTs from The Cancer Genome Atlas (TCGA), we further investigated CNAs, RNA expression, somatic mutations and fusion transcripts of these genes. Among primary ECs, ploidy ranged between 2.3 and 5.0, and the most common aberrations were DNA copy number gains at chromosome (arm) 7, 8, 12p, and 17, losses at 4, 10, 11, and 18, replicating known TGCT genome characteristics. Gain of whole or parts of 12p was found in all samples, including a highly amplified 100 kbp segment at 12p13.31, containing SLC2A3. Gain at 7p21, encompassing ETV1, was the second most frequent aberration. In conclusion, we present novel CNAs and the genes located within these regions, where the copy number gain of SLC2A3 and ETV1 are of interest, and which copy number levels also correlate with expression in TGCTs.

High Concordance and Negative Prognostic Impact of RAS/BRAF/PIK3CA Mutations in Multiple Resected Colorectal Liver Metastases.

BACKGROUND: The prevalence and clinical implications of genetic heterogeneity in patients with multiple colorectal liver metastases remain largely unknown. In a prospective series of patients undergoing resection of colorectal liver metastases, the aim was to investigate the inter-metastatic and primary-to-metastatic heterogeneity of mutations in KRAS, NRAS, BRAF, and PIK3CA and their prognostic impact. PATIENTS AND METHODS: We analyzed the mutation status among 372 liver metastases and 78 primary tumors from 106 patients by methods used in clinical routine testing, by Sanger sequencing, by next-generation sequencing (NGS), and/or by droplet digital polymerase chain reaction. The 3-year cancer-specific survival (CSS) was analyzed using the Kaplan-Meier method. RESULTS: Although Sanger sequencing indicated inter-metastatic mutation heterogeneity in 14 of 97 patients (14%), almost all cases were refuted by high-sensitive NGS. Also, heterogeneity among metastatic deposits was concluded only for PIK3CA in 2 patients. Similarly, primary-to-metastatic heterogeneity was indicated in 8 of 78 patients (10%) using Sanger sequencing but for only 2 patients after NGS, showing the emergence of 1 KRAS and 1 PIK3CA mutation in the metastatic lesions. KRAS mutations were present in 53 of 106 patients (50%) and were associated with poorer 3-year CSS after liver resection (37% vs. 61% for KRAS wild-type; P = .004). Poor prognostic associations were found also for the combination of KRAS/NRAS/BRAF mutations compared with triple wild-type (P = .002). CONCLUSION: Intra-patient mutation heterogeneity was virtually undetected, both between the primary tumor and the liver metastases and among the metastatic deposits. KRAS mutations separately, and KRAS/NRAS/BRAF mutations combined, were associated with poor patient survival after partial liver resection.

Technical differences between sequencing and microarray platforms impact transcriptomic subtyping of colorectal cancer.

Gene expression profiling has increasing relevance in the molecular screening of patients with colorectal cancer (CRC). We investigated potential platform-specific effects on transcriptomic subtyping according to established frameworks by comparisons of expression profiles from RNA sequencing and exon-resolution microarrays in 126 primary microsatellite stable CRCs. There was a strong platform correspondence in global gene expression levels, albeit with systematic technical bias likely attributed to few sequencing reads covering short (<2000 nucleotides) and/or lowly expressed genes (<1 FPKM), as well as over-saturation of highly expressed genes on microarrays. Classification concordances according to both the consensus molecular subtypes and CRC intrinsic subtypes (CRIS) were also strong, but with disproportionate subtype distributions between platforms caused by frequent disagreements in adherence to sample classification thresholds. Subtypes defined largely by genes expressed at low levels, including the CRIS-D subtype and the estimated level of tumor-infiltrating cytotoxic lymphocytes, had a weaker correspondence in classification metrics between platforms. In conclusion, even subtle differences between platforms suggest that clinical translation of transcriptomic CRC subtyping frameworks is dependent on assay standardization, and systematic technical biases reinforce the need for careful selection of classifier genes.

GREM1 is associated with metastasis and predicts poor prognosis in ER-negative breast cancer patients.

BACKGROUND: In breast cancer, activation of bone morphogenetic protein (BMP) signaling and elevated levels of BMP-antagonists have been linked to tumor progression and metastasis. However, the simultaneous upregulation of BMPs and their antagonist, and the fact that both promote tumor aggressiveness seems contradictory and is not fully understood. METHODS: We analyzed the transcriptomes of the metastatic 66cl4 and the non-metastatic 67NR cell lines of the 4T1 mouse mammary tumor model to search for factors that promote metastasis. CRISPR/Cas9 gene editing was used for mechanistic studies in the same cell lines. Furthermore, we analyzed gene expression patterns in human breast cancer biopsies obtained from public datasets to evaluate co-expression and possible relations to clinical outcome. RESULTS: We found that mRNA levels of the BMP-antagonist Grem1, encoding gremlin1, and the ligand Bmp4 were both significantly upregulated in cells and primary tumors of 66cl4 compared to 67NR. Depletion of gremlin1 in 66cl4 could impair metastasis to the lungs in this model. Furthermore, we found that expression of Grem1 correlated with upregulation of several stem cell markers in 66cl4 cells compared to 67NR cells. Both in the mouse model and in patients, expression of GREM1 associated with extracellular matrix organization, and formation, biosynthesis and modification of collagen. Importantly, high expression of GREM1 predicted poor prognosis in estrogen receptor negative breast cancer patients. Analyses of large patient cohorts revealed that amplification of genes encoding BMP-antagonists and elevation of the corresponding transcripts is evident in biopsies from more than half of the patients and much more frequent for the secreted BMP-antagonists than the intracellular inhibitors of SMAD signaling. CONCLUSION: In conclusion, our results show that GREM1 is associated with metastasis and predicts poor prognosis in ER-negative breast cancer patients. Gremlin1 could represent a novel target for therapy.

Transcriptional and functional consequences of TP53 splice mutations in colorectal cancer.

TP53 mutations are common in colorectal cancer (CRC). Most TP53 sequencing studies have been restricted to coding regions, but recent studies have revealed that splice mutations can generate transcript variants with distinct tumorigenic and prognostic properties. Here, we performed unrestricted sequencing of all coding sequences and splice regions of TP53 in a single-hospital series of 401 primary CRCs. TP53 splice mutations were detected in 4% of the cases (N = 16), considerably more frequent than reported in major databases, and they were mutually exclusive to exon mutations. RNA sequencing revealed high-level expression of aberrant transcript variants in the majority of splice mutated tumors (75%). Most variants were predicted to produce truncated TP53 proteins, including one sample expressing the potentially oncogenic and druggable p53ψ isoform. Despite heterogeneous transcript structures, downstream transcriptional profiling revealed that TP53 splice mutations had similar effects on TP53 target gene expression and pathway activity as exonic mutations. Intriguingly, TP53 splice mutations were associated with worse 5-year relapse-free survival in stage II disease, compared to both TP53 wild-type and exon mutations (P = 0.007). These data highlight the importance of including splice regions when examining the biological and clinical consequences of TP53 mutations in CRC.

Multifocal Primary Prostate Cancer Exhibits High Degree of Genomic Heterogeneity.

BACKGROUND: Most primary prostate cancers are multifocal with individual tumors harboring different aggressiveness; however, the genomic heterogeneity among these tumors is poorly understood. OBJECTIVE: To better understand the biological basis for clinical variability among different lesions, we sought to comprehensively characterize the heterogeneity of somatic gene mutations in multifocal prostate cancer. DESIGN, SETTING, AND PARTICIPANTS: High-coverage whole-exome sequencing of 153 frozen tissue samples, taken from two to three distinct tumor foci and one non-cancerous area from each of 41 patients, covering a total of 89 tumor foci. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: State-of-the-art bioinformatics tools for mutation calling and copy number determination from whole-exome sequencing data. RESULTS AND LIMITATIONS: We found a very high degree of interfocal heterogeneity among tumors, that is, 76% of pairwise-compared tumor foci from the same prostatectomy specimen had no point mutations in common and DNA copy number changes were rarely shared across cancer foci. The few point mutations shared across tumor foci were seldom in cancer-critical genes. CONCLUSIONS: In this first large genomic heterogeneity study of primary prostate cancer, we observe that different tumor foci within the same patient are genetically distinct, only rarely sharing any somatic gene mutations, including those in cancer driver genes. This heterogeneity affects how genomics-based management of prostate cancer can be implemented, as information from all tumor foci is necessary to draw valid conclusions about the cancer's genomic alterations. PATIENT SUMMARY: Most primary prostate cancers consist of multiple tumors within the same organ, but little is known about their relationships. We have compared the sets of gene mutations among such tumors and found that they only exceptionally have any in common. This will influence treatment decisions in the future as each tumor's mutations will render it unique and have to be considered to gain the best treatment results.

Involvement of DPP9 in gene fusions in serous ovarian carcinoma.

BACKGROUND: A fusion gene is a hybrid gene consisting of parts from two previously independent genes. Chromosomal rearrangements leading to gene breakage are frequent in high-grade serous ovarian carcinomas and have been reported as a common mechanism for inactivating tumor suppressor genes. However, no fusion genes have been repeatedly reported to be recurrent driver events in ovarian carcinogenesis. We combined genomic and transcriptomic information to identify novel fusion gene candidates and aberrantly expressed genes in ovarian carcinomas. METHODS: Examined were 19 previously karyotyped ovarian carcinomas (18 of the serous histotype and one undifferentiated). First, karyotypic aberrations were compared to fusion gene candidates identified by RNA sequencing (RNA-seq). In addition, we used exon-level gene expression microarrays as a screening tool to identify aberrantly expressed genes possibly involved in gene fusion events, and compared the findings to the RNA-seq data. RESULTS: We found a DPP9-PPP6R3 fusion transcript in one tumor showing a matching genomic 11;19-translocation. Another tumor had a rearrangement of DPP9 with PLIN3. Both rearrangements were associated with diminished expression of the 3' end of DPP9 corresponding to the breakpoints identified by RNA-seq. For the exon-level expression analysis, candidate fusion partner genes were ranked according to deviating expression compared to the median of the sample set. The results were collated with data obtained from the RNA-seq analysis. Several fusion candidates were identified, among them TMEM123-MMP27, ZBTB46-WFDC13, and PLXNB1-PRKAR2A, all of which led to stronger expression of the 3' genes. In view of our previous findings of nonrandom rearrangements of chromosome 19 in this cancer type, particular emphasis was given to changes of this chromosome and a DDA1-FAM129C fusion event was identified. CONCLUSIONS: We have identified novel fusion gene candidates in high-grade serous ovarian carcinoma. DPP9 was involved in two different fusion transcripts that both resulted in deregulated expression of the 3' end of the transcript and thus possible loss of the active domains in the DPP9 protein. The identified rearrangements might play a role in tumorigenesis or tumor progression.

Multi-omics of 34 colorectal cancer cell lines - a resource for biomedical studies.

BACKGROUND: Colorectal cancer (CRC) cell lines are widely used pre-clinical model systems. Comprehensive insights into their molecular characteristics may improve model selection for biomedical studies. METHODS: We have performed DNA, RNA and protein profiling of 34 cell lines, including (i) targeted deep sequencing (n = 612 genes) to detect single nucleotide variants and insertions/deletions; (ii) high resolution DNA copy number profiling; (iii) gene expression profiling at exon resolution; (iv) small RNA expression profiling by deep sequencing; and (v) protein expression analysis (n = 297 proteins) by reverse phase protein microarrays. RESULTS: The cell lines were stratified according to the key molecular subtypes of CRC and data were integrated at two or more levels by computational analyses. We confirm that the frequencies and patterns of DNA aberrations are associated with genomic instability phenotypes and that the cell lines recapitulate the genomic profiles of primary carcinomas. Intrinsic expression subgroups are distinct from genomic subtypes, but consistent at the gene-, microRNA- and protein-level and dominated by two distinct clusters; colon-like cell lines characterized by expression of gastro-intestinal differentiation markers and undifferentiated cell lines showing upregulation of epithelial-mesenchymal transition and TGFß signatures. This sample split was concordant with the gene expression-based consensus molecular subtypes of primary tumors. Approximately » of the genes had consistent regulation at the DNA copy number and gene expression level, while expression of gene-protein pairs in general was strongly correlated. Consistent high-level DNA copy number amplification and outlier gene- and protein- expression was found for several oncogenes in individual cell lines, including MYC and ERBB2. CONCLUSIONS: This study expands the view of CRC cell lines as accurate molecular models of primary carcinomas, and we present integrated multi-level molecular data of 34 widely used cell lines in easily accessible formats, providing a resource for preclinical studies in CRC.

chimeraviz: a tool for visualizing chimeric RNA.

SUMMARY: Advances in high-throughput RNA sequencing have enabled more efficient detection of fusion transcripts, but the technology and associated software used for fusion detection from sequencing data often yield a high false discovery rate. Good prioritization of the results is important, and this can be helped by a visualization framework that automatically integrates RNA data with known genomic features. Here we present chimeraviz , a Bioconductor package that automates the creation of chimeric RNA visualizations. The package supports input from nine different fusion-finder tools: deFuse, EricScript, InFusion, JAFFA, FusionCatcher, FusionMap, PRADA, SOAPfuse and STAR-FUSION. AVAILABILITY AND IMPLEMENTATION: chimeraviz is an R package available via Bioconductor ( https://bioconductor.org/packages/release/bioc/html/chimeraviz.html ) under Artistic-2.0. Source code and support is available at GitHub ( https://github.com/stianlagstad/chimeraviz ). CONTACT: rolf.i.skotheim@rr-research.no. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Multilevel genomics of colorectal cancers with microsatellite instability-clinical impact of JAK1 mutations and consensus molecular subtype 1.

BACKGROUND: Approximately 15% of primary colorectal cancers have DNA mismatch repair deficiency, causing a complex genome with thousands of small mutations-the microsatellite instability (MSI) phenotype. We investigated molecular heterogeneity and tumor immunogenicity in relation to clinical endpoints within this distinct subtype of colorectal cancers. METHODS: A total of 333 primary MSI+ colorectal tumors from multiple cohorts were analyzed by multilevel genomics and computational modeling-including mutation profiling, clonality modeling, and neoantigen prediction in a subset of the tumors, as well as gene expression profiling for consensus molecular subtypes (CMS) and immune cell infiltration. RESULTS: Novel, frequent frameshift mutations in four cancer-critical genes were identified by deep exome sequencing, including in CRTC1, BCL9, JAK1, and PTCH1. JAK1 loss-of-function mutations were validated with an overall frequency of 20% in Norwegian and British patients, and mutated tumors had up-regulation of transcriptional signatures associated with resistance to anti-PD-1 treatment. Clonality analyses revealed a high level of intra-tumor heterogeneity; however, this was not associated with disease progression. Among the MSI+ tumors, the total mutation load correlated with the number of predicted neoantigens (P = 4 × 10-5), but not with immune cell infiltration-this was dependent on the CMS class; MSI+ tumors in CMS1 were highly immunogenic compared to MSI+ tumors in CMS2-4. Both JAK1 mutations and CMS1 were favorable prognostic factors (hazard ratios 0.2 [0.05-0.9] and 0.4 [0.2-0.9], respectively, P = 0.03 and 0.02). CONCLUSIONS: Multilevel genomic analyses of MSI+ colorectal cancer revealed molecular heterogeneity with clinical relevance, including tumor immunogenicity and a favorable patient outcome associated with JAK1 mutations and the transcriptomic subgroup CMS1, emphasizing the potential for prognostic stratification of this clinically important subtype. See related research highlight by Samstein and Chan 10.1186/s13073-017-0438-9.

Novel RNA variants in colorectal cancers.

With an annual estimated incidence of 1.4 million, and a five-year survival rate of 60%, colorectal cancer (CRC) is a major clinical burden. To identify novel RNA variants in CRC, we analyzed exon-level microarray expression data from a cohort of 202 CRCs. We nominated 25 genes with increased expression of their 3' parts in at least one cancer sample each. To efficiently investigate underlying transcript structures, we developed an approach using rapid amplification of cDNA ends followed by high throughput sequencing (RACE-seq). RACE products from the targeted genes in 23 CRC samples were pooled together and sequenced. We identified VWA2-TCF7L2, DHX35-BPIFA2 and CASZ1-MASP2 as private fusion events, and novel transcript structures for 17 of the 23 other candidate genes. The high-throughput approach facilitated identification of CRC specific RNA variants. These include a recurrent read-through fusion transcript between KLK8 and KLK7, and a splice variant of S100A2. Both of these were overrepresented in CRC tissue and cell lines from external RNA-seq datasets.