Patterns of somatic alterations between matched primary and metastatic colorectal tumors characterized by whole-genome sequencing

Whole exome sequencing reveals diverse genomic relatedness between paired concurrent endometrial and ovarian carcinomas

INTRODUCTION

Concurrent non-serous endometrial and ovarian tumours are often managed clinically as two separate primary tumours, but almost all exhibit evidence of a genomic relationship.

METHODOLOGY

This study investigates the extent of relatedness using whole exome sequencing, which was performed on paired non-serous endometrial and ovarian carcinomas from 27 patients with concurrent tumours in a cohort with detailed clinicopathological annotation. Four whole exome sequencing-derived parameters (mutation, mutational burden, mutational signatures and mutant allele tumour heterogeneity scores) were used to develop a novel unsupervised model for the assessment of genomic similarity to infer genomic relatedness of paired tumours.

RESULTS

This novel model demonstrated genomic relatedness across all four parameters in all tumour pairs. Mutations in PTEN, ARID1A, CTNNB1, KMT2D and PIK3CA occurred most frequently and 24 of 27 (89 %) tumour pairs shared identical mutations in at least one of these genes, with all pairs sharing mutations in a number of other genes. Ovarian endometriosis, CTNNB1 exon 3 mutation, and progression and death from disease were present across the similarity ranking. Mismatch repair deficiency was associated with less genomically similar pairs.

DISCUSSION

Although there was diversity across the cohort, the presence of genomic similarity in all paired tumours supports the hypothesis that concurrent non-serous endometrial and ovarian carcinomas are genomically related and may have arisen from a common origin.

Distinct molecular profiles drive multifaceted characteristics of colorectal cancer metastatic seeds

Metastasis of primary tumors remains a challenge for early diagnosis and prevention. The cellular properties and molecular drivers of metastatically competent clones within primary tumors remain unclear. Here, we generated 10-16 single cell-derived lines from each of three colorectal cancer (CRC) tumors to identify and characterize metastatic seeds. We found that intrinsic factors conferred clones with distinct metastatic potential and cellular communication capabilities, determining organ-specific metastasis. Poorly differentiated or highly metastatic clones, rather than drug-resistant clones, exhibited poor clinical prognostic impact. Personalized genetic alterations, instead of mutation burden, determined the occurrence of metastatic potential during clonal evolution. Additionally, we developed a gene signature for capturing metastatic potential of primary CRC tumors and demonstrated a strategy for identifying metastatic drivers using isogenic clones with distinct metastatic potential in primary tumors. This study provides insight into the origin and mechanisms of metastasis and will help develop potential anti-metastatic therapeutic targets for CRC patients.

Effect of aza-BODIPY-photodynamic therapy on the expression of carcinoma-associated genes and cell death mode

BACKGROUND

Breast cancer is the most common cancer affecting women worldwide.Photodynamic therapy(PDT) has now proven to be a promising form of cancer therapy due to its targeted and low cytotoxicity to healthy cells and tissues.PDT is a technique used to create cell death localized by light after application of a light-sensitive agent.Aza-BODIPY is a promising photosensitizer for use in PDT. Our results showed that aza-BODIPY-PDT induced apoptosis, probably through p53 and caspase3 in MCF-7 cells. Future studies should delineate the molecular mechanisms underlying aza-BODIPY-PDT-induced cell death for a better understanding of the signaling pathways modulated by the therapy so that this novel technology could be implemented in the clinic for treating breast cancer.

AIM

In this study,we aimed to determine the change in the expression levels of 88 carcinoma-associated genes induced by aza-BODIPY-PDT were analyzed so as to understand the specific pathways that are modulated by aza-BODIPY-PDT.

MATERIAL METHOD

In this study,the molecular basis of the anti-cancer activity of aza-BODIPY-PDT was investigated.Induction of apoptosis and necrosis in MCF-7 breast cancer cells after treatment with aza- BODIPY derivative with phthalonitrile substituents (aza-BODIPY) followed by light exposure was evaluated by Annexin V 7- Aminoactinomycin D (7-AAD) flow cytometry.

RESULTS

Aza-BODIPY-PDT induced cell death in MCF-7 cells treated with aza-BODIPY-PDT; flow cytometry revealed that 28 % of the cells died by apoptosis. Seven of the 88 carcinoma-associated genes that were assayed were differentially expressed -EGF, LEF1, WNT1, TCF7, and TGFBR2 were downregulated, and CASP3 and TP53 were upregulated - in cells subjected to aza-BODIPY-PDT.This made us think that the aza-BODIPY-PDT induced caspase 3 and p53-mediated apoptosis in MCF7 cells.

CONCLUSION

In our study,it was determined that the application of aza-BODIPY-PDT to MCF7 cells had a negative effect on cell connectivity and cell cycle.The fact that the same effect was not observed in control cells and MCF7 cells in the dark field of aza-BODIPY indicates that aza-BODIPY has a strong phodynamic anticancer effect.

Clonal and "Intrinsic" Heterogeneity of Somatic Variants in Microsatellite-Stable Colorectal Carcinomas and Their Metastases

To test the traditional model of tumor progression, Darwinian-type evolution, against the more recent Big Bang model, we selected 6 microsatellite-stable colorectal standard-type adenocarcinomas and their synchronous lymph node and liver metastases. Somatic genomic variants were identified by whole-exome sequencing (WES) of large tumor fragments from the primaries and 1 liver metastasis each, and used to design targeted resequencing next-generation sequencing (NGS) panels, 1 per case. Targeted deep resequencing (mean coverage, 2725; median, 2222) was performed with DNA from punch samples (1-mm tissue microarrayer needles) obtained from different regions of the primaries and their metastases. In total, 255 genomic variants were interrogated in 108 punch samples. Clonal heterogeneity was infrequent: a pattern of clonal heterogeneity consistent with a role in metastasis formation was observed only in 1 case in a single gene (p. Asp604Tyr of the PTPRT gene). However, when comparing variant allele frequencies (VAFs) of genomic variants in adjacent positions on chromosomes ("matched genomic variant loci") across punch samples, differences that exceeded 2 SD of the NGS assay variations (ad hoc dubbed VAF dysbalance) were observed in 7.1% of the punch samples (2.6%-12.0% per case), which indicates an intricate intermixing of mutated and nonmutated tumor cells ("intrinsic heterogeneity"). Additional OncoScan array analyses on a subset of the punch samples (31 in total) showed gross genomic aberrations as a possible explanation in only some (39.2%) of the matched genomic variant loci with VAF dysbalance. Our study provides a fairly direct (statistical model-free) view of the genomic states of microsatellite-stable colorectal carcinomas and their metastases, and suggests that Darwinian-type tumor evolution is not the key pathway of the metastasizing disease; instead, we recorded an "intrinsic" genomic heterogeneity, which may echo an initial Big Bang-like event.

CAR T-cells to treat brain tumors

Tremendous success using CAR T therapy in hematological malignancies has garnered significant interest in developing such treatments for solid tumors, including brain tumors. This success, however, has yet to be mirrored in solid organ neoplasms. CAR T function has shown limited efficacy against brain tumors due to several factors including the immunosuppressive tumor microenvironment, blood-brain barrier, and tumor-antigen heterogeneity. Despite these considerations, CAR T-cell therapy has the potential to be implemented as a treatment modality for brain tumors. Here, we review adult and pediatric brain tumors, including glioblastoma, diffuse midline gliomas, and medulloblastomas that continue to portend a grim prognosis. We describe insights gained from different preclinical models using CAR T therapy against various brain tumors and results gathered from ongoing clinical trials. Furthermore, we outline the challenges limiting CAR T therapy success against brain tumors and summarize advancements made to overcome these obstacles.

Uncovering novel mutational signatures by de novo extraction with SigProfilerExtractor

