Co-Evolution of Somatic Variation in Primary and Metastatic Colorectal Cancer May Expand Biopsy Indications in the Molecular Era

IntroductionMetastasis is thought to be a clonal event whereby a single cell initiates the development of a new tumor at a distant site. However the degree to which primary and metastatic tumors differ on a molecular level remains unclear. To further evaluate these concepts, we used next generation sequencing (NGS) to assess the molecular composition of paired primary and metastatic colorectal cancer tissue specimens.Methods468 colorectal tumor samples from a large personalized medicine initiative were assessed by targeted gene sequencing of 1,321 individual genes. Eighteen patients produced genomic profiles for 17 paired primary:metastatic (and 2 metastatic:metastatic) specimens.ResultsAn average of 33.3 mutations/tumor were concordant (shared) between matched samples, including common well-known genes (APC, KRAS, TP53). An average of 2.3 mutations/tumor were discordant (unshared) among paired sites. KRAS mutational status was always concordant. The overall concordance rate for mutations was 93.5%; however, nearly all (18/19 (94.7%)) paired tumors showed at least one mutational discordance. Mutations were seen in: TTN, the largest gene (5 discordant pairs), ADAMTS20, APC, MACF1, RASA1, TP53, and WNT2 (2 discordant pairs), SMAD2, SMAD3, SMAD4, FBXW7, and 66 others (1 discordant pair).ConclusionsWhereas primary and metastatic tumors displayed little variance overall, co-evolution produced incremental mutations in both. These results suggest that while biopsy of the primary tumor alone is likely sufficient in the chemotherapy-naïve patient, additional biopsies of primary or metastatic disease may be necessary to precisely tailor therapy following chemotherapy resistance or insensitivity in order to adequately account for tumor evolution.     

Discovery of biclonal origin and a novel oncogene SLC12A5 in colon cancer by single-cell sequencing

Single-cell sequencing is a powerful tool for delineating clonal relationship and identifying key driver genes for personalized cancer management. Here we performed single-cell sequencing analysis of a case of colon cancer. Population genetics analyses identified two independent clones in tumor cell population. The major tumor clone harbored APC and TP53 mutations as early oncogenic events, whereas the minor clone contained preponderant CDC27 and PABPC1 mutations. The absence of APC and TP53 mutations in the minor clone supports that these two clones were derived from two cellular origins. Examination of somatic mutation allele frequency spectra of additional 21 whole-tissue exome-sequenced cases revealed the heterogeneity of clonal origins in colon cancer. Next, we identified a mutated gene SLC12A5 that showed a high frequency of mutation at the single-cell level but exhibited low prevalence at the population level. Functional characterization of mutant SLC12A5 revealed its potential oncogenic effect in colon cancer. Our study provides the first exome-wide evidence at single-cell level supporting that colon cancer could be of a biclonal origin, and suggests that low-prevalence mutations in a cohort may also play important protumorigenic roles at the individual level.     

Intra- and inter-tumor heterogeneity of BRAF(V600E))mutations in primary and metastatic melanoma

The rationale for using small molecule inhibitors of oncogenic proteins as cancer therapies depends, at least in part, on the assumption that metastatic tumors are primarily clonal with respect to mutant oncogene. With the emergence of BRAF^V600E as a therapeutic target, we investigated intra- and inter-tumor heterogeneity in melanoma using detection of the BRAF^V600E mutation as a marker of clonality. BRAF mutant-specific PCR (MS-PCR) and conventional sequencing were performed on 112 tumors from 73 patients, including patients with matched primary and metastatic specimens (n = 18). Nineteen patients had tissues available from multiple metastatic sites. Mutations were detected in 36/112 (32%) melanomas using conventional sequencing, and 85/112 (76%) using MS-PCR. The better sensitivity of the MS-PCR to detect the mutant BRAF^V600E allele was not due to the presence of contaminating normal tissue, suggesting that the tumor was comprised of subclones of differing BRAF genotypes. To determine if tumor subclones were present in individual primary melanomas, we performed laser microdissection and mutation detection via sequencing and BRAF^V600E-specific SNaPshot analysis in 9 cases. Six of these cases demonstrated differing proportions of BRAF^V600Eand BRAF^wild-type cells in distinct microdissected regions within individual tumors. Additional analyses of multiple metastatic samples from individual patients using the highly sensitive MS-PCR without microdissection revealed that 5/19 (26%) patients had metastases that were discordant for the BRAF^V600E mutation. In conclusion, we used highly sensitive BRAF mutation detection methods and observed substantial evidence for heterogeneity of the BRAF^V600E mutation within individual melanoma tumor specimens, and among multiple specimens from individual patients. Given the varied clinical responses of patients to BRAF inhibitor therapy, these data suggest that additional studies to determine possible associations between clinical outcomes and intra- and inter-tumor heterogeneity could prove fruitful.     

