- Research
- Open access
- Published:
Metabolomics analysis of patients with Schistosoma japonicum infection based on UPLC-MS method
Parasites & Vectors volume 17, Article number: 350 (2024)
Abstract
Background
Schistosomiasis is still one of the most serious parasitic diseases. Evidence showed that the metabolite profile in serum can potentially act as a marker for parasitic disease diagnosis and evaluate disease progression and prognosis. However, the serum metabolome in patients with Schistosoma japonicum infection is not well defined. In this study, we investigated the metabolite profiles of patients with chronic and with advanced S. japonicum infection.
Methods
The sera of 33 chronic S. japonicum patients, 15 patients with advanced schistosomiasis and 17 healthy volunteers were collected. Samples were extracted for metabolites and analyzed with ultra-performance liquid chromatography-mass spectrometry (UPLC-MS).
Results
We observed significant differences in metabolite profiles in positive and negative ion modes between patients with advanced and chronic S. japonicum infection. In patients with chronic S. japonicum infection, 199 metabolites were significantly upregulated while 207 metabolites were downregulated in advanced infection. These differential metabolites were mainly concentrated in steroid hormone biosynthesis, cholesterol metabolism and bile secretion pathways. We also found that certain bile acid levels were significantly upregulated in the progression from chronic to advanced S. japonicum infection. In receiver operator characteristic (ROC) analysis, we identified three metabolites with area under the curve (AUC) > 0.8, including glycocholic (GCA), glycochenodeoxycholate (GCDCA) and taurochenodeoxycholic acid (TCDCA) concentrated in cholesterol metabolism, biliary secretion and primary bile acid biosynthesis.
Conclusions
This study provides evidence that GCA, GCDCA and TCDCA can potentially act as novel metabolite biomarkers to distinguish patients in different stages of S. japonicum infection. This study will contribute to the understanding of the metabolite mechanisms of the transition from chronic to advanced S. japonicum infection, although more studies are needed to validate this potential role and explore the underlying mechanisms.
Graphical Abstract
Background
Schistosomiasis is one of the most serious parasitic diseases, second only to malaria in the global list of neglected tropical diseases, and one of the most devastating tropical parasitic diseases [1]. It is estimated that 240 million people are currently infected worldwide, and > 700 million people in 76 countries and territories are affected, mainly in tropical and subtropical regions [2,3,4,5,6]. Schistosomiasis caused by Schistosoma japonicum is much more severe than the disease caused by other Schistosoma species [7]. After eggs have been deposited in a patient’s liver, intestine and spleen, they can persist for a long time, leading to chronic damage and liver fibrosis, a stage known as chronic S. japonicum infection. Patients diagnosed with chronic S. japonicum infection usually are slightly symptomatic or even asymptomatic. If patients have recurrent infections, even with treatment, after 5–10 years they will develop advanced S. japonicum infection, resulting in cirrhotic ascites, portal hypertension and a marked reduction in quality of life [8].
Currently, some advances have been achieved in diagnosing S. japonicum infection. Nonetheless, S. japonicum infection diagnosis relies on traditional methods such as the Kato-Katz (KK) technique, miracidial hatching method and antigen or antibody-based detection [9,10,11]. These conventional diagnostic strategies show low sensitivity and specificity and are closely related to sample collection [10]. Nucleic acid-based techniques, including PCR-based methods, cell-free DNA-based detection and CRISPR-Cas12/13-mediated detection, are characterized by high sensitivity and specificity but are limited because they require expensive equipment, reagents and experienced personnel [9, 12]. Thus, there is an urgent need to develop new diagnostic methods [13]. Recent studies showed significant changes in serum metabolites in different stages of S. japonicum infection in mice, and serum metabolites have potential to act as an indicator for disease burden and assist early diagnosis [14, 15]. However, the serum metabolome has not been explored in patients infected with S. japonicum. Animal studies indicated that further studies on metabolite responses in patients with advanced and chronic S. japonicum infection may help us better understand the mechanisms of disease progression and shed light on tools for early diagnosis of S. japonicum infection.
Metabolomics involve the study of multivariate metabolite responses between the host internal environment and complex parasite infection, which can reflect host regulation of parasite growth and development as well as disease progression [16]. In this study, ultra-performance liquid chromatography-mass spectrometry (UPLC-MS)-based metabolomic method was applied to analyze serum metabolites changes in S. japonicum infection. Alterations in serum metabolites in different stages of S. japonicum infection were observed. This study provided evidence that the altered metabolites seem to have potential to be novel metabolite biomarkers to distinguish patients in different stages of S. japonicum infection. Also, this study will contribute to the understanding of the metabolite mechanisms of the transition from chronic to advanced S. japonicum infection.
Methods
Sample collection and clinical information
Patients with S. japonicum infection are diagnosed according to the National Standardized Diagnostic Criteria for S. japonicum infection (WS261-2006) of the Ministry of Health of China [17]. The diagnostic criteria for chronic schistosomiasis japonica were as follows: (i) a history of exposure to endemic water in endemic areas of S. japonicum; (ii) medical history of diagnosis of schistosomiasis based on Kato-Katz or indirect hemagglutination test and received praziquantel treatment; (iii) asymptomatic or intermittent abdominal pain, diarrhea or purulent and bloody stool; (iv) positive ovum test in the patient stool or positive serum immunological test; (v) B-ultrasound examination indicated fibrosis, without ascites, portal hypertension and splenomegaly, etc. [18]. Patients with advanced schistosomiasis should meet the following conditions: (i) history of repeated or prolonged exposure to water in endemic areas; (ii) history of treatment for S. japonicum infection; (iii) positive ovum stool test or positive serum immunological test; (iv) liver fibrosis; (v) one or more complications of splenomegaly, hypersplenism, portal hypertension or ascites. In this study, all patients (≥ 18 years old) with S. japonicum infection who visited the Xiangyue Hospital of Hunan Province from January 2022 to August 2022 and agreed to participate in this study were included. Patients with the following diseases were excluded from research: hybrid liver diseases; combined infection with other viruses, such as HAV, HBV, HCV or HDV; primary liver cancer or other malignant tumors; medium and severe fatty liver; severe liver dysfunction. Healthy volunteers were from the same area as those who received health examinations in the same hospital and agreed to participate in this study. According to the inclusion and exclusion criteria, 33 patients with chronic S. japonicum infection and 15 patients with advanced S. japonicum infection and 17 healthy volunteers were included. The serum and clinical information on the patients with S. japonicum infection in the hospital was collected. The related clinical information is shown in Table 1. Blood (5 ml) was collected early in the morning from all participants who fasted for 10 h and was coagulated at room temperature. Serum was obtained by centrifugation at 3000 rpm for 5 min and stored at − 80 °C. We used the same standard for all samples. This study was approved by the Institutional Review Board (IRB) of the Third Xiangya Hospital of Central South University (21149), and written informed consent was obtained from all participants.
Serum sample analysis using UPLC-MS
Chemicals and reagents
Methanol (A454-4) and acetonitrile (A998-4) were all analytical grade (Thermo Fisher Scientific, USA); ammonia formate (17843-250G, Honeywell Fluka, USA) and formic acid (50144-50Â ml, DIMKA, USA) were used. Water was supplied by ultra-pure water.
Sample preparation
Serum samples were stored at − 80 °C until analyzed. Metabolites were isolated from serum using methanol as described previously [19]; 700 µl extractant containing internal standard 1 (methanol: acetonitrile: water = 4:2:1, V/V/V) was added to 100 µl serum including quality control (QC) sample. After vortexing for 1 min, the samples were incubated at − 20 ℃ for 2 h and then centrifuged for 15 min at 25,000 g, 4 ℃; 600 µl of the supernatant was collected and dried. Then, 180 µl methanol:pure water (1:1 v/v) was added to resuspend that sample. After vortexing for 10 min, it was then centrifuged for 15 min at 25,000 g, 4 ℃. The supernatant was collected for LC-MS/MS analysis.
UPLC-MS analysis
This experiment used Waters UPLC I-Class Plus (Waters, USA) tandom Q Exactive high-resolution mass spectrometer (Thermo Fisher Scientific, USA) for separation and detection of metabolites.
Chromatographic separation was performed on a Waters ACQUITY UPLC BEH C18 column (1.7 μm, 2.1 × 100 mm, Waters, USA), and the column temperature was maintained at 45 °C. The mobile phase consisted of 0.1% formic acid (A) and acetonitrile (B) in positive mode, and in negative mode, the mobile phase consisted of 10 mM ammonium formate (A) and acetonitrile (B). The gradient conditions were as follows: 0–1 min, 2% B; 1–9 min, 2–98% B; 9–12 min, 98% B; 12–12.1 min, 98% B to 2% B; 12.1–15 min, 2% B. The flow rate was 0.35 ml/min, and the injection volume was 5 μl.
Q Exactive (Thermo Fisher Scientific, USA) was used to perform primary and secondary mass spectrometry data acquisition. The full scan range was 70–1050 m/z with a resolution of 70000, and the automatic gain control (AGC) target for MS acquisitions was set to 3e6 with a maximum ion injection time of 100 ms. The top three precursors were selected for subsequent MSMS fragmentation with a maximum ion injection time of 50 ms and resolution of 17,500; the AGC was 1e5. The stepped normalized collision energy was set to 20, 40 and 60 eV. ESI parameters were: sheath gas flow rate of 40, auxiliary gas flow rate of 10, positive ion mode spray voltage (|KV|) of 3.80, negative ion mode spray voltage (|KV|) of 3.20, capillary temperature of 320 °C and auxiliary gas heater temperature of 350 °C.
Metabolite ion peak extraction and metabolite identification
After quality control and pre-processing, ion peaks were extracted. The off-line data of mass spectrometry were imported into Compound Discoverer 3.3 (Thermo Fisher Scientific, USA) software, and metabolites were identified by comparison with standard references in the Beijing Genomics Institute metabolome database (BMBD), mzcloud database and chemspider online database. Then, the data were further analyzed and processed.
Data preprocessing
The result files from Compound Discoverer were input into MetaXto MetaX. Data preprocessing included: (i) data normalization to obtain relative peak areas by probabilistic quotient normalization (PQN) [20]; (ii) quality control-based robust loess signal correction to correct batch effect [19]; (iii) removal of metabolites with a coefficient of variation > 30% on the relative peak area in QC samples.
Quality control
Principal component analysis (PCA) was used to observe the overall distribution of QC for each group of samples and stability of the entire analysis process. PCA is an unsupervised pattern recognition method for statistical analysis of multidimensional data. Through orthogonal transformation, a group of variables that may be correlated is converted into a group of linear unrelated variables, which are called principal components after the transformation. This method is used to study how a few principal components can reveal the internal structure between multiple variables while keeping the original variable information. Log transformation and Pareto scaling were mainly used to compute PCA. The PCA plot reflects the real distribution of samples and was mainly used to observe the separation trend between sample groups and whether there were abnormal samples and to reflect the variability between and within groups from the original data. For QC samples, the better the QC samples aggregate, the more stable the instrument is and the better the reproducibility of the collected data.
Overall analysis
The identified metabolites were classified and annotated using Human Metabolome Database (HMBD). The HMBD contains chemical, molecular biological/biochemical and clinical information on metabolites and supports metabolic pathway and spectral searches. The function annotation of pathway was performed using Kyoto Encyclopedia of Genes and Genomes (KEGG) through the KEGG PATHWAY database (core of KEGG database) to determine the main biochemical metabolic pathways involved in the metabolites. Numerous metabolic pathways and the relationships between them can be found there. To investigate potential metabolite pathways involved in the progression of S. japonicum infection, we conducted KEGG pathway analysis with the differently expressed metabolites in different stages of S. japonicum infection. The differential abundance (DA) score was calculated according a previous publication [21]. DA = (no. of increased metabolites − no. of decreased metabolites)/no. of measured metabolites in a given pathway. The DA score can indicate the average and gross alterations of all measured metabolites in a pathway. A DA score of 1 indicates all metabolites in a pathway increased in abundance, while a DA score of − 1 denotes all metabolites decreased.
Screening the differences between groups
Partial least squares discriminant analysis (PLS-DA) and orthogonal partial least squares discriminant analysis (OPLS-DA) performed on preprocessed data in the R programming language version 4.0.3 (R Foundation for Statistical Computing, 2020) were used to evaluate the metabolite alterations between advanced and chronic S. japonicum infection. Unlike PCA, PLS-DA is a supervised statistical method. It can reflect the differences between classification groups better. This method uses partial least squares regression to establish a model between metabolite expression and sample categories to predict sample categories. After log 2 transformation of the data, a PLS-DA model was established between the comparative analysis groups (two groups of samples), the scaling method was Pareto, and a sevenfold cross validation was used for validation when building the model.
OPLS-DA is a combination of OSC and PLS-DA; it is an extension of PLS-DA, which can decompose X matrix information into two types of information related to Y and unrelated to Y, remove information irrelevant to classification and effectively reduce the complexity of the model without reducing the predictive ability of the model; this enhances the explanatory power of the model. Two hundred iterations of response permutation tests were performed to validate the OPLS-DA model and prevent model overfitting. Also, variable importance in projection (VIP) was used to measure the impact strength and explanatory power of each metabolite expression pattern on the classification and discrimination of each group of samples, which helped to screen the metabolic biomarkers. Generally, a > 1 VIP value could indicate metabolites that had a significant effect on distinct sample categories. VIP values, indicating the importance of variables for group discrimination, were obtained from the OPLS-DA model. Fold change analysis (FC analysis) and t-test were performed on the data. Fold change (FC) was obtained by FC analysis; P-value was obtained by t-test. P-value was used to assess statistical significance between two analyzed groups. Metabolites with VIP ≥ 1, fold change ≥ 1.2 or ≤ 0.83 and P-value < 0.05 were considered differential metabolites.
ROC curve and statistical analyses
ROC curve analysis was performed using SPSS software version 25.0 (SPSS Inc., Chicago, IL) to identify potential candidate biomarkers (area under the curve > 0.8) in the discovery set and to assess the predictive power of markers in metabolomics studies [22]. Pathway analyses of all potential candidate biomarkers in the discovery set were generated using MetaboAnalyst 5.0 software. A two-tailed P < 0.05 was considered statistically significant.
Patient demographic data were examined using SPSS software version 25.0 (SPSS Inc., Chicago, IL). Differences in gender were determined by using Fisher’s exact test. One-way analysis of variance (ANOVA) and Mann-Whitney U test were performed for multiple comparisons. Receiver operator characteristic (ROC) analyses were used to predict the diagnostic efficiency of the different metabolites. A P-value < 0.05 was considered a criterion for statistically significant differences.
Result
Clinical characteristics of participants
Sixty-five patients met the inclusion criteria and were divided into three different groups. No significant differences in age or sex were observed among the three groups (Table 1). The total bile acid level in patients with S. japonicum infection was higher than that in healthy individuals. Furthermore, patients with advanced S. japonicum infection exhibited even higher levels of bile acids compared to those with chronic infection (Table 1). Compared to the healthy control group, patients with S. japonicum infection showed an increased level of blood glucose, total cholesterol and low-density lipoprotein (Table 1). Nevertheless, there were no significant differences between ASJ and CSJ. The clinical characteristics of the participants are shown in Table 1.
Metabolomic analysis of S. japonicum infection sera
The overlaid base peak intensity chromatograms of the QC samples in positive and negative ion modes verified the reliability and stability of the system (Fig. 1a). PCA indicates variability among CON, CSJ and ASJ (Fig. 1b). The observed differences among the samples, as depicted in the figure, can be attributed to variations in metabolite levels. Then, PLS-DA was applied to these three groups of biological samples to establish the relationship model between metabolite expression and sample class, which can realize the modeling and prediction of sample class (Fig. 1c). OPLS-DA was performed on the collected data to further identify variations in the three groups. Two hundred iterations of the response permutation test showed negative Q2 intercept values (Fig. 1d), indicating that the OPLS-DA model was not overfit. OPLS-DA score plots showed great separation, confirming significant differences among the three groups (Fig. 1e).
Differential metabolites in different stages of S. japonicum infection
Herein, we conducted a differential analysis of high-quality metabolites to gain quantitative insights into their relative levels among samples from the different groups. Differential metabolites were identified based on specific criteria: FC ≥ 1.2 or FC ≤ 0.83, P-value < 0.05, and VIP ≥ 1.0. Comparing samples of CSJ with healthy controls, we found 185 upregulated metabolites and 297 downregulated metabolites. The significant metabolites mainly belonged to categories such as benzene and derivatives, organoheterocyclic compounds, glycerophospholipids, bile acids, alcohols and derivatives, and fatty acyls (Fig. 2a and S1). Furthermore, comparing the ASJ to the CSJ, we identified 199 upregulated metabolites and 207 downregulated metabolites in ASJ. These differentially abundant metabolites were categorized as steroids and derivatives, glycerophospholipids, amino acids, peptides and analogs, organoheterocyclic compounds, organic acids and fatty acyls (Fig. 2b and S2). Comparing the ASJ and healthy control groups, we found 187 upregulated metabolites and 383 downregulated metabolites in ASJ. These metabolites were classified in categories such as benzene and derivatives, steroids and derivatives, organoheterocyclic compounds, fatty acyls, glycerophospholipids, terpenoids and carbohydrates (Fig. 2c and S3).
Functional enrichment analysis of the differentially expressed metabolites
Using KEGG, we expressed metabolites in different stages of S. japonicum infection and found that 10 pathways are significantly enriched in CSJ and ASJ, including steroid hormone biosynthesis, cholesterol metabolism, bile secretion and primary bile acid biosynthesis (Table 2 and Fig. 3a). DA score is a pathway-based metabolite change analysis method, which can help to understand the overall changes in metabolite pathways [21]. Compared with the CSJ, DA score analysis indicated that the steroid hormone biosynthesis pathway was significantly downregulated in ASJ; pathways such as cholesterol metabolism, bile secretion, ATP-binding cassette (ABC) transporters, amino acid biosynthesis and primary bile acid biosynthesis were significantly upregulated (Fig. 3b).
Metabolites associated with the progression of S. japonicum infection
To identify metabolites that have potential to assist diagnosis of different stages of S. japonicum infection, we performed further analysis on the different expressed metabolites. A Venn diagram analysis revealed that there were 133 different metabolites when comparing CON and CSJ samples, 233 when comparing CSJ and ASJ samples and 174 when comparing ASJ and CON samples (Figure S4 and Table S1–3). These findings suggest that there are relatively fewer different metabolites associated with progression to advanced schistosomiasis. Specifically, we noticed that glycocholic acid (GCA), glycochenodeoxycholate (GCDCA) and taurochenodeoxycholic acid (TCDCA) exhibited a gradual rise in concentration as the disease progressed (Fig. 4a). To evaluate the predictive ability of potential biomarkers for advanced S. japonicum infection, we performed ROC analysis on different expressed metabolites. Our study identified three metabolites with area under the curve (AUC) > 0.8, including GCA, GCDCA and TCDCA, concentrated in cholesterol metabolism, biliary secretion and primary bile acid biosynthesis (Table 2). Thus, ROC analysis indicated that GCA, GCDCA and TCDCA seemed to have potential to act as candidate biomarkers for diagnosis of advanced S. japonicum infection (Fig. 4b).
Discussion
Schistosomiasis remains a public health problem in the world. Although there have been advances in the diagnosis of schistosomiasis, the methods used in the clinic are mainly conventional diagnostic technologies, characterized by low sensitivity and specificity. The main goal of controlling schistosomiasis is not only the development of novel diagnostic tools but also the implementation of the WASH measurements among mass drug administrations, which are stated in WHO Roadmap 2030 [23]. Moreover, the mechanism of the progression of S. japonicum infection from chronic to advanced stage is still not well defined. The metabolome has been reported to be important in various disease diagnoses and prognoses, and few studies have also addressed the metabolite status of clinical S. japonicum infection patients [13], but the role of metabolite profiles in the progression of S. japonicum infection in the clinic is not clear. In this study, we used UPLC-MS-based method to detect the serum metabolite profiles in patients with S. japonicum infection. Our results demonstrated that patients in different stages of S. japonicum infection had different metabolite signatures. In addition, our study revealed that steroid hormone biosynthesis, cholesterol metabolism and amino acid metabolism possibly become involved in the progression from chronic to advanced S. japonicum infection. Furthermore, ROC analysis indicated that GCA, GCDCA and TCDCA have potential to be novel metabolite biomarkers to differentiate patients with chronic from advanced S. japonicum infection.
Fluctuations in steroid hormone levels can greatly affect the diagnosis and prognosis of liver disease [24]. Schistosoma japonicum infection leads to reprogramming of glucose and lipid metabolism in the mouse colon [25]. The major differentially expressed metabolite molecules observed in this study were lipids, consistent with previous studies [26, 27]. Schistosoma japonicum expresses steroid hormone receptors on the surface, displaying the ability to process steroid hormones to anabolic derivatives [28]. The genome of Schistosoma encodes putative enzymes that convert progesterone and pregnenolone to estriol, estrone, androsterone and testosterone [28]. Therefore, Schistosoma may exploit this pathway during the developmental stages of their parasitic hosts. Our results indicated a significant decrease in the steroid hormone biosynthesis pathway in patients with advanced S. japonicum infection, particularly in the levels of sex hormones, such as androsterone, etiocholanolone and its derivatives, as well as dehydroepiandrosterone. Research on the direct effects of androsterone levels on liver disease was limited. However, studies have established a correlation between the decline of sexual hormone levels and late stages of fibrosis in non-alcoholic fatty liver disease (NAFLD), including androsterone sulfate, epiandrosterone sulfate and dehydroepiandrosterone sulfate [29]. Epiandrosterone is a natural byproduct of dehydroepiandrosterone (DHEA) that is formed through 5α-reductase enzyme. Typically, DHEA is present in the form of dehydroepiandrosterone sulfate (DHEA-S). A previous study found that patients with advanced nonalcoholic fatty liver disease (NASH) had lower levels of DHEA-S [30]. Earlier studies have suggested that administering DHEA to experimental animals can alleviate liver damage by hindering the activity of various inflammatory mediators, including tumor necrosis factor-α and macrophage mitogen inhibitory factor, and by averting the elevation of serum ALT levels [31]. According to the aforementioned research, the decrease in steroid hormone synthesis may diminish the protective effects of these hormones on the liver, consequently aggravating the progression of liver damage related to S. japonicum infection.
The clinical data showed that advanced S. japonicum infection patients displayed notably elevated serum bile acid levels. Meanwhile, the metabolomics investigation further underscored the substantial upregulation of pathways associated with cholesterol metabolism, primary bile acid synthesis and bile secretion in these advanced S. japonicum infection patients. However, it remains unclear whether the increased serum bile acids were the result of liver damage induced by S. japonicum eggs or the causes of progressed liver disease; thus, further study is needed. Earlier studies have already indicated that disruptions in cholesterol metabolism can lead directly to liver damage and potentially trigger the development of metabolite liver diseases [32]. Another consequence of heightened cholesterol metabolism is the escalation of bile acid synthesis [33]. Conditions like cholestasis and metabolite disorders have been linked to the accumulation of bile acids at elevated concentrations, which in turn could exert cytotoxic and pro-inflammatory effects on liver tissues [33]. Therefore, the interplay of cholesterol metabolism, bile acid synthesis and bile secretion pathways could indicate crucial mechanisms underlying the progression from chronic to advanced stages of S. japonicum infection-related liver disease. Of note, recent animal studies demonstrated that S. japonicum infection of mice lacking Farnesoid X receptor (FXR) in hepatocytes increased bile acid toxicity and inhibited autophagy in host hepatocytes, thus possibly accelerating the progression of S. japonicum infection by promoting hepatocyte injury [34]. Therefore, it is feasible to treat advanced S. japonicum infection by inhibiting cholesterol metabolism, bile acid synthesis and bile secretion-related pathways; however, it requires further animal studies and needs to be validated in patients.
In the pathways that were upregulated, we observed an unexpected surge in the expression of ABC transporters among patients with advanced S. japonicum infection. Previous study has indicated that ABC transporters play a pivotal role in modulating the sensitivity to drugs like praziquantel (PZQ) in Schistosoma infection as well as other parasitic worms [35]. Thus, ABC transporter inhibitors seemed to have potential benefits for the treatment of S. japonicum infection, although this warrants further study. In addition, our analysis revealed increased expression of specific amino acids including l-asparagine, l-threonine, l-histidine and l-methionine in ASJ (Table S2). Similar disruptions in amino acid metabolism have been identified not only in hosts infected with S. japonicum infection but also in patients suffering from a range of metabolite liver diseases, along with their corresponding disease models [36,37,38,39,40,41,42]. This suggests that targeting these amino acid metabolite abnormalities seemed have potential to provide some evidence for developing novel therapeutic strategies to mitigate liver damage in advanced S. japonicum infection patients.
The cytotoxic effect of bile acids is dependent on the composition of the bile acid pool, particularly hydrophobic bile acids exhibiting greater toxicity [43]. Bile acids currently hold significant promise as predictive markers in the progression trajectory of numerous liver diseases [44]. For instance, alterations in bile acid composition are strongly linked with conditions like cirrhosis [45]. Patients with cirrhosis have a lower proportion of bile acids compared to other patients, while the proportion of primary bile acids is higher in patients with advanced liver cirrhosis [46,47,48,49]. Recent research on serum metabolite groups has identified 16 bile acids, such as GCA, GCDCA and TCDCA, which are increased in patients with severe hepatitis B (Child-Pugh score) [50]. Besides, a recent study found that the bile acid levels of NASH patients were higher than in healthy controls [50, 51]. Moreover, in patients with cirrhosis, the levels of primary conjugated bile acids such as GCA, GCDCA and TCDCA were significantly elevated compared with non-cirrhotic patients [51]. These results suggested that an increase in primary conjugated bile acids is a characteristic feature of cirrhosis, irrespective of the presence or absence of hepatocellular carcinoma (HCC) [51]. In addition, the study also discovered the expression of GCA, GCDCA and TCDCA induced TGFβ1 in vitro hepatic stellate cells (HSCs), which indicated that primary conjugated bile acid is related to HSC activation and fibrosis in liver cirrhosis [51]. Consistent with previous studies, our study found that the robust area AUC values of primary conjugated bile acids like GCA, GCDCA and TCDCA were > 0.8. This suggests that variations in these bile acid levels hold promising diagnostic potential for assessing the extent of liver damage in both advanced and chronic S. japonicum infection. However, to pinpoint the most pivotal differential metabolites in the progression from chronic to advanced S. japonicum infection, further studies are required and need to be validated by large sample size clinical investigation. This will allow for a more in-depth understanding of the metabolite mechanisms underlying the transition from chronic to advanced S. japonicum infection.
Conclusions
UPLC-MS technology represents a novel approach for detecting changes in the serum metabolite profile of individuals with chronic or advanced S. japonicum infection. Our study revealed that patients in different stages of S. japonicum infection had different metabolome signatures. Furthermore, we found that the steroid hormone biosynthesis, cholesterol metabolism and partial amino acid metabolism correlated with the progression from chronic to advanced S. japonicum infection. Notably, our results indicated that GCA, GCDCA and TCDCA seemed to have potential to be novel metabolite biomarkers to differentiate patients with chronic S. japonicum infection from those with advanced S. japonicum infection, and this research also shed light on the mechanism of the progression from chronic to advanced S. japonicum infection in terms of metabolism. However, further studies in larger groups of patients with S. japonicum infection are needed to validate the potential role of the metabolites in the progression of S. japonicum infection. Overall, our study demonstrates the potential utility of UPLC-MS technology for identifying biomarkers associated with S. japonicum infection progression, which could have important implications for clinical management and the development of more effective treatments.
Availability of data and materials
The untargeted metabolomic data used in this publication are available from the corresponding author upon reasonable request.
Abbreviations
- S. japonicum :
-
Schistosoma japonicum
- UPLC-MS:
-
Ultra-performance liquid chromatography-mass spectrometry
- ROC:
-
Receiver operator characteristic
- GCA:
-
Glycocholic acid
- GCDCA:
-
Glycochenodeoxycholate
- TCDCA:
-
Taurochenodeoxycholic acid
- KK:
-
Kato-Katz
- PCA:
-
Principal component analysis
- PLS-DA:
-
Partial least squares discriminant analysis
- OPLS-DA:
-
Orthogonal partial least squares discriminant analysis
- VIP:
-
Variable importance for projected
- FC:
-
Fold change
- KEGG:
-
Kyoto Encyclopedia of Genes and Genomes
- ABC:
-
ATP-binding cassette
- AUC:
-
Area under the curve
References
Utzinger J, Becker SL, Knopp S, Blum J, Neumayr AL, Keiser J. Hatz CF Neglected tropical diseases: diagnosis, clinical management, treatment and control. Swiss Med Week. 2012;142:w13727.
Adenowo AF, Oyinloye BE, Ogunyinka BI, Kappo AP. Impact of human schistosomiasis in sub-Saharan Africa. Braz J Infect Dis. 2015;19:196–205.
Gray DJ, McManus DP, Li Y, Williams GM, Bergquist R, Ross AG. Schistosomiasis elimination: lessons from the past guide the future. Lancet Infect Dis. 2010;10:733–6.
Schistosomiasis: number of people receiving preventive chemotherapy in 2012. Releve epidemiologique hebdomadaire. 2014;89:21–8.
Freer JB, Bourke CD, Durhuus GH, Kjetland EF, Prendergast AJ. Schistosomiasis in the first 1000 days. Lancet Infect Dis. 2018;18:e193–203.
Huang Y, Lu J, Xu Y, Xiong C, Tong D, Hu N, et al. Xiaochaihu decorction relieves liver fibrosis caused by Schistosoma japonicum infection via the HSP47/TGF-β pathway. Parasit Vectors. 2020;13:254.
Von Lichtenberg F, Erickson DG, Sadun EH. Comparative histopathology of Schistosome granulomas in the hamster. Am J Pathol. 1973;72:149–78.
Zheng L, Wang L, Hu Y, Yi J, Wan L, Shen Y, et al. Higher frequency of circulating Vδ1 γδT cells in patients with advanced Schistosomiasis. Parasite Immunol. 2021;43:e12871.
MacGregor SR, McManus DP, Sivakumaran H, Egwang TG, Adriko M, Cai P, et al. Development of CRISPR/Cas13a-based assays for the diagnosis of Schistosomiasis. EBioMedicine. 2023;94:104730.
Amoah AS, Hoekstra PT, Casacuberta-Partal M, Coffeng LE, Corstjens PL, Greco B, et al. Sensitive diagnostic tools and targeted drug administration strategies are needed to eliminate schistosomiasis. Lancet Infect Dis. 2020;20:e165–72.
He P, Gordon CA, Williams GM, Li Y, Wang Y, Hu J, et al. Real-time PCR diagnosis of Schistosoma japonicum in low transmission areas of China. Infect Dis Poverty. 2018;7:8.
Weerakoon KG, McManus DP. Cell-free DNA as a diagnostic tool for human parasitic infections. Trends Parasitol. 2016;32:378–91.
Huang Y, Xu Y, Huang Y, Sun F, Tian H, Hu N, et al. Identification of newly developed advanced schistosomiasis with MALDI-TOF mass spectrometry and ClinProTools analysis. Parasite. 2019;26:33.
Huang Y, Wu Q, Zhao L, Xiong C, Xu Y, Dong X, et al. UHPLC-MS-based metabolomics analysis reveals the process of schistosomiasis in mice. Front Microbiol. 2020;11:1517.
Wu J, Xu W, Ming Z, Dong H, Tang H, Wang Y. Metabolic changes reveal the development of schistosomiasis in mice. PLoS Negl Trop Dis. 2010;4(8):e807.
Tounta V, Liu Y, Cheyne A, Larrouy-Maumus G. Metabolomics in infectious diseases and drug discovery. Mol Omics. 2021;17:376–93.
Chinese Center for Disease Control and Prevention ea. Diagnostic Criteria for Schistosomiasis (WS261–2006). 2006.
Li Y, Mei L, Qiang J, Ju S, Zhao S. Magnetic resonance spectroscopy for evaluating portal-systemic encephalopathy in patients with chronic hepatic Schistosomiasis Japonicum. PLoS Negl Trop Dis. 2016;10:e0005232.
Dunn WB, Broadhurst D, Begley P, Zelena E, Francis-McIntyre S, Anderson N, et al. Procedures for large-scale metabolic profiling of serum and plasma using gas chromatography and liquid chromatography coupled to mass spectrometry. Nat Protoc. 2011;6:1060–83.
Di Guida R, Engel J, Allwood JW, Weber RJ, Jones MR, Sommer U, et al. Non-targeted UHPLC-MS metabolomic data processing methods: a comparative investigation of normalisation, missing value imputation, transformation and scaling. Metabolomics. 2016;12:93.
Hakimi AA, Reznik ED, Lee CH, Creighton CJ, Brannon AR, Luna A, et al. An integrated metabolic atlas of clear cell renal cell carcinoma. Cancer Cell. 2016;29:104–16.
Li R, Yang L, Guan S, Lin M, Lai H, Liu K, et al. UPLC-MS-based serum metabolic profiling reveals potential biomarkers for predicting propofol responsiveness in females. J Proteome Res. 2021;20:4578–88.
Venkatesan P. Defeating meningitis by 2030: the WHO roadmap. Lancet Infect Dis. 2021;21:1635.
Kur P, Kolasa-Wołosiuk A, Misiakiewicz-Has K, Wiszniewska B. Sex hormone-dependent physiology and diseases of liver. Int J Environ Res Public Health. 2020;17(8):2620.
Yang X, Ding W, Qian X, Jiang P, Chen Q, Zhang X, et al. Schistosoma japonicum infection leads to the reprogramming of glucose and lipid metabolism in the colon of mice. Front Vet Sci. 2021;8:645807.
Liu Y, Zhang PP, Li JH, Li H, Zhou C, Zhang Y, et al. Correction: Association between serum lipid profile and liver fibrosis in patients infected with Schistosoma japonicum. Parasit Vectors. 2022;15:330.
Yang H, Li H, Chen W, Mei Z, Yuan Y, Wang X, et al. Therapeutic effect of Schistosoma japonicum cystatin on atherosclerotic renal damage. Front Cell Dev Biol. 2021;9:760980.
Wu W, LoVerde PT. Nuclear hormone receptors in parasitic platyhelminths. Mol Biochem Parasitol. 2019;233:111218.
Caussy C, Ajmera VH, Puri P, Hsu CL-S, Bassirian S, Mgdsyan M, et al. Serum metabolites detect the presence of advanced fibrosis in derivation and validation cohorts of patients with non-alcoholic fatty liver disease. Gut. 2019;68:1884–92.
Charlton M, Angulo P, Chalasani N, Merriman R, Viker K, Charatcharoenwitthaya P, et al. Low circulating levels of dehydroepiandrosterone in histologically advanced nonalcoholic fatty liver disease. Hepatology. 2008;47:484–92.
Abd-Ella EM, El-Kott AF, El-Kenawy AE, Khalifa HS, Bin-Meferij MM, Alramlawy AM. Dehydroepiandrosterone protects against acetaminophen-induced liver damage in rats by upregulation of Bcl-2 and activation of sirt signalling. J Physiol Pharmacol. 2020;71(6):781–92.
Goicoechea L, de la Conde Rosa L, Torres S, GarcÃa-Ruiz C, Fernández-Checa JC. Mitochondrial cholesterol: metabolism and impact on redox biology and disease. Redox Biol. 2023;61:102643.
Mizuochi T, Takei H, Nittono H, Kimura A. Inborn errors of bile acid metabolism in Japan. Pediatr int. 2023;65:e15490.
Zhang B, Li J, Zong X, Wang J, Xin L, Song H, et al. FXR deficiency in hepatocytes disrupts the bile acid homeostasis and inhibits autophagy to promote liver injury in Schistosoma japonicum-infected mice. PLoS Negl Trop Dis. 2022;16:e0010651.
Greenberg RM. Schistosome ABC multidrug transporters: from pharmacology to physiology. Int J Parasitol Drugs Drug Resist. 2014;4:301–9.
Li JV, Holmes E, Saric J, Keiser J, Dirnhofer S, Utzinger J, et al. Metabolic profiling of a Schistosoma mansoni infection in mouse tissues using magic angle spinning-nuclear magnetic resonance spectroscopy. Int J Parasitol. 2009;39:547–58.
Watanabe A, Hayashi S, Higashi T, Obata T, Sakata T, Takei N, et al. Characteristics change in serum amino acid levels in different types of hepatic encephalopathy. Gastroenterol Jpn. 1982;17:218–23.
An Z, Hu T, Lv Y, Li P, Liu L. Targeted amino acid and related amines analysis based on iTRAQ -LCMS/ MS for discovering potential hepatotoxicity biomarkers. J Pharm Biomed Anal. 2020;178:112812.
Locasale JW. Serine, glycine and one-carbon units: cancer metabolism in full circle. Nat Rev Cancer. 2013;13:572–83.
Meiser J, Tumanov S, Maddocks O, Labuschagne CF, Athineos D, Van Den Broek N, et al. Serine one-carbon catabolism with formate overflow. Sci Adv. 2016;2:e1601273.
Li Z, Wang F, Liang B, Su Y, Sun S, Xia S, et al. Methionine metabolism in chronic liver diseases: an update on molecular mechanism and therapeutic implication. Signal Transd Targeted Ther. 2020;5:280.
del Maria Bas J, RodrÃguez B, Puiggròs F, Mariné S, RodrÃguez MA, Morina D, et al. Hepatic accumulation of S-adenosylmethionine in hamsters with non-alcoholic fatty liver disease associated with metabolic syndrome under selenium and vitamin E deficiency. Clin Sci. 2019;133:409–23.
Perez MJ, Briz O. Bile-acid-induced cell injury and protection. World J Gastroenterol. 2009;15:1677–89.
Lai J, Luo L, Zhou T, Feng X, Ye J, Zhong B. Alterations in circulating bile acids in metabolic dysfunction-associated steatotic liver disease: a systematic review and meta-analysis. Biomolecules. 2023;13(9):1356.
Islam KS, Fukiya S, Hagio M, Fujii N, Ishizuka S, Ooka T, et al. Bile acid is a host factor that regulates the composition of the cecal microbiota in rats. Gastroenterology. 2011;141:1773–81.
Vlahcevic ZR, Buhac I, Bell CC Jr, Swell L. Abnormal metabolism of secondary bile acids in patients with cirrhosis. Gut. 1970;11:420–2.
Yin P, Wan D, Zhao C, Chen J, Zhao X, Wang W, et al. A metabonomic study of hepatitis B-induced liver cirrhosis and hepatocellular carcinoma by using RP-LC and HILIC coupled with mass spectrometry. Mol Biosyst. 2009;5:868–76.
Wang X, Wang X, Xie G, Zhou M, Yu H, Lin Y, et al. Urinary metabolite variation is associated with pathological progression of the post-hepatitis B cirrhosis patients. J Proteome Res. 2012;11:3838–47.
Chen T, Xie G, Wang X, Fan J, Qiu Y, Zheng X, et al. Serum and urine metabolite profiling reveals potential biomarkers of human hepatocellular carcinoma. Mol Cell Proteomics: MCP. 2011;10:M110004945.
Wang X, Xie G, Zhao A, Zheng X, Huang F, Wang Y, et al. Serum bile acids are associated with pathological progression of hepatitis B-induced cirrhosis. J Proteome Res. 2016;15:1126–34.
Sydor S, Best J, Messerschmidt I, Manka P, Vilchez-Vargas R, Brodesser S, et al. Altered microbiota diversity and bile acid signaling in cirrhotic and noncirrhotic NASH-HCC. Clin Transl Gastroenterol. 2020;11:e00131.
Acknowledgements
We thank Prof. Guanghui Ren for his insightful and thoughtful suggestions regarding this manuscript.
Funding
This study was supported by National Natural Science Foundation of China (grant no. 81771722 to Y Ming; grant no. 81901630 to J Li), Natural Science Foundation of Hunan Province of China (grant no. 2021JJ40959 to Y Zhu) and the Key Research and Development Plan of Hunan Province (grant no. 2021SK2032 to Y Ming).
Author information
Authors and Affiliations
Contributions
The study was conceived and designed by JL, JJ and YM. The patient recruitment and screening were completed by JJ, YZ and JZ. Data were collected and analyzed by YZ and JJ. The manuscript was written by JL, JJ and YM. YM supervised the whole project. YM, JL and YZ provided funding support for the study. All the authors have read and approved the submitted manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Institutional Review Board (IRB) of The Third Xiangya Hospital of Central South University (21149), and written informed consent was obtained from all participants.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
13071_2024_6429_MOESM4_ESM.pdf
Supplementary Material 4: Fig. S4. Venn diagram results among the three groups, CON vs. CSJ, CSJ vs. ASJ and ASJ vs. CON.
13071_2024_6429_MOESM5_ESM.xlsx
Supplementary Material 5: Table S1.Molecular information of differential metabolites between healthy_control and CSJ groups_ESM.
13071_2024_6429_MOESM7_ESM.xlsx
Supplementary Material 7: Table S3.Molecular information of differential metabolites between healthy_control and ASJ groups_ESM.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Li, J., Jiang, J., Zhu, Y. et al. Metabolomics analysis of patients with Schistosoma japonicum infection based on UPLC-MS method. Parasites Vectors 17, 350 (2024). https://doi.org/10.1186/s13071-024-06429-9
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s13071-024-06429-9