Key Points
IDH2R172 mutations define a unique subgroup with distinct TFH-like gene expression signatures in AITL.
IDH2R172 mutations can induce DNA and repressive histone hypermethylation in AITL.
Abstract
Angioimmunoblastic T-cell lymphoma (AITL) is a common subtype of peripheral T-cell lymphoma (PTCL) with a poor prognosis. We performed targeted resequencing on 92 cases of PTCL and identified frequent mutations affecting RHOA, TET2, DNMT3A, and isocitrate dehydrogenase 2 (IDH2). Although IDH2 mutations are largely confined to AITL, mutations of the other 3 can be found in other types of PTCL, although at lower frequencies. These findings indicate a key role of epigenetic regulation in the pathogenesis of AITL. However, the epigenetic alterations induced by these mutations and their role in AITL pathogenesis are still largely unknown. We correlated mutational status with gene expression and global DNA methylation changes in AITL. Strikingly, AITL cases with IDH2R172 mutations demonstrated a distinct gene expression signature characterized by downregulation of genes associated with TH1 differentiation (eg, STAT1 and IFNG) and a striking enrichment of an interleukin 12-induced gene signature. Ectopic expression of IDH2R172K in the Jurkat cell line and CD4+ T cells led to markedly increased levels of 2-hydroxyglutarate, histone-3 lysine methylation, and 5-methylcytosine and a decrease of 5-hydroxymethylcytosine. Correspondingly, clinical samples with IDH2 mutations displayed a prominent increase in H3K27me3 and DNA hypermethylation of gene promoters. Integrative analysis of gene expression and promoter methylation revealed recurrently hypermethylated genes involved in T-cell receptor signaling and T-cell differentiation that likely contribute to lymphomagenesis in AITL.
Introduction
Peripheral T-cell lymphoma (PTCL) is a heterogeneous group of generally aggressive lymphoid malignancies, accounting for 10% to 15% of all non-Hodgkin lymphomas.1 Angioimmunoblastic T-cell lymphoma (AITL) represents 20% to 25% of all PTCLs and is recognized as a distinct entity.1,2 Currently used chemotherapy is rarely curative with a 5-year overall survival of <30%.1,3,4 Gene expression profiling (GEP) and pathologic analysis suggest that the follicular helper T-cell (TFH cell) is the cell of origin for AITL.5-7
We and others have identified frequent mutations affecting isocitrate dehydrogenase 2 (IDH2) at arginine-172 (R172), ten-eleven translocation 2 (TET2), DNA methyltransferase 3A (DNMT3A), CD28, and RHOA in AITL.8-12 IDH2, TET2, and DNMT3A mutations may all affect epigenetic modifications. TET2 and DNMT3A mutations seem to occur at an early stage of hematopoietic cell differentiation, as these mutations have been found in nonmalignant hematopoietic cells of several PTCL cases and in normal elderly individuals.11-14 Loss-of-function mutations of DNMT3A have been observed in several hematologic malignancies, leading to DNA hypomethylation.15 Mitochondrial IDH2 protein and cytosolic IDH1 protein are homologs. Surprisingly, unlike other malignancies,16,17 AITL harbors mutations only at IDH2R172, but not at IDH2R140 or in IDH1. In the tricarboxylic acid cycle, IDH3 catalyzes the oxidative decarboxylation of isocitrate to α-ketoglutarate (α-KG). Under anabolic conditions, IDH2 primarily catalyzes the reverse reaction: reductive carboxylation of α-KG to isocitrate. The IDH2R172 mutation confers a neomorphic activity; namely, conversion of α-KG to the R-enantiomer of 2-hydroxyglutarate (R-2-HG),18 which, in normal cells, is maintained at low concentrations by 2-HG dehydrogenase. R-2-HG acts as a competitive antagonist of the α-KG–dependent dioxygenases,19 which notably include the TET family of 5-methylcytosine hydroxylases and jumonji-C domain-containing histone demethylases.19-21 Thus, mutant IDH2R172 is predicted to impair DNA and histone demethylation and leads to abnormal regulation of gene transcription, which may promote lymphomagenesis.
The epigenetic alterations induced by IDH2, TET2, and DNMT3A mutations and their roles in AITL pathogenesis are still largely unknown. The biochemical and functional consequences of IDH2R172 mutations in T cells have not yet been elucidated. In this study, we performed targeted resequencing of these 3 epigenetic regulators in molecularly defined PTCL cases. We also analyzed the biochemical changes associated with IDH2R172 mutations and investigated alterations in GEP, DNA methylation, and histone modification to improve our understanding of their roles in the pathogenesis of AITL.
Materials and methods
Patient samples, T-cell line, and normal CD4+ T cells
We included 90 molecularly defined PTCL cases with GEP and clinical outcome data, as well as 2 pathologically diagnosed cases (Table 1).22 This study was approved by the Institutional Review Board of the University of Nebraska Medical Center and City of Hope Medical Center. Jurkat cells (Clone E6-1, ATCC TIB-152) were cultured in 10 mM HEPES-buffered RPMI 1640 (Invitrogen) supplemented with 10% fetal bovine serum, penicillin G (100 U/mL), and streptomycin (100 μg/mL) at 37°C in a 5% CO2 atmosphere. Normal human CD4+ T cells were isolated from peripheral blood lymphocytes of healthy donors (Miltenyi Biotec) and were cultured in medium as mentioned above with interleukin (IL)-2 (30 U/mL; R&D Systems) and anti-CD3/CD28 Dynabeads (Invitrogen).
Sample type . | Targeted resequencing of TET2 and DNMT3A . | Sanger sequencing of IDH2* . | GEP . | RRBS . |
---|---|---|---|---|
Tumor specimens | ||||
AITL | 39 | 58 | 36 | 18 |
PTCL-NOS | 41 | 17 | NA | NA |
ALK(−) ALCL | 12 | NA | NA | NA |
Primary CD4+T cells | NA | NA | 5† | 4‡ |
Tonsils | NA | NA | NA | 3 |
Sample type . | Targeted resequencing of TET2 and DNMT3A . | Sanger sequencing of IDH2* . | GEP . | RRBS . |
---|---|---|---|---|
Tumor specimens | ||||
AITL | 39 | 58 | 36 | 18 |
PTCL-NOS | 41 | 17 | NA | NA |
ALK(−) ALCL | 12 | NA | NA | NA |
Primary CD4+T cells | NA | NA | 5† | 4‡ |
Tonsils | NA | NA | NA | 3 |
NA, not applicable.
Low tumor content cases were confirmed by pyrosequencing; 39 of 58 AITL cases have corresponding targeted resequencing data.
Treated with IL-12 for 0, 2, 8, 24, or 48 hours.
Two sets were freshly isolated from donor PBMCs, and the other 2 were cultured with anti-CD3/CD28 beads for 2 weeks before DNA extraction.
Targeted resequencing of TET2, DNMT3A, and RHOA and mutational analysis of IDH2
We used TruSeq Custom Amplicon design (Illumina) to interrogate the entire coding regions of TET2, DNMT3A, and RHOA in genomic DNA. Twenty-five cycles of amplification were used for 250 ng genomic DNA per sample. We used the Mutascope algorithm for sequence analysis.23 Mutational analysis of the IDH2 gene was performed by bi-directional Sanger sequencing8 and pyrosequecing (supplemental Methods, available on the Blood Web site) using primers that specifically amplify R140 or R172 regions.
GEP and pathway analysis
GEP was performed and reported in a previous study22 using HG-U133-plus 2.0 arrays (Affymetrix). BRB Array Tools (http://linus.nci.nih.gov/BRB-ArrayTools.html) were used for unsupervised hierarchical clustering and for determining significant differences in gene expression among predefined groups: IDH2/TET2 double-mutant, TET2 single-mutant, and wild-type AITL using normalized log2-transformed signals. We used Ingenuity pathway analysis (IPA; Qiagen) to perform functional annotation. The Gene Set Enrichment Analysis (GSEA) algorithm (http://www.broad.mit.edu/gsea/)24 was used to identify a leading edge set of the genes that best distinguished IDH2/TET2 double-mutant and TET2 single-mutant AITL.
Reduced representation bisulfite sequencing for methylation analysis
Reduced representation bisulfite sequencing (RRBS) was performed as previously described25 with some modifications. Alignment rates and coverage are shown in supplemental Figure 1. Full details are given in the supplemental Methods.
Histone trimethylation analysis by immunohistochemistry
We performed immunohistochemical staining on tissue microarray sections of AITL with 1:100 rabbit polyclonal anti-H3K27me3 (07-449; Millipore) and 1:400 rabbit polyclonal anti-H3K9me3 (ab8898; Abcam) as previously described.26 The stained sections were scanned by Ventana’s iScan Coreo slide scanner (Ventana Medical Systems) using a 40 × 0.25 objective and exported into ImageScope analysis software (Leica Microsystems) for quantification. The signal threshold for nuclei in the negative control was determined and applied to all scanned images. Areas above the signal threshold were automatically measured and divided into weak-positive, positive, and strong-positive based on the signal intensities.
Protein histone extraction and western blotting
Whole cell lysates were prepared as previously described,27 and histones were isolated using a histone extraction kit (ab113476; Abcam). Protein was separated by 10% sodium dodecyl sulfate-polyacrylamide gel electrophoresis and transferred to polyvinylidene difluoride membrane (Bio-Rad).27 After incubation with primary antibody (listed in supplemental Table 3) overnight at 4°C, the membrane was washed and incubated with anti-mouse immunoglobulin (Ig)G horseradish peroxidase (HRP)-linked secondary antibody (7076S; Cell Signaling) at 1:10 000 or anti-rabbit IgG HRP-linked secondary antibody (7074S; Cell Signaling) at 1:10 000 according to the primary antibody used. A Western enhanced chemiluminescence substrate (170-5060; Bio-Rad) was used for detection of HRP.
Liquid chromatography coupled to tandem mass spectrometry analysis of 2-HG, 5-methylcytosine, and 5-hydroxymethyl cytosine in genomic DNA
2-HG was quantified using reverse-phase liquid chromatography coupled to tandem mass spectrometry (LC-MS/MS) as previously described.28 DNA hydrolysis was performed as described previously29 with slight modifications. 5-Methylcytosine (5mC) and 5-hydroxymethyl cytosine (5hmC) levels were determined by mass spectrometry, and details are described in the supplemental Methods.
Statistical analyses
Details are given in the supplemental Methods.
Results
Mutations affecting epigenetic regulators are more frequent in AITL than other subtypes of PTCL
Consistent with our previous report,8 mutations in AITL were observed in IDH2 at R172 but not at R140, unlike in other cancer types. The overall frequency of IDH2R172 mutations in AITL was 32.8% (19 of 58), including cases from our previous cohort8 (Figure 1A). In contrast, the IDH2R172 mutation occurred in only 1 of 24 PTCL, not otherwise specified (PTCL-NOS) cases, which was molecularly classified into the TBX21 subgroup by GEP analysis.26 The most common mutant variants were arginine to serine or lysine, similar to the IDH2 mutational spectrum in acute myelogenous leukemia (AML). There was no difference in IDH2 mRNA expression levels between the IDH2R172 mutant and wild-type cases (supplemental Figure 2A). Unlike the findings in AML, where IDH2R172 mutations are associated with a poor prognosis,30,31 no association of IDH2R172 mutations with overall survival was observed (supplemental Figure 3A).
TET2 mutations were found in 82.1% (32 of 39) of AITL, 46.3% (19 of 41) of PTCL-NOS, and 33.3% (4 of 12) of ALK(−) anaplastic large cell lymphoma (ALCL) (Figure 1A; supplemental Table 1). In PTCL-NOS, both the TBX21 (57%, 10 of 18) and GATA3 subgroups (41.7%, 5 of 12) harbored TET2 mutations. The TET2 mutations were distributed throughout the entire coding region; however, the majority of mutations were observed in the catalytic domains, including the Cys-rich and double-stranded β-helix (DSBH) domains,32 with a few recurrently targeted codons (Figure 1B). We observed DNMT3A mutations at similar frequencies in AITL (38.5%, 15 of 39) and PTCL-NOS (36.6%, 15 of 41), but less frequently in ALK(−) ALCL (16.7%, 2 of 12). In PTCL-NOS, 33.3% (6 of 18) harbored DNMT3A mutations in the TBX21 subgroup and 25% (3 of 12) in the GATA3 subgroup (Figure 1A). Most DNMT3A mutations were clustered in the MTase-binding domain (Figure 1B), presumably leading to loss of function. Of note, the most frequent mutation site was residue Arg-882 (R882) located in the MTase-binding domain. This DNMT3AR882 mutant has been reported to inhibit wild-type DNMT3A activity by disrupting homotetramerization.33
The majority (26 of 29) of RHOA mutations observed in AITL occurred at the same codon, resulting in a Gly17Val substitution in the GTP-binding domain. Consistent with previous reports,11,34 RHOA mutations were very common in AITL (28 of 39), also present in PTCL-NOS (11 of 41), but absent in ALCL (0 of 12) (Figure 1A), and tended to co-occur with TET2 mutations (2-tailed Fisher’s exact test, P = .012) and IDH2 mutations (2-tailed Fisher’s exact test, P = .003). All IDH2/TET2 double-mutant AITL cases had a RHOA mutation.
Interestingly, IDH2 and TET2 mutations were not mutually exclusive, unlike in AML or glioblastoma, but showed significant (P < .001) co-occurrence in AITL (Figure 1A). Most IDH2 (68%) and DNMT3A (80%) mutant cases also had TET2 mutations.
IDH2R172 mutations define an AITL subgroup with a distinct gene expression signature
To gain an overview of the functional consequences of the mutations affecting epigenetic regulators, we first performed unsupervised hierarchical clustering of 37 AITL samples with GEP data and known mutational status of IDH2, TET2, and DNMT3A. Interestingly, most IDH2R172 mutant cases formed a single cluster in unsupervised hierarchical clustering analysis (Figure 2A). We used Fisher's exact test to determine the significance of the tendency of neighboring cases on the dendrogram to have identical mutation status. The fraction of IDH2R172 mutant cases that had mutant neighbors was significantly higher than the fraction of IDH2 wild-type cases that had mutant neighbors (P = .011). In contrast, neither TET2 mutant cases (P = .073) nor DNMT3A mutant cases (P = .727) showed significant clustering (Figure 2A), indicating that IDH2R172 mutations are associated with a unique gene expression profile in AITL.
To identify the genes that are differentially expressed between IDH2 wild-type and mutant groups, we performed supervised hierarchical clustering using a univariate permutation test for the significance of individual genes (P = .001). Among the 344 differentially expressed genes, ∼40% (136 of 344) were downregulated in IDH2R172 mutant cases compared with IDH2 wild-type cases, an observation distinct from AML (supplemental Figure 4A).35 Only a small subset of genes was downregulated in IDH2 mutant cases in both AITL and AML (supplemental Figure 4B), suggesting different functional consequences of IDH2 mutations in these 2 malignancies.
Because mutant IDH2 generates 2-HG and consequently can inhibit TET2 enzyme activity, IDH2 mutation is not expected to co-occur with TET2 mutation; nevertheless, most IDH2R172 mutant cases in AITL also harbor the TET2 mutation (30 of 36). Therefore, we performed class comparison among IDH2/TET2 double-mutant, TET2 single-mutant, and wild-type AITL cases to identify the differences among these groups and found 1035 differentially expressed genes. The tumor cell content of TET2 single-mutant AITL and IDH2/TET2 double-mutant AITL was similar (supplemental Figure 5). Wild-type cases and TET2 single-mutant cases are similar, but markedly different from IDH2 /TET2 double-mutant cases. Double-mutant cases showed upregulation of TFH-associated genes (IL21 and ICOS) and downregulation of genes associated with TH1 (IL2, STAT1, CXCR3), TH2 (IL10RA), and TH17 (IL17RA), indicating a more TFH cell-like phenotype (Figure 2B; supplemental Figure 6). To gain further biological insight, we performed IPA and identified several important clusters of genes associated with proliferation, apoptosis, cell cycle, and T-cell development and differentiation. Pathways involved in T helper cell differentiation, IL-17 signaling (IL17RA, AKT3, CXCL10, and CCL2), and interferon (IFN) signaling (IFNG, STAT1, IFITM3, IRF1, and IFIT3) were dysregulated in IDH2/TET2 double-mutant cases (Figure 3A; supplemental Figure 7). Strikingly, IPA upstream analysis predicted activation of the vascular endothelial growth factor pathway and inhibition of IFN signaling in the double-mutant group (Figure 3B). This observation was further supported by GSEA using the Broad Institute and lymphoid signature databases, which showed significant enrichment of IFN response genes (P < .01) in the TET2 single-mutant group. However, the double-mutant cases showed enrichment of EZH2 targets, MYC targets, FOXP3 targets, downregulated FOXO3 target genes, and genes associated with proliferation and cell cycle transitions or in apoptosis due to CDKN1 via TP53 (supplemental Table 2; Figure 3C).
We also observed a striking correlation of GEP of IDH2/TET2 double-mutant cases with that of primary CD4+ T cells treated with IL-12 (Figure 2B). IL-12 is a major driver of TFH cell differentiation, and it can induce naïve CD4+ T cells to differentiate into TFH-like cells in vitro.36 We analyzed 6245 genes that were highly expressed (top quartile) in primary CD4+ T cells either after 48-hour IL-12 treatment and/or no treatment. Among these 6245 genes, the differences in GEP between individual double-mutant cases and the median tumor value strongly correlated with the changes in GEP induced in normal CD4+ T-cells by IL-12 (Figure 3D).
To investigate the similarity between double-mutant AITL and primary TFH cells, we used transcriptome sequencing data of TFH cells (n = 3), naïve CD4+ T cells (n = 4), and non-TFH effector cells (CD4+CD45RA−TCRβ+PD-1loCXCR5loPSGL-1hi) (n = 3) from a published dataset.37 By identifying genes with significantly different expression (P ≤ .05, 1-sided Student t test) between TFH cells and both naïve CD4+ T cells and non-TFH effector cells, we defined 1045 “TFH-up genes” and 699 “TFH-down genes.” Similar numbers of “TFH-up genes” were upregulated (comparing individual cases to the median expression level) in all 3 mutation groups, but significantly more genes were downregulated in IDH2/TET2 double-mutant AITL than wild-type and TET2 single-mutant cases (Figure 3E). This suggests that mutant IDH2 is associated with the repression of non-TFH differentiation genes, presumably through epigenetic mechanisms.
IDH2R172 mutations are associated with DNA hypermethylation in AITL and T cells
To investigate the effect of IDH2R172 mutations on the T-cell methylome, we transduced Jurkat cells and primary CD4+ T cells with vectors expressing IDH2R172K or wild-type IDH2. As anticipated, ectopic expression of mutant IDH2 resulted in increased intracellular 2-HG in Jurkat cells (see Figure 7A). Remarkably, a global increase in 5-methylcytosine (5mC) and a global decrease in 5-hydroxymethyl cytosine (5hmC) were observed in the DNA of Jurkat cells expressing the IDH2R172K mutant compared with cells expressing wild-type IDH2 or empty vector (Figure 4A-B), indicating that IDH2R172K mutation leads to epigenetic changes in DNA through 2-HG. Similar results for 2-HG, 5mC, and 5hmC were observed in primary CD4+ T cells (Figure 4C-E).
Because an AITL-specific DNA methylation pattern is not known, we conducted RRBS to investigate the genome-wide alterations in DNA methylation in 18 AITL cases with relatively high tumor content and compared them with T-cell samples isolated from peripheral blood of four donors. Most of AITL cases with high tumor content harbor TET2 mutations (17 of 18), and 10 of 18 cases harbor IDH2R172 mutations. Because AITL tumors have a considerable number and variety of nonmalignant cells, we included 3 tonsils to serve as additional controls to filter out differences resulting from cell type-specific methylation of nonmalignant cells.
Studies in several cancers have shown that gain-of-function mutations of IDH2 are associated with a DNA hypermethylation phenotype,18 possibly through inhibition of TET2. Compared with normal tissues, IDH2 mutant AITL showed global alterations of the DNA methylation landscape including hypermethylation of regulatory regions, such as promoters and CpG islands (CGIs) (Figure 4F) that are known to be correlated with downregulation of gene transcription, whereas gene body and intergenic regions are hypomethylated.38 CpG island shores (CGSs) are ±2-kb regions flanking CGI, and CGI and CGS are normally unmethylated for most genes. We classified promoters into CGI-associated promoters, CGS-associated promoters, and non–CGI/CGS-associated promoters. We observed global hypermethylation of CGI-associated promoters and global hypomethylation of non–CGI/CGS-associated promoters in both IDH2 mutant and wild-type AITL (Figure 4G). We observed that IDH2 wild-type AITL had hypermethylation of many genes that are normally not expressed in CD4+ T cells. However, among genes normally expressed in CD4+ T cells, a significantly higher number were hypermethylated in the IDH2 mutant group (Figure 4H), indicating a more significant functional effect of IDH2R172 mutations in AITL. Additionally, using RRBS and GEP data, we identified genes most significantly hypermethylated in promoter regions and downregulated in IDH2 mutant compared with IDH2 wild-type AITL, including a subset of helper T-cell differentiation genes (SMAD7 and SMURF2) (Table 2).
Gene . | EntrezID . | Downregulation (IDH2 WT > IDH2 mutant) . | Hypermethylation (IDH2 mutant > normal tissue) . | Correlation between gene expression and methylation . |
---|---|---|---|---|
ABL1 | 25 | P = .042 | P = .017 | R = −0.308 |
AFF1 | 4299 | P = .009 | P = .109 | R = −0.559 |
AMIGO2 | 347902 | P = .010 | P = .029 | R = −0.704 |
ANKRD53 | 79998 | P = .006 | P = .065 | R = −0.499 |
ANTXR2 | 118429 | P = .017 | P = .035 | R = −0.462 |
APOBR | 55911 | P = .016 | P = .053 | R = −0.449 |
ATHL1 | 80162 | P = .009 | P = .063 | R = −0.661 |
ATP8B2 | 57198 | P = .008 | P = .007 | R = −0.301 |
CD7 | 924 | P = .019 | P = .037 | R = −0.559 |
CIDEB | 27141 | P = .019 | P = .067 | R = −0.311 |
CPLX1 | 10815 | P = .048 | P = .001 | R = −0.320 |
DGKQ | 1609 | P = .004 | P = .079 | R = −0.427 |
EVC2 | 132884 | P = .037 | P = .045 | R = −0.432 |
FAM149A | 25854 | P = .035 | P = .024 | R = −0.342 |
FXYD5 | 53827 | P = .046 | P = .002 | R = −0.400 |
FZD6 | 8323 | P = .049 | P = .005 | R = −0.411 |
GIMAP8 | 155038 | P = .008 | P = .021 | R = −0.392 |
GPR68 | 8111 | P = .049 | P = .064 | R = −0.439 |
HDAC11 | 79885 | P = .037 | P = .010 | R = −0.462 |
HKR1 | 284459 | P = .024 | P = .067 | R = −0.533 |
IL10RA | 3587 | P = .001 | P = .014 | R = −0.451 |
ISYNA1 | 51477 | P = .018 | P = .025 | R = −0.602 |
KBTBD11 | 9920 | P = .024 | P = .002 | R = −0.402 |
KCNJ4 | 3761 | P = .032 | P < .001 | R = −0.490 |
KLF3 | 51274 | P = .008 | P = .018 | R = −0.634 |
KLHL22 | 84861 | P = .004 | P = .008 | R = −0.534 |
LPAR2 | 9170 | P = .005 | P = .034 | R = −0.427 |
NOVA2 | 4858 | P = .025 | P = .013 | R = −0.341 |
NTN1 | 9423 | P = .047 | P = .266 | R = −0.309 |
PDLIM2 | 64236 | P = .001 | P = .005 | R = −0.385 |
PTGDR | 5729 | P = .044 | P = .016 | R = −0.320 |
PTPN3 | 5774 | P = .042 | P < .001 | R = −0.321 |
SATB1 | 6304 | P = .013 | P = .087 | R = −0.555 |
SECTM1 | 6398 | P = .002 | P = .083 | R = −0.444 |
SLC12A7 | 10723 | P = .009 | P = .006 | R = −0.572 |
SLC1A6 | 6511 | P = .043 | P = .117 | R = −0.426 |
SLC35D2 | 11046 | P = .002 | P = .019 | R = −0.380 |
SMAD7 | 4092 | P = .025 | P = .009 | R = −0.335 |
SMO | 6608 | P = .015 | P = .004 | R = −0.300 |
SMURF2 | 64750 | P = .006 | P = .003 | R = −0.378 |
SPNS2 | 124976 | P = .045 | P = .018 | R = −0.338 |
TCEA3 | 6920 | P = .004 | P = .029 | R = −0.462 |
XKR8 | 55113 | P = .012 | P < .001 | R = −0.666 |
ZBTB4 | 57659 | P = .001 | P = .001 | R = −0.774 |
ZFP3 | 124961 | P = .012 | P = .003 | R = −0.357 |
ZNF704 | 619279 | P = .024 | P = .059 | R = −0.363 |
Gene . | EntrezID . | Downregulation (IDH2 WT > IDH2 mutant) . | Hypermethylation (IDH2 mutant > normal tissue) . | Correlation between gene expression and methylation . |
---|---|---|---|---|
ABL1 | 25 | P = .042 | P = .017 | R = −0.308 |
AFF1 | 4299 | P = .009 | P = .109 | R = −0.559 |
AMIGO2 | 347902 | P = .010 | P = .029 | R = −0.704 |
ANKRD53 | 79998 | P = .006 | P = .065 | R = −0.499 |
ANTXR2 | 118429 | P = .017 | P = .035 | R = −0.462 |
APOBR | 55911 | P = .016 | P = .053 | R = −0.449 |
ATHL1 | 80162 | P = .009 | P = .063 | R = −0.661 |
ATP8B2 | 57198 | P = .008 | P = .007 | R = −0.301 |
CD7 | 924 | P = .019 | P = .037 | R = −0.559 |
CIDEB | 27141 | P = .019 | P = .067 | R = −0.311 |
CPLX1 | 10815 | P = .048 | P = .001 | R = −0.320 |
DGKQ | 1609 | P = .004 | P = .079 | R = −0.427 |
EVC2 | 132884 | P = .037 | P = .045 | R = −0.432 |
FAM149A | 25854 | P = .035 | P = .024 | R = −0.342 |
FXYD5 | 53827 | P = .046 | P = .002 | R = −0.400 |
FZD6 | 8323 | P = .049 | P = .005 | R = −0.411 |
GIMAP8 | 155038 | P = .008 | P = .021 | R = −0.392 |
GPR68 | 8111 | P = .049 | P = .064 | R = −0.439 |
HDAC11 | 79885 | P = .037 | P = .010 | R = −0.462 |
HKR1 | 284459 | P = .024 | P = .067 | R = −0.533 |
IL10RA | 3587 | P = .001 | P = .014 | R = −0.451 |
ISYNA1 | 51477 | P = .018 | P = .025 | R = −0.602 |
KBTBD11 | 9920 | P = .024 | P = .002 | R = −0.402 |
KCNJ4 | 3761 | P = .032 | P < .001 | R = −0.490 |
KLF3 | 51274 | P = .008 | P = .018 | R = −0.634 |
KLHL22 | 84861 | P = .004 | P = .008 | R = −0.534 |
LPAR2 | 9170 | P = .005 | P = .034 | R = −0.427 |
NOVA2 | 4858 | P = .025 | P = .013 | R = −0.341 |
NTN1 | 9423 | P = .047 | P = .266 | R = −0.309 |
PDLIM2 | 64236 | P = .001 | P = .005 | R = −0.385 |
PTGDR | 5729 | P = .044 | P = .016 | R = −0.320 |
PTPN3 | 5774 | P = .042 | P < .001 | R = −0.321 |
SATB1 | 6304 | P = .013 | P = .087 | R = −0.555 |
SECTM1 | 6398 | P = .002 | P = .083 | R = −0.444 |
SLC12A7 | 10723 | P = .009 | P = .006 | R = −0.572 |
SLC1A6 | 6511 | P = .043 | P = .117 | R = −0.426 |
SLC35D2 | 11046 | P = .002 | P = .019 | R = −0.380 |
SMAD7 | 4092 | P = .025 | P = .009 | R = −0.335 |
SMO | 6608 | P = .015 | P = .004 | R = −0.300 |
SMURF2 | 64750 | P = .006 | P = .003 | R = −0.378 |
SPNS2 | 124976 | P = .045 | P = .018 | R = −0.338 |
TCEA3 | 6920 | P = .004 | P = .029 | R = −0.462 |
XKR8 | 55113 | P = .012 | P < .001 | R = −0.666 |
ZBTB4 | 57659 | P = .001 | P = .001 | R = −0.774 |
ZFP3 | 124961 | P = .012 | P = .003 | R = −0.357 |
ZNF704 | 619279 | P = .024 | P = .059 | R = −0.363 |
To identify genes whose expression is likely downregulated by promoter methylation, we selected differentially methylated promoters in AITL cases compared with normal tissues. A total of 46 hypermethylated genes and 9 hypomethylated genes in AITL were identified having strong correlation (R < −0.4) with transcriptional changes (Figure 5). The epigenetically downregulated genes included negative regulators of T-cell receptor signaling (PTPN7, SIT1, and DGKA), positive regulators of TH1 activity (MATK and PHF11), and other cell surface marker genes (CD7 and IL10RA) (supplemental Figure 8).
One of the most striking findings from GSEA analysis of the genes differentially expressed among wild-type, TET2 single-mutant, and IDH2/TET2 double-mutant groups was the enrichment of genes associated with the repressive mark histone 3 trimethyllysine-27 (H3K27me3) (Figure 6A).We analyzed genes with both H3K4me3 and H3K27me3 in human embryonic stem cells (hESCs) (poised genes) as these genes are particularly prone to repression by H3K27me3. Indeed, hESC-poised genes expressed in CD4+ T cells are downregulated in AITL with IDH2R172 mutations (Figure 6B). Initially H3K27me3-repressed genes are prone to cancer-specific promoter DNA hypermethylation, definitively silencing them.39 Thus, we compared DNA methylation between hESC-poised genes and non–hESC-poised genes and found that poised genes are significantly more likely to be hypermethylated than nonpoised genes (Figure 6C).
IDH2R172 mutations are associated with enhanced H3K27me3 in AITL and T cells
To determine whether the mutant IDH2R172 can alter histone lysine methylation in vitro, we ectopically expressed mutant IDH2R172K, wild-type IDH2, and empty vector in Jurkat cells (Figure 7A). H3K27me3, H3K9me3, and H3K4me3 were significantly increased in Jurkat cells expressing IDH2R172K (Figure 7B). Immunohistochemical staining revealed significantly enhanced H3K27me3 in the IDH2R172 mutant AITL cases group (n = 3; average 42.2% positive cells) compared with cases in the IDH2 wild-type group (n = 3; average 16.2% positive cells; Figure 7C). Strongly positive areas averaged only 2.9% in the IDH2 wild-type cases, but constituted 18.8% in the IDH2R172 mutant cases (Figure 7D). Overall, these data demonstrate increased H3K27me3 both in vitro and in tumor samples in agreement with our GEP findings and strongly suggest that IDH2R172K mutations can dysregulate the posttranslational modification of histones in T cells.
Discussion
The aberrant DNA and histone methylation identified in various malignancies suggests an important role for these epigenetic changes in tumorigenesis, although the profiles vary among tumors, indicating a possible dependency on cellular context.40 Mutations of epigenetic regulators (IDH2R172, TET2, and DNMT3A) are common in AITL, emphasizing the importance of aberrant epigenetic alterations in lymphomagenesis. Given the presence of TET2 and DNMT3A mutations in other subtypes of PTCL, but not IDH2R172 mutations, it is highly likely that IDH2R172 mutations play a unique role in the pathogenesis of AITL. Indeed, unsupervised hierarchical clustering identified a distinct GEP associated with IDH2R172 mutations in AITL. The changes in GEP in IDH2R172 mutant cases likely resulted in part from promoter hypermethylation, but also from disruption of histone lysine demethylation.
One of the most interesting characteristics in the IDH2/TET2 double-mutant group is enrichment of the TFH-like phenotype and downregulation of features of other types of helper T-cells. IL-6, IL-12, and IL-23 have been shown to induce IL-21 expression in human naïve CD4+ T cells. IL-12 not only can induce IL-21 production but also can induce other TFH phenotypes, including upregulation of CXCR5, ICOS, and BCL6 and the capacity to facilitate B-cell differentiation into antibody-secreting cells.41 Treatment of normal human CD4+ T cells by IL-12 for 48 hours induced a steady increase of IL-21 and CXCR5 (supplemental Figure 9). The GEP of IDH2/TET2 double-mutant AITL was highly enriched in the signature of CD4+ T cells stimulated by IL-12, suggesting a more polarized TFH phenotype. In addition, IDH2/TET2 double-mutant AITL showed similarity of gene expression with primary TFH cells. These observations suggest that limiting the plasticity of differentiation to other TH lineages is an important consequence of IDH2R172 mutations. Of the signaling pathways attributable to differentiation of helper T-cells, IFN signaling (STAT1, IRF1, and IFIT3) and IL-17 signaling (IL17RA, CXCL10, and CCL2) were downregulated in the IDH2/TET2 double-mutant group. Conversely, transforming growth factor-β signaling, which contributes to the differentiation of TFH cells,42 was enriched. In keeping with this observation, we identified a subset of helper T-cell differentiation genes (SMAD7 and SMURF2) that were hypermethylated in patients with an IDH2R172 mutation (Table 2).
Unlike in AML, IDH2 and TET2 mutations are not mutually exclusive in AITL. Despite the significant differences between IDH2/TET2 double-mutant AITL and normal tissues, and between TET2 single-mutant AITL and normal tissues, there is only a moderate difference in genome-wide methylation between IDH2/TET2 double-mutant and TET2 single-mutant AITL, which seems insufficient to explain the much greater differences on GEP. However, GEP includes secondary changes from genetic and epigenetic alterations. The significant changes in gene expression in the IDH2R172 mutant group may also represent consequences of histone modification.
Mutant IDH2 inhibits various histone demethylases. Nearly all histone demethylation is brought about by the Jumonji family of histone demethylases, with 27 distinct enzymes in humans. H3K4, H3K9, H3K27, and H3K79 methylation has been shown to be elevated in glioma cells expressing mutant IDH1/2, both in vitro and in vivo.43 Consistent with this finding, we observed similar alterations in vitro and significantly elevated levels of H3K27me3 in AITL cases with IDH2R172 mutation, suggesting that 2-HG produced by IDH2 mutants indeed interferes with histone demethylases. In ESCs, bivalent chromatin modification is thought to silence genes controlling cell lineage through the repressive histone mark H3K27me3 while rending them poised for subsequent activation during differentiation through H3K4me3. Histone lysine demethylases (KDM) 6A/B are responsible for demethylating H3K27me3. Inactivating mutations of KDM6A have frequently been identified in hematopoietic malignancies and solid tumors.44,45 Genes that are poised in ESCs are especially likely to be persistently turned off by H3K27me3 due to inhibition of KDM6A/B by mutant IDH2. In the context of AITL, our results suggest that IDH2R172 mutation further promotes lymphomagenesis by altering gene expression through histone lysine methylation. This hypothesis is supported by the observation that hESC-poised genes were significantly downregulated in the IDH2R172 mutant group. Aberrant repression by H3K27me3 may be induced, albeit less readily, in other genes that are downregulated in IDH2R172 mutant AITL. In both IDH2R172 mutant and wild-type cases, genes with little or no expression in T cells are especially likely to undergo promoter DNA hypermethylation. Thus, some of the genes aberrantly repressed by H3K27me3 may later become hypermethylated in their promoters, locking in their transcriptional silencing. Thus, loss of function of TET2 and histone demethylases may cooperate in epigenetic dysregulation.
In AML, both R140 and R172 neomorphic mutants of IDH2 were identified. In contrast, AITL harbors mutations only at IDH2R172, but not at IDH1R132 or IDH2R140. This exclusivity is highly interesting, but the reason is currently not clear. Studies by others have shown that mutations differ quantitatively in 2-HG production33,46 ; this raises the possibility that different tumor types may favor certain ranges of 2-HG concentrations. In previous studies,33,46 IDH2R140 was shown to be a weaker 2-HG–producing mutant than IDH2R172. In addition, IDH2R172K mutant 10T mesenchymal progenitor cells, but not IDH2R140Q mutant cells, can grow as xenografts,33 indicating different biological consequences resulting from the 2 mutants. In AML, IDH2R140 mutant AML cells frequently have normal cytogenetics and NPM1 mutations, but IDH2R172 is frequently the only mutation detected,47 again suggesting biological differences. Understanding the mechanistic difference between the IDH2R140 and IDH2R172 mutants is certainly important and should be investigated further in the future.
Both IDH2R172 and TET2 mutations can lead to promoter hypermethylation, giving a strong rationale for considering the use of hypomethylating agents in AITL treatment. TET2-mutated myelodysplastic syndrome patients have been shown to response to hypomethylating agents.48 Moreover, a case report using 5-azacytidine to treat a patient with chronic myelomonocytic leukemia, who also had TET2-mutant AITL, led to a complete remission of both the AITL and chronic myelomonocytic leukemia.49 Vascular endothelial growth factor-A has been shown to be overexpressed in malignant cells and endothelial cells from AITL and associated with an adverse prognosis,50 suggesting that it is another potential therapeutic target. A small molecule inhibitor targeting mutant IDH2, AG-221, has also shown promising results in treating IDH2-mutant AML in animal models and in clinical trials,51,52 suggesting that inhibition of mutant IDH2 could be a useful therapeutic strategy that specifically targets an oncogenic mutation in AITL. However, because most IDH2-mutant cases also have TET2 mutations, it is possible that combination with a DNA hypomethylating agent or antiangiogenic drugs could improve efficacy.
Presented in part, at the 56th American Society of Hematology annual meeting, San Francisco, CA, December 6-9, 2014.
The online version of this article contains a data supplement.
The publication costs of this article were defrayed in part by page charge payment. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.
Acknowledgments
The authors acknowledge helpful discussions with Dr Katharine Yen and her help in measuring the 2-HG levels at Agios Pharmaceuticals (Cambridge, MA). The authors are grateful for Rhishikesh Thakare’s and Dr Yazen Alnouti’s help in measuring global 5hmC and 5mC and Victoria B. Smith’s and Dr Charles A Kuszynski’s technical assistance with fluorescence activated cell sorting at University of Nebraska Medical Center. The authors thank Dr David L. Klinkebiel and Dr James Eudy for constructive discussion, Dr Francoise Berger for contribution to sample collection and characterization, and Dr Xin Wang for support and suggestions.
This work was supported in part by the Lymphoma Research Foundation (F-263549, to J.I.), the translational research program of the Leukemia and Lymphoma Society (6129-14, to J.I.), the University of Nebraska Medical Center Clinical-Translational Research Scholars Program, the National Institutes of Health, National Cancer Institute Eppley Cancer Center support grant P30CA036727, City of Hope Cancer Center support grant P30 CA033572, and Specialized Programs of Research Excellence (1 P50 CA136411-01 01A1 PP-4, to W.C.C.). C.W. and B.J. were recipients of Chinese Scholarship Council funding.
Authorship
Contribution: C.W., T.W.M., J.I., and W.C.C. conceived and designed the project; C.W., A.B., J.W., Z.M., B.J., J.R., A.C., C.S., Y.L., and E.G.R. performed experiments; C.W., Q.G., W.Z., X.W., S.H., L.S., and T.W.M. analyzed the data; A.R., R.D.G., K.F., D.D.W., T.C.G., G.O., L.M.S., J.V., and W.C.C. provided clinical samples, reviewed pathology results and immunostaining data, or contributed to the clinical data; and C.W., T.W.M., J.I., and W.C.C. wrote and finalized the paper.
Conflict-of-interest disclosure: The authors declare no competing financial interests.
Correspondence: Javeed Iqbal, Department of Pathology and Microbiology, Nebraska Medical Center, Omaha, NE 68198; e-mail: jiqbal@unmc.edu; or Wing C. (John) Chan, Department of Pathology, City of Hope Medical Center, 1500 East Duarte Rd, Duarte, CA 91010; e-mail: jochan@coh.org.