Mutational signature analysis is commonly performed in cancer genomic studies. Here, we present SigProfilerExtractor, an automated tool for de novo extraction of mutational signatures, and benchmark it against another 13 bioinformatics tools by using 34 scenarios encompassing 2,500 simulated signatures found in 60,000 synthetic genomes and 20,000 synthetic exomes. For simulations with 5% noise, reflecting high-quality datasets, SigProfilerExtractor outperforms other approaches by elucidating between 20% and 50% more true-positive signatures while yielding 5-fold less false-positive signatures. Applying SigProfilerExtractor to 4,643 whole-genome- and 19,184 whole-exome-sequenced cancers reveals four novel signatures. Two of the signatures are confirmed in independent cohorts, and one of these signatures is associated with tobacco smoking. In summary, this report provides a reference tool for analysis of mutational signatures, a comprehensive benchmarking of bioinformatics tools for extracting signatures, and several novel mutational signatures, including one putatively attributed to direct tobacco smoking mutagenesis in bladder tissues.

FAT2 mutation is associated with better prognosis and responsiveness to immunotherapy in uterine corpus endometrial carcinoma

BACKGROUND

Uterine corpus endometrial carcinoma (UCEC) ranks sixth among malignant tumors in women and the mortality is still rising. FAT2 gene has been considered to be related to the survival and prognosis of some certain diseases in previous studies, but the FAT2 mutation status in UCEC and its prognostic value has been rarely studied. Hence, the purpose of our study was to explore the role of FAT2 mutations for predicting prognosis and responsiveness to immunotherapy in patients with UCEC.

METHODS

UCEC samples from the Cancer Genome Atlas database were analyzed. We evaluated the impact of FAT2 gene mutation status and clinicopathological characteristics on the prognosis of UCEC patients and used univariate and multivariate Cox analysis risk scores to independently predict patient overall survival (OS). Tumor mutation burden (TMB) values of the FAT2 mutant and non-mutant groups were computed by Wilcoxon rank sum test. The correlation of FAT2 mutation and half maximal inhibitory concentration (IC50) values of various anticancer drugs was analyzed. Gene Ontology data and Gene Set Enrichment Analysis (GSEA) were employed to examine the differential expression of genes between the two groups. Finally, a single-sample GSEA arithmetic was utilized to measure the abundance of tumor-infiltrating immune cells in UCEC patients.

RESULTS

FAT2 mutations suggested better OS (p < 0.001) and disease-free survival (DFS) (p = 0.007) in UCEC. The IC50 values of 18 anticancer drugs were upregulated in FAT2 mutation patients (p < 0.05). The TMB and microsatellite instability values of patients with FAT2 mutations were significantly higher (p < 0.001). Next, the Kyoto Encyclopedia of Genes and Genomes functional analysis and GSEA revealed the potential mechanism of FAT2 mutation on the tumorigenesis and progression of UCEC. In addition, in reference to the UCEC microenvironment, the infiltration levels of activated CD4/CD8 T cells (p < 0.001/p = 0.001) and plasmacytoid dendritic cells (p = 0.006) were upregulated in the non-FAT2 mutation group, and Type 2 T helper cells (p = 0.001) were downregulated in the FAT2 mutation group.

CONCLUSIONS

UCEC patients with FAT2 mutations have better prognosis and are more likely to respond to immunotherapy. FAT2 mutation may be a valuable predictor for prognosis and responsiveness to immunotherapy in UCEC patients.

Somatic Mutation of FAT Family Genes Implicated Superior Prognosis in Patients With Stomach Adenocarcinoma

FAT family genes encode protocadherin, which regulates tumor cell proliferation and migration. Although transcriptional levels of FAT family members had been reported in multiple malignant tumors, the association between mutation and prognosis of the FAT family in stomach adenocarcinoma (STAD) has not been investigated. Herein, we performed a multi-omics integrative bioinformatics analysis using genomic and mRNA expression data to explore the role of gene mutations across the FAT family on clinical outcomes of STAD. The results showed that FAT mutations occurred in 174 of 435 (40%) of the samples. Patients with FAT mutations possessed significantly better progression-free survival (P = 0.019) and overall survival (P = 0.034) than those with non-FAT mutations, and FAT mutations exhibited significantly higher tumor mutational burden (TMB) and microsatellite instability. Notably, FAT mutations had a greater effect on somatic single-nucleotide variation than copy number variation and resulted in more abundant DNA damage repair (DDR) mutations. Further investigation demonstrated that FAT mutations contributed to an inflammatory tumor microenvironment (TME), as indicated by significantly increased numbers of activated CD4 and CD8 T cells, and significantly decreased numbers of mast cell, plasmacytoid dendritic cell, type 2 T helper cell, and high expression of immune-promoting genes. Moreover, biological process antigen processing and presentation, DNA replication, and DDR-related pathways were significantly upregulated in patients with FAT mutations. Collectively, FAT mutations significantly improved the survival of patients with STAD by enhancing tumor immunogenicity (e.g., TMB and DDR mutations) and an inflamed TME, indicating that the FAT family might be a potential prognostic and therapeutic biomarker for STAD.

Discovering Innate Driver Variants for Risk Assessment of Early Colorectal Cancer Metastasis

