The origin and evolution of mutations in acute myeloid leukemia.

Most mutations in cancer genomes are thought to be acquired after the initiating event, which may cause genomic instability and drive clonal evolution. However, for acute myeloid leukemia (AML), normal karyotypes are common, and genomic instability is unusual. To better understand clonal evolution in AML, we sequenced the genomes of M3-AML samples with a known initiating event (PML-RARA) versus the genomes of normal karyotype M1-AML samples and the exomes of hematopoietic stem/progenitor cells (HSPCs) from healthy people. Collectively, the data suggest that most of the mutations found in AML genomes are actually random events that occurred in HSPCs before they acquired the initiating mutation; the mutational history of that cell is "captured" as the clone expands. In many cases, only one or two additional, cooperating mutations are needed to generate the malignant founding clone. Cells from the founding clone can acquire additional cooperating mutations, yielding subclones that can contribute to disease progression and/or relapse.

Immunological ignorance is an enabling feature of the oligo-clonal T cell response to melanoma neoantigens.

The impact of intratumoral heterogeneity (ITH) and the resultant neoantigen landscape on T cell immunity are poorly understood. ITH is a widely recognized feature of solid tumors and poses distinct challenges related to the development of effective therapeutic strategies, including cancer neoantigen vaccines. Here, we performed deep targeted DNA sequencing of multiple metastases from melanoma patients and observed ubiquitous sharing of clonal and subclonal single nucleotide variants (SNVs) encoding putative HLA class I-restricted neoantigen epitopes. However, spontaneous antitumor CD8+ T cell immunity in peripheral blood and tumors was restricted to a few clonal neoantigens featuring an oligo-/monoclonal T cell-receptor (TCR) repertoire. Moreover, in various tumors of the 4 patients examined, no neoantigen-specific TCR clonotypes were identified despite clonal neoantigen expression. Mature dendritic cell (mDC) vaccination with tumor-encoded amino acid-substituted (AAS) peptides revealed diverse neoantigen-specific CD8+ T responses, each composed of multiple TCR clonotypes. Isolation of T cell clones by limiting dilution from tumor-infiltrating lymphocytes (TILs) permitted functional validation regarding neoantigen specificity. Gene transfer of TCRαß heterodimers specific for clonal neoantigens confirmed correct TCR clonotype assignments based on high-throughput TCRBV CDR3 sequencing. Our findings implicate immunological ignorance of clonal neoantigens as the basis for ineffective T cell immunity to melanoma and support the concept that therapeutic vaccination, as an adjunct to checkpoint inhibitor treatment, is required to increase the breadth and diversity of neoantigen-specific CD8+ T cells.

Rapid expansion of preexisting nonleukemic hematopoietic clones frequently follows induction therapy for de novo AML.

There is interest in using leukemia-gene panels and next-generation sequencing to assess acute myelogenous leukemia (AML) response to induction chemotherapy. Studies have shown that patients with AML in morphologic remission may continue to have clonal hematopoiesis with populations closely related to the founding AML clone and that this confers an increased risk of relapse. However, it remains unknown how induction chemotherapy influences the clonal evolution of a patient's nonleukemic hematopoietic population. Here, we report that 5 of 15 patients with genetic clearance of their founding AML clone after induction chemotherapy had a concomitant expansion of a hematopoietic population unrelated to the initial AML. These populations frequently harbored somatic mutations in genes recurrently mutated in AML or myelodysplastic syndromes and were detectable at very low frequencies at the time of AML diagnosis. These results suggest that nonleukemic hematopoietic stem and progenitor cells, harboring specific aging-acquired mutations, may have a competitive fitness advantage after induction chemotherapy, expand, and persist long after the completion of chemotherapy. Although the clinical importance of these "rising" clones remains to be determined, it will be important to distinguish them from leukemia-related populations when assessing for molecular responses to induction chemotherapy.

Cancer exome analysis reveals a T-cell-dependent mechanism of cancer immunoediting.

