- 1Image and Graphic Methods Unit, University Hospital, Federal Sergipe University, Aracaju, Brazil
- 2Health Science Graduate Program, Federal University of Sergipe, Aracaju, Brazil
- 3Laboratory of Molecular Biology, University Hospital, Federal Sergipe University, Aracaju, Brazil
- 4Department of Medicine, University Hospital, Federal University of Sergipe, Aracaju, Brazil
- 5Hepatology Service, University Hospital, Federal University of Sergipe, Aracaju, Brazil
- 6Department of Biochemistry and Immunology, Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, Brazil
- 7Division of Pathology, University Hospital, Federal University of Sergipe, Aracaju, Brazil
- 8Infectology Service, University Hospital, Federal University of Sergipe, Sergipe, Brazil
- 9Immunology Institute of Investigation (III), National Institute of Science and Technology (INCT), Brazilian Research and Technology Council (CNPq), São Paulo, Brazil
- 10Department of Morphology, Federal University of Sergipe, São Cristóvão, Brazil
This is a case series study to evaluate immunological markers associated with schistosomiasis advanced fibrosis, including 69 patients from an endemic area from the State of Sergipe and from the Hepatology Service of the University Hospital in Sergipe, Brazil. Hepatic fibrosis was classified based on Niamey protocol for ultrasonography (US). Immune response to Schistosoma mansoni antigens was evaluated by stimulating peripheral blood mononuclear cells (PBMCs) from these patients with either adult worm (SWAP—10 μg/ml) or egg (SEA—10 μg/ml) antigens or purified protein derivative of turberculin (PPD—10 μg/ml) or phytohemagglutinin (PHA—1 μg/ml) for 72 h. The levels of IFN-γ, TNF-α, IL-5, IL-10, and IL-17 were measured in these supernatants by ELISA and IL-9 by Luminex. Single nucleotide polymorphisms in IL-17, IL10, and CD209 genes were genotyped using TaqMan probe by qPCR. Higher levels of IL-9, IL-10, and IL-17 were found in PBMC supernatants of patients with advanced hepatic fibrosis. Direct correlations were detected between IL-9 and IL-17 levels with US spleen sizes, portal vein diameters, and periportal thickening. The CD209 rs2287886 AG polymorphism patients produce higher IL-17 levels. Together, these data suggest a role of these cytokines in the immunopathogenesis of advanced fibrosis in human schistosomiasis.
Introduction
Schistosomiasis is an infectious and parasitic disease which affects about 240 million people in 78 countries and around 700 million people live in areas at risk from this infection (1). In Brazil, schistosomiasis is mainly caused by Schistosoma mansoni (S. mansoni). Most infected individuals develop mild gastrointestinal disease, but 5% to 10% develop hepatosplenic schistosomiasis, characterized by severe liver fibrosis, splenomegaly, and portal hypertension (2). About 22 million people are currently suffering from complications of chronic S. mansoni infection and up to 42% of infected people have been diagnosed with periportal fibrosis. About 0.2 million deaths are attributed to chronic schistosomiasis every year (3). Severe liver fibrosis and portal hypertension hold the main responsibility for morbidity and mortality related to schistosomiasis and high cost for the healthcare and control programs of countries, placing this disease as an important public health concern (4, 5).
The mechanisms related to the development of severe forms of the disease in humans have not been clearly defined yet (6). Schistosomiasis morbidity has varied directly with egg count, but the nature of the immune response of the host is what shapes the magnitude of the immunopathology in each case (6). Immune response mediated by CD4+ T cells from the host to S. mansoni eggs laid in the liver and intestine, with granuloma formation, leads to fibrosis, with consequent portal hypertension (7). Studies in both, experimental models and in human, show that the initial events of granuloma formation are associated with a Th2 response (8). Additionally, in an experimental model with an advanced fibrosis resembling advanced hepatic fibrosis from human schistosomiasis, a role of Th17 response is demonstrated, whereas in human schistosomiasis, there is no clear understanding of what occurs in the later stages of fibrosis (9). Furthermore, the genetic background of the host might be involved in the modulation of immune response and pathogenesis development, as seen in many infectious diseases (10–15). Thus, functional single nucleotide polymorphisms (SNPs) may be important as they may influence the gene expression or the structure of the corresponding protein (16). The understanding of the molecular mechanisms involved in the immunopathogenesis of this disease can have important implications for intervention strategies, mainly in immunoprophylaxis, not only in schistosomiasis but also in other parasitic diseases (5, 17).
he present study evaluates the cytokine profiles of Th1 (IFN-γ, TNF-α), Th2/Treg (IL-5, IL-10), Th9 (IL-9), and Th17 (IL-17) in peripheral blood mononuclear cell (PBMC) supernatants stimulated with soluble egg (SEA) and adult worm (SWAP) antigens from patients with schistosomiasis with different stages of hepatic fibrosis and investigates whether SNPs in IL17A, IL10, and CD209 receptor [dendritic cell-specific ICAM3-grabbing non-integrin (DC-SIGN)] are associated with this pathogenesis and influence the expression of IL-10 and IL-17 in this population.
Materials and Methods
Ethics Statement
This study was carried out in accordance with the recommendations of the Declaration of Helsinki. All donors or their legal guardians from whom information and biologic samples were obtained or clinical procedures were performed gave written informed consent for their participation in the study. The study was approved by the Ethics and Research Committee from the University Hospital of the Federal University of Sergipe, Brazil (CAAE 0022.0.107.000-08).
Study Design and Patient Selection
A case series analytical study was performed in patients with schistosomiasis, comparing patients classified according to the degree of fibrosis by ultrasonography (US). A total of 69 patients were included. Most patients were recruited from the city of Ilha das Flores, an important endemic area for schistosomiasis, located in the north of the State of Sergipe, Brazil, on the banks of the São Francisco River. According to the 2007 Schistosomiasis Control Program data, there is a high frequency (46.45%) in the positivity of parasitological stool test in this municipality, provided by constant and obligatory contacts with natural sources of water, which are abundant in the municipality, whether in work activities or in daily household chores. As only a mild hepatosplenic form was present in about 2% to 7% of infected patients from this endemic area (18), patients treated at the Hepatology Service of the University Hospital (UH) at the Federal University of Sergipe were recruited. These patients were from other endemic areas of the state or even from neighboring states, such as Alagoas and Bahia, and presented severe hepatosplenic schistosomiasis, but did not have any other associated diseases. Therefore, it was possible to analyze the immune response associated to more advanced stages of the disease.
The inclusion criteria for the study were as follows: positive parasitological diagnosis based on the Kato-Katz method for S. mansoni, or previous treatment background, classified using abdominal ultrasonography with different degrees of fibrosis compatible with S. mansoni etiology. Exclusion criteria were age <5 and >70 years old, pregnant or people with diseases which affect the immune response (HIV, HTLV-1, and diabetes mellitus), people with other liver diseases associated with portal hypertension (such as hepatitis A, B, or C; cirrhosis; Budd–Chiari; portal vein thrombosis) and acute viral or bacterial infections in the moment of blood collection, and patients who have not completed the immunological evaluation. After applying the exclusion criteria, the selected patients were divided into two groups: no fibrosis (n = 53), all of them from the endemic area of Ilha das Flores, and with advanced fibrosis (n = 16), all of them from the Hepatology Service of the University Hospital (UH). The small sample size in the advanced fibrosis group is due to the rigidus selection criteria, which include only patients with confirmed hepatic fibrosis and portal hypertension, also confirmed by esophageal varices in the endoscopy, and exclude patients with comorbidities and coinfections that are very frequent in these group of patients. The small sample might reduce the power and interfere in the negative results, but the significant data are still valid.
After all, a genetic polymorphism study was performed in both groups. Unfortunately, due to problems of DNA recovery from the blood samples and low quantification of DNA, only a total of 37 of these individuals were enrolled in the SNP analysis: 28 patients from the no fibrosis group and 9 patients from the fibrosis group.
Ultrasonography
Ultrasound exams (US) were done following the method defined by the World Health Organization at Niamey meeting, 2000, for the classification of hepatic fibrosis. The US was performed by the same sonographer on a high-resolution device, Medison brand, model: Sonoace 8000 SE. The Niamey classification consists of three stages, which result in three scores that assess intraparenchymal (IP), periportal (PT), and portal hypertension (PH) aspects. The assessment of periportal fibrosis in this protocol is analyzed by two distinct methods: a qualitative one, which considers the ultrasound aspect of the liver, evaluated by characteristic patterns, the IP score, associated with a quantitative method, which results from the measurement of two or three second-order portal branches, the PT score. In addition to these two criteria, the presence of portal hypertension (PH score) is also analyzed, which is given by the analysis of indicators of its existence, as the measure of the portal vein diameter, the presence of ascites, and the presence of collateral circulation. These scores (IP, PT, and PH) are interpreted using a table that classifies the periportal fibrosis as follows: group 1 (no sign of periportal fibrosis), group 2 (incipient periportal fibrosis not excluded), group 3 (possible periportal fibrosis), group 4 (probable periportal fibrosis), group 5 (periportal fibrosis), group 6 (advanced periportal fibrosis), and group 7 (advanced periportal fibrosis + portal hypertension). This method has been used in several previous studies (19–23).
Immunological Evaluation
The specific antigens of S. mansoni used were soluble extracts of adult worm (SWAP) and soluble egg antigen (SEA) which were kindly provided by Evan Secor from the Division of Parasitic Diseases and Malaria, Centers for Disease Control and Prevention, Atlanta, Georgia. For the isolation of PBMC from heparinized blood, it was subjected to a density gradient using a Ficoll Hypaque gradient (Ficoll-Paque™ PLUS; GE Healthcare, USA), as previously described (24, 25). In summary, 3 × 106 cells in 1 ml of RPMI 1640 (Gibco™, Thermo Fisher Scientific, Brazil) supplemented with 10% AB Rh-positive serum were either not stimulated or stimulated with specific S. mansoni antigens. The antigen concentration of 10 μg/ml for both SEA and SWAP was chosen because it was able to induce a response in PBMC from patients with schistosomiasis, but not from healthy controls from the same endemic area. The mitogen phytohemagglutinin (PHA) at 1 μg/ml and the purified protein derivative of turberculin (PPD) at 10 μg/ml were used as positive controls. After a 72-h incubation at 37°C in 5% CO2, the plates were centrifuged and the supernatants collected, aliquoted, and stored at −80°C for later cytokine dosage.
Determination of Cytokine Levels
The levels of IFN-γ and TNF-α (Th1 profile), IL-5 and IL-10 (Th2/Treg profiles), and IL-17 (Th17 profile) cytokines were measured in these supernatants by the ELISA technique using commercially available kits (DuoSet, R&D Systems Inc., USA) on microtiter plates (Kartell SPA, Italy). The techniques described by the manufacturers were followed, except that the diluent for establishing the standard curves was the same (RPMI + 10% AB-positive sera) as for the tested samples. IL-9 (Th9) levels were measured by Luminex (Milliplex kit, Panomics Inc., Affymetrix, USA). IL-9 was only analyzed in SEA-stimulated supernatants because of the limitation of the number of samples in the kit available. The measurement techniques were also described in previous publications of our group (5, 26).The results were expressed in pg/ml based on a standard curve using the recombinant cytokines. It was not possible to measure all cytokines in all patients because we run out of some of the PBMC supernatants.
SNP Genotyping
SNPs in IL17A, IL10, and CD209 (DC-SIGN) gene were genotyped by allelic discriminations. In brief, genomic DNA was extracted using PureLink® Genomic DNA Mini Kit (Invitrogen™, USA). The concentration and purity of DNA were quantified using NanoDrop™ (Thermo Scientific™, USA). We genotyped DNA samples using IL17A rs2275913, IL10 rs1800871 and rs1800872, and CD209 rs2287886 and rs4804803 TaqMan® probe by qPCR using 7500 Real-Time PCR (Applied Biosystems, USA) following the instructions of the manufacturer. The results were assessed using TaqMan® Genotyper software version 1.6.0 (Applied Biosystems, USA). Information about the analyzed SNPs are presented in Supplementary Table 1.
Statistical Analyses
D’Agostino & Pearson (D’Agostino & Pearson omnibus normality test) and Shapiro–Wilk normality tests were applied to the numerical data of the groups. Differences between cytokine levels in PBMC supernatants were compared between groups with different degrees of fibrosis using the Kruskal–Wallis non-parametric test, with Dunn’s post-test. Correlations were also made between numerical variables such as cytokine levels with measurements of periportal thickening, portal vein diameter, and spleen size by Spearman’s correlation test. Allelic and genotype frequency associations in the SNPs were analyzed by Fisher exact test or Chi-square test. Odds ratio value was used to compare the groups. These analyses were performed using the GraphPad Prism program, version 6.0. A value of α ≤5% (p ≤ 0.05) was considered for statistical significance. The Hardy–Weinberg equilibrium test was performed by GENEPOP online 4.2 available at http://genepop.curtin.edu.au/.
Results
Demographic and Clinical Aspects of Individuals From the Endemic Area and From the University Hospital
The demographic and clinical data of patients recruited from the endemic area and from the UH, classified according to their clinical forms, are shown in Table 1.
Table 1 Demographic, clinical, and ultrasound data of individuals from the endemic area and hepatosplenic ones from the University Hospital (UH) classified according to their clinical forms.
It is observed that hepatosplenic patients (HE) from the UH are significantly older than the ones with intestinal (I) and hepatointestinal (HI) forms from the endemic area. All patients from the Hepatology Service of UH (n = 16) had advanced fibrosis on ultrasound, associated with portal hypertension (Niamey 7), mainly verified by the measurements of the spleen size, the portal vein diameter, and the periportal space and showed portal hypertension signs, confirmed by detection of esophageal varices in upper digestive endoscopy. Patients from the endemic area (n = 53) did not have advanced fibrosis, being classified as Niamey 1.
Cellular Immune Response in Individuals Classified According to the Degrees of Fibrosis Determined by Ultrasound
The levels of cytokines IFN-γ, TNF-α, IL-5, IL-9, IL-10, and IL-17 were compared in the groups of patients without fibrosis (classified as Niamey 1) (n = 53) and with advanced fibrosis (classified as Niamey 7) (n = 16) (Figure 1). IL-17 levels were significantly higher in PBMC supernatants stimulated with SWAP and SEA in patients with fibrosis (Figures 1B, C). IL-9 and IL-10 levels were significantly higher in the PBMC supernatants stimulated with SEA in patients with fibrosis (Figures 1A, K). There were no significant differences in the levels of IFN-γ, TNF-α, and IL-5 in the PBMC supernatants stimulated with SEA and SWAP (Figures 1D–I). All patients presented a response to PPD and PHA stimulations.
Figure 1 Peripheral blood mononuclear cell (PBMC)-specific cytokines response to Schistosoma mansoni antigens in patients with different degrees of fibrosis. Comparisons between the levels of IL-9 (A), IL-17 (B, C), IFN-γ (D, E), TNF-α (F, G), IL-5 (H, I), and IL-10 (J, K) in supernatants from PBMC cultures stimulated with adult worm (SWAP) and egg (SEA) antigens in individuals with different degrees of fibrosis by the Niamey criteria (no fibrosis and advanced fibrosis). IL-9 was only analyzed in SEA-stimulated supernatants because of the limitation of the number of samples in the kit available. The data represent individual values and the averages are shown. Comparisons were analyzed using the Mann–Whitney U test.
Correlations Between Cellular Immune Response and Ultrasound Data Associated With Fibrosis and Portal Hypertension
Correlations were made between the levels of IFN-γ, TNF-α, IL-5, IL-9, IL-10, and IL-17 cytokines and the average of the measurements of the periportal space, the portal vein diameter, and the spleen, taken by ultrasound (Figure 2). There was a direct correlation between these measurements and the levels of IL-9 and IL-17 cytokines in response to SEA (Figures 2A–F). There was also a direct correlation between IL-10 levels in response to SEA and the average measurements of periportal space (Figure 2G). IFN-γ in response to SEA showed an inverse correlation with periportal thickening and spleen size (Figures 2H, I). There were no significant correlations between these measurements and the levels of TNF-α and IL-5 (Supplementary Table 2).
Figure 2 Correlations between the response of peripheral blood mononuclear cell (PBMC)-specific cytokines to S. mansoni antigens and the ultrasound data associated with fibrosis and portal hypertension. Correlations between the levels of IL-9 (A–C), IL-17 (D–F), IL-10 (G), and IFN-γ (H, I) in the supernatants from PBMC cultures stimulated with soluble egg antigen (SEA) and the periportal thickening, portal vein diameter, and the spleen size measured on ultrasound. Correlations were analyzed using Spearman’s correlation test.
Genotyping of IL17A, IL10, and CD209 SNPs
The SNPs in IL17A (rs2275913), IL10 (rs1800871 and rs1800872), and CD209 (rs2287886 and rs4804803) genes were genotyped in the no fibrosis patients (classified as Niamey 1) (n = 28) and fibrosis patients (classified as Niamey 7) (n = 9). No deviation was found in the Hardy–Weinberg equilibrium test. There were no significant differences in genotype and allelic comparisons between schistosomiasis patients with fibrosis and without fibrosis, despite the differences in genotype frequency distributions of IL-10 (rs1800871 and rs1800872) and CD209 (rs2287886 and rs4804803) SNPs (Supplementary Table 3).
CD209 rs2287886 SNP Is Associated With Differences in IL-17 Levels in Schistosomiasis Patients
When we compared the amount of IL-17 quantified in the PBMC supernatants of all patients included in the genetic analysis together, stratified according to genotypes, we found an association between higher levels of IL-17 in the AG genotype to CD209 rs2287886 when compared with the GG genotype in non-stimulated PBMC supernatants (p = 0.051) but not in the SWAP (p = 0.063) nor SEA (p = 0.964)-stimulated PBMC supernatants (Mann–Whitney test—Supplementary Figure 1). There were no other differences between the genotype of other SNPs assessed and the cytokine measurement found in the PBMC supernatants in this population when all patients included in the genetic analysis were evaluated together.
Complementary to this, we analyzed the levels of IL-17 and IL-10 in non-stimulated and stimulated with SWAP or SEA PBMC supernatants according to different genotypes for the SNPs studied in each group of patients separately (no fibrosis and advanced fibrosis) (Supplementary Figure 2). We observed significant differences between cytokine dosages when advanced fibrosis and no fibrosis patients with the same genotype for the SNPs were compared (Supplementary Figure 2).
Discussion
The present study detects a predominance of IL-17 (Th17) and IL-9 (Th9) production in hepatosplenic patients with advanced fibrosis in PBMC cultures stimulated with S. mansoni antigens. In addition, there is a direct correlation between these cytokine levels and data related to the severity of fibrosis, which suggests an important association of the Th17 and Th9 response with the severe form of human schistosomiasis. Besides, the study demonstrates a predominance of IL-10 production in hepatosplenic patients with advanced fibrosis in cultures stimulated with SEA and a correlation between this response and the average measurements of periportal space, indicating a possible Th2/Treg response remaining from the initial periportal fibrosis or an attempt to downmodulate the fibrosis.
Many studies show the importance of the immune response to parasite antigens in the pathogenesis of schistosomiasis. Studies in experimental models have shown that cytokines are important factors in the formation and modulation of the immune response to S. mansoni egg. Most experimental studies demonstrate the role of Th2 cytokines (IL-4, IL-5, and IL-13) in the formation of granuloma (27–29) and of the Th17 response in advanced fibrosis (9, 30). In human schistosomiasis, the influence of the Th2 response and the TGF-β cytokine in periportal fibrosis has been demonstrated (5), but they have failed to identify the immune response profile associated with advanced fibrosis, because the studied patients did not present advanced fibrosis. To our knowledge, this is the first study that demonstrates the association of these cytokines with advanced fibrosis in human schistosomiasis.
The present study observes a predominance of IL-17 (Th17) production in hepatosplenic patients with advanced fibrosis in cultures stimulated with S. mansoni antigens. In addition, there is a direct correlation between these cytokine levels with the spleen size, portal vein diameter, and periportal thickening, data related to the severity of the disease, which suggests an important association of the Th17 response with the pathogenesis of the severe form of human schistosomiasis. These findings are in accordance with literature data in an experimental model that develops severe fibrosis, corresponding to advanced human fibrosis, known as Symmers’ fibrosis. In these mice, the IL-23 and the induction of IL-17 have been shown to play a critical role in the immunopathology of schistosomiasis (31). Kalantari et al. (32) demonstrated that CD209a-deficient CBA mice display decreased Th17 responses and are protected from severe immunopathology. A recent study found that acute schistosome infection induced a transient Th17 response to cathepsin B1 cysteine proteases secreted by the worms and that this early Th17 response may determine the pathogenic progression of the infection (7, 33). The Th17 response has also been associated with the pathogenesis of other human helminth infections (34).
The present study also demonstrates a predominance of IL-9 production in hepatosplenic patients with advanced fibrosis in SEA-stimulated cultures, in addition to correlations between this response and the spleen size, portal vein diameter, and periportal thickening, thus suggesting a role of the Th9 response in the pathogenesis of the severe form of human schistosomiasis. In mice infected with S. japonicum, proliferation of Th9 cells and elevated levels of IL-9 have been demonstrated. In addition, dynamic changes of Th9 and IL-9 were synchronous with the developmental trend of hepatic egg granulomatous inflammation, suggesting that Th9 cells might be a new subset in the pathogenesis of schistosomiasis (35). Several studies demonstrate the role of this cytokine in several classes of inflammatory diseases. Treatment with IL-9 neutralizing antibodies effectively attenuates liver inflammation and necrosis (36). In allergic airway disease, allergic patients have more circulating Th9 cells than non-allergic control subjects (37). Neutralization of IL-9 was able to reduce allergen-induced inflammation (38).
The only study carried out in humans by Barreto et al. (39), evaluating the association of IL-17 and IL-9 in sera from patients with several clinical forms of schistosomiasis, did not show significant differences in the levels of these cytokines in patients with different degrees of fibrosis. However, in this study, the cytokines were measured in serum samples of these patients, and they did not evaluate the specific immune response to S. mansoni antigens. More importantly, it is noteworthy that the study only included two patients with advanced fibrosis.
The present study shows that there is a difference in the IL-10 (Th2/Treg) levels in patients with advanced fibrosis in the supernatants of PBMC cultures when stimulated with SEA. There are no differences in the response to the SWAP, demonstrating that this response is specific to the egg antigen. Previous studies in experimental models demonstrate the IL-10 induction mediated by carbohydrates of the Schistosoma eggshell (40, 41). We also observed a direct correlation between these cytokine levels and periportal thickening, with no correlations between IL-10 and the US data most associated with portal hypertension (portal vein diameter and spleen size). IL-10 is known for its regulatory function, and some experimental studies show an association between its absence and accelerated liver fibrosis in schistosomiasis (42). A study in mice showed that IL-10 is an important immunoregulatory cytokine in acute schistosomiasis, but it plays no essential role in the process of immune downmodulation in chronic S. mansoni infection in mice (43). Jesus et al. (5) also observed an increase in IL-10 in patients with degree III hepatic fibrosis, representing the initial stages of periportal fibrosis, but in that study, there were no patients with advanced fibrosis. Our hypothesis is that the Th2/Treg response shown to occur in the early stages of liver fibrosis and reported in experimental models (44, 45) and the high IL-10 levels found in the present study in patients with advanced fibrosis represent an attempt to modulate the injury-inducing immune response but are insufficient to prevent the progression of fibrosis.
Considering previous reports of functional SNPs in these cytokine genes, in this study, the IL17A rs2275913, IL10 rs1800871, IL10 rs1800872, CD209 rs2287886, and CD209 rs4804803 were genotyped. Although no statistical differences were in the allelic and genotype frequencies comparisons, probably due to the small sample size, differences in the genotype frequency distributions of IL-10 (rs1800871 and rs1800872) and CD209 (rs2287886 and rs4804803) SNPs were observed between schistosomiasis patients with fibrosis and without fibrosis.
Here, we found an association between the presence of genotype AG for the CD209 rs2287886 with higher levels of IL-17 in non-stimulated PBMC supernatant when compared with the GG genotype in schistosomiasis patients.
The CD209 receptor, also known by DC-SIGN, is part of a human C-type lectin receptor family and has as primary function in the recognition of conserved carbohydrate domains that are present in many pathogens, including S. mansoni (46). Its recognition leads DCs to produce IL-1β and IL-23 and induces Th17 responses (9). Experimental studies demonstrated that, despite similar infection protocols, there is a variation in disease severity among different mouse strains regarding schistosomiasis. While CBA mice develop large egg granulomas and severe disease due to CD209a activation to produce IL-1β and IL-23 and to induce Th17 cells, BL/6 mice develop smaller lesions, explained by a disability to induce Th17 response by the CD209a recognition pathway (47, 48).
The A allele in human CD209 rs2287886 is associated with a higher expression of CD209 in DC from healthy donors, when compared with the G allele. Interestingly, the CDs that expressed significantly higher levels of CD209 were also more susceptible to cytomegalovirus infection (13). No patient from the fibrosis group presented the AA genotype, but a higher frequency of the AG genotype was observed in this group. Our findings that the genotype AG is associated with higher levels of IL-17 agree with a possible explanation derived from previous studies in a murine model, since the increased expression of CD209 receptor in CDs is associated with induction of Th17 response (9, 47, 48).
The lack of associations between the IL17A rs2275913, IL10 rs1800871, IL10 rs1800872, CD209 rs2287886, and CD209 rs4804803 alleles or genotype frequencies with the schistosomiasis clinical manifestations in this population could be due to small sample size, and further studies are necessary to rule out the influence of these SNPs in schistosomiasis pathogenesis. In spite of this, we could observe significant differences between cytokine dosages when advanced fibrosis and no fibrosis patients with the same genotype for the SNPs were compared, suggesting that this modulation may occur as a consequence of disease severity variations and not due to genetic background of the SNPs assessed.
There is no association between the levels of IL-5, TNF-α, and INF-γ and the advanced degree of fibrosis in the present study. Most studies in mice clearly show the role of Th2-type cytokines such as IL-4 (27) and IL-13 (29) in determining granulomatous lesion and hepatic fibrosis in schistosomiasis, with IL-5 being important for eosinophilia present in these granulomas (49). Shainheit et al. (50) evidenced that there was no fibrosis progression after blocking the production of INF-γ in mice. However, these mice do not have the classic Symmers’ fibrosis and portal hypertension from the schistosomiasis observed in the human species. In another mouse model, which presents fibrosis more like that which occurs in humans, it was evidenced that cytokines such as TNF-α (51) and IL-17 (31, 52) have an important role in this fibrosis process. In the human species, it has been observed that, similar to mice, the immune response of the Th2 profile (IL-5, IL-13) in response to egg antigen (SEA) is associated with the early stages of fibrosis, using the classification of Abdel-Wahab (5). However, in this study, there were no patients with more advanced hepatosplenic form. A study performed in Africa evidenced higher levels of TNF-α and IFN-γ and lower levels of IL-5 in PBMC supernatants stimulated with S. mansoni egg antigen in hepatosplenic individuals when compared with non-hepatosplenic ones (53). However, the patients in this study were from a malaria-endemic area and were classified only by clinical examination, and the hepatosplenomegaly may also be due to malaria. Malaria infection, in turn, could also justify the high levels of TNF-α and IFN-γ found in these patients (54).
In conclusion, the severity of fibrosis observed by ultrasonography in human schistosomiasis is associated with a Th9 and Th17 response specific to S. mansoni antigens. The direct correlation between the levels of IL-9 and IL-17 cytokines and the spleen size, portal vein diameter, and thickening of the periportal space reinforces the association of these cytokines with the immunopathogenesis of human schistosomiasis. The observation that IL-10 is also increased in patients with advanced fibrosis and the correlation found between IL-10 levels only with thickening of the periportal space suggest the hypothesis of an attempt of fibrosis modulation by this cytokine.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, under a reasonable request.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics and Research Committee from the University Hospital of the Federal University of Sergipe, Brazil. Written informed consent to participate in this study was provided by the legal guardian/next of kin of the participants.
Author Contributions
AJ, AS, MS, CR, and FA recruited and clinically characterized the schistosomiasis patients from the endemic area of Ilha das Flores. AF, AS, and FL recruited and clinically characterized the patients from the University Hospital of Federal University of Sergipe. KF and DF performed the ultrasound exams and analyzed and interpreted the data related to the classification of hepatic fibrosis. FA, CR, FO, CS, and LM performed the immunological evaluation experiments. CS and RC performed the SNP genotyping experiments and analyzed and interpreted the data. AJ and RA conceived and designed the study. AJ, MS, and VC supervised the research and interpreted the data. KF, HB, LB, and AJ drafted the manuscript. AJ, RA, and JS revised the manuscript. All the authors discussed the results and commented on the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by Edital MS/CNPq/FAPITEC/SE/SES N° 06/2007 – PPSUS (Grant number 19.203-00775/2007-3), CAPES-Basic Parasitology (032/2010), and EDITAL CAPES/FAPITEC N° 11/2016 – PROEF (Grant number 88881.157436/2017-01). AJ is a scientist supported by the Brazilian Research and Technology Council (CNPq).
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors thank the population of Ilha das Flores and the patients from University Hospital, Sergipe, Brazil, who accepted to participate in this study. This study would not have been possible without the participation of the healthcare agents Jorge Feitosa and Edmilson Nicolaus dos Santos and the medical students Cybele Santos, Debora Barreto, Cinthia Andrade, and Karina Pessoa who helped the research team in the clinical, epidemiological, and parasitological survey.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.779534/full#supplementary-material
Supplementary Figure 1 | CD209 rs2287886 SNP is associated with differences IL-17 in schistosomiasis patients. Levels of IL-17 in 37 schistosomiasis patients stratified by CD209 rs2287886 genotype. *p < 0.05 (Mann-Whitney test).
Supplementary Figure 2 | Cytokine levels in no fibrosis and advanced fibrosis schistosomiasis patients by different genotypes to IL17A rs2275913, CD209 rs2287886, CD209 rs4804803 and IL10 rs1800871 SNP. Comparisons between the levels of IL-17 and IL-10 cytokines in supernatants from PBMC cultures not stimulated (Medium) or stimulated with adult worm (SWAP) and egg (SEA) S. mansoni antigens from schistosomiasis patient stratified by groups of advanced fibrosis (n=9) and no fibrosis (n=28) and by genotypes to CD209 rs2287886 (A), CD209 rs4804803 (B), IL-17A rs2275913 (C), and IL-10 rs1800871(D). *p < 0.05 (Mann-Whitney test).
References
1. World Health Organization. Schistosomiasis (Bilharzia). (2021). Available at: https://www.who.int/health-topics/schistosomiasis#tab=tab_1 (Accessed August 13, 2021).
2. Ponichtera HE, Shainheit MG, Liu BC, Raychowdhury R, Larkin BM, Russo JM, et al. CD209a Expression on Dendritic Cells Is Critical for the Development of Pathogenic Th17 Cell Responses in Murine Schistosomiasis. J Immunol (2014) 192:4655–65. doi: 10.4049/jimmunol.1400121
3. Gunda DW, Kilonzo SB, Manyiri PM, Peck RN, Mazigo HD. Morbidity and Mortality Due to Schistosoma Mansoni Related Periportal Fibrosis : Could Early Diagnosis of Varices Improve the Outcome Following Available Treatment Modalities in Sub Saharan Africa ? A Scoping Review. Trop Med Infect Dis (2020) 5:20. doi: 10.3390/tropicalmed5010020
4. Bergquist NR. Schistosomiasis Consortium for Operational Research and Evaluation : Mission Accomplished. Am J Trop Med Hyg (2020) 103:1–4. doi: 10.4269/ajtmh.19-0838
5. de Jesus AR, Magalhães A, Miranda DG, Miranda RG, Araújo MI, de Jesus AA, et al. Association of Type 2 Cytokines With Hepatic Fibrosis in Human Schistosoma Mansoni Infection. Infect Immun (2004) 72:3391–7. doi: 10.1128/IAI.72.6.3391-3397.2004
6. Kamdem SD, Moyou-Somo R, Brombacher F, Nono JK. Host Regulators of Liver Fibrosis During Human Schistosomiasis. Front Immunol (2018) 9:2781. doi: 10.3389/fimmu.2018.02781
7. Zheng B, Zhang J, Chen H, Nie H, Miller H. Gong Q. T Lymphocyte-Mediated Liver Immunopathology of Schistosomiasis. Front Immunol (2020) 11:61. doi: 10.3389/fimmu.2020.00061
8. Schwartz C, Fallon PG. Schistosoma “Eggs-Iting” the Host: Granuloma Formation and Egg Excretion. Front Immunol (2018) 9:2492. doi: 10.3389/fimmu.2018.02492
9. Kalantari P, Bunnell SC, Stadecker MJ. The C-Type Lectin Receptor-Driven , Pathology in Schistosomiasis: Not All Immune Responses to Helminth Parasites Are Th2 Dominated. Front Immunol (2019) 10:26. doi: 10.3389/fimmu.2019.00026
10. Castellucci LC, Jamieson SE, Miller EN, Menezes E, Oliveira J, Magalhães A, et al. CXCR1 and SLC11A1polymorphisms Affect Susceptibility to Cutaneous Leishmaniasis in Brazil: A Case-Control and Family-Based Study. BMC Med Genet (2010) 11:10. doi: 10.1186/1471-2350-11-10
11. Castellucci L, Menezes E, Oliveira J, Magalhaes A, Guimaraes LH, Lessa M, et al. IL6 -174 G/C Promoter Polymorphism Influences Susceptibility to Mucosal But Not Localized Cutaneous Leishmaniasis in Brazil. J Infect Dis (2006) 194:519–27. doi: 10.1086/505504
12. Morenikeji OB, Metelski JL, Hawkes ME, Capria AL, Seamans BN, Falade CO, et al. CD209 and Not CD28 or STAT6 Polymorphism Mediates Clinical Malaria and Parasitemia Among Children From Nigeria. Microorganisms (2020) 8:158. doi: 10.3390/microorganisms8020158
13. Mezger M, Steffens M, Semmler C, Arlt EM, Zimmer M, Kristjanson GI, et al. Investigation of Promoter Variations in Dendritic Cell-Specific ICAM3-Grabbing Non-Integrin (DC-SIGN) (CD209) and Their Relevance for Human Cytomegalovirus Reactivation and Disease After Allogeneic Stem-Cell Transplantation. Clin Microbiol Infect (2008) 14:228–34. doi: 10.1111/j.1469-0691.2007.01902.x
14. Silva PCV, Gomes AV, De Souza TKG, Coêlho MRCD, Cahu GGDOM, Muniz MTC, et al. Association of SNP (-G1082A) IL-10 With Increase in Severity of Periportal Fibrosis in Schistosomiasis, in the Northeast of Brazil. Genet Test Mol Biomarkers (2014) 18:646–52. doi: 10.1089/gtmb.2014.0098
15. Li Y-P, Wang M-Q, Liu C-R, Deng H-L, Wu Y, Dang S-S, et al. Polymorphisms in the DC-SIGN Gene and Their Association With the Severity of Hand, Foot, and Mouth Disease Caused by Enterovirus 71. Arch Virol (2021) 166:1133–40. doi: 10.1007/s00705-021-04991-6
16. Sachidanandam R, Weissman D, Schmidt SC, Kakol JM, Stein LD, Marth G, et al. A Map of Human Genome Sequence Variation Containing 1.42 Million Single Nucleotide Polymorphisms. Nature (2001) 409:928–33. doi: 10.1038/35057149
17. Mbow M, Larkin BM, Meurs L, Wammes LJ, de Jong SE, Labuda LA, et al. T-Helper 17 Cells Are Associated With Pathology in Human Schistosomiasis. J Infect Dis (2013) 207:186–95. doi: 10.1093/infdis/jis654
18. Mesquita GV, Brandt CT, Dos Santos SM, Maricevich JPBR. Avaliação De Pacientes Esquistossomóticos Quanto Ao Conteúdo Mineral Ósseo, Nível De Magnésio Corporal Total E Atividade Da Enzima Lecitina: Colesterol Acetiltransferase Em Pós-Esplenectomizados Com Auto-Implante Esplênico. Rev Bras Ortop (2004) 38(1/2):49–58.
19. Fernandes DA, Chagas ACP, de Jesus AR, França AVC, de Lima FS, da Silva ÂM, et al. Aspectos Ultrassonográficos Associados À Morbidade De Formas Clínicas Crônicas De Esquistossomose Mansônica, Utilizando-Se Protocolo Proposto Pela Organização Mundial Da Saúde. Radiol Bras (2013) 46:1–6. doi: 10.1590/S0100-39842013000100006
20. Barreto na VMS, Alecrim VM, De Medeiros TB, Domingues ALC, Lopes EP, Martins JRM, et al. New Index for the Diagnosis of Liver Fibrosis in Schistosomiasis Mansoni. Arq Gastroenterol (2017) 54:51–6. doi: 10.1590/s0004-2803.2017v54n1-10
21. Davi M. Schistosomiasis Mansoni: Ultrasound-Evaluated Hepatic Fibrosis and Serum Concentrations of Hyaluronic Acid. Ann Trop Med Parasitol (2011) 105:233–9. doi: 10.1179/136485911X12987676649629
22. Akpata R, Neumayr A, Holtfreter MC, Krantz I, Singh DD, Mota R, et al. The WHO Ultrasonography Protocol for Assessing Morbidity Due to Schistosoma Haematobium . Acceptance and Evolution Over 14 Years . Systematic Review. Parasitol Res (2015) 114:1279–89. doi: 10.1007/s00436-015-4389-z
23. Richter J, Botelho MC, Holtfreter MC, Akpata R, El Scheich T, Neumayr A, et al. Ultrasound Assessment of Schistosomiasis. Z Gastroenterol (2016) 54:653–60. doi: 10.1055/s-0042-107359
24. de Jesus AR, Silva A, Santana LB, Magalhães A, de Jesus AA, de Almeida RP, et al. Clinical and Immunologic Evaluation of 31 Patients With Acute Schistosomiasis Mansoni. J Infect Dis (2002) 185:98–105. doi: 10.1086/324668
25. de Oliveira FA, Barreto AS, Bomfim LGS, Leite TRS, dos Santos PL, de Almeida RP, et al. Soluble CD40 Ligand in Sera of Subjects Exposed to Leishmania Infantum Infection Reduces the Parasite Load in Macrophages. PloS One (2015) 10:e0141265. doi: 10.1371/journal.pone.0141265
26. Santos MB, De Oliveira DT, Cazzaniga RA, Varj CS, Santos PL, Duthie MS, et al. Distinct Roles of Th17 and Th1 Cells in Inflammatory Responses Associated With the Presentation of Paucibacillary Leprosy and Leprosy Reactions. Scand J Immunol (2017) 86:40–9. doi: 10.1111/sji.12558
27. Cheever AW, Williams ME, Wynn TA, Finkelman FD, Seder RA, Cox TM, et al. Anti-IL-4 Treatment of Schistosoma Mansoni-Infected Mice Inhibits Development of T Cells and Non-B, Non-T Cells Expressing Th2 Cytokines While Decreasing Egg-Induced Hepatic Fibrosis. J Immunol (1994) 153:753–9.
28. Metwali A, Elliott D, Blum AM, Weinstock JV. Granuloma Eosinophils Enhance IL-5 Production by Lymphocytes From Mice Infected With Schistosoma Mansoni. J Immunol (1993) 151:7048–56.
29. Chiaramonte MG, Mentink-Kane M, Jacobson BA, Cheever AW, Whitters MJ, Goad MEP, et al. Regulation and Function of the Interleukin 13 Receptor Alpha 2 During a T Helper Cell Type 2-Dominant Immune Response. J Exp Med (2003) 197:687–701. doi: 10.1084/jem.20020903
30. Smith PM, Shainheit MG, Bazzone LE, Rutitzky LI, Stadecker MJ. Genetic Control of Severe Egg-Induced Immunopathology and IL-17 Production in Murine Schistosomiasis. J Immunol (2009) 183:3317–23. doi: 10.4049/jimmunol.0901504
31. Shainheit MG, Lasocki KW, Finger E, Larkin BM, Patrick M, Sharpe AH, et al. The Pathogenic Th17 Cell Response to Major Schistosome Egg Antigen Is Sequentially Dependent on IL-23 and IL-1 β. J Immunol (2011) 187:5328–35. doi: 10.4049/jimmunol.1101445
32. Kalantari P, Morales Y, Miller EA, Jaramillo LD, Holly E, Wuethrich MA, et al. CD209a Synergizes With Dectin-2 and Mincle to Drive Severe Th17 Cell-Mediated Schistosome Egg-Induced Immunopathology. Cell Rep (2018) 22:1288–300. doi: 10.1016/j.celrep.2018.01.001.CD209a
33. Soloviova K, Fox EC, Dalton JP, Caffrey CR, Davies SJ. A Secreted Schistosome Cathepsin B1 Cysteine Protease and Acute Schistosome Infection Induce a Transient T Helper 17 Response. PloS Negl Trop Dis (2019) 13:1–22. doi: 10.1371/journal.pntd.0007070
34. Babu S, Bhat SQ, Kumar NP, Lipira AB, Kumar S, Karthik C, et al. Filarial Lymphedema Is Characterized by Antigen- Specific Th1 and Th17 Proinflammatory Responses and a Lack of Regulatory T Cells. PloS Negl Trop Dis (2009) 3(4):e420. doi: 10.1371/journal.pntd.0000420
35. Zhan T, Zhang T, Wang Y, Wang X, Lin C, Ma H, et al. Dynamics of Th9 Cells and Their Potential Role in Immunopathogenesis of Murine Schistosomiasis. Parasit Vectors (2017) 10:305. doi: 10.1186/s13071-017-2242-1
36. Qin S, Lu D, Guo X, Luo W, Hu B, Huang X. A Deleterious Role for Th9 / IL-9 in Hepatic Fibrogenesis. Sci Rep (2016) 6(18694):1–12. doi: 10.1038/srep18694
37. Jones CP, Gregory LG, Causton B, Campbell GA, Lloyd CM. Activin a and TGF- β Promote T H 9 Cell – Mediated Pulmonary Allergic Pathology. J Allergy Clin Immunol (2012) 129:1000–10. doi: 10.1016/j.jaci.2011.12.965
38. Kearley J, Erjefalt JS, Andersson C, Benjamin E, Carla P, Robichaud A, et al. IL-9 Governs Allergen-Induced Mast Cell Numbers in the Lung and Chronic Remodeling of the Airways. Am J Respir Crit Care Med (2011) 183:865–75. doi: 10.1164/rccm.200909-1462OC.IL-9
39. Barreto AVMS, Lacerda GAN, Figueiredo AL de C, Diniz GTN, Gomes ECS, Domingues ALC, et al. Evaluation of Serum Levels of IL-9 and IL-17 in Human Schistosoma Mansoni Infection and Their Relationship With Periportal Fibrosis. Immunobiology (2016) 221:1351–4. doi: 10.1016/j.imbio.2016.07.014
40. Okano M, Satoskar AR, Nishizaki K, Abe M, Harn DA. Induction of Th2 Responses and Ige Is Largely Due to Carbohydrates Functioning as Adjuvants on Schistosoma Mansoni Egg Antigens. J Immunol (1999) 163:6712–7.
41. Silveira AMS, Gazzinelli G, Alves-Oliveira LF, Bethony J, Gazzinelli A, Carvalho-Queiroz C, et al. Human Schistosomiasis Mansoni: Intensity of Infection Differentially Affects the Production of Interleukin-10, Interferon-Gamma and Interleukin-13 by Soluble Egg Antigen or Adult Worm Antigen Stimulated Cultures. Trans R Soc Trop Med Hyg (2004) 98:514–9. doi: 10.1016/j.trstmh.2003.11.009
42. Il I, Mentink-kane MM, Cheever AW, Wilson MS, Madala SK, Beers M, et al. Accelerated and Progressive and Lethal Liver Fibrosis in Mice That Lack Interleukin (IL)-10, IL-12p40, and IL-13R α2. Gastroenterology (2011) 141:2200–9. doi: 10.1053/j.gastro.2011.08.008
43. Wynn TA, Cheever AW, Williams ME, Caspar P, Kühn R, Müller W, et al. IL-10 Regulates Liver Pathology in Acute Murine Schistosomiasis Mansoni But Is Not Required for Immune Down-Modulation of Chronic Disease. J Immunol (1998) 160:4473–80.
44. Wilson MS, Mentink-Kane MM, Pesce JT, Ramalingam TR, Wynn TA. Immunopathology of Schistosomiasis. Immunol Cell Biol (2007) 85:148–54. doi: 10.1038/sj.icb.7100014
45. Gieseck RL III, Wilson MS, Wynn TA. Type 2 Immunity in Tissue Repair and Fibrosis. Nat Rev Immunol (2018) 18:62–76. doi: 10.1038/nri.2017.90
47. Possomato-Vieira JS, Khalil RA. Dendritic Cell Expression of the C-Type Lectin Receptor CD209a: A Novel Innate Parasite-Sensing Mechanism Inducing Th17 Cells That Drive Severe Immunopathology in Murine Schistosome Infection. Physiol Behav (2016) 176:139–48. doi: 10.1016/j.exppara.2015.04.006
48. Logan PGB, GD. CD209a Expression on Dendritic Cells Is Critical for the Development of Pathogenic Th17 Cell Responses in Murine Schistosomiasis. Bone (2014) 23:1–7. doi: 10.4049/jimmunol.1400121.CD209a
49. Chiaramonte MG, Cheever AW, Malley JD, Donaldson DD, Wynn TA. Studies of Murine Schistosomiasis Reveal Interleukin-13 Blockade as a Treatment for Established and Progressive Liver Fibrosis. Hepatology (2001) 34:273–82. doi: 10.1053/jhep.2001.26376
50. Shainheit MG, Saraceno R, Bazzone LE, Rutitzky LI, Stadecker MJ. Disruption of Interleukin-27 Signaling Results in Impaired Gamma Interferon Production But Does Not Significantly Affect Immunopathology in Murine Schistosome Infection. Infect Immun (2007) 75:3169–77. doi: 10.1128/IAI.01053-06
51. Adewusi OI, Nix NA, Lu X, Colley DG, Secor WE. Schistosoma Mansoni: Relationship of Tumor Necrosis Factor-a to Morbidity and Collagen Deposition in Chronic Experimental Infection. Exp Parasitol (1996) 84:115–23. doi: 10.1006/expr.1996.0097
52. Rutitzky LI, Lopes JR, Miguel J, Rutitzky LI, Lopes JR, Stadecker MJ. Severe CD4 T Cell-Mediated Immunopathology in Murine Schistosomiasis Is Dependent on IL-12p40 and Correlates With High Levels of IL-17. J Immunol (2021) 175:3920–6. doi: 10.4049/jimmunol.175.6.3920
53. Mwatha JK, Kimani G, Kamau T, Mbugua GG, Ouma JH, Mumo J, et al. High Levels of TNF, Soluble TNF Receptors, Soluble ICAM-1, and IFN-Gamma, But Low Levels of IL-5, Are Associated With Hepatosplenic Disease in Human Schistosomiasis Mansoni. J Immunol (1998) 160:1992–9.
Keywords: schistosomiasis, IL-9, IL-10, IL-17, hepatic fibrosis
Citation: Franco KGS, de Amorim FJR, Santos MA, Rollemberg CVV, de Oliveira FA, França AVC, Santos CNO, Magalhães LS, Cazzaniga RA, de Lima FS, Benevides L, Carregaro V, Silva JS, Brito HLdF, Fernandes DA, da Silva ÂM, de Almeida RP, Bezerra-Santos M and de Jesus AR (2021) Association of IL-9, IL-10, and IL-17 Cytokines With Hepatic Fibrosis in Human Schistosoma mansoni Infection. Front. Immunol. 12:779534. doi: 10.3389/fimmu.2021.779534
Received: 19 September 2021; Accepted: 15 November 2021;
Published: 14 December 2021.
Edited by:
Fabrizio Bruschi, University of Pisa, ItalyReviewed by:
Rashika El Ridi, Cairo University, EgyptVicente de Paulo Martins, University of Brasilia, Brazil
Copyright © 2021 Franco, de Amorim, Santos, Rollemberg, de Oliveira, França, Santos, Magalhães, Cazzaniga, de Lima, Benevides, Carregaro, Silva, Brito, Fernandes, da Silva, de Almeida, Bezerra-Santos and de Jesus. 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: Amélia Ribeiro de Jesus, ameliaribeirodejesus@gmail.com
†These authors have contributed equally to this work and share first authorship
‡These authors have contributed equally to this work and share last authorship