- Department of Dermatology and Venerology, The First Hospital of Jilin University, Changchun, China
Sporotrichosis is a subcutaneous mycotic infection, and Sporothrix globosa is one of the causative agents with a worldwide distribution, notably in Asia. However, the immune profile in human sporotrichosis caused by S. globosa still remains obscure. Here, we demonstrated enhanced Th2 response in circulation with significant increases in Th2 frequency, Th2/Tregs as well as IL-4 seretion in patients. Elevated IL-17A+Th17 percentage was accompanied with reduced IL-17A level in serum, which may imply a dysfunction of this CD4+T subset in S. globosa infection. In addition, Th2 percentage, the ratios of Th2/Tregs and Th17/Tregs were all raised in patients with fixed cutaneous form, while only Th2/Tregs displayed increment in lymphocutaneous form. Meanwhile, the percentage of double negative B cells was significantly increased and positively correlated with Th2 and Tregs in whole patients. Except naïve B cells, all memory B cells together with Th2 cells increased in patients with short duration (less than 6 months), which may suggest a collaboration of T cells with altered B cell profile in human sporotrichosis caused by S. globosa. In consistent with the changes of IFN-γ+Th1, IL-4+Th2 and IL-17A+Th17 in patients with short duration, the percentages of these effector T cells all expanded when cocultured with S. globosa yeast cells in vitro. These data shed light on the potential involvement of peripheral T and B cell immunity against this mycotic infection and indicated that different immune responses existed in different stages of sporotrichosis; meanwhile different immune profile may contribute to different clinical manifestations of this disease.
Introduction
Sporotrichosis is a subcutaneous mycotic infection caused by Sporothrix complex, among which Sporothrix schenckii sensu stricto (S. schenckii), Sporothrix brasiliensis (S. brasiliensis), and Sporothrix globosa (S. globosa) are most clinically relevant. S. globosa has a global distribution evidenced by isolates from Europe (the United Kingdom, Spain, and Italy), the United States, South America (Mexico, Guatemala, and Columbia), and Asia (India, China, and Japan) (1–6). Compared with other species of Sporothrix complex, the clinical manifestations caused by S. globosa are relatively mild owing to its low to moderate virulence (1, 7–9), which indicates it may trigger different immune responses in the host.
Up to now, the defense mechanisms against Sporothrix complex in human are very limited. In mice, previous studies have acquired a consensus that T cell-mediated immunity particularly Th1 response is considered critical to host defense against S. schenckii infection (10, 11) which participates in the granuloma formation in skin lesions (12) and enhances phagocytosis of the pathogen by macrophages (13, 14). Besides Th1, evidence suggested that control of S. schenckii is Th17 dependent (12, 15). Th17 response is found crucial for optimal S. schenckii clearance (16) and the host tends to use a polarized Th17 response when phagocytosis capacity is impaired in vivo (17). So far, only one published literature mentioned that great reduction in Th1 and Th17 was found in patients with sporotrichosis caused by S. globosa (18).
In addition to cellular immunity, humoral response in experimental sporotrichosis in mice has also been evaluated (19). IgG antibodies against Sporothrix antigens could be detected throughout the infection in mice (20). Anti-Gp70 is one of the antibodies that has been widely recognized as playing protective role against S. schenckii and S. brasiliensis infection (21, 22). Recently, anti-enolase antibodies were also characterized involved in antifungal responses in S. schenckii and S. brasilensis infected mice (23, 24). In human sporotrichosis, the IgG, IgM and IgA antibodies were detected in all of the clinical forms (25).
Despite the progress in the knowledge of humoral response in mice, little is known about the humoral immunity against clinically relevant Sporothrix species, especially S. globosa in human (1, 19). Follicular helper T cells (Tfh), a newly discovered CD4+ T cell subgroup, bridges the cellular and humoral immunity by assisting B cells’ differentiation and antibody production (26). However, in the host defense against Sporothrix complex, the change of B cells and Tfh cells and the mechanism driving the production of antibodies are still not clear.
To elucidate these questions, we investigated the immune profile in circulation by application of phenotyping circulating cells and cytokine detection in sera as well as cells-pathogen coculture in vitro. We showed a Th2 predominant response together with altered B cell profile in whole patients. Moreover, double negative B cells (IgD-CD27- B cells) were positively correlated with Th2 and Tregs which indicated that B and T cells collaborated in the immunity against this pathogen. An elevated ratio of Th17/Tregs was noticed in patients with fixed form, not in patients with lymphocutaneous type. Except the decrease in IgG2, we did not find any change in the level of IgE, total IgG and IgM as well as the proportion of Tfh cells in patients, though CD27+ memory B subsets responsible for antibodies production were greatly changed.
Materials and Methods
Subjects
Sixty-eight patients with sporotrichosis (SP) and 45 age and sex matched healthy controls (HC) were recruited in this study. The definitive diagnosis of sporotrichosis was established based on clinical manifestation, skin pathology, and positive fungal culture. Isolates from patients were identified as S. globosa based on morphological characteristics and calmodulin-encoding gene (CAL) sequencing. All patients had not been treated with any drugs and had no history of any auto-immune or infectious diseases by asking questions about their symptoms and medical history. Patients enrolled were divided into subgroups according to their duration (till the moment of diagnosis) and clinical types respectively: shorter duration (SD) (< 6 months) and longer duration (LD) (> 6 months); fixed cutaneous form (FF) and lymphocutaneous form (LF). The detail information was listed in Tables 1, 2, and S1. Written informed consents from all subjects were obtained. The study was approved by the First Hospital of Jilin University Research Ethics Committee.
Multiparametric Flow Cytometry Analysis
Peripheral blood mononuclear cells (PBMCs) were freshly isolated from heparinized venous blood by density gradient centrifugation on Ficoll-Lympholyte (Cedarlane Laboratories Limited, Ontario, Canada). Cells were stained with antibodies specific to surface markers at 4°C for 30 min. For intracellular staining, cells were further fixed, permeabilized and incubated with antibodies according to the manufacturer’s protocol. For intracellular detection of IFN-γ, IL-4, and IL-17A, freshly isolated PBMCs were first stimulated with Protein Transport Inhibitor (BD GolgiPlug containing Brefeldin A, BD Biosciences, USA) for 4 h in vitro. Unlabeled cells were used as autofluorescence control and IgG isotypes of corresponding antibodies were used as non-specific fluorescence control. FACS analysis was performed with a FACS Calibur flow cytometer (BD Biosciences, US). A total of 3×104-5×104 events for a selected gate (PBMCs gate) were collected for each sample and the data were analyzed by the FlowJo software (version 7.6). All antibodies used were listed in Table 3.
ELISA
Serum levels of total IgG, IgG subclasses (IgG1, IgG2, IgG3 and IgG4), IgM, IgE, and T cell cytokines (IFN-γ, IL-4, IL-17A, and TGF-β1), respectively were determined using Invitrogen human uncoated ELISA Kit (Bender MedSystems GmbH, Campus Vienna BIocenter 2,1030 Vienna, Austria). For Ig measurement, we prediluted the samples in Assay Buffer according to the manufacturer’s instructions. The values of OD were detected by spectrophotometric instrument at a wave length of 450 nm.
Strain and Co-Culture of PBMCs With S. globosa
The S. globosa strain was Sporothrix globosa ATCC4912, and yeast cells were grown in Brain Heart Infusion Agar (BD Bioxon) for 5–7 days at 37°C. After centrifugation at 3,000 r/min for 10 min, the yeasts were washed twice with phosphate-buffered saline (PBS), and then were resuspended in complete RPMI 1640 medium (supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin). 2×105 freshly isolated PBMCs were seeded in one well of 96-well plate and challenged with S. globosa yeast cells at different ratios (1:1; 1:2.5; 1:5; and 1:10). The coculture at each ratio had triplicates. 3 days later, PBMCs were restimulated with Protein Transport Inhibitor for additional 4 h. After restimulation, cells were harvested and stained with surface and intracellular mAbs as mentioned above.
Statistical Analysis
All statistical analysis was performed using GraphPad Prism 5.01 (GraphPad software Inc, San Diego, CA, USA). The difference among three groups was analyzed with one-way ANOVA with Bonferroni correction. Unpaired t-test (Welch corrected) was utilized to analyze the data between two groups. The data was presented as mean±SD and P value <0.05 was considered significant.
Results
Th2 and Th17 Inclined Responses in Peripheral Blood of Sporotrichosis Patients Infected With S. globosa
Previous studies have suggested that T cell-mediated immunity plays an essential role in the immunity against infection with S. schenckii and S. brasiliensis (19). To make clear T cells’ profile in S. globosa infection in human, effector CD4+ T (Teff) cells were evaluated in SP (patients with sporotrichosis) and HC (Figures 1, S1). We analyzed the expressions of lineage specific transcription factors (T-bet, GATA-3, and ROR-γt) and cytokines (IFN-γ, IL-4, and IL-17A) in CD4+ T cells for the detection of Th1, Th2, and Th17 cells, respectively (Figures 1A, S1A). Interestingly, the percentage of CD4+GATA-3+ Th2 cells was significantly higher in whole patients than that of HC (P = 0.01, n = 49) (Figure 1B) which was further confirmed by the greatly increased frequency of CD4+IL-4+ Th2 cells (P < 0.0001, n = 18) (Figure S1B). Moreover, the proportion of Th2 cells remained at higher level throughout the course of disease as evidenced by its great increase both in the patients with short duration (SD) and long duration (LD) (P = 0.04, n = 24 and P = 0.02, n = 25 for GATA3+Th2; P < 0.0001, n = 7 and P = 0.0018, n = 11 for IL-4+Th2) (Figures 1C, S1C). In addition, similar increment was also observed in patients with fixed cutaneous form (FF) (P = 0.005, n = 16) rather than with lymphocutaneous form (LF) (Figure 1D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/570888/fimmu-11-570888-HTML/image_m/fimmu-11-570888-g001.jpg)
Figure 1 Altered distribution of Th1, Th2, and Th17 in PBMCs of patients. (A) PBMCs were stained intracellularly with T-bet, GATA-3, and ROR-γt mAbs after surface staining of CD4 mAb. According to CD4 staining and SSC, CD4+ T cells were gated. The parameters shown in quadrants of the representative graphs are mean frequency of GATA-3 for each group. (B–D) The average percentages of CD4+T-bet+ Th1, CD4+GATA-3+ Th2, CD4+ROR-γt+ Th17 are compared between HC (n = 24) and whole patients (n = 50) as well as subgroups of patients (SD, n = 24; LD, n = 26; FF, n = 33; LF, n = 17). Error bars represent mean±SD. **P < 0.01, *P < 0.05, and NS P ≥ 0.05.
It had been reported that Th1 response dominated in the immunity against infection with S. schenckii and S. brasiliensis (11). However, neither assessed by CD4+T-bet+ nor CD4+IFN-γ+ as definition of Th1 cells did we find statistical change in the proportion of this subset in whole patients compared with HC (Figures 1B, S1B). Meanwhile, we observed the dynamic change of Th1 cells in patients with different duration and found that only CD4+IFN-γ+ Th1 frequency had a significant elevation in SD (n = 7, P = 0.0006) while it dropped down in LD compared with HC (Figures 1C, S1C). No statistical difference was found among HC and patients with FF and LF (Figure 1D).
Besides, evidence suggested that control of S. schenckii is Th1/Th17 dependent (27) and optimal fungal clearance depends on an intact Th17 response during the S. schenckii systemic infection in mice (16). Similarly, in our study, we found that the frequency of CD4+IL-17A+ Th17 cells was greatly elevated in whole patients (P < 0.0001, n = 18) and remained high level in both patients with SD (P < 0.0001, n = 7) and LD (P < 0.0001, n = 11) (Figures S1B, C). However, CD4+RORγt+ Th17 cells had no change in whole patients and in patients with different duration or clinical types (Figures 1B–D).
Imbalance Between Effector T Cells and Immunosuppressive T Cells in PBMCs of Sporotrichosis Patients
Tregs are potent immunosuppressive CD4+ T cells which keep the balance of immunity. The ratios of Teff cells (Th1, Th2, Th17)/Tregs can better reflect the inclination of T cell-mediated immunity. We next analyzed CD4+CD25+CD127- Tregs in the subjects’ peripheral blood (Figure 2A). As shown in Figure 2B, Tregs frequency showed no statistical change in whole patients and patients with different clinical types. However, a significantly decreased percentage of Tregs was found in SD group (P = 0.03, n = 20) compared with HC; while in LD, it showed a trend towards normalization.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/570888/fimmu-11-570888-HTML/image_m/fimmu-11-570888-g002.jpg)
Figure 2 Tregs’ frequency and the ratio of effector T cells/Tregs in PBMCs of patients and HC. (A) PBMCs were stained with CD4, CD25 and CD127 mAbs. CD4+ T cells were gated according to CD4 staining and SSC and the proportion of CD4+CD25+CD127- Tregs in CD4+ T cells were determined. The graphs are representative for patients, indicated subgroups and HC. Mean value of Tregs’ percentage is shown in each target quadrant. (B–E) Statistical graphs for comparing Tregs’ percentage and Th/Tregs ratios of the whole patients and indicated subgroups with HC. (SP, n = 46; SD, n = 20; LD, n = 26, FF, n = 30; LF, n = 16; HC, n = 25). Error bars represent mean±SD. **P < 0.01, *P < 0.05, and NS P ≥ 0.05.
Remarkably, Th2 inclined immune response was evidenced again by significant elevation in the ratio of Th2/Tregs not only in whole patients (P = 0.001, n = 46) but also in patients with different duration and clinical types (SD, P = 0.005, n = 20; LD, P = 0.006, n = 26; LF, P = 0.03, n = 16; and FF, P = 0.001, n = 30) (Figures 2C–E). In contrast to Th2 response, the ratio of Th1/Tregs kept statistically comparable to HC in patients regardless of the difference in duration and clinical types (Figures 2C–E). Besides, higher ratio of Th17/Tregs was found in whole patients (P = 0.04, n = 46) and in FF subgroup (P = 0.01, n = 30) (Figures 2C, E), which may suggest the participation of Th17 reaction in the formation of FF.
Altered Distribution of CD19+ B Cell Subsets in Circulation of Sporotrichosis Patients
Previous reports indicated that humoral immunity also participate in the antifungal response in S. schenckii and S. brasiliensis infection (28). To figure out whether B cell populations are involved in the defense against S. globosa infection, we determined the profile of circulating CD19+ B cells in patients and HC. According to the expression of CD27 and IgD, B cells were divided into four subsets (Figure 3A): IgD+CD27- naïve B cells (NB cells), IgD+CD27+ unswitched memory B cells (USM B cells), IgD-CD27+ switched memory B cells (SM B cells), and IgD-CD27- double negative B cells (DN B cells). Although the frequency of total CD19+ B cells was unchanged in patients (Figure 3B), remarkably altered distribution of B cell subsets in whole patients (n = 48) was found (Figure 3C): a great reduction of USM B cells’ percentage (P = 0.03), significant increase of SM B cells (P = 0.03) and DN B cells (P < 0.001). Meanwhile, we noticed a similar variation in patients with SD (P = 0.03, P = 0.007, P = 0.008, respectively) (Figure 3D). However, in patients with LD, only DN B subset presented a significant elevation (P = 0.001). Patients with FF and LF both demonstrated enhancement in DN B cells (P = 0.007, P = 0.001, respectively) (Figure 3E). The frequency of naïve B cells remained comparable to that of HC in whole patients and patients with different duration and clinical types.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/570888/fimmu-11-570888-HTML/image_m/fimmu-11-570888-g003.jpg)
Figure 3 Phenotypic characterization of circulating B cell subsets in patients and HC. (A) Gated by CD19+ B cells, IgD+CD27- NB, IgD+CD27+ USM B cells, IgD-CD27+ SM B cells, and IgD-CD27- DN B cells were identified. The graphs are representative for HC and patients with different duration. Mean value of each B cell subset’s percentage is shown in the quadrants. (B) Statistical graphs for comparison of CD19+ B cells’ percentages between patient groups and HC. (C–E) Statistical graphs for distinct CD19+ B cells subsets between patients (SP, n = 48; SD, n = 23; LD, n = 25; FF, n = 33; LF, n = 15.) and HC (n = 25). Error bars represent mean±SD. *P < 0.05, **P < 0.01, ***P < 0.001, and NS P ≥ 0.05.
Unchanged Tfh Frequency and Decreased IgG2 in Sporotrichosis Patients
It was reported that Tfh cells play a critical role in helping B cells differentiation and antibody production (29). To clarify whether the alteration of B cell subsets was associated with Tfh cells, we next investigated Tfh cells in circulation (cTfh). Compared with HC, the percentage of cTfh cells showed no statistical change in patients (Figure 4A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/570888/fimmu-11-570888-HTML/image_m/fimmu-11-570888-g004.jpg)
Figure 4 Unchanged cTfh cells and altered Ig profile in sporotrichosis patients. (A) Statistical graphs in the upper row were for CD4+CXCR5+Tfh while the lower row were for CD4+CXCR5+PD1+Tfh. The comparison of cTfh percentages were between HC (n = 24) and whole patients (n = 50), patients with different duration (< 6 mon, n = 24; > 6 mon, n = 26) and clinical types (FF, n = 33; LF, n = 17). (B) Comparison of serum levels of total IgG, IgM, and IgE between patients (n = 46) and HC (n = 24). (C, D) Distribution of IgG subtypes (IgG1, IgG2, and IgG3 and IgG4) in patients (in whole: n = 20; SD: n = 7; LD: n = 13) and HC (n = 19). Error bars represent mean±SD. *P < 0.05, **P < 0.01, and NS P ≥ 0.05.
Early studies have revealed that USM B and SM B cells are mainly responsible for IgM and IgG secretion respectively (30) and Th2 cells are widely recognized as IgE inducer. Considering significant changes in B cells and Th2 cells in our study, we next measured the concentration of IgG, IgM, and IgE in sera of the subjects. However, compared with HC, the levels of total IgG, IgM, and IgE had no significant changes in whole patients (Figure 4B) and patients with SD, LD, FF, and LF (data not shown). We further assessed subclasses of IgG to explore differences in small subpopulations. In comparison with HC, a decreased level of IgG2 (P = 0.02, n = 20) was found in whole patients. Both of the IgG2 and IgG4 were found significantly decreased in SD patients (P = 0.002, P = 0.04, respectively, n = 7) (Figures 4C, D).
Correlation of B Cells With CD4+ T Cells in Sporotrichosis Patients
Researchers discovered that in addition to dendritic cells, B cells also take part in the differentiation of CD4+ T cells (30). To find out if there is any association between altered B cell subsets and CD4+ T cells, we further performed correlation analysis of all the altered B cell subsets (i.e. SM B cells, SM B cells and DN B cells) with Th1, Th2, Th17, Tfh, and Tregs, respectively. Interestingly, we found that the frequency of DN B cells was positively correlated with that of Th2 and Tregs (P = 0.01, P = 0.04, respectively) in patients, whereas there was no such correlation in HC (Figure 5). No correlation was found between the rest of altered B cell subsets and T cells in patients (data not shown).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/570888/fimmu-11-570888-HTML/image_m/fimmu-11-570888-g005.jpg)
Figure 5 Correlation analysis between the percentage of DN B cells and that of Th2 or Tregs in patients and HC. The graphs showed linear fit in patients. The data were analyzed with Pearsoncorrelation analysis. The dotted line demonstrated a 95% confidence interval. The upper panel, the graphs for DN B cells and Th2; the lower panel, the graphs for DN B cells and Tregs. Pearson’s correlation coefficient (r) and associated P values are shown on each graph.
Altered Serum Levels of IFN-γ, IL-4, IL-17A, and TGF-β1 in Sporotrichosis Patients
To evaluate the impact of the infection in the pattern of cytokines in sera, we next performed ELISA tests of IL-4, IL-17A, IFN-γ, and TGF-β1 in patients and HC (Figure S2A). Compared with HC, the levels of IFN-γ (P = 0.04, n = 21) and IL-4 (P = 0.008, n = 43) were significantly increased in patients, while that of IL-17A (P = 0.001, n = 34) and TGF-β1 (P = 0.001, n = 16) were greatly reduced. However, the concentrations of IFN-γ, IL-4, and IL-17A in patients with different duration had no statistical variation from those of HC (Figure S2B). With regard to TGF-β1, significant lower level was found both in patients with SD (P = 0.04, n = 5) and those with LD (P = 0.005, n = 11) (Figure S2C).
In Vitro Expansion of Specific CD4+ Teff Cells by S. globosa
To further understand T cell immunity against S. globosa infection, we established an in vitro co-culture assay to mimic in vivo interaction between PBMCs and this pathogen. PBMCs from healthy donors were cultured with/without live S. globosa yeast cells at different ratios (Yeasts: PBMCs, i.e. Y:P 1:1/2.5:1/5:1/10:1) for 3 days. The percentages of IFN-γ, IL-4, and IL-17 positive CD4+ T cells (i.e. Th1, Th2, and Th17 cells, respectively) were determined by flow cytometry (Figure S3A). Except the coculture at the highest Y/P ratio (10:1), all of these Teff cells’ frequencies in PBMCs were significantly elevated under the stimulation of yeast cells in a dose-dependent manner. When the load of pathogen to PBMCs reached 10:1, the proportions of these Teff cells were all comparable to that of control, however, Th1(P = 0.001, n = 3) and Th17 cells (P = 0.0005, n = 3) both notably dropped down nearly by half (Figure S3B).
Discussion
We herein described a previously unrecognized Th2 inclined response in circulation of patients with sporotrichosis caused by S. globosa: great increases in the percentage of Th2, the ratio of Th2/Tregs as well as IL-4 level in serum in whole patients. Except for the classical role in allergic diseases, Th2 cells are also found involved in the immunity to fungal infection. Th2 response was reported to begin in later phase in systemic S. schenckii infected mice (27), whereas it was significantly enhanced throughout the duration in our patients. We speculated that it might be due to differences in species and strains. Till now, only one published study mentioned that Th2 response was unchanged in S. globosa infection in human (18). Different number of subjects and the additional investigation of Tregs in our study may contribute to the discrepancy.
So far, most studies have shown that Th2-type immunity to fungi may increase susceptibility and prolong the survival of fungal pathogens partially by inhibiting Th1 and Th17 responses in the host (31). In our study, significantly increased IL-17A+ Th17 in circulation was accompanied with greatly reduced serum IL-17A, indicating a defect in Th17 function. In addition, it was reported that shift of the Th1-Th2 balance toward the Th2-dominant condition induces alternatively activated macrophages and results in worsened fungal infection (32). Similarly, we also found that the ratio of IFN-γ+ Th1/IL-4+ Th2 greatly decreased in patients suggesting the shift of Th1 to Th2 pattern in S. globosa infection. Taken together, the enhanced Th2 response may promote pathogenesis of sporotrichosis caused by S. globosa partly by suppressing defensive responses of Th1 and Th17. The underlying mechanism still needs further investigations.
Tregs, the potent immunosuppressive cells, play a pivotal role in maintenance of the homeostasis of immunity. Experimental evidence reveals a complex interaction between Tregs and Teff cells in fungal infections. In vitro expanded Tregs inhibited Th1 and Th2 responses but promoted Th17 cell response to C. albicans (33, 34). In S. schenckii-infected mice, an increased Tregs response was present at early period, but declined after 21 days post infection, while S. brasiiensis-infected mice showed a long lasting elevated Tregs response (35). Tregs depletion in mice was accompanied with enhanced Th1/Th17 response in schenckii-infection (36). However, Tregs profile in human sporotrichosis is still unclear. Our data demonstrated that Tregs were significantly reduced at early stage (in 6 months) and returned to normal level in later course (more than 6 months). Whereas, TGF-β1, the most important cytokine produced by Tregs, was remarkably decreased in serum regardless of duration suggesting a dysfunction of Tregs in LD patients infected with S. globosa. Taken together, in early stage of S. globosa infection, the immune status of the host skewed to inflammatory immunity, which was further supported by enhanced Th2/Tregs and Th17/Tregs in whole patients and increases of Teff cells.
Besides Tregs, different Teff cell profiles also seem to correspond to the inoculated strain with different virulence. Among the Sporothrix complex, S. brasiliensis is the most virulent one, followed by S. schenckii and then S. globosa (7–9). In vitro, the peptides of S. brasiliensis could elicit a mixed Th1/Th17 cell response and the cutaneous origin of S. schenckii induced potent Th1 reaction (15, 37, 38). In contrast, we found unchanged Th1 and increased dysfunctional IL-17A+Th17 cells in addition to prominent Th2 reaction while Chen reported great decrease in both Th1 and Th17 reactions in S. globosa infection (18). We supposed that S. globosa as a low virulent fungus may favor a Th2 environment instead of intense Th1 and Th17 response in the host.
For the first time, we found greatly altered distribution of B cell subpopulations in patients infected with S. globosa. Our data showed substantial increased frequencies of SM B cells and DN B cells and declined USM B cells along with unchanged naïve B cells in whole patients. Similar B cell profile was found in SD compared with HC, while only DN B subset’s percentage sustained high level at later time point. Previous studies revealed that USM B cells are responsible for IgM production and SM B cells are the main subtype producing IgG (30, 39, 40). However, only levels of IgG2 and IgG4 were reduced in patients with SD, with comparable levels of total IgG and IgM in patients to those of HC. Similar to our results, Esteban found that T. cruzi-infected individuals had significantly lower frequencies of SM B cells with no change in the level of total IgG (41). We speculated that decreased USM B cells may be responsible for decreased IgG2 and/or IgG4 production which need further verification. Nevertheless, since we did not measure the S. globosa specific antibodies, we could not rule out that these altered B cell subsets may affect the production of the pathogen specific antibodies. Meanwhile, considering the unchanged Tfh cells, we presumed that there may be other factors that interact with B cell subsets in this disease.
It has been shown that antigen activated B cells can influence CD4+ T cell differentiation (30). In vivo and in vitro, B cells promoted Th1 and Th17 cell responses against C. albicans infection in IL-6 dependent manner (42, 43). In addition, Th1 differentiation was suppressed by human B cells partially via IL-10 in systemic lupus erythematosus (SLE) patients (44). Therefore, we speculated that the remodeling of B cell subsets may be associated with the disturbed T cell reaction in S. globosa infection. Interestingly, for the first time, we found that the greatly increased DN B cells was positively correlated with Th2 and Tregs in all patients, while no such correlation was found in HC. DN B cells were also found expanded in SLE, HIV and malaria patients, respectively (45–47), but there is no literature mentioned the relation between DN B cells and T cells till now. Although the origin and function of DN B cells remain unclear, it was suggested that DN B cells may be exhausted memory B cells (48), or precursor memory B cells that have not yet to upregulate CD27 (49), or a product which CD27 fails to upregulate appropriately (50) and the increased DN B cells in SLE patients were considered as precursors of autoantibody producing plasma cells (51). The correlation found in our study implied that DN B cells collaborated with Th2 and Tregs in the immunity against this mycotic infection in human. However, besides their function, whether DN B cells regulated Th2 and Tregs differentiation or they were expanded by Th2 and Tregs in the immunity against S. globosa infection still needs further investigation.
According to clinical presentations, sporotrichosis are classified into fixed cutaneous, lymphocutaneous, disseminated cutaneous, and extracutaneous sporotrichosis, among which the first two types are the most common in clinic (52). However, little is known about the underlying immunological mechanism that may contribute to various clinical manifestations. In our study, in both FF and LF patients, we found similar changes in B cell compartment and elevated ratio of Th2/Tregs; however, additionally enhanced Th2 and Th17/Tregs were found in FF patients. Involvement of Th2 immune response to fungal infection usually results in susceptibility to infection partially owing to the inhibition of protective Th1 response by IL-4 and IL-10 (53), while immunity mediated by Th17 cells is protective as evidenced by its anti-fungal role in the defense of several fungi (54–56). Although the role of Th17 function seemed to be discounted in S. globosa infection, the different Th17/Tregs profile in the two clinical types indicated that lack of the help of Th17 response may contribute to the fungi spreading along lymphatic vessels in lymphocutaneous form rather than besieged in a single erythematous nodule as fixed type presented.
In summary, we first demonstrated the symphony of Th2 and remarkably altered B cell subsets in patients with sporotrichosis caused by S. globosa and these data will help us to further understand the complex immunity against this pathogen and may provide us new strategy for clinical diagnosis and therapy of this disease and similar fungal infection.
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
The studies involving human participants were reviewed and approved by the First Hospital of Jilin University Research Ethics Committee. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
JZ performed the FACS analysis and analyzed the data. LY, YS, and YC performed the collection of samples. MG and RC collected the patient information. YZ and SL supervised research and reviewed data and revised the paper. JZ and YZ wrote the paper with contributions from all authors. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by the grants from the National Natural Science Foundation of China (No. 81401351 and No.81773317) and Science and Technology Development Programs of Jilin Province (20180101110JC).
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.
Acknowledgments
The authors thank all patients and healthy donors who participated in current study.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2020.570888/full#supplementary-material
Supplementary Figure 1 | Immunophenotyping of circulating CD4+ Teff cells (Th1, Th2, Th17) in sporotrichosis patients and HC. Freshly isolated PBMCs were first stimulated with Protein Transport Inhibitor for 4 h in vitro. Cells were then subjected to CD4 staining followed by intracellular staining of IFN-γ, IL-4, IL-17A. (A) The subsets of CD4+ T cells were examined by flow cytometry. The numbers displayed are mean values of each quadrant. (B) Comparsion of the frequency of Teff cells in patients with HC. (C) Teff cell subsets are also compared between HC and patients with different duration. (D) The ratio of Th1/Th2 in whole patients, patients with different duration and HC. Error bars represent mean±SD. *P < 0.05, **P < 0.01, ***P < 0.001, and NS P ≥ 0.05.
Supplementary Figure 2 | Comparison of concentration of IFN-γ, IL-4, IL-17A, and TGF-β1 in sera of sporotrichosis patients and HC. (A) The statistical graphs for comparison of the cytokines between patients and HC. (B, C) The cytokine levels are analyzed between HC and subgroups of patients with different duration. Error bars represent mean±SD. *P < 0.05, **P < 0.01, and NS P ≥ 0.05.
Supplementary Figure 3 | Coculture of PBMCs and S. globosa in vitro. Freshly isolated PBMCs from healthy donors were co-cultured with live S. globosa yeast cells at different ratios [Yeasts:PBMCs (Y:P), 1:1/2.5:1/5:1/10:1]. After 3 days co-culture, cells were restimulated with Protein Transport Inhibitor for additional 4 h and then harvested and stained as indicated. (A) PBMCs were fixed, permeabilized and intracellularly stained for cytokines (IFN-γ, IL-4, IL-17A) after surface staining of CD4 mAb. The flow cytometry plot presents the average proportion of Th1, Th2, and Th17 subsets in CD4+ T cells of infected or control. Each graph is representative of three independent experiments (n = 3). (B) The frequencies of CD4+IFN-γ+, CD4+IL-4+, and CD4+IL-17A+ T cells are compared between the groups with yeasts and control. Error bars represent mean±SD. *P < 0.05, **P < 0.01, ***P < 0.001, and NS P ≥ 0.05.
References
1. Lopes-Bezerra LM, Mora-Montes HM, Zhang Y, Nino-Vega G, Rodrigues AM, de Camargo ZP, et al. Sporotrichosis between 1898 and 2017: The evolution of knowledge on a changeable disease and on emerging etiological agents. Med Mycol (2018) 56:126–43. doi: 10.1093/mmy/myx103
2. Chakrabarti A, Bonifaz A, Gutierrez-Galhardo MC, Mochizuki T, Li S. Global epidemiology of sporotrichosis. Med Mycol (2015) 53:3–14. doi: 10.1093/mmy/myu062
3. Moussa TAA, Kadasa NMS, Al Zahrani HS, Ahmed SA, Feng P, Gerrits van den Ende AHG, et al. Origin and distribution of Sporothrix globosa causing sapronoses in Asia. J Med Microbiol (2017) 66:560–69. doi: 10.1099/jmm.0.000451
4. Song Y, Li SS, Zhong SX, Liu YY, Yao L, Huo SS. Report of 457 sporotrichosis cases from Jilin province, northeast China, a serious endemic region. J Eur Acad Dermatol Venereol (2013) 27:313–8. doi: 10.1111/j.1468-3083.2011.04389.x
5. Yao L, Song Y, Cui Y, Zhou JF, Zhong SX, Zhao DY, et al. Pediatric Sporotrichosis in Jilin Province of China (2010-2016): a Retrospective Study of 704 Cases. J Pediatr Infect Dis Soc (2019) 9:342–45. doi: 10.1093/jpids/piz052
6. Zhao L, Cui Y, Zhen Y, Yao L, Shi Y, Song Y, et al. Genetic variation of Sporothrix globosa isolates from diverse geographic and clinical origins in China. Emerg Microbes Infect (2017) 6:e88. doi: 10.1038/emi.2017.75
7. Almeida-Paes R, de Oliveira LC, Oliveira MM, Gutierrez-Galhardo MC, Nosanchuk JD, Zancope-Oliveira RM. Phenotypic characteristics associated with virulence of clinical isolates from the Sporothrix complex. BioMed Res Int (2015) 2015:212308. doi: 10.1155/2015/212308
8. Della Terra PP, Rodrigues AM, Fernandes GF, Nishikaku AS, Burger E, de Camargo ZP. Exploring virulence and immunogenicity in the emerging pathogen Sporothrix brasiliensis. PloS Negl Trop Dis (2017) 11:e0005903. doi: 10.1371/journal.pntd.0005903
9. Mario DN, Schaffer LF, Peroza LR, Jesus FPK, Denardi LB, Fachinetto R, et al. Sporothrix brasiliensis produces the highest levels of oxidative stress in a murine model among the species of the Sporothrix schenckii complex. Rev Soc Bras Med Trop (2017) 50:554–57. doi: 10.1590/0037-8682-0171-2016
10. Tachibana T, Matsuyama T, Mitsuyama M. Involvement of CD4+ T cells and macrophages in acquired protection against infection with sporothrix schenckii in mice. Med Mycol (1999) 37:397–404. doi: 10.1046/j.1365-280x.1999.00239.x
11. Dickerson CL, Taylor RL, Drutz DJ. Susceptibility of congenitally athymic (nude) mice to sporotrichosis. Infect Immun (1983) 40:417–20. doi: 10.1128/IAI.40.1.417-420.1983
12. Verdan FF, Faleiros JC, Ferreira LS, Monnazzi LG, Maia DC, Tansine A, et al. Dendritic cell are able to differentially recognize Sporothrix schenckii antigens and promote Th1/Th17 response in vitro. Immunobiology (2012) 217:788–94. doi: 10.1016/j.imbio.2012.04.006
13. Koga T, Duan H, Urabe K, Furue M. Immunohistochemical localization of activated and mature CD83+ dendritic cells in granulomas of sporotrichosis. Eur J Dermatol (2001) 11:527–29. doi: 10.1080/mmy.40.2.111.114
14. Koga T, Duan H, Furue M. Immunohistochemical detection of interferon-gamma-producing cells in granuloma formation of sporotrichosis. Med Mycol (2002) 40:111–14. doi: 10.1080/mmy.40.2.111.114
15. de Almeida JRF, Jannuzzi GP, Kaihami GH, Breda LCD, Ferreira KS, de Almeida SR. An immunoproteomic approach revealing peptides from Sporothrix brasiliensis that induce a cellular immune response in subcutaneous sporotrichosis. Sci Rep (2018) 8:4192. doi: 10.1038/s41598-018-22709-8
16. Ferreira LS, Goncalves AC, Portuondo DL, Maia DC, Placeres MC, Batista-Duharte A, et al. Optimal clearance of Sporothrix schenckii requires an intact Th17 response in a mouse model of systemic infection. Immunobiology (2015) 220:985–92. doi: 10.1016/j.imbio.2015.02.009
17. Rossato L, Silvana Dos Santos S, Ferreira LG, Rogerio de Almeida S. The impact of the absence of Toll-like receptor-2 during Sporothrix brasiliensis infection. J Med Microbiol (2019) 68:87–94. doi: 10.1099/jmm.0.000876
18. Chen F, Yan BL, An L, Jiang RH. Is there a relationship between the percentage of T helper cells and prognosis for sporotrichosis? Br J Dermatol (2016) 174:1398–400. doi: 10.1111/bjd.14377
19. Garcia Carnero LC, Lozoya Perez NE, Gonzalez Hernandez SE, Martinez Alvarez JA. Immunity and Treatment of Sporotrichosis. J Fungi (Basel) (2018) 4:100. doi: 10.3390/jof4030100
20. Nascimento RC, Almeida SR. Humoral immune response against soluble and fractionate antigens in experimental sporotrichosis. FEMS Immunol Med Microbiol (2005) 43:241–7. doi: 10.1016/j.femsim.2004.08.004
21. Almeida SR. Therapeutic monoclonal antibody for sporotrichosis. Front Microbiol (2012) 3:409:409. doi: 10.3389/fmicb.2012.00409
22. de Almeida JR, Santiago KL, Kaihami GH, Maranhao AQ, de Macedo Brigido M, de Almeida SR. The Efficacy of Humanized Antibody against the Sporothrix Antigen, gp70, in Promoting Phagocytosis and Reducing Disease Burden. Front Microbiol (2017) 8:345. doi: 10.3389/fmicb.2017.00345
23. Portuondo DL, Batista-Duharte A, Ferreira LS, Martínez DT, Polesi MC, Duarte RA, et al. A cell wall protein-based vaccine candidate induce protective immune response against Sporothrix schenckii infection. Immunobiology (2016) 221:300–9. doi: 10.1016/j.imbio.2015.10.005
24. Téllez-Martínez D, Leandro Portuondo D, Loesch ML, Batista-Duharte A, Zeppone Carlos I. A Recombinant Enolase-Montanide™ PetGel A Vaccine Promotes a Protective Th1 Immune Response against a Highly Virulent Sporothrix schenckii by Toluene Exposure. Pharmaceutics (2019) 11:144. doi: 10.3390/pharmaceutics11030144
25. Almeida-Paes R, Pimenta MA, Monteiro PCF, Nosanchuk JD, Zancope-Oliveira RM. Immunoglobulins G, M, and A against Sporothrix schenckii Exoantigens in Patients with Sporotrichosis before and during Treatment with Itraconazole. Clin Vaccine Immunol (2007) 14:1149–57. doi: 10.1128/cvi.00149-07
26. Vinuesa CG, Linterman MA, Yu D, MacLennan IC. Follicular Helper T Cells. Annu Rev Immunol (2016) 34:335–68. doi: 10.1146/annurev-immunol-041015-055605
27. Maia DC, Sassa MF, Placeres MC, Carlos IZ. Influence of Th1/Th2 cytokines and nitric oxide in murine systemic infection induced by Sporothrix schenckii. Mycopathologia (2006) 161:11–9. doi: 10.1007/s11046-005-0142-y
28. Nascimento RC, Espíndola NM, Castro RA, Teixeira PA, Loureiro y Penha CV, Lopes-Bezerra LM, et al. Passive immunization with monoclonal antibody against a 70-kDa putative adhesin of Sporothrix schenckii induces protection in murine sporotrichosis. Eur J Immunol (2008) 38:3080–9. doi: 10.1002/eji.200838513
29. Kong F-Y, Feng B, Zhang H-H, Rao H-Y, Wang J-H, Cong X, et al. CD4+CXCR5+ T cells activate CD27+IgG+ B cells via IL-21 in patients with hepatitis C virus infection. Hepatobiliary Pancreatic Dis Int (2016) 15:55–64. doi: 10.1016/s1499-3872(16)60054-1
30. Agematsu K, Nagumo H, Shinozaki K, Hokibara S, Yasui K, Terada K, et al. Absence of IgD-CD27(+) memory B cell population in X-linked hyper-IgM syndrome. J Clin Invest (1998) 15:853–60. doi: 10.1172/JCI3409
31. Hansakon A, Mutthakalin P, Ngamskulrungroj P, Chayakulkeeree M, Angkasekwinai P. Cryptococcus neoformans and Cryptococcus gattii clinical isolates from Thailand display diverse phenotypic interactions with macrophages. Virulence (2019) 10:26–36. doi: 10.1080/21505594.2018.1556150
32. Arora S, Hernandez Y, Erb-Downward JR, McDonald RA, Toews GB, Huffnagle GB. Role of IFN-gamma in regulating T2 immunity and the development of alternatively activated macrophages during allergic bronchopulmonary mycosis. J Immunol (2005) 174:6346–56. doi: 10.4049/jimmunol.174.10.6346
33. Pandiyan P, Conti HR, Zheng L, Peterson AC, Mathern DR, Hernandez-Santos N, et al. CD4(+)CD25(+)Foxp3(+) regulatory T cells promote Th17 cells in vitro and enhance host resistance in mouse Candida albicans Th17 cell infection model. Immunity (2011) 34:422–34. doi: 10.1016/j.immuni.2011.03.002
34. Whibley N, Maccallum DM, Vickers MA, Zafreen S, Waldmann H, Hori S, et al. Expansion of Foxp3(+) T-cell populations by Candida albicans enhances both Th17-cell responses and fungal dissemination after intravenous challenge. Eur J Immunol (2014) 44:1069–83. doi: 10.1002/eji.201343604
35. Batista-Duharte A, Tellez-Martinez D, Roberto de Andrade C, Portuondo DL, Jellmayer JA, Polesi MC, et al. Sporothrix brasiliensis induces a more severe disease associated with sustained Th17 and regulatory T cells responses than Sporothrix schenckii sensu stricto in mice. Fungal Biol (2018) 122:1163–70. doi: 10.1016/j.funbio.2018.08.004
36. Batista-Duharte A, Téllez-Martínez D, de Andrade CR, Polesi MC, Portuondo DL, Carlos IZ. Transient Foxp3(+) regulatory T-cell depletion enhances protective Th1/Th17 immune response in murine sporotrichosis caused by Sporothrix schenckii. Immunobiology (2020) 225(5):151993. doi: 10.1016/j.imbio.2020.151993
37. Uenotsuchi T, Takeuchi S, Matsuda T, Urabe K, Koga T, Uchi H, et al. Differential induction of Th1-prone immunity by human dendritic cells activated with Sporothrix schenckii of cutaneous and visceral origins to determine their different virulence. Int Immunol (2006) 18:1637–46. doi: 10.1093/intimm/dxl097
38. Rodrigues AM, de Melo Teixeira M, de Hoog GS, Schubach TM, Pereira SA, Fernandes GF, et al. Phylogenetic analysis reveals a high prevalence of Sporothrix brasiliensis in feline sporotrichosis outbreaks. PloS Negl Trop Dis (2013) 7:e2281. doi: 10.1371/journal.pntd.0002281
39. Maurer D, Fischer GF, Fae I, Majdic O, Stuhlmeier K, Von Jeney N, et al. IgM and IgG but not cytokine secretion is restricted to the CD27+ B lymphocyte subset. J Immunol (1992) 148:3700–05.
40. Hu F, Zhang W, Shi L, Liu X, Jia Y, Xu L, et al. Impaired CD27(+)IgD(+) B Cells With Altered Gene Signature in Rheumatoid Arthritis. Front Immunol (2018) 9:626. doi: 10.3389/fimmu.2018.00626
41. Rodrigues MM, Fernández ER, Olivera GC, Quebrada Palacio LP, González MN, Hernandez-Vasquez Y, et al. Altered Distribution of Peripheral Blood Memory B Cells in Humans Chronically Infected with Trypanosoma cruzi. PloS One (2014) 9:e104951. doi: 10.1371/journal.pone.0104951
42. Shen P, Fillatreau S. Antibody-independent functions of B cells: a focus on cytokines. Nat Rev Immunol (2015) 15:441–51. doi: 10.1038/nri3857
43. Li R, Rezk A, Li H, Gommerman JL, Prat A, Bar-Or A, et al. Antibody-Independent Function of Human B Cells Contributes to Antifungal T Cell Responses. J Immunol (2017) 198:3245–54. doi: 10.4049/jimmunol.1601572
44. Blair PA, Noreña LY, Flores-Borja F, Rawlings DJ, Isenberg DA, Ehrenstein MR, et al. CD19+CD24hiCD38hi B Cells Exhibit Regulatory Capacity in Healthy Individuals but Are Functionally Impaired in Systemic Lupus Erythematosus Patients. Immunity (2010) 32:129–40. doi: 10.1016/j.immuni.2009.11.009
45. Wei C, Anolik J, Cappione A, Zheng B, Pugh-Bernard A, Brooks J, et al. A new population of cells lacking expression of CD27 represents a notable component of the B cell memory compartment in systemic lupus erythematosus. J Immunol (2007) 178:6624–33. doi: 10.4049/jimmunol.178.10.6624
46. Moir S, Ho J, Malaspina A, Wang W, DiPoto AC, O’Shea MA, et al. Evidence for HIV-associated B cell exhaustion in a dysfunctional memory B cell compartment in HIV-infected viremic individuals. J Exp Med (2008) 205:1797–805. doi: 10.1084/jem.20072683
47. Isnardi I, Ng Y-S, Menard L, Meyers G, Saadoun D, Srdanovic I, et al. Complement receptor 2/CD21- human naive B cells contain mostly autoreactive unresponsive clones. Blood (2010) 115:5026–36. doi: 10.1182/blood-2009-09-243071
48. Centuori SM, Gomes CJ, Kim SS, Putnam CW, Larsen BT, Garland LL, et al. Double-negative (CD27(-)IgD(-)) B cells are expanded in NSCLC and inversely correlate with affinity-matured B cell populations. J Transl Med (2018) 16:30. doi: 10.1186/s12967-018-1404-z
49. Fraussen J, Marquez S, Takata K, Beckers L, Montes Diaz G, Zografou C, et al. Phenotypic and Ig Repertoire Analyses Indicate a Common Origin of IgD(-)CD27(-) Double Negative B Cells in Healthy Individuals and Multiple Sclerosis Patients. J Immunol (2019) 203:1650–64. doi: 10.4049/jimmunol.1801236
50. Agematsu K. Memory B cells and CD27. Histol Histopathol (2000) 15:573–76. doi: 10.14670/HH-15.573
51. Jenks SA, Cashman KS, Zumaquero E, Marigorta UM, Patel AV, Wang X, et al. Distinct Effector B Cells Induced by Unregulated Toll-like Receptor 7 Contribute to Pathogenic Responses in Systemic Lupus Erythematosus. Immunity (2018) 49:203. doi: 10.1016/j.immuni.2018.08.015
52. Morgado FN, Schubach AO, Barros MB, Conceicao-Silva F. The in situ inflammatory profile of lymphocutaneous and fixed forms of human sporotrichosis. Med Mycol (2011) 49:612–20. doi: 10.3109/13693786.2011.552532
53. Blanco JL, Garcia ME. Immune response to fungal infections. Vet Immunol Immunopathol (2008) 125:47–70. doi: 10.1016/j.vetimm.2008.04.020
54. Silva AA, Criado PR, Nunes RS, da Silva WL, Kanashiro-Galo L, Duarte MI, et al. In situ immune response in human chromoblastomycosis–a possible role for regulatory and Th17 T cells. PloS Negl Trop Dis (2014) 8:e3162. doi: 10.1371/journal.pntd.0003162
55. Bozza S, Gaziano R, Spreca A, Bacci A, Montagnoli C, di Francesco P, et al. Dendritic cells transport conidia and hyphae of Aspergillus fumigatus from the airways to the draining lymph nodes and initiate disparate Th responses to the fungus. J Immunol (2002) 168:1362–71. doi: 10.4049/jimmunol.168.3.1362
Keywords: sporotrichosis, Sporothrix globosa, immunity, T cells, B cells
Citation: Zu J, Yao L, Song Y, Cui Y, Guan M, Chen R, Zhen Y and Li S (2020) Th2 Biased Immunity With Altered B Cell Profiles in Circulation of Patients With Sporotrichosis Caused by Sporothrix globosa. Front. Immunol. 11:570888. doi: 10.3389/fimmu.2020.570888
Received: 10 June 2020; Accepted: 16 October 2020;
Published: 13 November 2020.
Edited by:
Juarez Antonio Simões Quaresma, Evandro Chagas Institute, BrazilReviewed by:
Alexander Batista Duharte, Universidade Estadual Paulista, BrazilMarcos Vinícius Da Silva, Universidade Federal do Triângulo Mineiro, Brazil
Copyright © 2020 Zu, Yao, Song, Cui, Guan, Chen, Zhen and Li. 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: Shanshan Li, shansalee@163.com; Yu Zhen, zyu2007@hotmail.com