Cancer immunoediting, the process by which the immune system controls tumour outgrowth and shapes tumour immunogenicity, is comprised of three phases: elimination, equilibrium and escape. Although many immune components that participate in this process are known, its underlying mechanisms remain poorly defined. A central tenet of cancer immunoediting is that T-cell recognition of tumour antigens drives the immunological destruction or sculpting of a developing cancer. However, our current understanding of tumour antigens comes largely from analyses of cancers that develop in immunocompetent hosts and thus may have already been edited. Little is known about the antigens expressed in nascent tumour cells, whether they are sufficient to induce protective antitumour immune responses or whether their expression is modulated by the immune system. Here, using massively parallel sequencing, we characterize expressed mutations in highly immunogenic methylcholanthrene-induced sarcomas derived from immunodeficient Rag2(-/-) mice that phenotypically resemble nascent primary tumour cells. Using class I prediction algorithms, we identify mutant spectrin-ß2 as a potential rejection antigen of the d42m1 sarcoma and validate this prediction by conventional antigen expression cloning and detection. We also demonstrate that cancer immunoediting of d42m1 occurs via a T-cell-dependent immunoselection process that promotes outgrowth of pre-existing tumour cell clones lacking highly antigenic mutant spectrin-ß2 and other potential strong antigens. These results demonstrate that the strong immunogenicity of an unedited tumour can be ascribed to expression of highly antigenic mutant proteins and show that outgrowth of tumour cells that lack these strong antigens via a T-cell-dependent immunoselection process represents one mechanism of cancer immunoediting.

Whole-genome analysis informs breast cancer response to aromatase inhibition.

To correlate the variable clinical features of oestrogen-receptor-positive breast cancer with somatic alterations, we studied pretreatment tumour biopsies accrued from patients in two studies of neoadjuvant aromatase inhibitor therapy by massively parallel sequencing and analysis. Eighteen significantly mutated genes were identified, including five genes (RUNX1, CBFB, MYH9, MLL3 and SF3B1) previously linked to haematopoietic disorders. Mutant MAP3K1 was associated with luminal A status, low-grade histology and low proliferation rates, whereas mutant TP53 was associated with the opposite pattern. Moreover, mutant GATA3 correlated with suppression of proliferation upon aromatase inhibitor treatment. Pathway analysis demonstrated that mutations in MAP2K4, a MAP3K1 substrate, produced similar perturbations as MAP3K1 loss. Distinct phenotypes in oestrogen-receptor-positive breast cancer are associated with specific patterns of somatic mutations that map into cellular pathways linked to tumour biology, but most recurrent mutations are relatively infrequent. Prospective clinical trials based on these findings will require comprehensive genome sequencing.

Correction: Tumor Evolution in Two Patients with Basal-like Breast Cancer: A Retrospective Genomics Study of Multiple Metastases.

[This corrects the article DOI: 10.1371/journal.pmed.1002174.].

Clonal architecture of secondary acute myeloid leukemia defined by single-cell sequencing.

Next-generation sequencing has been used to infer the clonality of heterogeneous tumor samples. These analyses yield specific predictions-the population frequency of individual clones, their genetic composition, and their evolutionary relationships-which we set out to test by sequencing individual cells from three subjects diagnosed with secondary acute myeloid leukemia, each of whom had been previously characterized by whole genome sequencing of unfractionated tumor samples. Single-cell mutation profiling strongly supported the clonal architecture implied by the analysis of bulk material. In addition, it resolved the clonal assignment of single nucleotide variants that had been initially ambiguous and identified areas of previously unappreciated complexity. Accordingly, we find that many of the key assumptions underlying the analysis of tumor clonality by deep sequencing of unfractionated material are valid. Furthermore, we illustrate a single-cell sequencing strategy for interrogating the clonal relationships among known variants that is cost-effective, scalable, and adaptable to the analysis of both hematopoietic and solid tumors, or any heterogeneous population of cells.

Tumor Evolution in Two Patients with Basal-like Breast Cancer: A Retrospective Genomics Study of Multiple Metastases.

BACKGROUND: Metastasis is the main cause of cancer patient deaths and remains a poorly characterized process. It is still unclear when in tumor progression the ability to metastasize arises and whether this ability is inherent to the primary tumor or is acquired well after primary tumor formation. Next-generation sequencing and analytical methods to define clonal heterogeneity provide a means for identifying genetic events and the temporal relationships between these events in the primary and metastatic tumors within an individual. METHODS AND FINDINGS: We performed DNA whole genome and mRNA sequencing on two primary tumors, each with either four or five distinct tissue site-specific metastases, from two individuals with triple-negative/basal-like breast cancers. As evidenced by their case histories, each patient had an aggressive disease course with abbreviated survival. In each patient, the overall gene expression signatures, DNA copy number patterns, and somatic mutation patterns were highly similar across each primary tumor and its associated metastases. Almost every mutation found in the primary was found in a metastasis (for the two patients, 52/54 and 75/75). Many of these mutations were found in every tumor (11/54 and 65/75, respectively). In addition, each metastasis had fewer metastatic-specific events and shared at least 50% of its somatic mutation repertoire with the primary tumor, and all samples from each patient grouped together by gene expression clustering analysis. TP53 was the only mutated gene in common between both patients and was present in every tumor in this study. Strikingly, each metastasis resulted from multiclonal seeding instead of from a single cell of origin, and few of the new mutations, present only in the metastases, were expressed in mRNAs. Because of the clinical differences between these two patients and the small sample size of our study, the generalizability of these findings will need to be further examined in larger cohorts of patients. CONCLUSIONS: Our findings suggest that multiclonal seeding may be common amongst basal-like breast cancers. In these two patients, mutations and DNA copy number changes in the primary tumors appear to have had a biologic impact on metastatic potential, whereas mutations arising in the metastases were much more likely to be passengers.

