CpG Island Hypermethylation Mediated by DNMT3A Is a Consequence of AML Progression.

DNMT3A mutations occur in ∼25% of acute myeloid leukemia (AML) patients. The most common mutation, DNMT3AR882H, has dominant negative activity that reduces DNA methylation activity by ∼80% in vitro. To understand the contribution of DNMT3A-dependent methylation to leukemogenesis, we performed whole-genome bisulfite sequencing of primary leukemic and non-leukemic cells in patients with or without DNMT3AR882 mutations. Non-leukemic hematopoietic cells with DNMT3AR882H displayed focal methylation loss, suggesting that hypomethylation antedates AML. Although virtually all AMLs with wild-type DNMT3A displayed CpG island hypermethylation, this change was not associated with gene silencing and was essentially absent in AMLs with DNMT3AR882 mutations. Primary hematopoietic stem cells expanded with cytokines were hypermethylated in a DNMT3A-dependent manner, suggesting that hypermethylation may be a response to, rather than a cause of, cellular proliferation. Our findings suggest that hypomethylation is an initiating phenotype in AMLs with DNMT3AR882, while DNMT3A-dependent CpG island hypermethylation is a consequence of AML progression.

PML::RARA and GATA2 proteins interact via DNA templates to induce aberrant self-renewal in mouse and human hematopoietic cells.

The underlying mechanism(s) by which the PML::RARA fusion protein initiates acute promyelocytic leukemia is not yet clear. We defined the genomic binding sites of PML::RARA in primary mouse and human hematopoietic progenitor cells with V5-tagged PML::RARA, using anti-V5-PML::RARA chromatin immunoprecipitation sequencing and CUT&RUN approaches. Most genomic PML::RARA binding sites were found in regions that were already chromatin-accessible (defined by ATAC-seq) in unmanipulated, wild-type promyelocytes, suggesting that these regions are "open" prior to PML::RARA expression. We found that GATA binding motifs, and the direct binding of the chromatin "pioneering factor" GATA2, were significantly enriched near PML::RARA binding sites. Proximity labeling studies revealed that PML::RARA interacts with ~250 proteins in primary mouse hematopoietic cells; GATA2 and 33 others require PML::RARA binding to DNA for the interaction to occur, suggesting that binding to their cognate DNA target motifs may stabilize their interactions. In the absence of PML::RARA, Gata2 overexpression induces many of the same epigenetic and transcriptional changes as PML::RARA. These findings suggested that PML::RARA may indirectly initiate its transcriptional program by activating Gata2 expression: Indeed, we demonstrated that inactivation of Gata2 prior to PML::RARA expression prevented its ability to induce self-renewal. These data suggested that GATA2 binding creates accessible chromatin regions enriched for both GATA and Retinoic Acid Receptor Element motifs, where GATA2 and PML::RARA can potentially bind and interact with each other. In turn, PML::RARA binding to DNA promotes a feed-forward transcriptional program by positively regulating Gata2 expression. Gata2 may therefore be required for PML::RARA to establish its transcriptional program.

IL-1ß expression in bone marrow dendritic cells is induced by TLR2 agonists and regulates HSC function.

Hematopoietic stem/progenitor cells (HSPCs) reside in localized microenvironments, or niches, in the bone marrow that provide key signals regulating their activity. A fundamental property of hematopoiesis is the ability to respond to environmental cues such as inflammation. How these cues are transmitted to HSPCs within hematopoietic niches is not well established. Here, we show that perivascular bone marrow dendritic cells (DCs) express a high basal level of Toll-like receptor-1 (TLR1) and TLR2. Systemic treatment with a TLR1/2 agonist induces HSPC expansion and mobilization. It also induces marked alterations in the bone marrow microenvironment, including a decrease in osteoblast activity and sinusoidal endothelial cell numbers. TLR1/2 agonist treatment of mice in which Myd88 is deleted specifically in DCs using Zbtb46-Cre show that the TLR1/2-induced expansion of multipotent HPSCs, but not HSPC mobilization or alterations in the bone marrow microenvironment, is dependent on TLR1/2 signaling in DCs. Interleukin-1ß (IL-1ß) is constitutively expressed in both murine and human DCs and is further induced after TLR1/2 stimulation. Systemic TLR1/2 agonist treatment of Il1r1-/- mice show that TLR1/2-induced HSPC expansion is dependent on IL-1ß signaling. Single-cell RNA-sequencing of low-risk myelodysplastic syndrome bone marrow revealed that IL1B and TLR1 expression is increased in DCs. Collectively, these data suggest a model in which TLR1/2 stimulation of DCs induces secretion of IL-1ß and other inflammatory cytokines into the perivascular niche, which in turn, regulates multipotent HSPCs. Increased DC TLR1/2 signaling may contribute to altered HSPC function in myelodysplastic syndrome by increasing local IL-1ß expression.

