- 1Jiangsu Province Key Laboratory of Immunity and Metabolism, Xuzhou Medical University, Xuzhou, China
- 2Department of Pathogenic Biology and Immunology, Xuzhou Medical University, Xuzhou, China
- 3The Department of Hematology, Fukushima Medical University, Fukushima, Japan
- 4Department of Hematology, YUASA Foundation Jusendo General Hospital, Koriyama, Japan
- 5Department of Hematology, Iwaki City Medical Center, Iwaki, Japan
A substitution mutation of valine to phenylalanine at codon encoding position 617 of the Janus kinase 2 (JAK2) gene (JAK2V617F) has been detected in myeloid cells of some individuals with higher levels of proinflammatory cytokine production such as interleukin (IL)-6. However, the mechanisms by which JAK2V617F mutation mediating those cytokines remain unclear. We, therefore, established JAK2V617F-expressing murine macrophages (JAK2V617F macrophages) and found that the levels of p-STAT3 were markedly elevated in JAK2V617F macrophages in association with an increase in IL-6 production. However, inhibition of STAT3 by C188-9 significantly decreased the production of IL-6. Furthermore, the JAK2V617F mutation endowed macrophages with an elevated glycolytic phenotype in parallel with aberrant expression of PKM1. Interestingly, silencing of PKM1 inactivated STAT3 in parallel with reduced IL-6 production. In contrast, ectopic expression of PKM1 elevated IL-6 production via STAT3 activation. Importantly, the JAK2V617F mutation contributed to PKM1 protein stabilization via blockade of lysosomal-dependent degradation via chaperone-mediated autophagy (CMA), indicating that the JAK2V617F mutation could protect PKM1 from CMA-mediated degradation, leading to activation of STAT3 and promoting IL-6 production.
Introduction
JAK2 is a non-receptor tyrosine kinase that activates signal transduction and activator of transcription (STAT) and mediates cytokine signaling. For example, erythropoietin binds to its receptor which in turn phosphorylates JAK2/STAT3 and regulates erythropoiesis (1). The JAK2V617F mutation was identified in the majority of individuals with polycythemia vera (PV) and more than 50% of those with essential thrombocythemia (ET) and primary myelofibrosis (2). JAK2V617F mutation confers ligand-independent activation of STATs (3). Additionally, transplantation of JAK2V617F -expressing hematopoietic cells resulted in erythrocytosis in a murine model (4).
The clinical features of PV and ET include thromboembolism, which limits the survival of these patients. For example, the JAK2V617F mutation was detected in individuals who developed Budd-Chiari syndrome, a potentially lethal thrombotic disease, even though they did not have any thrombotic risk factors or complete blood cell count abnormalities (5). Also, the JAK2V617F mutation was found in some individuals with age-related hematopoiesis, which is intimately linked to a high incidence of cardiovascular events (5). Moreover, veno-occlusive disease was more frequently noted in individuals with myelodysplastic syndromes harboring the JAK2V617F mutation than in those without the mutation (6).
Previous studies showed the aberrant production of cytokines in patients with myeloproliferative neoplasms (MPN). For example, multiplex cytokine analyses examined the serum levels of a variety of cytokines and chemokines in more than 400 patients with MPN and found the aberrant production of epidermal growth factor and GRO-α in these patients (7). Importantly, these chemokines were produced by CD56+CD14+ proinflammatory monocytes and their serum levels were predictable of disease progression (7). Intriguingly, monocytes isolated from MPN patients were insensitive to IL-10, a negative regulator of pro-inflammatory cytokines, resulting in the sustained production of tumor necrosis factor-α (8). Another study found that monocytes with JAK2V617F expressed a greater amount of CD25, an activate marker of monocytes than those with wild-type JAK2. Also, monocytes isolated from MPN patients with thrombotic history expressed a greater amount of CD25 compared with those isolated from MPN patients without thrombotic history (9). Collectively these results suggest the possible roles of cytokine produced by monocytes in the development of thrombosis in MPN patients.
A large number of studies suggest that the relationship between activated cytokine signaling, including that of interleukin (IL)-1β and IL-6, which are predominantly produced by macrophages, is involved in the development of atherothrombosis (10–12). In response to stimuli, the main macrophage functions are under metabolic control (13). Upon treatment with lipopolysaccharide (LPS), a specific Toll-like receptor 4 (TLR4) ligand, macrophages undergo metabolic reprogramming, leading to increased glycolysis and a disrupted Krebs cycle to supply cell metabolic adaptations and cytokine production (13). Pyruvate kinase M2 (PKM2), an indispensable member of the pyruvate kinase family, has been shown to influence glycolytic reprogramming in activated immune cells and tumor cells and to function as a proinflammatory mediator (14, 15). Upon LPS stimulation, PKM2 interacts with hypoxia-inducible factor 1-α (HIF-1α) in macrophages (16, 17). Additionally, PKM2 can be shuttled into nuclei and phosphorylate STAT3 (18), leading to an increase in IL-6 production (16). These observations prompted us to hypothesize that JAK2V617F-expressing cells, especially macrophages, produce large amounts of cytokines, which could play a role in the development of thrombosis in individuals possessing this mutation. To test this hypothesis, we established JAK2V617F-expressing macrophages and explored their function.
Materials and Methods
Cells
Murine RAW264.7 and J774A.1 cells were commercially obtained from FuHeng (FuHeng Cell Center, Shanghai, China). These cells were maintained in Dulbecco’s modified Eagle’s medium (DMEM) containing 10% heat-inactivated FBS, 100 mg/L streptomycins, and 100 U/ml penicillin at 37°C in a 5% CO2 environment.
Gene Transfection and Silencing
For small interfering RNA (siRNA) transfections, cells were transfected with 20 nmol/L siRNA (Jima, Shanghai, China) by jetPRIME transfection reagent (Polyplus, France). For plasmid transfections, cells were transfected using jetPRIME (Polyplus, France) according to the manufacturer’s instructions. The siRNA sequences were as follows: siPKM1: GUGGAGGCCUCUUAUAAGUTT, siPPKM1#2 (the siRNA against PKM1 and PKM2): GCCACAGAAAGCUUUGCAUTT, siHSC70: CCAGGCCAGUAUUGAGAUUTT, and siLAMP2: GCCGUUCAGUCCAAUGCAUTT.
Plasmids and Virus Infection
Human JAK2 cDNA was synthesized and cloned into the pLenti-CMV-GFP-2A-Puro lentivirus vector (Abm, Canada) at NheI/XbaI sites. The JAK2V617F mutation was cloned using a TransStart FastPfu Fly DNA Polymerase Kit (TransGen Biotech). Lentiviruses were generated by cotransfecting 293T cells with the other two packing vectors pMD2G and psPAX and concentrated as described previously (19). Stable cell lines were generated by infection with lentivirus and selected with 2 mg/ml puromycin (Vicmed, Xuzhou, Jiangsu, China) for approximately 4 weeks. Murine pkm1 cDNA (NM_001253883.1) was synthesized and cloned into the pLenti-GIII-CMV-CBH-GFP-2A-Puro vector (Abcam, Canada).
Colony Formation Assay
For the colony formation assay, 200 infected cells were cultured in six-well plates at 37°C for 7–10 days, and visible colonies were washed twice with PBS, fixed with 4% paraformaldehyde, and stained with crystal violet. Pictures of the colonies were taken, and the number of colonies was counted by ImageJ software.
RNA Isolation and Real-Time Reverse Transcription-Polymerase Chain Reaction
We measured the expression of actin for normalization. Real-time RT-PCR was performed using TB Green PCR Master Mix (TaKaRa, Japan). The primer sets for the PCR are shown in Table 1. The PCR conditions for all the genes were as follows: initial activation at 95°C for 30 s, followed by 40 cycles at 95°C for 5 s and at 60°C for 20 s, and fluorescence determination was performed at the melting temperature of the product for 15 s on a LightCycler 480 (Roche).
Determination of mRNA Half-Life
To measure the half-life of endogenous HK1 or PKM1 messenger RNA (mRNA), cells were cultured in the presence of either control diluent or actinomycin D (2 µg/ml). Total RNA was extracted at the indicated time points and subjected to real-time RT-PCR. mRNA levels were normalized to actin and plotted as a percentage of the value at time 0 (set at 100%).
Immunoblotting
Immunoblotting was performed as previously described (20). For examination of the protein involved in pro-survival signaling pathway, cell lysis from empty-, JAK2-, or JAK2V617F-expressing macrophages were extracted and subpackaged equally for three groups (each group included empty, JAK2 and JAK2V617F lysates), followed by separation on 8% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) gel. After semi-dry transfer, the membranes were sequentially probed with the indicated antibodies. Anti-p-JAK2 (T1007/1008) (#3776), anti-JAK2 (D2E12) (#3230), anti-p-STAT1 (T701) (#7649), anti-STAT1 (D1K9Y) (#14994), anti-p-STAT3 (Tyr705, D3A7) (#9145), anti-STAT3 (D3Z2G) (#12640), anti-p-STAT5 (Tyr694, C11C5) (#9359), anti-STAT5 (D206Y) -(#94205), anti-p-p38 (Thr180/Tyr182, D3F9) (#4511), anti-p38 (D13E1) (#8690), anti-p-AKT (Ser473, D9E) (#10019), anti-AKT (#9272), anti-p-JNK (Thr183/Tyr185, 81E11) (#4668), and anti-JNK (#9252) antibodies were purchased from Cell Signaling Technology. Anti-HK1 (19662-1-AP), anti-HK2 (22029-1-AP), anti-PKM1 (15821-1-AP), anti-flag (Sigma), and anti-β-actin (66009-1-Ig) antibodies were purchased from ProteinTech.
Half-Life Assay
Cells were treated with 200 μg/ml cycloheximide and harvested using sodium dodecyl sulfate (SDS) lysis buffer at the indicated time points. The levels of HK2 and β-actin were analyzed by Western blotting. HK2 bands were quantified after normalization to those of β-actin, and the quantification was plotted as the relative amount of protein remaining compared to that at the 0 min treatment time. The bands were compared quantitatively using ImageJ software.
Immunoprecipitation
Immunoprecipitation was performed as described previously (21). Briefly, RAW264.7 lysates were prepared in immunoprecipitation lysis buffer (20 mmol/L Tris-Cl, pH 8.0, 10 mmol/L NaCl, 1 mmol/L EDTA, 0.1% NP-40) containing a protease inhibitor cocktail (Sigma). Two micrograms of cell extracts were precleared with 50 µl of protein A/G-agarose (Thermo Fisher) at 4°C for 2 h, and the supernatant was incubated with the corresponding antibodies with gentle shaking at 4°C overnight, followed by the addition of 20 µl of protein A/G-agarose for another 1 h. The beads were washed and then resuspended in 30 µl of loading buffer and boiled for 5 min, followed by Western blot detection.
Seahorse Analyzer
The extracellular acidification rate (ECAR) was measured with an XF 24 extracellular flux analyzer (Seahorse Bioscience). Briefly, 8x104 cells were seeded in each well of Seahorse XF 24 plates with 250 μl of DMEM and incubated overnight. ECAR measurements were normalized to the cell number. Cells were initially plated in XF Seahorse media with 2 mM glutamine in ECAR tests using the following concentrations of injected compounds, as indicated in the text: oligomycin, 2 µM; 2-DG, 100 mM; and glucose, 30 mM.
Glucose Uptake Assay
Glucose uptake was examined by a Glucose Assay Kit (Sigma) according to the manufacturer’s instructions.
Measurement of Serum Levels of IL-6 in MPN Patients
Serum was drawn from patients with PV or ET after obtaining written informed consent. Levels of IL-6 were measured using a human IL-6 ELISA kit (ProteinTech, Rosemont, IL, USA). This study was approved by the ethical committee of Fukushima Medical University.
Statistical Analysis
Data are expressed as means ± SEMs. Data shown are representative of one of two or three independent experiments. The comparison was performed by GraphPad Prism V5 (GraphPad Software). Mann-Whitney U test was used to analyze the comparison of IL-6 concentration in MPN patients with/without the JAK2V617F mutation. The student’s t-test was used to analyze the comparison of IL-6 concentrations in the knockdown or overexpression system. One-way analysis of variance (ANOVA) was used to analyze the comparison of empty-, JAK2-, or JAK2V617F-expressing macrophages. Upon stimuli, two-way ANOVA was used for comparison of empty-, JAK2-, or JAK2V617F-expressing macrophages.
Results
JAK2V617F Mutation Confers STAT Activation and Production of IL-6
JAK2V617F mutation is associated with the development of thrombosis, in which cytokines may be key orchestrators. To investigate whether and how JAK2V617F -expressing macrophages can produce a large amount of cytokines, such as IL-1β and IL-6, murine macrophage RAW264.7 cells were transduced with either an empty vector or the JAK2 or JAK2V617F gene by using lentiviruses, and sublines stably expressing these genes were established (Figure 1A). JAK2V617F macrophages aberrantly expressed the phosphorylated forms of JAK2 (p-JAK2) and STAT3 (p-STAT3) (Figure 1A). Additionally, the levels of p-JNK were also dramatically increased in JAK2V617F macrophages (Figure 1A). In contrast, the levels of p-STAT5 were modestly increased in JAK2V617F macrophages (Figure 1A). Interestingly, p-STAT1, the key mediator of interferon signaling involved in thrombus resolution (22), was undetectable in JAK2 and JAK2V617F macrophages (Figure 1A). Unexpectedly, the levels of p-AKT were decreased in JAK2V617F macrophages compared to those in empty vector-transduced macrophages, although the levels of p-AKT were elevated in JAK2 -expressing macrophages (Figure 1A). Similarly, JAK2V617F mutation endowed murine J774A.1 macrophages with higher levels of p-STAT3, but not p-STAT1, p-p38, p-AKT, or p-JNK, compared with that in JAK2-expressing J774A.1 macrophages (Figure S1A). The levels of p-STAT5 were slightly increased in JAK2V617F-expressing J774A.1 macrophages, compared with that in JAK2-expressing J774A.1 macrophages (Figure S1A).
Figure 1 IL-6 production by JAK2V617F-expressing cells in vitro and in vivo. (A) Western blot analysis. Proteins were extracted from the indicated cell lines and sub packaged equally for three groups (each group included the dame amount of empty, JAK2 and JAK2V617F lysates), followed by separation on 8% sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) gel. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (B) The relative expression of Il6 in cells treated with or without 100 ng/ml lipopolysaccharide (LPS) for the indicated times. The data shown are representative of one of three independent experiments performed in duplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. *P < 0.05. (C) The relative expression of Il1b in cells treated with or without 100 ng/ml LPS for the indicated times. The data shown are representative of one of three independent experiments performed in duplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01; *P < 0.05; n.s., not significant. (D) The IL-6 concentration of the cell culture medium was detected by ELISA. The data shown are representative of one of two independent experiments performed in duplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01; n.s., not significant. (E) Measurement of IL-6 in MPN patients with/without the JAK2V617F mutation. Blood was withdrawn from patients, and the levels of IL-6 were measured by ELISA. The data shown are representative of one of two independent experiments. The results represent the mean ± SEM. Statistical significance was assessed by the Mann-Whitney U test. **, P < 0.01. (F) Colony assay. The number of colonies was counted by ImageJ software after crystal violet staining. The data shown are representative of three independent experiments performed in duplicate. The results represent the mean ± SEM. One-way analysis of ANOVA was used for comparison. **P < 0.01; *P < 0.05; n.s., not significant.
Exposure of macrophages to LPS increased the expression levels of IL-6 in a time-dependent manner up to 6 h (Figure 1B). JAK2V617F macrophages expressed greater amounts of IL-6 mRNA than empty vector- or JAK2 expression vector-transduced RAW264.7 macrophages (empty or JAK2 macrophages) after exposure to LPS (Figure 1B). On the other hand, the levels of IL-1β mRNA were lower in JAK2V617F macrophages at 12 h after exposure to LPS than in empty macrophages (Figure 1C). The levels of IL-1β mRNA after exposure to LPS were comparable between JAK2V617F and JAK2 macrophages (Figure 1C). An ELISA showed that the production of IL-6 by JAK2V617F macrophages was greater than that by empty and JAK2 macrophages at baseline and after exposure to LPS for 6 h (Figure 1D).
We next measured the serum levels of IL-6 in MPN patients with/without the JAK2V617F mutation. Importantly, an increase in the levels of IL-6 was noted in PV and ET patients harboring JAK2V617F compared with those in patients without JAK2V617F (Figure 1E). White blood cells (WBCs) and values of C-reactive proteins (CRPs) were comparable between the two groups (Table 2). Of note, five thrombotic events were noted in MPN patients with JAK2V617F but not in patients without JAK2 mutations (Table 2).
Because the JAK/STAT pathway is involved in cell proliferation, we performed a colony assay and found that JAK2V617F macrophages formed more colonies than empty and JAK2 macrophages (Figures 1F and S1B). All of these observations suggested that the JAK2V617F mutation endowed macrophages with activation of STAT along with abundant production of IL-6.
Activation of STAT3 Was Required for IL-6 Production in JAK2V617F Macrophages
To determine whether STATs were critical for IL-6 production in JAK2V617F macrophages, we treated these cells with the STAT3 inhibitor C188-9. As shown in Figure 2A, upon LPS treatment, the levels of p-STAT3 were dramatically increased. By contrast, the levels of p-STAT3 were inhibited in the presence of C188-9 (Figure 2A). As expected, inhibition of STAT3 attenuated the IL-6 production induced by LPS (Figure 2B). Since STAT5 and JNK were activated in the presence of the JAK2V617F mutation, to test whether STAT5 and JNK were also involved in mediating IL-6 production in JAK2V617F macrophages, we blocked the STAT5 and JNK signal transduction pathways by specific inhibitors. We found that the inhibition of these pathways did not modulate the levels of IL-6 (data not shown). These observations suggested that activation of STAT3 was critical for IL-6 production.
Figure 2 Activation of STAT3 mediated by the JAK2V617F mutation was critical for IL-6 production. (A) Western blot analysis. Cells were cultured with 100 ng/ml LPS and/or 5 µM C188-9, and proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (B) The IL-6 concentration of the cell culture medium was detected by ELISA. The data shown are representative of one of two independent experiments performed in duplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01; *P < 0.0; n.s., not significant.
JAK2V617F Mutation Increases Glycolysis in Association with Constitutive Expression of PKM1 in Macrophages
In response to stimulation, macrophages undergo dynamic metabolic reprogramming (23). Activation of macrophages upon stimulation by LPS occurs in association with a metabolic switch from oxidative phosphorylation to glycolysis (24). To investigate whether an elevated glycolytic phenotype occurs in JAK2V617F macrophages, we compared the glycolytic capacities among empty, JAK2, and JAK2V617F macrophages. A Seahorse analyzer measured the ECAR in macrophages after sequential addition of glucose, oligomycin, and 2-DG (Figures 3A and S2A). Both glycolysis and glycolytic capacity were increased in JAK2V617F macrophages compared to those in empty and JAK2 macrophages, and these capacities were further enhanced in the presence of LPS (Figures 3B, C). Similarly, both glycolysis and glycolytic capacity were increased in JAK2V617F-expressing J774A.1 cells compared to those in empty-expressing J774A.1 cells, and these capacities were further enhanced in the presence of LPS (Figures S2B, C). Additionally, there’s a tendence that both glycolysis and glycolytic capacity were increased in JAK2V617F-expressing J774A.1 macrophages compared to those in JAK2-expressing J774A.1 cells, which were further enhanced in the presence of LPS (Figures S2B, C). On the other hand, glucose intake was almost identical in each subline (Figure 3D).
Figure 3 The JAK2V617F mutation endowed cells with an increase in glycolysis. (A) Cells treated with or without LPS were seeded into wells, and the ECAR was determined by extracellular flux analysis. A representative plot of the ECAR over time of these cells with the addition of glucose (30 mM), oligomycin (2 µM), and 2-DG (50 mM), as indicated. (B) Glycolysis in the assay shown in Figure 3A was quantified. The data shown are representative of one of two independent experiments performed in triplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. *P < 0.05. (C) The glycolytic capacity in the assay shown in Figure 3A was quantified. The data shown are representative of one of two independent experiments performed in triplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. *P < 0.05. (D) The levels of glucose uptake in the indicated macrophages. The data shown are representative of one of two independent experiments. The results represent the mean ± SEM. One-way analysis of ANOVA was used for comparison. n.s., not significant. (E) Western blot analysis. Proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (F) Western blot analysis. RAW264.7 cells were transfected with empty or Myc-PKM1 expression vectors. After 24 h, proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (G) The IL-6 concentration of the cell culture medium was detected by ELISA. The data shown are representative of one of two independent experiments performed in triplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01. (H) Western blot analysis. JAK2-expressing RAW264.7 cells were transfected with a negative control siRNA or an siRNA against PKM1. After 24 h, proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (I) The IL-6 concentration of the cell culture medium was detected by ELISA. The data shown are representative of one of two independent experiments performed in triplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01.
To explore the underlying mechanisms causing the elevated glycolytic phenotype noted in JAK2V617F macrophages, we next measured the expression levels of enzymes involved in glycolysis and found that pyruvate kinase M 1 (PKM1) but not hexokinase 2 (HK2) or PKM2 was aberrantly expressed in JAK2V617F macrophages compared to those in the other two sublines (Figure 3E). HK1 was slightly induced in JAK2V617F macrophages compared to in the other two sublines (Figure 3E). As expected, the expression of PKM1, but not HK1/2 or PKM2 was dramatically increased in JAK2V617F-expressing J774A.1 macrophages compared to those in the other two sublines (Figure S2D).
JAK2V617F Mutation-Induced Phosphorylation of STAT3 by PKM1
A previous study showed that PKM2 promotes proinflammatory cytokine production by mediating the activation of STAT3 (25). To further investigate the role of elevated PKM1 observed in JAK2V617F macrophages, we transfected RAW264.7 macrophages with PKM1 plasmid DNA and found that overexpression of PKM1 caused an increase in IL-6 production (Figures 3F, G). Additionally, overexpression of PKM1 increased the levels of p-STAT3 but not those of p-STAT5 and p-JNK (Figure 3F). Silencing of PKM1 decreased the levels of p-STAT3 (Figure 3H) in parallel with a decrease in IL-6 production (Figure 3I). Similarly, ectopic expression of PKM1 in J774A.1 macrophages contributed to higher levels of p-STAT3 in parallel with an increase in IL-6 production (Figures S2E, F). Whereas, knocking down of PKM1 reduced the expression of p-STAT3 along with a decrease in IL-6 production (Figures S2G, H), suggesting that the JAK2V617F mutation could endow macrophages with elevated IL-6 production by PKM1-mediated activation of STAT3.
JAK2V617F Mutation Contributed to PKM1 Protein Stabilization
The aberrant expression of PKM1 in JAK2V617F macrophages prompted us to ask whether the JAK2V617F mutation induced an increase in pkm1 gene expression. We found that in JAK2V617F macrophages, the level of PKM1 mRNA was increased by 2.5-fold compared to that in empty macrophages (Figure 4A). However, the level of PKM1 mRNA in JAK2 macrophages was almost identical to that in empty macrophages (Figure 4A). Similarly, the level of HK1 mRNA was increased in JAK2V617F macrophages compared to that in empty or JAK2 macrophages (Figure 4B). Next, actinomycin D was used to block de novo RNA synthesis (26), and the persistence of either PKM1 or HK1 mRNA was measured by real-time RT-PCR. The JAK2V617F mutation did not change PKM1 (Figure 4C) and HK1 mRNA degradation (Figure 4D), which excludes the possibility that the JAK2V617F mutation induced PKM1 expression via mainly stabilization of PKM1 mRNA. In addition, these results prompted us to investigate whether protein degradation was critical for JAK2V617F mutation-mediated PKM1 expression.
Figure 4 The JAK2V617F mutation facilitated PKM1 protein stabilization. (A) The relative expression of HK1 messenger RNA (mRNA) in empty vector-transfected RAW264.7, JAK2-expressing RAW264.7, or JAK2V617F-expressing RAW264.7 cells. The data shown are representative of one of three independent experiments performed in duplicate. The results represent the mean ± SEM. One-way analysis of ANOVA was used for comparison. *P < 0.05; n.s., not significant. (B) The relative expression of PKM1 mRNA in empty vector-transfected RAW264.7, JAK2-expressing RAW264.7, or JAK2V617F-expressing RAW264.7 cells. The data shown are representative of one of three independent experiments performed in duplicate. The results represent the mean ± SEM. One-way analysis of ANOVA was used for comparison. *P < 0.05; n.s., not significant. (C) Actinomycin D was added for the indicated period to block RNA synthesis, and HK1 mRNA was analyzed using real-time RT-PCR. The data shown are representative of one of two independent experiments performed in duplicate. Two-way analysis of ANOVA was used for comparison. (D) Actinomycin D was added for the indicated period to block RNA synthesis, and PKM1 mRNA was analyzed using real-time RT-PCR. The data shown are representative of one of two independent experiments performed in duplicate. Two-way analysis of ANOVA was used for comparison. (E–F) The levels of HK1, PKM1, and actin were analyzed by Western blotting. The amounts of HK1 and PKM1 in empty vector-transfected, JAK2-expressing, or JAK2V617F-expressing RAW264.7 cells were quantified by densitometry, normalized to the level of actin, and plotted. The data shown are representative of one of two independent experiments.
To elucidate whether the JAK2V617F mutation influences PKM1 protein degradation, we performed a protein half-life assay. In macrophages infected with an empty lentivirus or the JAK2 lentivirus, PKM1 had a modest half-life of approximately 6 h (Figures 4E, F). By contrast, in the presence of the JAK2V617F mutation, PKM1 degradation was hampered (Figures 4E, F). After 12 h of tracing, the protein levels of PKM1 were only slightly reduced in JAK2V617F macrophages (Figures 4E, F). However, HK1 was relatively stable, and its half-life appeared to be unaffected (Figure 4E). Together, these data strongly indicated that the JAK2V617F mutation played an important role in PKM1 protein degradation.
The JAK2V617F Mutation Facilitated HSC70 Binding
In search of the mechanisms by which the JAK2V617F mutation inhibited PKM1 degradation, we treated the three macrophage sublines with either the proteasome inhibitor MG-132 or the lysosomal protease inhibitor leupeptin and found that upon MG-132 treatment, the levels of PKM1 were not modulated (Figure 5A). On the other hand, the levels of PKM1 were increased in the presence of leupeptin in JAK2 and JAK2V617F macrophages (Figure 5B). A recent study indicated that acetylation of PKM2 on lysine promoted its lysosomal-dependent degradation via chaperone-mediated autophagy (CMA) (27), a pathway for selective degradation of individual proteins that is mediated through binding to heat shock cognate 71 kDa protein (HSC70) followed by unfolding and translocation of proteins through the lysosomal membrane by lysosome-associated membrane protein type 2A (LAMP2A) (28). Given that JAK2V617F macrophages possess higher levels of PKM1 (Figures 2–4) than empty macrophages, we hypothesized that lysosomal degradation of PKM1 occurred through CMA. To test this hypothesis, we transfected JAK2 macrophages with siRNA against either HSC70 or LAMP2A. As expected, silencing of HSC70 or LAMP2A induced an increase in the levels of PKM1 in JAK2V617F macrophages (Figures 5C, D).
Figure 5 The JAK2V617F mutation protected PKM1 protein from CMA-mediated degradation. (A) Western blot analysis. Cells were cultured with 10 µM MG-152 for 5 h, and proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (B) Western blot analysis. Cells were cultured with 10 µM leupeptin for 18 h, and proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of three independent experiments. (C) Western blot analysis. JAK2-expressing RAW264.7 cells were transfected with a negative control siRNA or an siRNA against HSC70. After 24 h, proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of three independent experiments. (D) Western blot analysis. JAK2-expressing RAW264.7 cells were transfected with a negative control siRNA or an siRNA against LAMP2A. After 24 h, proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of three independent experiments. (E) IP. Cell lysates were harvested from empty vector-transfected, JAK2-expressing and JAK2V617F-expressing RAW264.7 cells, submitted to IP for either HSC70 or PKM1, and immunoblotted with the indicated antibodies. The loading represents 5% of the total cell lysate utilized for IP. The results shown are representative of one of three independent experiments. (F) Signaling pathways in normal and JAK2V617F cells. In normal cells, HSC70 interacts with PKM1 and promotes PKM1 degradation (left panel). The JAK2V617F mutation sequesters PKM1 from binding with HSC70, resulting in PKM1 accumulation, which causes upregulation of the transcriptional activity of STAT3 and subsequently leads to enhancement of IL-6 production (right panel).
Since JAK2 interacts with HSC70 (29), we hypothesize that the JAK2V617F mutant protein could be a preferred form for HSC70 binding compared with wild-type JAK2, protecting PKM1 from CMA-mediated degradation. To test this hypothesis, we performed an immunoprecipitation (IP) analysis. IP of HSC70 pulled down a significant amount of flag in JAK2V617F macrophage lysates (Figure 5E, IP lane 3). In comparison, the same amount of HSC70 pulled down less flag in JAK2 macrophage lysates (Figure 5E, IP lane 2). Although it has been shown that HSC70 binds PKM2 (27), unfortunately, the IP of HSC70 did not pull down PKM1 in vivo (data not shown), perhaps because the anti-HSC70 antibody interfered with binding. IP of PKM1 pulled down a significant amount of HSC70 in JAK2 macrophage lysates (Figure 5E, IP PKM1, lane 2). In comparison, a higher amount of PKM1 pulled down less HSC70 in JAK2V617F macrophage lysates (Figure 5E, IP PKM1, lane 3).
Together, these data demonstrated that the JAK2V617F mutant protein bound HSC70 more efficiently than the wild-type JAK2 protein. They also suggested that the JAK2V617F mutant protein could protect PKM1 from CMA-mediated degradation through abrogating HSC70-PKM1 binding.
Discussion
Recent studies have demonstrated that the JAK2V617F mutation triggers the activation of STAT1 in myeloid cells. As a consequence, exposure of JAK2V617F-expressing myeloid cells to LPS produced greater amounts of IL-6 than that produced by control cells, which was consistent with the current study results (Figure 1). Interestingly, transplantation of JAK2V617F-harboring hematopoietic stem cells into mice accelerated cardiac hypertrophy and fibrosis, which were accompanied by a large number of infiltrating macrophages and elevated IL-6 expression in the heart (30). Another study also showed that the transplantation of JAK2V617F-expressing hematopoietic stem cells enhanced the formation of atherosclerotic lesions in lipoprotein receptor knockout mice. The authors also found that JAK2V617F macrophages produced large amounts of cytokines and chemokines (31). Collectively, these observations suggested that JAK2V617F macrophages acquire activated inflammatory properties and contribute to the development of cardiovascular diseases.
The present study found that forced expression of JAK2V617F resulted in an increase in the production of IL-6 but not IL-1β in macrophages even in the absence of LPS (Figure 1D and data not shown). Additionally, elevated serum levels of IL-6 were noted in MPN patients harboring JAK2V617F (Figure 1E, Table 2). To explore the mechanisms by which JAK2V617F mutant proteins cause aberrant expression of inflammatory cytokines in macrophages, we focused on metabolic reprogramming in these cells. Surprisingly, we found that JAK2V617F macrophages possessed increased levels of PKM1 in parallel with an increase in the glycolytic phenotype (Figures 3 and S2). Previous studies showed that PKM2 phosphorylated STAT3 (18), contributing to elevated IL-6 production (16). Similar to the role of PKM2 in immune cells (25), PKM1 also activates STAT3. However, in contrast to PKM2, forced expression of PKM1 endowed macrophages with an increase in the levels of IL-6 but not IL-1β (Figure 3), suggesting that PKM1 specifically regulates the expression of IL-6. To support this idea, silencing of PKM1 hampered IL-6 production but not IL-1β production (data not shown) in association with a decrease in the levels of p-STAT3 (Figure 3). Furthermore, a specific inhibitor of STAT3 downregulated endogenous expression of IL-6 in JAK2V617F macrophages (Figure 2). These observations suggested that an enhanced glycolytic phenotype resulted in activation of STAT3, leading to elevated production of IL-6 mediated by JAK2V617F; aberrant expression of PKM1 mediated by JAK2V617F could play an important role in regulating this process in macrophages.
Unlike PKM2, PKM1 forms a stable, constitutively active tetramer with high PK activity and is aberrantly expressed in differentiated cells, such as those of the muscle and brain (32). Whether tetrameric PKM1 could be converted into a PKM1 dimer following posttranslational modifications remains obscure. Interestingly, in JAK2V617F macrophages, PKM1 protein half-life was prolonged compared to that in their parental cell lines (Figures 4E, F), although the levels of PKM1 mRNA were also increased in JAK2V617F macrophages (Figure 4), indicating that the JAK2V617F protein could be involved in controlling PKM1 protein stabilization, at least partially. In a further study, we found that JAK2V617F interacted with HSC70, releasing PKM1 from HSC70-mediated degradation, which in turn could enhance PKM1 protein stabilization (Figure 5).
Importantly, JAK2/STAT3 is not the only signaling that governs the cytokine production in macrophages. The use of ruxolitinib, an inhibitor of JAK2 was not able to decrease cytokine production in MPN patients. Ex vivo studies found that inhibition of nuclear factor κB and MAPK signaling is also required to shut down the production of cytokines in MPN patients (33).
Altogether, our results demonstrate that the JAK2V617F mutation stimulates IL-6 production through activation of STAT3 signaling, which is at least in part mediated by PKM1-regulated glycolysis. Targeting the PKM1/STAT3/IL-6 axis may be a promising treatment strategy to prevent thrombosis and cardiovascular events in individuals harboring the JAK2V617F mutation.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
Conception and design: JY and TI. Development of methodology: RL, NS, and JY. Acquisition of data (provided animals, enrolled and managed patients, provided facilities, etc.): TI. Analysis and interpretation of data (e.g., statistical analysis, biostatistics, and computational analysis): JY and TI. Writing, review, and/or revision of the manuscript: JY and TI. Administrative, technical, or material support (i.e., reporting or organizing data and constructing databases): RL, NS, XC, XL, JZ, HH, WC, MFukat, HM, HT, HO, MFukam, MO, YH, KZ. Study supervision: TI. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the Jiangsu Distinguished Professorship Program to JY and KAKENHI (18H02844) to TI.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2020.589048/full#supplementary-material
Supplementary Figure 1 | Related to Figure 1. (A) Western blot analysis. Proteins were extracted from the indicated cell lines and sub packaged equally for three groups (each group included the same amount of empty, JAK2 and JAK2V617F lysates), followed by separation on 8% SDS-PAGE gel. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (B) Colony assay. The number of colonies was counted by ImageJ software after crystal violet staining. The data shown are representative of three independent experiments performed in duplicate. The results represent the mean ± SEM. One-way analysis of ANOVA was used for comparison. **P < 0.01; *P < 0.05; n.s., not significant.
Supplementary Figure 2 | Related to Figures 3. (A) J774A.1 macrophages treated with or without LPS were seeded into wells, and the ECAR was determined by extracellular flux analysis. A representative plot of the ECAR over time of these cells with the addition of glucose (30 mM), oligomycin (2 mM), and 2-DG (50 mM), as indicated. (B) Glycolysis in the assay shown in Supplementary Figure 2A was quantified. The data shown are representative of one of two independent experiments. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. *P < 0.05. (C) The glycolytic capacity in the assay shown in Figures 3A was quantified. The data shown are representative of one of two independent experiments. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. *P < 0.05. (D) Proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (E) J774A.1 macrophages were transfected with empty or Myc-PKM1 expression vectors. After 24 hrs, proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (F) The IL-6 concentration of the cell culture medium was detected by ELISA. The data shown are representative of one of two independent experiments performed in triplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01. (G) JAK2-expressing J774A.1 cells were transfected with a negative control siRNA or a siRNA against PKM1. After 24 hrs, proteins were extracted and subjected to Western blot analysis. The membrane was sequentially probed with the indicated antibodies. The results shown are representative of one of two independent experiments. (H) The IL-6 concentration of the cell culture medium was detected by ELISA. The data shown are representative of one of two independent experiments performed in triplicate. The results represent the mean ± SEM. Two-way analysis of ANOVA was used for comparison. **P < 0.01.
References
1. Hodges VM, Rainey S, Lappin TR, Maxwell AP. Pathophysiology of anemia and erythrocytosis. Crit Rev Oncol Hematol (2007) 64(2):139–58. doi: 10.1016/j.critrevonc.2007.06.006
2. Godfrey AL, Chen E, Pagano F, Silber Y, Campbell PJ, Green AR, et al. Clonal analyses reveal associations of JAK2V617F homozygosity with hematologic features, age and gender in polycythemia vera and essential thrombocythemia. Haematologica (2013) 98(5):718–21. doi: 10.3324/haematol.2012.079129
3. Grisouard J, Shimizu T, Duek A, Kubovcakova L, Hao-Shen H, Dirnhofer S, et al. Deletion of Stat3 in hematopoietic cells enhances thrombocytosis and shortens survival in a JAK2-V617F mouse model of MPN. Blood (2015) 125(13):2131–40. doi: 10.1182/blood-2014-08-594572
4. Wernig G, Mercher T, Okabe R, Levine RL, Lee BH, Gilliland DG. Expression of Jak2V617F causes a polycythemia vera-like disease with associated myelofibrosis in a murine bone marrow transplant model. Blood (2006) 107(11):4274–81. doi: 10.1182/blood-2005-12-4824
5. Foucar CE, Stein BL. JAK2 V617F Mutation Testing in Patients Presenting With Hepatic and Portal Vein Thrombosis. JAMA (2017) 317(21):2228–9. doi: 10.1001/jama.2017.2329
6. Naqvi K, Sasaki K, Montalban-Bravo G, Alfonso Pierola A, Yilmaz M, Short N, et al. Clonal hematopoiesis of indeterminate potential-associated mutations and risk of comorbidities in patients with myelodysplastic syndrome. Cancer (2019) 125(13):2233–41. doi: 10.1002/cncr.32056
7. Obro NF, Grinfeld J, Belmonte M, Irvine M, Shepherd MS, Rao TN, et al. Longitudinal Cytokine Profiling Identifies GRO-alpha and EGF as Potential Biomarkers of Disease Progression in Essential Thrombocythemia. Hemasphere (2020) 4(3):e371. doi: 10.1097/HS9.0000000000000371
8. Lai HY, Brooks SA, Craver BM, Morse SJ, Nguyen TK, Haghighi N, et al. Defective negative regulation of Toll-like receptor signaling leads to excessive TNF-alpha in myeloproliferative neoplasm. Blood Adv (2019) 3(2):122–31. doi: 10.1182/bloodadvances.2018026450
9. Goette NP, Lev PR, Heller PG, Kornblihtt LI, Korin L, Molinas FC, et al. Monocyte IL-2Ralpha expression is associated with thrombosis and the JAK2V617F mutation in myeloproliferative neoplasms. Cytokine (2010) 51(1):67–72. doi: 10.1016/j.cyto.2010.04.011
10. Moss JW, Ramji DP. Cytokines: roles in atherosclerosis disease progression and potential therapeutic targets. Future Med Chem (2016) 8(11):1317–30. doi: 10.4155/fmc-2016-0072
11. Dinarello CA. Interleukin-1 in the pathogenesis and treatment of inflammatory diseases. Blood (2011) 117(14):3720–32. doi: 10.1182/blood-2010-07-273417
12. Schuett H, Oestreich R, Waetzig GH, Annema W, Luchtefeld M, Hillmer A, et al. Transsignaling of interleukin-6 crucially contributes to atherosclerosis in mice. Arterioscler Thromb Vasc Biol (2012) 32(2):281–90. doi: 10.1161/ATVBAHA.111.229435
13. Stunault MI, Bories G, Guinamard RR, Ivanov S. Metabolism Plays a Key Role during Macrophage Activation. Mediators Inflammation (2018) 2018:2426138. doi: 10.1155/2018/2426138
14. Christofk HR, Vander Heiden MG, Harris MH, Ramanathan A, Gerszten RE, Wei R, et al. The M2 splice isoform of pyruvate kinase is important for cancer metabolism and tumour growth. Nature (2008) 452(7184):230–3. doi: 10.1038/nature06734
15. Alves-Filho JC, Palsson-McDermott EM. Pyruvate Kinase M2: A Potential Target for Regulating Inflammation. Front Immunol (2016) 7:145. doi: 10.3389/fimmu.2016.00145
16. Shirai T, Nazarewicz RR, Wallis BB, Yanes RE, Watanabe R, Hilhorst M, et al. The glycolytic enzyme PKM2 bridges metabolic and inflammatory dysfunction in coronary artery disease. J Exp Med (2016) 213(3):337–54. doi: 10.1084/jem.20150900
17. Palsson-McDermott EM, Curtis AM, Goel G, Lauterbach MA, Sheedy FJ, Gleeson LE, et al. Pyruvate kinase M2 regulates Hif-1alpha activity and IL-1beta induction and is a critical determinant of the warburg effect in LPS-activated macrophages. Cell Metab (2015) 21(1):65–80. doi: 10.1016/j.cmet.2014.12.005
18. Demaria M, Poli V. PKM2, STAT3 and HIF-1alpha: The Warburg’s vicious circle. JAKSTAT (2012) 1(3):194–6. doi: 10.4161/jkst.20662
19. Hou P, Li L, Chen F, Chen Y, Liu H, Li J, et al. PTBP3-Mediated Regulation of ZEB1 mRNA Stability Promotes Epithelial-Mesenchymal Transition in Breast Cancer. Cancer Res (2018) 78(2):387–98. doi: 10.1158/0008-5472.CAN-17-0883
20. Zhang B, Gu X, Han X, Gao Q, Liu J, Guo T, et al. Crosstalk between DNA methylation and histone acetylation triggers GDNF high transcription in glioblastoma cells. Clin Epigenet (2020) 12(1):47. doi: 10.1186/s13148-020-00835-3
21. Ayanlaja AA, Ji G, Wang J, Gao Y, Cheng B, Kanwore K, et al. Doublecortin undergo nucleocytoplasmic transport via the RanGTPase signaling to promote glioma progression. Cell Commun Signal (2020) 18(1):24. doi: 10.1186/s12964-019-0485-5
22. Nosaka M, Ishida Y, Kimura A, Kuninaka Y, Inui M, Mukaida N, et al. Absence of IFN-gamma accelerates thrombus resolution through enhanced MMP-9 and VEGF expression in mice. J Clin Invest (2011) 121(7):2911–20. doi: 10.1172/JCI40782
23. Seim GL, Britt EC, John SV, Yeo FJ, Johnson AR, Eisenstein RS, et al. Two-stage metabolic remodelling in macrophages in response to lipopolysaccharide and interferon-gamma stimulation. Nat Metab (2019) 1(7):731–42. doi: 10.1038/s42255-019-0083-2
24. Kelly B, O’Neill LA. Metabolic reprogramming in macrophages and dendritic cells in innate immunity. Cell Res (2015) 25(7):771–84. doi: 10.1038/cr.2015.68
25. van de Wetering C, Aboushousha R, Manuel AM, Chia SB, Erickson C, MacPherson MB, et al. Pyruvate Kinase M2 Promotes Expression of Proinflammatory Mediators in House Dust Mite-Induced Allergic Airways Disease. J Immunol (2020) 204(4):763–74. doi: 10.4049/jimmunol.1901086
26. An MX, Li S, Yao HB, Li C, Wang JM, Sun J, et al. BAG3 directly stabilizes Hexokinase 2 mRNA and promotes aerobic glycolysis in pancreatic cancer cells. J Cell Biol (2017) 216(12):4091–105. doi: 10.1083/jcb.201701064
27. Lv L, Li D, Zhao D, Lin R, Chu Y, Zhang H, et al. Acetylation targets the M2 isoform of pyruvate kinase for degradation through chaperone-mediated autophagy and promotes tumor growth. Mol Cell (2011) 42(6):719–30. doi: 10.1016/j.molcel.2011.04.025
28. Hubbi ME, Hu H, Kshitiz, Ahmed I, Levchenko A, Semenza GL. Chaperone-mediated autophagy targets hypoxia-inducible factor-1alpha (HIF-1alpha) for lysosomal degradation. J Biol Chem (2013) 288(15):10703–14. doi: 10.1074/jbc.M112.414771
29. Sarkar S, Pollack BP, Lin KT, Kotenko SV, Cook JR, Lewis A, et al. hTid-1, a human DnaJ protein, modulates the interferon signaling pathway. J Biol Chem (2001) 276(52):49034–42. doi: 10.1074/jbc.M103683200
30. Sano S, Wang Y, Yura Y, Sano M, Oshima K, Yang Y, et al. JAK2 (V617F) -Mediated Clonal Hematopoiesis Accelerates Pathological Remodeling in Murine Heart Failure. JACC Basic Transl Sci (2019) 4(6):684–97. doi: 10.1016/j.jacbts.2019.05.013
31. Wang W, Liu W, Fidler T, Wang Y, Tang Y, Woods B, et al. Macrophage Inflammation, Erythrophagocytosis, and Accelerated Atherosclerosis in Jak2 (V617F) Mice. Circ Res (2018) 123(11):e35–47. doi: 10.1161/CIRCRESAHA.118.313283
32. Mazurek S. Pyruvate kinase type M2: a key regulator of the metabolic budget system in tumor cells. Int J Biochem Cell Biol (2011) 43(7):969–80. doi: 10.1016/j.biocel.2010.02.005
Keywords: JAK2V617F, glycolysis, PKM1, STAT3, IL-6
Citation: Li R, Sun N, Chen X, Li X, Zhao J, Cheng W, Hua H, Fukatsu M, Mori H, Takahashi H, Ohkawara H, Fukami M, Okamoto M, Hamazaki Y, Zheng K, Yang J and Ikezoe T (2021) JAK2V617F Mutation Promoted IL-6 Production and Glycolysis via Mediating PKM1 Stabilization in Macrophages. Front. Immunol. 11:589048. doi: 10.3389/fimmu.2020.589048
Received: 01 September 2020; Accepted: 23 December 2020;
Published: 08 February 2021.
Edited by:
Joanna Cichy, Jagiellonian University, PolandReviewed by:
Jonas Samuel Jutzi, Harvard Medical School, United StatesChristoph Garbers, University Hospital Magdeburg, Germany
Copyright © 2021 Li, Sun, Chen, Li, Zhao, Cheng, Hua, Fukatsu, Mori, Takahashi, Ohkawara, Fukami, Okamoto, Hamazaki, Zheng, Yang and Ikezoe. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Takayuki Ikezoe, ikezoet@fmu.ac.jp; Jing Yang, jingyang@xzhmu.edu.cn
†These authors have contributed equally to this work