Somatic neurofibromatosis type 1 (NF1) inactivation characterizes NF1-associated pilocytic astrocytoma.

Low-grade brain tumors (pilocytic astrocytomas) arising in the neurofibromatosis type 1 (NF1) inherited cancer predisposition syndrome are hypothesized to result from a combination of germline and acquired somatic NF1 tumor suppressor gene mutations. However, genetically engineered mice (GEM) in which mono-allelic germline Nf1 gene loss is coupled with bi-allelic somatic (glial progenitor cell) Nf1 gene inactivation develop brain tumors that do not fully recapitulate the neuropathological features of the human condition. These observations raise the intriguing possibility that, while loss of neurofibromin function is necessary for NF1-associated low-grade astrocytoma development, additional genetic changes may be required for full penetrance of the human brain tumor phenotype. To identify these potential cooperating genetic mutations, we performed whole-genome sequencing (WGS) analysis of three NF1-associated pilocytic astrocytoma (PA) tumors. We found that the mechanism of somatic NF1 loss was different in each tumor (frameshift mutation, loss of heterozygosity, and methylation). In addition, tumor purity analysis revealed that these tumors had a high proportion of stromal cells, such that only 50%-60% of cells in the tumor mass exhibited somatic NF1 loss. Importantly, we identified no additional recurrent pathogenic somatic mutations, supporting a model in which neuroglial progenitor cell NF1 loss is likely sufficient for PA formation in cooperation with a proper stromal environment.

Association Between Mutation Clearance After Induction Therapy and Outcomes in Acute Myeloid Leukemia.

IMPORTANCE: Tests that predict outcomes for patients with acute myeloid leukemia (AML) are imprecise, especially for those with intermediate risk AML. OBJECTIVES: To determine whether genomic approaches can provide novel prognostic information for adult patients with de novo AML. DESIGN, SETTING, AND PARTICIPANTS: Whole-genome or exome sequencing was performed on samples obtained at disease presentation from 71 patients with AML (mean age, 50.8 years) treated with standard induction chemotherapy at a single site starting in March 2002, with follow-up through January 2015. In addition, deep digital sequencing was performed on paired diagnosis and remission samples from 50 patients (including 32 with intermediate-risk AML), approximately 30 days after successful induction therapy. Twenty-five of the 50 were from the cohort of 71 patients, and 25 were new, additional cases. EXPOSURES: Whole-genome or exome sequencing and targeted deep sequencing. Risk of identification based on genetic data. MAIN OUTCOMES AND MEASURES: Mutation patterns (including clearance of leukemia-associated variants after chemotherapy) and their association with event-free survival and overall survival. RESULTS: Analysis of comprehensive genomic data from the 71 patients did not improve outcome assessment over current standard-of-care metrics. In an analysis of 50 patients with both presentation and documented remission samples, 24 (48%) had persistent leukemia-associated mutations in at least 5% of bone marrow cells at remission. The 24 with persistent mutations had significantly reduced event-free and overall survival vs the 26 who cleared all mutations. Patients with intermediate cytogenetic risk profiles had similar findings. [table: see text]. CONCLUSIONS AND RELEVANCE: The detection of persistent leukemia-associated mutations in at least 5% of bone marrow cells in day 30 remission samples was associated with a significantly increased risk of relapse, and reduced overall survival. These data suggest that this genomic approach may improve risk stratification for patients with AML.

PML-RARA requires DNA methyltransferase 3A to initiate acute promyelocytic leukemia.

