Next Article in Journal
Identification of the Main Metabolites of a Marine-Derived Strain of Penicillium brevicompactum Using LC and GC MS Techniques
Previous Article in Journal
Toward a Standardized Strategy of Clinical Metabolomics for the Advancement of Precision Medicine
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 10
Issue 2
10.3390/metabo10020053
metabolites-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Metabolomic Profiling of Wildtype and Transgenic Giardia lamblia Strains by 1H HR-MAS NMR Spectroscopy

by
Joachim Müller
1,*,
Martina Vermathen
2,
David Leitsch
3,
Peter Vermathen
4 and
Norbert Müller
1
1
Institute of Parasitology, Vetsuisse Faculty, University of Bern, Länggass-Strasse 122, CH-3012 Bern, Switzerland
2
Department of Chemistry and Biochemistry, University of Bern, Freiestrasse 3, CH-3012 Bern, Switzerland
3
Institute of Specific Prophylaxis and Tropical Medicine, Medical University of Vienna, Kinderspitalgasse 15, A-1090 Vienna, Austria
4
Departments of BioMedical Research and Radiology, University and Inselspital Bern, sitem-insel AG Freiburgstr. 3, CH-3010 Bern, Switzerland
*
Author to whom correspondence should be addressed.
Metabolites 2020, 10(2), 53; https://doi.org/10.3390/metabo10020053
Submission received: 10 December 2019 / Revised: 21 January 2020 / Accepted: 28 January 2020 / Published: 30 January 2020
Browse Figures
Versions Notes

Abstract

:
Giardia lamblia, a causative agent of persistent diarrhea in humans, domestic animals, and cattle, is usually treated with nitro compounds. Consequently, enzymes involved in anaerobic nitro reduction have been investigated in detail as potential targets. Their role within the normal metabolic context is, however, not understood. Using 1H high-resolution magic angle spinning (HR-MAS) NMR spectroscopy, we analyzed the metabolomes of G. lamblia trophozoites overexpressing three nitroreductases (NR1–NR3) and thioredoxin reductase (TrxR), most likely a scavenger of reactive oxygen species, as suggested by the results published in this study. We compared the patterns to convenient controls and to the situation in the nitro drug resistant strain C4 where NR1 is downregulated. We identified 27 metabolites in G. lamblia trophozoites. Excluding metabolites of high variability among different wildtype populations, only trophozoites overexpressing NR1 presented a distinct pattern of nine metabolites, in particular arginine catabolites, differing from the respective controls. This pattern matched a differential pattern between wildtype and strain C4. This suggests that NR1 interferes with arginine and thus energy metabolism. The exact metabolic function of NR1 (and the other nitroreductases) remains to be elucidated.

Graphical Abstract

1. Introduction

Giardia lamblia is an anaerobic or microaerophilic unicellular diplomonad causing persistent diarrhea in humans and various domestic animals. The life cycle comprises two stages, namely, the pear-shaped, flagellated trophozoite proliferating in the upper part of the small intestine, and the non-proliferative cyst representing the infectious stage of the parasite [1,2].
Analysis of the genome [3] has revealed prokaryote-like features such as short promoter sequences, a nearly complete absence of introns, and high similarities of some key enzymes of energy and intermediate metabolism, with prokaryotic enzymes most likely acquired by lateral transfer [4,5]. Thus, it is not surprising that the metabolism of G. lamblia trophozoites has unusual features in comparison to host cells. The anaerobic energy metabolism [6] is chemoheterotrophic yielding ethanol, acetic acid, and alanine as fermentation end products. Moreover, arginine is used as the main energy supply via arginine deiminase (ADI), ornithine carbamoyl transferase, and carbamate kinase, as detailed in earlier studies [7,8,9].
The treatment against giardiasis comprises nitro compounds such as metronidazole (MET), other 5-nitroimidazole compounds, or nitazoxanide (NTZ) [10,11]. Moreover, G. lamblia is susceptible to a variety of antibiotics because of its prokaryote-like transcription and translation machineries, as well as to albendazole [12]. According to a commonly accepted model, nitro compounds are activated by reduction yielding toxic intermediates causing nitrosative stress [13,14]. Pyruvate ferredoxin oxidoreductase (PFOR) is a main supplier of electrons [6,15]. This explains why knock-down of this enzyme is correlated with increased resistance to MET [16]. Studies with resistant strains have revealed, however, that resistance does not always correlate with reduced PFOR activity [17,18]. Other enzymes potentially involved in the reduction of nitro drugs are the nitroreductases (NR)1 [19,20] and NR2 [21,22]. Recently, a third homologue of NR1 and NR2 with weak nitroreductase activity, NR3, has been described [23]. All three NRs have higher quinone reductase than nitroreductase activities in vitro [22,23]. Overexpression of NR1 results in an increased susceptibility to nitro compounds in Escherichia coli [19,22,23] and in G. lamblia [19,23], whereas overexpression of NR2 results in an increased resistance in E. coli [21,22,23] in independent studies. Their function within the intermediary metabolism under physiological conditions (i.e., in the absence of nitro drugs) is completely unknown.
Thioredoxin-reductase (TrxR) has nitroreductase activity in vitro [15]. Overexpression of the wildtype allele in the G. lamblia strain WBC6 does not alter the susceptibility to NTZ and only slightly alters the susceptibility to MET, and overexpression of a dominant negative mutant allele has no effects on drug resistance [24]. In resistant strains of a different genetic background, however, TrxR is overexpressed [25]. Specific functions are not known, to date.
In order to get an insight into potential functions of these enzymes, the first step would be to compare metabolomic patterns of strains with altered expression levels of the respective genes. In a previous study, we introduced high-resolution magic angle spinning (HR-MAS) NMR spectroscopy as a suitable tool to study the metabolome of G. lamblia [26]. HR-MAS yields well-resolved liquid-like 1H-NMR spectra from semi-solid materials, reducing resonance line broadening effects caused by dipolar interactions and magnetic susceptibility differences inherent to most semi-solid materials [27]. Therefore, HR-MAS NMR has become a very attractive tool in the analysis of small metabolites from biological materials such as cells and tissue samples [28,29] or food samples [30,31,32]. Moreover, this method has proven to be suitable for the determination of biomarkers for the metabolic response to anticancer drugs [33,34,35].
Here, we present a metabolomic study using 1H HR-MAS NMR spectroscopy comparing transgenic G. lamblia strains overexpressing genes encoding for thioredoxin reductase and nitroreductases. Unfortunately, knock-down of the expression of the respective genes is not possible in our hands. An exception is TrxR, where a dominant negative mutant allele exists, the gene product of which forms an inactive complex with the wildtype TrxR. Another exception is the nitro drug resistant strain G. lamblia C4 derived from the wildtype WBC6 [17], which has significantly reduced NR1 protein levels when grown in the presence of NTZ [19,25].
In order to identify potential metabolomic patterns related to the overexpression of the respective genes, we eliminated metabolites from the discussion that showed strong differences between independent wildtype populations. Moreover, as we identified a distinct pattern in NR1 overexpressors only, we tested the observed pattern by a comparison between the wildtype WBC6 and the derived nitro drug-resistant strain C4. The strains used in this study are listed in Table 1.