Metastasis is the main fatal cause of colorectal cancer (CRC). Although enormous efforts have been made to date to identify biomarkers associated with metastasis, there is still a huge gap to translate these efforts into effective clinical applications due to the poor consistency of biomarkers in dealing with the genetic heterogeneity of CRCs. In this study, a small cohort of eight CRC patients was recruited, from whom we collected cancer, paracancer, and normal tissues simultaneously and performed whole-exome sequencing. Given the exomes, a novel statistical parameter LIP was introduced to quantitatively measure the local invasion power for every somatic and germline mutation, whereby we affirmed that the innate germline mutations instead of somatic mutations might serve as the major driving force in promoting local invasion. Furthermore, via bioinformatic analyses of big data derived from the public zone, we identified ten potential driver variants that likely urged the local invasion of tumor cells into nearby tissue. Of them, six corresponding genes were new to CRC metastasis. In addition, a metastasis resister variant was also identified. Based on these eleven variants, we constructed a logistic regression model for rapid risk assessment of early metastasis, which was also deployed as an online server, AmetaRisk (http://www.bio-add.org/AmetaRisk). In summary, we made a valuable attempt in this study to exome-wide explore the genetic driving force to local invasion, which provides new insights into the mechanistic understanding of metastasis. Furthermore, the risk assessment model can assist in prioritizing therapeutic regimens in clinics and discovering new drug targets, and thus substantially increase the survival rate of CRC patients.

Salivary gland papillary adenocarcinoma with intestinal-like features: Clinicopathologic, immunohistochemical, and genetic study of six cases

BACKGROUND

Salivary gland tumors with papillary architecture and intestinal-like mucinous cytologic features are rare. Their clinicopathologic and genetic features are not fully understood, and whether they represent one separate entity remains unclear.

METHODS

Six salivary adenocarcinomas with papillary architecture and intestinal-like mucinous cytologic features were reported. Immunostaining was done for CK7, CK20, CDX2, SOX10, S100, MUC1, MUC2 and MUC5AC. Tumor DNA samples were extracted for Sanger sequencing. Previously reported morphology-analogous cases were reviewed.

RESULTS

Six cases involved the palate (2), retromolar region (1), submandibular region (1), tongue (1), and mandible (1). Five cases were followed up, with one case of recurrence one year after surgery, one death from cerebral infarction seven days after surgery, and three cases without signs of recurrence or metastasis over five years. All cases had abundant mucinous production and presented a typical immunophenotype common to salivary primaries, CK7 & MUC1 positive, CK20 & CDX-2 negative. Sanger sequencing demonstrated recurrent AKT1 E17K mutations in four cases (4/6, 66.7%). A review of reported salivary intestinal-like tumors revealed 3 out of 13 cases presented with papillary morphology and CDX-2 negative. Some salivary papillary neoplasms with mucinous cytologic features termed as intraductal papillary neoplasms or mucinous adenocarcinomas were also reported with AKT1 E17K mutations.

CONCLUSION

We describe 6 cases of salivary gland papillary adenocarcinoma with intestinal-like mucinous cytologic features, which are different from conventional intestinal-type adenocarcinoma, presenting a consistent immunophenotype of CK7 & MUC1 positive, CK20 & CDX-2 negative and exhibiting recurrent AKT1 E17K mutations.

Evolving Landscape of Long Non-coding RNAs in Cerebrospinal Fluid: A Key Role From Diagnosis to Therapy in Brain Tumors

Long non-coding RNAs (lncRNAs) are a type of non-coding RNAs that act as molecular fingerprints and modulators of many pathophysiological processes, particularly in cancer. Specifically, lncRNAs can be involved in the pathogenesis and progression of brain tumors, affecting stemness/differentiation, replication, invasion, survival, DNA damage response, and chromatin dynamics. Furthermore, the aberrations in the expressions of these transcripts can promote treatment resistance, leading to tumor recurrence. The development of next-generation sequencing technologies and the creation of lncRNA-specific microarrays have boosted the study of lncRNA etiology. Cerebrospinal fluid (CSF) directly mirrors the biological fluid of biochemical processes in the brain. It can be enriched for small molecules, peptides, or proteins released by the neurons of the central nervous system (CNS) or immune cells. Therefore, strategies that identify and target CSF lncRNAs may be attractive as early diagnostic and therapeutic options. In this review, we have reviewed the studies on CSF lncRNAs in the context of brain tumor pathogenesis and progression and discuss their potential as biomarkers and therapeutic targets.

Identification of the RP11-21C4.1/SVEP1 gene pair associated with FAT2 mutations as a potential biomarker in gastric cancer

Gastric cancer (GC) is one of the most common malignancies worldwide. Despite rapid advances in systemic therapy, GC remains the third leading cause of cancer-related deaths. We aimed to identify a novel prognostic signature associated with FAT2 mutations in GC. We analyzed the expression levels of FAT2-mutant and FAT2-wildtype GC samples obtained from The Cancer Genome Atlas (TCGA). The Kaplan–Meier survival curve showed that patients with FAT2 mutations showed better prognosis than those without the mutation. Sixteen long non-coding RNAs (lncRNAs) and 62 messenger RNAs (mRNAs) associated with FAT2 mutations were correlated with the prognosis of GC. We then constructed a 4-mRNA signature and a 5-lncRNA signature for GC. Finally, we identified the most relevant RP11-21 C4.1/SVEP1 gene pair as a prognostic signature of GC that exhibited superior predictive performance in comparison with the 4-mRNA or 5-lncRNA signature by weighted gene correlation network analysis (WGCNA) and Cox proportional hazards regression analysis. In this study, we constructed a prognostic signature of GC by integrative genomics analysis, which also provided insights into the molecular mechanisms linked to FAT2 mutations in GC.

Comparative Analysis of Patient-Matched PDOs Revealed a Reduction in OLFM4-Associated Clusters in Metastatic Lesions in Colorectal Cancer

Metastasis is the major cause of cancer-related death, but whether metastatic lesions exhibit the same cellular composition as primary tumors has yet to be elucidated. To investigate the cellular heterogeneity of metastatic colorectal cancer (CRC), we established 72 patient-derived organoids (PDOs) from 21 patients. Combined bulk transcriptomic and single-cell RNA-sequencing analysis revealed decreased gene expression of markers for differentiated cells in PDOs derived from metastatic lesions. Paradoxically, expression of potential intestinal stem cell markers was also decreased. We identified OLFM4 as the gene most strongly correlating with a stem-like cell cluster, and found OLFM4⁺ cells to be capable of initiating organoid culture growth and differentiation capacity in primary PDOs. These cells were required for the efficient growth of primary PDOs but dispensable for metastatic PDOs. These observations demonstrate that metastatic lesions have a cellular composition distinct from that of primary tumors; patient-matched PDOs are a useful resource for analyzing metastatic CRC.

•

Seventy-two PDOs were established from 21 stage IV CRC patients

•

Forty-one DEGs were identified between primary and corresponding metastatic PODs

•

scRNA-seq analysis identified OLFM4 as a potential cancer stem cell marker

•

Different roles of OLFM4⁺ cells in primary and metastatic PDOs were demonstrated

A genome-wide study of the relationship between chromosomal abnormalities and gene expression in colorectal tumors

The role of somatic copy number alterations (SCNAs) that occur in colorectal tumors is poorly understood. SCNAs are correlated with corresponding gene expression changes that may contribute to neoplastic progression. Thus, we examined SCNAs and the expression of messenger RNAs (mRNAs) located at corresponding loci in colorectal neoplasia, a progression model of human neoplasm. We used 42 colorectal neoplastic samples, including adenomas, intramucosal cancers (IMC) and invasive colorectal cancers (CRC) that were microsatellite stable (MSS) using a genome-wide SNP array and gene expression array (first cohort). In addition, validation analyses were examined (37 colorectal neoplasias). None of the mRNAs with a corresponding SCNA was found in the adenomas. However, three mRNAs, including ARFGEF2 at 20q13.13, N4BP2L2 at 13q13.1 and OLFM4 at 13q14.3 with a copy number (CN) gain at the corresponding locus were upregulated in IMCs of the first cohort. Moreover, upregulated expression of ARFGEF2 and OLFM4 was upregulated in the validation analysis. Finally, 28 mRNAs with gains of corresponding loci were pooled in invasive CRC of the first cohort. The mRNAs, including ACSS2 (20q11.22), DDX27 (20q13.13), MAPRE1 (20q11.21), OSBPL2 (20q11.22) and PHF20 (20q11.22-q11.23) with CN gains of the corresponding loci were identified in 28 mRNAs. Four of these mRNAs (DDX27, MAPRE1, OSBPL2 and PHF20) were upregulated in the invasive CRC in the validation analysis. We conclude that specific 13q and 22q CN gains with gene expression changes in the corresponding loci may play an important role in IMC cells' progression into invasive CRC.

Feature selection algorithm based on dual correlation filters for cancer-associated somatic variants

Background

Since the development of sequencing technology, an enormous amount of genetic information has been generated, and human cancer analysis using this information is drawing attention. As the effects of variants on human cancer become known, it is important to find cancer-associated variants among countless variants.

Results

We propose a new filter-based feature selection method applicable for extracting cancer-associated somatic variants considering correlations of data. Both variants associated with the activation and deactivation of cancer’s characteristics are analyzed using dual correlation filters. The multiobjective optimization is utilized to consider two types of variants simultaneously without redundancy. To overcome high computational complexity problem, we calculate the correlation-based weight to select significant variants instead of directly searching for the optimal subset of variants. The proposed algorithm is applied to the identification of melanoma metastasis or breast cancer stage, and the classification results of the proposed method are compared with those of conventional single correlation filter-based method.

Conclusions

We verified that the proposed dual correlation filter-based method can extract cancer-associated variants related to the characteristics of human cancer.

Genomic profiling of sporadic liver metastatic colorectal cancer

Sporadic colorectal cancer (sCRC) is the third leading cause of cancer death in the Western world. Approximately, a quarter of sCRC patients present metastatic dissemination at the moment of diagnosis, the liver being the most frequently affected organ. Additionally, this group of CRC patients is characterized by a worse prognosis. In the last decades, significant technological developments for genome analysis have fostered the identification and characterization of genetic alterations involved in the pathogenesis of sCRC. However, genetic alterations involved in the metastatic process through which tumor cells are able to colonize other tissues with a different microenvironment, still remain to be fully identified. Here, we review current knowledge about the most relevant genomic alterations involved in the liver metastatic process of sCRC, including detailed information about the genetic profile of primary colorectal tumors vs. their paired liver metastases.

Identification of FBXW7α-regulated genes in M1-polarized macrophages in colorectal cancer by RNA sequencing

OBJECTIVES

To determine the FBXW7α-regulated genes in tumor-polarized macrophages in colorectal cancer.

METHODS

This experimental study was performed between June 2017 and March 2019. FBXW7α siRNA transfected RAW264.7 cells, together with the control group, were co-cultured with the colon cancer cell line, Colon-26. M1  marker production from the macrophages was determined by ELISA and quantitative reverse transcription-polymerase chain reaction. Whole genomic differential expression between the FBXW7α siRNA group and the control group were determined by RNA-sequencing analysis. The target site of the microRNA-205 gene was predicted using Targetscan and was verified by the luciferase assay. By transfecting mimics or inhibitors of microRNA-205, we explored the role of FBXW7α/microRNA-205 axis in regulating the polarization of tumor-associated macrophages (TAM).

RESULTS

FBXW7α knockdown in RAW264.7 enhanced the expression of cyclooxigenase (COX)-2 and inducible nitric oxide synthase (iNOS), mRNA expression and IL6, IL12, p40, and tumor necrosis factor-α (TNFα) production upon co-culture with Colon-26 cells in vitro. Further, compared with the control group, 648 genes in total were enhanced and 416 targets were downregulated in FBXW7α siRNA transfected cells, among which miR-205 was the most significantly upregulated. SMAD1 was identified as an miR-205 target. The FBXW7α/miR-205 axis might regulate TAM polarization by affecting SMAD1 expression.

CONCLUSION

These results prove that the FBXW7α/miR-205 axis plays an important role in TAM polarization and could facilitate further exploration of its molecular mechanism.

Chromosomal Instability in Tumor Initiation and Development

Chromosomal instability (CIN) is one of the major forms of genomic instability in various human cancers and is recognized as a common hallmark of tumorigenesis and heterogeneity. However, some malignant tumors show a paucity of chromosomal alterations, suggesting that tumor progression and evolution can occur in the absence of CIN. It is unclear whether CIN is stable between precursor lesions, primary tumor, and metastases or if it evolves during these steps. In this review, we describe the influence of CIN on the various steps in tumor initiation and development. Given the recognized significant effects of CIN in cancer, CIN-targeted therapeutics could have a major impact on improving clinical outcomes.

HCMDB: the human cancer metastasis database

Metastasis is the main event leading to death in cancer patients. Over the past decade, high-throughput technologies have provided genome-wide view of transcriptomic changes associated with cancer metastases. Many microarray and RNA sequencing studies have addressed metastases-related expression patterns in various types of cancer, and the number of relevant works continues to increase rapidly. These works have characterized genes that orchestrate the metastatic phenotype of cancer cells. However, these expression data have been deposited in various repositories, and efficiently analyzing these data is still difficult because of the lack of an integrated data mining platform. To facilitate the in-depth analyses of transcriptome data on metastasis, it is quite important to make a comprehensive integration of these metastases-related expression data. Here, we presented a database, HCMDB (the human cancer metastasis database, http://hcmdb.i-sanger.com/index), which is freely accessible to the research community query cross-platform transcriptome data on metastases. HCMDB is developed and maintained as a useful resource for building the systems-biology understanding of metastasis.

The landscape of genomic copy number alterations in colorectal cancer and their consequences on gene expression levels and disease outcome

Aneuploidy, the unbalanced state of the chromosome content, represents a hallmark of most solid tumors, including colorectal cancer. Such aneuploidies result in tumor specific genomic imbalances, which emerge in premalignant precursor lesions. Moreover, increasing levels of chromosomal instability have been observed in adenocarcinomas and are maintained in distant metastases. A number of studies have systematically integrated copy number alterations with gene expression changes in primary carcinomas, cell lines, and experimental models of aneuploidy. In fact, chromosomal aneuploidies target a number of genes conferring a selective advantage for the metabolism of the cancer cell. Copy number alterations not only have a positive correlation with expression changes of the majority of genes on the altered genomic segment, but also have effects on the transcriptional levels of genes genome-wide. Finally, copy number alterations have been associated with disease outcome; nevertheless, the translational applicability in clinical practice requires further studies. Here, we (i) review the spectrum of genetic alterations that lead to colorectal cancer, (ii) describe the most frequent copy number alterations at different stages of colorectal carcinogenesis, (iii) exemplify their positive correlation with gene expression levels, and (iv) discuss copy number alterations that are potentially involved in disease outcome of individual patients.

Quantitative evidence for early metastatic seeding in colorectal cancer

Both the timing and molecular determinants of metastasis are unknown, hindering treatment and prevention efforts. Here we characterize the evolutionary dynamics of this lethal process by analyzing exome-sequencing data from 118 biopsies from 23 patients with colorectal cancer with metastases to the liver or brain. The data show that the genomic divergence between the primary tumor and metastasis is low and that canonical driver genes were acquired early. Analysis within a spatial tumor growth model and statistical inference framework indicates that early disseminated cells commonly (81%, 17 out of 21 evaluable patients) seed metastases while the carcinoma is clinically undetectable (typically, less than 0.01 cm3). We validated the association between early drivers and metastasis in an independent cohort of 2,751 colorectal cancers, demonstrating their utility as biomarkers of metastasis. This conceptual and analytical framework provides quantitative in vivo evidence that systemic spread can occur early in colorectal cancer and illuminates strategies for patient stratification and therapeutic targeting of the canonical drivers of tumorigenesis.

Engineering Multidimensional Evolutionary Forces to Combat Cancer

With advances in technology and bioinformatics, we are now positioned to view and manage cancer through an evolutionary lens. This perspective is critical as our appreciation for the role of tumor heterogeneity, tumor immune compartment, and tumor microenvironment on cancer pathogenesis and evolution grows. Here, we explore recent knowledge on the evolutionary basis of cancer pathogenesis and progression, viewing tumors as multilineage, multicomponent organisms whose growth is regulated by subcomponent fitness relationships. We propose reconsidering some current tenets of the cancer management paradigm in order to take better advantage of crucial fitness relationships to improve outcomes of patients with cancer. SIGNIFICANCE: Tumor and tumor immune compartment and microenvironment heterogeneity, and their evolution, are critical disease features that affect treatment response. The impact and interplay of these components during treatment are viable targets to improve clinical response. In this article, we consider how tumor cells, the tumor immune compartment and microenvironment, and epigenetic factors interact and also evolve during treatment. We evaluate the convergence of these factors and suggest innovative treatment concepts that leverage evolutionary relationships to limit tumor growth and drug resistance.

High-resolution clonal mapping of multi-organ metastasis in triple negative breast cancer

Most triple negative breast cancers (TNBCs) are aggressively metastatic with a high degree of intra-tumoral heterogeneity (ITH), but how ITH contributes to metastasis is unclear. Here, clonal dynamics during metastasis were studied in vivo using two patient-derived xenograft (PDX) models established from the treatment-naive primary breast tumors of TNBC patients diagnosed with synchronous metastasis. Genomic sequencing and high-complexity barcode-mediated clonal tracking reveal robust alterations in clonal architecture between primary tumors and corresponding metastases. Polyclonal seeding and maintenance of heterogeneous populations of low-abundance subclones is observed in each metastasis. However, lung, liver, and brain metastases are enriched for an identical population of high-abundance subclones, demonstrating that primary tumor clones harbor properties enabling them to seed and thrive in multiple organ sites. Further, clones that dominate multi-organ metastases share a genomic lineage. Thus, intrinsic properties of rare primary tumor subclones enable the seeding and colonization of metastases in secondary organs in these models.

It is unclear how intra-tumoral heterogeneity contributes to metastasis. Here the authors study the clonal dynamics of triple negative breast cancer metastasis using patient derived xenografts and demonstrate that primary tumor clones harbor properties that support seeding and colonization of multiple organs.

Whole genome sequencing puts forward hypotheses on metastasis evolution and therapy in colorectal cancer

Incomplete understanding of the metastatic process hinders personalized therapy. Here we report the most comprehensive whole-genome study of colorectal metastases vs. matched primary tumors. 65% of somatic mutations originate from a common progenitor, with 15% being tumor- and 19% metastasis-specific, implicating a higher mutation rate in metastases. Tumor- and metastasis-specific mutations harbor elevated levels of BRCAness. We confirm multistage progression with new components ARHGEF7/ARHGEF33. Recurrently mutated non-coding elements include ncRNAs RP11-594N15.3, AC010091, SNHG14, 3’ UTRs of FOXP2, DACH2, TRPM3, XKR4, ANO5, CBL, CBLB, the latter four potentially dual protagonists in metastasis and efferocytosis-/PD-L1 mediated immunosuppression. Actionable metastasis-specific lesions include FAT1, FGF1, BRCA2, KDR, and AKT2-, AKT3-, and PDGFRA-3’ UTRs. Metastasis specific mutations are enriched in PI3K-Akt signaling, cell adhesion, ECM and hepatic stellate activation genes, suggesting genetic programs for site-specific colonization. Our results put forward hypotheses on tumor and metastasis evolution, and evidence for metastasis-specific events relevant for personalized therapy.

The evolution and genetic nature of metastatic lesions is not completely characterized. Here the authors perform a comprehensive whole-genome study of colorectal metastases in comparison to matched primary tumors and define a multistage progression model and metastasis-specific changes that, in part, are therapeutically actionable.

Comparison of EpCAMhighCD44+ cancer stem cells with EpCAMhighCD44- tumor cells in colon cancer by single-cell sequencing

Cancer stem cells (CSCs) are considered to be responsible for tumorigenesis and cancer relapse. EpCAM^highCD44⁺ tumor cells are putative colorectal CSCs that express high levels of stem cell genes, while the EpCAM^highCD44^- population mostly contains differentiated tumor cells (DTCs). This study aims to determine whether single CSC (EpCAM^highCD44⁺) and DTC (EpCAM^highCD44^-) can be distinguished in terms of somatic copy number alterations (SCNAs). We applied fluorescence-activated cell sorting to isolate the CD45^-EpCAM^highCD44⁺ and CD45^-EpCAM^highCD44^- populations from two primary colon tumors, on which low-coverage single-cell whole-genome sequencing (WGS) was then performed ∼0.1x depth. We compared the SCNAs of the CSCs and DTCs at single-cell resolution. In total, 47 qualified single cells of the two populations underwent WGS. The single-cell SCNA profiles showed that there were obvious SCNAs in both the CSCs and DTCs of each patient, and each patient had a specific copy number alteration pattern. Hierarchical clustering and correlation analysis both showed that the SCNA profiles of CSCs and DTCs from the same patient had similar SCNA pattern, while there were regional differences in the CSCs and DTCs in certain patient. SCNAs of CSCs in the same patient were highly reproducible. Our data suggest that major SCNAs occurred at an early stage and were inherited steadily. The similarity of ubiquitous SCNAs between the CSCs and DTCs might have arisen from lineage differentiation. CSCs from the same patient had reproducible SCNA profiles, indicating that gain or loss in certain chromosome is required for colon cancer development.

Genetic insights into the morass of metastatic heterogeneity

Tumour heterogeneity poses a substantial problem for the clinical management of cancer. Somatic evolution of the cancer genome results in genetically distinct subclones in the primary tumour with different biological properties and therapeutic sensitivities. The problem of heterogeneity is compounded in metastatic disease owing to the complexity of the metastatic process and the multiple biological hurdles that the tumour cell must overcome to establish a clinically overt metastatic lesion. New advances in sequencing technology and clinical sample acquisition are providing insights into the phylogenetic relationship of metastases and primary tumours at the level of somatic tumour genetics while also illuminating fundamental mechanisms of the metastatic process. In addition to somatically acquired genetic heterogeneity in the tumour cells, inherited population-based genetic heterogeneity can profoundly modify metastatic biology and further complicate the development of effective, broadly applicable antimetastatic therapies. Here, we examine how genetic heterogeneity impacts metastatic disease and the implications of current knowledge for future research endeavours and therapeutic interventions.

Genomic heterogeneity in primary colorectal carcinomas and their metastases: born bad or brought up a villain?

Progression of solid cancers, colorectal carcinomas among them, from their primaries to metastatic lesions traditionally is thought to proceed by a stepwise acquisition of and selection for genomic aberrations. To test if patterns of genomic aberrations would be consistent with this model, we studied 10 colorectal carcinoma primary-metastasis pairs, 9 with 1 liver metastasis each and 1 with 2 metastases. Next-generation targeted sequencing (50-gene panel) with samples obtained from different regions of the primaries and their metastases demonstrated 1-11 gene mutations per lesion. But only in 2 tumors were there seen mutations in all samples from the metastasis and not any of the primaries (BRAF^D594N and SMARCB1^R377C mutation, respectively). However, allelotyping the multiregional samples with polymorphous microsatellite markers (17p13.1, D9S942, D9S1748, D5S346, D5S1385) and DNA methylation studies with a marker panel (MLH1, CDNK2A, NEUROG1, CRABP1, CACNA1G, IGF2, RUNX3, SOCS1) showed remarkably "insular" genomic aberrations in all cases for at least some of the analyses. The marked preponderance of mutations shared by the primaries and their metastases throughout the lesions over mutations private to metastases suggests that, at least in many cases, colorectal carcinomas might be endowed with a mutational load sufficient for fully fledged metastases even at a very early stage ("born bad"). But the very focal allelic imbalances and methylations observed here, hypothetically, could play a role in clinically metastasizing disease, a process of years rather than months and very much a matter of tumor-host interactions when tumor cells adapt to the host microenvironment.

Exome sequencing of synchronously resected primary colorectal tumours and colorectal liver metastases to inform oncosurgical management

BACKGROUND

Next generation sequencing technology has facilitated mapping of the colorectal cancer genotype and furthered our understanding of metastogenesis. The aim of this study was to investigate for conserved and different mutations in the exomes of synchronously resected primary colorectal tumour and liver metastases. This information could potentially be utilised to guide the treatment of advanced disease with the help of biological information from the primary tumour.

METHODS

We performed exome sequencing of synchronously resected primary colorectal cancer and colorectal liver metastases as well as normal colonic mucosa and liver parenchyma, from four patients who had received neo-adjuvant chemotherapy, at a depth of 50X using the Ion Proton platform. Raw data was mapped to the reference genome prior to variant calling, annotation and downstream analysis.

RESULTS

Exome sequencing identified 585 non-synonymous missense single nucleotide variants (SNVs), of which 215 (36.8%) were unique to the primary tumour, 226 (38.6%) unique to the metastasis and 81 (13.8%) present in patient matched pairs. SNVs identified in the ErbB pathway appear to be concordant between primary and metastatic tumours.

CONCLUSION

Only 13.8% of the metastatic exome can be predicted by the genotype of the primary tumour. We have demonstrated concordance of a number of SNVs in the ErbB pathway, which may inform selection of therapeutic agents in advanced colorectal cancer.

Molecular signatures reflecting microenvironmental metabolism and chemotherapy-induced immunogenic cell death in colorectal liver metastases

Background

Metastatic colorectal cancer (CRC) is associated with highly variable clinical outcome and response to therapy. The recently identified consensus molecular subtypes (CMS1-4) have prognostic and therapeutic implications in primary CRC, but whether these subtypes are valid for metastatic disease is unclear. We performed multi-level analyses of resectable CRC liver metastases (CLM) to identify molecular characteristics of metastatic disease and evaluate the clinical relevance.

Methods

In this ancillary study to the Oslo-CoMet trial, CLM and tumor-adjacent liver tissue from 46 patients were analyzed by profiling mutations (targeted sequencing), genome-wide copy number alteration (CNAs), and gene expression.

Results

Somatic mutations and CNAs detected in CLM were similar to reported primary CRC profiles, while CNA profiles of eight metastatic pairs suggested intra-patient divergence. A CMS classifier tool applied to gene expression data, revealed the cohort to be highly enriched for CMS2. Hierarchical clustering of genes with highly variable expression identified two subgroups separated by high or low expression of 55 genes with immune-related and metabolic functions. Importantly, induction of genes and pathways associated with immunogenic cell death (ICD) was identified in metastases exposed to neoadjuvant chemotherapy (NACT).

Conclusions

The uniform classification of CLM by CMS subtyping may indicate that novel class discovery approaches need to be explored to uncover clinically useful stratification of CLM. Detected gene expression signatures support the role of metabolism and chemotherapy in shaping the immune microenvironment of CLM. Furthermore, the results point to rational exploration of immune modulating strategies in CLM, particularly by exploiting NACT-induced ICD.

Sequencing brain metastases and opportunities for targeted therapies

CNS metastases have long been recognized as a common and late complication of systemic malignancies. They represent the most common tumor of the brain. As outcomes and overall survival improve with better tolerated and more durable responses from therapies for systemic cancers, the incidence and prevalence of brain metastases is likely to increase. Among the most common systemic cancers leading to brain metastases include lung, melanoma, breast (triple-negative histology) and renal cell cancers. To date, there has been infrequent involvement of gastrointestinal and gynecologic malignancies; however, this may also change, reflecting improvement in overall survival and therapeutic regimens. Traditional therapy of brain metastases has focused on surgery, radiation therapy or best supportive/palliative care. The advent of modern genomic techniques, including next-generation and whole-exome sequencing, has allowed for the identification of unique markers and potential drivers of metastatic pathways. This review aims to discuss and highlight the known drivers of disease and the opportunities for ultimate development of targeted therapies.

DNA copy number changes define spatial patterns of heterogeneity in colorectal cancer

Genetic heterogeneity between and within tumours is a major factor determining cancer progression and therapy response. Here we examined DNA sequence and DNA copy-number heterogeneity in colorectal cancer (CRC) by targeted high-depth sequencing of 100 most frequently altered genes. In 97 samples, with primary tumours and matched metastases from 27 patients, we observe inter-tumour concordance for coding mutations; in contrast, gene copy numbers are highly discordant between primary tumours and metastases as validated by fluorescent in situ hybridization. To further investigate intra-tumour heterogeneity, we dissected a single tumour into 68 spatially defined samples and sequenced them separately. We identify evenly distributed coding mutations in APC and TP53 in all tumour areas, yet highly variable gene copy numbers in numerous genes. 3D morpho-molecular reconstruction reveals two clusters with divergent copy number aberrations along the proximal-distal axis indicating that DNA copy number variations are a major source of tumour heterogeneity in CRC.

Molecular evolution of colorectal cancer: from multistep carcinogenesis to the big bang

Colorectal cancer is characterized by exquisite genomic instability either in the form of microsatellite instability or chromosomal instability. Microsatellite instability is the result of mutation of mismatch repair genes or their silencing through promoter methylation as a consequence of the CpG island methylator phenotype. The molecular causes of chromosomal instability are less well characterized. Genomic instability and field cancerization lead to a high degree of intratumoral heterogeneity and determine the formation of cancer stem cells and epithelial-mesenchymal transition mediated by the TGF-β and APC pathways. Recent analyses using integrated genomics reveal different phases of colorectal cancer evolution. An initial phase of genomic instability that yields many clones with different mutations (big bang) is followed by an important, previously not detected phase of cancer evolution that consists in the stabilization of several clones and a relatively flat outgrowth. The big bang model can best explain the coexistence of several stable clones and is compatible with the fact that the analysis of the bulk of the primary tumor yields prognostic information.

Unexpected effects of different genetic backgrounds on identification of genomic rearrangements via whole-genome next generation sequencing

Background

Whole genome next generation sequencing (NGS) is increasingly employed to detect genomic rearrangements in cancer genomes, especially in lymphoid malignancies. We recently established a unique mouse model by specifically deleting a key non-homologous end-joining DNA repair gene, Xrcc4, and a cell cycle checkpoint gene, Trp53, in germinal center B cells. This mouse model spontaneously develops mature B cell lymphomas (termed G1XP lymphomas).

Results

Here, we attempt to employ whole genome NGS to identify novel structural rearrangements, in particular inter-chromosomal translocations (CTXs), in these G1XP lymphomas. We sequenced six lymphoma samples, aligned our NGS data with mouse reference genome (in C57BL/6J (B6) background) and identified CTXs using CREST algorithm. Surprisingly, we detected widespread CTXs in both lymphomas and wildtype control samples, majority of which were false positive and attributable to different genetic backgrounds. In addition, we validated our NGS pipeline by sequencing multiple control samples from distinct tissues of different genetic backgrounds of mouse (B6 vs non-B6). Lastly, our studies showed that widespread false positive CTXs can be generated by simply aligning sequences from different genetic backgrounds of mouse.

Conclusions

We conclude that mapping and alignment with reference genome might not be a preferred method for analyzing whole-genome NGS data obtained from a genetic background different from reference genome. Given the complex genetic background of different mouse strains or the heterogeneity of cancer genomes in human patients, in order to minimize such systematic artifacts and uncover novel CTXs, a preferred method might be de novo assembly of personalized normal control genome and cancer cell genome, instead of mapping and aligning NGS data to mouse or human reference genome. Thus, our studies have critical impact on the manner of data analysis for cancer genomics.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3153-9) contains supplementary material, which is available to authorized users.

Spatio-temporal mutation profiles of case-matched colorectal carcinomas and their metastases reveal unique de novo mutations in metachronous lung metastases by targeted next generation sequencing

Background

Targeted next generation sequencing (tNGS) has become part of molecular pathology diagnostics for determining RAS mutation status in colorectal cancer (CRC) patients as predictive tool for decision on EGFR-targeted therapy. Here, we investigated mutation profiles of case-matched tissue specimens throughout the disease course of CRC, to further specify RAS-status dynamics and to identify de novo mutations associated with distant metastases.

Methods

Case-matched formalin-fixed and paraffin-embedded (FFPE) resection specimens (n = 70; primary tumours, synchronous and/or metachronous liver and/or lung metastases) of 14 CRC cases were subjected to microdissection of normal colonic epithelial, primary and metastatic tumour cells, their DNA extraction and an adapted library protocol for limited DNA using the 48 gene TruSeq Amplicon Cancer Panel^TM, MiSeq sequencing and data analyses (Illumina).

Results

By tNGS primary tumours were RAS wildtype in 5/14 and mutated in 9/14 (8/9 KRAS exon 2; 1/9 NRAS Exon 3) of cases. RAS mutation status was maintained in case-matched metastases throughout the disease course, albeit with altered allele frequencies. Case-matched analyses further identified a maximum of three sequence variants (mainly in APC, KRAS, NRAS, TP53) shared by all tumour specimens throughout the disease course per individual case. In addition, further case-matched de novo mutations were detected in synchronous and/or metachronous liver and/or lung metastases (e.g. in APC, ATM, FBXW7, FGFR3, GNAQ, KIT, PIK3CA, PTEN, SMAD4, SMO, STK11, TP53, VHL). Moreover, several de novo mutations were more frequent in synchronous (e.g. ATM, KIT, PIK3CA, SMAD4) or metachronous (e.g. FBXW7, SMO, STK11) lung metastases. Finally, some de novo mutations occurred only in metachronous lung metastases (CDKN2A, FGFR2, GNAS, JAK3, SRC).

Conclusion

Together, this study employs an adapted FFPE-based tNGS approach to confirm conservation of RAS mutation status in primary and metastatic tissue specimens of CRC patients. Moreover, it identifies genes preferentially mutated de novo in late disease stages of metachronous CRC lung metastases, several of which might be actionable by targeted therapies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12943-016-0549-8) contains supplementary material, which is available to authorized users.