Genomic Analyses Reveal Mutational Signatures and Frequently Altered Genes in Esophageal Squamous Cell Carcinoma

Esophageal squamous cell carcinoma (ESCC) is one of the most common cancers worldwide and the fourth most lethal cancer in China. However, although genomic studies have identified some mutations associated with ESCC, we know little of the mutational processes responsible. To identify genome-wide mutational signatures, we performed either whole-genome sequencing (WGS) or whole-exome sequencing (WES) on 104 ESCC individuals and combined our data with those of 88 previously reported samples. An APOBEC-mediated mutational signature in 47% of 192 tumors suggests that APOBEC-catalyzed deamination provides a source of DNA damage in ESCC. Moreover, PIK3CA hotspot mutations (c.1624G>A [p.Glu542Lys] and c.1633G>A [p.Glu545Lys]) were enriched in APOBEC-signature tumors, and no smoking-associated signature was observed in ESCC. In the samples analyzed by WGS, we identified focal (<100 kb) amplifications of CBX4 and CBX8. In our combined cohort, we identified frequent inactivating mutations in AJUBA, ZNF750, and PTCH1 and the chromatin-remodeling genes CREBBP and BAP1, in addition to known mutations. Functional analyses suggest roles for several genes (CBX4, CBX8, AJUBA, and ZNF750) in ESCC. Notably, high activity of hedgehog signaling and the PI3K pathway in approximately 60% of 104 ESCC tumors indicates that therapies targeting these pathways might be particularly promising strategies for ESCC. Collectively, our data provide comprehensive insights into the mutational signatures of ESCC and identify markers for early diagnosis and potential therapeutic targets.     

Exome Sequencing of Cell-Free DNA from Metastatic Cancer Patients Identifies Clinically Actionable Mutations Distinct from Primary Disease

The identification of the molecular drivers of cancer by sequencing is the backbone of precision medicine and the basis of personalized therapy; however, biopsies of primary tumors provide only a snapshot of the evolution of the disease and may miss potential therapeutic targets, especially in the metastatic setting. A liquid biopsy, in the form of cell-free DNA (cfDNA) sequencing, has the potential to capture the inter- and intra-tumoral heterogeneity present in metastatic disease, and, through serial blood draws, track the evolution of the tumor genome.In order to determine the clinical utility of cfDNA sequencing we performed whole-exome sequencing on cfDNA and tumor DNA from two patients with metastatic disease; only minor modifications to our sequencing and analysis pipelines were required for sequencing and mutation calling of cfDNA. The first patient had metastatic sarcoma and 47 of 48 mutations present in the primary tumor were also found in the cell-free DNA. The second patient had metastatic breast cancer and sequencing identified an ESR1 mutation in the cfDNA and metastatic site, but not in the primary tumor. This likely explains tumor progression on Anastrozole. Significant heterogeneity between the primary and metastatic tumors, with cfDNA reflecting the metastases, suggested separation from the primary lesion early in tumor evolution. This is best illustrated by an activating PIK3CA mutation (H1047R) which was clonal in the primary tumor, but completely absent from either the metastasis or cfDNA. Here we show that cfDNA sequencing supplies clinically actionable information with minimal risks compared to metastatic biopsies. This study demonstrates the utility of whole-exome sequencing of cell-free DNA from patients with metastatic disease. cfDNA sequencing identified an ESR1 mutation, potentially explaining a patient’s resistance to aromatase inhibition, and gave insight into how metastatic lesions differ from the primary tumor.     