2. Results and Discussion

2.1. Relative Amounts of Identified Metabolites

Using 1H-HR-MAS NMR, we identified 27 metabolites in G. lamblia WBC6 trophozoites, as listed in Table 2. The values of all individual integrals are available in the Supplementary Materials section (Table S1). In order to discriminate “natural” metabolite differences between wildtype populations from those related to transgenic modifications, we compared, in a first step, the metabolomes of two completely independent populations of WBC6 wildtype trophozoites (Table 2).
Nine of the metabolites identified, namely, glucose-1-P (Glc-1-P), seven amino acids (alanine, asparagine, cysteine, glutamine, leucine, threonine, valine), and trimethylamine (TMA) were significantly different between independent populations of WBC6 wildtype populations (Table 2, columns on the left). The variability of these nine metabolites may be due to rapid turn-over rates due to their pivotal roles in energy and intermediary metabolism. Furthermore, the variability may reflect the potential of the Giardia strain WBC6 to adapt to changing environmental conditions. Consequently, in order to identify patterns discriminating between strains, these metabolites had to be excluded. For an overview, principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (oPLS-DA) plots containing all Giardia strains given in Table 2 are displayed in Figure S1.

2.2. Phenotype and Metabolite Patterns Related to Thioredoxin Reductase Overexpression

Because overexpression of either allele of TrxR had little effect on nitro drug susceptibility, we investigated whether oxygen differentially affected the growth of strains overexpressing either the wildtype allele (TrxR) or a dominant negative allele (TrxR DN). There was no marked difference between the three strains tested under anaerobic growth conditions (Figure 1A). Under semi-aerobic conditions, however, the untransformed wildtype strain grew only in the presence of more than 10 mM D-cysteine (D-Cys) added to the medium. The strain overexpressing the dominant negative TrxR allele started growing already at 2.5 mM D-Cys. The strain overexpressing the wildtype TrxR grew even in the absence of D-Cys, reaching higher densities after 3 days than the two other strains (Figure 1B).
Our results suggest that TrxR acts as a scavenger of reactive oxygen species, as suggested by previous studies [5,36]. Because D-Cys is non-proteinogenic, we attribute the protective effect to a direct inactivation of oxygen, not to an indirect effect such as an increased synthesis of Cys-rich proteins. The TrxR DN strain may need less D-Cys than the wildtype trophozoites (WT) to resume growth because the protein itself may act as a scavenger independently of enzyme activities.
The comparison of metabolite patterns in G. lamblia trophozoites overexpressing either the wildtype thioredoxin reductase (TrxR) or a dominant negative allele (TrxR DN) revealed clear clustering of each group in principal component analysis (PCA; Figure 2A). Orthogonal partial least squares discriminant analysis (oPLS-DA) showed complete and significant separation between the groups along the first, predictive PLS component latent variable (LV) 1 (Figure 2B).
A detailed comparison of single metabolites revealed that these differences could be attributed to acetate, which was significantly increased, and lysine, which was significantly decreased in TrxR DN as compared to TrxR (Table 2, columns 3 and 4). Other separating contributions only derived from metabolites with pool sizes varying between wildtype populations. Because acetate is an indirect product of PFOR, overexpressors of the wildtype allele may reduce it more efficiently to ethanol than the overexpressors of the dominant negative mutant allele. As the catabolism of lysine yields acetyl-coenzymeA (and thus acetate) as a final product (see [37] and references therein), the increase of acetate may also be directly linked to an increased catabolism of lysine, thereby causing the observed changes in the pool sizes of both metabolites.

2.3. Metabolite Patterns Related to Nitroreductase Overexpression

When comparing trophozoites overexpressing the nitro reductases NR1, NR2, and NR3 to a control strain overexpressing glucuronidase A (see Table 1 for phenotypes), PCA resulted in complete separation of all classes, and oPLS-DA revealed significant differences between all these strains (Figure 3). More specifically, the NR3 and GusA groups were clustered closer together, whereas the NR1 group was the furthest away and was already distinguished in PCA from all other groups solely by its positive scores along principal component (PC) 1 (Figure 3A).
Upon exclusion of metabolites differing in wildtype population, the most striking differences were observed in NR1 overexpressing trophozoites, with nine metabolites significantly changed as compared to the GusA control, namely, acetate, glutamine, phenylalanine, and tyrosine with significantly lower levels, and cystine, methionine, ornithine, citrulline, and pipecolic acid with significantly higher levels in NR1 overexpressors as compared to GusA control. Acetate was significantly decreased in all NR overexpressors as compared to the GusA control (Table 2, columns on the right). In the strain overexpressing NR2, besides acetate, only ornithine and cysteine followed the same pattern as in the NR1 overexpressor. In the strain overexpressing NR3, only acetate was significantly affected.
In particular, amino acids generated during the catabolism of the basic amino acids arginine (ornithine and citrulline) and lysine (pipecolic acid) were significantly increased in NR1 expressing trophozoites as compared to the GusA control (Table 2).
A typical NMR spectrum representing peaks for ornithine and citrulline from NR1 and GusA trophozoites is shown in Figure 4.
Thus, taken together, we could identify a distinct pattern of nine metabolites differential between NR1 and GusA control trophozoites.

2.4. Comparison of the NR1 Pattern with the Situation in Nitro Drug-Resistant Trophozoites

These findings prompted us to test whether the pattern obtained by comparing NR1 overexpressing trophozoites (with an increased susceptibility to NTZ) to GusA control trophozoites could be linked to the nitro drug susceptibility. For this, we compared the metabolomes of WBC6 wildtype trophozoites with trophozoites of the nitro drug-resistant strain C4 with reduced NR1 activity. Individual normalized peak integrals and a comparison of all metabolites between WBC6 and C4 trophozoites are given in the Supplementary Materials section (Table S3 and Figure S2).
Subsequently, an oPLS-DA model calculated for discrimination between NR1 and GusA (Figure 5A) was applied to the nitro drug-resistant C4 trophozoites and wildtype trophozoites (WT), as shown in Figure 5B. The predictive value of the model led to complete separation of C4 versus WT, that is, suggesting a correlation between nitroreductase activity and the metabolome. In a next step, we plotted the relative changes observed in NR1 compared to GusA of the nine differential metabolites listed above versus the corresponding relative changes obtained in WT compared to C4 (Figure 5C).
As shown in Figure 5C, the pattern of the nine identified metabolites inversely correlated with a pattern of the same metabolites in the nitro drug-resistant strain C4, which has lower NR1 levels than its corresponding wildtype. The correlation was highly significant.
The catabolism of arginine, a major energy source for G. lamblia [8,26], is of particular interest within this context. Citrulline and ornithine, the first and the second catabolites of arginine, are positively correlated with NR1 levels, whereas glutamine, an indirect catabolite of ornithine [38], is negatively correlated with NR1 levels. This suggests that NR1 interferes with the catabolism of arginine. Because the ATP levels are unaltered in NR1, as well as in all other transgenic trophozoites, a direct interference with energy supply is unlikely. In G. lamblia, arginine catabolism is not only a major source for ATP, but has also other functions. Within this context, arginine deiminase (ADI) is of particular interest. ADI does not only transform free arginine into citrulline [39], but also arginine residues in a variety of proteins including variant surface proteins [40]. Accumulation of inhibitory amounts of citrulline of this enzyme may thus affect regulatory pathways depending on ADI, a “moonlighting” protein, as evidenced by recent studies [41,42]. Because NR1 overexpressed and purified from E. coli has a much higher quinone than nitroreductase activity in vitro [22,23], it is possible that it interferes with susceptibility/resistance to nitro compounds, not (only) by direct reduction of these compounds, but rather indirectly via reduction of an unknown quinone substrate or by interfering with coenzymes involved in electron transfer [18]. The exact function of this enzyme remains to be elucidated.