Genome Analysis of Osteosarcoma Progression Samples Identifies FGFR1 Overexpression as a Potential Treatment Target and CHM as a Candidate Tumor Suppressor Gene

Osteosarcoma (OS) is the most common primary malignant tumor of bone, showing complex chromosomal rearrangements but with few known consistent changes. Deeper biological understanding is crucial to find new therapies to improve patient survival. We have sequenced the whole exome of two primary tumors (before and after chemotherapy), one metastatic tumor and a matched normal sample from two OS patients, to identify mutations involved in cancer biology. The metastatic samples were also RNA sequenced. By RNA sequencing we identified dysregulated expression levels of drug resistance- and apoptosis-related genes. Two fusion transcripts were identified in one patient (OS111); the first resulted in p53 inactivation by fusing the first exon of TP53 to the fifth exon of FAM45A. The second fusion joined the two first exons of FGFR1 to the second exon of ZNF343. Furthermore, FGFR1 was amplified and highly expressed, representing a potential treatment target in this patient. Whole exome sequencing revealed large intertumor heterogeneity, with surprisingly few shared mutations. Careful evaluation and validation of the data sets revealed a number of artefacts, but one recurrent mutation was validated, a nonsense mutation in CHM (patient OS106), which also was the mutation with the highest expression frequency (53%). The second patient (OS111) had wild-type CHM, but a downregulated expression level. In a panel of 71 clinical samples, we confirmed significant low expression of CHM compared to the controls (p = 0.003). Furthermore, by analyzing public datasets, we identified a significant association between low expression and poor survival in two other cancer types. Together, these results suggest CHM as a candidate tumor suppressor gene that warrants further investigation.

