- 1Institute of Immunology, PLA, Third Military Medical University, Chongqing, China
- 2National Clinical Research Center of Kidney Diseases, Jinling Hospital, Nanjing, China
- 3Cancer Center, The General Hospital of Western Theater Command, Chengdu, China
- 4Chongqing Public Health Medical Center, Chongqing, China
Follicular helper T cells (TFH cells), known as the primary “helpers” of the germinal center (GC) reaction, promote the humoral immune response to defend against various pathogens. Under conditions of infection by different types of pathogens, many shared transcription factors (TFs), such as Bcl-6, TCF-1, and Maf, are selectively enriched in pathogen-specific TFH cells, orchestrating TFH cell differentiation and function. In addition, TFH cells also coexpress environmentally associated TFs as their conventional T cell counterparts (such as T-bet, GATA-3, or ROR-γt, which are expressed in Th1, Th2, or Th17 cells, respectively). These features likely indicate both the lineage-specificity and environmental adaption of the TFH cell responses. However, the extent to which the TFH cell response relies on these environmentally specific TFs is not completely understood. Here, we found that T-bet was specifically expressed in Type I TFH cells but not Type II TFH cells. While dispensable for the early fate commitment of TFH cells, T-bet was essential for the maintenance of differentiated TFH cells, promoting their proliferation, and inhibiting their apoptosis during acute viral infection. Microarray analysis showed both similarities and differences in transcriptome dependency on T-bet in TFH and TH1 cells, suggesting the distinctive role of T-bet in TFH cells. Collectively, our findings reveal an important and specific supporting role for T-bet in type I TFH cell response, which can help us gain a deeper understanding of TFH cell subsets.
Introduction
Because of the complexity and diversity of pathogens, organisms have developed highly organized and well-adapted immune systems to eliminate invaders. To defend against different microorganisms, the immune system elicits optimal responses according to the species of invader (1). For example, intracellular microbes induce type I immune response, which consists of IFN-γ-producing group1 innate lymphoid cell (ILC) lineages (including natural killer cells and ILC1s) (2–4), CD4+ type I T helper cells (TH1) (5, 6), and CD8+ type I cytotoxic T cells (TC1) mediated responses (7, 8); venoms or helminthes induce type II immune response, which includes IL-4-producing ILC2s (9–11), TH2 cells (6, 12), and TC2 cells (7, 13); and extracellular fungi or bacteria induce type III immune response, which comprises IL-17-producing ILC3s (14, 15), TH17 cells (16–18), and TC17 cells (19, 20). This phenomenon reflects the high plasticity and environmental dependency of immune cells.
As a CD4+ helper T cell subset specialized to “help” the germinal center (GC) reaction (21, 22), follicular helper T cells (TFH cells) have been reported to play important roles in type I immune response (23–27), type II immune response (28–32), and type III immune response (33–35). TFH cells express high levels of CXCR5, which is required for their localization in lymphoid follicles (21, 22, 36–40). In the light zone of GCs, they provide crucial signals to antigen-specific B cells and promote somatic hyper mutation, class switch recombination (CSR), and affinity maturation of GC B cells through cellular interactions and cytokine secretion (41–45). In addition, TFH cells also facilitate the differentiation of memory B cells and long-lived plasma cells from GC B cells (21, 22, 46).
TFH cells share similar differentiation processes during different types of immune responses; during the initiation phase of TFH cell differentiation, the expression of some TFs (such as Bcl-6, Ascl2, Maf, and TCF-1) is regulated in certain activated CD4+ T cells, which promotes CXCR5 expression (47–51). Next, CXCR5+Bcl-6+TFH precursor cells migrate to the T-B border zone, where they receive more differentiation signals from activated B cells (52). After this engagement, the reinforced expression of Bcl-6 regulates surface markers, which promote the migration of the TFH cells into GCs, where they provide helper signals to B cells (53, 54). Despite these similarities, TFH cells are also endowed with some unique characteristics for responding to distinct microenvironment associated with different types of microbial infection. Previous studies showed that TFH cells also express lineage-specific TFs like their conventional counterparts when defending against different types of pathogens, such as T-bet, GATA-3, or ROR-γt in type I, II, or III immune responses, respectively (28, 55–60). The production of IFN-γ, IL4, or IL17 driven by these specific TFs in TFH cells can help B cells switch to the optimal class of antibody to clear the microbes (25, 29, 33, 61–66). However, the extent to which TFH cells rely on these TFs for their differentiation or maintenance is not clear.
The transcription factor T-bet was originally discovered as a lineage marker of TH1 cells because it can establish TH1 differentiation and inhibit polarization of other CD4+ T cell subsets such as TH2 or TH17 cells (67–69). Later, it was also found to be extensively expressed by multiple different lymphocyte lineages during type I immune response, including both innate and adaptive immune cell subsets (70–72). For example, T-bet has been reported to promote the early differentiation and terminal maturation of NK cells (73–75). It has also been found that T-bet can promote IFN-γ production in ILCs and γδ T cells (4, 76, 77). Moreover, T-bet expressed by DCs can enhance their TH1-priming capacity (78). In NKT cells, T-bet can upregulate CD122 levels and promote survival (74, 79). In addition, T-bet is required for optimal terminal differentiation and granzyme B secretion in CD8+ T cells (80, 81). Moreover, T-bet expressed by B cells can promote the survival of memory B cells and enhance IgG2 switching (66, 82). Together, these facts highlight T-bet as the master regulator of type I immune response. In type I immune response, activated antigen-specific CD4+ T helper cells mainly differentiate into TH1 and TFH cell subsets (27, 83). Most studies have focused on the role of T-bet in TH1 differentiation and have generally considered T-bet to be a suppressor of TFH differentiation (72, 84–86). However, the exact role of T-bet in the TFH cell response is not well-understood.
In this study, using a combined conditional/inducible knockout system, we investigated the putative role of T-bet in regulating the response of virus-specific TFH cell in acute viral infection. We found the constitutive expression of T-bet in TFH cells during acute viral infection. A great reduction in the magnitude of the TFH cell response was observed when T-bet expression was deficient. Furthermore, microarray analysis showed significant differences in function- and proliferation-related genes between WT and Tbx21−/− TFH cells. In addition, TFH and TH1 cells showed different levels of T-bet dependency in their lineage-specific expression patterns. Thus, our findings demonstrate the crucial and specific role of T-bet in type I TFH cell responses, which suggests that modulation of T-bet expression in TFH cells may be a powerful therapeutic method for the treatment of infectious diseases and autoimmune diseases.
Materials and Methods
Mice and Treatment
C57BL/6J (CD45.1+ and CD45.2+), CD4cre transgenic, Ifng−/− and Tbx21fl/fl mice were purchased from Jackson Laboratory. ERT2cre transgenic mice were kindly provided by Yisong Wan (University of North Carolina). SMARTA (CD45.1+) mice were a kind gift from Rafi Ahmed (Emory University). All these strains had a C57BL/6J background. All mice were housed and bred under specific-pathogen-free (SPF) conditions. All mouse experiments were performed following the guidelines of the Institutional Animal Care and Use Committees of Army Medical University. All mice were infected/immunized at 6–10 weeks of age. Lymphocytic choriomeningitis virus (LCMV, Armstrong strain) was provided by Rafi Ahmed (Emory University). A total of 2 × 105 plaque-forming units of LCMV (Armstrong strain) were injected intraperitoneally to establish an acute viral infection model in mice. The Listeria monocytogenes-expressing LCMV-gp61-80 was created from vector strain 1. A total of 1 × 107 colony-forming units of recombinant bacteria were injected intravenously to establish a mouse bacterial infection model. NP-KLH (100 μg; N-5060-25; Biosearch Technology) was emulsified 1:1 with Aluminum hydroxide gel (Alum) (21645-51-2; InvivoGen) and was injected subcutaneously to establish a protein immunization model in mice. Tamoxifen (1 mg; T5648; Sigma-Aldrich) was diluted with sunflower oil and injected intraperitoneally into ERT2cre-Tbx21fl/fl or ERT2cre-Tbx21fl/fl mice to induce gene deletion at the indicated timepoints.
Flow Cytometry and Antibodies
Stained cells were analyzed by flow cytometry with a FACS Canto II flow cytometer (BD Bioscience). Flow cytometry data were analyzed with FlowJo software (Tree Star). LCMV-GP66 tetramer staining was described previously (51). CXCR5 staining has also been described previously (47). Surface staining was performed in PBS containing 2% fetal bovine serum (weight/volume). Staining for intracellular IgG2c, Bcl2, and IFN-γ was performed using a Cytofix/Cytoperm Fixation/Permeabilization Kit (554714; BD Bioscience). Staining for intranuclear TCF-1, T-bet and FOXP3 was performed with a Foxp3/Transcription Factor Staining Buffer Set (00-5523; eBioscience). For intracellular cytokine production analysis, before surface and intracellular staining, cells were stimulated with GP61-80 peptide for 5 h at 37°C, 5% CO2 in the presence of GolgiPlug (BD Bioscience), GolgiStop (BD Bioscience), and DNase I (Sigma-Aldrich). For in vivo incorporation of BrdU, mice were given BrdU (1.5 mg of BrdU in 0.5 ml of DPBS) intraperitoneally 3 h before staining. BrdU staining was performed with a BrdU Flow Kit (559619; BD Bioscience) according to the manufacturer's instructions. Annexin V staining was performed with an Annexin V Apoptosis Detection Kit I (559763; BD Bioscience) according to the manufacturer's instructions. The antibodies and reagents used in flow cytometry staining are listed in Table S1.
Enzyme-Linked Immunosorbent Assay
LCMV-specific IgG and IgG2c were titrated with LCMV lysates and the secondary antibodies HRP-conjugated goat anti-mouse IgG (1036-05; SouthernBiotech) and HRP-conjugated goat anti-mouse IgG2c (1078-05; SouthernBiotech) as previously described (87).
Adoptive Cell Transfer
To examine the LCMV-specific TFH cell response, 1 × 106 (for analysis before day 3 or after day 30) or 2 × 105 (for analysis between day 3 and day 30) sorted naïve or retrovirus-transduced CD45.1+ SMARTA cells (WT or Tbx21−/−) were adoptively transferred into naïve or infection-matched CD45.2+ mice (WT or Tbx21−/−) according to the requirements of the experiments. After being allowed to rest for one day, the cell-transferred hosts were infected intravenously with 2 × 106 plaque-forming units (for analysis at day 3 or earlier) or infected intraperitoneally with 2 × 105 plaque-forming units (for analysis at day 5 or later).
Mixed Bone Marrow Chimera
To determine the intrinsic role of T-bet, bone marrow cells collected from CD45.2+ Tbx21−/− mice and CD45.1+ WT mice were mixed at a ratio of 3:7 and transferred intravenously into lethally irradiated (5.5 Gy, twice) naïve WT CD45.1+ mice (5 × 106 cells/mouse). After at least 8 weeks of bone marrow reconstitution, the recipients were infected with LCMV.
Quantitative RT-PCR
To compare gene expression in LCMV-specific TH1 cells and TFH cells differentiated from naïve WT and Tbx21−/− SMARTA cells, SLAMhiCXCR5− and SLAMlowCXCR5+ SMARTA cells were sorted from recipient mice and directly lysed with TRIzol LS reagent (10296; Life Technologies). Total RNA was extracted with isopropyl ethanol and reverse-transcribed with a RevertAid H Minus First Strand cDNA Synthesis Kit (K1632; Thermo Scientific). Quantitative PCR of cDNA was carried out with a QuantiNova SYBR Green PCR Kit (208054; Qiagen) on a CFX96 Touch Real-Time System (Bio-Rad). The sequences of Tbx21 primers used in RT-qPCR are listed here: Tbx21 (F)-5′ CAATGTGACCCAGATGATCG 3′; Tbx21(R)-5′ CAATGTGACCCAGATGATCG 3′. Expression was calculated normalized to Hprt.
Microarray and Analysis
For the isolation of LCMV-specific TH1 cells and TFH cells, SLAMhiCXCR5− and SLAMlowCXCR5+ SMARTA cells were sorted from recipient mice adoptively transferred with WT or Tbx21−/− SMARTA cells at day 6 post LCMV infection. For the isolation of naïve CD4+ T cells, CD44−CD62L+CD4+ T cells were sorted from naïve WT C57BL/6J mice. The cells were sorted directly into TRIzol LS reagent (10296; Life Technologies). Total RNA was extracted with isopropyl ethanol and submitted to the CapitalBio Corporation for microarray analysis. Gene set enrichment analysis (GSEA) was performed as described previously (88). Clustering analysis was performed and heat maps were constructed using Cluster 3.0 with a hierarchical average linkage method, and the results were visualized using Java TreeView software. Pathway enrichment analysis was performed using KOBAS 3.0 (89).
Immunofluorescence Staining
Spleen tissues were snap frozen in O.C.T. compound (4583; SAKURA) and stored at −80°C until frozen sectioning. The tissues were cut into 10 μm-thick cryosections and fixed with ice-cold acetone. The sections were rehydrated and blocked with 5% rat serum and 3% BSA with 0.1% Tween and stained with fluorescent-labeled antibodies and reagents, including CD4 (RM4-5; BioLegend), IgD (11-26c.2a; BioLegend), GL7 (GL7; BD Bioscience), and DAPI (R37606; Invitrogen). Images were obtained with an EVOS FL Imaging System (ThermoFisher).
Statistical Analysis
Statistical analysis was performed with Prism 6.0. Differences between groups were analyzed with paired (for bone marrow chimera experiments) or unpaired two-tailed t-tests. A p-value <0.05 was considered significant.
Results
The Transcription Factor T-Bet Is Selectively Expressed in Type I but not Type II TFH Cells
Previously, it has been reported that T-bet is expressed in mouse TFH cells in LCMV infection model (27), which belongs to type I immune response. However, whether TFH cells express T-bet in other type I immune response models or in type II immune responses is not clear. Thus, we first examined the expression of T-bet in Listeria monocytogenes (LM) infection, NP-KLH immunization and LCMV infection models. Based on the expression of CD44 and CXCR5, FOXP3−CD4+ T cells were divided into three subsets, CD44+CXCR5+, CD44+CXCR5−, and CD44−CXCR5− cells, which were referred to as TFH, non-TFH and naïve CD4+ T cells, respectively (Figures 1A–C). At day 8 after immunization, we observed that TFH cells and non-TFH cells generated from the LCMV/LM infection model expressed much higher levels of T-bet than naïve CD4+ T cells (Figures 1D,E). In addition, we noticed that TFH cells expressed less T-bet than non-TFH cells in the LCMV/LM infection model (Figures 1D,E), which is consistent with published data (27). However, in the NP-KLH immunization model, there is nearly no detectable T-bet expression in both TFH cells and non-TFH cells (Figure 1F). These data demonstrated that T-bet is selectively expressed in TFH cells derived from type I rather than type II immune responses, suggesting that unlike common transcription factors such as TCF1 or Bcl6, T-bet may be an immune response type-dependent feature of TFH cells.
Figure 1. Transcription factor T-bet is selectively expressed in Type I but not Type II TFH cells. (A–F) WT C57BL/6 mice were infected with LCMV, LM, or immunized with NP-KLH. Lymphocytes of Spleen (for LCMV and LM infection) or draining lymph nodes (for NP-KLH immunization in alum) were isolated and analyzed for T-bet expression in TFH cells at day 8 post infection/immunization. (A–C) Representative flow cytometry of TFH cells (CD44+CXCR5+), Non-TFH cells (CD44+CXCR5−) and Naïve CD4+ T cells (CD44−CXCR5−) in LCMV (A), LM (B), infection or NP-KLH (C) immunization model. Numbers adjacent to outlined areas indicate percent of each subset in parent subset. (D–F) Representative Flow cytometry of T-bet expression in TFH cells, Non-TFH cells and Naïve CD4+ T cells (left) and the summary of T-bet expression by calculating the mean fluorescence intensity (MFI) of T-bet in each cell subsets (right) during LCMV (D), LM (E) infection, or NP-KLH (F) immunization in alum. (G,H) WT SMARTA cells (CD45.1+) were transferred into WT naïve C57BL/6 mice (CD45.2+) and the splenocytes were analyzed for T-bet expression of virus-specific TFH and TH1 cells at day3, 6, 10, 40, 90, and 160 post LCMV infection. (G) Representative Flow cytometry of T-bet expression in TFH cells (CD45.1+CXCR5+), TH1 cells (CD45.1+CXCR5−) and Naïve CD4+ T cells (CD45.2+CD44−). (H) Kinetics of T-bet expression of TFH and TH1 cells by calculating the MFI of T-bet in each subset followed by normalization to the T-bet MFI of Naïve CD4+ T cells. ns, not significant; *P < 0.05, **P < 0.01,****P < 0.0001 (unpaired two-tailed t-test). Data are representative of two independent experiments with 3–5 mice per group (error bars, SEM).
Next, focusing on the expression of T-bet in type I TFH cells, we investigated the expression kinetics of T-bet in LCMV-specific TFH cells using a SMARTA cell adoptive transfer system. SMARTA cells express LCMV-gp66-specific TCRs, so they can recognize and respond to LCMV and other LCMV-gp66 epitope-carrying microbes (90). We purified naïve ly5.1 SMARTA cells from spleen tissue, transferred them into naïve C57BL/6J recipient mice, and infected the recipient mice with the LCMV strain Armstrong. At different time points post infection, we measured the expression of T-bet in donor SMARTA TH1 and TFH cells. At day 3 post infection, T-bet expression in TFH and TH1 cells was ~2- and 4-fold higher, respectively, than that in naïve CD4+ T cells (Figures 1G,H). At day 6 post infection, T-bet expression in TFH and TH1 cells was upregulated to nearly 8- and 24-fold, respectively, compared to that in naïve CD4+ T cells (Figures 1G,H). At day 10 post infection, T-bet expression in TFH and TH1 cells had decreased back to levels ~6- and 13-fold higher than those in naïve CD4+ T cells (Figures 1G,H). At day 40, 90, and 160 post infection, T-bet expression in TFH and TH1 cells had further decreased and remained ~4- and 6-fold higher, respectively, than that in naive CD4+ T cells (Figures 1G,H). Taken together, these data indicate that TFH and TH1 cells share a dynamic similarity in their T-bet expression patterns: T-bet expression sharply increases in the early effect phase, gradually falls back in the contraction phase, and is stably maintained at a certain level in the memory phase. Meanwhile, consistently lower levels of T-bet were observed in TFH cells than in TH1 cells throughout the entire response.
T-Bet Is Required for TFH Cell Expansion During Acute Viral Infection
To investigate whether T-bet is required for optimal TFH cell responses during acute viral infection, we generated a CD4cre-Tbx21fl/fl strain of mice (called Tbx21−/− mice here) by crossing transgenic CD4cre mice with Tbx21fl/fl mice to selectively knock out the Tbx21 gene (encoding the T-bet protein) in T cells. These Tbx21−/− mice showed normal T cell development in vivo (Figures S1A,B). At day 8 post LCMV Armstrong infection, the RT-qPCR and flow cytometry results showed that TH1 and TFH cells from Tbx21−/− mice did not express T-bet (Figures 2A,B). In addition, we observed a significant decrease in the frequency and number of polyclonal CD44+CXCR5+ TFH cells in the spleens of Tbx21−/− mice (Figure 2C). To determine whether this phenotype was caused by deficient clonal expansion or abnormal differentiation, we used the gp66 tetramer to measure antigen-specific CD4+ T cells. The results showed that fewer gp66-specific CD4+ T cells were present inTbx21−/− mice (Figure 2D), indicating that clonal expansion of gp66-specific Tbx21−/− CD4+ T cells was heavily affected. Consistent with the decreased number of gp66-specific CD4+ T cells, the number of gp66-specific TFH cells was also greatly decreased in Tbx21−/− mice (Figure 2E). Furthermore, we observed a mild increase in the frequency of gp66-specific Tbx21−/− TFH cells (Figure 2E), which was consistent with a report that T-bet inhibits TFH cell differentiation in vitro (84). Together, these data suggest that T-bet is required for TFH cell response mainly by promoting clonal expansion during acute viral infection.
Figure 2. T-bet is required for TFH cell expansion during acute viral infection. WT and Tbx21−/− mice were infected with LCMV. Splenocytes were isolated and analyzed by flow cytometry and RT-qPCR at day 8 post infection. (A) Summary of T-bet expression (showed as MFI) of TFH cells (CD44+CXCR5+) and TH1 cells (CD44+CXCR5−) in WT and Tbx21−/− mice. (B) T-bet gene expression relative to Hprt mRNA of sorted TFH and TH1 cells from WT and Tbx21−/− mice was assessed by RT-qPCR. (C–E) Representative Flow cytometry of TFH cells (CD44+CXCR5+) (C), tetramer-positive CD4+ T cells (D), and tetramer-positive TFH cells (E) in WT and Tbx21−/− mice (left), and the summary of percentages and number of these cell subsets (right). Numbers adjacent to outlined areas indicate percent of each cell subset in parent subset. ns, not significant; *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 (unpaired two-tailed t-test). Data are representative of three independent experiments with 3–5 mice per group (error bars, SEM).
Optimal Germinal Center Response Requires T-Bet Expression in TFH Cells
Based on the critical role of TFH cells in “helping” GC response, we next investigated whether T-bet deficiency in TFH cells would influence the germinal center response. At day 8 post LCMV Armstrong infection, we observed severely affected GC formation in the spleens of Tbx21−/− mice (Figure 3A). Loss of T-bet expression in TFH cells strongly reduced the frequency of GC B cells and plasma cells (Figures 3B,C). In addition, Tbx21−/− mice showed less IgG2c class switching in plasma cells (Figure 3D) as well as much lower LCMV-specific IgG and subtype IgG2c titers in serum (Figures 3E,F) than WT mice. These data further verified the vital role of T-bet in promoting the TFH cell response and antibody IgG2 class switching.
Figure 3. Optimal germinal center response requires T-bet expressed in TFH cell. WT and Tbx21−/− mice were infected with LCMV. Spleens were harvested at day 8 post infection. (A) Representative Immunofluorescent staining of splenic B cell follicular with GL-7 (green), anti-IgD (blue), and anti-CD4 (red). (B,C) Representative flow cytometry of GC B cells (FAS+PNA+) (B) and plasma cells (CD138+B220low) (C) (left) with its percentages (right) in WT and Tbx21−/− mice. (D) Representative flow cytometry of IgG2c+ plasma cells (left) with its percentages (right) in WT and Tbx21−/− mice. Numbers adjacent to outlined areas indicate percent of each cell subset in parent subset. (E,F) Serum of WT and Tbx21−/− mice were collected and tested for anti-LCMV IgG (E) and anti-LCMV IgG2c (F) by enzyme-linked immunosorbent assay (ELISA). *P < 0.05, ***P < 0.001, ****P < 0.0001 (unpaired two-tailed t-test). Data are representative of two independent experiments with 3–5 mice per group (error bars, SEM).
T-Bet Is not Required for Type II TFH Cell Response
The observation of a compromised type I TFH response in Tbx21−/− mice during acute viral infection led us to investigate whether T-bet plays an important role in type II TFH cell response. Thus, we tested the TFH and GC responses of Tbx21−/− mice in a protein immunization model. At day 8 post NP-KLH immunization, we observed similar frequencies and numbers of TFH cells in Tbx21−/− mice and control mice (Figure S2A). In addition, we did not find any reductions in the frequency or number of GC B cells as well as plasma cells in Tbx21−/− mice (Figures S2B,C). These results are consistent with the observation that T-bet is not expressed in CD4+ T cells during type II immune response. Together with the crucial role of T-bet in regulating type I TFH cell response during acute viral infection, these results confirm that T-bet is an environmentally specific regulator of type I TFH cell response.
T-Bet Promotes TFH Cell Expansion in a T Cell Intrinsic Manner
In the CD4cre-induced Tbx21 knockout system, both CD4+ and CD8+ T cells lost their capacity to express T-bet. In addition, a deficient GC response might reciprocally amplify the impairment of the TFH cell response in Tbx21−/− mice. To further clarify the cell-intrinsic role of T-bet in regulating the TFH cell response, we set up bone marrow chimeras by reconstituting lethally irradiated WT (ly5.1+) recipient mice with a 3:7 ratio mixture of bone marrow cells from Tbx21−/− (ly5.2+) and WT (ly5.1+) donor mice, respectively (Figure 4A). Chimera mice were infected with LCMV Armstrong after successful bone marrow reconstitution (Figures S3A,B). At day 8 post infection, we still observed a largely decreased frequency of gp66-specific CD4+ T cells and polyclonal TFH cells in Tbx21−/− (ly5.2+) mice compared to control mice (Figures 4B,C). Similar to what we found in Tbx21−/− mice, the frequency of gp66-specific TFH cells was slightly increased in Tbx21−/− cells of chimera mice (Figure 4D). Moreover, the Tbx21−/−: WT ratios of polyclonal TFH cells, gp66-specific CD4+ T cells and gp66-specific TFH cells were markedly decreased relative to that of total CD4+ T cells (Figures 4B–D). These data confirmed the intrinsic role of T-bet in regulating the TFH cell response during acute viral infection.
Figure 4. T-bet promotes TFH cell expansion in a T cell intrinsic manner. (A) Setup of Bone marrow chimera mice. Bone marrow cells from WT (CD45.1+) and Tbx21−/− (CD45.2+) mice were mixed in 7:3 ratio and adoptively transferred into lethally irradiated recipient WT (CD45.1+) mice. After bone marrow reconstitution, recipients were infected with LCMV (Armstrong strain). Spleens were harvested at day8 post infection. (B–D) Representative flow cytometry of TFH cells (CD44+CXCR5+) (B), tetramer-positive CD4+ T cells (C) and tetramer-positive TFH cells (D) in WT and Tbx21−/−. Bone marrow derived cells (left), and the summary of percentages and ratio (by comparing the T-bet−/−: WT ratio in target cells and total CD4+ T cells) of these cell subsets (right). Numbers adjacent to outlined areas indicate percent of each cell subset in parent subset. **P < 0.01, ***P < 0.001, ****P < 0.0001 (paired two-tailed t-test). Data are representative of two independent experiments with 6 or 9 chimera mice.
T-Bet Promotes TFH Cell Maintenance by Regulating Proliferation and Apoptosis
It was clear that the deficiency of the TFH cell response in Tbx21−/− mice was mainly caused by the greatly reduced magnitude of the TFH cell response. To investigate the kinetics of virus-specific TFH cell expansion in Tbx21−/− mice, we transferred the same number of WT or Tbx21−/− SMARTA cells into naïve recipient mice and then infected host mice with LCMV Armstrong. From day 2 post infection, we detected a continuous slightly higher frequency of TFH cells in the Tbx21−/− SMARTA group than in the WT group (Figures 5A,B). At day 2 and day 5 post infection, we did not observe any differences in the numbers of Tbx21−/− and WT SMARTA TFH cells (Figures 5A,B). To our surprise, the number of Tbx21−/− SMARTA TFH cells decreased sharply at day8 post infection (Figures 5A,B). Besides, the reduction in the virus-specific TFH cell population in Tbx21−/− mice might not have been caused by impaired early activation of CD4+ T cells (Figures S4A,B). Taken together, the results suggest the possibility that the loss of TFH cells in Tbx21−/− mice was mainly caused by reduced maintenance at the late phase of the anti-viral immune response.
Figure 5. T-bet promotes TFH cell maintenance. (A,B) WT or Tbx21−/− SMARTA cells (CD45.1+) were transferred into naïve WT mice (CD45.2+), which would be infected with LCMV. Splenocytes were isolated and analyzed for TFH response at day2, 5 and 8 post infection. Flow cytometry of TFH cells derived from transferred SMARTA cells (A), and the summary of the percentages and numbers of those cells (B). (C,D) WT and iT-bet−/− mice were treated with Tamoxifen for 3 days before or after LCMV infection. Spleens were harvested at day 9 post infection. (C) Flow cytometry of TFH cells (left) with its number (right) in mice treat with Tamoxifen before infection (day−3 to−1). (D) Flow cytometry of TFH cells (left) with its number (right) in mice treat with Tamoxifen after infection (day 5–7). (E) Setup of CD4+ T cell transfer experiment. CD4+ T cells were purified form spleens of iT-bet−/− mice (CD45.2+) without Tamoxifen treatment and adoptively transferred into infection-matched recipient mice (CD45.1+). At day 8–10 post infection, recipient mice were treated with Tamoxifen or vehicle. Spleens were harvested at day 14 post infection. (F) Representative flow cytometry of TFH cells (left) and summary of TFH cell number (right) in Tamoxifen treated mice (iTbx21−/−) and vehicle treated mice (WT) was showed here. (G) Setup of TFH co-transfer experiments. WT SMARTA cells (CD45.1+CD45.2+) and Tbx21−/− SMARTA cells (CD45.1−CD45.2+) were transferred into naïve B6 mice (CD45.1+CD45.2−) separately before infecting recipients with LCMV. At day 6 post infection, WT SMARTA TFH cells and Tbx21−/− SMARTA TFH cells were FACS sorted and mixed (about 1:1), then co-transferred into infection-matched B6 mice (CD45.1+CD45.2−). Spleens were harvested and analyzed for transferred TFH cells at day 9 post infection. (H) Flow cytometry of sorted SMARTA TFH cells before (day 6 p.i.) and after (day 9 p.i.) co-transfer (left). The percentage of Tbx21−/− SMARTA TFH cells in total donor TFH cells and the number of WT and Tbx21−/− SMARTA TFH cells at day 9 post infection were summarized (right). Numbers adjacent to outlined areas in (A,C,D,F,H) indicate percent of each cell subset in parent subset. ns, not significant; *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 (unpaired (B–D,F) or paired (H) two-tailed t-test). Data in (C,D,F) are representative of two independent experiments with 3–5 mice per group (error bars, SEM).
To more carefully investigate the influence of T-bet on TFH cell maintenance at the late phase of infection, we generated ERT2cre-Tbx21fl/fl mice (iTbx21−/−) by crossing Tbx21fl/fl mice with ERT2cre transgenic mice, in which Tbx21 gene knockout could be induced by tamoxifen treatment. We treated mice with tamoxifen at 1–3 days before or 5–7 days after LCMV infection to induce T-bet deletion before or after TFH cell commitment, respectively. Under both of these circumstances, we observed a lower abundance of TFH cells in iTbx21−/− mice at day 9 post infection (Figures 5C,D). In addition, iTbx21−/− CD4+ T cells were purified and adoptively transferred into infection-matched recipient mice at day 7 post LCMV infection (Figure 5E). After 3 days of tamoxifen or vehicle administration (days 8–10), we observed a decreased number of donor TFH cells in mice treated with tamoxifen than in control mice at day 14 post infection (Figure 5F). These results suggest that T-bet is required for the TFH cell response even after TFH commitment.
Furthermore, to investigate the sharp decrease in the number of Tbx21−/− SMARTA TFH cells at the late effector phase, we sorted WT and Tbx21−/− SMARTA TFH cells from recipient mice at day 6 post LCMV Armstrong infection and adoptively transferred a 1:1 ratio mixture of WT and Tbx21−/− SMARTA TFH cells into infection-matched mice (Figure 5G). At day 9 post infection, we observed that the ratio of WT and Tbx21−/− SMARTA TFH cells had changed to ~4:1 (Figure 5H). Taken together, these results indicated that intrinsic expression of T-bet is essential for TFH maintenance at the late effector phase.
Next, to gain insight into the reason for the reduction in the TFH cell population, we measured the proliferation and apoptosis of SMARTA TFH cells. At day 2 post infection, we observed even higher proliferation and expression of the survival marker Bcl2 but comparable apoptosis in Tbx21−/− SMARTA TFH cells compared to WT cells (Figures 6A–C). However, at day 5 post infection, the proliferation rate of Tbx21−/− SMARTA TFH cells dropped quickly to a significantly lower level than that of WT cells (Figure 6A), which was consistent with the lower Bcl2 expression in Tbx21−/− SMARTA TFH cells (Figure 6C). At day 8 post infection, higher apoptosis rates and lower Bcl2 expression were found in Tbx21−/− SMARTA TFH cells than in WT cells (Figures 6B,C). These results suggested that T-bet controls the TFH cell maintenance ultimately by promoting proliferation at the mid phase and inhibiting apoptosis at the late effector phase.
Figure 6. T-bet regulates TFH cell proliferation and apoptosis. WT or Tbx21−/− SMARTA cells (CD45.1+) were transferred into naïve WT mice (CD45.2+), which would be infected with LCMV. Splenocytes were isolated and analyzed for TFH response at day 2, 5, and 8 post infection. (A) Flow cytometry of BrdU incorporation in TFH cells (left), and summary of BrdU-incorporated cells percentage (right). Numbers adjacent to outlined areas indicate percent of each cell subset in parent subset. (B) Flow cytometry of Annexin-V and 7-AAD staining in TFH cells (left), and summary of Annexin-V+7-AAD− percentage (right). Numbers in quadrants indicate percent cells in each. (C) Flow cytometry of BCL-2 expression in TFH cells (left), and the summary of Bcl-2 expression (showed as MFI) of TFH cells. ns, not significant; **P < 0.01, ***P < 0.001, ****P < 0.0001 (unpaired two-tailed t-test). Data are representative of two independent experiments with 3–5 mice per group (error bars, SEM).
T-Bet Dependency of the TFH and TH1 Cell Transcriptomes
To investigate the molecular mechanisms regulated by T-bet in TFH and TH1 cell response, we sorted WT and Tbx21−/− SMARTA TFH and TH1 cells from recipient mice at day 6 post LCMV Armstrong infection after adoptive transfer, as well as naïve mice-derived CD4+ T cells, for gene expression profile analysis. The gene expression patterns differed greatly between WT and Tbx21−/− cell population at the genome-wide level (Figures 7A,B). We observed 822 upregulated and 899 downregulated genes in Tbx21−/− TH1 cells relative to WT TH1 cells (Figure S5A). Accordingly, we identified 151 upregulated and 367 downregulated genes in Tbx21−/− TFH cells (Figure S5B). Among these differentially expressed genes, 103-up and 223-down regulated genes were shared by Tbx21−/− TFH and TH1 cells (Figure 7C). Besides, PANTHER pathway enrichment analysis of the changed genes in Tbx21−/− TFH and TH1 cells also showed similar enrichment in many important pathways, such as DNA replication, Apoptosis, and Interferon-gamma signaling pathway, which account for the impaired maintenance of Tbx21−/− TFH cells (Figures S5C,D). To further figure out the downstream factors involved in the maintenance of differentiated TFH and TH1 cells controlled by T-bet, we did the gene set enrichment analysis focusing on cell proliferation and survival and found a reduction to a similar extent in an array of proliferation and survival relevant genes, such as Ccna2, Ccnb2, Aurkb, E2f1, E2f7, and E2f8 in both Tbx21−/− TH1 and TFH cells compared to their WT counterpart (Figure 7D), highlighting the shared regulatory pathway important for both TFH and TH1 proliferation and survival that is likely imprinted by the same Type-I microenvironment (58).
Figure 7. T-bet dependency of the TFH and TH1 cell transcriptomes. (A) Clustering of genes based on their expression in WT TH1, WT TFH, Tbx21−/− TH1, and Tbx21−/− TFH cells. (B) Principal component analysis (PCA) of genes in the five cell populations listed in (A). (C) Venn diagram comparing the differentially expressed genes in TH1 and TFH cells. The red numbers represent upregulated genes, and the green numbers represent downregulated genes. (D) Heat map of genes related to cell proliferation in WT TH1, Tbx21−/− TH1, WT TFH, and Tbx21−/− TFH cells. (E) GSEA of WT and Tbx21−/− TH1 genes showing gene sets upregulated in TH1 cells relative to TFH cells (left) and gene sets upregulated in TFH cells relative to TH1 cells (right). (F) GSEA of WT and Tbx21−/− TFH genes showing gene sets upregulated in TH1 cells relative to TFH cells (left) and gene sets upregulated in TFH cells relative to TH1 cells (right). The gene sets used in GSEA include published data (GEO accession code GSE21379) and our microarray data (GEO accession code GSE122931). (G) Heat map of genes related to T cell differentiation and function in WT TH1, Tbx21−/− TH1, WT TFH, and Tbx21−/− TFH cells.
Despite these similarities, bioinformatics analysis of the microarray data also observed many differences in the dependency of T-bet in TFH and TH1 cell development. Principal component analysis (PCA) showed that after T-bet deletion, the gene expression profile of Tbx21−/− TH1 cells shifted toward that of TFH cells to some extent (Figure 7B). Then, we selected sets of genes that were upregulated in TFH cells (TFH-UP) compared with non-TFH cells or upregulated in TH1 cells (TH1-UP) compared with TFH cells based on both published data (91) (GEO accession code GSE21379) and our microarray data (GEO accession code GSE122931). Gene set enrichment analysis (GSEA) of TH1 cells showed that the TH1-UP gene set was enriched in WT TH1 cells, whereas the TFH-UP gene set was enriched in Tbx21−/− TH1 cells, which suggested to some extent that the gene-expression pattern of Tbx21−/− TH1 cells lost TH1 signature and change to TFH signature (Figure 7E). Similar to that of TH1 cells, GSEA of TFH cells showed that the TH1-UP gene set was enriched in WT TFH cells (Figure 7F). However, different from what we have observed in TH1 cells, genes in the TFH-UP gene set were enriched in WT TFH cells much more than in Tbx21−/− TFH cells, which suggested that the gene-expression pattern of Tbx21−/− TFH cells lost TH1 signature but also lost part of its TFH signature (Figure 7F). Through further assessed 48 genes known to be associated with T cell differentiation or function, we found a distinct T-bet dependency of the TFH and TH1 cell gene-expression pattern (Figure 7G). In TH1 cells, We found that the transcription levels of a set of important Th1 lineage associated genes like Il12rb2 (92) Cxcr3, Prdm1, Gzmb, Gzmk, Gzma, Cx3cr1, and Ifng were significantly decreased in Tbx21−/− Th1 cells than in WT TH1 cells. We also observed Prdm1 expression was lower in Tbx21−/− TH1 cells than in WT TH1 cells. In addition, the expression of Foxp3, Gata3, and Rorc, which are essential for Treg, TH2, and TH17 cell differentiation, respectively, was higher in Tbx21−/− TH1 cells than in WT TH1 cells (Figure 7G). On the other hand, the abundances of TFH lineage-specification associated genes, including Tox2, Id3, Bhlhe40, Il6st, Il6ra, and Tcf7 were much increased in Tbx21−/− TFH cells (Figure 7G), at least in part explaining the differential role of T-bet in regulating the program of TH1 and TFH differentiation. We compared the expression of Batf4 (93) and Irf4 in WT and Tbx21−/− TFH cells. As expected, we observed the dramatically decrease expression of Batf and Irf4 in Tbx21−/− TFH cells. Besides, the significant lower expression of Icos (94) was also detected in Tbx21−/− TFH cells compared to that in WT counterparts (Figure 7G). Whereas, expression of other TFH cell-relevant genes like Prdm1, Bcl6, and Cxcr5 was not significantly influenced by T-bet deletion (Figure 7G). The imbalanced impact of T-bet deletion on TH1 and TFH cells may interpret the mildly higher frequency of TFH cells than TH1 cells in antigen-specific CD4+ T cells that we observed (Figures 2E, 4D, 5A). Together, these data indicate both similarities and differences in transcriptome dependency on T-bet in TFH and TH1 cells during acute viral infection.
IFN-γ as a Candidate Downstream Target of T-Bet in Regulating TFH Cell Expansion
As a direct target of T-bet, it has been reported that IFN-γ could promote clonal expansion and survival of CD4+ T cells (95). Microarray analysis also showed that the interferon-gamma signaling pathway was severely impaired in TFH cells after T-bet deletion. In addition, we observed largely decreased IFN-γ production in ex vivo recalled Tbx21−/− TFH cells (Figures 8A,B). To investigate if IFN-γ could regulate TFH cell expansion, we compared the TFH cell responses in Tbx21−/− and Ifng−/− mice during LCMV Armstrong infection. At day 8 post infection, we found that the frequency and number of TFH cells were decreased in both Tbx21−/− and Ifng−/− mice (Figure 8C). In addition, T-bet expression was not affected by IFN-γ deletion (Figure 8D). Although not sufficient, these results suggest that IFN-γ might be a candidate target of T-bet in regulating TFH cell expansion.
Figure 8. IFN-γ as a candidate downstream target of T-bet in regulating TFH cell expansion. (A) Flow cytometry of IFN-γ production in TFH cells (left), with the summary of IFN-γ percentages (right) derived from WT and Tbx21−/− mice at day 8 post LCMV infection. (B) Flow cytometry of IFN-γ production in WT and Tbx21−/− TFH cells (left), with the summary of IFN-γ percentages (right) derived from same bone marrow chimera mice at day 8 post LCMV infection. (C,D) WT, Tbx21−/− and Ifng−/− mice were infected with LCMV. Spleens were harvested at day 8 post infection. (C) Flow cytometry of TFH cells (left) with the summary of percentages and numbers (right) in these mice. (D) Flow cytometry of T-bet expression (left), and the summary of T-bet expression (showed as MFI) in TFH cells of WT, Tbx21−/−, and Ifng−/− mice. Numbers adjacent to outlined areas in (A–C) indicate percent of each cell subset in parent subset. ns, not significant; *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001 (unpaired two-tailed t-test). Data are representative of two independent experiments with 3–5 mice per group (error bars, SEM).
Discussion
Follicular helper T cells (TFH cells) play critical roles in type I, II, and III immune responses, but how TFH cells adapt to different environments has remained largely unknown. In this study, we identified a key link between the transcription factor T-bet and type I TFH cell response during acute viral infection. We found that T-bet is specifically expressed in TFH cells originating from type I but not type II immune response. Tbx21−/− mice exhibited significant deficiency in the TFH cell response during acute viral infection. We observed a greatly decreased magnitude of the TFH cell response in Tbx21−/− mice compared to WT mice, although a slightly increased ratio of TFH cells was observed. Based on these results, we concluded that T-bet is required for optimal type I TFH cell response.
In addition to LCMV infection, type I TFH cells have also been discovered as TH1-biased TFH cells in simian immunodeficiency virus infection and TH1-polarized TFH cells in malaria infection (96, 97). Type I TFH cells not only express Bcl6, CXCR5, IL21, and PD-1 but also coexpress CXCR3 and IFN-γ. The secreted IFN-γ can help antibodies from B cells switch to the IgG2a/c class, which is essential for efficient elimination of viruses and other pathogens. Similar to the case in TH1 cells, the expression of TH1-associated molecules in type I TFH cells is also induced by the transcription factor T-bet. In addition, our research demonstrated the crucial role of T-bet in promoting TFH cell proliferation and maintenance. Based on the specialized role of type I TFH cells in defending against intracellular pathogens, we propose that this Th1-like effector TFH population be named TFH1.
During TH1 differentiation, it has been reported that T-bet attenuates the TFH cell-like phenotype in the late phase of TH1 specification by repressing the expression of Bcl6 and other molecules associated with TFH cell development (84). In addition, T-bet has been found to inhibit Tcf7 expression by directly binding with the Tcf7 gene promoter and suppress Bcl6 function by physically interacting with the Bcl6 protein (85, 86, 98). These in vitro results suggest that T-bet uses multiple mechanisms to inhibit TFH differentiation. However, evidence supporting the role of T-bet in regulating the TFH phenotype is not sufficient at the in vivo level. In our study, we clearly showed that the maintenance of TFH at the later effector phase is sharply impaired after T-bet deletion even though mildly increased early TFH differentiation was observed. Additionally, the constitutive expression of T-bet in type I TFH cells may suppress early TFH generation but sustain the clonal expansion of TFH cells at the late stage. Thorough understanding of the distinct role of T-bet in type I TFH cells at different stages necessitates further investigation in the future.
Notably, we did not observe any T-bet expression in type II TFH cells during protein immunization. Accordingly, Tbx21−/− mice showed normal TFH and GC responses during protein immunization. These results remind us that, unlike Bcl6 or Blimp1, the transcription factor T-bet is not a fundamental regulator of “all-weather” TFH cell responses under natural conditions. In other words, T-bet is a type I TFH cell-specific regulator, suggesting that diverse transcription factors are required for optimal TFH cell responses in different environments.
In addition, our microarray analysis results showed that there are not only many similar but also many different important changes in the TFH and TH1 cell transcriptomes that occur in a T-bet-dependent manner. On the one hand, Tbx21−/− TFH cells share many altered genes with TH1 cells, including genes enriched in signaling pathways involved in DNA replication, apoptosis and interferon-gamma signaling. On the other hand, many TFH differentiation-related genes were altered in different directions and to different degrees in Tbx21−/− TH1 and TFH cells. Four possible mechanisms might be involved in this scenario. First, some genes are regulated by T-bet in a redundant way, which means that the regulatory role of T-bet may be unnecessary if these genes have already been up- or downregulated by other transcription factors. Second, differences in chromatin accessibility between TFH and TH1 cells would lead to differences in the binding affinity of T-bet. Third, by interacting with different transcription factors, T-bet could differentially regulate gene expression in TFH and TH1 cells. Fourth, the post transcriptional modification of T-bet is different in TFH and TH1 cells, which might result in different or even opposite regulatory functions at the same gene loci. Further studies are needed to explore the exact mechanism underlying the contradictory effects of T-bet in regulating the development of TH1 and TFH cells.
Overall, this study revealed that T-bet, although slightly inhibiting TFH differentiation, mainly supports type I TFH cell response by promoting cell proliferation and apoptotic intervention to maintain the TFH cell response at the late effector phase during acute viral infection. These findings provide important insights into the transcription factor-mediated regulation of the environmental suitability of TFH cells.
Data Availability
The datasets generated for this study can be found in Gene Expression Omnibus, GSE122931.
Ethics Statement
All mouse experiments were performed following the guidelines of the Institutional Animal Care and Use Committees (IACUCs) of Army Medical University. The protocols were approved by the IACUCs.
Author Contributions
PW, YC, YuW, and LY designed and supervised the study. PW, YW, LuX, MX, JW, QH, BL, XC, LiX, SY, and YH performed experiments. PW wrote the manuscript with YW. QB and RH helped with analysis. All authors read and approved the manuscript.
Funding
This study was supported by grants from the National Natural Science Foundation of China (No. 31825011 to LY; No. 31700774 to LiX) and Open research fund of State Key Laboratory of Veterinary Biotechnology (SKLVBF2018XX to LY) and the National Key Research Development Plan (No. 2016YFA0502202 to LY) and the National Science and Technology Major Project of China during the 13th Five-year Plan (2018ZX10302104 to YC).
Conflict of Interest Statement
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.
Acknowledgments
We thank Yisong Wan (University of North Carolina) for kindly providing us ERT2cre transgenic mice; We thank Rafi Ahmed (Emory University) for generous providing us LCMV Armstrong virus and SMARTA transgenic mice. We also thank CapitalBio Corporation (Beijing, China) for performing microarray experiments and data analysis, and we thank the core facility center of Army Medical University for cell sorting.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2019.00606/full#supplementary-material
References
1. Annunziato F, Romagnani C, Romagnani S. The 3 major types of innate and adaptive cell-mediated effector immunity. J Allergy Clin Immunol. (2015) 135:626–35. doi: 10.1016/j.jaci.2014.11.001
2. Klose CSN, Flach M, Mohle L, Rogell L, Hoyler T, Ebert K, et al. Differentiation of type 1 ILCs from a common progenitor to all helper-like innate lymphoid cell lineages. Cell. (2014) 157:340–56. doi: 10.1016/j.cell.2014.03.030
3. Sun JC, Lanier LL. NK cell development, homeostasis and function: parallels with CD8(+) T cells. Nat Rev Immunol. (2011) 11:645–57. doi: 10.1038/nri3044
4. Bernink JH, Peters CP, Munneke M, te Velde AA, Meijer SL, Weijer K, et al. Human type 1 innate lymphoid cells accumulate in inflamed mucosal tissues. Nat Immunol. (2013) 14:221–9. doi: 10.1038/ni.2534
5. Mosmann TR, Cherwinski H, Bond MW, Giedlin MA, Coffman RL. Two types of murine helper T cell clone. I. Definition according to profiles of lymphokine activities and secreted proteins. J Immunol. (1986) 136:2348–57.
6. Farrar JD, Asnagli H, Murphy KM. T helper subset development: roles of instruction, selection, and transcription. J Clin Invest. (2002) 109:431–5. doi: 10.1172/JCI0215093
7. Mosmann TR, Sad S, Krishnan L, Wegmann TG, Guilbert LJ, Belosevic M. Differentiation of subsets of CD4+ and CD8+ T cells. Ciba Found Symp. (1995) 195:42–50. discussion-50–4.
8. Carter LL, Murphy KM. Lineage-specific requirement for signal transducer and activator of transcription (Stat)4 in interferon gamma production from CD4(+) versus CD8(+) T cells. J Exp Med. (1999) 189:1355–60. doi: 10.1084/jem.189.8.1355
9. Price AE, Liang HE, Sullivan BM, Reinhardt RL, Eisley CJ, Erle DJ, et al. Systemically dispersed innate IL-13-expressing cells in type 2 immunity. Proc Natl Acad Sci USA. (2010) 107:11489–94. doi: 10.1073/pnas.1003988107
10. Fort MM, Cheung J, Yen D, Li J, Zurawski SM, Lo S, et al. IL-25 induces IL-4, IL-5, and IL-13 and Th2-associated pathologies in vivo. Immunity. (2001) 15:985–95. doi: 10.1016/S1074–7613(01)00243–6
11. Walker JA, McKenzie AN. Development and function of group 2 innate lymphoid cells. Curr Opin Immunol. (2013) 25:148–55. doi: 10.1016/j.coi.2013.02.010
12. Romagnani S. Lymphokine production by human T cells in disease states. Annu Rev Immunol. (1994) 12:227–57. doi: 10.1146/annurev.iy.12.040194.001303
13. Maggi E, Giudizi MG, Biagiotti R, Annunziato F, Manetti R, Piccinni MP, et al. Th2-like CD8+ T cells showing B cell helper function and reduced cytolytic activity in human immunodeficiency virus type 1 infection. J Exp Med. (1994) 180:489–95. doi: 10.1084/jem.180.2.489
14. Spits H, Artis D, Colonna M, Diefenbach A, Di Santo JP, Eberl G, et al. Innate lymphoid cells–a proposal for uniform nomenclature. Nat Rev Immunol. (2013) 13:145–9. doi: 10.1038/nri3365
15. Eberl G, Marmon S, Sunshine MJ, Rennert PD, Choi Y, Littman DR. An essential function for the nuclear receptor RORgamma(t) in the generation of fetal lymphoid tissue inducer cells. Nat Immunol. (2004) 5:64–73. doi: 10.1038/ni1022
16. Harrington LE, Hatton RD, Mangan PR, Turner H, Murphy TL, Murphy KM, et al. Interleukin 17-producing CD4+ effector T cells develop via a lineage distinct from the T helper type 1 and 2 lineages. Nat Immunol. (2005) 6:1123–32. doi: 10.1038/ni1254
17. Acosta-Rodriguez EV, Rivino L, Geginat J, Jarrossay D, Gattorno M, Lanzavecchia A, et al. Surface phenotype and antigenic specificity of human interleukin 17-producing T helper memory cells. Nat Immunol. (2007) 8:639–46. doi: 10.1038/ni1467
18. Annunziato F, Cosmi L, Santarlasci V, Maggi L, Liotta F, Mazzinghi B, et al. Phenotypic and functional features of human Th17 cells. J Exp Med. (2007) 204:1849–61. doi: 10.1084/jem.20070663
19. Intlekofer AM, Banerjee A, Takemoto N, Gordon SM, Dejong CS, Shin H, et al. Anomalous type 17 response to viral infection by CD8+ T cells lacking T-bet and eomesodermin. Science. (2008) 321:408–11. doi: 10.1126/science.1159806
20. Yen HR, Harris TJ, Wada S, Grosso JF, Getnet D, Goldberg MV, et al. Tc17 CD8 T cells: functional plasticity and subset diversity. J Immunol. (2009) 183:7161–8. doi: 10.4049/jimmunol.0900368
21. Breitfeld D, Ohl L, Kremmer E, Ellwart J, Sallusto F, Lipp M, et al. Follicular B helper T cells express CXC chemokine receptor 5, localize to B cell follicles, and support immunoglobulin production. J Exp Med. (2000) 192:1545–52. doi: 10.1084/jem.192.11.1545
22. Schaerli P, Willimann K, Lang AB, Lipp M, Loetscher P, Moser B. CXC chemokine receptor 5 expression defines follicular homing T cells with B cell helper function. J Exp Med. (2000) 192:1553–62. doi: 10.1084/jem.192.11.1553
23. Fahey LM, Wilson EB, Elsaesser H, Fistonich CD, McGavern DB, Brooks DG. Viral persistence redirects CD4 T cell differentiation toward T follicular helper cells. J Exp Med. (2011) 208:987–99. doi: 10.1084/jem.20101773
24. Harker JA, Lewis GM, Mack L, Zuniga EI. Late interleukin-6 escalates T follicular helper cell responses and controls a chronic viral infection. Science. (2011) 334:825–9. doi: 10.1126/science.1208421
25. Boyden AW, Legge KL, Waldschmidt TJ. Pulmonary infection with influenza A virus induces site-specific germinal center and T follicular helper cell responses. PLoS ONE. (2012) 7:e40733. doi: 10.1371/journal.pone.0040733
26. Feng J, Hu X, Guo H, Sun X, Wang J, Xu L, et al. Patients with chronic hepatitis C express a high percentage of CD4(+)CXCR5(+) T follicular helper cells. J Gastroenterol. (2012) 47:1048–56. doi: 10.1007/s00535–012-0568–1
27. Hale JS, Youngblood B, Latner DR, Mohammed AU, Ye L, Akondy RS, et al. Distinct memory CD4+ T cells with commitment to T follicular helper- and T helper 1-cell lineages are generated after acute viral infection. Immunity. (2013) 38:805–17. doi: 10.1016/j.immuni.2013.02.020
28. Glatman Zaretsky A, Taylor JJ, King IL, Marshall FA, Mohrs M, Pearce EJ. T follicular helper cells differentiate from Th2 cells in response to helminth antigens. J Exp Med. (2009) 206:991–9. doi: 10.1084/jem.20090303
29. King IL, Mohrs M. IL-4-producing CD4+ T cells in reactive lymph nodes during helminth infection are T follicular helper cells. J Exp Med. (2009) 206:1001–7. doi: 10.1084/jem.20090313
30. Kemeny DM. The role of the T follicular helper cells in allergic disease. Cell Mol Immunol. (2012) 9:386–9. doi: 10.1038/cmi.2012.31
31. Kamekura R, Shigehara K, Miyajima S, Jitsukawa S, Kawata K, Yamashita K, et al. Alteration of circulating type 2 follicular helper T cells and regulatory B cells underlies the comorbid association of allergic rhinitis with bronchial asthma. Clin Immunol. (2015) 158:204–11. doi: 10.1016/j.clim.2015.02.016
32. Meli AP, Fontes G, Leung Soo C, King IL. T Follicular helper cell-derived IL-4 is required for IgE production during intestinal helminth infection. J Immunol. (2017) 199:244–52. doi: 10.4049/jimmunol.1700141
33. Ding R, Gao W, He Z, Wang H, Conrad CD, Dinney CM, et al. Overrepresentation of Th1- and Th17-like follicular helper T cells in coronary artery disease. J Cardiovasc Transl Res. (2015) 8:503–5. doi: 10.1007/s12265–015-9662–0
34. Smith C, Buhlmann JE, Wang X, Bartlett A, Lim B, Barrington RA. CD275-Independent IL-17-producing T follicular helper-like cells in lymphopenic autoimmune-prone mice. J Immunol. (2016) 196:4935–46. doi: 10.4049/jimmunol.1402193
35. Yang Y, Zhang M, Ye Y, Ma S, Fan L, Li Z. High frequencies of circulating Tfh-Th17 cells in myasthenia gravis patients. Neurol Sci. (2017) 38:1599–1608. doi: 10.1007/s10072–017-3009–3
36. Moser B, Ebert L. Lymphocyte traffic control by chemokines: follicular B helper T cells. Immunol Lett. (2003) 85:105–12. doi: 10.1016/S0165–2478(02)00233-X
37. Hardtke S, Ohl L, Forster R. Balanced expression of CXCR5 and CCR7 on follicular T helper cells determines their transient positioning to lymph node follicles and is essential for efficient B-cell help. Blood. (2005) 106:1924–31. doi: 10.1182/blood-2004–11-4494
38. Rasheed AU, Rahn HP, Sallusto F, Lipp M, Muller G. Follicular B helper T cell activity is confined to CXCR5(hi)ICOS(hi) CD4 T cells and is independent of CD57 expression. Eur J Immunol. (2006) 36:1892–903. doi: 10.1002/eji.200636136
39. Fazilleau N, Mark L, McHeyzer-Williams LJ, McHeyzer-Williams MG. Follicular helper T cells: lineage and location. Immunity. (2009) 30:324–35. doi: 10.1016/j.immuni.2009.03.003
40. Moser B. CXCR5, the defining marker for follicular B helper T (TFH) Cells. Front Immunol. (2015) 6:296. doi: 10.3389/fimmu.2015.00296
41. Vinuesa CG, Tangye SG, Moser B, Mackay CR. Follicular B helper T cells in antibody responses and autoimmunity. Nat Rev Immunol. (2005) 5:853–65. doi: 10.1038/nri1714
42. Bryant VL, Ma CS, Avery DT, Li Y, Good KL, Corcoran LM, et al. Cytokine-mediated regulation of human B cell differentiation into Ig-secreting cells: predominant role of IL-21 produced by CXCR5+ T follicular helper cells. J Immunol. (2007) 179:8180–90. doi: 10.4049/jimmunol.179.12.8180
43. King C, Tangye SG, Mackay CR. T follicular helper (TFH) cells in normal and dysregulated immune responses. Annu Rev Immunol. (2008) 26:741–66. doi: 10.1146/annurev.immunol.26.021607.090344
44. Nurieva RI, Chung Y, Hwang D, Yang XO, Kang HS, Ma L, et al. Generation of T follicular helper cells is mediated by interleukin-21 but independent of T helper 1, 2, or 17 cell lineages. Immunity. (2008) 29:138–49. doi: 10.1016/j.immuni.2008.05.009
45. Vogelzang A, McGuire HM, Yu D, Sprent J, Mackay CR, King C. A fundamental role for interleukin-21 in the generation of T follicular helper cells. Immunity. (2008) 29:127–37. doi: 10.1016/j.immuni.2008.06.001
46. Haynes NM. Follicular associated T cells and their B-cell helper qualities. Tissue Antigens. (2008) 71:97–104. doi: 10.1111/j.1399–0039.2007.00995.x
47. Johnston RJ, Poholek AC, DiToro D, Yusuf I, Eto D, Barnett B, et al. Bcl6 and Blimp-1 are reciprocal and antagonistic regulators of T follicular helper cell differentiation. Science. (2009) 325:1006–10. doi: 10.1126/science.1175870
48. Kroenke MA, Eto D, Locci M, Cho M, Davidson T, Haddad EK, et al. Bcl6 and Maf cooperate to instruct human follicular helper CD4 T cell differentiation. J Immunol. (2012) 188:3734–44. doi: 10.4049/jimmunol.1103246
49. Liu X, Chen X, Zhong B, Wang A, Wang X, Chu F, et al. Transcription factor achaete-scute homologue 2 initiates follicular T-helper-cell development. Nature. (2014) 507:513–8. doi: 10.1038/nature12910
50. Andris F, Denanglaire S, Anciaux M, Hercor M, Hussein H, Leo O. The Transcription Factor c-Maf promotes the differentiation of follicular helper T cells. Front Immunol. (2017) 8:480. doi: 10.3389/fimmu.2017.00480
51. Xu L, Cao Y, Xie Z, Huang Q, Bai Q, Yang X, et al. The transcription factor TCF-1 initiates the differentiation of T(FH) cells during acute viral infection. Nat Immunol. (2015) 16:991–9. doi: 10.1038/ni.3229
52. Nurieva RI, Chung Y, Martinez GJ, Yang XO, Tanaka S, Matskevitch TD, et al. Bcl6 mediates the development of T follicular helper cells. Science. (2009) 325:1001–5. doi: 10.1126/science.1176676
53. Yu D, Rao S, Tsai LM, Lee SK, He Y, Sutcliffe EL, et al. The transcriptional repressor Bcl-6 directs T follicular helper cell lineage commitment. Immunity. (2009) 31:457–68. doi: 10.1016/j.immuni.2009.07.002
54. Baumjohann D, Okada T, Ansel KM. Cutting Edge: distinct waves of BCL6 expression during T follicular helper cell development. J Immunol. (2011) 187:2089–92. doi: 10.4049/jimmunol.1101393
55. Ryg-Cornejo V, Ioannidis LJ, Ly A, Chiu CY, Tellier J, Hill DL, et al. Severe malaria infections impair germinal center responses by inhibiting t follicular helper cell differentiation. Cell Rep. (2016) 14:68–81. doi: 10.1016/j.celrep.2015.12.006
56. Fang D, Cui K, Mao K, Hu G, Li R, Zheng M, et al. Transient T-bet expression functionally specifies a distinct T follicular helper subset. J Exp Med. (2018) 215:2705–15. doi: 10.1084/jem.20180927
57. Ma X, Nakayamada S, Kubo S, Sakata K, Yamagata K, Miyazaki Y, et al. Expansion of T follicular helper-T helper 1 like cells through epigenetic regulation by signal transducer and activator of transcription factors. Ann Rheum Dis. (2018) 77:1354–61. doi: 10.1136/annrheumdis-2017–212652
58. Weinstein JS, Laidlaw BJ, Lu Y, Wang JK, Schulz VP, Li N, et al. STAT4 and T-bet control follicular helper T cell development in viral infections. J Exp Med. (2018) 215:337–55. doi: 10.1084/jem.20170457
59. Hercor M, Anciaux M, Denanglaire S, Debuisson D, Leo O, Andris F. Antigen-presenting cell-derived IL-6 restricts the expression of GATA3 and IL-4 by follicular helper T cells. J Leukoc Biol. (2017) 101:5–14. doi: 10.1189/jlb.1HI1115–511R
60. Schmitt N, Liu Y, Bentebibel SE, Munagala I, Bourdery L, Venuprasad K, et al. The cytokine TGF-beta co-opts signaling via STAT3-STAT4 to promote the differentiation of human TFH cells. Nat Immunol. (2014) 15:856–65. doi: 10.1038/ni.2947
61. Reinhardt RL, Liang HE, Locksley RM. Cytokine-secreting follicular T cells shape the antibody repertoire. Nat Immunol. (2009) 10:385–93. doi: 10.1038/ni.1715
62. Lee SK, Silva DG, Martin JL, Pratama A, Hu X, Chang PP, et al. Interferon-gamma excess leads to pathogenic accumulation of follicular helper T cells and germinal centers. Immunity. (2012) 37:880–92. doi: 10.1016/j.immuni.2012.10.010
63. Luthje K, Kallies A, Shimohakamada Y, Belz GT, Light A, Tarlinton DM, et al. The development and fate of follicular helper T cells defined by an IL-21 reporter mouse. Nat Immunol. (2012) 13:491–8. doi: 10.1038/ni.2261
64. Wong MT, Chen J, Narayanan S, Lin W, Anicete R, Kiaang HT, et al. Mapping the Diversity of follicular helper T cells in human blood and tonsils using high-dimensional mass cytometry analysis. Cell Rep. (2015) 11:1822–33. doi: 10.1016/j.celrep.2015.05.022
65. Vijayanand P, Seumois G, Simpson LJ, Abdul-Wajid S, Baumjohann D, Panduro M, et al. Interleukin-4 production by follicular helper T cells requires the conserved Il4 enhancer hypersensitivity site V. Immunity. (2012) 36:175–87. doi: 10.1016/j.immuni.2011.12.014
66. Peng SL, Szabo SJ, Glimcher LH. T-bet regulates IgG class switching and pathogenic autoantibody production. Proc Natl Acad Sci USA. (2002) 99:5545–50. doi: 10.1073/pnas.082114899
67. Szabo SJ, Kim ST, Costa GL, Zhang X, Fathman CG, Glimcher LH. A novel transcription factor, T-bet, directs Th1 lineage commitment. Cell. (2000) 100:655–69. doi: 10.1016/S0092–8674(00)80702–3
68. Hwang ES, Szabo SJ, Schwartzberg PL, Glimcher LH. T helper cell fate specified by kinase-mediated interaction of T-bet with GATA-3. Science. (2005) 307:430–3. doi: 10.1126/science.1103336
69. Oestreich KJ, Weinmann AS. T-bet employs diverse regulatory mechanisms to repress transcription. Trends Immunol. (2012) 33:78–83. doi: 10.1016/j.it.2011.10.005
70. Kallies A, Good-Jacobson KL. Transcription factor T-bet orchestrates lineage development and function in the immune system. Trends Immunol. (2017) 38:287–97. doi: 10.1016/j.it.2017.02.003
71. Lazarevic V, Glimcher LH. T-bet in disease. Nat Immunol. (2011) 12:597–606. doi: 10.1038/ni.2059
72. Lazarevic V, Glimcher LH, Lord GM. T-bet: a bridge between innate and adaptive immunity. Nat Rev Immunol. (2013) 13:777–89. doi: 10.1038/nri3536
73. Gordon SM, Chaix J, Rupp LJ, Wu J, Madera S, Sun JC, et al. The transcription factors T-bet and Eomes control key checkpoints of natural killer cell maturation. Immunity. (2012) 36:55–67. doi: 10.1016/j.immuni.2011.11.016
74. Townsend MJ, Weinmann AS, Matsuda JL, Salomon R, Farnham PJ, Biron CA, et al. T-bet regulates the terminal maturation and homeostasis of NK and Valpha14i NKT cells. Immunity. (2004) 20:477–94. doi: 10.1016/S1074–7613(04)00076–7
75. Werneck MB, Lugo-Villarino G, Hwang ES, Cantor H, Glimcher LH. T-bet plays a key role in NK-mediated control of melanoma metastatic disease. J Immunol. (2008) 180:8004–10. doi: 10.4049/jimmunol.180.12.8004
76. Klose CS, Kiss EA, Schwierzeck V, Ebert K, Hoyler T, d'Hargues Y, et al. A T-bet gradient controls the fate and function of CCR6-RORgammat+ innate lymphoid cells. Nature. (2013) 494:261–5. doi: 10.1038/nature11813
77. Chen L, He W, Kim ST, Tao J, Gao Y, Chi H, et al. Epigenetic and transcriptional programs lead to default IFN-gamma production by gammadelta T cells. J Immunol. (2007) 178:2730–6. doi: 10.4049/jimmunol.178.5.2730
78. Lugo-Villarino G, Maldonado-Lopez R, Possemato R, Penaranda C, Glimcher LH. T-bet is required for optimal production of IFN-gamma and antigen-specific T cell activation by dendritic cells. Proc Natl Acad Sci USA. (2003) 100:7749–54. doi: 10.1073/pnas.1332767100
79. Matsuda JL, Zhang Q, Ndonye R, Richardson SK, Howell AR, Gapin L. T-bet concomitantly controls migration, survival, and effector functions during the development of Valpha14i NKT cells. Blood. (2006) 107:2797–805. doi: 10.1182/blood-2005–08-3103
80. Cruz-Guilloty F, Pipkin ME, Djuretic IM, Levanon D, Lotem J, Lichtenheld MG, et al. Runx3 and T-box proteins cooperate to establish the transcriptional program of effector CTLs. J Exp Med. (2009) 206:51–9. doi: 10.1084/jem.20081242
81. Takemoto N, Intlekofer AM, Northrup JT, Wherry EJ, Reiner SL. Cutting Edge: IL-12 inversely regulates T-bet and eomesodermin expression during pathogen-induced CD8+ T cell differentiation. J Immunol. (2006) 177:7515–9. doi: 10.4049/jimmunol.177.11.7515
82. Wang NS, McHeyzer-Williams LJ, Okitsu SL, Burris TP, Reiner SL, McHeyzer-Williams MG. Divergent transcriptional programming of class-specific B cell memory by T-bet and RORalpha. Nat Immunol. (2012) 13:604–11. doi: 10.1038/ni.2294
83. Marshall HD, Chandele A, Jung YW, Meng H, Poholek AC, Parish IA, et al. Differential expression of Ly6C and T-bet distinguish effector and memory Th1 CD4(+) cell properties during viral infection. Immunity. (2011) 35:633–46. doi: 10.1016/j.immuni.2011.08.016
84. Nakayamada S, Kanno Y, Takahashi H, Jankovic D, Lu KT, Johnson TA, et al. Early Th1 cell differentiation is marked by a Tfh cell-like transition. Immunity. (2011) 35:919–31. doi: 10.1016/j.immuni.2011.11.012
85. Oestreich KJ, Huang AC, Weinmann AS. The lineage-defining factors T-bet and Bcl-6 collaborate to regulate Th1 gene expression patterns. J Exp Med. (2011) 208:1001–13. doi: 10.1084/jem.20102144
86. Oestreich KJ, Mohn SE, Weinmann AS. Molecular mechanisms that control the expression and activity of Bcl-6 in TH1 cells to regulate flexibility with a TFH-like gene profile. Nat Immunol. (2012) 13:405–11. doi: 10.1038/ni.2242
87. Rasheed MA, Latner DR, Aubert RD, Gourley T, Spolski R, Davis CW, et al. Interleukin-21 is a critical cytokine for the generation of virus-specific long-lived plasma cells. J Virol. (2013) 87:7737–46. doi: 10.1128/JVI.00063–13
88. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. (2005) 102:15545–50. doi: 10.1073/pnas.0506580102
89. Wu J, Mao X, Cai T, Luo J, Wei L. KOBAS server: a web-based platform for automated annotation and pathway identification. Nucleic Acids Res. (2006) 34:W720–4. doi: 10.1093/nar/gkl167
90. Oxenius A, Bachmann MF, Zinkernagel RM, Hengartner H. Virus-specific MHC-class II-restricted TCR-transgenic mice: effects on humoral and cellular immune responses after viral infection. Eur J Immunol. (1998) 28:390–400. doi: 10.1002/(SICI)1521–4141(199801)28:01<390::AID-IMMU390>3.0.CO;2-O
91. Yusuf I, Kageyama R, Monticelli L, Johnston RJ, Ditoro D, Hansen K, et al. Germinal center T follicular helper cell IL-4 production is dependent on signaling lymphocytic activation molecule receptor (CD150). J Immunol. (2010) 185:190–202. doi: 10.4049/jimmunol.0903505
92. Xiao L, Jia L, Zhang Y, Yu S, Wu X, Yang B, et al. Human IL-21+IFN-gamma+CD4+ T cells in nasal polyps are regulated by IL-12. Sci Rep. (2015) 5:12781. doi: 10.1038/srep12781
93. Ise W, Kohyama M, Schraml BU, Zhang T, Schwer B, Basu U, et al. The transcription factor BATF controls the global regulators of class-switch recombination in both B cells and T cells. Nat Immunol. (2011) 12:536–43. doi: 10.1038/ni.2037
94. Stone EL, Pepper M, Katayama CD, Kerdiles YM, Lai CY, Emslie E, et al. ICOS coreceptor signaling inactivates the transcription factor FOXO1 to promote Tfh cell differentiation. Immunity. (2015) 42:239–51. doi: 10.1016/j.immuni.2015.01.017
95. Reed JM, Branigan PJ, Bamezai A. Interferon gamma enhances clonal expansion and survival of CD4+ T cells. J Interferon Cytokine Res. (2008) 28:611–22. doi: 10.1089/jir.2007.0145
96. Velu V, Mylvaganam GH, Gangadhara S, Hong JJ, Iyer SS, Gumber S, et al. Induction of Th1-Biased T follicular helper (Tfh) cells in lymphoid tissues during chronic simian immunodeficiency virus infection defines functionally distinct germinal center Tfh cells. J Immunol. (2016) 197:1832–42. doi: 10.4049/jimmunol.1600143
97. Obeng-Adjei N, Portugal S, Tran TM, Yazew TB, Skinner J, Li S, et al. Circulating Th1-Cell-type Tfh cells that exhibit impaired B cell help are preferentially activated during acute malaria in children. Cell Rep. (2015) 13:425–39. doi: 10.1016/j.celrep.2015.09.004
Keywords: T-bet, follicular helper T cells, type I immune response, humoral response, T cell differentiation, transcriptional regulation
Citation: Wang P, Wang Y, Xie L, Xiao M, Wu J, Xu L, Bai Q, Hao Y, Huang Q, Chen X, He R, Li B, Yang S, Chen Y, Wu Y and Ye L (2019) The Transcription Factor T-Bet Is Required for Optimal Type I Follicular Helper T Cell Maintenance During Acute Viral Infection. Front. Immunol. 10:606. doi: 10.3389/fimmu.2019.00606
Received: 28 November 2018; Accepted: 07 March 2019;
Published: 29 March 2019.
Edited by:
Georgia Fousteri, San Raffaele Hospital (IRCCS), ItalyReviewed by:
Masato Kubo, Tokyo University of Science, JapanColleen Jean Winstead, Merck Sharp & Dohme Corp, United States
Copyright © 2019 Wang, Wang, Xie, Xiao, Wu, Xu, Bai, Hao, Huang, Chen, He, Li, Yang, Chen, Wu and Ye. 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: Yaokai Chen, yaokaichen@hotmail.com
Yuzhang Wu, wuyuzhang@tmmu.edu.cn
Lilin Ye, yelilinlcmv@tmmu.edu.cn
†These authors have contributed equally to this work