3. Conclusions

In contrast to other organisms, the metabolome plasticity of the early diverging eukaryont G. lamblia is far from being elucidated, especially in response to genetic transformation and/or drug-based selection. The present study could therefore be regarded as a paradigm extendable to studies on other extracellular parasites such as, for example, amebae [2], trichomonads, or even helminths [43]. These investigations may yield novel antiparasitic drug targets.

4. Materials and Methods

4.1. Biochemicals

If not otherwise stated, all biochemical reagents were from Sigma (St Louis, MO, USA). D2O-based 50 mM phosphate-buffered saline (PBS) was prepared from mixing adequate amounts of 50 mM K2HPO4 and NaH2PO4 solutions in D2O containing 0.9% NaCl to give a pH of 7.0 (corresponding to a pD of 7.4).

4.2. Growth of Thioredoxin Reductase Overexpressors

G. lamblia trophozoites overexpressing a wildtype allele of thioredoxin reductase (TrxR), a dominant negative mutant allele (TrxR DN; see [24]), or wildtype trophozoites were seeded into 96-well plates (103 trophozoites per well) containing normal culture medium or medium supplemented with increasing amounts of D-cysteine. D-cysteine was used instead of L-cysteine because it is not proteinogenic. The plates were incubated under anaerobic (85% N2, 10% H2, 5% CO2) or semi-aerobic (85% N2, 10% CO2, 5% O2) conditions. After 3 days, cells were quantified using the Alamar blue assay, as described elsewhere [44,45]. The assays were run in quadruplicate.

4.3. Preparation of Cell Samples for HR-MAS NMR

G. lamblia trophozoites from the strains listed in Table 1 were grown to confluence in modified TYI-S33 medium and harvested by chilling on ice, as previously described [18,23,24]. After centrifugation, trophozoites were washed once in ice-cold PBS, followed by two washes with ice-cold D2O-based PBS, counted in a Neubauer chamber, and shock frozen in liquid nitrogen. All samples were stored in aliquots corresponding to 3 × 106 cells at −70 °C until further analysis [26]. The nitroreductase overexpressors were compared to a GusA overexpressor and not directly to an untransfected wildtype strain because all strains had to undergo a selection procedure of transfectants, which was shown to influence expression of various genes [46].

4.4. HR-MAS NMR Spectroscopy and Processing of Spectra

1H HR-MAS NMR experiments of cell samples were performed on a 500 MHz Bruker Avance II spectrometer (Bruker BioSpin, Fällanden, Switzerland) with a 4 mm HR-MAS dual inverse 1H/13C probe equipped with a z-gradient directed along the magic angle axis. Each sample was thawed during 2 min immediately before HR-MAS NMR measurement. Suspensions corresponding to 3 × 106 cells were filled into ZrO2-HR-MAS rotors with PTFE (Teflon) inserts providing a sample volume of 50 μL. The sample amount inside the rotor was determined by weight control. The rotor was spun at a MAS rate of 3 kHz during measurements. The temperature was set to 279 K (nominal), allowing the sample to equilibrate during 5 min. 1H spectra were recorded using a 1D PROJECT (periodic refocusing of J evolution by coherence transfer) pulse sequence with water pre-saturation and a rotor-synchronized T2 filter of 115 ms to suppress broad components with short T2 relaxation times. Each spectrum was acquired with 1024 scans, a spectral width of 6010 Hz, a data size of 16 k points, an acquisition time of 1.36 s, and a relaxation delay of 4 s. 1H NMR resonance assignments were based on the spectral analysis previously published [26]. To confirm resonance assignments, additional 2D 1H1H-TOCSY spectra were acquired for at least one sample within each sample group using the DIPSI2 sequence (“dipsi2phpr” from the Bruker pulse program library) with water pre-saturation during relaxation delay.
Data acquisition and spectral processing were performed with the Bruker TopSpin v 3.2 Pl5 software. All spectra were processed identically—the summed free induction decays (FIDs) were exponentially multiplied with a line broadening factor of 1 Hz and Fourier transformed. The resulting spectra were phased and calibrated to the γ-CH3 resonance of threonine (Thr) at 1.326 ppm. Baseline correction for the spectral region of 0–5.7 ppm was performed by applying a cubic spline function with individually set baseline points. The spectral region between 5.7 and 10 ppm was baseline-corrected using the underground removal tool of the Bruker AMIX 3.9.14 software applying a filter of 30 Hz.

4.5. Bioinformatics and Statistical Methods