Proteomic and phosphoproteomic landscapes of acute myeloid leukemia.

We have developed a deep-scale proteome and phosphoproteome database from 44 representative acute myeloid leukemia (AML) patients from the LAML TCGA dataset and 6 healthy bone marrow-derived controls. After confirming data quality, we orthogonally validated several previously undescribed features of AML revealed by the proteomic data. We identified examples of posttranscriptionally regulated proteins both globally (ie, in all AML samples) and also in patients with recurrent AML driver mutations. For example, samples with IDH1/2 mutations displayed elevated levels of the 2-oxoglutarate-dependent histone demethylases KDM4A/B/C, despite no changes in messenger RNA levels for these genes; we confirmed this finding in vitro. In samples with NPMc mutations, we identified several nuclear importins with posttranscriptionally increased protein abundance and showed that they interact with NPMc but not wild-type NPM1. We identified 2 cell surface proteins (CD180 and MRC1/CD206) expressed on AML blasts of many patients (but not healthy CD34+ stem/progenitor cells) that could represent novel targets for immunologic therapies and confirmed these targets via flow cytometry. Finally, we detected nearly 30 000 phosphosites in these samples; globally, AML samples were associated with the abnormal phosphorylation of specific residues in PTPN11, STAT3, AKT1, and PRKCD. FLT3-TKD samples were associated with increased phosphorylation of activating tyrosines on the cytoplasmic Src-family tyrosine kinases FGR and HCK and related signaling proteins. PML-RARA-initiated AML samples displayed a unique phosphorylation signature, and TP53-mutant samples showed abundant phosphorylation of serine-183 on TP53 itself. This publicly available database will serve as a foundation for further investigations of protein dysregulation in AML pathogenesis.

Dnmt3a deficiency in the skin causes focal, canonical DNA hypomethylation and a cellular proliferation phenotype.

DNA hypomethylation is a feature of epidermal cells from aged and sun-exposed skin, but the mechanisms responsible for this methylation loss are not known. Dnmt3a is the dominant de novo DNA methyltransferase in the skin; while epidermal Dnmt3a deficiency creates a premalignant state in which keratinocytes are more easily transformed by topical mutagens, the conditions responsible for this increased susceptibility to transformation are not well understood. Using whole genome bisulfite sequencing, we identified a focal, canonical DNA hypomethylation phenotype in the epidermal cells of Dnmt3a-deficient mice. Single-cell transcriptomic analysis revealed an increased proportion of cells with a proliferative gene expression signature, while other populations in the skin were relatively unchanged. Although total DNMT3A deficiency has not been described in human disease states, rare patients with an overgrowth syndrome associated with behavioral abnormalities and an increased risk of cancer often have heterozygous, germline mutations in DNMT3A that reduce its function (Tatton-Brown Rahman syndrome [TBRS]). We evaluated the DNA methylation phenotype of the skin from a TBRS patient with a germline DNMT3AR882H mutation, which encodes a dominant-negative protein that reduces its methyltransferase function by ∼80%. We detected a focal, canonical hypomethylation phenotype that revealed considerable overlap with hypomethylated regions found in Dnmt3a-deficient mouse skin. Together, these data suggest that DNMT3A loss creates a premalignant epigenetic state associated with a hyperproliferative phenotype in the skin and further suggest that DNMT3A acts as a tumor suppressor in the skin.

Immunosuppression and outcomes in adult patients with de novo acute myeloid leukemia with normal karyotypes.