Genetic evolution of nevus of Ota reveals clonal heterogeneity acquiring BAP1 and TP53 mutations

Melanoma presents molecular alterations based on its anatomical location and exposure to environmental factors. Due to its intrinsic genetic heterogeneity, a simple snapshot of a tumor's genetic alterations does not reflect the tumor clonal complexity or specific gene–gene cooperation. Here, we studied the genetic alterations and clonal evolution of a unique patient with a Nevus of Ota that developed into a recurring uveal‐like dermal melanoma. The Nevus of Ota and ulterior lesions contained GNAQ mutations were c‐KIT positive, and tumors showed an increased RAS pathway activity during progression. Whole‐exome sequencing of these lesions revealed the acquisition of BAP1 and TP53 mutations during tumor evolution, thereby unmasking clonal heterogeneity and allowing the identification of cooperating genes within the same tumor. Our results highlight the importance of studying tumor genetic evolution to identify cooperating mechanisms and delineate effective therapies.     

Biological processes, properties and molecular wiring diagrams of candidate low-penetrance breast cancer susceptibility genes

Background

The process of malignant transformation, progression and metastasis of melanoma is poorly understood. Gene expression profiling of human cancer has allowed for a unique insight into the genes that are involved in these processes. Thus, we have attempted to utilize this approach through the analysis of a series of primary, non-metastatic cutaneous tumors and metastatic melanoma samples.

Methods

We have utilized gene microarray analysis and a variety of molecular techniques to compare 40 metastatic melanoma (MM) samples, composed of 22 bulky, macroscopic (replaced) lymph node metastases, 16 subcutaneous and 2 distant metastases (adrenal and brain), to 42 primary cutaneous cancers, comprised of 16 melanoma, 11 squamous cell, 15 basal cell skin cancers. A Human Genome U133 Plus 2.0 array from Affymetrix, Inc. was utilized for each sample. A variety of statistical software, including the Affymetrix MAS 5.0 analysis software, was utilized to compare primary cancers to metastatic melanomas. Separate analyses were performed to directly compare only primary melanoma to metastatic melanoma samples. The expression levels of putative oncogenes and tumor suppressor genes were analyzed by semi- and real-time quantitative RT-PCR (qPCR) and Western blot analysis was performed on select genes.

Results

We find that primary basal cell carcinomas, squamous cell carcinomas and thin melanomas express dramatically higher levels of many genes, including SPRR1A/B, KRT16/17, CD24, LOR, GATA3, MUC15, and TMPRSS4, than metastatic melanoma. In contrast, the metastatic melanomas express higher levels of genes such as MAGE, GPR19, BCL2A1, MMP14, SOX5, BUB1, RGS20, and more. The transition from non-metastatic expression levels to metastatic expression levels occurs as melanoma tumors thicken. We further evaluated primary melanomas of varying Breslow's tumor thickness to determine that the transition in expression occurs at different thicknesses for different genes suggesting that the "transition zone" represents a critical time for the emergence of the metastatic phenotype. Several putative tumor oncogenes (SPP-1, MITF, CITED-1, GDF-15, c-Met, HOX loci) and suppressor genes (PITX-1, CST-6, PDGFRL, DSC-3, POU2F3, CLCA2, ST7L), were identified and validated by quantitative PCR as changing expression during this transition period. These are strong candidates for genes involved in the progression or suppression of the metastatic phenotype.

Conclusion

The gene expression profiling of primary, non-metastatic cutaneous tumors and metastatic melanoma has resulted in the identification of several genes that may be centrally involved in the progression and metastatic potential of melanoma. This has very important implications as we continue to develop an improved understanding of the metastatic process, allowing us to identify specific genes for prognostic markers and possibly for targeted therapeutic approaches.     

Comparative sequencing analysis reveals high genomic concordance between matched primary and metastatic colorectal cancer lesions

Background

