Clonal evolution is a key feature of cancer progression and relapse. We studied intratumoral heterogeneity in 149 chronic lymphocytic leukemia (CLL) cases by integrating whole-exome sequence and copy ...number to measure the fraction of cancer cells harboring each somatic mutation. We identified driver mutations as predominantly clonal (e.g., MYD88, trisomy 12, and del(13q)) or subclonal (e.g., SF3B1 and TP53), corresponding to earlier and later events in CLL evolution. We sampled leukemia cells from 18 patients at two time points. Ten of twelve CLL cases treated with chemotherapy (but only one of six without treatment) underwent clonal evolution, predominantly involving subclones with driver mutations (e.g., SF3B1 and TP53) that expanded over time. Furthermore, presence of a subclonal driver mutation was an independent risk factor for rapid disease progression. Our study thus uncovers patterns of clonal evolution in CLL, providing insights into its stepwise transformation, and links the presence of subclones with adverse clinical outcomes.
Display omitted
► Whole-exome analysis of clonal heterogeneity in 149 chronic lymphocytic leukemias ► Earlier and later mutations in the temporal evolution of CLL are identified ► Clonal evolution is commonly seen with treatment, typically in a branched pattern ► A subclonal driver in a pretreatment sample is associated with adverse outcome
The intratumoral heterogeneity in 149 chronic lymphocytic leukemia (CLL) cases was evaluated by whole-exome sequencing. The evolutionary patterns of distinct clones enabled a temporal ordering of mutations in CLL, revealed the association of clonal evolution with chemotherapy, and linked the presence of subclonal driver mutations with adverse clinical outcomes.
Detection of somatic point substitutions is a key step in characterizing the cancer genome. However, existing methods typically miss low-allelic-fraction mutations that occur in only a subset of the ...sequenced cells owing to either tumor heterogeneity or contamination by normal cells. Here we present MuTect, a method that applies a Bayesian classifier to detect somatic mutations with very low allele fractions, requiring only a few supporting reads, followed by carefully tuned filters that ensure high specificity. We also describe benchmarking approaches that use real, rather than simulated, sequencing data to evaluate the sensitivity and specificity as a function of sequencing depth, base quality and allelic fraction. Compared with other methods, MuTect has higher sensitivity with similar specificity, especially for mutations with allelic fractions as low as 0.1 and below, making MuTect particularly useful for studying cancer subclones and their evolution in standard exome and genome sequencing data.
Small cell lung carcinoma (SCLC) is a highly lethal, smoking-associated cancer with few known targetable genetic alterations. Using genome sequencing, we characterized the somatic evolution of a ...genetically engineered mouse model (GEMM) of SCLC initiated by loss of Trp53 and Rb1. We identified alterations in DNA copy number and complex genomic rearrangements and demonstrated a low somatic point mutation frequency in the absence of tobacco mutagens. Alterations targeting the tumor suppressor Pten occurred in the majority of murine SCLC studied, and engineered Pten deletion accelerated murine SCLC and abrogated loss of Chr19 in Trp53; Rb1; Pten compound mutant tumors. Finally, we found evidence for polyclonal and sequential metastatic spread of murine SCLC by comparative sequencing of families of related primary tumors and metastases. We propose a temporal model of SCLC tumorigenesis with implications for human SCLC therapeutics and the nature of cancer-genome evolution in GEMMs.
Display omitted
•Murine SCLCs acquire few point mutations in the absence of tobacco mutagens•Pten is recurrently mutated, and engineered deletion accelerates tumor progression•Mycl1 amplifications play an early, central role in mSCLC tumorigenesis•Analysis of related primary and metastatic mSCLC suggests complex clonal evolution
Comprehensive genomic analysis in a mouse model of small-cell lung carcinoma delineates metastatic progression in this tumor type, showing that selection for specific cancer mutations can drive large genomic rearrangements and that distant metastases likely arise from a common metastatic seeding step involving the lymph nodes.
We performed massively parallel sequencing of paired tumor/normal samples from 203 multiple myeloma (MM) patients and identified significantly mutated genes and copy number alterations and discovered ...putative tumor suppressor genes by determining homozygous deletions and loss of heterozygosity. We observed frequent mutations in KRAS (particularly in previously treated patients), NRAS, BRAF, FAM46C, TP53, and DIS3 (particularly in nonhyperdiploid MM). Mutations were often present in subclonal populations, and multiple mutations within the same pathway (e.g., KRAS, NRAS, and BRAF) were observed in the same patient. In vitro modeling predicts only partial treatment efficacy of targeting subclonal mutations, and even growth promotion of nonmutated subclones in some cases. These results emphasize the importance of heterogeneity analysis for treatment decisions.
•Deep sequencing reveals significant genetic events in multiple myeloma•Intratumor genetic heterogeneity is common in multiple myeloma•Recurrent mutations can occur either early or late in the evolution of a tumor•Genetic heterogeneity in cancer may limit effectiveness of targeted therapy
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.
Determining how somatic copy number alterations (SCNAs) promote cancer is an important goal. We characterized SCNA patterns in 4,934 cancers from The Cancer Genome Atlas Pan-Cancer data set. ...Whole-genome doubling, observed in 37% of cancers, was associated with higher rates of every other type of SCNA, TP53 mutations, CCNE1 amplifications and alterations of the PPP2R complex. SCNAs that were internal to chromosomes tended to be shorter than telomere-bounded SCNAs, suggesting different mechanisms underlying their generation. Significantly recurrent focal SCNAs were observed in 140 regions, including 102 without known oncogene or tumor suppressor gene targets and 50 with significantly mutated genes. Amplified regions without known oncogenes were enriched for genes involved in epigenetic regulation. When levels of genomic disruption were accounted for, 7% of region pairs were anticorrelated, and these regions tended to encompass genes whose proteins physically interact, suggesting related functions. These results provide insights into mechanisms of generation and functional consequences of cancer-related SCNAs.
Response to immune checkpoint blockade in mesenchymal tumors is poorly characterized, but immunogenomic dissection of these cancers could inform immunotherapy mediators. We identified a ...treatment-naive patient who has metastatic uterine leiomyosarcoma and has experienced complete tumor remission for >2 years on anti-PD-1 (pembrolizumab) monotherapy. We analyzed the primary tumor, the sole treatment-resistant metastasis, and germline tissue to explore mechanisms of immunotherapy sensitivity and resistance. Both tumors stained diffusely for PD-L2 and showed sparse PD-L1 staining. PD-1+ cell infiltration significantly decreased in the resistant tumor (p = 0.039). Genomically, the treatment-resistant tumor uniquely harbored biallelic PTEN loss and had reduced expression of two neoantigens that demonstrated strong immunoreactivity with patient T cells in vitro, suggesting long-lasting immunological memory. In this near-complete response to PD-1 blockade in a mesenchymal tumor, we identified PTEN mutations and reduced expression of genes encoding neoantigens as potential mediators of resistance to immune checkpoint therapy.
Display omitted
•Anti-PD-1 monotherapy can induce complete tumor remission in uterine leiomyosarcoma•Longitudinal profiling of patient tumors reveals response and resistance mechanisms•Tumor neoantigens are patient specific and can induce long-term immunoreactivity•Biallelic PTEN loss is associated with resistance to PD-1 blockade in sarcoma
George et al. report an exceptional responder to anti-PD-1 monotherapy in uterine leiomyosarcoma and propose mediators of treatment sensitivity and resistance. Neoantigen-directed immunoreactivity was associated with sensitivity to anti-PD-1 monotherapy, whereas resistance was associated with reduction in neoantigen expression consistent with immune evasion, and biallelic PTEN loss was associated with induction of an immunosuppressive microenvironment.
Infiltrating stromal and immune cells form the major fraction of normal cells in tumour tissue and not only perturb the tumour signal in molecular studies but also have an important role in cancer ...biology. Here we describe 'Estimation of STromal and Immune cells in MAlignant Tumours using Expression data' (ESTIMATE)--a method that uses gene expression signatures to infer the fraction of stromal and immune cells in tumour samples. ESTIMATE scores correlate with DNA copy number-based tumour purity across samples from 11 different tumour types, profiled on Agilent, Affymetrix platforms or based on RNA sequencing and available through The Cancer Genome Atlas. The prediction accuracy is further corroborated using 3,809 transcriptional profiles available elsewhere in the public domain. The ESTIMATE method allows consideration of tumour-associated normal cells in genomic and transcriptomic studies. An R-library is available on https://sourceforge.net/projects/estimateproject/.
Everolimus, an inhibitor of the mammalian target of rapamycin (mTOR), is effective in treating tumors harboring alterations in the mTOR pathway. Mechanisms of resistance to everolimus remain ...undefined. Resistance developed in a patient with metastatic anaplastic thyroid carcinoma after an extraordinary 18-month response. Whole-exome sequencing of pretreatment and drug-resistant tumors revealed a nonsense mutation in TSC2, a negative regulator of mTOR, suggesting a mechanism for exquisite sensitivity to everolimus. The resistant tumor also harbored a mutation in MTOR that confers resistance to allosteric mTOR inhibition. The mutation remains sensitive to mTOR kinase inhibitors.
We developed a computational method to characterize aneuploidy in tumor samples based on coordinated aberrations in expression of genes localized to each chromosomal region. We summarized the total ...level of chromosomal aberration in a given tumor in a univariate measure termed total functional aneuploidy. We identified a signature of chromosomal instability from specific genes whose expression was consistently correlated with total functional aneuploidy in several cancer types. Net overexpression of this signature was predictive of poor clinical outcome in 12 cancer data sets representing six cancer types. Also, the signature of chromosomal instability was higher in metastasis samples than in primary tumors and was able to stratify grade 1 and grade 2 breast tumors according to clinical outcome. These results provide a means to assess the potential role of chromosomal instability in determining malignant potential over a broad range of tumors.