Statistical analysis of the inhibition tests was performed with suitable tools from the open source software package R [47]. The results were compared to the controls separately for each strain and analysed by two-sided paired t-tests. The p-values were corrected for multiple comparisons.
Multivariate statistical data analysis was performed using MATLAB R2012a (Mathworks) and PLS-Toolbox 7.5.2 (Eigenvector Research, Inc., Manson, WA, USA). The processed NMR spectra were subdivided into 138 individually sized buckets (integral regions) to cover single resonances as much as possible in the spectral range between 0 and 10 ppm. Bucket selection was guided by an overlay of summed spectra from each sample group. The spectral region of the residual water resonance (4.66–5.2 ppm) and noise regions were excluded. Both bucket selection and integration of bucket regions were performed with the Bruker TopSpin v 3.2 Pl5 software. The numeric bucket integrals were exported into the Excel (Microsoft Office Professional Plus 2016) software, and the buckets were assigned as far as possible to contributing metabolites. The subsequent steps were performed using the PLS-Toolbox software. Probabilistic quotient normalization (PQN) was applied to the integral regions to account for different sample amounts [48]. Mean centering was applied to the buckets, as well as pareto scaling to reduce the relative importance of large values [49]. The samples (spectra) were grouped according to the corresponding G. lamblia strain, as given in Table 2, and the number of samples per group was n = 4–10 (Table 2). Principal component analysis (PCA) was applied to probe for clustering. The first two principal components (PC 1, PC 2) were plotted to generate the PC score plots to visualize unsupervised variation of the datasets within each group. For testing, if the groups can be distinguished on the basis of their metabolite spectra, orthogonal partial least squares discriminant analysis (oPLS-DA) was performed using sample group assignment as y-variables. One to three latent variables (LV) were calculated, and LV1 and LV2 were plotted to obtain the oPLS-score plots. When only one LV was calculated (as suggested by the software in order not to overfit the data), this value was plotted against the Q residuals. For both PCA and oPLS, 83.4% confidence ellipses, indicating a statistical significant separation of p < 0.05 when not overlapping, were superimposed on the score plots to illustrate the goodness of class separation [50]. The venetian blinds method was used for cross validation of the oPLS models. Parameters of the oPLS including the goodness of fit (R2) and the predictive ability (Q2) are summarized in Table S2.
For univariate statistical analysis, integrals of representative spectral regions for each of the 27 metabolites were evaluated. These integrals were selected such that either single or averages of multiple isolated resonances from each metabolite were included, minimizing contributions from other metabolite resonances as much as possible. The selection of integrals is visualized in Figure S3 on representative spectra. The resonance integrals or corresponding buckets were extracted from the PQN-normalized data matrix that was created for multivariate statistics, as described above. For each integral, the mean value and standard deviation was calculated using the Excel software. Single factor ANOVA and t-test (Excel) were used to probe for metabolite differences between the G. lamblia strains. All p-values were multiplied with a factor of 27 to correct for multiple comparisons (Bonferroni correction).

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-1989/10/2/53/s1, Figure S1: PCA and oPLS-DA of all integral regions (138 buckets) obtained by 1H-HR-MAS NMR analysis of TrxR, TrxR DN, WT1, NR1, NR2, NR3, GusA, and WT2 trophozoites. The strains are detailed in Table 1. A, PCA; B, oPLS-DA, Figure S2: Metabolites identified in G. lamblia WBC6 wildtype trophozoites (WT, n = 4) and in trophozoites of the nitro drug-resistant strain C4 (n = 10); mean values ± SD. Metabolites marked with ** and * were significantly different with p < 0.001 and p < 0.005 respectively according to t-test and correction for multiple comparisons. Significant levels in gray were excluded from discussion due to high variability between different WT-batches, Figure S3. Representative 1H HR-MAS NMR spectra of Giardia trophozoites indicating selection of spectral regions / buckets and their metabolite assignments that were used for integration in univariate analysis. Ac (acetate) and Ile are better visible in the upper spectrum (A) than in the lower spectrum (A’). For metabolites marked with * the average of multiple regions as indicated were used. A-D, contiguous spectral regions; Cyn, cysteine, Table S1: Peak integrals obtained by 1H-HR-MAS NMR analysis of G. lamblia WBC6 wildtype trophozoites (WT; two independent experiments), in trophozoites overexpressing thioredoxin reductase (TrxR) or a dominant negative mutation of TrxR (TrxR DN) as a control, and in trophozoites overexpressing the nitroreductases NR1, NR2, NR3, or GusA as a control, Table S2: Summary of parameters and statistical results of the oPLS-DA shown in Figure 2, Figure 3 and Figure 5. Table S3: Peak integrals obtained by 1H-HR-MAS NMR analysis of G. lamblia WBC6 wildtype and nitro drug-resistant C4 trophozoites.

Author Contributions

Conceptualization, J.M. and N.M.; data curation, P.V.; formal analysis, M.V., D.L., and P.V.; funding acquisition, N.M.; investigation, N.M.; methodology, M.V., D.L., and P.V.; project administration, N.M.; resources, P.V.; software, M.V. and P.V.; writing—original draft, J.M.; writing—review and editing, J.M., M.V., D.L., P.V., and N.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Swiss National Science Foundation (grant no. 31003A_163230 and grant no. 206021-128736 for the purchase of the HR-MAS probe) and by the Austrian Science Fund (FWF; grant J3492).

Acknowledgments