Colorectal cancer is the second leading cause of cancer death in the United States, with over 50,000 deaths estimated in 2014. Molecular profiling for somatic mutations that predict absence of response to anti-EGFR therapy has become standard practice in the treatment of metastatic colorectal cancer; however, the quantity and type of tissue available for testing is frequently limited. Further, the degree to which the primary tumor is a faithful representation of metastatic disease has been questioned. As next-generation sequencing technology becomes more widely available for clinical use and additional molecularly targeted agents are considered as treatment options in colorectal cancer, it is important to characterize the extent of tumor heterogeneity between primary and metastatic tumors.

Results

We performed deep coverage, targeted next-generation sequencing of 230 key cancer-associated genes for 69 matched primary and metastatic tumors and normal tissue. Mutation profiles were 100% concordant for KRAS, NRAS, and BRAF, and were highly concordant for recurrent alterations in colorectal cancer. Additionally, whole genome sequencing of four patient trios did not reveal any additional site-specific targetable alterations.

Conclusions

Colorectal cancer primary tumors and metastases exhibit high genomic concordance. As current clinical practices in colorectal cancer revolve around KRAS, NRAS, and BRAF mutation status, diagnostic sequencing of either primary or metastatic tissue as available is acceptable for most patients. Additionally, consistency between targeted sequencing and whole genome sequencing results suggests that targeted sequencing may be a suitable strategy for clinical diagnostic applications.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0454-7) contains supplementary material, which is available to authorized users.     

Whole Genome Sequencing versus Traditional Genotyping for Investigation of a Mycobacterium tuberculosis Outbreak: A Longitudinal Molecular Epidemiological Study

Background

Understanding Mycobacterium tuberculosis (Mtb) transmission is essential to guide efficient tuberculosis control strategies. Traditional strain typing lacks sufficient discriminatory power to resolve large outbreaks. Here, we tested the potential of using next generation genome sequencing for identification of outbreak-related transmission chains.

Methods and Findings

During long-term (1997 to 2010) prospective population-based molecular epidemiological surveillance comprising a total of 2,301 patients, we identified a large outbreak caused by an Mtb strain of the Haarlem lineage. The main performance outcome measure of whole genome sequencing (WGS) analyses was the degree of correlation of the WGS analyses with contact tracing data and the spatio-temporal distribution of the outbreak cases. WGS analyses of the 86 isolates revealed 85 single nucleotide polymorphisms (SNPs), subdividing the outbreak into seven genome clusters (two to 24 isolates each), plus 36 unique SNP profiles. WGS results showed that the first outbreak isolates detected in 1997 were falsely clustered by classical genotyping. In 1998, one clone (termed “Hamburg clone”) started expanding, apparently independently from differences in the social environment of early cases. Genome-based clustering patterns were in better accordance with contact tracing data and the geographical distribution of the cases than clustering patterns based on classical genotyping. A maximum of three SNPs were identified in eight confirmed human-to-human transmission chains, involving 31 patients. We estimated the Mtb genome evolutionary rate at 0.4 mutations per genome per year. This rate suggests that Mtb grows in its natural host with a doubling time of approximately 22 h (400 generations per year). Based on the genome variation discovered, emergence of the Hamburg clone was dated back to a period between 1993 and 1997, hence shortly before the discovery of the outbreak through epidemiological surveillance.

Conclusions

Our findings suggest that WGS is superior to conventional genotyping for Mtb pathogen tracing and investigating micro-epidemics. WGS provides a measure of Mtb genome evolution over time in its natural host context. Please see later in the article for the Editors'' Summary     

Distinct mutational backgrounds and clonal architectures implicated prognostic discrepancies in small-cell carcinomas of the esophagus and lung

