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.     

Absolute quantification of somatic DNA alterations in human cancer

We describe a computational method that infers tumor purity and malignant cell ploidy directly from analysis of somatic DNA alterations. The method, named ABSOLUTE, can detect subclonal heterogeneity and somatic homozygosity, and it can calculate statistical sensitivity for detection of specific aberrations. We used ABSOLUTE to analyze exome sequencing data from 214 ovarian carcinoma tumor-normal pairs. This analysis identified both pervasive subclonal somatic point-mutations and a small subset of predominantly clonal and homozygous mutations, which were overrepresented in the tumor suppressor genes TP53 and NF1 and in a candidate tumor suppressor gene CDK12. We also used ABSOLUTE to infer absolute allelic copy-number profiles from 3,155 diverse cancer specimens, revealing that genome-doubling events are common in human cancer, likely occur in cells that are already aneuploid, and influence pathways of tumor progression (for example, with recessive inactivation of NF1 being less common after genome doubling). ABSOLUTE will facilitate the design of clinical sequencing studies and studies of cancer genome evolution and intra-tumor heterogeneity.     

Mutational Profiling Can Establish Clonal or Independent Origin in Synchronous Bilateral Breast and Other Tumors

Background

Synchronous tumors can be independent primary tumors or a primary-metastatic (clonal) pair, which may have clinical implications. Mutational profiling of tumor DNA is increasingly common in the clinic. We investigated whether mutational profiling can distinguish independent from clonal tumors in breast and other cancers, using a carefully defined test based on the Clonal Likelihood Score (CLS = 100 x # shared high confidence (HC) mutations/ # total HC mutations).

Methods

Statistical properties of a formal test using the CLS were investigated. A high CLS is evidence in favor of clonality; the test is implemented as a one-sided binomial test of proportions. Test parameters were empirically determined using 16,422 independent breast tumor pairs and 15 primary-metastatic tumor pairs from 10 cancer types using The Cancer Genome Atlas.

Results

We validated performance of the test with its established parameters, using five published data sets comprising 15,758 known independent tumor pairs (maximum CLS = 4.1%, minimum p-value = 0.48) and 283 known tumor clonal pairs (minimum CLS 13%, maximum p-value <0.01), across renal cell, testicular, and colorectal cancer. The CLS test correctly classified all validation samples but one, which it appears may have been incorrectly classified in the published data. As proof-of-concept we then applied the CLS test to two new cases of invasive synchronous bilateral breast cancer at our institution, each with one hormone receptor positive (ER+/PR+/HER2-) lobular and one triple negative ductal carcinoma. High confidence mutations were identified by exome sequencing and results were validated using deep targeted sequencing. The first tumor pair had CLS of 81% (p-value < 10–15), supporting clonality. In the second pair, no common mutations of 184 variants were validated (p-value >0.99), supporting independence. A plausible molecular mechanism for the shift from hormone receptor positive to triple negative was identified in the clonal pair.

Conclusion

We have developed the statistical properties of a carefully defined Clonal Likelihood Score test from mutational profiling of tumor DNA. Under identified conditions, the test appears to reliably distinguish between synchronous tumors of clonal and of independent origin in several cancer types. This approach may have scientific and clinical utility.     

SciClone: Inferring Clonal Architecture and Tracking the Spatial and Temporal Patterns of Tumor Evolution

The sensitivity of massively-parallel sequencing has confirmed that most cancers are oligoclonal, with subpopulations of neoplastic cells harboring distinct mutations. A fine resolution view of this clonal architecture provides insight into tumor heterogeneity, evolution, and treatment response, all of which may have clinical implications. Single tumor analysis already contributes to understanding these phenomena. However, cryptic subclones are frequently revealed by additional patient samples (e.g., collected at relapse or following treatment), indicating that accurately characterizing a tumor requires analyzing multiple samples from the same patient. To address this need, we present SciClone, a computational method that identifies the number and genetic composition of subclones by analyzing the variant allele frequencies of somatic mutations. We use it to detect subclones in acute myeloid leukemia and breast cancer samples that, though present at disease onset, are not evident from a single primary tumor sample. By doing so, we can track tumor evolution and identify the spatial origins of cells resisting therapy.     

VDJ-Seq: Deep Sequencing Analysis of Rearranged Immunoglobulin Heavy Chain Gene to Reveal Clonal Evolution Patterns of B Cell Lymphoma

Understanding tumor clonality is critical to understanding the mechanisms involved in tumorigenesis and disease progression. In addition, understanding the clonal composition changes that occur within a tumor in response to certain micro-environment or treatments may lead to the design of more sophisticated and effective approaches to eradicate tumor cells. However, tracking tumor clonal sub-populations has been challenging due to the lack of distinguishable markers. To address this problem, a VDJ-seq protocol was created to trace the clonal evolution patterns of diffuse large B cell lymphoma (DLBCL) relapse by exploiting VDJ recombination and somatic hypermutation (SHM), two unique features of B cell lymphomas.In this protocol, Next-Generation sequencing (NGS) libraries with indexing potential were constructed from amplified rearranged immunoglobulin heavy chain (IgH) VDJ region from pairs of primary diagnosis and relapse DLBCL samples. On average more than half million VDJ sequences per sample were obtained after sequencing, which contain both VDJ rearrangement and SHM information. In addition, customized bioinformatics pipelines were developed to fully utilize sequence information for the characterization of IgH-VDJ repertoire within these samples. Furthermore, the pipeline allows the reconstruction and comparison of the clonal architecture of individual tumors, which enables the examination of the clonal heterogeneity within the diagnosis tumors and deduction of clonal evolution patterns between diagnosis and relapse tumor pairs. When applying this analysis to several diagnosis-relapse pairs, we uncovered key evidence that multiple distinctive tumor evolutionary patterns could lead to DLBCL relapse. Additionally, this approach can be expanded into other clinical aspects, such as identification of minimal residual disease, monitoring relapse progress and treatment response, and investigation of immune repertoires in non-lymphoma contexts.     

Single-cell exome sequencing and monoclonal evolution of a JAK2-negative myeloproliferative neoplasm

Tumor heterogeneity presents a challenge for inferring clonal evolution and driver gene identification. Here, we describe a method for analyzing the cancer genome at a single-cell nucleotide level. To perform our analyses, we first devised and validated a high-throughput whole-genome single-cell sequencing method using two lymphoblastoid cell line single cells. We then carried out whole-exome single-cell sequencing of 90 cells from a JAK2-negative myeloproliferative neoplasm patient. The sequencing data from 58 cells passed our quality control criteria, and these data indicated that this neoplasm represented a monoclonal evolution. We further identified essential thrombocythemia (ET)-related candidate mutations such as SESN2 and NTRK1, which may be involved in neoplasm progression. This pilot study allowed the initial characterization of the disease-related genetic architecture at the single-cell nucleotide level. Further, we established a single-cell sequencing method that opens the way for detailed analyses of a variety of tumor types, including those with high genetic complex between patients.     

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.     

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.     

Implications of genetic heterogeneity in cancer

DNA sequencing studies have established that many cancers contain tens of thousands of clonal mutations throughout their genomes, which is difficult to reconcile with the very low rate of mutation in normal human cells. This observation provides strong evidence for the mutator phenotype hypothesis, which proposes that a genome-wide elevation in the spontaneous mutation rate is an early step in carcinogenesis. An elevated mutation rate implies that cancers undergo continuous evolution, generating multiple subpopulations of cells that differ from one another in DNA sequence. The extensive heterogeneity in DNA sequence and continual tumor evolution that would occur in the context of a mutator phenotype have important implications for cancer diagnosis and therapy.     

Hormone resistance and neuroendocrine differentiation due to accumulation of genetic lesions during clonal evolution of prostate cancer

Progression of malignant tumors is largely due to clonal evolution of the primary tumor, clones acquiring different sets of molecular genetic lesions. Lesions can confer a selective advantage in proliferation rate or metastasis on the tumor cell population, especially if developing resistance to anticancer therapy. Prostate cancer (PCa) provides an illustrative example of clinically significant clonal evolution. The review considers the genetic alterations that occur in primary PCa and the mechanism whereby hormone-refractory PCa develops on hormone therapy, including mutations and alternative splicing of the androgen receptor gene (AR) and intratumoral androgen synthesis. Certain molecular genetic lesions determine resistance to new generation inhibitors (AR mutations that block the antagonist effect or allow other hormones to activate the receptor) or lead to neuroendocrine differentiation (repression of the AR signaling pathway, TP53 mutations, and amplification of the AURKA or MYCN oncogene). Multistep therapy based on the data about somatic mutations associated with progression and metastasis of the primary tumor can be expected to significantly improve the survival of patients with advanced PCa in the nearest future.     

Inferring Clonal Composition from Multiple Sections of a Breast Cancer

Cancers arise from successive rounds of mutation and selection, generating clonal populations that vary in size, mutational content and drug responsiveness. Ascertaining the clonal composition of a tumor is therefore important both for prognosis and therapy. Mutation counts and frequencies resulting from next-generation sequencing (NGS) potentially reflect a tumor''s clonal composition; however, deconvolving NGS data to infer a tumor''s clonal structure presents a major challenge. We propose a generative model for NGS data derived from multiple subsections of a single tumor, and we describe an expectation-maximization procedure for estimating the clonal genotypes and relative frequencies using this model. We demonstrate, via simulation, the validity of the approach, and then use our algorithm to assess the clonal composition of a primary breast cancer and associated metastatic lymph node. After dividing the tumor into subsections, we perform exome sequencing for each subsection to assess mutational content, followed by deep sequencing to precisely count normal and variant alleles within each subsection. By quantifying the frequencies of 17 somatic variants, we demonstrate that our algorithm predicts clonal relationships that are both phylogenetically and spatially plausible. Applying this method to larger numbers of tumors should cast light on the clonal evolution of cancers in space and time.     

Deep sequencing reveals clonal evolution patterns and mutation events associated with relapse in B-cell lymphomas

Background

Molecular mechanisms associated with frequent relapse of diffuse large B-cell lymphoma (DLBCL) are poorly defined. It is especially unclear how primary tumor clonal heterogeneity contributes to relapse. Here, we explore unique features of B-cell lymphomas - VDJ recombination and somatic hypermutation - to address this question.

Results

We performed high-throughput sequencing of rearranged VDJ junctions in 14 pairs of matched diagnosis-relapse tumors, among which 7 pairs were further characterized by exome sequencing. We identify two distinctive modes of clonal evolution of DLBCL relapse: an early-divergent mode in which clonally related diagnosis and relapse tumors diverged early and developed in parallel; and a late-divergent mode in which relapse tumors developed directly from diagnosis tumors with minor divergence. By examining mutation patterns in the context of phylogenetic information provided by VDJ junctions, we identified mutations in epigenetic modifiers such as KMT2D as potential early driving events in lymphomagenesis and immune escape alterations as relapse-associated events.

Conclusions

Altogether, our study for the first time provides important evidence that DLBCL relapse may result from multiple, distinct tumor evolutionary mechanisms, providing rationale for therapies for each mechanism. Moreover, this study highlights the urgent need to understand the driving roles of epigenetic modifier mutations in lymphomagenesis, and immune surveillance factor genetic lesions in relapse.

Electronic supplementary material

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

scDPN for High-throughput Single-cell CNV Detection to Uncover Clonal Evolution During HCC Recurrence

Single-cell genomics provides substantial resources for dissecting cellular heterogeneity and cancer evolution. Unfortunately, classical DNA amplification-based methods have low throughput and introduce coverage bias during sample preamplification. We developed a single-cell DNA library preparation method without preamplification in nanolitre scale (scDPN) to address these issues. The method achieved a throughput of up to 1800 cells per run for copy number variation (CNV) detection. Also, our approach demonstrated a lower level of amplification bias and noise than the multiple displacement amplification (MDA) method and showed high sensitivity and accuracy for cell line and tumor tissue evaluation. We used this approach to profile the tumor clones in paired primary and relapsed tumor samples of hepato-cellular carcinoma (HCC). We identified three clonal subpopulations with a multitude of aneuploid alterations across the genome. Furthermore, we observed that a minor clone of the primary tumor containing additional alterations in chro-mosomes 1q, 10q, and 14q developed into the dominant clone in the recurrent tumor, indicating clonal selection during recurrence in HCC. Overall, this approach provides a comprehensive and scalable solution to understand genome hetero-geneity and evolution.     

Novel molecular and computational methods improve the accuracy of insertion site analysis in Sleeping Beauty-induced tumors

The recent development of the Sleeping Beauty (SB) system has led to the development of novel mouse models of cancer. Unlike spontaneous models, SB causes cancer through the action of mutagenic transposons that are mobilized in the genomes of somatic cells to induce mutations in cancer genes. While previous methods have successfully identified many transposon-tagged mutations in SB-induced tumors, limitations in DNA sequencing technology have prevented a comprehensive analysis of large tumor cohorts. Here we describe a novel method for producing genetic profiles of SB-induced tumors using Illumina sequencing. This method has dramatically increased the number of transposon-induced mutations identified in each tumor sample to reveal a level of genetic complexity much greater than previously appreciated. In addition, Illumina sequencing has allowed us to more precisely determine the depth of sequencing required to obtain a reproducible signature of transposon-induced mutations within tumor samples. The use of Illumina sequencing to characterize SB-induced tumors should significantly reduce sampling error that undoubtedly occurs using previous sequencing methods. As a consequence, the improved accuracy and precision provided by this method will allow candidate cancer genes to be identified with greater confidence. Overall, this method will facilitate ongoing efforts to decipher the genetic complexity of the human cancer genome by providing more accurate comparative information from Sleeping Beauty models of cancer.     

Characterizing mutational heterogeneity in a glioblastoma patient with double recurrence

Human cancers are driven by the acquisition of somatic mutations. Separating the driving mutations from those that are random consequences of general genomic instability remains a challenge. New sequencing technology makes it possible to detect mutations that are present in only a minority of cells in a heterogeneous tumor population. We sought to leverage the power of ultra-deep sequencing to study various levels of tumor heterogeneity in the serial recurrences of a single glioblastoma multiforme patient. Our goal was to gain insight into the temporal succession of DNA base-level lesions by querying intra- and inter-tumoral cell populations in the same patient over time. We performed targeted "next-generation" sequencing on seven samples from the same patient: two foci within the primary tumor, two foci within an initial recurrence, two foci within a second recurrence, and normal blood. Our study reveals multiple levels of mutational heterogeneity. We found variable frequencies of specific EGFR, PIK3CA, PTEN, and TP53 base substitutions within individual tumor regions and across distinct regions within the same tumor. In addition, specific mutations emerge and disappear along the temporal spectrum from tumor at the time of diagnosis to second recurrence, demonstrating evolution during tumor progression. Our results shed light on the spatial and temporal complexity of brain tumors. As sequencing costs continue to decline and deep sequencing technology eventually moves into the clinic, this approach may provide guidance for treatment choices as we embark on the path to personalized cancer medicine.     

Mitochondrial genome variations in advanced stage breast cancer: A case–control study

Entire mitochondrial DNA (mtDNA) sequencing was carried out in 101 primary breast cancer patients and 90 controls of south Indian origin. We identified 69 novel mutations in breast cancer patients and 637 reported polymorphisms in patients and/or controls. PolyPhen-2 analysis predicted 5 out of 14 novel missense mutations as ‘probably damaging variants’. Haplogrouping analysis identified a significant association between haplogroup M5 and breast cancer risk. Microsatellite instability and tumor specific large scale mtDNA deletions were not observed in tumor tissues from the patients. In conclusion, mtDNA mutations and haplogroups may constitute an inheritable risk factor for pathogenesis of breast cancer.