The authors thank Vreni Balmer (Institute of Parasitology, Berne) for technical assistance. Open Access Funding by the Austrian Science Fund (FWF).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cernikova, L.; Faso, C.; Hehl, A.B. Five facts about Giardia lamblia. PLoS Pathog. 2018, 14, e1007250. [Google Scholar] [CrossRef] [Green Version]
  2. Hemphill, A.; Müller, N.; Müller, J. Comparative pathobiology of the intestinal protozoan parasites Giardia lamblia, Entamoeba histolytica, and Cryptosporidium parvum. Pathogens 2019, 8, 116. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Morrison, H.G.; McArthur, A.G.; Gillin, F.D.; Aley, S.B.; Adam, R.D.; Olsen, G.J.; Best, A.A.; Cande, W.Z.; Chen, F.; Cipriano, M.J.; et al. Genomic minimalism in the early diverging intestinal parasite Giardia lamblia. Science 2007, 317, 1921–1926. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Nixon, J.E.; Wang, A.; Field, J.; Morrison, H.G.; McArthur, A.G.; Sogin, M.L.; Loftus, B.J.; Samuelson, J. Evidence for lateral transfer of genes encoding ferredoxins, nitroreductases, NADH oxidase, and alcohol dehydrogenase 3 from anaerobic prokaryotes to Giardia lamblia and Entamoeba histolytica. Eukaryot. Cell 2002, 1, 181–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ma’ayeh, S.Y.; Knorr, L.; Svard, S.G. Transcriptional profiling of Giardia intestinalis in response to oxidative stress. Int. J. Parasitol. 2015, 45, 925–938. [Google Scholar] [CrossRef]
  6. Brown, D.M.; Upcroft, J.A.; Edwards, M.R.; Upcroft, P. Anaerobic bacterial metabolism in the ancient eukaryote Giardia duodenalis. Int. J. Parasitol. 1998, 28, 149–164. [Google Scholar] [CrossRef]
  7. Schofield, P.J.; Costello, M.; Edwards, M.R.; O’Sullivan, W.J. The arginine dihydrolase pathway is present in Giardia intestinalis. Int. J. Parasitol. 1990, 20, 697–699. [Google Scholar] [CrossRef]
  8. Edwards, M.R.; Schofield, P.J.; O’Sullivan, W.J.; Costello, M. Arginine metabolism during culture of Giardia intestinalis. Mol. Biochem. Parasitol. 1992, 53, 97–103. [Google Scholar] [CrossRef]
  9. Schofield, P.J.; Edwards, M.R.; Matthews, J.; Wilson, J.R. The pathway of arginine catabolism in Giardia intestinalis. Mol. Biochem. Parasitol. 1992, 51, 29–36. [Google Scholar] [CrossRef]
  10. Nash, T.E. Treatment of Giardia lamblia infections. Pediatr. Infect. Dis. J. 2001, 20, 193–195. [Google Scholar] [CrossRef]
  11. Leitsch, D. A review on metronidazole: An old warhorse in antimicrobial chemotherapy. Parasitology 2017, 146, 1167–1178. [Google Scholar] [CrossRef] [PubMed]
  12. Müller, J.; Hemphill, A. New approaches for the identification of drug targets in protozoan parasites. Int Rev. Cell Mol. Biol. 2013, 301, 359–401. [Google Scholar] [CrossRef] [PubMed]
  13. Lloyd, D.; Harris, J.C.; Maroulis, S.; Mitchell, A.; Hughes, M.N.; Wadley, R.B.; Edwards, M.R. Nitrosative stress induced cytotoxicity in Giardia intestinalis. J. Appl. Microbiol. 2003, 95, 576–583. [Google Scholar] [CrossRef] [PubMed]
  14. Müller, N.; Müller, J. Giardia . In Molecular Parasitology; Walochnik, J., Duchêne, M., Eds.; Springer: Vienna, Austria, 2016; pp. 93–114. [Google Scholar]
  15. Leitsch, D.; Burgess, A.G.; Dunn, L.A.; Krauer, K.G.; Tan, K.; Duchêne, M.; Upcroft, P.; Eckmann, L.; Upcroft, J.A. Pyruvate:ferredoxin oxidoreductase and thioredoxin reductase are involved in 5-nitroimidazole activation while flavin metabolism is linked to 5-nitroimidazole resistance in Giardia lamblia. J. Antimicrob. Chemother. 2011, 66, 1756–1765. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Dan, M.; Wang, A.L.; Wang, C.C. Inhibition of pyruvate-ferredoxin oxidoreductase gene expression in Giardia lamblia by a virus-mediated hammerhead ribozyme. Mol. Microbiol. 2000, 36, 447–456. [Google Scholar] [CrossRef]
  17. Müller, J.; Sterk, M.; Hemphill, A.; Müller, N. Characterization of Giardia lamblia WB C6 clones resistant to nitazoxanide and to metronidazole. J. Antimicrob. Chemother. 2007, 60, 280–287. [Google Scholar] [CrossRef] [Green Version]
  18. Müller, J.; Hemphill, A.; Müller, N. Physiological aspects of nitro drug resistance in Giardia lamblia. Int, J. Parasitol. Drugs Drug Resist. 2018, 8, 271–277. [Google Scholar] [CrossRef]
  19. Nillius, D.; Müller, J.; Müller, N. Nitroreductase (GlNR1) increases susceptibility of Giardia lamblia and Escherichia coli to nitro drugs. J. Antimicrob. Chemother. 2011, 66, 1029–1035. [Google Scholar] [CrossRef] [Green Version]
  20. Ansell, B.R.; Baker, L.; Emery, S.J.; McConville, M.J.; Svard, S.G.; Gasser, R.B.; Jex, A.R. Transcriptomics indicates active and passive metronidazole resistance mechanisms in three seminal Giardia lines. Front. Microbiol. 2017, 8, 398. [Google Scholar] [CrossRef] [Green Version]
  21. Müller, J.; Schildknecht, P.; Müller, N. Metabolism of nitro drugs metronidazole and nitazoxanide in Giardia lamblia: Characterization of a novel nitroreductase (GlNR2). J. Antimicrob. Chemother. 2013, 68, 1781–1789. [Google Scholar] [CrossRef] [Green Version]
  22. Müller, J.; Rout, S.; Leitsch, D.; Vaithilingam, J.; Hehl, A.; Müller, N. Comparative characterisation of two nitroreductases from Giardia lamblia as potential activators of nitro compounds. Int. J. Parasitol. Drugs Drug Resist. 2015, 5, 37–43. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Müller, J.; Müller, N. Nitroreductases of bacterial origin in Giardia lamblia: Potential role in detoxification of xenobiotics. MicrobiologyOpen 2019, 8, e904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Leitsch, D.; Müller, J.; Müller, N. Evaluation of Giardia lamblia thioredoxin reductase as drug activating enzyme and as drug target. Int. J. Parasitol. Drugs Drug Resist. 2016, 6, 148–153. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Müller, J.; Braga, S.; Heller, M.; Müller, N. Resistance formation to nitro drugs in Giardia lamblia: No common markers identified by comparative proteomics. Int. J. Parasitol. Drugs Drug Resist. 2019, 9, 112–119. [Google Scholar] [CrossRef]
  26. Vermathen, M.; Müller, J.; Furrer, J.; Müller, N.; Vermathen, P. 1H HR-MAS NMR spectroscopy to study the metabolome of the protozoan parasite Giardia lamblia. Talanta 2018, 188, 429–441. [Google Scholar] [CrossRef] [PubMed]
  27. Power, W.P. High-Resolution Magic Angle Spinning-enabling applications of NMR spectroscopy to semi-solid phases. Annu Rep. NMR Spectro. 2010, 72, 111–156. [Google Scholar] [CrossRef]
  28. Dietz, C.; Ehret, F.; Palmas, F.; Vandergrift, L.A.; Jiang, Y.; Schmitt, V.; Dufner, V.; Habbel, P.; Nowak, J.; Cheng, L.L. Applications of high-resolution magic angle spinning MRS in biomedical studies II-Human diseases. NMR Biomed. 2017, 30, e3784. [Google Scholar] [CrossRef]
  29. Kaebisch, E.; Fuss, T.L.; Vandergrift, L.A.; Toews, K.; Habbel, P.; Cheng, L.L. Applications of high-resolution magic angle spinning MRS in biomedical studies I-cell line and animal models. NMR Biomed. 2017, 30, e3700. [Google Scholar] [CrossRef]
  30. Santos, A.D.C.; Fonseca, F.A.; Liao, L.M.; Alcantara, G.B.; Barison, A. High-resolution magic angle spinning nuclear magnetic resonance in foodstuff analysis. Trac-Trend Anal. Chem. 2015, 73, 10–18. [Google Scholar] [CrossRef]
  31. Vermathen, M.; Marzorati, M.; Baumgartner, D.; Good, C.; Vermathen, P. Investigation of different apple cultivars by High Resolution Magic Angle Spinning NMR. A feasibility study. J. Agr. Food Chem. 2011, 59, 12784–12793. [Google Scholar] [CrossRef]
  32. Vermathen, M.; Marzorati, M.; Diserens, G.; Baumgartner, D.; Good, C.; Gasser, F.; Vermathen, P. Metabolic profiling of apples from different production systems before and after controlled atmosphere (CA) storage studied by 1H high resolution-magic angle spinning (HR-MAS) NMR. Food Chem. 2017, 233, 391–400. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Duarte, I.F.; Ladeirinha, A.F.; Lamego, I.; Gil, A.M.; Carvalho, L.; Carreira, I.M.; Melo, J.B. Potential markers of cisplatin treatment response unveiled by NMR metabolomics of human lung cells. Mol. Pharm. 2013, 10, 4242–4251. [Google Scholar] [CrossRef] [PubMed]
  34. Lamego, I.; Duarte, I.F.; Marques, M.P.; Gil, A.M. Metabolic markers of MG-63 osteosarcoma cell line response to doxorubicin and methotrexate treatment: Comparison to cisplatin. J. Proteome Res. 2014, 13, 6033–6045. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Vermathen, M.; Paul, L.E.; Diserens, G.; Vermathen, P.; Furrer, J. 1H HR-MAS NMR Based Metabolic Profiling of Cells in Response to Treatment with a Hexacationic Ruthenium Metallaprism as Potential Anticancer Drug. PLoS ONE 2015, 10, e0128478. [Google Scholar] [CrossRef] [Green Version]
  36. Brogi, S.; Fiorillo, A.; Chemi, G.; Butini, S.; Lalle, M.; Ilari, A.; Gemma, S.; Campiani, G. Structural characterization of Giardia duodenalis thioredoxin reductase (gTrxR) and computational analysis of its interaction with NBDHEX. Eur. J. Med. Chem. 2017, 135, 479–490. [Google Scholar] [CrossRef]
  37. Takei, M.; Ando, Y.; Saitoh, W.; Tanimoto, T.; Kiyosawa, N.; Manabe, S.; Sanbuissho, A.; Okazaki, O.; Iwabuchi, H.; Yamoto, T.; et al. Ethylene glycol monomethyl ether-induced toxicity is mediated through the inhibition of flavoprotein dehydrogenase enzyme family. Toxicol. Sci. 2010, 118, 643–652. [Google Scholar] [CrossRef] [Green Version]
  38. Morris, S.M., Jr. Recent advances in arginine metabolism. Curr. Opin. Clin. Nutr. Metab. Care. 2004, 7, 45–51. [Google Scholar] [CrossRef]
  39. Knodler, L.A.; Sekyere, E.O.; Stewart, T.S.; Schofield, P.J.; Edwards, M.R. Cloning and expression of a prokaryotic enzyme, arginine deiminase, from a primitive eukaryote Giardia intestinalis. J. Biol. Chem. 1998, 273, 4470–4477. [Google Scholar] [CrossRef] [Green Version]
  40. Touz, M.C.; Rópolo, A.S.; Rivero, M.R.; Vranych, C.V.; Conrad, J.T.; Svard, S.G.; Nash, T.E. Arginine deiminase has multiple regulatory roles in the biology of Giardia lamblia. J. Cell Sci. 2008, 121, 2930–2938. [Google Scholar] [CrossRef] [Green Version]
  41. Ropolo, A.S.; Feliziani, C.; Touz, M.C. Unusual proteins in Giardia duodenalis and their role in survival. Adv. Parasitol. 2019, 106, 1–50. [Google Scholar] [CrossRef]
  42. Ortega-Pierres, M.G.; Arguello-Garcia, R. Giardia duodenalis: Role of secreted molecules as virulent factors in the cytotoxic effect on epithelial cells. Adv. Parasitol. 2019, 106, 129–169. [Google Scholar] [CrossRef] [PubMed]
  43. Bringaud, F.; Ebikeme, C.; Boshart, M. Acetate and succinate production in amoebae, helminths, diplomonads, trichomonads and trypanosomatids: Common and diverse metabolic strategies used by parasitic lower eukaryotes. Parasitology 2010, 137, 1315–1331. [Google Scholar] [CrossRef] [PubMed]
  44. Müller, J.; Nillius, D.; Hehl, A.; Hemphill, A.; Müller, N. Stable expression of Escherichia coli beta-glucuronidase A (GusA) in Giardia lamblia: Application to high-throughput drug susceptibility testing. J. Antimicrob. Chemother 2009, 64, 1187–1191. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Bénéré, E.; da Luz, R.A.; Vermeersch, M.; Cos, P.; Maes, L. A new quantitative in vitro microculture method for Giardia duodenalis trophozoites. J. Microbiol. Methods 2007, 71, 101–106. [Google Scholar] [CrossRef] [PubMed]
  46. Su, L.H.; Lee, G.A.; Huang, Y.C.; Chen, Y.H.; Sun, C.H. Neomycin and puromycin affect gene expression in Giardia lamblia stable transfection. Mol. Biochem Parasitol 2007, 156, 124–135. [Google Scholar] [CrossRef]
  47. R_Core_Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2012. [Google Scholar]
  48. 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–4290. [Google Scholar] [CrossRef]
  49. Van den Berg, R.A.; Hoefsloot, H.C.; Westerhuis, J.A.; Smilde, A.K.; van der Werf, M.J. Centering, scaling, and transformations: Improving the biological information content of metabolomics data. BMC Genom. 2006, 7, 142. [Google Scholar] [CrossRef] [Green Version]
  50. Knol, M.J.; Pestman, W.R.; Grobbee, D.E. The (mis)use of overlap of confidence intervals to assess effect modification. Eur. J. Epidemiol 2011, 26, 253–254. [Google Scholar] [CrossRef] [Green Version]