Acute myeloid leukemia (AML) patients rarely have long first remissions (LFRs; >5 y) after standard-of-care chemotherapy, unless classified as favorable risk at presentation. Identification of the mechanisms responsible for long vs. more typical, standard remissions may help to define prognostic determinants for chemotherapy responses. Using exome sequencing, RNA-sequencing, and functional immunologic studies, we characterized 28 normal karyotype (NK)-AML patients with >5 y first remissions after chemotherapy (LFRs) and compared them to a well-matched group of 31 NK-AML patients who relapsed within 2 y (standard first remissions [SFRs]). Our combined analyses indicated that genetic-risk profiling at presentation (as defined by European LeukemiaNet [ELN] 2017 criteria) was not sufficient to explain the outcomes of many SFR cases. Single-cell RNA-sequencing studies of 15 AML samples showed that SFR AML cells differentially expressed many genes associated with immune suppression. The bone marrow of SFR cases had significantly fewer CD4+ Th1 cells; these T cells expressed an exhaustion signature and were resistant to activation by T cell receptor stimulation in the presence of autologous AML cells. T cell activation could be restored by removing the AML cells or blocking the inhibitory major histocompatibility complex class II receptor, LAG3. Most LFR cases did not display these features, suggesting that their AML cells were not as immunosuppressive. These findings were confirmed and extended in an independent set of 50 AML cases representing all ELN 2017 risk groups. AML cell-mediated suppression of CD4+ T cell activation at presentation is strongly associated with unfavorable outcomes in AML patients treated with standard chemotherapy.

Tumor suppressor function of Gata2 in acute promyelocytic leukemia.

Most patients with acute promyelocytic leukemia (APL) can be cured with combined all-trans retinoic acid (ATRA) and arsenic trioxide therapy, which induces the destruction of PML-RARA, the initiating fusion protein for this disease. However, the underlying mechanisms by which PML-RARA initiates and maintains APL cells are still not clear. Therefore, we identified genes that are dysregulated by PML-RARA in mouse and human APL cells and prioritized GATA2 for functional studies because it is highly expressed in preleukemic cells expressing PML-RARA, its high expression persists in transformed APL cells, and spontaneous somatic mutations of GATA2 occur during APL progression in mice and humans. These and other findings suggested that GATA2 may be upregulated to thwart the proliferative signal generated by PML-RARA and that its inactivation by mutation (and/or epigenetic silencing) may accelerate disease progression in APL and other forms of acute myeloid leukemia (AML). Indeed, biallelic knockout of Gata2 with CRISPR/Cas9-mediated gene editing increased the serial replating efficiency of PML-RARA-expressing myeloid progenitors (as well as progenitors expressing RUNX1-RUNX1T1, or deficient for Cebpa), increased mouse APL penetrance, and decreased latency. Restoration of Gata2 expression suppressed PML-RARA-driven aberrant self-renewal and leukemogenesis. Conversely, addback of a mutant GATA2R362G protein associated with APL and AML minimally suppressed PML-RARA-induced aberrant self-renewal, suggesting that it is a loss-of-function mutation. These studies reveal a potential role for Gata2 as a tumor suppressor in AML and suggest that restoration of its function (when inactivated) may provide benefit for AML patients.

Remethylation of Dnmt3a-/- hematopoietic cells is associated with partial correction of gene dysregulation and reduced myeloid skewing.

Mutations in the DNA methyltransferase 3A (DNMT3A) gene are the most common cause of age-related clonal hematopoiesis (ARCH) in older individuals, and are among the most common initiating events for acute myeloid leukemia (AML). The most frequent DNMT3A mutation in AML patients (R882H) encodes a dominant-negative protein that reduces methyltransferase activity by ∼80% in cells with heterozygous mutations, causing a focal, canonical DNA hypomethylation phenotype; this phenotype is partially recapitulated in murine Dnmt3a-/- bone marrow cells. To determine whether the hypomethylation phenotype of Dnmt3a-/- hematopoietic cells is reversible, we developed an inducible transgene to restore expression of DNMT3A in transplanted bone marrow cells from Dnmt3a-/- mice. Partial remethylation was detected within 1 wk, but near-complete remethylation required 6 mo. Remethylation was accurate, dynamic, and highly ordered, suggesting that differentially methylated regions have unique properties that may be relevant for their functions. Importantly, 22 wk of DNMT3A addback partially corrected dysregulated gene expression, and mitigated the expansion of myeloid cells. These data show that restoring DNMT3A expression can alter the epigenetic "state" created by loss of Dnmt3a activity; this genetic proof-of-concept experiment suggests that this approach could be relevant for patients with ARCH or AML caused by loss-of-function DNMT3A mutations.

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.

Rapid and accurate remethylation of DNA in Dnmt3a-deficient hematopoietic cells with restoration of DNMT3A activity.

