- Research Article
- Open Access
Gut microbiota signatures in Schistosoma japonicum infection-induced liver cirrhosis patients: a case–control study
Infectious Diseases of Poverty volume 10, Article number: 43 (2021)
Several studies have assessed the role of gut microbiota in various cirrhosis etiologies, however, none has done so in the context of Schistosoma japonicum infection in humans. We, therefore, sought to determine whether gut microbiota is associated with S. japonicum infection-induced liver cirrhosis.
From December 2017 to November 2019, 24 patients with S. japonicum infection-induced liver cirrhosis, as well as 25 age- and sex-matched controls from the Zhejiang Province, China, were enrolled. Fecal samples were collected and used for 16S rRNA gene sequencing (particularly, the hypervariable V4 region) using the Illumina MiSeq system. Wilcoxon Rank-Sum and PERMANOVA tests were used for analysis.
Eight hundred and seven operational taxonomic units (OTUs) were detected, of which, 491 were common between the two groups, whereas 123 and 193 were unique to the control and cirrhosis groups, respectively. Observed species, Chao, ACE, Shannon, Simpson, and Good’s coverage indexes, used for alpha diversity analysis, showed values of 173.4 ± 63.8, 197.7 ± 73.0, 196.3 ± 68.9, 2.96 ± 0.57, 0.13 ± 0.09, and 1.00 ± 0.00, respectively, in the control group and 154.0 ± 68.1, 178.6 ± 75.1, 179.9 ± 72.4, 2.68 ± 0.76, 0.19 ± 0.18, and 1.00 ± 0.00, respectively, in the cirrhosis group, with no significant differences observed between the groups. Beta diversity was evaluated by weighted UniFrac distances, with values of 0.40 ± 0.13 and 0.40 ± 0.11 in the control and cirrhosis groups, respectively (P > 0.05). PCA data also confirmed this similarity (P > 0.05). Meanwhile, the relative abundance of species belonging to the Bacilli class was higher in cirrhosis patients [median: 2.74%, interquartile range (IQR): 0.18–7.81%] than healthy individuals (median: 0.15%, IQR: 0.47–0.73%; P < 0.01), and that of Lactobacillales order was also higher in cirrhosis patients (median: 2.73%, IQR: 0.16–7.80%) than in healthy individuals (median: 0.12%, IQR: 0.03–0.70%; P < 0.05).
Cumulatively, our results suggest that the gut microbiota of S. japonicum infection-induced liver cirrhosis patients is similar to that of healthy individuals, indicating that bacterial taxa cannot be used as non-invasive biomarkers for S. japonicum infection-induced liver cirrhosis.
Schistosomiasis is one of the most prevalent parasitic diseases in the world . According to the World Health Organization (2020), an estimated 240 million people are affected by this disease, together with an additional 700 million at risk of infection . Worryingly, schistosomiasis alone is associated with a high disease burden, of as much as 3.3 million disability-adjusted life years, accounting for more than 10% of the global burden of neglected tropical diseases [3, 4]. Five main etiologic agents have been described for schistosomiasis, all of which are trematode parasites of Schistosoma genus, including S. japonicum, S. mansoni, S. haematobium, S. intercalatum, and S. mekongi . However, only S. japonicum exists in China, and is predominant in 12 provinces south of the Yangtze River . In the 1950s, an estimated 11.6 million people were infected, with more than 100 million at risk of infection . From the mid-1950s, the number of infected cases was reduced by over 99%; however, schistosomiasis remains endemic in 140 counties in China .
One of the serious complications of S. japonicum infection is the development of liver cirrhosis . Most egg-laying adult Schistosoma trematodes have liver tropism, and the local (hepatic) host immune responses to the deposited eggs result in granuloma formation, which ultimately evolves to progressive liver fibrosis and portal hypertension . Both in vivo and clinical studies focusing on S. japonicum infections revealed that hepatic stellate cells are the key drivers of hepatic fibrosis [8, 10], via collagen production and deposition. Importantly, liver fibrosis presentations have the potential to evolve into liver cirrhosis or liver cancer . Several studies have assessed the role of the gut microbiota in cirrhosis of different etiologies [12,13,14,15,16] with microbiota dysbiosis having been reported in cirrhosis patients as follows: (i) autochthonous/non-autochthonous taxa ratio reduction, (ii) Firmicutes/Bacteroidetes ratio inversion, or (iii) increased prevalence of potentially pathogenic bacteria [17,18,19]. Interestingly, the liver-gut axis is involved in the physiopathology of cirrhosis; specifically, overgrown intestinal bacteria have the potential to translocate across the leaky gut barrier and reach the liver, leading to cirrhosis development or progression .
Recent studies using not only animal models, but also human clinical samples have demonstrated a clear relationship between gut microbiota alterations and Schistosoma spp. infection [21,22,23,24,25,26]. However, no study has explored the relationship between the gut microbiota and S. japonicum infection-induced liver cirrhosis. Accordingly, the potential role of gut microbiota alterations in disease onset and progression remains unclear.
Here, we hypothesized that an association between gut microbiota and S. japonicum infection-induced liver cirrhosis exists. To explore this hypothesis, a case-control study was conducted. The main objectives of this study were as follows: (1) to characterize the potential gut microbiota alterations associated with S. japonicum infection-induced liver cirrhosis; (2) to identify whether there are potential bacterial taxa that can be used as non-invasive biomarkers of liver cirrhosis secondary to schistosomiasis.
This study was conducted according to the principles expressed in the Declaration of Helsinki. All experiments were approved by the First Affiliated Hospital Clinical Research Ethics Committee, School of Medicine, Zhejiang University (ref: 2017411-1), and by the Jiaxing Wangdian People's Hospital Ethics Committee (ref: 2017002). All participants provided written informed consent.
The S. japonicum infection-induced liver cirrhosis patient inclusion criteria were as follows: (i) diagnosis of schistosomiasis performed using the Kato-Katz method to detect S. japonicum eggs in stool samples; (ii) diagnosis of liver cirrhosis according to clinical and biochemical data, as well as imaging examination [computed tomography (CT) or B-ultrasound] [27, 28]; (iii) age ≥ 18 years. Eligibility was defined by the application of the following exclusion criteria: (i) history of digestive disorders, including other causes of liver disease, irritable bowel syndrome, inflammatory bowel disease, or chronic diarrhea; (ii) history of previous gastrointestinal surgery; (iii) history of probiotic, prebiotic, synbiotic, antibiotic, acid suppressor, metformin, or gastrointestinal motility targeting-drug use less than eight weeks before enrollment; (iv) history of bowel preparation, less than four weeks before enrollment; (v) diagnosis of major diseases that can affect the gut microbiota, such as malignant tumors, cardiovascular, respiratory, autoimmune, or allergic diseases, neurological disorders, renal insufficiency, diabetes, uncontrolled hypertension (blood pressure ≥ 150/90 mmHg), dyslipidemia, depression, mania, and bipolar affective disorder; (vi) hepatitis B virus, hepatitis C virus, or HIV infection.
Due to the lack of data on effect size, power analysis could not be performed. However, a recent similar study on S. japonicum infection in which each group consisted of 24 animals supports the appropriateness of our sample size . From December 2017 to November 2019, 24 patients with S. japonicum infection-induced liver cirrhosis were enrolled from Hangzhou and Jiaxing, northern cities of the Zhejiang Province, China, a schistosomiasis japonicum endemic area. Liver cirrhosis patients (either inpatients or outpatient) conforming to all inclusion criteria and not meeting any of the exclusion criteria were enrolled by specified physicians. All 24 patients included in this study were previously infected (neither currently infected nor recurrently infected). Eleven patients were treated with praziquantel in the acute infection period, whereas eight were not treated due to patient refusal of treatment for no or mild symptoms; for the remaining five patients, no information was available with respect to schistosomiasis treatment. In addition, 26 age- and sex-matched healthy individuals not meeting any of the exclusion criteria were enrolled from Hangzhou and Jiaxing. Since the stool specimen of one healthy individual was deemed unsatisfactory (stool mixed with blood), the healthy control group was ultimately composed of 25 individuals. The demographic information, laboratory tests, CT, and B-ultrasound results from all participants were collected from the hospital electronic medical records system. Additionally, questionnaires were filled by all subjects to obtain their sex, age, marriage status, occupation, education level, diet, eating behavior, exercise habits, and smoking or alcohol use status. Questionnaire data quality was controlled through special training of investigators on objectives, filling, and administration methods. Liver function was evaluated using the Child–Pugh classification system  and was based on five parameters, including ascites, hepatic encephalopathy, serum albumin, total bilirubin (TB), and prothrombin time (PT). Child–Pugh scores were calculated and defined as follows: grade A, 5–6 points; grade B, 7–9 points; grade C, ≥ 10 points.
Fecal sample collection
The participants were trained by a researcher to obtain a complete stool sample. Fresh stool samples (self-collected by patients and controls) were collected into anaerobic bags (Mitsubishi Gas Chemical, Tokyo, Japan). Samples were divided into five 200-mg aliquots, within the first 30 min after collection, and immediately stored at − 80 °C until further analysis.
Genomic DNA extraction
Total genomic DNA was extracted per sample using the MagPure Stool DNA KF kit B (Magen, Guangzhou, China) according to the manufacturer’s instructions. DNA was quantified using the Qubit Fluorometer (Qubit 2.0, Invitrogen, CA, USA) and the Qubit® dsDNA BR Assay kit (Invitrogen). DNA quality was assessed using gel electrophoresis (1% agarose gels). Negative controls containing elution buffer only were included for DNA extraction and quantification to monitor possible contamination.
Library construction and sequencing
The bacterial 16S rRNA gene hypervariable region V4 was amplified using degenerate polymerase chain reaction (PCR) with the primers 515Fw (5′-GTGCCAGCMGCCGCGGTAA-3′) and 806Rv (5′-GGACTACHVGGGTWTCTAAT-3′), which contained the Illumina adapter sequence, a pad, and a linker. Fifty-microliter reaction volumes, containing 30 ng of DNA templates together with the primers (10 μmol/L) and the PCR master mix, were used. Amplification was performed under the following conditions: 95 °C for 3 min; 30 cycles of 95 °C for 45 s, 56 °C for 45 s, and 72 °C for 45 s; a final extension step at 72 °C for 10 min. PCR products were purified using Agencourt AMPure XP beads (Beckman Coulter, CA, USA) according to the manufacturer’s instructions. Library quality was evaluated using the 2100 Bioanalyzer (HT DNA 1K Labchip 760517) instrument (Agilent Technologies, Santa Clara, CA, USA). Sequencing was performed using the Illumina HiSeq 2500 platform (BGI, Shenzhen, China) following the Illumina standard pipelines to obtain 2 × 250 bp paired-end reads.
Raw reads were filtered to remove adaptors and low-quality and/or ambiguous bases, and paired-end reads were merged using the Fast Length Adjustment of Short reads (FLASH, v1.2.11) software . Sequences were clustered into operational taxonomic units (OTUs) with a cutoff value of 97% using the UPARSE software (v7.0.1090) . Chimeric sequences were detected (compared with the Gold database) using UCHIME software (v4.2.40)  and removed from the analysis. Sequences of representative OTUs were then taxonomically classified using the Ribosomal Database Project Classifier (v.2.2) with a minimum confidence threshold of 0.6 and trained on the Greengenes database (v201305) using the QIIME pipeline (v1.8.0) . The USEARCH_global command  was used to quantify the abundance of OTUs in each sample.
Bacterial diversity analysis
Alpha and beta diversities were estimated using MOTHUR software (v1.31.2)  and the QIIME pipeline (v1.8.0) , respectively. Sample clustering was conducted using the QIIME pipeline (v1.8.0)  and the unweighted pairwise grouping method with averaging. Barplots and heatmaps of different classification levels were plotted with the R package (v3.4.1) “heatmap” and “gplots” tools, respectively. Venn diagrams and accumulation curves were obtained using the R package “VennDiagram” (v 3.1.1) tool. A principal component analysis (PCA) was performed using the R package “ade4” tool. Taxonomic characterization was performed using a linear discriminant analysis (LDA) with the calculation of the LDA effect size (LEfSe).
Continuous variables are presented as the mean ± standard deviation (SD) or median [interquartile range (IQR)], while categorical variables are presented as percentages (%). The Student’s t-test or Wilcoxon Rank-Sum test and the χ2 test were used to evaluate differences between two groups for continuous and categorical variables, respectively. Whenever applicable, the false discovery rate (FDR) was used to correct the calculated P-value. Statistical analyses were performed using SPSS software (v 22.0; SPSS, IL, USA). Statistical significance was given by P < 0.05, FDR < 0.05, or LDA scores > 2.
The study population included a total of 24 patients with S. japonicum infection-induced liver cirrhosis (12 males and 12 females), and 25 age- and sex-matched healthy individuals (12 males and 13 females). The average ages of cirrhosis and control groups were comparable, 82.7 ± 7.3 and 82.3 ± 7.1 years, respectively (Table 1); in fact, since this was a case-control study with matched controls, no differences in terms of age or sex were expected. Furthermore, no significant differences were observed concerning body mass index, or smoking and drinking habits between the two groups. In addition, the platelet counts, and serum levels of TB, direct bilirubin, and alanine aminotransferase were comparable between the groups. However, whereas the serum levels of aspartate aminotransferase (AST) and glutamyl transpeptidase (GGT), together with the PT were significantly increased in the S. japonicum infection-induced liver cirrhosis group, the total protein, and albumin serum levels were significantly decreased in the same group, in comparison with those in the control group (Table 1; at least P < 0.05). Regarding disease severity stratification, performed only for the cirrhosis group, 13 patients (54.2%) were classified as Child–Pugh A, seven (29.2%) as Child–Pugh B, and four (16.7%) as Child–Pugh C (Table 1).
Gut microbial alpha and beta diversity analysis between the two groups
The obtained accumulation curve reached a plateau, indicating that most of the gut microbial populations had been detected at the employed sequencing depth (Fig. 1). A total of 807 OTUs were obtained, of which 491 were equally detected in samples from cirrhosis individuals and healthy controls (Fig. 2). Notably, as represented in the Venn diagram, 123 and 193 OTUs were uniquely detected in healthy individuals and patients with S. japonicum infection-induced liver cirrhosis, respectively.
Concerning alpha diversity, we used six different metrics for analysis; specifically, the observed species, Chao, ACE, Shannon, Simpson, and Good’s coverage indexes were calculated. The resulting values were 173.4 ± 63.8, 197.7 ± 73.0, 196.3 ± 68.9, 2.96 ± 0.57, 0.13 ± 0.09, and 1.00 ± 0.00, respectively, in the control group and 154.0 ± 68.1, 178.6 ± 75.1, 179.9 ± 72.4, 2.68 ± 0.76, 0.19 ± 0.18, and 1.00 ± 0.00, respectively, in the cirrhosis group. Overall, we did not detect any significant changes when comparing healthy individuals to patients with S. japonicum infection-induced liver cirrhosis. However, a slight tendency toward a lower bacterial diversity was observed in the cirrhosis group (particularly for observed species, Chao, ACE, and Shannon indexes; Fig. 3). In addition, even considering the small number of cases showing disease severity in the cirrhosis group (13 Child–Pugh A, seven Child–Pugh B, and four Child–Pugh C cases), we further explored the variance in microbiota data by disease severity. The results had no statistical power likely due to the smaller sample size (particularly of the Child–Pugh C cohort); however, major trends could still be observed in Child–Pugh C patients (decreased Shannon index and increased Simpson index).
We further analyzed beta diversity (based on weighted UniFrac distances), to evaluate differences between samples, and again, no major differences were observed between control and cirrhosis groups (0.40 ± 0.13 vs 0.40 ± 0.11; Fig. 4a). PCA data (Fig. 4b) confirmed this similarity. Samples from each group showed heterogeneous distributions that were not species-specific, indicating that the structure of gut bacterial communities is similar between healthy individuals and patients with S. japonicum infection-induced liver cirrhosis (Fig. 4b).
Taxonomic composition of gut microbiota between the two groups
We evaluated bacterial relative abundance at both phylum and genus levels (Fig. 5a and b, respectively). Analysis of the most abundant bacterial phyla in the fecal samples collected from the two groups revealed a similar pattern; specifically, Bacteroidetes, Firmicutes, and Proteobacteria were the three most prevalent phyla in both groups, accounting for 97.1% and 94.4% of all bacterial populations in healthy controls and cirrhosis patients, respectively (Fig. 5a). However, different patterns were observed at the genus level; specifically, Bacteroides, Prevotella, and Faecalibacterium were the three most prevalent genera in the guts of healthy individuals, whereas in cirrhosis patients, the third most prevalent genus was Escherichia (Fig. 5b).
Additionally, to explore statistical relevance, we performed a LEfSe analysis, wherein P < 0.05 (Wilcoxon Rank-Sum test) and LDA scores > 2 were considered statistically significant (Figs. 6 and 7). Of note, we also considered FDR values in this analysis to minimize the number of false-positive results. Overall, we found no significant differences between the two groups at the phylum or genus levels, the most relevant categories in the context of taxonomic composition. Regarding the taxonomic composition of gut microbiota at other levels (class–order–family–species), only the Bacilli class and the Lactobacillales order were significantly increased in cirrhosis patients compared to those in healthy individuals (Fig. 8).
Liver cirrhosis is a common (end-stage) consequence of several liver diseases and a major cause of mortality worldwide . Previous studies on gut-liver axis have demonstrated that the gut microbiota (and bacterial translocation) plays an important role in the pathogenesis of cirrhosis and its complications [37, 38]. However, few previous studies have explored the relationship between gut microbiota and S. japonicum infection-induced liver cirrhosis. Hence, this study explored the possible association between the gut microbiota and S. japonicum infection-induced liver cirrhosis in humans.
In this study, patients with liver cirrhosis showed increased levels of serum PT, AST, and GGT and decreased total protein, and albumin; of note, liver cirrhosis led to all of these changes. However, the globulin serum levels, defined previously as a good biomarker of inflammation and the immune status , did not change in the context of S. japonicum infection-induced liver cirrhosis. This result suggests that S. japonicum infection-induced liver cirrhosis might have a limited impact on systemic immune functions.
Surprisingly, we found that the microbiota of S. japonicum infection-induced liver cirrhosis patients was as rich and diverse as that of matched healthy individuals. Overall, these results are not consistent with findings reported in the context of in vivo studies. For example, in a mouse model of S. japonicum ova-induced granulomas, a significantly higher beta diversity was detected versus that in controls . Similarly, an overall reduction in alpha and a significant increase in beta diversities were found in the guts of S. mansoni-infected mice compared to those in uninfected controls . In our study, compared to healthy individuals, Child–Pugh C patients exhibited a decreased trend in the Shannon index and increased trend in the Simpson index. As such, an increase in the sample size is essential in future investigations to improve the statistical relevance and determine potential differences among groups.
Furthermore, the taxonomic composition at the phylum level was similar between cirrhosis patients and healthy individuals (prevalence of Bacteroidetes, Firmicutes, and Proteobacteria), aligning with in vivo data on S. japonicum infections . Additionally, even at the genus level, no major changes were detected, indicating that in the context of our study, there is no gut dysbiosis in S. japonicum infection-induced liver cirrhosis patients. This result was unexpected and differed from that of other studies reporting gut dysbiosis in the context of human schistosomiasis (e.g., those caused by S. mansoni and S. haematobium) [21,22,23,24,25]. Additionally, in vivo studies have reported Schistosoma infection-related alterations in gut microbial composition. For example, S. japonicum infection in mice was found to be associated with significant changes in the gut prevalence of Firmicutes (relative decrease), Bacteroidetes, and Proteobacteria (relative increase) .
Moreover, no differences were observed between the two groups at the phylum or genus levels. However, when observed at other levels (class–order–family–species) differences were noted. Specifically, Bacilli (at the class level) and Lactobacillales (at the order level) were significantly increased in patients with S. japonicum infection-induced liver cirrhosis compared to levels in healthy individuals. Importantly, our results are consistent with those of a previous study, which might suggest that the increased prevalence of Lactobacillales in the gut microbiota is a hallmark of human liver cirrhosis . Additionally, a different study showed that Lactobacillales and Bacilli were the key contributors to gut dysbiosis in primary liver cancer and proposed that they are potential diagnostic markers . Taken together, these findings suggest a strong association between gut Bacilli and Lactobacillales and liver disease that warrants further exploration in different pathological contexts.
Overall, our data do not support the idea that S. japonicum infection-induced cirrhosis (in elderly individuals) is associated with significant alterations to gut microbiota profiles, opposing hypotheses proposed by others concerning schistosomiasis. A few reasons might justify our opposing results. For instance, studies have focused on different Schistosoma spp. with different organ tropisms; for example, S. mansoni and S. japonicum target the gut and liver, whereas S. haematobium targets the bladder and urogenital system [21,22,23, 25]. Importantly, the notion that S. haematobium infection is equally associated with human gut dysbiosis suggests that potential alterations in intestinal microbial communities in the context of schistosomiasis are likely a consequence of infection-induced systemic immune responses, not local events [21,22,23]. For example, mesenteric lymph node cells isolated from mice infected with S. mansoni (dysbiotic) respond differently than those isolated from controls, including higher interleukin (IL)-4, IL-10, and IL-17 secretion, together with lower interferon-γ production . Appropriately, aging is accompanied by a deterioration of immune responses, recently named immunosenescence or immune-aging . In our study, the enrolled patients were of advanced age (average of 82.7 ± 7.3 years), which might justify the absence of this phenotype; in other words, their potentially weaker immune responses to S. japonicum infection were not strong enough to impact their gut microbial communities. Of note, the serum globulin levels in liver cirrhosis patients did not change significantly, which could indirectly reflect the weaker effect of S. japonicum infection on systemic immunity in elderly individuals.
Nevertheless, our geographic target-region dictated the age range of our study population, which is a limitation of our study. Since 1995, due to the application of effective schistosomiasis control measures, there are almost no newly reported cases in Zhejiang . Therefore, patients with S. japonicum infection-induced liver cirrhosis are not prevalent and tend to be of advanced ages . Hence, we cannot exclude a possible age selection bias that reflects the reality of the population studied; specifically, among the cirrhosis patients, 17 and 7 individuals were 80 or 90 years of age or older, respectively, and none of the others was less than 60 years of age. We can, however, also hypothesize that our population has a survival bias. Living patients could be those with better outcomes due to stable gut microbiota profiles, whereas prior fatalities might have been a consequence of gut microbiota disorders. This is, however, only speculative. Therefore, further studies focusing on much larger and heterogeneous populations (different geographies, with a wider age range) are needed to fully elucidate whether S. japonicum infection-induced liver cirrhosis is accompanied by gut microbiota alterations. In addition, the effective schistosomiasis control in Zhejiang and advanced age of participants might affect our conclusion because it was almost impossible to include a cohort of age-matched S. japonicum acutely infected patients without liver cirrhosis.
An additional factor that should be considered in future studies is the timing of infection, particularly if the study is designed to identify biomarkers for the diagnosis and/or prognosis of cirrhosis. Indeed, the infection status might also determine the results. Human studies have primarily targeted adolescents  or children , and consequently, the early phase of infection. Additionally, the most commonly used murine model of schistosomiasis in microbiome studies can be used to evaluate acute infection, specifically, less than eight weeks post-challenge according to the model characterization [43,44,45]. The in vivo published data refer to stool samples collected at 28 days , 42 days , and 8 weeks  after infection. Hence, we cannot exclude the possibility that our conflicting findings are simply a consequence of our study context. Cirrhosis is the end-stage of chronic schistosomiasis that occurs significantly later than the early phase of infection. Both systemic immune responses and Schistosoma infection were found to be associated with gut microbiota alterations . However, in S. japonicum infection-induced liver cirrhosis, both are expected to decrease in magnitude, which might partially explain our findings.
Overall, the results show that the gut microbial taxonomic profiles and diversity in elderly S. japonicum infection-induced liver cirrhosis patients appear to be similar to those from age-matched healthy individuals. Therefore, in our study, patients with S. japonicum infection-induced liver cirrhosis present a “healthy” gut structure, which might partially explain the better prognosis of cirrhosis secondary to this infection, as compared to that with cirrhosis of different etiologies in China. Future studies with a much larger sample size and considering a wider age range are, however, needed to further explore the relationship between S. japonicum infection-induced liver cirrhosis and gut microbiota, and potentially disclose the possible mechanisms by which gut microbiota alterations influence cirrhosis and/or vice versa. Importantly, these studies should also evaluate immunity and inflammation to finally determine if the gut-liver axis is implicated in S. japonicum infection-induced liver cirrhosis development and progression. Such studies should be further complemented by in vivo works with the transplantation of human gut microbiota samples to ultimately establish (or exclude) a causal relationship between dysbiosis and cirrhosis secondary to schistosomiasis.
Availability of data and materials
The data supporting the conclusions of this article are included in the manuscript. The datasets analyzed during the present study are available from the corresponding author upon reasonable request.
Operational taxonomic units
Polymerase chain reaction
Principal component analysis
Linear discriminant analysis
LDA effect size
- SD :
False discovery rate
S. japonicum Infection-induced liver cirrhosis
Cao ZG, Zhao YE, Lee Willingham A, Wang TP. Towards the elimination of schistosomiasis japonica through control of the disease in domestic animals in the People’s Republic of China: a tale of over 60 years. Adv Parasitol. 2016;92:269–306.
World Health Organization. World Health Organisation Schistosomiasis (Bilharzia). World Health Organization. 2021. https://www.who.int/health-topics/schistosomiasis#tab=tab_1. Accessed 24 Jan 2021.
Steinmann P, Keiser J, Bos R, Tanner M, Utzinger J. Schistosomiasis and water resources development: systematic review, meta-analysis, and estimates of people at risk. Lancet Infect Dis. 2006;6:411–25.
Murray CJ, Vos T, Lozano R, Naghavi M, Flaxman AD, Michaud C, et al. Disability-adjusted life years (DALYs) for 291 diseases and injuries in 21 regions, 1990–2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2197–223.
Chuah C, Gobert GN, Latif B, Heo CC, Leow CY. Schistosomiasis in Malaysia: a review. Acta Trop. 2019;190:137–43.
Zhang SQ, Sun CS, Wang M, Lin DD, Zhou XN, Wang TP. Epidemiological features and effectiveness of schistosomiasis control programme in lake and marshland region in the People’s Republic of China. Adv Parasitol. 2016;92:39–71.
Li ZJ, Ge J, Dai JR, Wen LY, Lin DD, Madsen H, et al. Biology and control of snail intermediate host of Schistosoma japonicum in the People’s Republic of China. Adv Parasitol. 2016;92:197–236.
Carson JP, Ramm GA, Robinson MW, McManus DP, Gobert GN. Schistosome-induced fibrotic disease: the role of hepatic stellate cells. Trends Parasitol. 2018;34:524–40.
Anthony B, Allen JT, Li YS, McManus DP. Hepatic stellate cells and parasite-induced liver fibrosis. Parasit Vectors. 2010;3:60.
Seki E, Brenner DA. Recent advancement of molecular mechanisms of liver fibrosis. J Hepatobiliary Pancreat Sci. 2015;22:512–8.
Kamdem SD, Moyou-Somo R, Brombacher F, Nono JK. Host regulators of liver fibrosis during human schistosomiasis. Front Immunol. 2018;9:2781.
Wei X, Jiang S, Chen Y, Zhao X, Li H, Lin W, et al. Cirrhosis related functionality characteristic of the fecal microbiota as revealed by a metaproteomic approach. BMC Gastroenterol. 2016;16:121.
Schwimmer JB, Johnson JS, Angeles JE, Behling C, Belt PH, Borecki I, et al. Microbiome signatures associated with steatohepatitis and moderate to severe fibrosis in children with nonalcoholic fatty liver disease. Gastroenterology. 2019;157:1109–22.
Bajaj JS, Thacker LR, Fagan A, White MB, Gavis EA, Hylemon PB, et al. Gut microbial RNA and DNA analysis predicts hospitalizations in cirrhosis. JCI Insight. 2018;3:e98019.
Heidrich B, Vital M, Plumeier I, Doscher N, Kahl S, Kirschner J, et al. Intestinal microbiota in patients with chronic hepatitis C with and without cirrhosis compared with healthy controls. Liver Int. 2018;38:50–8.
Loomba R, Seguritan V, Li W, Long T, Klitgord N, Bhatt A, et al. Gut microbiome-based metagenomic signature for non-invasive detection of advanced fibrosis in human nonalcoholic fatty liver disease. Cell Metab. 2017;25(1054–62):e5.
Bajaj JS, Heuman DM, Hylemon PB, Sanyal AJ, White MB, Monteith P, et al. Altered profile of human gut microbiome is associated with cirrhosis and its complications. J Hepatol. 2014;60:940–7.
Bajaj JS, Betrapally NS, Hylemon PB, Heuman DM, Daita K, White MB, et al. Salivary microbiota reflects changes in gut microbiota in cirrhosis with hepatic encephalopathy. Hepatology. 2015;62:1260–71.
Ling Z, Liu X, Cheng Y, Jiang X, Jiang H, Wang Y, et al. Decreased diversity of the oral microbiota of patients with hepatitis b virus-induced chronic liver disease: a pilot project. Sci Rep. 2015;5:17098.
Oikonomou T, Papatheodoridis GV, Samarkos M, Goulis I, Cholongitas E. Clinical impact of microbiome in patients with decompensated cirrhosis. World J Gastroenterol. 2018;24:3813–20.
Ajibola O, Rowan AD, Ogedengbe CO, Mshelia MB, Cabral DJ, Eze AA, et al. Urogenital schistosomiasis is associated with signatures of microbiome dysbiosis in Nigerian adolescents. Sci Rep. 2019;9:829.
Jenkins TP, Peachey LE, Ajami NJ, MacDonald AS, Hsieh MH, Brindley PJ, et al. Schistosoma mansoni infection is associated with quantitative and qualitative modifications of the mammalian intestinal microbiota. Sci Rep. 2018;8:12072.
Floudas A, Aviello G, Schwartz C, Jeffery IB, O’Toole PW, Fallon PG. Schistosoma mansoni worm infection regulates the intestinal microbiota and susceptibility to colitis. Infect Immun. 2019. https://doi.org/10.1128/IAI.00275-19.
Zhao Y, Yang S, Li B, Li W, Wang J, Chen Z, et al. Alterations of the mice gut microbiome via Schistosoma japonicum ova-induced granuloma. Front Microbiol. 2019;10:352.
Kay GL, Millard A, Sergeant MJ, Midzi N, Gwisai R, Mduluza T, et al. Differences in the faecal microbiome in Schistosoma haematobium infected children vs. uninfected children. PLoS Negl Trop Dis. 2015;9:e0003861.
Schneeberger PHH, Coulibaly JT, Panic G, Daubenberger C, Gueuning M, Frey JE, et al. Investigations on the interplays between Schistosoma mansoni, praziquantel and the gut microbiome. Parasit Vectors. 2018;11:168.
Wu W, Feng A, Huang Y. Research and control of advanced schistosomiasis japonica in China. Parasitol Res. 2015;114:17–27.
Song LG, Wu XY, Sacko M, Wu ZD. History of schistosomiasis epidemiology, current status, and challenges in China: on the road to schistosomiasis elimination. Parasitol Res. 2016;115:4071–81.
Cholongitas E, Papatheodoridis GV, Vangeli M, Terreni N, Patch D, Burroughs AK. Systematic review: the model for end-stage liver disease–should it replace Child-Pugh’s classification for assessing prognosis in cirrhosis? Aliment Pharmacol Ther. 2005;22:1079–89.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.
Lozano R, Naghavi M, Foreman K, Lim S, Shibuya K, Aboyans V, et al. Global and regional mortality from 235 causes of death for 20 age groups in 1990 and 2010: a systematic analysis for the Global Burden of Disease Study 2010. Lancet. 2012;380:2095–128.
Arab JP, Martin-Mateos RM, Shah VH. Gut-liver axis, cirrhosis and portal hypertension: the chicken and the egg. Hepatol Int. 2018;12:24–33.
Acharya C, Bajaj JS. Gut Microbiota and complications of liver disease. Gastroenterol Clin North Am. 2017;46:155–69.
Zhou Q, Cao H, Xu Z, Lan R, Chen X, Wang D, et al. Baseline serum globulin as a predictor of the recurrence of lone atrial fibrillation after radiofrequency catheter ablation. Anatol J Cardiol. 2017;17:381–5.
Bajaj JS, Ahluwalia V, Steinberg JL, Hobgood S, Boling PA, Godschalk M, et al. Elderly patients have an altered gut-brain axis regardless of the presence of cirrhosis. Sci Rep. 2016;6:38481.
Zhang L, Wu YN, Chen T, Ren CH, Li X, Liu GX. Relationship between intestinal microbial dysbiosis and primary liver cancer. Hepatobiliary Pancreat Dis Int. 2019;18:149–57.
Walford RL. The immunologic theory of aging. Gerontologist. 1964;4:195–7.
Ferrari TC, Moreira PR. Neuroschistosomiasis: clinical symptoms and pathogenesis. Lancet Neurol. 2011;10:853–64.
Seki T, Kumagai T, Kwansa-Bentum B, Furushima-Shimogawara R, Anyan WK, Miyazawa Y, et al. Interleukin-4 (IL-4) and IL-13 suppress excessive neutrophil infiltration and hepatocyte damage during acute murine schistosomiasis japonica. Infect Immun. 2012;80:159–68.
Chen L, Chen Q, Hou W, He L. High-throughput dynamic analysis of differentially expressed genes in splenic dendritic cells from mice infected with Schistosoma japonicum. Immunol Lett. 2017;184:15–22.
We would like to thank Editage (www.editage.cn) for English language editing and BGI-Shenzhen for genome-sequencing services.
This study was supported by the National Key Research and Development Program of China (2018YFC2000500, 2018YFC2000301), the Key Disciplines Construction Plan of the Zhejiang Province for Traditional Chinese Medicine (2017-XK-A31), and the Chinese Medicine Research Program of the Zhejiang Province (2020ZB148).
Ethics approval and consent to participate
The study was approved by the First Affiliated Hospital Clinical Research Ethics Committee, School of Medicine, Zhejiang University (ref: 2017411-1), and by the Jiaxing Wangdian People’s Hospital Ethics Committee (ref: 2017002). All participants provided written informed consent in Chinese. Participation was entirely voluntary and only participants that provided the signed informed consent forms were enrolled in the study. The rights of these subjects were protected.
Consent for publication
The authors declare that they have no competing interests.
About this article
Cite this article
Gui, QF., Jin, HL., Zhu, F. et al. Gut microbiota signatures in Schistosoma japonicum infection-induced liver cirrhosis patients: a case–control study. Infect Dis Poverty 10, 43 (2021). https://doi.org/10.1186/s40249-021-00821-8
- Gut microbiota
- 16S rRNA
- Schistosoma japonicum