Metabolites 10 00053 g001 550
Figure 1. Thioredoxin reductase protected Giardia lamblia trophozoites from oxygen. Growth of G. lamblia wildtype trophozoites (WT) and of trophozoites overexpressing either wildtype thioredoxin reductase (TrxR) or a dominant negative mutant allele (TrxR DN) under anaerobic (A) or semi-aerobic (B) conditions in medium supplemented with different amounts of D-cysteine (D-Cys).
Figure 1. Thioredoxin reductase protected Giardia lamblia trophozoites from oxygen. Growth of G. lamblia wildtype trophozoites (WT) and of trophozoites overexpressing either wildtype thioredoxin reductase (TrxR) or a dominant negative mutant allele (TrxR DN) under anaerobic (A) or semi-aerobic (B) conditions in medium supplemented with different amounts of D-cysteine (D-Cys).
Metabolites 10 00053 g001
Metabolites 10 00053 g002 550
Figure 2. Principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (oPLS-DA) of all integral regions (138 buckets) obtained by 1H high-resolution magic angle spinning (HR-MAS) NMR analysis of TrxR and TrxR DN trophozoites. The ellipses show 83.4% confidence intervals. The strains are detailed in Table 1. (A) PCA; (B) oPLS-DA. R2 and Q2 values for the oPLS model are given in Table S2. PC: principal component; LV: latent variable.
Figure 2. Principal component analysis (PCA) and orthogonal partial least squares discriminant analysis (oPLS-DA) of all integral regions (138 buckets) obtained by 1H high-resolution magic angle spinning (HR-MAS) NMR analysis of TrxR and TrxR DN trophozoites. The ellipses show 83.4% confidence intervals. The strains are detailed in Table 1. (A) PCA; (B) oPLS-DA. R2 and Q2 values for the oPLS model are given in Table S2. PC: principal component; LV: latent variable.
Metabolites 10 00053 g002
Metabolites 10 00053 g003 550
Figure 3. PCA and oPLS-DA of all integral regions (138 buckets) obtained by 1H-HR-MAS NMR analysis of GusA, NR1, NR2, and NR3 trophozoites. The ellipses show 83.4% confidence intervals. R2 and Q2 values for the oPLS model are given in Table S2. The strains are detailed in Table 1. (A) PCA; (B) oPLS-DA.
Figure 3. PCA and oPLS-DA of all integral regions (138 buckets) obtained by 1H-HR-MAS NMR analysis of GusA, NR1, NR2, and NR3 trophozoites. The ellipses show 83.4% confidence intervals. R2 and Q2 values for the oPLS model are given in Table S2. The strains are detailed in Table 1. (A) PCA; (B) oPLS-DA.
Metabolites 10 00053 g003
Metabolites 10 00053 g004 550
Figure 4. Representative 1H HR-MAS NMR spectra of suspensions of NR1 and GusA trophozoites in phosphate buffered saline (PBS). The resonances of citrulline (Ctl) and Ornithine (Orn) between 3.0 and 3.2 ppm are highlighted by the gray box. (A) Spectral region from 1.5 ppm to 4 ppm, (B) highlighted region in detail.
Figure 4. Representative 1H HR-MAS NMR spectra of suspensions of NR1 and GusA trophozoites in phosphate buffered saline (PBS). The resonances of citrulline (Ctl) and Ornithine (Orn) between 3.0 and 3.2 ppm are highlighted by the gray box. (A) Spectral region from 1.5 ppm to 4 ppm, (B) highlighted region in detail.
Metabolites 10 00053 g004
Metabolites 10 00053 g005 550
Figure 5. The metabolic pattern obtained by comparing NR1 versus GusA can be linked to the nitro drug susceptibility. (A) oPLS-DA score plots displaying the model calculated for distinguishing NR1 from GusA along the latent variable LV 1. (B) Class prediction based on the model shown in (A) applied to C4 and WT strains. The ellipses show 83.4% confidence intervals. (C) Ratios of the peak integrals of the eight differential metabolites listed above from NR1/GusA vs. WT/C4. The correlation between metabolite ratios in G. lamblia NR1 vs. GusA and WT vs. C4 trophozoites (see Table 1) was highly significant (Pearsons correlation coefficient 0.942; p < 0.0005). R2 and Q2 values for the oPLS model are given in Table S2.
Figure 5. The metabolic pattern obtained by comparing NR1 versus GusA can be linked to the nitro drug susceptibility. (A) oPLS-DA score plots displaying the model calculated for distinguishing NR1 from GusA along the latent variable LV 1. (B) Class prediction based on the model shown in (A) applied to C4 and WT strains. The ellipses show 83.4% confidence intervals. (C) Ratios of the peak integrals of the eight differential metabolites listed above from NR1/GusA vs. WT/C4. The correlation between metabolite ratios in G. lamblia NR1 vs. GusA and WT vs. C4 trophozoites (see Table 1) was highly significant (Pearsons correlation coefficient 0.942; p < 0.0005). R2 and Q2 values for the oPLS model are given in Table S2.
Metabolites 10 00053 g005
Table
Table 1. List of the transgenic Giardia lamblia strains used in this study. All transgenes were fused to the strong, constitutive arginine deiminase (ADI)-promotor. CDS, coding sequence; NR, nitroreductase; TrxR, thioredoxin reductase; TrxR DN, dominant negative mutation of TrxR. Please note that we have kept our original designations NR1, NR2, and NR3 to be in frame with our previous publications.
Table 1. List of the transgenic Giardia lamblia strains used in this study. All transgenes were fused to the strong, constitutive arginine deiminase (ADI)-promotor. CDS, coding sequence; NR, nitroreductase; TrxR, thioredoxin reductase; TrxR DN, dominant negative mutation of TrxR. Please note that we have kept our original designations NR1, NR2, and NR3 to be in frame with our previous publications.
StrainCharacteristicsAccessionReferences
WBC6Wildtype strain (WT)No transgene
C4Nitro drug-resistant strain derived from WTNo transgene
GusAWT overexpressing E. coli glucuronidase AE. coli K-12 CDS 1785 074 –1785 074
NR1WT overexpressing G. lamblia Fd-NR2 (“NR1”)
Higher susceptibility to nitro drugs in G. lamblia and E. coli		Giardia DB 22677[19,23]
NR2WT overexpressing G. lamblia Fd-NR1 (“NR2”)
Better resistance to nitro drugs in E. coli		Giardia DB 6175[21,23]
NR3WT overexpressing G. lamblia NR family protein (“NR3”)Giardia DB 15307[23]
TrxRWT overexpressing G. lamblia TrxR wildtypeGiardia DB 9827[24]
TrxR DNWT overexpressing G. lamblia TrxR-dominant negative allele
Higher susceptibility to oxygen		Giardia DB 9827[24]
This work
Table
Table 2. Metabolites identified in G. lamblia WBC6 wildtype trophozoites (WT; two independent experiments), in trophozoites overexpressing thioredoxin reductase (TrxR) or a dominant negative mutation of TrxR (TrxR DN) as a control, and in trophozoites overexpressing the nitroreductases NR1, NR2, NR3, or GusA as a control. The values represent normalized integrated peak areas obtained by integration of 1H-NMR resonances. Values are given as mean values ± SD of arbitrary units for n replicates as indicated. Metabolites with significant differences from the respective controls (i.e., TrxR DN for TrxR and GusA for the NRs) after Bonferroni correction for multiple comparisons are printed in bold (ANOVA followed by t-tests; +, p < 0.05; °, p < 0.01; *, p < 0.001). TMA: trimethylamine.
Table 2. Metabolites identified in G. lamblia WBC6 wildtype trophozoites (WT; two independent experiments), in trophozoites overexpressing thioredoxin reductase (TrxR) or a dominant negative mutation of TrxR (TrxR DN) as a control, and in trophozoites overexpressing the nitroreductases NR1, NR2, NR3, or GusA as a control. The values represent normalized integrated peak areas obtained by integration of 1H-NMR resonances. Values are given as mean values ± SD of arbitrary units for n replicates as indicated. Metabolites with significant differences from the respective controls (i.e., TrxR DN for TrxR and GusA for the NRs) after Bonferroni correction for multiple comparisons are printed in bold (ANOVA followed by t-tests; +, p < 0.05; °, p < 0.01; *, p < 0.001). TMA: trimethylamine.
MetaboliteWT 1
n = 10		WT 2
n = 4		TrxR
n = 10		TrxR_DN
n = 10		GusA
n = 6		NR1
n = 5		NR2
n = 6		NR3
n = 5
Carbohydrates and related
Acetate0.77 ± 0.050.85 ± 0.030.57 ± 0.030.65 ± 0.04 *0.65 ± 0.020.49 ± 0.02 *0.58 ± 0.02 +0.58 ± 0.02 +
Citrate0.17 ± 0.070.20 ± 0.010.18 ± 0.090.24 ± 0.060.11 ± 0.050.32 ± 0.100.18 ± 0.080.16 ± 0.05
Glucose-1-phosphate4.46 ± 0.651.65 ± 0.1 *4.59 ± 0.405.16 ± 0.661.69 ± 0.192.31 ± 0.13 °1.71 ± 0.241.88 ± 0.14
Amino acids
Alanine6.02 ± 0.777.63 ± 0.33 *4.88 ± 0.885.20 ± 0.796.95 ± 0.305.40 ± 0.256.00 ± 0.486.80 ± 0.36
Asparagine0.18 ± 0.020.26 ± 0.03 *0.19 ± 0.040.22 ± 0.030.22 ± 0.020.18 ± 0.070.21 ± 0.030.26 ± 0.02
Citrulline0.19 ± 0.090.26 ± 0.090.14 ± 0.110.18 ± 0.050.26 ± 0.071.16 ± 0.14 *0.46 ± 0.180.24 ± 0.13
Cysteine0.56 ± 0.120.29 ± 0.07 +0.40 ± 0.130.32 ± 0.040.22 ± 0.090.38 ± 0.070.29 ± 0.150.29 ± 0.14
Cystine (Cys-Cys)0.12 ± 0.040.13 ± 0.030.12 ± 0.050.15 ± 0.040.24 ± 0.050.47 ± 0.04 +0.42 ± 0.06 +0.27 ± 0.04
Glutamate7.03 ± 0.284.94 ± 0.31 *7.27 ± 0.257.12 ± 0.245.32 ± 0.155.35 ± 0.155.14 ± 0.125.41 ± 0.20
Glutamine0.48 ± 0.120.72 ± 0.240.45 ± 0.140.48 ± 0.090.94 ± 0.100.34 ± 0.09 *0.66 ± 0.320.91 ± 0.11
Glycine2.18 ± 0.102.85 ± 0.032.38 ± 0.092.21 ± 0.122.86 ± 0.062.40 ± 0.112.65 ± 0.032.96 ± 0.07
Histidine0.19 ± 0.040.28 ± 0.020.14 ± 0.040.16 ± 0.040.26 ± 0.020.14 ± 0.050.21 ± 0.030.22 ± 0.04
Isoleucine0.42 ± 0.180.68 ± 0.150.35 ± 0.100.35 ± 0.100.72 ± 0.120.68 ± 0.040.88 ± 0.110.77 ± 0.08
Leucine4.27 ± 0.185.89 ± 0.10 *3.92 ± 0.163.76 ± 0.175.95 ± 0.064.09 ± 0.185.52 ± 0.155.19 ± 0.14
Lysine4.67 ± 0.364.26 ± 0.185.45 ± 0.245.02 ± 0.24 °4.37 ± 0.274.68 ± 0.184.08 ± 0.244.33 ± 0.25
Methionine0.95 ± 0.081.06 ± 0.071.02 ± 0.050.98 ± 0.051.02 ± 0.051.30 ± 0.04 °1.15 ± 0.051.09 ± 0.05
Ornithine1.84 ± 0.131.71 ± 0.031.90 ± 0.111.72 ± 0.131.66 ± 0.022.65 ± 0.17 *1.92 ± 0.12 °1.81 ± 0.08
Phenylalanine0.30 ± 0.020.33 ± 0.020.28 ± 0.020.27 ± 0.040.29 ± 0.010.20 ± 0.04 +0.25 ± 0.040.28 ± 0.02
Pipecolic acid0.64 ± 0.080.69 ± 0.120.73 ± 0.080.72 ± 0.111.00 ± 0.041.23 ± 0.06 *1.07 ± 0.071.00 ± 0.10
Proline0.56 ± 0.120.40 ± 0.220.50 ± 0.120.50 ± 0.120.68 ± 0.080.54 ± 0.050.64 ± 0.170.66 ± 0.09
Threonine2.96 ± 0.122.61 ± 0.07 +3.27 ± 0.123.03 ± 0.142.75 ± 0.112.50 ± 0.072.83 ± 0.143.00 ± 0.07
Tryptophan0.11 ± 0.030.07 ± 0.030.09 ± 0.020.10 ± 0.030.05 ± 0.020.04 ± 0.020.08 ± 0.040.03 ± 0.01
Tyrosine0.38 ± 0.040.37 ± 0.020.37 ± 0.020.34 ± 0.020.36 ± 0.020.26 ± 0.02 *0.32 ± 0.070.36 ± 0.02
Valine4.24 ± 0.195.85 ± 0.08 *4.72 ± 0.204.46 ± 0.265.87 ± 0.135.27 ± 0.195.94 ± 0.255.93 ± 0.16
Cofactors
Adenosine triphosphate0.25 ± 0.030.22 ± 0.020.27 ± 0.040.25 ± 0.030.24 ± 0.020.26 ± 0.020.30 ± 0.050.25 ± 0.02
Nicotinamide-adenine-dinucleotide 0.08 ± 0.030.06 ± 0.010.08 ± 0.010.08 ± 0.020.07 ± 0.030.07 ± 0.020.08 ± 0.030.08 ± 0.01
Other
Trimethylamine1.18 ± 0.041.09 ± 0.02 +1.16 ± 0.061.16 ± 0.061.12 ± 0.030.71 ± 0.03 *0.87 ± 0.08 °1.04 ± 0.03