Genome-wide mutation profiles of colorectal tumors and associated liver metastases at the exome and transcriptome levels

To characterize the mutation profiles of colorectal cancer (CRC) primary tumors (PTs) and liver metastases (CLMs), we performed both whole-exome and RNA sequencing. Ten significantly mutated genes, including BMI1, CARD11, and NRG1, were found in 34 CRCs with CLMs. We defined three mutation classes (Class 1 to 3) based on the absence or presence of mutations during liver metastasis. Most mutations were classified into Class 1 (shared between PTs and CLMs), suggesting the common clonal origin of PTs and CLMs. Class 1 was more strongly associated with the clinical characteristics of advanced cancer and was more frequently superimposed with chromosomal deletions in CLMs than Class 2 (PT-specific). The integration of exome and RNA sequencing revealed that variant-allele frequencies (VAFs) of mutations in the transcriptome tended to have stronger functional implications than those in the exome. For instance, VAFs of the TP53 and APC mutations in the transcriptome significantly correlated with the expression level of their target genes. Additionally, mutations with high functional impact were enriched with high VAFs in the CLM transcriptomes. We identified 11 mutation-associated splicing events in the CRC transcriptomes. Thus, the integration of the exome and the transcriptome may elucidate the underlying molecular events responsible for CLMs.

