- Open Access
Deciphering the metabolic perturbation in hepatic alveolar echinococcosis: a 1H NMR-based metabolomics study
Parasites & Vectors volume 12, Article number: 300 (2019)
Hepatic alveolar echinococcosis (HAE) is caused by the growth of Echinococcus multilocularis larvae in the liver. It is a chronic and potentially lethal parasitic disease. Early stage diagnosis for this disease is currently not available due to its long asymptomatic incubation period. In this study, a proton nuclear magnetic resonance (1H NMR)-based metabolomics approach was applied in conjunction with multivariate statistical analysis to investigate the altered metabolic profiles in blood serum and urine samples obtained from HAE patients. The aim of the study was to identify the metabolic signatures associated with HAE.
A total of 21 distinct metabolic differences between HAE patients and healthy individuals were identified, and they are associated with perturbations in amino acid metabolism, energy metabolism, glyoxylate and dicarboxylate metabolism. Furthermore, the present results showed that the Fischer ratio, which is the molar ratio of branched-chain amino acids to aromatic amino acids, was significantly lower (P < 0.001) in the blood serum obtained from the HAE patients than it was in the healthy patient group.
The altered Fischer ratio, together with perturbations in metabolic pathways identified in the present study, may provide new insights into the mechanistic understanding of HAE pathogenesis and potential therapeutic interventions.
Echinococcosis, or hydatid disease, is a near-cosmopolitan zoonotic parasitic infection caused by the larval stage of pathogenic cestode parasites in the genus Echinococcus . Different species of Echinococcus cause different diseases. The main types of echinococcosis include cystic echinococcosis (CE) and alveolar echinococcosis (AE), which are caused by Echinococcus granulosus and Echinococcus multilocularis, respectively . Compared with CE, AE is associated with a lower occurrence, but higher lethality, if not treated with timely and proper treatments .
AE has been reported in its definitive hosts (canids) as well as in humans . In humans, AE affects the liver in approximately 95% of diagnosed cases and is known as hepatic alveolar echinococcosis (HAE) . HAE leads to liver-tissue injury or hepatic failure primarily through an infiltrative behaviour. It resembles the infiltrative proliferation of tumour growth and is clinically known as “worm cancer” or “parasite liver cancer”. After a long period of latent and asymptomatic stage, HAE can progress into a cirrhotic stage , or it can metastasize to other organs (e.g. the lungs and brain) [7,8,9] and cause local organ-function impairment and metastatic infiltration .
HAE is endemic to the northern hemisphere, including North America (Canada ), central Europe (France , Germany , Austria , Poland ), and central Asia (northern Iran , Mongolia  and western China ). Its global occurrence is mainly attributed to chronic anthropogenic influences, including increased globalization of animals and animal products, and altered human-animal interactions .
HAE patients often suffer from symptoms including fatigue, abdominal pain, hepatomegaly, nausea and vomiting. Notably, its long asymptomatic period of the early infection stage leads to difficulty in determining the time-point or place of infection . In general, an HAE diagnosis is based on the integrated information including medical history, contact history with potential livestock hosts, clinical and pathologic findings, imaging examination, nucleic acid detection and serologic tests . Modern imaging techniques, e.g. ultrasonography (US), magnetic resonance imaging (MRI), fluorodeoxyglucose positron emission tomography (FDG-PET), and computed tomography (CT), have been applied for the clinical diagnosis of HAE disease and have been used for long-term follow-up after therapeutic interventions . However, these imaging methods have several limitations: CT is unable to delineate disease-induced perihepatic extension; MRI provides good definition of hydatid lesions but fails to assess hydatid fertility (viability); and FDG-PET provides non-invasive and accurate evaluation of metabolic activity in HAE at the complexity as well as the high expense of administering a radiotracer synthesised with an on-site cyclotron. Additionally, the effective diagnosis of early stage HAE remains difficult even with these advanced imaging techniques. This situation leads to late diagnosis at the middle or advanced stages, causing the patients not to undergo timely curative resection . Thus, an effective early diagnosis and treatment for HAE may help to prevent complications, reduce postoperative reoccurrence, and improve the recovery rate and prognosis.
In the last decade, metabolomics approaches have been broadly used to study various liver diseases to determine potential early biomarkers and perturbed metabolic pathways . However, there are limited reports on the applications of the proton nuclear magnetic resonance (1H NMR)-based metabolomics method on echinococcosis. In the present study, we used high-resolution 1H NMR based metabolomics and multivariate statistical analyses to (i) determine distinct metabolic patterns to differentiate HAE patients from healthy individuals and establish a metabolic fingerprint of biofluids obtained from HAE patients; (ii) provide some clues for the research on the molecular mechanism of HAE disease; (iii) demonstrate the feasibility and effectiveness of using metabolomics in echinococcosis studies; and (iv) provide a diagnostic reference for HAE patients.
Recruitment of participants
The recruited HAE patients were diagnosed at the Department of Hepatopancreatobiliary Surgery of the Affiliated Hospital of Qinghai University between July and September 2016. The diagnoses of HAE relied on the integrated use of imaging examination, serologic test, nucleic acid detection and pathologic observation (Additional file 1: Figure S1). In addition, patients’ medical history and contact with potential animal host of the disease were also taken into consideration. All HAE patients were classified according to the WHO/IWGE classification : P1, peripheral lesions without proximal vascular and/or biliary involvement; P2, central lesions with proximal vascular and/or biliary involvement of one lobe; P3, central lesions with hilar vascular or biliary involvement of both lobes and/or with involvement of two hepatic veins; P4, any liver lesion with extension along the vessels and the biliary tree; N0, no regional involvement; N1, regional involvement of contiguous organs or tissues; M0, no metastasis; M1, metastasis.
To prevent confounding the metabolic effect from other diseases, patients with the following diseases were excluded from the study: diabetes, nephrosis, autoimmunity disease, malignant hepatic tumour, severe hepatorenal dysfunction, and post-transplantation immunoreaction. Prior to this study, no patient underwent surgery or received anti-inflammatory drugs (e.g. non-steroidal anti-inflammatory drugs and corticosteroids). Healthy volunteers from the patients’ families were recruited as the control group.
The sample size used in this study was estimated based on the prior power test with some knowledge from our previous preliminary study. The power analysis was performed using MetaboAnalyst v.4.0 software (http://www.metaboanalyst.ca) [24, 25], and the result indicated that a sample size of 18 was adequate to provide sufficiently high statistical power (Additional file 2: Figure S2).
Sample collection and preparation
All blood and urine samples were collected in the morning, 12 h after the last meal of the previous day (fasting conditions). Blood samples were collected into tubes without anticoagulant and centrifuged at 6700×g at 4 °C for 15 min to obtain blood serum. The first morning urine samples were collected and centrifuged at 8500×g at 4 °C for 15 min and the supernatants were transferred into tubes. The blood serum and urine samples were aliquoted, snap-frozen in liquid nitrogen and stored at − 80 °C until further analysis.
Prior to analysis, an aliquot of 400 μl of blood serum was mixed with 200 μl of phosphate buffer solution (90 mM K2HPO4/NaH2PO4, pH 7.4, 0.9% NaCl, 99.9% D2O). Additionally, 300 μl of urine samples were mixed with a different phosphate buffer solution (300 μl, 1.5 M K2HPO4/NaH2PO4, pH 7.4, 99.9% D2O containing 0.3 mM 3-trimethylsilyl-propionic-2,2,3,3-d4 acid (TSP) as a chemical-shift reference for 0 ppm). D2O was used to provide the NMR spectrometer with a field frequency for locking. Buffered serum and urine samples were then centrifuged at 6700×g at 4 °C for 10 min to remove debris, and 500 μl of supernatant from each 600 μl mixture was transferred to 5-mm NMR tubes. In total, 36 serum and urine samples in NMR tubes were prepared and stored at 4 °C before NMR analysis.
1H NMR experiments were performed using a Bruker NMR system (Bruker Biospin, Karlsruhe, Germany) operating at the proton frequency of 600 MHz. The operating temperature was set at 298 K. The Carr–Purcell–Meiboom–Gill (CPMG) sequence (waiting time ~ π/2 ~ [τ ~ π ~ τ]n ~ acquisition) was used to acquire spectra of blood serum samples with an echo time (τ) of 250 μs and a free relaxation duration (2nτ) of 100 ms. For urine samples, nuclear overhauser effect spectroscopy (NOESY, waiting time ~ π/2 ~ t1 ~ π/2 ~ tm ~ π/2 ~ acquisition) was implemented with a 2 s water suppression and mixing time (tm) of 120 ms. Free induction decays (FIDs) were recorded with 64 scans at a spectral width of 10 kHz. The FIDs were zero-padding to 32 k data points prior to fast Fourier transformation.
Data processing of 1H-NMR spectra
Data pre-processing for the acquired 1H NMR spectra (including Fourier transformation, baseline correction and phase correction) was performed using MestReNova v.8.1.2 software (Mestrelab Research S.L., La Coruña, Spain). For peak alignment, the TSP signal was set as δ 0.00 for urine samples, and left split of the doublet of lactate signals was set as δ 1.336 for serum samples. Residual water signals (serum: δ 4.65–5.15; urine: δ 4.75–5.15), urea resonances (δ 5.70–6.40) and peak-free regions were selectively excluded for further analyses. The remaining spectra over the ranges of δ 0.8–8.5 for blood serum and δ 0.8–9.5 for urine were segmented into bucketed data using self-adaptive integration , and the results were exported as Microsoft Excel files. The data were normalized using the probabilistic quotient normalization (PQN)  method to compensate for overall concentration variations. The peaks in the acquired 1H NMR spectra were assigned based on previously published studies [28, 29], the KEGG database  and the HMDB database . The relative concentrations for all assigned metabolites were evaluated based on their normalized peak integral area. For metabolites that gave rise to multiple peaks, those in the least overlapping spectral region were chosen for quantification.
Multivariate and univariate statistical data analyses
Multivariate analyses of the pre-processed data were performed using SIMCA v.14.1 software (Umetrics, Umeå, Sweden). The data were examined by the non-supervised principal components analysis (PCA, unit variance scaling) and orthogonal partial least squares-discrimination analysis (OPLS-DA, unit variance scaling). A seven-fold cross-validation (i.e. 6/7 samples were used as training data and the remaining 1/7 samples as validation data) and permutation test (200 permutations) were performed and the obtained values of R2 (total explained variation) and Q2 (model predictability) were applied to validate the constructed models [32, 33]. An additional CV-ANOVA analysis was performed and the obtained P-value provided another venue of further model validation. The variable importance projection (VIP) and absolute correlation coefficient (|r|) calculated based on the OPLS-DA analysis were used as parameters for differential selection of metabolites. The VIP parameter, which was denoted as a unitless number, delineates the contribution of each predictor variable to the model and presents the influence of each predictor on the response variables. A larger VIP corresponds to a greater discriminatory power for the metabolites. In addition, metabolite sets were also analysed using Student’s t-test, which constitutes a simple statistical analysis to determine statistically significant metabolic variations in univariate analysis with the transformed P-value . For a particular metabolite, its fold change was calculated based on the ratio of average concentrations between the HAE group and the control group.
In the present study, the results of multivariate statistical analysis were visualized and integrated with volcano plots to identify metabolites with significant difference between the HAE and control groups [35, 36]. In particular, the interactive volcano plot denotes − log10 (P-value) on the vertical axis and log2 (fold change) on the horizontal axis with circles of different sizes and colours for displaying VIP and |r| values, respectively.
Additionally, a post-hoc power analysis with specified significance level α, sample size, and effect size was executed with the online statistics software G*power v.3.1 (http://www.gpower.hhu.de/) [37, 38] to demonstrate the statistical power of the reported results. Moreover, receiver operator characteristic (ROC) curve analysis was applied to the resulting differential metabolites in HAE patients on Metaboanalyst v.4.0 (http://www.metaboanalyst.ca/) . The area under the curve (AUC) of ROC was used to determine the specificity and the sensitivity of altered metabolites to discriminate HAE patients from healthy individuals.
Clinical characteristics of subjects
The metadata of the participants (including age, gender, diameters of lesions, PNM classification and liver function test results) are summarized in Table 1. Comparing the control and the HAE groups, the age and gender of the participants showed no significant difference (P > 0.05). All HAE patients were classified into six categories according to the PNM classification by WHO/IWGE, as stated in methods. Among the participants, patients in the P3N1M1 and P4N1M1 classes showed lung metastasis.
1H-NMR spectra of biological samples
Typical 1H NMR spectra of biofluids collected from the HAE patients and healthy volunteers are shown in Fig. 1. In general, several metabolite classes were detected in the serum and urine samples, including metabolites involved in energy metabolism, neurotransmission, membrane metabolism and osmoregulation. In addition, serum metabolome comprised of a few glycoproteins, while a number of gut microbiota-related metabolites were detected in urine samples.
Pattern recognition analysis
To identify HAE-induced metabolic changes, we performed multivariate analyses on the processed NMR data. First, the data were analysed by principal component analysis (PCA), of which the resulting score plots are shown in Fig. 2. Moderate groupings could be observed between the control and the HAE groups in the PCA score plots. In contrast, a distinct group separation could be achieved when the dataset was analysed using the supervised orthogonal partial least squares discrimination analysis (OPLS-DA). From the OPLS-DA score plots (Fig. 3), the control samples primarily distributed in the left hemisphere, while the HAE group clustered in the right hemisphere. The group separation indicated that HAE induced distinct metabolic changes in human serum and urine. The corresponding probability P-values by CV-ANOVA were also used to assess the statistical significance of metabolic differences between the control and the HAE groups. Moreover, the OPLS-DA models were found to be robust with good explained variations R2Y and predictive powers Q2 following permutation tests (200 permutations).
Determination of differentiating metabolites for HAE
To identify metabolic markers that differentiate the HAE patients from the healthy controls, we next presented the data using a four-dimensional enhanced volcano plot which may also offer intuitive data visualization , as shown in Fig. 4. In the volcano plot, the importance and significance of the metabolic changes were determined using the following criteria: variable importance projection (VIP) > top 30%, absolute correlation coefficient values (|r|) > 0.5, − log10 (P-value) > 2 (i.e. P < 0.01), and absolute log2 (fold change) > 0.20. Generally, identified metabolites with significant changes (those with a significant P-value, high fold change, VIP and |r|) by multivariate statistical analyses tended to locate at the upper-left or upper-right zones of the enhanced volcano plot with larger circle shapes and warmer colours, as marked with abbreviations in Fig. 4.
For blood serum samples, the volcano plot showed that HAE caused significant increases in the concentrations of phenylalanine (Phe), glutamate (Glu), tyrosine (Tyr), formate (For), and lactate (Lac), together with decreases of valine (Val), leucine (Leu), isoleucine (Ile), lysine (Lys), serine (Ser), glutamine (Gln), betaine (Bet), creatine (Cr) and 1-methylhistidine (1-MH) (Fig. 4a). Using a similar approach, urine samples from the HAE group were characterized by significantly higher levels of glutamate (Glu), malate (Mal), glycolate (Gla), 3-hydroxybutyrate (3-HB), γ-aminobutyrate (GABA), 3-methylhistidine (3-MH), and p-hydroxyphenylacetate (p-HPA), as well as lower concentration of creatinine (Cn) (Fig. 4b).
Assessment of distinct metabolites for HAE
The statistical power analyses of the selected differential metabolites are presented in Table 2. First, the absolute difference between two independent means (the HAE and control groups) and the intra-group standard deviation were used to calculate the effect size. The effect size was considered to be large at a value of 0.80 according to the report by Cohen . In addition, post-hoc analysis of the achieved statistical power calculation yielded an average value of 0.92 (SD = 0.09) under the given condition of significance level α = 0.05, sample size n = 18 and specified effective sizes (Additional file 3: Figure S3). The reliability and effectiveness of the present results were validated based on sufficiently high effect sizes and statistical power. Classic univariate ROC curve analyses of the total selected differential metabolites in serum and urine are listed in Table 3.
The sensitivity and specificity were evaluated for each potential metabolite with area under the curve (AUC) values. In principle, an AUC value close to 1 indicates an excellent diagnostic power, and an AUC of 0.5 suggests no diagnostic power. Therefore, metabolites with high AUC values have a better predictive power for identifying HAE patients. In serum samples, ROC curves showed that glutamate (AUC = 0.924), leucine (AUC = 0.920), valine (AUC = 0.920) and phenylalanine (AUC = 0.914) have a high discriminatory ability for HAE. Overall, 10 metabolites had an AUC > 0.8. In urine samples, glutamate, creatinine and p-hydroxyphenylacetate exhibited an AUC > 0.8 (Additional file 4: Figure S4 and Additional file 5: Figure S5). As individual metabolites may not be sufficient to predict the correct phenotype, we performed multivariate ROC curve analyses to investigate combinations of metabolites as potential differential metabolic signatures in serum and urine samples (Fig. 5). ROC curves were generated by Monte-Carlo cross validations. In each Monte-Carlo cross validation, two thirds of the samples were used to evaluate the feature importance. Top important features were used to build classification models, which were validated on the remaining one third of the samples. In serum samples, the best model consisting of the top 10 metabolites (glutamate, leucine, valine and phenylalanine, etc.) showed an AUC value of 0.984 (95% CI: 0.889–1), which did not increase further even when more metabolites were added into the analyses. The results suggest that the combination of the 10 differential metabolites has the potential of distinguishing the HAE patients from the healthy controls (Fig. 5a). In urine samples, metabolite combination increased the AUC value from 0.792 to 0.835 (Fig. 5b). These results suggested the potential roles of these differential metabolites in the HAE pathogenesis.
In addition, we used the 21 characteristic metabolites to build a PLS-DA model (Additional file 6: Figure S6) for HAE identification. The misclassification table of the HAE PLS-DA model with six-fold cross-validation suggests that the combination of 21 characteristic metabolites can serve as a good identifier for HAE samples in the present dataset (Table 4).
Metabolic pathway analysis
Next, metabolic pathway analysis of the identified differential metabolites was performed using MetaboAnalyst v.4.0 software to investigate the most perturbed metabolic pathways in HAE. The analysis suggested perturbations in multiple metabolic pathways including alanine, aspartate and glutamate, arginine and proline metabolism, d-glutamine and d-glutamate metabolism, glycine, serine and threonine metabolism, glyoxylate and dicarboxylate metabolism, lysine metabolism, methane metabolism, phenylalanine metabolism, tyrosine metabolism and valine, leucine and isoleucine metabolism (Fig. 6). These pathways were selected based on an impact value ≥ 0.02 and − log(P) ≥ 5 and considered as the most relevant pathways in HAE (Table 5).
The present study reveals metabolic perturbation in response to HAE disease using 1H NMR based metabolomics. First, 1H NMR spectra and pattern recognition analyses of serum and urine showed substantial changes in the human metabolome due to HAE. Secondly, a total number of 21 metabolites were identified that were able to discriminate the HAE patients from the healthy individuals. Thirdly, metabolic network analyses revealed HAE-induced modulations occur in several metabolic pathways, i.e. amino acid metabolism, energy metabolism, glyoxylate and dicarboxylate metabolism, and methane metabolism (Figs. 6, 7). These metabolic changes may provide new perspectives into the understanding of biological mechanisms that occur during HAE infection.
In the present study, HAE was found to result in significant changes in amino acid metabolism. This is evidenced by decreased levels in lysine, serine, glutamine and branched-chain amino acids (BCAA; i.e. valine, leucine and isoleucine), together with increased levels of glutamate and aromatic amino acids (AAA; e.g. tyrosine and phenylalanine) (Figs. 4, 7). Tyrosine is the first product in phenylalanine catabolism. Previously, it was reported that directionalities of the plasma-concentration shift of these two metabolites are the same and conversion of phenylalanine to tyrosine is an exclusive function of the liver [42, 43]. Thus, the parallel accumulations of AAA (i.e. tyrosine and phenylalanine) in blood serum suggested the possible loss of liver function during HAE infection. As essential amino acids, BCAA account for approximately 20% of our dietary protein intake and are key regulators of protein synthesis in animals and humans. BCAA disposal is tightly mediated by phosphorylation and dephosphorylation of the branched-chain α-keto acid dehydrogenase complex (BCKDC) . The BCKDC is the rate-limiting enzyme in BCAA oxidative catabolism; the possible activation of BCKDC by HAE increases the BCAA catabolism to promote protein synthesis in case of increased protein degradation due to cell damage or necrosis during HAE infection. Moreover, the role of BCAA in the immune function has gained research interest in the last decade . BCAA are oxidized by immune cells as fuel sources and incorporated as precursors for the synthesis of new immune cells, effector molecules and protective molecules . The host immune response induced by HAE may lead to a higher consumption of BCAA to maintain immune function.
In addition, increased urinary concentrations of GABA and 3-methylhistidine were found in the HAE patients. GABA is known to have relaxant effects on muscle tone while 3-methylhistidine is an indicator of muscle catabolism . Taken together, the data indicates that HAE may have metabolic effect on muscle function. In general, 3-methylhistidine is released during muscle-protein degradation and is not reutilized for muscle protein synthesis. This observation suggests that the immune response to HAE requires not only redirection of abundant energy but also metabolic resources from other tissues, especially from skeletal muscles. Hence, the consumption of BCAA in promoting protein synthesis may be due to muscle-protein synthesis. On the other hand, 98% of body creatine is found in skeletal muscle and a constant fraction of the body creatine pool is converted each day to creatinine. Hence, it has long been recognized that urinary creatinine excretion constitutes a good reflection of the skeletal muscle mass . The decline of urine creatinine reflects a corresponding muscle mass loss in the HAE patients. The symptoms of fatigue and weight loss for HAE patients may be attributed to muscle mass loss and altered muscle function.
Being large neutral amino acids, BCAA and AAA can travel across the blood-brain barrier into the brain through the same type of neutral amino acid transporter LAT-1 . Therefore, the raised AAA (tyrosine and phenylalanine) and reduced BCAA (valine, leucine and isoleucine) levels in the blood serum of the HAE patients may contribute to a subsequent increased influx of tyrosine and phenylalanine in the brain during HAE infection. However, the production of important bioamine neurotransmitters (i.e. dopamine, epinephrine and norepinephrine) is relatively insensitive to the cerebral levels of their precursors (i.e. tyrosine and phenylalanine). By contrast, the production of false neurotransmitters (e.g. octopamine or phenylethanolamine) can be stimulated by the cerebral elevated levels of tyrosine and phenylalanine, leading to an increased number of false neurotransmitters. These false neurotransmitters are structurally similar to “authentic” neurotransmitters, but they lack the ability to propagate neural excitation. This may result in imbalances of neurotransmitter synthesis. Homeostasis between excitatory (e.g. glutamate) and inhibitory (e.g. GABA) neurotransmitters is essential for maintaining the normal functioning of the central nervous system (CNS) . In our study, elevations of excitatory glutamate and inhibitory GABA along with glutamine reduction indicated that neurotransmitter recycling disorder occurs in the CNS during HAE infection due to the imbalance between excitation and inhibition. Moreover, alterations of GABA and glutamate in the same direction have been reported to contribute to the pathophysiology of depression based on behavioural observations [51, 52].
Previously, Fischer et al.  reported decreased BCAA and increased AAA levels in the blood during hepatic failure in patients with liver disease. The ratio between the three BCAA (valine, leucine and isoleucine) and two AAA (tyrosine and phenylalanine) has been termed as the Fischer ratio for a quick assessment of liver disease. With a simultaneous consideration over variations of BCAA and AAA, the amino acid imbalance will lead to a lower Fischer ratio in cases of hepatic failures [54,55,56]. In this study, a decreased Fischer ratio of serum samples for the HAE group (P < 0.001) indicated that HAE may have triggered hepatic failure in these patients. Therefore, the Fischer ratio may emerge as a reliable index for diagnosing HAE. Although the present 1H NMR data do not provide absolute concentrations of AAA and BCAA, the Fischer ratio calculated using normalized data was found to provide good discrimination between the HAE patients and the healthy individuals.
Alveoli, host metastasis and massive immune response due to HAE are known to lead to the generation of a large amount of unstable reactive oxygen species (ROS) [11, 57]. Over-generation of ROS induces an oxidative-stress response with upregulation of the tricarboxylic acid (TCA) cycle and glycolysis to produce more ATP for cell maintenance and proliferation. In the present study, the increased levels of malate and lactate in the HAE group support the notion of enhanced TCA cycle and glycolysis. In addition, the present results also showed decreased creatine concentration in the HAE group, which provides further support for the perturbed energy metabolism in liver mitochondria during HAE infection. Creatine acts as an intracellular high-energy phosphate shuttle and plays a crucial role in maintaining cellular energy homeostasis. It can be catabolised to creatinine to provide immediate adequate levels of ATP in cases of energy deficit . Thus, a decreased level of creatine in blood serum may reflect its consumption to alleviate HAE-induced energy deficit. Apart from being an ergogenic aid, creatine has a cytoprotective property as a direct radical scavenger against ROS . Therefore, the lower concentration of creatine in the HAE group can partly be attributed to the antioxidant reactions during HAE infections.
It should be noted that the present study has several limitations. First, a larger number of participants are needed to further increase the reliability and accuracy of the present results. Secondly, an integrated application of NMR with LC–MS or GC–MS will extend the metabolite coverage, and multiple -omic technologies (e.g. genomics and proteomics) can cross-validate and better support the experimental results. Thirdly, data for early stage HAE are currently not available due to the limitation of sample number. Data from different stages of HAE are important for the early detection, prediction and diagnosis of HAE. Finally, a discriminating method that can distinguish HAE from other types of hepatic disease is not available in the present study. If developing such a method, it is crucial to prevent clinical misdiagnoses as HAE exhibits similarities (e.g. lower Fischer ratio) with other hepatic diseases. The possibility of using the Fischer ratio and a unique metabolic signature for HAE to distinguish between different hepatic diseases warrants further study.
The attractive properties of the NMR-based metabolomics approach, i.e. excellent repeatability, short detection duration and multiple metabolite coverage in a single measurement, makes it suitable for the HAE studies. A systematic exploration of multiple metabolites of small molecules helps to reveal the comprehensive metabolic variations induced by HAE infection. The early stage HAE diagnosis may be possible with a combination of both metabolomics and imaging techniques. For instance, chemical-exchange weighted magnetic resonance imaging techniques can be utilized to yield maps weighted by metabolites of interest [60,61,62]. The maps are obtained through an exchange effect with water and are characterized by significant signal enhancement to observe small metabolic variations unnoticeable with other imaging modalities. A general issue for these techniques is the lack of specificity, and modulations from adjacent resonances are inevitably and unfavourably introduced. NMR-based metabolomics with high specificity can thus be combined for improved analyses and better interpretation of the results.
In this study, a 1H NMR-based metabolomics approach was used to investigate metabolic perturbations of HAE. Multivariate statistical analyses have highlighted several characteristic metabolites for HAE, e.g. decreased BCAA (valine, leucine and isoleucine), increased AAA (tyrosine and phenylalanine), increased lactate and decreased creatine, to list a few. The altered Fischer ratio together with specific metabolic changes identified in the present study may provide new insights into the molecular mechanism of HAE, and provide some clues to therapeutic intervention for HAE.
Availability of data and materials
The data supporting the findings of this article are included within the article and its additional files. The 1H NMR spectral raw data have been submitted to the MetaboLights repository under study identifier MTBLS981.
aromatic amino acids
area under the curve
branched-chain amino acids
branched-chain α-keto acid dehydrogenase complex
central nervous system
ANOVA of the cross-validated residuals
fluorodeoxyglucose positron emission tomography
free induction decays
gas chromatography–mass spectrometry
hepatic alveolar echinococcosis
human metabolome database
informal working group on echinococcosis
Kyoto Encyclopedia of Genes and Genomes
liquid chromatography–mass spectrometry
magnetic resonance imaging
nuclear magnetic resonance
orthogonal partial least squares-discrimination analysis
principal components analysis
probabilistic quotient normalization
receiver operator characteristic
reactive oxygen species
variable importance projection
World Health Organization
Taxy JB, Gibson WE, Kaufman MW. Echinococcosis unexpected occurrence and the diagnostic contribution of routine histopathology. Am J Surg Pathol. 2017;41:94–100.
McManus DP, Zhang WB, Li J, Bartley PB. Echinococcosis. Lancet. 2003;362:1295–304.
Budke CM, Casulli A, Kern P, Vuitton DA. Cystic and alveolar echinococcosis: successes and continuing challenges. PLoS Negl Trop Dis. 2017;11:e0005477.
Davidson RK, Romig T, Jenkins E, Tryland M, Robertson LJ. The impact of globalisation on the distribution of Echinococcus multilocularis. Trends Parasitol. 2012;28:239–47.
Farrokh D, Zandi B, Rad MP, Tavakoli M. Hepatic alveolar echinococcosis. Arch Iran Med. 2015;18:199–202.
Sekiya C. Parasitic cirrhosis of the liver. Jpn J Clin Med. 1994;52:234–9.
Tappe D, Weise D, Ziegler U, Muller A, Mullges W, Stich A. Brain and lung metastasis of alveolar echinococcosis in a refugee from a hyperendemic area. J Med Microbiol. 2008;57:1420–3.
Takci E, Sengul G, Akar A, Uslu H, Alper F, Erdogan F, et al. Alveolar echinococcosis of the brain in five patients. J Clin Neurosci. 2008;15:1105–9.
Ohsaki Y, Sasaki T, Shibukawa K, Takahashi T, Osanai S. Radiological findings of alveolar hydatid disease of the lung caused by Echinococcus multilocularis. Respirology. 2007;12:458–61.
Wang H, Lu C, Liu X, Zhang W. Metastatic and prognostic factors in patients with alveolar Echinococcosis. Int J Clin Exp Pathol. 2015;8:11192–8.
Nunnari G, Pinzone MR, Gruttadauria S, Celesia BM, Madeddu G, Malaguarnera G, et al. Hepatic echinococcosis: clinical and therapeutic aspects. World J Gastroenterol. 2012;18:1448–58.
Piarroux M, Piarroux R, Giorgi R, Knapp J, Bardonnet K, Sudre B, et al. Clinical features and evolution of alveolar echinococcosis in France from 1982 to 2007: results of a survey in 387 patients. J Hepatol. 2011;55:1025–33.
Romig T, Kratzer W, Kimmig P, Frosch M, Gaus W, Flegel WA, et al. An epidemiologic survey of human alveolar echinococcosis in southwestern Germany. Am J Trop Med Hyg. 1999;61:566–73.
Kreidl P, Allerberger F, Judmaier G, Auer H, Aspock H, Hall AJ. Domestic pets as risk factors for alveolar hydatid disease in Austria. Am J Epidemiol. 1998;147:978–81.
Nahorski WL, Knap JP, Pawlowski ZS, Krawczyk M, Polanski J, Stefaniak J, et al. Human alveolar echinococcosis in Poland: 1990–2011. PLoS Negl Trop Dis. 2013;7:e1986.
Geramizadeh B, Baghernezhad M. Hepatic alveolar hydatid cyst: a brief review of published cases from Iran in the last 20 years. Hepat Mon. 2016;16:4.
Ito A, Budke CM. The present situation of echinococcoses in Mongolia. J Helminthol. 2015;89:680–8.
Zhang WB, Zhang ZZ, Wu WP, Shi BX, Li J, Zhou XN, et al. Epidemiology and control of echinococcosis in central Asia, with particular reference to the People’s Republic of China. Acta Trop. 2015;141:235–43.
Romig T, Dinkel A, Mackenstedt U. The present situation of echinococcosis in Europe. Parasitol Int. 2006;55:S187–91.
Liu WY, Delabrousse E, Blagosklonov O, Wang J, Zeng HC, Jiang Y, et al. Innovation in hepatic alveolar echinococcosis imaging: best use of old tools, and necessary evaluation of new ones. Parasite. 2014;21:14.
Qu B, Guo L, Sheng G, Yu F, Chen G, Wang YP, et al. Management of advanced hepatic alveolar echinococcosis: report of 42 cases. Am J Trop Med Hyg. 2017;96:680–5.
Amathieu R, Triba MN, Goossens C, Bouchemal N, Nahon P, Savarin P, et al. Nuclear magnetic resonance based metabolomics and liver diseases: recent advances and future clinical applications. World J Gastroenterol. 2016;22:417–26.
Kern P, Wen H, Sato N, Vuitton DA, Gruener B, Shao YM, et al. WHO classification of alveolar echinococcosis: principles and application. Parasitol Int. 2006;55:S283–7.
Xia JG, Psychogios N, Young N, Wishart DS. MetaboAnalyst: a web server for metabolomic data analysis and interpretation. Nucleic Acids Res. 2009;37:W652–60.
Xia JG, Sinelnikov IV, Han B, Wishart DS. MetaboAnalyst 3.0-making metabolomics more meaningful. Nucleic Acids Res. 2015;43:W251–7.
Anderson PE, Mahle DA, Doom TE, Reo NV, DelRaso NJ, Raymer ML. Dynamic adaptive binning: an improved quantification technique for NMR spectroscopic data. Metabolomics. 2011;7:179–90.
Dieterle F, Ross A, Schlotterbeck G, Senn H. Probabilistic quotient normalization as robust method to account for dilution of complex biological mixtures. Application in 1H NMR metabonomics. Anal Chem. 2006;78:4281–90.
Psychogios N, Hau DD, Peng J, Guo AC, Mandal R, Bouatra S, et al. The human serum metabolome. PLoS ONE. 2011;6:e16957.
Bouatra S, Aziat F, Mandal R, Guo AC, Wilson MR, Knox C, et al. The human urine metabolome. PLoS ONE. 2013;8:e73076.
Ogata H, Goto S, Sato K, Fujibuchi W, Bono H, Kanehisa M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999;27:29–34.
Wishart DS, Knox C, Guo AC, Eisner R, Young N, Gautam B, et al. HMDB: a knowledgebase for the human metabolome. Nucleic Acids Res. 2009;37:D603–10.
Triba MN, Le Moyec L, Amathieu R, Goossens C, Bouchemal N, Nahon P, et al. PLS/OPLS models in metabolomics: the impact of permutation of dataset rows on the K-fold cross-validation quality parameters. Mol Biosyst. 2015;11:13–9.
McLoughlin GA, Ma D, Tsang TM, Jones DNC, Cilia J, Hill MD, et al. Analyzing the effects of psychotropic drugs on metabolite profiles in rat brain using 1H NMR spectroscopy. J Proteome Res. 2009;8:1943–52.
Cui XQ, Churchill GA. Statistical tests for differential expression in cDNA microarray experiments. Genome Biol. 2003;4:210.
De Hertogh B, De Meulder B, Berger F, Pierre M, Bareke E, Gaigneaux A, et al. A benchmark for statistical microarray data analysis that preserves actual biological and technical variance. BMC Bioinform. 2010;11:17.
Hur M, Campbell AA, Almeida-de-Macedo M, Li L, Ransom N, Jose A, et al. A global approach to analysis and interpretation of metabolic data for plant natural product discovery. Nat Prod Rep. 2013;30:565–83.
Faul F, Erdfelder E, Buchner A, Lang AG. Statistical power analyses using G*Power 3.1: tests for correlation and regression analyses. Behav Res Methods. 2009;41:1149–60.
Faul F, Erdfelder E, Lang AG, Buchner A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods. 2007;39:175–91.
Xia JG, Broadhurst DI, Wilson M, Wishart DS. Translational biomarker discovery in clinical metabolomics: an introductory tutorial. Metabolomics. 2013;9:280–99.
Lin CG, Wei ZL, Cheng KK, Xu JJ, Shen GP, She C, et al. 1H NMR-based investigation of metabolic response to electro-acupuncture stimulation. Sci Rep. 2017;7:13.
Cohen J. A powr primer. Psychol Bull. 1992;112:155–9.
Dejong CHC, van de Poll MCG, Soeters PB, Jalan R, Damink S. Aromatic amino acid metabolism during liver failure. J Nutr. 2007;137:1579S–85S.
Abu-Omar MM, Loaiza A, Hontzeas N. Reaction mechanisms of mononuclear non-heme iron oxygenases. Chem Rev. 2005;105:2227–52.
Harris RA, Joshi M, Jeoung NH. Mechanisms responsible for regulation of branched-chain amino acid catabolism. Biochem Biophys Res Commun. 2004;313:391–6.
Zhang SH, Zeng XF, Ren M, Mao XB, Qiao SY. Novel metabolic and physiological functions of branched chain amino acids: a review. J Anim Sci Biotechnol. 2017;8:12.
Negro M, Giardina S, Marzani B, Marzatico F. Branched-chain amino acid supplementation does not enhance athletic performance but affects muscle recovery and the immune system. J Sport Med Phys Fit. 2008;48:347–51.
Allen PJ, Wise D, Greenway T, Khoo L, Griffin MJ, Jablonsky M. Using 1-D 1H and 2-D 1H J-resolved NMR metabolomics to understand the effects of anemia in channel catfish (Ictalurus punctatus). Metabolomics. 2015;11:1131–43.
Davies KM, Heaney RP, Rafferty K. Decline in muscle mass with age in women: a longitudinal study using an indirect measure. Metabolism. 2002;51:935–9.
Fernstrom JD. Large neutral amino acids: dietary effects on brain neurochemistry and function. Amino Acids. 2013;45:419–30.
Benitez-Diaz P, Miranda-Contreras L, Pena-Contreras Z, Davila-Vera D, Mendoza-Briceno RV, Palacios-Pru E. Histotypic mouse parietal cortex cultures: excitation/inhibition ratio and ultrastructural analysis. J Neurosci Methods. 2006;156:64–70.
Cryan JF, Kaupmann K. Don’t worry ‘B’ happy!: a role for GABAB receptors in anxiety and depression. Trends Pharmacol Sci. 2005;26:36–43.
Merali Z, Du LS, Hrdina P, Palkovits M, Faludi G, Poulter MO, et al. Dysregulation in the suicide brain: mRNA expression of corticotropin-releasing hormone receptors and GABAA receptorsubunits in frontal cortical brain region. J Neurosci. 2004;24:1478–85.
Fischer JE, Yoshimura N, Aguirre A, James JH, Cummings MG, Abel RM, et al. Plasma amino acids in patients with hepatic encephalopathy. Effects of amino acid infusions. Am J Surg. 1974;127:40–7.
Hagenfeldt L, Eriksson LS, Wahren J. Amino acids in liver disease. Proc Nutr Soc. 1983;42:497–506.
Benjamin IS, Engelbrecht GHC, Saunders SJ, Vanhoornhickman R. Amino acid imbalance following portal diversion in the rat. The relevance of nutrition and of hepatic function. J Hepatol. 1988;7:208–14.
Blonde-Cynober F, Aussel C, Cynober L. Abnormalities in branched-chain amino acid metabolism in cirrhosis: influence of hormonal and nutritional factors and directions for future research. Clin Nutr. 1999;18:5–13.
Cabrera G, Cabrejos ME, Morassutti AL, Cabezon C, Orellana J, Hellman U, et al. DNA damage, RAD9 and fertility/infertility of echinococcus granulosus hydatid cysts. J Cell Physiol. 2008;216:498–506.
Snow RJ, Murphy RM. Creatine and the creatine transporter: a review. Mol Cell Biochem. 2001;224:169–81.
Sestili P, Martinelli C, Bravi G, Piccoli G, Curci R, Battistelli M, et al. Creatine supplementation affords cytoprotection in oxidatively injured cultured mammalian cells via direct antioxidant activity. Free Radic Biol Med. 2006;40:837–49.
Cai KJ, Haris M, Singh A, Kogan F, Greenberg JH, Hariharan H, et al. Magnetic resonance imaging of glutamate. Nat Med. 2012;18:302–6.
Zhou JY, Tryggestad E, Wen ZB, Lal B, Zhou TT, Grossman R, et al. Differentiation between glioma and radiation necrosis using molecular magnetic resonance imaging of endogenous proteins and peptides. Nat Med. 2011;17:U130–308.
van Zijl PCM, Lam WW, Xu J, Knutsson L, Stanisz GJ. Magnetization transfer contrast and chemical exchange saturation transfer MRI. Features and analysis of the field-dependent saturation spectrum. Neuroimage. 2017;168:222–41.
This study was funded by the National Natural Science Foundation of China (grant nos. 81871445 and 81801788), and the Qinghai Basic Research Plan Project (grant no. 2018-ZJ-954Q).
Ethics approval and consent to participate
All procedures involving human subjects were approved by the Ethics Committee of the Affiliated Hospital of Qinghai University in Xining, Qinghai, China. Written informed consent was obtained from all participants involved in this study. All institutional and national guidelines for the care and use of biological samples were followed. The data acquired were in accordance with appropriate ethical requirements.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Figure S1.
Typical abdominal CT images of echinococcosis lesions in patients. Echinococcosis lesions are indicated by arrowheads. a Hepatic alveolar echinococcosis; CT image reveals dense calcification with irregular and indistinct margins in the lesion. b Hepatic cystic echinococcosis; CT image shows a sharply defined homogeneous cyst with multiple daughter cyst.
Additional file 2: Figure S2.
Predicted power profile with sample size per group (false discovery rate of 0.005).
Additional file 3: Figure S3.
The main window of the post-hoc power analysis specification in Gpower v.3.1 and the “effect size” drawer. Take as an example the metabolite of 1-methylhistidine which is different between control and HAE groups.
Additional file 4: Figure S4.
Univariate ROC curve analyses of metabolites in serum for discrimination of HAE patients from healthy individuals.
Additional file 5: Figure S5.
Univariate ROC curve analyses of metabolites in urine for discrimination of HAE patients from healthy individuals.
Additional file 6: Figure S6.
PLS-DA scores plot (left panel) and permutation test (right panel) of PLS-DA model comprising 21 identified characteristic metabolites.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
About this article
Cite this article
Lin, C., Chen, Z., Zhang, L. et al. Deciphering the metabolic perturbation in hepatic alveolar echinococcosis: a 1H NMR-based metabolomics study. Parasites Vectors 12, 300 (2019). https://doi.org/10.1186/s13071-019-3554-0