Share and Cite

MDPI and ACS Style

Müller, J.; Vermathen, M.; Leitsch, D.; Vermathen, P.; Müller, N. Metabolomic Profiling of Wildtype and Transgenic Giardia lamblia Strains by 1H HR-MAS NMR Spectroscopy. Metabolites 2020, 10, 53. https://doi.org/10.3390/metabo10020053

AMA Style

Müller J, Vermathen M, Leitsch D, Vermathen P, Müller N. Metabolomic Profiling of Wildtype and Transgenic Giardia lamblia Strains by 1H HR-MAS NMR Spectroscopy. Metabolites. 2020; 10(2):53. https://doi.org/10.3390/metabo10020053

Chicago/Turabian Style

Müller, Joachim, Martina Vermathen, David Leitsch, Peter Vermathen, and Norbert Müller. 2020. "Metabolomic Profiling of Wildtype and Transgenic Giardia lamblia Strains by 1H HR-MAS NMR Spectroscopy" Metabolites 10, no. 2: 53. https://doi.org/10.3390/metabo10020053

APA Style

Müller, J., Vermathen, M., Leitsch, D., Vermathen, P., & Müller, N. (2020). Metabolomic Profiling of Wildtype and Transgenic Giardia lamblia Strains by 1H HR-MAS NMR Spectroscopy. Metabolites, 10(2), 53. https://doi.org/10.3390/metabo10020053

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop