Jin Shao

Clonal Architecture of Secondary Acute Myeloid Leukemia

BACKGROUND The myelodysplastic syndromes are a group of hematologic disorders that often evolve into secondary acute myeloid leukemia (AML). The genetic changes that underlie progression from the myelodysplastic syndromes to secondary AML are not well understood. METHODS We performed whole-genome sequencing of seven paired samples of skin and bone marrow in seven subjects with secondary AML to identify somatic mutations specific to secondary AML. We then genotyped a bone marrow sample obtained during the antecedent myelodysplastic-syndrome stage from each subject to determine the presence or absence of the specific somatic mutations. We identified recurrent mutations in coding genes and defined the clonal architecture of each pair of samples from the myelodysplastic-syndrome stage and the secondary-AML stage, using the allele burden of hundreds of mutations. RESULTS Approximately 85% of bone marrow cells were clonal in the myelodysplastic-syndrome and secondary-AML samples, regardless of the myeloblast count. The secondary-AML samples contained mutations in 11 recurrently mutated genes, including 4 genes that have not been previously implicated in the myelodysplastic syndromes or AML. In every case, progression to acute leukemia was defined by the persistence of an antecedent founding clone containing 182 to 660 somatic mutations and the outgrowth or emergence of at least one subclone, harboring dozens to hundreds of new mutations. All founding clones and subclones contained at least one mutation in a coding gene. CONCLUSIONS Nearly all the bone marrow cells in patients with myelodysplastic syndromes and secondary AML are clonally derived. Genetic evolution of secondary AML is a dynamic process shaped by multiple cycles of mutation acquisition and clonal selection. Recurrent gene mutations are found in both founding clones and daughter subclones. (Funded by the National Institutes of Health and others.).

Recurrent DNMT3A mutations in patients with myelodysplastic syndromes

Alterations in DNA methylation have been implicated in the pathogenesis of myelodysplastic syndromes (MDS), although the underlying mechanism remains largely unknown. Methylation of CpG dinucleotides is mediated by DNA methyltransferases, including DNMT1, DNMT3A and DNMT3B. DNMT3A mutations have recently been reported in patients with de novo acute myeloid leukemia (AML), providing a rationale for examining the status of DNMT3A in MDS samples. In this study, we report the frequency of DNMT3A mutations in patients with de novo MDS, and their association with secondary AML. We sequenced all coding exons of DNMT3A using DNA from bone marrow and paired normal cells from 150 patients with MDS and identified 13 heterozygous mutations with predicted translational consequences in 12/150 patients (8.0%). Amino acid R882, located in the methyltransferase domain of DNMT3A, was the most common mutation site, accounting for 4/13 mutations. DNMT3A mutations were expressed in the majority of cells in all tested mutant samples regardless of myeloblast counts, suggesting that DNMT3A mutations occur early in the course of MDS. Patients with DNMT3A mutations had worse overall survival compared with patients without DNMT3A mutations (P=0.005) and more rapid progression to AML (P=0.007), suggesting that DNMT3A mutation status may have prognostic value in de novo MDS.

RECURRENT MUTATIONS IN THE U2AF1 SPLICING FACTOR IN MYELODYSPLASTIC SYNDROMES

Myelodysplastic syndromes (MDS) are hematopoietic stem cell disorders that often progress to chemotherapy-resistant secondary acute myeloid leukemia (sAML). We used whole-genome sequencing to perform an unbiased comprehensive screen to discover the somatic mutations in a sample from an individual with sAML and genotyped the loci containing these mutations in the matched MDS sample. Here we show that a missense mutation affecting the serine at codon 34 (Ser34) in U2AF1 was recurrently present in 13 out of 150 (8.7%) subjects with de novo MDS, and we found suggestive evidence of an increased risk of progression to sAML associated with this mutation. U2AF1 is a U2 auxiliary factor protein that recognizes the AG splice acceptor dinucleotide at the 3′ end of introns, and the alterations in U2AF1 are located in highly conserved zinc fingers of this protein. Mutant U2AF1 promotes enhanced splicing and exon skipping in reporter assays in vitro. This previously unidentified, recurrent mutation in U2AF1 implicates altered pre-mRNA splicing as a potential mechanism for MDS pathogenesis.

