The molecular etiology of myeloproliferative neoplasms (MPNs) remains incompletely understood, despite recent advances incurred through the discovery of several different mutations in MPN patients. We have recently described overexpression of the transcription factor NF-E2 in MPN patients and shown that elevated NF-E2 levels in vivo cause an MPN phenotype and predispose to leukemic transformation in transgenic mice. We report the presence of acquired insertion and deletion mutations in the NF-E2 gene in MPN patients. These result in truncated NF-E2 proteins that enhance wild-type (WT) NF-E2 function and cause erythrocytosis and thrombocytosis in a murine model. NF-E2 mutant cells acquire a proliferative advantage, witnessed by clonal dominance over WT NF-E2 cells in MPN patients. Our data underscore the role of increased NF-E2 activity in the pathophysiology of MPNs.
Myeloproliferative neoplasms (MPNs) constitute a group of clonal, malignant hematopoietic disorders, which includes polycythemia vera (PV), essential thrombocythemia (ET), and primary myelofibrosis (PMF). Patients present clinically with erythrocytosis and/or thrombocytosis (PV and ET, respectively) or with pan-cytopenia (PMF) and display a propensity to transform to acute leukemia. Despite the discovery of a variety of mutations in MPN patients (Vainchenker et al., 2011), most notably an activating point mutation in the JAK2 kinase (JAK2V617F; Baxter et al., 2005; James et al., 2005; Kralovics et al., 2005), the pathophysiology of these disorders remains incompletely understood. Consequently, at present, no curative treatment exists besides bone marrow transplantation, which in this patient population is associated with considerable morbidity and mortality.
We have shown that expression of the hematopoietic transcription factor nuclear factor erythroid-2 (NF-E2) is elevated in MPN patients irrespective of the presence or absence of the JAK2V617F mutation (Goerttler et al., 2005; Wang et al., 2010). Based on this observation, we have generated a novel murine model for MPN by overexpressing NF-E2 in transgenic mice (Kaufmann et al., 2012). NF-E2–overexpressing mice display many MPN features including thrombocytosis, leukocytosis, Epo-independent colony formation, characteristic bone marrow histology, expansion of stem and progenitor compartments, and spontaneous transformation to acute myeloid leukemia (AML; Kaufmann et al., 2012). During AML transformation, NF-E2 transgenic mice acquire recurrent genetic alterations, particularly a gain of genetic material corresponding to human chromosome 8, which is also frequently observed in MPN patients and is associated with a higher rate of transformation to acute leukemia (Hussein et al., 2009). These data have established a role for aberrant NF-E2 expression in the pathophysiology of MPN. Here we test the hypothesis that alterations in the NF-E2 gene may contribute to the development of MPN.
To further investigate the role of NF-E2 in the pathophysiology of MPN, we sequenced the NF-E2 gene in a cohort of MPN patients (144 PV, 120 ET and 192 MF patients), as well as in patients with MDS (n = 57), CMML (n = 67), secondary, post-MPN AML (n = 39), and healthy controls (n = 65; Table S1). Seven different insertions and deletions leading to frameshift mutations were detected in the NF-E2 coding sequence in eight patients with MPNs, three with PV and five with MF, either PMF or secondary post-MPN MF (Table 1 and Fig. S1). Two mutations, c.782-785delAGAG and c.662_663insG, were found in two different patients each, and one patient harbored two separate mutations (Table 1). The frameshifts introduce premature stop codons in the open reading frame, leading to truncations in the NF-E2 protein (Fig. 1 A). One 12 bp deletion, c.889-900del, causes an in-frame deletion of four amino acids within the leucine zipper heterodimerization domain, causing the two proximal leucine to lose the seven amino acid spacing typical of a leucine zipper (O’Shea et al., 1989).
Insertion and deletion mutations in NF-E2 were detected exclusively in PV and MF patients (3 of 144 patients, 2.1%; and 5 of 192 patients, 2.6%, respectively). These mutations were not observed in patients with ET, MPN-U, MDS, secondary post-MPN AML and CMML, or in healthy controls. Constitutive DNA, obtained from buccal swabs or T cells, was available from five patients with insertions or deletions, and in all cases, we were able to demonstrate that the mutations were acquired (Table 1).
Because the truncated NF-E2 proteins contain neither the DNA binding domain nor the leucine zipper required for dimerization to small Maf proteins or, in one case, contain a deletion in the leucine zipper, we investigated whether NF-E2 mutants retain DNA binding activity in an electrophoretic mobility shift assay (EMSA). As previously demonstrated (Igarashi et al., 1994), NF-E2 requires interaction with a small Maf protein, here MafG, to bind DNA (Fig. 1 B, compare lane 2 and lane 4). Specificity of the NF-E2/MafG heterodimer binding to its cognate DNA sequence was verified by competition and super-shift experiments (Fig. 1 B, lanes 5–8).
In contrast to WT NF-E2, two truncated NF-E2 mutants (p.L245VfsX5, here called 248aa, and p.E261AfsX3, here called 262aa) were unable to bind DNA even in the presence of MafG (Fig. 1 C, compare lanes 1 and 2–3). Furthermore, the NF-E2 mutant carrying the 4 aa deletion, here called Δ297-300, was likewise unable to bind DNA (Fig. 1 C, compare lanes 1 and 4). Protein expression of all mutants was verified by Western blotting (Fig. 1 D). Importantly, analysis of protein extracts from primary cells of MPN patients also demonstrates expression of the truncated, mutant protein in MPN patient cells (Fig. 1 E).
Subsequently, we investigated whether the NF-E2 mutants retained transactivation potential in reporter gene assays. Two different reporter gene constructs were used, one containing a known NF-E2 binding site from the β-globin promoter (Blank et al., 1997), and a second containing an NF-E2 binding site in the β1-tubulin promoter, recently characterized by us (unpublished data). Both reporter constructs yielded similar results (Figs. 2, A and B); all three NF-E2 mutants tested were no longer able to evoke reporter gene activity above the background levels observed with MafG alone.
As transcription off a plasmid DNA does not reflect physiological conditions within the chromatin-bound DNA, we used CB3 cells to test the mutants under physiological conditions. CB3 cells carry a homozygous viral insertion in the NF-E2 locus and therefore express neither NF-E2 nor its target β-globin (Ben-David et al., 1992). Upon reintroduction of NF-E2, β-globin is robustly expressed (Li et al., 2001). We introduced WT NF-E2 or the three NF-E2 mutants tested above into CB3 cells and measured β-globin mRNA expression by qRT-PCR. Although transduction of WT NF-E2 caused a 100-fold increase in β-globin mRNA expression, all three NF-E2 mutants were unable to activate β-globin expression (Fig. 2 C).
Because all NF-E2 mutations were observed in a heterozygous state, hence unfolding their effect in the presence of WT NF-E2, we investigated the effect of coexpression of WT NF-E2 and NF-E2 mutants in CB3 cells. CB3 cells stably transfected with WT NF-E2 were infected with an empty pLeGO-iT2 virus (Weber et al., 2008), or viruses expressing the two truncation mutants, 248aa and 262aa, or the 4aa deletion Δ297-300, and assayed for β-globin expression (Fig. 2, D and E). Although transduction with empty pLeGO-iT2 had no effect on β-globin mRNA expression in WT NF-E2–expressing CB3 cells (Fig. 2 E), transduction with the three mutants, which on their own retain no transactivation activity (Fig. 2 C), enhances WT NF-E2 activity between three- and ninefold, with the Δ297-300 mutant displaying the strongest effect (Fig. 2 E). Thus, although unable to transactivate transcription on their own, the NF-E2 mutants significantly enhance WT NF-E2 activity.
We investigated the molecular mechanism by which mutant NF-E2 enhances the activity of WT NF-E2. NF-E2 mRNA levels were determined by qRT-PCR in five patients carrying NF-E2 mutations, seven MPN patients with WT NF-E2, and three healthy controls (Fig. 3 A). As we have previously published, NF-E2 is overexpressed 2–10-fold in MPN patients expressing WT NF-E2 (Fig. 3 A; Goerttler et al., 2005). However, expression of NF-E2 mRNA is further significantly augmented in patients carrying NF-E2 mutations, reaching a median 2.5-fold of that observed in WT NF-E2 patients (Fig. 3 A). The presence of mutant NF-E2 therefore significantly increases NF-E2 mRNA levels.
The qRT-PCR assay used does not distinguish between RNA encoding WT or mutant, truncated NF-E2. We therefore investigated NF-E2 protein levels in primary cells from MPN patients carrying NF-E2 mutations. Corresponding to the observed increase in NF-E2 mRNA, we detected elevated NF-E2 protein levels in patients expressing NF-E2 mutants (Fig. 1 E). Interestingly, the expression of full-length, WT NF-E2 was strongly increased, whereas the levels of truncated, mutant protein were low. Thus, presence of mutant NF-E2 results in elevated levels of WT NF-E2. These data provide one possible mechanism by which mutant NF-E2 augments the activity of its WT counterpart.
In a second series of experiments, we compared the stability of WT NF-E2 protein in the presence or absence of mutant NF-E2. The CB3 cells expressing WT NF-E2 were cotransfected either with a virus expressing the 262aa NF-E2 truncation mutant or with an empty control virus. The levels of WT and truncated, mutant NF-E2 were determined by Western blotting and shown to be comparable in their relative ratio to those observed in MPN patients; mutant NF-E2 was present at much lower levels than WT NF-E2 (compare Fig. 3 C, right, lane 1 with Fig. 1 E). Cells were treated with 20 µg/ml cycloheximide to inhibit de novo protein biosynthesis and assayed for NF-E2 protein levels by Western blotting for up to 6 h after the start of treatment (Fig. 3 C). The half-life of WT NF-E2, which is 1.35 h in control cells, is extended 1.65-fold to 2.23 h by the presence of truncated NF-E2 (Fig. 3 D). Expression of mutant NF-E2 therefore not only increases mRNA levels but, in addition, increases protein stability. Both mechanisms are likely to contribute to the increase in NF-E2 protein levels seen in NF-E2 mutant carrying MPN patients (Fig. 1 E).
MPN patients carrying NF-E2 mutations were also tested for the presence of the JAK2V617F mutation as well as mutations in cMpl, ASXL1, TET2, and EZH2. Three patients displayed mutations in TET2, whereas all seven patients assayed carried both a NF-E2 mutation and the JAK2V617F mutation (Table 1). Analysis of individual hematopoietic colonies grown in methylcellulose revealed that in three patients, 202, 209, and 409, the NF-E2 mutation was acquired subsequent to the JAK2V617F mutation, as some colonies carrying only the JAK2V617F mutation but not the NF-E2 mutation were found and all colonies positive for the NF-E2 mutation also carried JAK2V617F (Fig. 4, A-C). Interestingly, these assays revealed a proliferative advantage of cells carrying the NF-E2 mutation over cells carrying only JAK2V617F (Fig. 4 A-C). In all three patients, cells carrying both a NF-E2 mutation and the JAK2V617F mutation represent the vast majority of the colonies analyzed (81, 94, and 79% of the erythroid colonies, respectively), demonstrating that the JAK2V617F-positive cell that acquired a NF-E2 mutation outcompeted cells carrying the JAK2V617F mutation alone.
Patient 409 has acquired two separate and mutually exclusive NF-E2 mutations, c.236delC and c.234delT (Fig. 4 C and Table 1). The vast majority of both myeloid and erythroid colonies carry a homozygous JAK2V617F mutation. However, 67% of the erythroid and 55% of the myeloid colonies carry the c.236delC, whereas 12% of the erythroid and 11% of the myeloid colonies carry a c.234delT mutation (Fig. 4 C). The remaining colonies are WT for NF-E2, demonstrating that both mutations are acquired and again suggesting that both NF-E2 mutations confer a growth advantage, as they outcompete the JAK2V617F NF-E2 WT cells. The presence of both heterozygous and homozygous JAK2V617F clones containing either the c.234delT or the c.236delC NF-E2 mutation implies that in this patient, uniparental disomy (loss of heterozygosity) has occurred in two independent clones, a rare situation previously reported in individual MPN patients (Schaub et al., 2010; Wang et al., 2011).
In patient 532, it cannot be determined for sure whether the NF-E2 or the JAK2V617F mutation was incurred first. However, again, the presence of both mutations clearly provided a proliferative advantage over either of the mutations alone, as all colonies carry both mutations (Fig. 4 D). In addition, in the context of an NF-E2 mutation, as in WT NF-E2 patients, homozygous JAK2V617F cells outgrow heterozygous JAK2V617F cells.
As acquisition of the NF-E2 mutation consistently conferred a selective advantage, so that cells carrying both mutations outnumbered those carrying JAK2V617F alone (Fig. 4, A–D), we sought to determine the mechanism mediating this effect. JAK2V617F has been shown to confer growth factor independence on the Epo- and IL-3–dependent cell lines FDCP-1-EpoR and BaF3, respectively (James et al., 2005). We transduced FDCP-1-EpoR and BaF3 cells carrying the JAK2V617F mutation either with a vector expressing the 262aa NF-E2 truncation mutant or with an empty control vector. The levels of endogenous, WT NF-E2 and truncated, mutant NF-E2 were determined by Western blotting and shown to be comparable in their relative ratio to those observed in MPN patients; mutant NF-E2 was present at much lower levels than WT NF-E2 (compare Fig. 4 E with Fig. 1 E).
FDCP-1-EpoR and BaF3 cells were starved of growth factors and cell numbers recorded in the absence of Epo or IL-3, respectively (Fig. 4, F and G, left). As published previously (James et al., 2005), cells carrying the JAK2V617F mutation alone display growth in the absence of cytokines. Consistently, in the absence of cytokines, cells expressing mutant NF-E2 in addition to the mutant JAK2V617F kinase grew faster than cells expressing only JAK2V617F, reaching statistically significantly higher cell numbers (Fig. 4, F and G, left). In the presence of Epo or IL-3, all cells proliferated equally (Fig. 4 F and G, middle and right). The presence of mutant NF-E2 therefore enhanced cytokine independence and conferred a proliferative advantage in the absence of cytokines beyond that conferred by JAK2V617F alone. This was true for both the Epo independence of FDCP-1-EpoR cells and the independence from IL-3 in BaF3 cells (Fig. 4, F and G).
The proliferative advantage conferred by mutant NF-E2 was further investigated by conducting cell cycle analysis on FDCP-1-EpoR cells expressing both JAK2V617F and the NF-E2 262aa mutant as well as FDCP-1-EpoR cells expressing JAK2V617F alone (Fig. 4 H). Cells carrying mutant NF-E2 displayed a significant increase in the proportion of cells in S-phase and a concomitant significant decrease in the proportion of cells in G1-phase (Fig. 4 I). Expression of mutant NF-E2 therefore enhances cell division, causing increased proliferation.
To investigate the molecular basis for the proliferative advantage conferred by mutant NF-E2, we conducted microarray analysis comparing FDCP-1-EpoR cells expressing JAK2V617F alone to those carrying both JAK2V617F and the NF-E2 262aa mutant. We combined these results with ChIP-seq data available in silico, depicting in vivo NF-E2 DNA binding in K562 cells (Consortium, 2011). A series of candidates was selected whose RNA expression was increased in cells carrying both JAK2V617F and the NF-E2 262aa mutant compared with cells expressing JAK2V617F alone and whose promoter or enhancer regions displayed NF-E2 binding, suggesting that this transcription factor is involved in regulating their expression. Several candidate genes encode cell cycle regulators. Protein expression of three candidates, cyclin D3, CDK4 and CDK6, was examined by Western blotting (Fig. 4, J and K). Cells carrying both JAK2V617F and the NF-E2 262aa mutant expressed higher levels of all three of these cell cycle regulators. These data provide a molecular mechanism for the enhanced proliferation conferred by mutant NF-E2: by increasing expression of cell cycle regulators that promote the G1-S transition (Sherr, 1995; Herzinger and Reed, 1998), mutant NF-E2 causes an increase in the number of cells entering S-phase, hence promoting proliferation.
To test whether the NF-E2 mutants also confer an MPN phenotype in vivo, a murine bone marrow transplant model was used. Donor bone marrow was transduced with empty lentiviruses or viruses encoding WT NF-E2, the 248aa, the 262aa, or the Δ297-300 mutants and transplanted into recipient mice (Fig. 5 A). Donor engraftment exceeded 90% in all mice.
We determined NF-E2 mRNA and protein expression in this novel BMT model. Bone marrow was harvested by aspiration and NF-E2 protein levels assayed by Western blotting. Similar to what we observed in MPN patients (Fig. 1 E), truncated, mutant NF-E2 is expressed at low levels in our murine model, equivalent to those observed in patients (Fig. 5, B and C; in Fig. 5 C, mutant NF-E2 is 26% of WT in the transplanted mouse and 29% of WT in patient 202), demonstrating that the mouse model faithfully represents the pathophysiologic situation in MPN patients.
NF-E2 mRNA levels were determined by qRT-PCR. As noted above (Fig. 3 A), MPN patients carrying NF-E2 mutations display significantly elevated NF-E2 mRNA levels compared with patients with WT NF-E2. A similar effect was observed in mice carrying NF-E2 mutations. Two of the three mutants likewise express significantly elevated levels of murine NF-E2 mRNA compared with control mice expressing WT NF-E2 (Fig. 3 B).
To determine the effect of mutant NF-E2 on the transcription of a well characterized target gene, we measured β-globin expression in cells from control and NF-E2 mutant mice. Cells at equivalent stages of erythroid differentiation were isolated by FACS sorting of CD71, Ter-119 double-positive cells. β-globin mRNA expression was determined by qRT-PCR in this sorted population. As depicted in Fig. 2 F, β-globin expression is significantly elevated in NF-E2 mutant carrying mice compared with control mice. Moreover, the degree of β-globin overexpression in mice carrying the three different NF-E2 mutations precisely mirrors that previously observed in transiently transfected cells (Fig. 2, compare E and F). Therefore, the truncated mutant NF-E2 protein is expressed and functional in the BMT model.
After transplantation, peripheral blood was analyzed at monthly intervals. Compared with mice transplanted with empty vector or with vector expressing WT NF-E2, which, similar to NF-E2 transgenic mice in the first year of life, display normal platelet values (Kaufmann et al., 2012), mice transplanted with all three mutants displayed a statistically significant elevation in platelet number beginning 1–2 mo after transplantation (Fig. 5 D). Moreover, although WT NF-E2 transplanted mice, like NF-E2 transgenic mice, display normal hematocrit and normal RBCs (Kaufmann et al., 2012), mice transplanted with NF-E2 mutants displayed an elevation in hematocrit and in RBC number (Fig. 5, E and F). Furthermore, the absolute neutrophil count (ANC) is elevated in NF-E2 mutant mice beginning as early as 1 mo after transplantation (Fig. 5 G). These data demonstrate that presence of mutant NF-E2 is sufficient to cause thrombocytosis, erythrocytosis, and neutrophilia in a murine model.
We investigated whether expression of mutant NF-E2 alters the proportion of myeloid and lymphoid cells in the peripheral blood. FACS analysis of peripheral blood leukocytes demonstrates that the proportion of Gr-1+-Mac-1+ myeloid cells is statistically significantly increased in mice carrying NF-E2 mutations (Fig. 6 A and Fig. S1). Consequently, the proportion of lymphocytes is statistically significantly reduced (Fig. 6 B). Among the lymphocytes, the proportion of B cells remained constant, whereas the proportion of T cells was significantly reduced (Fig. 6, C and D; and Fig. S1). In MPN patients, the B cells are frequently part of the malignant clone (Raskind et al., 1985), whereas the T cells often are not (el-Kassar et al., 1997). Therefore, mutant NF-E2 causes a lineage bias and increases the production of those lineages affected in MPN patients.
To determine the origin of the observed lineage bias, we analyzed the stem and progenitor compartments in NF-E2 mutant-bearing mice. Compared with control mice, NF-E2 mutant-expressing mice show a significant expansion of KSL cells (Kit+, Sca+, lin− cells; Fig. 6 E and Fig. S2). Moreover, the size of the KL compartment (Kit+, Sca−, lin− cells) is more than doubled in NF-E2 mutant mice, compared with empty control transduced mice (Fig. 6 F). This expansion is the result of an increase in the numbers of all three types of KL progenitors, megakaryocyte/erythroid progenitor (MEP), granulocyte/monocyte progenitor (GMP), and common myeloid progenitor (CMP; Fig. 6, G–I; and Fig. S2). An increase of the latter three cell types in PV patients has previously been reported by Jamieson et al. (2006).
We have recently shown that transgenic mice overexpressing WT NF-E2 contain elevated numbers of KSL and KL cells (Kaufmann et al., 2012). In addition, WT NF-E2 transgenic mice display increased numbers of CMP and MEP but not GMP (Kaufmann et al., 2012). In contrast, NF-E2 mutant-bearing mice contain an expansion of all three progenitor types, CMP, MEP, and GMP, and show a highly significant increase in the number of KL cells, CMP, and GMP over WT NF-E2–expressing mice (Fig. 6, F–I). Therefore, mutant, truncated NF-E2 increases and augments the myeloid lineage bias already evident in NF-E2 WT transgenic mice. In the murine model, as in the cell lines, mutant NF-E2 enhances and extends the effect of WT NF-E2.
Histopathological analysis was conducted in mice sacrificed at 6–7 mo of age. As shown in Fig. 5, mice expressing WT NF-E2, similar to NF-E2 transgenic mice (Kaufmann et al., 2012), did not show changes in blood counts in the first 6 mo after transplant (Fig. 5, D–G). Accordingly, bone marrow of mice expressing WT NF-E2 is histologically normal and similar to control bone marrow (Fig. 7, F and G). Mice expressing mutant NF-E2, in contrast, display an increased ANC (Fig. 5 G) as well as a myeloid lineage bias (Fig. 6). Correspondingly, bone marrow histologies of mutant NF-E2–expressing mice show increased granulopoiesis and display a right shift, manifested by the increased presence of mature forms, which can be found extending to the trabecula (Fig. 7, H–J).
Similar to what we have shown in WT NF-E2 transgenic animals (Kaufmann et al., 2012), animals transplanted with WT NF-E2 show an increase in siderin deposits in the spleen (Fig. 7 B), as well as an expanded red pulp, compared with control mice (Fig. 7 A). Mice expressing mutant NF-E2 show changes in splenic architecture very similar to those observed in WT NF-E2 transduced mice (Fig. 7, C–E).
The livers and lungs of NF-E2 WT and mutant NF-E2–expressing mice display normal histology in all animals, notably an absence of myeloid infiltrates (Fig. 7, K–N and O–R). We wished to determine whether the cooperation between JAK2V617F and mutant NF-E2 observed in MPN patients and in vitro models (Fig. 4) can be recapitulated in vivo.
To correctly model the order of acquisition of the two mutations, the following experiment was conducted (Fig. 8 A). Because the patients first acquire JAK2V617F and subsequently mutant NF-E2 (Fig. 4, A–C), we infected murine bone marrow from balb/c mice with JAK2V617F and transplanted recipient mice. This JAK2V617F model is well documented in the literature (Bumm et al., 2006; Wernig et al., 2006; Zaleskas et al., 2006). After 10 wk, the transplanted bone marrow was harvested and infected either with a lentivirus expressing mutant NF-E2 (pLeGO-iC2-NF-E2-262aa) or an empty LeGO-iC2 control. Cells expressing both alleles, JAK2V617F and either mutant NF-E2 or empty control, were sorted by FACS. Subsequently, 50,000 double-positive cells were transplanted, together with 300,000 bone marrow support cells, into lethally irradiated recipient mice (Fig. 8 A).
As previously published, mice expressing JAK2V617F alone (with the empty pLeGo-iC2 control) display an increased hemoglobin of 17 g/dl (Bumm et al., 2006; Wernig et al., 2006; Zaleskas et al., 2006). Mice expressing both JAK2V617F and mutant NF-E2 display a further, statistically significant increase in hemoglobin (Fig. 8 B). Furthermore, addition of mutant NF-E2 significantly increased the MCV, also previously shown to be elevated in the JAK2V617F BMT model (Lacout et al., 2006), and increased MCH, as well as the WBC count, above that affected by JAK2V617F alone (Fig. 8, C–E). These data demonstrate that in vivo mutant NF-E2 cooperates with JAK2V617F to augment the MPN phenotype.
We evaluated bone marrow histologies of mice expressing JAK2V617F alone (with the empty pLeGo-iC2 control) and mice expressing both JAK2V617F and mutant NF-E2. Compared with JAK2V617F-expressing mice, mice expressing both mutations show increased granulopoiesis and display a right shift, similar to mice expressing mutant NF-E2 alone (Fig. 8, F–I). Histological analysis thus supports the data from peripheral blood counts in that the presence of mutant NF-E2 confers a phenotype beyond the changes affected by the expression of JAK2V617F alone.
Based on our previous data documenting a role for aberrant NF-E2 expression in the pathophysiology of MPNs (Kaufmann et al., 2012), we hypothesized that MPN patients may harbor alterations in the NF-E2 gene. Here, we report the presence of insertion and deletion mutations in patients with PV and PMF, which lead to the expression of truncated NF-E2 proteins (Fig. 1). Although these truncation mutants demonstrate neither DNA binding nor transcriptional transactivating activity on their own (Fig. 1 C; and Fig. 2, A–C), in the presence of WT NF-E2, they augment the effect of this transcription factor, both in vitro and in vivo (Fig. 2, E and F; Fig. 3, A and B; Fig. 5, D–G; and Fig. 6).
We demonstrate two distinct effects of truncated NF-E2 on the expression of its WT counterpart. First, mutant NF-E2 increases the steady-state levels of WT NF-E2 mRNA (Fig. 3, A and B). Second, the presence of mutant NF-E2 increases protein stability of WT NF-E2 (Fig. 3, C and D). How the truncated protein causes these effects mechanistically remains to be investigated. Increased amounts of steady-state mRNA can arise from an increased rate of transcription, from delayed mRNA degradation, or from a combination of both effects. Run-on transcription assays, as well as the determination of mRNA decay, will be used to investigate by which mechanisms truncated NF-E2 proteins increase WT NF-E2 mRNA levels.
The presence of mutant NF-E2 augments the levels and activity of WT NF-E2 protein (Fig. 2, E and F; Fig. 3 C; Fig. 5, D–G; and Fig. 6). Because the truncated proteins lack a DNA binding domain, they most likely exert their effect via protein–protein interactions. However, as the leucine zipper domain required for interactions with other β-zip proteins is also deleted or inactivated, these protein–protein interactions must take place within the remaining portion of the protein, the N-terminal transactivating domain, which extends from amino acids 1 to 206. The shortest truncating mutant detected in our cohort retains only the first 79 amino acids in frame (Fig. 1 A and Table 1). Interestingly, amino acids 1–79 constitute the proline-rich domain, which contains two heme binding sites (Moore et al., 2006).
NF-E2 has been reported to interact with several partners through sites within its transactivating domain. These include the ubiquitin ligases ITCH and NEDD4, which bind the NF-E2 PY domain (aa 79–83) through their WW domains (Gavva et al., 1997; Chen et al., 2001; Lee et al., 2008). This interaction domain is retained in all truncation mutants except for the shortest p.P79LfsX32 protein. In addition, NF-E2 binds the transcriptional coactivators YAP and CBP (Gavva et al., 1997; Hung et al., 2001) as well as a member of the TFIID transcription complex, TAFII130 (Amrolia et al., 1997).
Revealing a potential role as both an activator and a repressor, NF-E2 also binds the transcriptional repressor NAPP2, which acts epigenetically by recruiting HDAC1 and causing histone deacetylation (Gavva et al., 2002). NF-E2 is known to exert further epigenetic functions. It has been shown to remodel chromatin at the β-globin locus, establishing two chromatin marks, histone 3 trimethylated lysine 4 and trimethylated lysine 79 (Kiekhaefer et al., 2004). NF-E2 is required for recruitment of the MLL2 (mixed-lineage leukemia 2) and the G9a methyltransferase complexes to the β-globin gene and interacts physically with several of their members (Demers et al., 2007; Chaturvedi et al., 2009).
Presence of truncated NF-E2 proteins, which may still be able to bind their interaction partners but lack additional domains required for full function, could perturb any of the protein–protein interactions detailed above. We are currently investigating which interactions are retained by the truncated NF-E2 mutants and whether these alter the functions of the resulting complexes. Alternatively or in addition, mutant protein could serve as a decoy, occupying binding sites in regulatory proteins, thereby preventing the degradation of WT NF-E2 leading to its accumulation. These hypotheses are currently being addressed experimentally.
In summary, our data identify NF-E2 as a novel mutational target in MPN patients. Based on our studies, these alterations are rare in MPN patients, affecting ∼2% of cases. We report truncating mutations of NF-E2 that enhance WT NF-E2 function, promote erythrocytosis and thrombocytosis in a murine model, and cooperate with JAK2V617F to augment the MPN phenotype. Our data underscore the role of elevated NF-E2 activity in the pathophysiology of MPNs.
MATERIALS AND METHODS
MPN patients and healthy controls.
Peripheral blood samples were obtained from PV (n = 144), ET (n = 120), MF (n = 192, including 139 PMF), 53 post-PV or post-ET secondary MF patients), MPN-U (n = 49), MDS (n = 57), post-MPN-sAML (n = 39), and CMML (n = 67) patients, all fulfilling the WHO criteria for diagnosis (Swerdlow et al., 2008). Patient characteristics are summarized in Table S1. Buffy coats of healthy volunteer blood donors were obtained from the University Hospital Freiburg Center for Blood Transfusion and from Sanquin Bloodsupply, Nijmegen. Isolated granulocytes, obtained by dextran sedimentation and Ficoll gradient centrifugation as previously described (Temerinac et al., 2000), were used for DNA preparation and subsequent NF-E2 sequence analysis. Buccal swabs or isolated T cells (Nikoloski et al., 2010) were used to obtain DNA for germline analysis. The study protocol was approved by the local ethics committees (University Hospital Freiburg, University Hospital Ulm, Radboud University Nijmegen Medical Centre, as well as the Member Institutions of the MPD-RC) and informed consent was obtained from all patients. Each patient was assigned a unique patient number (UPN), which was used thereafter for the protection of privacy. All patients save one were tested for the presence of the JAK2V617F mutation by qRT-PCR or LNA PCR as previously described (Steimle et al., 2007).
Exons 2 and 3 of the NF-E2 locus were amplified by PCR and subjected to direct sequencing according to standard protocols. Primer sequences are available upon request. Sequence variations were confirmed by an independent PCR amplification of the original genomic DNA and subsequent sequence analysis. Constitutive DNA, obtained either from buccal swabs (patients 532 and 980) or isolated T cells (patient 209), was subjected to PCR amplification and sequencing to determine the acquired nature of the detected mutations. Detected mutations were described according to the consensus nomenclature (http://www.hgvs.org/mutnomen; den Dunnen and Antonarakis, 2001).
HEK-293 nuclear extracts were prepared as previously described (Schreiber et al., 1989). For EMSA analysis, 1–2 µg of nuclear extracts were added to a binding reaction containing 2 µg poly dI-dC, 2 µl 10× binding buffer (10 mM Hepes, pH 7.9, 5 mM MgCl2, 30 mM KCl, 1 mM EDTA, 1 mM dithiothreitol, 12% glycerol), and 0.5 ng 32P-labeled oligonucleotide. The oligonucleotide containing the NF-E2 consensus binding site at bp −160 of the human porphobilinogen deaminase gene promoter (Mignotte et al., 1989) was synthesized by Eurofins MWG Operon and the NF-κB consensus oligonucleotide was purchased from Promega. When indicated, a 100-fold excess of nonradioactive oligonucleotide, an NF-E2 antibody (Sigma-Aldrich), or an NF-κB p65 antibody (Santa Cruz Biotechnologies, Inc.) were added. The reaction was incubated at room temperature for 15 min. For supershift assays, extracts and antibodies were preincubated for 10 min at room temperature before addition of the radioactive nucleotide.
Murine bone marrow, peripheral blood MNCs from MPN patients, CB3 cells, or HEK-293 cells, transduced with expression vectors encoding WT NF-E2 or mutants, were lysed in SC buffer (50 mM Tris/HCl, pH 8.0, 170 mM NaCl, 0.1% vol/vol NP-40, 50 mM NaF, 2 mM Na3VO4, 0.2 mM DTT, 20% vol/vol glycerol, and 1 µg/ml BSA) supplemented with 1× Complete (Roche). Cell lysates were subjected to SDS-PAGE-gel-electrophoresis and Western blotting. Primary antibodies against NF-E2 (Goerttler et al., 2005) were used as described. FDCP-1-EpoR cells expressing JAK2V617F as well as the 262aa NF-E2 truncation mutant or carrying an empty control vector were analyzed for CDK4 (clone DCS 156; Cell Signaling Technology), CDK6 (clone DCS 83; Cell Signaling Technology), and Cyclin D3 (clone DCS 22; Cell Signaling Technology) expression.
The blots were stripped and reprobed against β-actin (Clone AC-15; Sigma-Aldrich) or GAPDH (Clone 71.1; Sigma-Aldrich) to control for equal loading. Immune complexes were detected using chemiluminescence Western Blotting Reagents (GE Healthcare).
Plasmids: NF-E2 reporter constructs and lentiviral vectors.
A 5.2 kb region, spanning bp −1 to bp −5221 of the β1 tubulin gene, was amplified by PCR from human genomic DNA and inserted into the pGL3basic luciferase reporter vector (Promega), resulting in the β1-tubulin −5.2 kb-luc construct. The pRBGP2-Luciferase (pRBGP2-Luc) reporter plasmid (Igarashi et al., 1994), which contains three copies of the NF-E2 binding site from the chicken β-globin enhancer in a TATA-luciferase reporter vector, was a gift of M. Yamamoto (Tsukuba University, Tsukuba, Japan).
Insertion and deletion mutations detected in MPN patients were inserted into the NF-E2 cDNA by site-directed mutagenesis using the GeneTailor site-directed mutagenesis system (Invitrogen). WT and mutant NF-E2 cDNAs were cloned into pLeGO-iG, pLeGO-iT2, pLeGo-iC2 (Weber et al., 2008), and pRc/CMV vectors by restriction enzyme digestion. The MafG-pCMV6-XL4 expression vector was purchased from OriGene.
Transient transfections and luciferase assays.
HEK-293 cells were transiently transfected with 0.2 µg of either the pRBGP2-Luc or the β1-tubulin-Luciferase reporter gene construct, together with 1.36 µg of either pRc/CMV-NF-E2 WT or a pRc/CMV-NF-E2 mutant expression vector or the pRc/CMV empty control vector as well as 0.34 µg of a MafG-pCMV6-XL4 expression vector (OriGene). In addition, 0.1 µg of a TK-Renilla plasmid (Promega) were cotransfected as an internal control.
Cells were harvested 16 h after transfection, and luciferase activity was determined using the Dual Luciferase Reporter Assay System (Promega). Luciferase activity was normalized to the Renilla internal control to compensate for variations in transfection efficiency.
CB3 cells, a kind gift of V. Blank (McGill University, Montreal, Quebec, Canada) and E. Bresnick (University of Wisconsin, Madison, WI), were transduced with empty pLeGO-iG, pLeGO-iT, or pLeGO-iC2 containing WT NF-E2 or NF-E2 mutants as previously described (Roelz et al., 2010) using an MOI between 1 and 8. GFP-expressing cells were sorted to obtain pure populations. CB3-NF-E2wt cells were obtained by transducing CB3 cells with pLeGO-iG-WT-NF-E2 and subsequently cloning single transduced cells by limiting dilution. A stable clone was selected and used for coexpression of NF-E2 mutants and variants. CB3-NF-E2wt cells were transduced with either pLeGO-iT2 or pLeGO-iC2 containing mutant and variant NF-E2 cDNAs or the empty control vector using an MOI between 0.5 and 9. Tomato or Cherry and GFP-expressing cells were sorted to obtain pure populations.
FDCP-1-Epo-R and FDCP-1-EpoR-JAK2V617F cells were transduced either with pLeGO-iC2 expressing the 262aa NF-E2 truncation mutant or the empty control vector. BaF3 cells were transduced with both the pMIG retrovirus expressing JAK2V617F as well as with pLeGO-iC2 expressing the 262aa NF-E2 truncation mutant or the empty control pLeGo-iC2 vector. All cells were sorted to obtain pure populations. Cells were starved of growth factors and subsequently counted daily to assess the cumulated cell number. As a control, FDCP-1-Epo-R cells were cultured in 1 IU/ml Epo and both BaF3 and FDCP-1-Epo-R cells in 5 ng/ml mIL-3.
Measurement of NF-E2 protein stability.
CB3 cells stably transduced with WT NF-E2 were subsequently transduced with either pLeGO-iC2 expressing the 262aa NF-E2 truncation mutant or the empty control vector and sorted to obtain pure populations. Cells were treated with 20 µg/ml cycloheximide to arrest protein synthesis and protein extracts collected at 1, 2, 4, and 6 h thereafter. NF-E2 and GAPDH levels were detected by Western blotting as described above and quantitated densitometrically. A nonlinear best-fit regression model was used to calculate the half-life of WT NF-E2 protein in the presence or absence of the 262aa truncation mutant.
Cell cycle analysis.
FDCP-1-EpoR cells expressing JAK2V617F, as well as the 262aa NF-E2 truncation mutant, or carrying an empty control vector were analyzed for cell cycle stage by staining with Hoechst 33342 (Molecular Probes) and 1 µg/ml LIVE/DEAD Fixable Near-IR Dead Cell Stain (Molecular Probes) at the manufacturer’s recommendation and analyzed by FACS using FlowJo software (Tree Star).
Quantitative RT-PCR assays.
Quantitative RT-PCR measurements were performed using an Assay on Demand (Mm01611268_g1 Hbb-b1; Applied Biosystems) for analysis of murine β-globin as well as the following primer and probe sequences for murine β-2-Microglobulin: forward primer, 5′-CTTTCTGGTGCTTGTCTCACTGAC-3′; reverse primer, 5′-GGTGGCGTGAGTATACTTGAATTTG-3′; TaqMan Probe, 6-FAM-ATCCAGAAAACCCC-MGBNFQ. NF-E2 mRNA was quantitated as previously described (Kaufmann et al., 2012). Reverse transcription of total CB3 cell RNA was performed using the Multiscribe Reverse transcription (Applied Biosystems). Q-PCR assays were performed in duplicate in an ABI PRISM 7000 Cycler and analyzed using the ABI PRISM 7000 software (Applied Biosystems). β-globin expression was determined relative to expression of the β-2-Microglobulin house-keeping gene using the ΔΔCt method.
Colony assays and mutation analysis.
Cells were seeded in methylcellulose media (1,000 cells/ml) containing SCF, IL-3, GM-CSF, and EPO (H4434; Stem Cell Technologies) and incubated for 14 d at 37°C with 5% CO2. Individual colonies (erythroid burst-forming units), as well as granulocytic/macrophage colony-forming units, were harvested on day 14 and assayed for the presence of the JAK2V617F mutation as well as the NF-E2 mutation. Presence of the JAK2V617F mutation was assessed by PCR amplification and BsaXI digestion, as well as by a JAK2V617F mutant-specific PCR, whereas presence of the NF-E2 mutation was detected by GeneScan analysis (patients 202, 209, and 532) and by direct sequencing (patients 209 and 409).
GeneScan analysis of NF-E2 mutations.
GeneScan analysis of genomic DNA extracted from single colonies was performed if the NF-E2 mutation comprised a deletion of three or more nucleotides (patients 202, 209, and 532). Two different primer pairs were used: forward primer, 5′-GACTTTAATGAGCTATTGGCAAGGTA-3′, and reverse primer, 5′-ACAATGGTTTCCAGCTTCCTCTT-3′ (coupled to 6FAM); or forward primer, 5′-CCCTACTCACTCATGCCCAACT-3′, and reverse primer, 5′-AAGCCGCTCCCGTTCATT-3′ (coupled to 6FAM). GeneScan analysis was performed on a 3730 DNA Analyzer using the GeneScan 500 LIZ Size Standard (Life Technologies). The generated peaks were analyzed using the Peak Scanner Software (Life Technologies).
Assessment of the JAK2V617F burden by PCR amplification and BsaXI digestion.
Detection of the JAK2V617F mutation by BsaXI digestion was performed essentially as previously described (Baxter et al., 2005). In brief, the following primers, located in exon 14 and intron 14 of the JAK2 gene, were used for PCR of genomic DNA: forward primer, 5′-AGCAGCAAGTATGATGAGCAAG-3′; reverse primer, 5′-GAGAAAGGCATTAGAAAGCCTGTAG-3′. The DNA fragment generated was subjected to BsaXI restriction digestion. The c.2343 G>T transversion, which generates the V617F mutation, destroys a BsaXI recognition sequence. Therefore, DNA containing only WT JAK2 alleles is cut by the enzyme, resulting in smaller fragments. Whereas DNA heterozygous for the JAK2V617F mutation generates smaller fragments, visible as one band, as well as the uncut PCR amplicon, the homozygous JAK2V617F DNA shows only the uncut PCR amplicon. Genomic DNA from the JAK2 WT EBV-DUZ cell line and from HEL cells, which harbor 10 copies of the JAK2V617F allele, was used as a negative and a positive control, respectively.
Assessment of the JAK2V617F burden by mutant-specific PCR.
Detection of the JAK2V617F mutation by a mutant-specific PCR assay was performed essentially as described previously (Sidon et al., 2006). Genomic DNA was amplified using the following primers: a forward primer recognizing both WT and mutant JAK2 alleles, 5′-TTATGGACAACAGTCAAACAACAATTC-3′; and a mutant-specific reverse primer, 5′-CTTACTCTCGTCTCCACAAAA-3′. The mutant-specific reverse primer contains one additional mutation to enhance specificity. Genomic DNA from the EBV-DUZ and HEL cell lines was used as a negative and a positive control, respectively.
Murine bone marrow transplantation.
FVB/N-45.2 donor mice (Guibal et al., 2009) were treated with 150 mg/kg body weight 5-fluorouracil (5-FU) 4 d before bone marrow harvest. Donor bone marrow was either infected with empty lentiviral particles (pLeGO-iG) or with lentiviruses expressing WT or mutant NF-E2 at an MOI of 10 as previously described (Roelz et al., 2010). FVB/N mice (45.1) were transplanted intrafemorally with lentivirally transduced FVB/N-45.2 bone marrow, as previously described (Roelz et al., 2010; Kaufmann et al., 2012). 12 wk after transplantation, engraftment exceeded 90% in all cases.
Bone marrow isolated from 5-FU–treated balb/c mice was infected with pMIG-JAK2V617F, as previously described (Bumm et al., 2006; Zaleskas et al., 2006) and transplanted into lethally irradiated recipients. 10 wk after transplantation, bone marrow was harvested and transduced with either the empty pLeGO-iC2 lentivirus or a lentivirus expressing the 262aa NF-E2-262 truncation mutant (pLEGO-iC2-NF-E2-262aa). Double-positive cells, expressing both GFP and Cherry, were FACS sorted. 50,000 double-positive cells were transplanted intrafemorally together with 300,000 support bone marrow cells into lethally irradiated recipients. Peripheral blood was obtained by retro-orbital puncture and analyzed on an Advia 120 system (Siemens).
FACS analysis of murine peripheral blood and bone marrow.
Murine peripheral blood was stained with antibodies against Gr1 (clone RB6-8C5; BioLegend), Mac-1 (clone M1/70; BioLegend), B220 (clone RA3-6B2,103212; BioLegend), and CD3ε (clone 145-2C11, 100308; BioLegend) and analyzed by FACS. For FACS sorting of populations to determine β-globin expression, antibodies against Ter-119 (clone Ter-119, 22155234; ImmunoTools) and CD71 (clone R17217; Leinco) were used.
Bone marrow was stained with a cocktail for lineage markers (CD3, Mac-1/CD11b, B220, Ter-119, and Gr1/Ly6G/6C; BioLegend) and lineage-negative cells (Lin−) were analyzed for c-Kit/CD117 (clone 2B8; eBioscience) and Sca-1 (clone D7; BioLegend) expression as previously described (Akashi et al., 2000; Christensen and Weissman, 2001; Passegué et al., 2003). Progenitor and stem cell subpopulations were identified as described by staining with the following additional markers: CD34 (clone MEC14.7; BioLegend), Fc-gamma-II/III-R (clone 93; eBioscience), as well as Flt3-biotin (clone A2F10; BioLegend) detected with streptavidin-PerCP (405213, BioLegend). Positively staining cells, as well as negative cells, were determined by fluorescence-minus-one (FMO) stainings, as previously published (Kaufmann et al., 2012).
Organs were fixed in 4% formalin, femurs subsequently decalcified in 10% buffered ethylene-diamine tetra-acetic acid (EDTA), pH 7.2, and all organs paraffin embedded. Sections were stained with H&E for histological analysis.
Paired or unpaired Student’s t tests (two sided or one sided) and the Mann-Whitney Rank sum test were used to determine whether a significant (P < 0.05) difference existed between two groups. When comparing more than two groups, a one- or two-way ANOVA was used. These analyses were performed using SigmaPlot 11.0 (Systat) or Prism 5.0 (GraphPad Software) software.
Online supplemental material.
Fig. S1 depicts the quantification of peripheral blood leukocytes (granulocytes, B and T cells), in PB of transplanted mice. Fig. S2 shows the quantification of hematopoietic stem and precursor (KL, KSL, CMP, GMP, and MEP) cells in BM of transplanted mice. Table S1 shows the clinical characteristics of the MPN, MDS, sAML, and CMML patients studied. Table S2 depicts the clinical features of MPN patients carrying NF-E2 mutations and Table S3 lists variants in NF-E2 that could not be verified as somatic.
The authors sincerely thank Martina de Groot for expert technical assistance, Sandra Kaiser and Titiksha Basu for lending a helping hand, and Prof. Dr. Torsten Haferlach, MLL, München for providing patient samples and gratefully acknowledge support from the MPD Research Consortium (MPD-RC) Tissue Bank, which also provided tissue samples. CB3 cells were kindly provided by Dr. Blank. The pRBGP2-Luciferase (pRBGP2-Luc) reporter plasmid was a kind gift of Dr. Masayuki Yamamoto (Tsukuba University, Tsukuba, Japan). FDCP-1-Epo-R and FDCP-1-EpoR-JAK2V617F cells were generous gifts of Dr. Stefan Constantinescu, Brussels. We gratefully acknowledge the gift of FVB/N-45.2 mice by Dr. Daniel G. Tenen, Harvard Medical School through Prof. Dr. Steffen Koschmieder, University Hospital Aachen. The pLeGO-iG, pLeGO-iT2, and pLeGO-iC2 vectors were generously provided by Prof. Dr. B. Fehse, Hamburg.
This work was supported by grants from the National Cancer Institute (PO1 CA108671 to H.L. Pahl), the Deutsche Forschungsgemeinschaft (Pa 611/6-1 to H.L. Pahl), and the Netherlands Institute for Regenerative Medicine to J.H. Jansen. H.L. Pahl is a member of the MPD Research Consortium. R. Bogeska was funded by the Excellence Initiative of the German Federal and States Governments (GSC-4, Spemann Graduate School).
The authors have no conflicting financial interests.
acute myeloid leukemia
absolute neutrophil count
common myeloid progenitor
electrophoretic mobility shift assay
myeloid erythroid progenitor
nuclear factor erythroid-2
unique patient number
white blood cell
R. Bogeska and G. Nikoloski contributed equally to this paper.