Small-cell carcinoma of the esophagus (SCCE) is a rare and aggressive cancer. Although several consistent genomic changes were observed previously between SCCE and small-cell lung cancer (SCLC), detailed mutational landscapes revealing discrepancies in genetic underpinnings of tumorigenesis between these two cancers are scarce, and little attention has been paid to answer whether these genetic alterations were related to the prognosis. Herein by performing whole-exome sequencing of 48 SCCE and 64 SCLC tumor samples, respectively we have shown that the number of driver mutations in SCCE was significantly lower than in SCLC (p = 0.0042). In SCCE, 46% of recurrent driver mutations were clonal, which occurred at an early stage during tumorigenesis, while 16 driver mutations were found clonal in SCLC. NOTCH1/3, PIK3CA, and ATM were specifically clonal in SCCE, while TP53 was clonal in SCLC. The total number of clonal mutations differed between two cancers and presented lower in SCCE compared to SCLC (p = 0.0036). Moreover, overall survival (OS) was shorter in patients with higher numbers of clonal mutations for both cancers. In summary, SCCE showed distinct mutational background and clonal architecture compared with SCLC. Organ-specific clonal events revealed different molecular mechanisms underlying tumorigenesis, tumor development, patients’ prognosis, and possible variations in therapeutic outcomes to candidate treatments.Subject terms: Oesophageal cancer, Small-cell lung cancer     

IDH1 mutation identified in one human melanoma metastasis, but not correlated with metastases to the brain

Isocitrate dehydrogenase 1 (IDH1) and isocitrate dehydrogenase 2 (IDH2) are enzymes which convert isocitrate to α-ketoglutarate while reducing nicotinamide adenine dinucleotide phosphate (NADP + to NADPH). IDH1/2 were recently identified as mutated in a large percentage of progressive gliomas. These mutations occur at IDH1^R132 or the homologous IDH2^R172. Melanomas share some genetic features with IDH1/2-mutated gliomas, such as frequent TP53 mutation. We sought to test whether melanoma is associated with IDH1/2 mutations. Seventy-eight human melanoma samples were analyzed for IDH1^R132 and IDH2^R172 mutation status. A somatic, heterozygous IDH1 c.C394T (p.R132C) mutation was identified in one human melanoma metastasis to the lung. Having identified this mutation in one metastasis, we sought to test the hypothesis that certain selective pressures in the brain environment may specifically favor the cell growth or survival of tumor cells with mutations in IDH1/2, regardless of primary tumor site. To address this, we analyzed IDH1^R132 and IDH2^R172 mutation status 53 metastatic brain tumors, including nine melanoma metastases. Results revealed no mutations in any samples. This lack of mutations would suggest that mutations in IDH1^R132 or IDH2^R172 may be necessary for the formation of tumors in a cell-lineage dependent manner, with a particularly strong selective pressure for mutations in progressive gliomas; this also suggests the lack of a particular selective pressure for growth in brain tissue in general. Studies on the cell-lineages of tumors with IDH1/2 mutations may help clarify the role of these mutations in the development of brain tumors.     

Chromothripsis is a common mechanism driving genomic rearrangements in primary and metastatic colorectal cancer

Background

Structural rearrangements form a major class of somatic variation in cancer genomes. Local chromosome shattering, termed chromothripsis, is a mechanism proposed to be the cause of clustered chromosomal rearrangements and was recently described to occur in a small percentage of tumors. The significance of these clusters for tumor development or metastatic spread is largely unclear.

Results

We used genome-wide long mate-pair sequencing and SNP array profiling to reveal that chromothripsis is a widespread phenomenon in primary colorectal cancer and metastases. We find large and small chromothripsis events in nearly every colorectal tumor sample and show that several breakpoints of chromothripsis clusters and isolated rearrangements affect cancer genes, including NOTCH2, EXO1 and MLL3. We complemented the structural variation studies by sequencing the coding regions of a cancer exome in all colorectal tumor samples and found somatic mutations in 24 genes, including APC, KRAS, SMAD4 and PIK3CA. A pairwise comparison of somatic variations in primary and metastatic samples indicated that many chromothripsis clusters, isolated rearrangements and point mutations are exclusively present in either the primary tumor or the metastasis and may affect cancer genes in a lesion-specific manner.

Conclusions

We conclude that chromothripsis is a prevalent mechanism driving structural rearrangements in colorectal cancer and show that a complex interplay between point mutations, simple copy number changes and chromothripsis events drive colorectal tumor development and metastasis.     

TrAp: a tree approach for fingerprinting subclonal tumor composition