Here, we characterize the DNA methylation phenotypes of bone marrow cells from mice with hematopoietic deficiency of Dnmt3a or Dnmt3b (or both enzymes) or expressing the dominant-negative Dnmt3aR878H mutation [R882H in humans; the most common DNMT3A mutation found in acute myeloid leukemia (AML)]. Using these cells as substrates, we defined DNA remethylation after overexpressing wild-type (WT) DNMT3A1, DNMT3B1, DNMT3B3 (an inactive splice isoform of DNMT3B), or DNMT3L (a catalytically inactive "chaperone" for DNMT3A and DNMT3B in early embryogenesis). Overexpression of DNMT3A for 2 weeks reverses the hypomethylation phenotype of Dnmt3a-deficient cells or cells expressing the R878H mutation. Overexpression of DNMT3L (which is minimally expressed in AML cells) also corrects the hypomethylation phenotype of Dnmt3aR878H/+ marrow, probably by augmenting the activity of WT DNMT3A encoded by the residual WT allele. DNMT3L reactivation may represent a previously unidentified approach for restoring DNMT3A activity in hematopoietic cells with reduced DNMT3A function.

Genomic landscape of TP53 -mutated myeloid malignancies.

TP53 -mutated myeloid malignancies are most frequently associated with complex cytogenetics. The presence of complex and extensive structural variants complicates detailed genomic analysis by conventional clinical techniques. We performed whole genome sequencing of 42 AML/MDS cases with paired normal tissue to characterize the genomic landscape of TP53 -mutated myeloid malignancies. The vast majority of cases had multi-hit involvement at the TP53 genetic locus (94%), as well as aneuploidy and chromothripsis. Chromosomal patterns of aneuploidy differed significantly from TP53 -mutated cancers arising in other tissues. Recurrent structural variants affected regions that include ETV6 on chr12p, RUNX1 on chr21, and NF1 on chr17q. Most notably for ETV6 , transcript expression was low in cases of TP53 -mutated myeloid malignancies both with and without structural rearrangements involving chromosome 12p. Telomeric content is increased in TP53 -mutated AML/MDS compared other AML subtypes, and telomeric content was detected adjacent to interstitial regions of chromosomes. The genomic landscape of TP53 -mutated myeloid malignancies reveals recurrent structural variants affecting key hematopoietic transcription factors and telomeric repeats that are generally not detected by panel sequencing or conventional cytogenetic analyses. Key Points: WGS comprehensively determines TP53 mutation status, resulting in the reclassification of 12% of cases from mono-allelic to multi-hit Chromothripsis is more frequent than previously appreciated, with a preference for specific chromosomes ETV6 is deleted in 45% of cases, with evidence for epigenetic suppression in non-deleted cases NF1 is mutated in 48% of cases, with multi-hit mutations in 17% of these cases TP53 -mutated AML/MDS is associated with altered telomere content compared with other AMLs.

Genomic landscape of TP53-mutated myeloid malignancies.

TP53-mutated myeloid malignancies are associated with complex cytogenetics and extensive structural variants, which complicates detailed genomic analysis by conventional clinical techniques. We performed whole-genome sequencing (WGS) of 42 acute myeloid leukemia (AML)/myelodysplastic syndromes (MDS) cases with paired normal tissue to better characterize the genomic landscape of TP53-mutated AML/MDS. WGS accurately determines TP53 allele status, a key prognostic factor, resulting in the reclassification of 12% of cases from monoallelic to multihit. Although aneuploidy and chromothripsis are shared with most TP53-mutated cancers, the specific chromosome abnormalities are distinct to each cancer type, suggesting a dependence on the tissue of origin. ETV6 expression is reduced in nearly all cases of TP53-mutated AML/MDS, either through gene deletion or presumed epigenetic silencing. Within the AML cohort, mutations of NF1 are highly enriched, with deletions of 1 copy of NF1 present in 45% of cases and biallelic mutations in 17%. Telomere content is increased in TP53-mutated AMLs compared with other AML subtypes, and abnormal telomeric sequences were detected in the interstitial regions of chromosomes. These data highlight the unique features of TP53-mutated myeloid malignancies, including the high frequency of chromothripsis and structural variation, the frequent involvement of unique genes (including NF1 and ETV6) as cooperating events, and evidence for altered telomere maintenance.

Focal disruption of DNA methylation dynamics at enhancers in IDH-mutant AML cells.