The DNA methyltransferases DNMT3A and DNMT3B are primarily responsible for de novo methylation of specific cytosine residues in CpG dinucleotides during mammalian development. While loss-of-function mutations in DNMT3A are highly recurrent in acute myeloid leukemia (AML), DNMT3A mutations are almost never found in AML patients with translocations that create oncogenic fusion genes such as PML-RARA, RUNX1-RUNX1T1, and MLL-AF9. Here, we explored how DNMT3A is involved in the function of these fusion genes. We used retroviral vectors to express PML-RARA, RUNX1-RUNX1T1, or MLL-AF9 in bone marrow cells derived from WT or DNMT3A-deficient mice. Additionally, we examined the phenotypes of hematopoietic cells from Ctsg-PML-RARA mice, which express PML-RARA in early hematopoietic progenitors and myeloid precursors, with or without DNMT3A. We determined that the methyltransferase activity of DNMT3A, but not DNMT3B, is required for aberrant PML-RARA-driven self-renewal ex vivo and that DNMT3A is dispensable for RUNX1-RUNX1T1- and MLL-AF9-driven self-renewal. Furthermore, both the PML-RARA-driven competitive transplantation advantage and development of acute promyelocytic leukemia (APL) required DNMT3A. Together, these findings suggest that PML-RARA requires DNMT3A to initiate APL in mice.

A common founding clone with TP53 and PTEN mutations gives rise to a concurrent germ cell tumor and acute megakaryoblastic leukemia.

We report the findings from a patient who presented with a concurrent mediastinal germ cell tumor (GCT) and acute myeloid leukemia (AML). Bone marrow pathology was consistent with a diagnosis of acute megakaryoblastic leukemia (AML M7), and biopsy of an anterior mediastinal mass was consistent with a nonseminomatous GCT. Prior studies have described associations between hematological malignancies, including AML M7 and nonseminomatous GCTs, and it was recently suggested that a common founding clone initiated both cancers. We performed enhanced exome sequencing on the GCT and the AML M7 from our patient to define the clonal relationship between the two cancers. We found that both samples contained somatic mutations in PTEN (C136R missense) and TP53 (R213 frameshift). The mutations in PTEN and TP53 were present at ∼100% variant allele frequency (VAF) in both tumors. In addition, we detected and validated five other shared somatic mutations. The copy-number analysis of the AML exome data revealed an amplification of Chromosome 12p. We also identified a heterozygous germline variant in FANCA (S858R), which is known to be associated with Fanconi anemia but is of uncertain significance here. In summary, our data not only support a common founding clone for these cancers but also suggest that a specific set of distinct genomic alterations (in PTEN and TP53) underlies the rare association between GCT and AML. This association is likely linked to the treatment resistance and extremely poor outcome of these patients. We cannot resolve the clonal evolution of these tumors given limitations of our data.

Identification of a novel fusion transcript between human relaxin-1 (RLN1) and human relaxin-2 (RLN2) in prostate cancer.

Simultaneous expression of highly homologous RLN1 and RLN2 genes in prostate impairs their accurate delineation. We used PacBio SMRT sequencing and RNA-Seq in LNCaP cells in order to dissect the expression of RLN1 and RLN2 variants. We identified a novel fusion transcript comprising the RLN1 and RLN2 genes and found evidence of its expression in the normal and prostate cancer tissues. The RLN1-RLN2 fusion putatively encodes RLN2 isoform with the deleted secretory signal peptide. The identification of the fusion transcript provided information to determine unique RLN1-RLN2 fusion and RLN1 regions. The RLN1-RLN2 fusion was co-expressed with RLN1 in LNCaP cells, but the two gene products were inversely regulated by androgens. We showed that RLN1 is underrepresented in common PCa cell lines in comparison to normal and PCa tissue. The current study brings a highly relevant update to the relaxin field, and will encourage further studies of RLN1 and RLN2 in PCa and broader.

Optimizing cancer genome sequencing and analysis.

Tumors are typically sequenced to depths of 75-100× (exome) or 30-50× (whole genome). We demonstrate that current sequencing paradigms are inadequate for tumors that are impure, aneuploid or clonally heterogeneous. To reassess optimal sequencing strategies, we performed ultra-deep (up to ~312×) whole genome sequencing (WGS) and exome capture (up to ~433×) of a primary acute myeloid leukemia, its subsequent relapse, and a matched normal skin sample. We tested multiple alignment and variant calling algorithms and validated ~200,000 putative SNVs by sequencing them to depths of ~1,000×. Additional targeted sequencing provided over 10,000× coverage and ddPCR assays provided up to ~250,000× sampling of selected sites. We evaluated the effects of different library generation approaches, depth of sequencing, and analysis strategies on the ability to effectively characterize a complex tumor. This dataset, representing the most comprehensively sequenced tumor described to date, will serve as an invaluable community resource (dbGaP accession id phs000159).

cDNA hybrid capture improves transcriptome analysis on low-input and archived samples.

The use of massively parallel sequencing for studying RNA expression has greatly enhanced our understanding of the transcriptome through the myriad ways these data can be characterized. In particular, clinical samples provide important insights about RNA expression in health and disease, yet these studies can be complicated by RNA degradation that results from the use of formalin as a clinical preservative and by the limited amounts of RNA often available from these precious samples. In this study we describe the combined use of RNA sequencing with an exome capture selection step to enhance the yield of on-exon sequencing read data when compared with RNA sequencing alone. In particular, the exome capture step preserves the dynamic range of expression, permitting differential comparisons and validation of expressed mutations from limited and FFPE preserved samples, while reducing the data generation requirement. We conclude that cDNA hybrid capture has the potential to significantly improve transcriptome analysis from low-yield FFPE material.

Functional heterogeneity of genetically defined subclones in acute myeloid leukemia.

The relationships between clonal architecture and functional heterogeneity in acute myeloid leukemia (AML) samples are not yet clear. We used targeted sequencing to track AML subclones identified by whole-genome sequencing using a variety of experimental approaches. We found that virtually all AML subclones trafficked from the marrow to the peripheral blood, but some were enriched in specific cell populations. Subclones showed variable engraftment potential in immunodeficient mice. Xenografts were predominantly comprised of a single genetically defined subclone, but there was no predictable relationship between the engrafting subclone and the evolutionary hierarchy of the leukemia. These data demonstrate the importance of integrating genetic and functional data in studies of primary cancer samples, both in xenograft models and in patients.

Background mutations in parental cells account for most of the genetic heterogeneity of induced pluripotent stem cells.

To assess the genetic consequences of induced pluripotent stem cell (iPSC) reprogramming, we sequenced the genomes of ten murine iPSC clones derived from three independent reprogramming experiments, and compared them to their parental cell genomes. We detected hundreds of single nucleotide variants (SNVs) in every clone, with an average of 11 in coding regions. In two experiments, all SNVs were unique for each clone and did not cluster in pathways, but in the third, all four iPSC clones contained 157 shared genetic variants, which could also be detected in rare cells (<1 in 500) within the parental MEF pool. These data suggest that most of the genetic variation in iPSC clones is not caused by reprogramming per se, but is rather a consequence of cloning individual cells, which "captures" their mutational history. These findings have implications for the development and therapeutic use of cells that are reprogrammed by any method.

De novo gene disruptions in children on the autistic spectrum.

Exome sequencing of 343 families, each with a single child on the autism spectrum and at least one unaffected sibling, reveal de novo small indels and point substitutions, which come mostly from the paternal line in an age-dependent manner. We do not see significantly greater numbers of de novo missense mutations in affected versus unaffected children, but gene-disrupting mutations (nonsense, splice site, and frame shifts) are twice as frequent, 59 to 28. Based on this differential and the number of recurrent and total targets of gene disruption found in our and similar studies, we estimate between 350 and 400 autism susceptibility genes. Many of the disrupted genes in these studies are associated with the fragile X protein, FMRP, reinforcing links between autism and synaptic plasticity. We find FMRP-associated genes are under greater purifying selection than the remainder of genes and suggest they are especially dosage-sensitive targets of cognitive disorders.

Hybrid capture and next-generation sequencing identify viral integration sites from formalin-fixed, paraffin-embedded tissue.

Although next-generation sequencing (NGS) has been the domain of large genome centers, it is quickly becoming more accessible to general pathology laboratories. In addition to finding single-base changes, NGS allows for the detection of larger structural variants, including insertions/deletions, translocations, and viral insertions. We describe the use of targeted NGS on DNA extracted from formalin-fixed, paraffin-embedded (FFPE) tissue, and show that the short read lengths of NGS are ideally suited to fragmented DNA obtained from FFPE tissue. Further, we describe a novel method for performing hybrid-capture target enrichment using PCR-generated capture probes. As a model, we captured the 5.3-kb Merkel cell polyomavirus (MCPyV) genome in FFPE cases of Merkel cell carcinoma using inexpensive, PCR-derived capture probes, and achieved up to 37,000-fold coverage of the MCPyV genome without prior virus-specific PCR amplification. This depth of coverage made it possible to reproducibly detect viral genome deletions and insertion sites anywhere within the human genome. Out of four cases sequenced, we identified the 5' insertion sites in four of four cases and the 3' sites in three of four cases. These findings demonstrate the potential for an inexpensive gene targeting and NGS method that can be easily adapted for use with FFPE tissue to identify large structural rearrangements, opening up the possibility for further discovery from archival tissue.