Pathophysiology of colorectal peritoneal carcinomatosis: Role of the peritoneum

Colorectal cancer (CRC) is the third most common cancer and the fourth most common cause of cancer-related death worldwide. Besides the lymphatic and haematogenous routes of dissemination, CRC frequently gives rise to transcoelomic spread of tumor cells in the peritoneal cavity, which ultimately leads to peritoneal carcinomatosis (PC). PC is associated with a poor prognosis and bad quality of life for these patients in their terminal stages of disease. A loco-regional treatment modality for PC combining cytoreductive surgery and hyperthermic intraperitoneal peroperative chemotherapy has resulted in promising clinical results. However, this novel approach is associated with significant morbidity and mortality. A comprehensive understanding of the molecular events involved in peritoneal disease spread is paramount in avoiding unnecessary toxicity. The emergence of PC is the result of a molecular crosstalk between cancer cells and host elements, involving several well-defined steps, together known as the peritoneal metastatic cascade. Individual or clumps of tumor cells detach from the primary tumor, gain access to the peritoneal cavity and become susceptible to the regular peritoneal transport. They attach to the distant peritoneum, subsequently invade the subperitoneal space, where angiogenesis sustains proliferation and enables further metastatic growth. These molecular events are not isolated events but rather a continuous and interdependent process. In this manuscript, we review current data regarding the molecular mechanisms underlying the development of colorectal PC, with a special focus on the peritoneum and the role of the surgeon in peritoneal disease spread.