Clonal diversity of recurrently mutated genes in myelodysplastic syndromes

Recent studies suggest that most cases of myelodysplastic syndrome (MDS) are clonally heterogeneous, with a founding clone and multiple subclones. It is not known whether specific gene mutations typically occur in founding clones or subclones. We screened a panel of 94 candidate genes in a cohort of 157 patients with MDS or secondary acute myeloid leukemia (sAML). This included 150 cases with samples obtained at MDS diagnosis and 15 cases with samples obtained at sAML transformation (8 were also analyzed at the MDS stage). We performed whole-genome sequencing (WGS) to define the clonal architecture in eight sAML genomes and identified the range of variant allele frequencies (VAFs) for founding clone mutations. At least one mutation or cytogenetic abnormality was detected in 83% of the 150 MDS patients and 17 genes were significantly mutated (false discovery rate ⩽0.05). Individual genes and patient samples displayed a wide range of VAFs for recurrently mutated genes, indicating that no single gene is exclusively mutated in the founding clone. The VAFs of recurrently mutated genes did not fully recapitulate the clonal architecture defined by WGS, suggesting that comprehensive sequencing may be required to accurately assess the clonal status of recurrently mutated genes in MDS.

Mutant U2AF1 Expression Alters Hematopoiesis and Pre-mRNA Splicing In Vivo

Heterozygous somatic mutations in the spliceosome gene U2AF1 occur in ∼ 11% of patients with myelodysplastic syndromes (MDS), the most common adult myeloid malignancy. It is unclear how these mutations contribute to disease. We examined in vivo hematopoietic consequences of the most common U2AF1 mutation using a doxycycline-inducible transgenic mouse model. Mice expressing mutant U2AF1(S34F) display altered hematopoiesis and changes in pre-mRNA splicing in hematopoietic progenitor cells by whole transcriptome analysis (RNA-seq). Integration with human RNA-seq datasets determined that common mutant U2AF1-induced splicing alterations are enriched in RNA processing genes, ribosomal genes, and recurrently mutated MDS and acute myeloid leukemia-associated genes. These findings support the hypothesis that mutant U2AF1 alters downstream gene isoform expression, thereby contributing to abnormal hematopoiesis in patients with MDS.

Genomic analysis of germ line and somatic variants in familial myelodysplasia/acute myeloid leukemia

Familial clustering of myelodysplastic syndromes (MDSs) and acute myeloid leukemia (AML) can be caused by inherited factors. We screened 59 individuals from 17 families with 2 or more biological relatives with MDS/AML for variants in 12 genes with established roles in predisposition to MDS/AML, and identified a pathogenic germ line variant in 5 families (29%). Extending the screen with a panel of 264 genes that are recurrently mutated in de novo AML, we identified rare, nonsynonymous germ line variants in 4 genes, each segregating with MDS/AML in 2 families. Somatic mutations are required for progression to MDS/AML in these familial cases. Using a combination of targeted and exome sequencing of tumor and matched normal samples from 26 familial MDS/AML cases and asymptomatic carriers, we identified recurrent frameshift mutations in the cohesin-associated factor PDS5B, co-occurrence of somatic ASXL1 mutations with germ line GATA2 mutations, and recurrent mutations in other known MDS/AML drivers. Mutations in genes that are recurrently mutated in de novo AML were underrepresented in the familial MDS/AML cases, although the total number of somatic mutations per exome was the same. Lastly, clonal skewing of hematopoiesis was detected in 67% of young, asymptomatic RUNX1 carriers, providing a potential biomarker that could be used for surveillance in these high-risk families.

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.

Genome-wide association study to identify novel loci associated with therapy-related myeloid leukemia susceptibility

Therapy-related acute myeloid leukemia (t-AML) is a rare but fatal complication of cytotoxic therapy. Whereas sporadic cancer results from interactions between complex exposures and low-penetrance alleles, t-AML results from an acute exposure to a limited number of potent genotoxins. Consequently, we hypothesized that the effect sizes of variants associated with t-AML would be greater than in sporadic cancer, and, therefore, that these variants could be detected even in a modest-sized cohort. To test this, we undertook an association study in 80 cases and 150 controls using Affymetrix Mapping 10K arrays. Even at nominal significance thresholds, we found a significant excess of associations over chance; for example, although 6 associations were expected at P less than .001, we found 15 (P(enrich) = .002). To replicate our findings, we genotyped the 10 most significantly associated single nucleotide polymorphisms (SNPs) in an independent t-AML cohort (n = 70) and obtained evidence of association with t-AML for 3 SNPs in the subset of patients with loss of chromosomes 5 or 7 or both, acquired abnormalities associated with prior exposure to alkylator chemotherapy. Thus, we conclude that the effect of genetic factors contributing to cancer risk is potentiated and more readily discernable in t-AML compared with sporadic cancer.

Integrated Genomic Analysis Implicates Haploinsufficiency of Multiple Chromosome 5q31.2 Genes in De Novo Myelodysplastic Syndromes Pathogenesis

Deletions spanning chromosome 5q31.2 are among the most common recurring cytogenetic abnormalities detectable in myelodysplastic syndromes (MDS). Prior genomic studies have suggested that haploinsufficiency of multiple 5q31.2 genes may contribute to MDS pathogenesis. However, this hypothesis has never been formally tested. Therefore, we designed this study to systematically and comprehensively evaluate all 28 chromosome 5q31.2 genes and directly test whether haploinsufficiency of a single 5q31.2 gene may result from a heterozygous nucleotide mutation or microdeletion. We selected paired tumor (bone marrow) and germline (skin) DNA samples from 46 de novo MDS patients (37 without a cytogenetic 5q31.2 deletion) and performed total exonic gene resequencing (479 amplicons) and array comparative genomic hybridization (CGH). We found no somatic nucleotide changes in the 46 MDS samples, and no cytogenetically silent 5q31.2 deletions in 20/20 samples analyzed by array CGH. Twelve novel single nucleotide polymorphisms were discovered. The mRNA levels of 7 genes in the commonly deleted interval were reduced by 50% in CD34+ cells from del(5q) MDS samples, and no gene showed complete loss of expression. Taken together, these data show that small deletions and/or point mutations in individual 5q31.2 genes are not common events in MDS, and implicate haploinsufficiency of multiple genes as the relevant genetic consequence of this common deletion.

Dynamic changes in the clonal structure of MDS and AML in response to epigenetic therapy

Traditional response criteria in myelodysplastic syndrome (MDS) and acute myeloid leukemia (AML) are based on bone marrow morphology and may not accurately reflect clonal tumor burden in patients treated with non-cytotoxic chemotherapy. We used next-generation sequencing of serial bone marrow samples to monitor MDS and AML tumor burden during treatment with epigenetic therapy (decitabine and panobinostat). Serial bone marrow samples (and skin as a source of normal DNA) from 25 MDS and AML patients were sequenced (exome or 285 gene panel). We observed that responders, including those in complete remission (CR), can have persistent measurable tumor burden (that is, mutations) for at least 1 year without disease progression. Using an ultrasensitive sequencing approach, we detected extremely rare mutations (equivalent to 1 heterozygous mutant cell in 2000 non-mutant cells) months to years before their expansion at disease relapse. While patients can live with persistent clonal hematopoiesis in a CR or stable disease, ultimately we find evidence that expansion of a rare subclone occurs at relapse or progression. Here we demonstrate that sequencing of serial samples provides an alternative measure of tumor burden in MDS or AML patients and augments traditional response criteria that rely on bone marrow blast percentage.