|Home | About | Journals | Submit | Contact Us | Français|
Users may view, print, copy, download and text and data- mine the content in such documents, for the purposes of academic research, subject always to the full Conditions of use: http://www.nature.com/authors/editorial_policies/license.html#terms
TET2 is a close relative of TET1, an enzyme that converts 5-methylcytosine (5-mC) to 5-hydroxymethylcytosine (5-hmC) in DNA1,2. The gene encoding TET2 resides at chromosome 4q24, in a region showing recurrent microdeletions and copy-neutral loss of heterozygosity (CN-LOH) in patients with diverse myeloid malignancies3. Somatic TET2 mutations are frequently observed in myelodysplastic syndromes (MDS), myeloproliferative neoplasms (MPN), MDS/MPN overlap syndromes including chronic myelomonocytic leukemia (CMML), acute myeloid leukemias (AML) and secondary AML (sAML)4–12. We show here that TET2 mutations associated with myeloid malignancies compromise TET2 catalytic activity. Bone marrow samples from patients with TET2 mutations displayed uniformly low levels of 5-hmC in genomic DNA compared to bone marrow samples from healthy controls. Moreover, small hairpin RNA (shRNA)-mediated depletion of Tet2 in mouse haematopoietic precursors skewed their differentiation towards monocyte/macrophage lineages in culture. There was no significant difference in DNA methylation between bone marrow samples from patients with high 5-hmC versus healthy controls, but samples from patients with low 5-hmC showed hypomethylation relative to controls at the majority of differentially-methylated CpG sites. Our results demonstrate that TET2 is important for normal myelopoiesis, and suggest that disruption of TET2 enzymatic activity favours myeloid tumorigenesis. Measurement of 5-hmC levels in myeloid malignancies may prove valuable as a diagnostic and prognostic tool, to tailor therapies and assess responses to anti-cancer drugs.
We transiently transfected HEK293T cells with Myc-tagged murine Tet2 and assessed 5-mC and 5-hmC levels by immunocytochemistry (Fig. 1, Suppl. Figs. 1–4). Myc-Tet2-expressing cells displayed a strong increase in 5-hmC staining and a concomitant decrease in 5-mC staining in the nucleus (Fig. 1b, c, quantified in Suppl. Fig. 4). In contrast, 5-hmC was undetectable or barely detected in nuclei of cells expressing mutant Tet2 with H1302Y, D1304A substitutions in the signature HxD motif1,12,17 involved in coordinating Fe2+, and there was no obvious decrease in nuclear 5-mC staining (Fig. 1b, c, Suppl. Fig. 4). These studies confirm13 that Tet2 is a catalytically active enzyme that converts 5-mC to 5-hmC in genomic DNA.
Mutations in TET2 residues H1881 and R1896, predicted to bind Fe2+ and 2OG, respectively, have been identified repeatedly in patients with myeloid malignancies4,5,7,10. HEK293T cells expressing Tet2 mutants H1802R and H1802Q (Fig. 1a, Suppl. Fig. 2) showed greatly diminished 5-hmC staining and no loss of 5-mC staining, consistent with participation of this residue in catalysis (Fig. 1b, c, Suppl. Fig. 4a, b). We analysed missense mutations identified in TET2 in our own (Suppl. Table S1) and other3–6, 11 studies (P1367S, W1291R, G1913D, E1318G and I1873T). HEK293T cells expressing Tet2 mutants P1287S, W1211R or C1834D (Suppl. Figs. 2, 3a) displayed low 5-hmC staining and strong 5-mC staining (Suppl. Figs. 3b, 3c, 4c, 4d), suggesting a role for these residues in the integrity of the catalytic or DNA binding domains. Cells expressing Tet2 R1817S/M (Fig. 1a, Suppl. Figs. 2, 3a) were positive for 5-hmC staining but changes in 5-mC staining could not be reliably assessed (Figs. 1b,c, Suppl. Figs. 3b, 3c, 4).
To quantify these findings, we developed dot blot assays to detect 5-hmC in genomic DNA (Suppl. Fig. 5). In the first assay format, the blot was developed with a specific antiserum to 5-hmC (Suppl. Fig. 5b, left), whose ability to recognize 5-hmC depended strongly on the density of 5-hmC in DNA (Suppl. Fig. 5c, top). We therefore developed a more sensitive and quantitative assay in which DNA was treated with bisulfite to convert 5-hmC to cytosine 5-methylenesulfonate (CMS) 14 (Suppl. Fig. 5a), after which CMS was measured with a specific anti-CMS antiserum (Suppl. Fig. 5b, right). Unlike anti-5-hmC which reacted efficiently only with DNA containing high densities of 5-hmC, the anti-CMS antiserum recognized DNA with an average of only a single 5-hmC per 201 bp (Suppl. Fig. 5c, bottom). This lack of density dependence allowed us to plot the signal obtained with 2-fold dilutions of a standard oligonucleotide containing a known amount of 5-hmC against the amount of CMS obtained after bisulfite conversion. We assumed 100% conversion efficiency15 and used the linear portion of the standard curve to compute the amount of CMS, and therefore 5-hmC, in the DNA samples (e.g. see Fig. 2a, right).
To assess 5-hmC levels, we obtained uniform populations of Tet2-expressing HEK293T cells by transfection with Tet2-IRES-CD25 plasmid followed by magnetic isolation of CD25-expressing cells1. Wild type and mutant Tet2 proteins were expressed at comparable levels (Fig. 1d, Suppl. Fig. 3d). Anti-5-hmC/ CMS dot blots of genomic DNA revealed, as expected, that 5-hmC was barely detectable in DNA from cells transfected with empty vector; DNA from cells expressing wild type Tet2 showed a substantial increase in 5-hmC and a corresponding decrease in 5-mC; and DNA from cells expressing the HxD mutant Tet2 protein had very low 5-hmC (Fig. 1e, Suppl. Figs. 3e, 6). DNA from cells expressing 7 of the 9 mutant Tet2 proteins tested -- H1802Q/R, R1817S/M, W1211R, P1287S and C1834D -- contained significantly less 5-hmC than DNA from cells expressing wild type Tet2 (Fig. 1e, Suppl. Figs. 3e, 6), confirming our previous conclusion that these mutations impair enzymatic activity.
We measured 5-hmC (CMS) levels in genomic DNA extracted from bone marrow or blood (with >20% immature myeloid cells) of 88 patients with myeloid malignancies and 17 healthy controls (Suppl. Table S1). In blinded experiments, DNA was treated with bisulfite and CMS levels were evaluated. TET2 mutations were strongly associated with low genomic 5-hmC (Fig. 2, Suppl. Fig. 7a). To confirm these conclusions in a statistically rigorous fashion, we tested samples for which a sufficient amount of DNA was available to make independent dilutions in triplicate, so that a median and standard deviation for 5-hmC (CMS) levels in each patient could be derived (Suppl. Fig. 7b). Analysis of DNA from 9 healthy donors and 41 patients (28 with wild type TET2 and 13 with TET2 mutations, Suppl. Table S1) revealed a strong, statistically significant correlation of TET2 mutations with low 5-hmC (Fig. 2c). In contrast, samples from patients with wild type TET2 showed a bimodal distribution, with 5-hmC levels ranging from ~0.4 to ~3.8 pmol/µg DNA (Fig. 2c; Suppl. Fig. 7; also see Fig. 4).
We examined Tet2 expression in haematopoietic cell subsets isolated from bone marrow and thymus of C57BL/6 mice (Suppl. Figs. 8, 9). Tet2 mRNA was highly expressed in lineage-negative (Lin−) Sca-1+c-Kithi multipotent progenitors (LSK), at levels similar to those in embryonic stem cells (ESC). Expression was maintained at high levels in myeloid progenitors (common myeloid progenitors, CMPs, and granulocyte-monocyte progenitors, GMPs), was low in mature granulocytes (Gr-1+Mac-1+) and high in monocytes (Gr-1− Mac-1+) (Suppl. Fig. 9a, middle panel).
To test the role of Tet2 in myelopoiesis, we transduced bone marrow stem/progenitor cells with Tet2 shRNA (Suppl. Fig. 10a), effectively downregulating Tet2 mRNA and protein relative to control cells transduced with empty vector or scrambled shRNA (Fig. 3a, b) (refer to Suppl Fig. 10 b for choice of Tet2 shRNA). Tet2 depletion promoted expansion of Mac-1+ F4/80+ and Mac1+ CD115+ monocyte/macrophage cells in the presence of G-CSF or GM-CSF, cytokines that support granulocyte and granulocyte/ monocyte development respectively, but not in the presence of M-CSF, which promotes growth of monocytic progenitors (Fig. 3c; Suppl. Fig. 10d). Simultaneous treatment with GM-CSF and M-CSF, or GM-CSF and G-CSF, also led to increased numbers of monocyte/macrophage cells (Fig. 3c). These results indicate that Tet2 plays an important role in normal myelopoiesis. However, Tet2 does not markedly influence short-term proliferation of myeloid-lineage cells: when shRNA-transduced Lin− cells were cultured in the presence of GM-CSF and pulse-labeled with BrdU, Tet2 depletion promoted monocyte/ macrophage expansion but CD115+ (M-CSFR+) cells from the two cultures showed no difference in acute BrdU incorporation (Suppl. Fig. 11).
We asked whether 5-hmC levels in tumour samples correlated with DNA methylation status. A histogram of normalised values from 88 patients and 17 healthy individuals showed the expected bimodal distribution (see Online Methods): healthy controls and most patient samples with wild type TET2 had high 5-hmC, while the majority of patient samples with mutant TET2 had low 5-hmC (Fig. 4b). The DNA methylation status of 62 samples was interrogated at 27,578 CpG sites. As expected16, the resulting histograms were strikingly bimodal, with sites within and outside CpG islands showing low and predominantly high methylation respectively (Fig. 4c). Comparison of 28 control samples with 24 high 5-hmC tumour samples (22 TET2 wild type, 2 TET2 mutant) showed no significant difference in DNA methylation; in contrast comparison of the control samples with 29 low 5-hmC tumour samples (7 TET2 wild type, 22 TET2 mutant) yielded 2512 differentially methylated sites, of which the majority (2510 sites) were hypomethylated compared to controls (Fig. 4d, Suppl. Table S2 online). Thus TET2 loss-of-function is predominantly associated with decreased methylation at CpG sites.
To summarise, our studies demonstrate a strong correlation between myeloid malignancies and loss of TET2 catalytic activity. The leukemia-associated missense mutations associated with diminished 5-hmC levels provide clues to the structure of the TET2 catalytic domain. The W1211R, P1287S and C1834D mutations affect positions that are highly conserved within the catalytic domain of the TET subfamily of dioxygenases2: W1211 is located at the beginning of the strand just N-terminal to the core of the DSBH, and is predicted to constitute part of the “mouth” of the active site pocket of the enzyme; P1287 is predicted to stabilize the conformation of the junction between the N-terminal helix and the first core strand of the DSBH; and G1913/ C1834 is predicted to be the N-terminal capping residue of a helix that lines the “mouth” of the DSBH and potentially interacts with substrate DNA2. The E1238G mutation had no detectable effect on 5-hmC production in our overexpression assays, however, the patient with this mutation also showed CN-LOH spanning 4q24, a feature that likely contributes to the significant reduction in 5-hmC levels observed in the bone marrow.
Low 5-hmC levels were observed in a subset of patients with apparently wild type TET2, whose clinical phenotypes resembled those of patients with mutant TET2. In several of these patients, TET2 mRNA expression was not significantly different from controls; mutations in other TET proteins have not been described (Suppl. Text). Some patients in the wild type TET2/ low 5-hmC category may harbour mutations in regulatory or partner proteins for TET2, or in cis-regulatory regions controlling TET2 mRNA expression. Alternatively, the primary event in some of these patients may be CpG hypomethylation, resulting in decreased 5-hmC secondary to depletion of the substrate, 5-mC.
There is little consensus on whether TET2 mutations correlate with clinical outcome. One study reported an association with decreased survival in AML4, whereas others report little prognostic value in MPN diseases7,10,12. Assays for 5-hmC may increase our options for the molecular classification of myeloid malignancies, making it possible to ask whether patients with high or low levels of genomic 5-hmC show differences in disease progression or therapeutic response. Notably, histone deacetylase and DNA methyltransferase inhibitors show clinical efficacy in patients with CMML and AML17; and genomic 5-hmC levels could potentially be a useful prognostic indicator or predictor of patient responses or refractoriness to “epigenetic” therapy with demethylating agents.
DNA methylation is highly aberrant in cancer18–20. Since TET operates on 5-mC, we were surprised to find that TET2 loss-of-function in myeloid tumours was associated with widespread hypomethylation rather than the expected hypermethylation at differentially-methylated CpG sites. Tumour samples with low 5-hmC may have expanded cells with localized hypomethylation at these sites; or TET2 may control DNA methylation indirectly, for instance by regulating the expression or recruitment of one or more DNA methyltransferases, perhaps via 5-hmC-binding proteins. Alternatively, if TET2 and 5-hmC are required for cells to exit the stem cell state, loss of TET2 function in myeloid neoplasms may reactivate a stem-like state characterized by generalized hypomethylation and consequent genomic instability21,22. Indeed, hypomorphic DNMT1 mutations associated with genome-wide DNA hypomethylation skew haematopoietic differentiation toward myelo-erythroid lineages23, and promote the development of aggressive T cell lymphomas due to activation and insertion of endogenous retroviruses24,25. Further studies of the role of TET2 in haematopoietic differentiation should uncover the relation between TET2 loss-of-function, DNA methylation changes and myeloid neoplasia.
Genomic DNA was extracted from bone marrow/ peripheral blood samples from healthy donors and patients with MDS, MDS/MPN, primary and secondary AMLs. Clinical features and other detailed information pertaining to the patient samples are summarized in Suppl. Table 1.
For CMS detection, genomic DNA was treated with sodium bisulfite using the EpiTect Bisulfite kit (QIAGEN). DNA samples were denatured and 2-fold serial dilutions were spotted on a nitrocellulose membrane in an assembled Bio-Dot apparatus (Bio-Rad). The blotted membrane was washed, air-dried, vacuum-baked, blocked, and incubated with anti-5-hmC or anti-CMS antibody (1:1,000) and HRP-conjugated anti-rabbit IgG secondary antibody. To ensure equal spotting of total DNA on the membrane, the same blot was stained with 0.02% methylene blue in 0.3 M sodium acetate (pH 5.2). To compare results obtained in different experiments, we used the normalisation procedure described in Online Methods (see Figs. 4a, b, which incorporate data from Fig. 2 and Suppl. Fig. 6).
The DNA methylation status of bisulfite-treated genomic DNA was probed at 27,578 CpG dinucleotides using the Illumina® Infinium® 27k array (Illumina, San Diego, CA)26. Methylation status was calculated from the ratio of methylation-specific and demethylation-specific fluorophores (β-value) using BeadStudio® Methylation Module (Illumina, San Diego, CA). We removed sites on the Y and X chromosomes from the analysis because of inconsistent methylation status with respect to gender (a known problem based on communication with Illumina. Calculations are based on β values, which correspond to the methylation status of a site ranging from 0 to 1, returned by Illumina's BeadStudio software. We tested sites for differential methylation using an empirical Bayes approach employing a modified t-test (limma). The false discovery rate (FDR) is controlled at a level of 5% by the Benjamini-Hochberg correction.
This work was supported by NIH grants R01 AI44432 and RC1 DA028422 (to A.R.), NIH grants K24 HL077522 and R01 HL098522, an Established Investigator award from the Aplastic Anemia & MDS Foundation, and an award from the Bob Duggan Memorial Research Fund (to J.P.M), NIH grant R01 HG4069 (to X.S.L.) and a pilot grant from Harvard Catalyst, The Harvard Clinical and Translational Science Center (NIH Grant #1 UL1 RR 025758-02, to S.A.). Y.H. was supported by postdoctoral fellowships from the GlaxoSmithKline-Immune Disease Institute (GSK-IDI) Alliance and the Leukemia and Lymphoma Society of America. H.S.B. is supported by a postdoctoral fellowship from the GSK-IDI Alliance.
Author contributionsM.K. analysed the biochemical effects of patient-associated TET2 mutations and performed the in vitro differentiation studies; Y.H. generated and characterized the anti-CMS antiserum, developed the quantitative dot-blot assay and quantified 5-hmC in DNA samples from patients and healthy controls. A.M.J., R.G. and J.P.M. provided patient and control DNA for 5hmC quantification, performed DNA methylation arrays and analysed TET2 mutational status in patients. U.J.P. and X.S.L. carried out the statistical analysis of 5-hmC levels and methylation data; M.T., H.S.B. and K.P.K. provided critical reagents; J.A. and E.D.L. contributed to molecular cloning and mouse maintenance respectively; and L.A. and S.A. provided essential intellectual input. A.R. set overall goals and coordinated the collaborations.
Competing Financial Interests
The authors declare no competing financial interests.