Defining the Influence of Germline Variation on Metastasis Using Systems Genetics Approaches

Cancer is estimated to be responsible for 8 million deaths worldwide and over half a million deaths every year in the United States. The majority of cancer-related deaths in solid tumors is directly associated with the effects of metastasis. While the influence of germline factors on cancer risk and development has long been recognized, the contribution of hereditary variation to tumor progression and metastasis has only gained acceptance more recently. A variety of approaches have been used to define how hereditary variation influences tumor progression and metastasis. One approach that garnered much early attention was epidemiological studies of cohorts of cancer patients, which demonstrated that specific loci within the human genome are associated with a differential propensity for aggressive tumor development. However, a powerful, and somewhat underutilized approach has been the use of systems genetics approaches in transgenic mouse models of human cancer. Such approaches are typically multifaceted, and involve integration of multiple lines of evidence derived, for example, from genetic and transcriptomic screens of genetically diverse mouse models of cancer, coupled with bioinformatics analysis of human cancer datasets, and functional analysis of candidate genes. These methodologies have allowed for the identification of multiple hereditary metastasis susceptibility genes, with wide-ranging cellular functions including regulation of gene transcription, cell proliferation, and cell-cell adhesion. In this chapter, we review how each of these approaches have facilitated the identification of these hereditary metastasis modifiers, the molecular functions of these metastasis-associated genes, and the implications of these findings upon patient survival.

TNFRSF10C copy number variation is associated with metastatic colorectal cancer

BACKGROUND

Genetic markers for distant metastatic disease in patients with colorectal cancer (CRC) are not well defined. Identification of genetic alterations associated with metastatic CRC could help to guide systemic and local treatment strategies. We evaluated the association of tumor necrosis factor receptor superfamily member 10C (TNFRSF10C) copy number variation (CNV) with distant metastatic disease in patients with CRC using The Cancer Genome Atlas (TCGA).

METHODS

Genetic sequencing data and clinical characteristics were obtained from TCGA for all available patients with CRC. There were 515 CRC patient samples with CNV and clinical outcome data, including a subset of 144 rectal adenocarcinoma patient samples. Using the TCGA CRC dataset, CNV of TNFRSF10C was evaluated for association with distant metastatic disease (M1 vs. M0). Multivariate logistic regression analysis with odds ratio (OR) using a 95% confidence interval (CI) was performed adjusting for age, T stage, N stage, adjuvant chemotherapy, gender, microsatellite instability (MSI), location, and surgical margin status.

RESULTS

TNFRSF10C CNV in patients with CRC was associated with distant metastatic disease [OR 4.81 (95% CI, 2.13-10.85) P<0.001] and positive lymph nodes [OR 18.83 (95% CI, 8.42-42.09)]; P<0.001) but not MSI (OR P=0.799). On multivariate analysis, after adjusting for pathologic T stage, N stage, adjuvant chemotherapy, gender, and MSI, TNFRSF10C CNV remained significantly associated with distant metastatic disease (OR P=0.018). Subset analysis revealed that TNFRSF10C CNV was also significantly associated with distant metastatic disease in patients with rectal adenocarcinoma (OR P=0.016).

CONCLUSIONS

TNFRSF10C CNV in patients with CRC is associated with distant metastatic disease. With further validation, such genetic profiles could be used clinically to support optimal systemic treatment strategies versus more aggressive local therapies in patients with CRC, including radiation therapy for rectal adenocarcinoma.

Whom to blame for metastasis, the epithelial-mesenchymal transition or the tumor microenvironment?

Changes in the tumor microenvironment (TME) can trigger the activation of otherwise non-malignant cells to become highly aggressive and motile. This is evident during initial tumor growth when the poor vascularization in tumors generates hypoxic regions that trigger the latent embryonic program, epithelial-to-mesenchymal transition (EMT), in epithelial carcinoma cells (e-cars) leading to highly motile mesenchymal-like carcinoma cells (m-cars), which also acquire cancer stem cell properties. After that, specific bidirectional interactions take place between m-cars and the cellular components of TME at different stages of metastasis. These interactions include several vicious positive feedback loops in which m-cars trigger a phenotypic switch, causing normal stromal cells to become pro-tumorigenic, which then further promote the survival, motility, and proliferation of m-cars. Accordingly, there is not a single culprit accounting for metastasis. Instead both m-cars and the TME dynamically interact, evolve and promote metastasis. In this review, we discuss the current status of the known interactions between m-cars and the TME during different stages of metastasis and how these interactions promote the metastatic activity of highly malignant m-cars by promoting their invasive mesenchymal phenotype and CSC properties.

Tumor evolution and intratumor heterogeneity in colorectal carcinoma: insights from comparative genomic profiling of primary tumors and matched metastases