Recurrent mutations in IDH1 or IDH2 in acute myeloid leukemia (AML) are associated with increased DNA methylation, but the genome-wide patterns of this hypermethylation phenotype have not been comprehensively studied in AML samples. We analyzed whole-genome bisulfite sequencing data from 15 primary AML samples with IDH1 or IDH2 mutations, which identified ~4000 focal regions that were uniquely hypermethylated in IDHmut samples vs. normal CD34+ cells and other AMLs. These regions had modest hypermethylation in AMLs with biallelic TET2 mutations, and levels of 5-hydroxymethylation that were diminished in IDH and TET-mutant samples, indicating that this hypermethylation results from inhibition of TET-mediated demethylation. Focal hypermethylation in IDHmut AMLs occurred at regions with low methylation in CD34+ cells, implying that DNA methylation and demethylation are active at these loci. AML samples containing IDH and DNMT3AR882 mutations were significantly less hypermethylated, suggesting that IDHmut-associated hypermethylation is mediated by DNMT3A. IDHmut-specific hypermethylation was highly enriched for enhancers that form direct interactions with genes involved in normal hematopoiesis and AML, including MYC and ETV6. These results suggest that focal hypermethylation in IDH-mutant AML occurs by altering the balance between DNA methylation and demethylation, and that disruption of these pathways at enhancers may contribute to AML pathogenesis.

Kdm6a deficiency restricted to mouse hematopoietic cells causes an age- and sex-dependent myelodysplastic syndrome-like phenotype.