Revealing the clonal composition of a single tumor is essential for identifying cell subpopulations with metastatic potential in primary tumors or with resistance to therapies in metastatic tumors. Sequencing technologies provide only an overview of the aggregate of numerous cells. Computational approaches to de-mix a collective signal composed of the aberrations of a mixed cell population of a tumor sample into its individual components are not available. We propose an evolutionary framework for deconvolving data from a single genome-wide experiment to infer the composition, abundance and evolutionary paths of the underlying cell subpopulations of a tumor. We have developed an algorithm (TrAp) for solving this mixture problem. In silico analyses show that TrAp correctly deconvolves mixed subpopulations when the number of subpopulations and the measurement errors are moderate. We demonstrate the applicability of the method using tumor karyotypes and somatic hypermutation data sets. We applied TrAp to Exome-Seq experiment of a renal cell carcinoma tumor sample and compared the mutational profile of the inferred subpopulations to the mutational profiles of single cells of the same tumor. Finally, we deconvolve sequencing data from eight acute myeloid leukemia patients and three distinct metastases of one melanoma patient to exhibit the evolutionary relationships of their subpopulations.     

Mutational spectra of human cancer

The purpose of this review is to summarize the evidence that can be used to reconstruct the etiology of human cancers from mutations found in tumors. Mutational spectra of the tumor suppressor gene p53 (TP53) are tumor specific. In several cases, these mutational spectra can be linked to exogenous carcinogens, most notably for sunlight-associated skin cancers, tobacco-associated lung cancers, and aristolochic acid-related urothelial tumors. In the TP53 gene, methylated CpG dinucleotides are sequences selectively targeted by endogenous and exogenous mutagenic processes. Recent high-throughput sequencing efforts analyzing a large number of genes in cancer genomes have so far, for the most part, produced mutational spectra similar to those in TP53 but have unveiled a previously unrecognized common G to C transversion mutation signature at GpA dinucleotides in breast cancers and several other cancers. Unraveling the origin of these G to C mutations will be of importance for understanding cancer etiology.     

Curcumin Intake Affects miRNA Signature in Murine Melanoma with mmu-miR-205-5p Most Significantly Altered

Melanoma is the most aggressive form of skin cancer with estimated 48,000 deaths per year worldwide. The polyphenol curcumin derived from the plant Curcuma longa is well known for its anti-inflammatory and anti-cancerogenic properties. Accordingly, dietary intake of this compound may be suitable for melanoma prevention. However, how this compound affects basic cellular mechanisms in developing melanoma still remains elusive. Therefore, the aim of this study was to investigate for the first time the impact of oral curcumin administration on the miRNA signature of engrafting melanoma. For this purpose, the effects of a 4% curcumin diet were tested on melanoma, which were established by injection of murine B78H1 cells in the flank of C57BL/6 mice. Curcumin diet or standard chow (control) was administered two weeks prior to injection of tumor cells until termination of the experiment. High throughput chip-based array analysis was deployed to detect alterations in the miRNA signature of the tumors. Curcumin treatment significantly reduced the growth of the flank tumors. Furthermore the miRNA expression signature in tumors was substantially altered by curcumin intake with mmu-miR-205-5p over 100 times higher expressed when compared to controls. The expression levels of identified key miRNAs in the tumor samples were confirmed by quantitative real-time polymerase chain reaction (qRT-PCR). A comparable expression pattern of these miRNAs was also detected in other curcumin-treated melanoma cell lines under in vitro conditions. Putative targets of curcumin-induced up-regulated miRNAs were enriched in ‘o-glycan biosynthesis’, ‘endoplasmatic reticulum protein processing’ and different cancer-related pathways. Western Blot analyses revealed that of these targets anti-apoptotic B-cell CLL/lymphoma 2 (Bcl-2) and proliferating cell nuclear antigen (PCNA) were significantly down-regulated in curcumin-treated tumors. These findings demonstrate a profound alteration of the miRNA expression signature in engrafting curcumin-treated melanoma with mmu-miR-205-5p being up-regulated most significantly.     

Tumor Homogeneity between Primary and Metastatic Sites for BRAF Status in Metastatic Melanoma Determined by Immunohistochemical and Molecular Testing