Metastatic colorectal cancer (CRC) is one of the leading causes of cancer-related mortality among men and women worldwide. Over the past few decades, advances in our understanding of the genetic and epigenetic underpinnings of CRC have led to important insights into the pathogenesis of invasive tumors and have identified different molecular subgroups. Nonetheless, the events that might facilitate dissemination of tumor cells to distant sites giving rise to metastatic disease are not well characterized. Furthermore, in contrast to intertumor heterogeneity the extent of intratumor heterogeneity in different types of CRC has not been fully defined. In this paper, we review studies that have compared the genetic profile of primary invasive carcinomas to that of matched metastases and discuss the implications of their findings for our understanding of tumor evolution and for the clinical management of patients with advanced CRC.

Whole Exome Sequencing of Rapid Autopsy Tumors and Xenograft Models Reveals Possible Driver Mutations Underlying Tumor Progression

Pancreatic Ductal Adenocarcinoma (PDAC) is a highly lethal malignancy due to its propensity to invade and rapidly metastasize and remains very difficult to manage clinically. One major hindrance towards a better understanding of PDAC is the lack of molecular data sets and models representative of end stage disease. Moreover, it remains unclear how molecularly similar patient-derived xenograft (PDX) models are to the primary tumor from which they were derived. To identify potential molecular drivers in metastatic pancreatic cancer progression, we obtained matched primary tumor, metastases and normal (peripheral blood) samples under a rapid autopsy program and performed whole exome sequencing (WES) on tumor as well as normal samples. PDX models were also generated, sequenced and compared to tumors. Across the matched data sets generated for three patients, there were on average approximately 160 single-nucleotide mutations in each sample. The majority of mutations in each patient were shared among the primary and metastatic samples and, importantly, were largely retained in the xenograft models. Based on the mutation prevalence in the primary and metastatic sites, we proposed possible clonal evolution patterns marked by functional mutations affecting cancer genes such as KRAS, TP53 and SMAD4 that may play an important role in tumor initiation, progression and metastasis. These results add to our understanding of pancreatic tumor biology, and demonstrate that PDX models derived from advanced or end-stage likely closely approximate the genetics of the disease in the clinic and thus represent a biologically and clinically relevant pre-clinical platform that may enable the development of effective targeted therapies for PDAC.

Genomic Characterization of Brain Metastases Reveals Branched Evolution and Potential Therapeutic Targets

Brain metastases are associated with a dismal prognosis. Whether brain metastases harbor distinct genetic alterations beyond those observed in primary tumors is unknown. We performed whole-exome sequencing of 86 matched brain metastases, primary tumors, and normal tissue. In all clonally related cancer samples, we observed branched evolution, where all metastatic and primary sites shared a common ancestor yet continued to evolve independently. In 53% of cases, we found potentially clinically informative alterations in the brain metastases not detected in the matched primary-tumor sample. In contrast, spatially and temporally separated brain metastasis sites were genetically homogenous. Distal extracranial and regional lymph node metastases were highly divergent from brain metastases. We detected alterations associated with sensitivity to PI3K/AKT/mTOR, CDK, and HER2/EGFR inhibitors in the brain metastases. Genomic analysis of brain metastases provides an opportunity to identify potentially clinically informative alterations not detected in clinically sampled primary tumors, regional lymph nodes, or extracranial metastases.

SIGNIFICANCE

Decisions for individualized therapies in patients with brain metastasis are often made from primary-tumor biopsies. We demonstrate that clinically actionable alterations present in brain metastases are frequently not detected in primary biopsies, suggesting that sequencing of primary biopsies alone may miss a substantial number of opportunities for targeted therapy.

Applying Genomics and Proteomics in Translational Surgical Oncology Research

BACKGROUND

Translational surgical research can range from investigating clinically oriented basic pathophysiological research via the correlation of molecular findings with a patient's clinical outcome to conducting treatment response studies. Herein, the specialist's opinion as a 'surgeon scientist' plays a pivotal role, e.g. in planning experimental setups focusing on the clinically most relevant needs. SUMMARY AND KEY MESSAGES: This review provides a survey of genomic and proteomic technologies and gives an up-to-date overview of those studies applying these technologies exemplarily in colorectal cancer-based biomarker research. Although current results are promising, future validation studies within multicenter networks are required to promote the translation of 'omics' from bench to bedside. In this endeavor, departments of surgery and institutes of experimental surgery together should play a fundamental role. The program for 'clinician scientists' recently launched by the German Research Society (DFG) could be one tool to promote interdisciplinary networks and surgeons' impact on translational research.

Somatic gene copy number alterations in colorectal cancer: new quest for cancer drivers and biomarkers

Colorectal cancer (CRC) results from the accumulation of genetic alterations, and somatic copy number alterations (CNAs) are crucial for the development of CRC. Genome-wide survey of CNAs provides opportunities for identifying cancer driver genes in an unbiased manner. The detection of aberrant CNAs may provide novel markers for the early diagnosis and personalized treatment of CRC. A major challenge in array-based profiling of CNAs is to distinguish the alterations that play causative roles from the random alterations that accumulate during colorectal carcinogenesis. In this view, we systematically discuss the frequent CNAs in CRC, focusing on functional genes that have potential diagnostic, prognostic and therapeutic significance.

Epidermal growth factor receptor immunohistochemistry: new opportunities in metastatic colorectal cancer

The treatment of cancer is becoming more precise, targeting specific oncogenic drivers with targeted molecular therapies. The epidermal growth factor receptor has been found to be over-expressed in a multitude of solid tumours. Immunohistochemistry is widely used in the fields of diagnostic and personalised medicine to localise and visualise disease specific proteins. To date the clinical utility of epidermal growth factor receptor immunohistochemistry in determining monoclonal antibody efficacy has remained somewhat inconclusive. The lack of an agreed reproducible scoring criteria for epidermal growth factor receptor immunohistochemistry has, in various clinical trials yielded conflicting results as to the use of epidermal growth factor receptor immunohistochemistry assay as a companion diagnostic. This has resulted in this test being removed from the licence for the drug panitumumab and not performed in clinical practice for cetuximab. In this review we explore the reasons behind this with a particular emphasis on colorectal cancer, and to suggest a way of resolving the situation through improving the precision of epidermal growth factor receptor immunohistochemistry with quantitative image analysis of digitised images complemented with companion molecular morphological techniques such as in situ hybridisation and section based gene mutation analysis.

Exome Sequencing Reveals AMER1 as a Frequently Mutated Gene in Colorectal Cancer

PURPOSE

Somatic mutations occur at early stages of adenoma and accumulate throughout colorectal cancer progression. The aim of this study was to characterize the mutational landscape of stage II tumors and to search for novel recurrent mutations likely implicated in colorectal cancer tumorigenesis.

EXPERIMENTAL DESIGN

The exomic DNA of 42 stage II, microsatellite-stable colon tumors and their paired mucosae were sequenced. Other molecular data available in the discovery dataset [gene expression, methylation, and copy number variations (CNV)] were used to further characterize these tumors. Additional datasets comprising 553 colorectal cancer samples were used to validate the discovered mutations.

RESULTS

As a result, 4,886 somatic single-nucleotide variants (SNV) were found. Almost all SNVs were private changes, with few mutations shared by more than one tumor, thus revealing tumor-specific mutational landscapes. Nevertheless, these diverse mutations converged into common cellular pathways, such as cell cycle or apoptosis. Among this mutational heterogeneity, variants resulting in early stop codons in the AMER1 (also known as FAM123B or WTX) gene emerged as recurrent mutations in colorectal cancer. Losses of AMER1 by other mechanisms apart from mutations such as methylation and copy number aberrations were also found. Tumors lacking this tumor suppressor gene exhibited a mesenchymal phenotype characterized by inhibition of the canonical Wnt pathway.

CONCLUSIONS

In silico and experimental validation in independent datasets confirmed the existence of functional mutations in AMER1 in approximately 10% of analyzed colorectal cancer tumors. Moreover, these tumors exhibited a characteristic phenotype.