Kdm6a/Utx, a gene on the X chromosome, encodes a histone H3K27me3 demethylase that has an orthologue on the Y chromosome (Uty) (Zheng et al. 2018). We previously identified inactivating mutations of Kdm6a in approximately 50% of mouse acute promyelocytic leukemia samples; however, somatic mutations of KDM6A are more rare in human AML samples, ranging in frequency from 2-15% in different series of patients, where their role in pathogenesis is not yet clear. In this study, we show that female Kdm6aflox/flox mice (with allele inactivation initiated by Vav1-Cre in hematopoietic stem and progenitor cells (HSPCs) have a sex-specific phenotype that emerges with aging, with features resembling a myelodysplastic syndrome (MDS). Female Kdm6a-knockout (KO) mice have an age-dependent expansion of their HSPCs with aberrant self-renewal, but they did not differentiate normally into downstream progeny. These mice became mildly anemic and thrombocytopenic, but did not develop overt leukemia, or die from these cytopenias. ChIP-seq and ATAC-seq studies showed only minor changes in H3K27me3, H3K27ac, H3K4me, H3K4me3 and chromatin accessibility between Kdm6a-WT and Kdm6a-KO mice. Utilizing scRNA-seq, Kdm6a loss was linked to the transcriptional repression of genes that mediate hematopoietic cell fate determination. These data demonstrate that Kdm6a plays an important role in normal hematopoiesis, and that its inactivation may contribute to AML pathogenesis.

Functional and epigenetic phenotypes of humans and mice with DNMT3A Overgrowth Syndrome.

Germline pathogenic variants in DNMT3A were recently described in patients with overgrowth, obesity, behavioral, and learning difficulties (DNMT3A Overgrowth Syndrome/DOS). Somatic mutations in the DNMT3A gene are also the most common cause of clonal hematopoiesis, and can initiate acute myeloid leukemia (AML). Using whole genome bisulfite sequencing, we studied DNA methylation in peripheral blood cells of 11 DOS patients and found a focal, canonical hypomethylation phenotype, which is most severe with the dominant negative DNMT3AR882H mutation. A germline mouse model expressing the homologous Dnmt3aR878H mutation phenocopies most aspects of the human DOS syndrome, including the methylation phenotype and an increased incidence of spontaneous hematopoietic malignancies, suggesting that all aspects of this syndrome are caused by this mutation.

Cellular stressors contribute to the expansion of hematopoietic clones of varying leukemic potential.

Hematopoietic clones harboring specific mutations may expand over time. However, it remains unclear how different cellular stressors influence this expansion. Here we characterize clonal hematopoiesis after two different cellular stressors: cytotoxic therapy and hematopoietic transplantation. Cytotoxic therapy results in the expansion of clones carrying mutations in DNA damage response genes, including TP53 and PPM1D. Analyses of sorted populations show that these clones are typically multilineage and myeloid-biased. Following autologous transplantation, most clones persist with stable chimerism. However, DNMT3A mutant clones often expand, while PPM1D mutant clones often decrease in size. To assess the leukemic potential of these expanded clones, we genotyped 134 t-AML/t-MDS samples. Mutations in non-TP53 DNA damage response genes are infrequent in t-AML/t-MDS despite several being commonly identified after cytotoxic therapy. These data suggest that different hematopoietic stressors promote the expansion of distinct long-lived clones, carrying specific mutations, whose leukemic potential depends partially on the mutations they harbor.

Haploinsufficiency for DNA methyltransferase 3A predisposes hematopoietic cells to myeloid malignancies.

The gene that encodes de novo DNA methyltransferase 3A (DNMT3A) is frequently mutated in acute myeloid leukemia genomes. Point mutations at position R882 have been shown to cause a dominant negative loss of DNMT3A methylation activity, but 15% of DNMT3A mutations are predicted to produce truncated proteins that could either have dominant negative activities or cause loss of function and haploinsufficiency. Here, we demonstrate that 3 of these mutants produce truncated, inactive proteins that do not dimerize with WT DNMT3A, strongly supporting the haploinsufficiency hypothesis. We therefore evaluated hematopoiesis in mice heterozygous for a constitutive null Dnmt3a mutation. With no other manipulations, Dnmt3a+/- mice developed myeloid skewing over time, and their hematopoietic stem/progenitor cells exhibited a long-term competitive transplantation advantage. Dnmt3a+/- mice also spontaneously developed transplantable myeloid malignancies after a long latent period, and 3 of 12 tumors tested had cooperating mutations in the Ras/MAPK pathway. The residual Dnmt3a allele was neither mutated nor downregulated in these tumors. The bone marrow cells of Dnmt3a+/- mice had a subtle but statistically significant DNA hypomethylation phenotype that was not associated with gene dysregulation. These data demonstrate that haploinsufficiency for Dnmt3a alters hematopoiesis and predisposes mice (and probably humans) to myeloid malignancies by a mechanism that is not yet clear.

Comprehensive discovery of noncoding RNAs in acute myeloid leukemia cell transcriptomes.

To detect diverse and novel RNA species comprehensively, we compared deep small RNA and RNA sequencing (RNA-seq) methods applied to a primary acute myeloid leukemia (AML) sample. We were able to discover previously unannotated small RNAs using deep sequencing of a library method using broader insert size selection. We analyzed the long noncoding RNA (lncRNA) landscape in AML by comparing deep sequencing from multiple RNA-seq library construction methods for the sample that we studied and then integrating RNA-seq data from 179 AML cases. This identified lncRNAs that are completely novel, differentially expressed, and associated with specific AML subtypes. Our study revealed the complexity of the noncoding RNA transcriptome through a combined strategy of strand-specific small RNA and total RNA-seq. This dataset will serve as an invaluable resource for future RNA-based analyses.

Genetic heterogeneity of induced pluripotent stem cells: results from 24 clones derived from a single C57BL/6 mouse.

Induced pluripotent stem cells (iPSCs) have tremendous potential as a tool for disease modeling, drug testing, and other applications. Since the generation of iPSCs "captures" the genetic history of the individual cell that was reprogrammed, iPSC clones (even those derived from the same individual) would be expected to demonstrate genetic heterogeneity. To assess the degree of genetic heterogeneity, and to determine whether some cells are more genetically "fit" for reprogramming, we performed exome sequencing on 24 mouse iPSC clones derived from skin fibroblasts obtained from two different sites of the same 8-week-old C57BL/6J male mouse. While no differences in the coding regions were detected in the two parental fibroblast pools, each clone had a unique genetic signature with a wide range of heterogeneity observed among the individual clones: a total of 383 iPSC variants were validated for the 24 clones (mean 16.0/clone, range 0-45). Since these variants were all present in the vast majority of the cells in each clone (variant allele frequencies of 40-60% for heterozygous variants), they most likely preexisted in the individual cells that were reprogrammed, rather than being acquired during reprogramming or cell passaging. We then tested whether this genetic heterogeneity had functional consequences for hematopoietic development by generating hematopoietic progenitors in vitro and enumerating colony forming units (CFUs). While there was a range of hematopoietic potentials among the 24 clones, only one clone failed to differentiate into hematopoietic cells; however, it was able to form a teratoma, proving its pluripotent nature. Further, no specific association was found between the mutational spectrum and the hematopoietic potential of each iPSC clone. These data clearly highlight the genetic heterogeneity present within individual fibroblasts that is captured by iPSC generation, and suggest that most of the changes are random, and functionally benign.