BRAF inhibitors have demonstrated improvement of overall survival in patients with metastatic melanoma and BRAF^V600 mutations. In order to evaluate BRAF tumor heterogeneity between primary and metastatic site, we have evaluated the performance of immunohistochemistry (IHC) with an anti-BRAF^V600E antibody in both localization by comparison with high resolution melting analysis followed by Sanger sequencing in a parallel blinded study. A total of 230 samples distributed as primary melanoma (n = 88) and different types of metastatic samples (n = 142) were studied in 99 patients with advanced or metastatic melanoma (stage III or IV). The prevalence of each BRAF mutation was c.1799T>A, BRAF^V600E (45.2%), c.1799_1800TG>AA, BRAF^V600E2 (3.0%), c.1798_1799GT>AA, BRAF^V600K (3.0%), c.1801 A>G, BRAF^K601E (1.3%), c.1789_1790CT>TC, BRAF^L597S (0.4%), c.1780G>A, BRAF^D594N (0.9%) respectively. IHC was positive in 109/112 samples harboring BRAF^V600E/E2 mutations and negative in other cases. The cytoplasmic staining was either strongly positive in tumor cells of BRAF^V600E mutated cases. It appeared strong brown, different from the vesicular grey cytoplasmic pigmentation of melanophages. Concordance between the two techniques was 96.4%. Sensitivity of IHC for detecting the BRAF^V600E/E2 mutations was 97.3%, while specificity was 100%. Both our IHC and molecular study demonstrated homogeneity between primary and metastatic sites for BRAF status in melanoma. This study also provides evidence that IHC may be a cost-effective first-line method for BRAF^V600E detection. Thereafter, molecular techniques should be used in negative, ambiguous or non-contributive cases.     

Colorectal Cancer Genetic Heterogeneity Delineated by Multi-Region Sequencing

Intratumor heterogeneity (ITH) leads to an underestimation of the mutational landscape portrayed by a single needle biopsy and consequently affects treatment precision. The extent of colorectal cancer (CRC) genetic ITH is not well understood in Chinese patients. Thus, we conducted deep sequencing by using the OncoGxOne^™ Plus panel, targeting 333 cancer-specific genes in multi-region biopsies of primary and liver metastatic tumors from three Chinese CRC patients. We determined that the extent of ITH varied among the three cases. On average, 65% of all the mutations detected were common within individual tumors. KMT2C aberrations and the NCOR1 mutation were the only ubiquitous events. Subsequent phylogenetic analysis showed that the tumors evolved in a branched manner. Comparison of the primary and metastatic tumors revealed that PPP2R1A (E370X), SETD2 (I1608V), SMAD4 (G382T), and AR splicing site mutations may be specific to liver metastatic cancer. These mutations might contribute to the initiation and progression of distant metastasis. Collectively, our analysis identified a substantial level of genetic ITH in CRC, which should be considered for personalized therapeutic strategies.     

Evolution under Spatially Heterogeneous Selection in Solid Tumors

Spatial genetic and phenotypic diversity within solid tumors has been well documented. Nevertheless, how this heterogeneity affects temporal dynamics of tumorigenesis has not been rigorously examined because solid tumors do not evolve as the standard population genetic model due to the spatial constraint. We therefore, propose a neutral spatial (NS) model whereby the mutation accumulation increases toward the periphery; the genealogical relationship is spatially determined and the selection efficacy is blunted (due to kin competition). In this model, neutral mutations are accrued and spatially distributed in manners different from those of advantageous mutations. Importantly, the distinctions could be blurred in the conventional model. To test the NS model, we performed a three-dimensional multiple microsampling of two hepatocellular carcinomas. Whole-genome sequencing (WGS) revealed a 2-fold increase in mutations going from the center to the periphery. The operation of natural selection can then be tested by examining the spatially determined clonal relationships and the clonal sizes. Due to limited migration, only the expansion of highly advantageous clones can sweep through a large part of the tumor to reveal the selective advantages. Hence, even multiregional sampling can only reveal a fraction of fitness differences in solid tumors. Our results suggest that the NS patterns are crucial for testing the influence of natural selection during tumorigenesis, especially for small solid tumors.