Next Article in Journal
Serum Metabolomic Analysis of Male Patients with Cannabis or Amphetamine Use Disorder
Previous Article in Journal
The Endogenous Metabolite Glycerophosphocholine Promotes Longevity and Fitness in Caenorhabditis elegans
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Metabolites
Volume 12
Issue 2
10.3390/metabo12020178
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

Combined Plasma and Urinary Metabolomics Uncover Metabolic Perturbations Associated with Severe Respiratory Syncytial Viral Infection and Future Development of Asthma in Infant Patients

by
Shao Thing Teoh
1,
Mara L. Leimanis-Laurens
2,3,
Sarah S. Comstock
4,
John W. Winters
2,3,
Nikita L. Vandenbosch
3,
Jeremy W. Prokop
2,
André S. Bachmann
2,
Sophia Y. Lunt
1,5,* and
Surender Rajasekaran
2,3,6,*
1
Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, MI 48824, USA
2
Department of Pediatrics and Human Development, College of Human Medicine, Michigan State University, Grand Rapids, MI 49503, USA
3
Pediatric Intensive Care Unit, Helen DeVos Children’s Hospital, Grand Rapids, MI 49503, USA
4
Department of Food Science and Human Nutrition, Michigan State University, East Lansing, MI 48824, USA
5
Department of Chemical Engineering and Materials Science, Michigan State University, East Lansing, MI 48824, USA
6
Office of Research, Spectrum Health, Grand Rapids, MI 49503, USA
*
Authors to whom correspondence should be addressed.
Metabolites 2022, 12(2), 178; https://doi.org/10.3390/metabo12020178
Submission received: 15 December 2021 / Revised: 21 January 2022 / Accepted: 11 February 2022 / Published: 14 February 2022
(This article belongs to the Section Endocrinology and Clinical Metabolic Research)
Browse Figures
Review Reports Versions Notes

Abstract

:
A large percentage of infants develop viral bronchiolitis needing medical intervention and often develop further airway disease such as asthma. To characterize metabolic perturbations in acute respiratory syncytial viral (RSV) bronchiolitis, we compared metabolomic profiles of moderate and severe RSV patients versus sedation controls. RSV patients were classified as moderate or severe based on the need for invasive mechanical ventilation. Whole blood and urine samples were collected at two time points (baseline and 72 h). Plasma and urinary metabolites were extracted in cold methanol and analyzed by liquid chromatography coupled with tandem mass spectrometry (LC-MS/MS), and data from the two biofluids were combined for multivariate data analysis. Metabolite profiles were clustered according to severity, characterized by unique metabolic changes in both plasma and urine. Plasma metabolites that correlated with severity included intermediates in the sialic acid biosynthesis, while urinary metabolites included citrate as well as multiple nucleotides. Furthermore, metabolomic profiles were predictive of future development of asthma, with urinary metabolites exhibiting higher predictive power than plasma. These metabolites may offer unique insights into the pathology of RSV bronchiolitis and may be useful in identifying patients at risk for developing asthma.

Graphical Abstract

1. Introduction

Respiratory syncytial virus (RSV) is a major cause of childhood morbidity across the world and mortality in low to middle income countries [1]. In the United States, RSV represents the most frequent cause of hospitalization in the first year of life [2]. Twenty percent of the annual birth cohort develops viral bronchiolitis that requires medical attention [3], and acutely ill patients have been reported to continue having other respiratory problems requiring medical attention throughout their childhood [4]. Improved understanding of the molecular mechanisms underlying development of severe symptoms in RSV infection is therefore crucial to identify potential risk factors for both severe RSV and other associated respiratory problems, as well as improve treatment strategies. We previously reported the results of a multidimensional RNA analysis, including gene, transcript, immune cell deconvolution, secondary infections, and immune repertoire comparison between RSV patients and controls, and discovered that infants with severe RSV had significant blood signatures connected to the immune system, interferon signaling, and cytokine activation [5]. Metabolomics has been used in a number of studies to identify associations between metabolic changes in biofluids and perturbations to the immune system; for example, altered levels of aromatic amino acids, fatty acids and energy metabolites have been reported to associate with serum inflammatory markers in a number of different medical conditions [6,7,8,9]. Here, we extend our study of the aforementioned cohort by using liquid chromatography–tandem mass spectrometry (LC-MS/MS) metabolomics to characterize metabolic perturbations associated with RSV severity in plasma and urine of these patients. LC-MS/MS has been employed in a number of metabolomic analyses of clinical samples including urinary and serum or plasma [10,11,12,13]; although less comprehensive in coverage than untargeted metabolomics, targeted LC-MS/MS has the advantage of increased sensitivity and quantitative performance [14,15], and our established method [16,17,18,19] is capable of profiling a wide range of metabolites including amino acids, phosphorylated compounds in glycolysis and pentose phosphate pathway, organic acids in the tricarboxylic acid (TCA) cycle, as well as intermediates in nucleotide biosynthesis and nucleotide sugar metabolism to detect perturbations within a variety of metabolic pathways. Since metabolomics are predictive of the development of disease [20], biomarkers identified in this study may also have prognostic utility in determining whether severe disease may develop after initial hospitalization.
An epidemiological link between bronchiolitis caused by RSV and the development of recurrent wheezing and/or asthma was first described more than 60 years ago [21]. Together with genetic predisposition to allergic sensitization (atopy), RSV infection in early life is now recognized as one of the most important risk factors for pediatric asthma [22,23,24]. Here, we also tested the hypothesis that RSV patient metabolomic profiles may be predictive for future development of asthma over the two years after pediatric ICU admission.

2. Results

2.1. Patient Demographics

A detailed description of the patient cohort can be found in our previous blood transcriptomics study [5] and are summarized in Table 1. RSV patients requiring invasive mechanical ventilation were considered severe, while those only requiring non-invasive respiratory support were considered moderate. From this patient cohort, we identified plasma and urinary metabolites correlated with severity of RSV infection, and additionally found urinary metabolites with strong predictive power for future asthma development in RSV patients.
Additional details of the patient demographics are included in Supplementary Table S1. Weight at time of study enrollment PICU admission was found to be statistically lower in RSV patients as compared to control (4.25 kg vs. 7.32 kg; p ≤ 0.005). It was a largely Caucasian population, and there were three premature (less than 36 weeks) infants among the RSV patients. All RSV patients presented with acute bronchiolitis, and a significant number (15/20) were also diagnosed with pneumonia; detailed diagnoses as well as the severity of illness scores for the 20 RSV patients are summarized in Supplementary Table S2. PELOD was included, given that it is a measure of the first 10 days in the ICU [25]. Both comprehensive and Day 0 PELOD scores were similar for both severe and moderate patient populations. The greatest variability was seen on day 3, whereby moderate patients had less than half the scores of the severe (2.75 vs. 7.93). Although included in the table, both Pediatric Index of Mortality 2 (PIM2) [26] and Pediatric Risk of Mortality III (PRISMIII) [27] are calculated during the first hours of ICU admission only, the time to consenting patients was within the first 24 h of admission, and their relevance is not clear in this context.

2.2. Principal Components Analysis Reveals Sample Clustering by RSV Severity

We quantified a total of 64 plasma and 70 urinary metabolite peaks over two time points. Unsupervised principal component analysis (PCA) was performed on the combined metabolite profiles to address sample clustering. The first two principal components (PCs) were each driven by one or two strong outliers—data points that are highly divergent from the majority of samples (Appendix A Figure A1). We therefore focused on PCs 3 and 4 and found that these two PCs cluster the samples into distinct control (non-RSV), severe RSV baseline, and severe RSV 72 h groups (Figure 1A), indicating that metabolic differences corresponding to RSV severity are reflected in the data. Control and severe RSV samples separated along PC 3, while the PC 4 axis separated baseline and 72 h samples. Most moderate RSV baseline samples fell in between the severe RSV baseline and control samples, which is consistent with these samples originating from patients exhibiting a moderate pathology.

2.3. Plasma Sialic Acid Metabolic Intermediates and Urinary Nucleotides Discriminate RSV Samples

To further investigate the differences between RSV and control patients, we performed PLS regression on combined plasma and urinary metabolites, using RSV severity as the response variable. Numerical values for RSV severity were assigned as follows: control, 0; moderate RSV, 1; severe RSV, 2. Since 72 h samples were not available for all RSV patients, nor control patients, we limited the analysis to baseline samples for the RSV patients. We first investigated the PLS score plot and confirmed that the model successfully identified metabolite trends that correspond to RSV severity, with control and severe RSV samples falling on opposite ends of the plot and moderate RSV samples in the middle (Figure 1B).
To identify plasma metabolites that correlate positively or negatively with RSV severity, we examined the PLS coefficients (Figure 2A, Supplementary Table S3). Metabolites with large positive coefficients were positively correlated with the RSV response variable. Plasma metabolites with the highest positive coefficients were xanthosine, phosphoserine, and N-acetylneuraminate. Box plots of the normalized metabolite levels showed that these metabolites were highest in severe RSV and lowest in control groups (xanthosine, p = 7.34 × 10−6; phosphoserine, p = 8.5 × 10−5; N-acetylneuraminate, p = 1.84 × 10−5) (Figure 2B). Meanwhile, metabolites with negative coefficients included N-acetylglucosamine-1-phosphate, glycerol 3-phosphate, alpha ketoglutarate, and CMP-N-acetylneuraminate. In the box plots, N-acetylglucosamine-1-phosphate showed a clear stepwise decrease from controls to moderate to severe RSV (control vs. moderate RSV, p = 0.04; control vs. severe RSV, p = 4.28 × 10−6), while glycerol 3-phosphate, alpha ketoglutarate and CMP-N-acetylneuraminate also significantly decreased between control and severe RSV samples (glycerol 3-phosphate, p = 0.0163; alpha ketoglutarate, p = 0.00603; CMP-N-acetylneuraminate, p = 0.0357) (Figure 2C).
Next, we examined the PLS coefficients for urinary metabolites (Figure 3A, Supplementary Table S4). Adenosine monophosphate (AMP), uridine monophosphate (UMP), xylulose-5-phosphate, and inosine monophosphate (IMP) had large positive coefficients indicating correlation with RSV severity. From boxplots of the normalized metabolite levels, we found that the positive correlation of these metabolites with RSV severity was driven by large increases in severe RSV samples (AMP, p = 0.02; UMP, p = 0.00872; xylulose-5-phosphate, p = 0.207; IMP, p = 0.0408), while the metabolite levels were not significantly different between control and moderate RSV samples (Figure 3B). Conversely, metabolites with large negative coefficients included citrate/isocitrate, guanosine diphosphate (GDP) and cytidine diphosphate (CDP). Citrate/isocitrate showed a stepwise decrease from control to moderate and further to severe RSV (control vs. moderate RSV, p = 0.0303; control vs. severe RSV, p = 0.000361), while GDP (control vs. moderate RSV, p = 0.0303; control vs. severe RSV, p = 0.00227) and CDP (control vs. moderate RSV, p = 0.00433; control vs. severe RSV, p = 0.00318) were decreased in RSV samples compared to controls but were not different between moderate and severe RSV (Figure 3C).

2.4. Patient Weight and Plasma Protein Levels Negatively Correlate with RSV Severity

Unsupervised and supervised multivariate analyses identified plasma and urinary metabolites altered in RSV relative to control. Using these multivariate analysis strategies also enabled investigation of the relationships between other patient data. We performed PCA on available clinical covariates including white blood cell measurements. Protein levels quantified in the biofluid samples following metabolite extraction were included as variables in this analysis, and categorical data such as gender, race, prematurity, diet, presence of ENT problems, or classification as asthmatic within a two-year follow-up period were converted to numerical variables (Appendix A Table A1).
Length on ventilator, RSV severity, and length of stay had positive loadings on PC 1 in this model, indicating that they were all strongly correlated with each other (Figure 4A, Supplementary Table S5). We noted that the future asthma classification covariate also had a positive loading, indicating that future classification as asthmatic correlated with severe RSV. Weight, plasma protein, and age had negative loadings, indicating that they were negatively correlated with RSV severity. Thus, severe RSV patients tended to have lower weight, lower levels of plasma protein, as well as younger age. These trends were confirmed from the PCA score plot, with the color coding of data points varied to visualize RSV severity (Figure 4B), weight (Figure 4C), plasma protein (Figure 4D), and asthma classification (Figure 4E). The PC 2 of this model seems to be driven by an inverse relationship between prematurity and ear nose throat (ENT) problems, as well as birth weight (Appendix A Figure A2). The relatively high positive loading of eosinophil counts on this PC also suggests a correlation between prematurity and eosinophils.

2.5. Urinary Metabolites Are Predictive of Future Asthma Development

As noted above, asthma classification showed a correlation with RSV severity. Only one out of 10 control patients were classified as asthmatic in the two-year follow-up period after hospital visit. A larger proportion of moderate and severe RSV patients (3 out of 5 moderate RSV, 7 out of 15 severe RSV; in total 10 out of 20 RSV patients, or 50%) were classified as asthmatic in the following two years (Figure 5A). To investigate whether metabolite profiles of samples taken at initial hospitalization is predictive of future diagnosis of asthma, we performed PLS regression on combined plasma and urinary metabolites, using asthma classification as the response variable. Since RSV infection exerts a strong effect on metabolite profiles of baseline samples (Figure 1B), we only used metabolite data from control and 72 h timepoint RSV samples for asthma prediction to minimize potential confounding effects from RSV-driven metabolic perturbations. This strategy resulted in two out of three moderate RSV patients who developed asthma being excluded from the analysis, since 72 h samples were not collected for these patients due to early discharge from intensive care. The asthma PLS model separated asthma and non-asthma samples as indicated by the score plot (Figure 5B). Since the metabolite profiles were sampled during initial hospitalization and the patient classification as asthmatic occurred at a later date, this suggests that metabolite profiles are predictive of future asthma development in infant patients. The score plot showed that non-asthma severe RSV samples, while still clustering together with non-asthma control and moderate RSV samples, tended to fall closer to the asthma group (Figure 5B) indicating a potential underlying predisposition of severe RSV patients to develop asthma. To exclude the possibility of overlap between high-contributing metabolites of the RSV and the asthma models, we plotted PLS coefficients from the RSV model versus the asthma model and determined that there was low correlation (Pearson’s r = 0.417) between the two sets of coefficients (Appendix A Figure A3).
Examination of the PLS coefficients in the asthma prediction model showed that many urinary metabolites had positive coefficients and therefore are positively correlated with asthma classification, including lysine, ornithine, UMP, IMP and GDP-mannose (Figure 5C, Supplementary Table S6). The metabolite with largest negative coefficient (i.e., negatively correlated with asthma classification) was also a urinary metabolite, xanthosine. Receiver–operator characteristic (ROC) curves using either top positive and negative urinary metabolites (lysine and xanthosine) or top positive and negative plasma metabolites (uracil and citrate/isocitrate) showed area under the curve (AUC) values close to 1, indicating good predictive power of both top urinary and plasma metabolites (Appendix A Figure A4). However, AUC was larger with urinary metabolites (0.956) than plasma metabolites (0.895). Thus, urinary metabolites have a higher predictive power toward future asthma development compared to plasma metabolites.

3. Discussion

In this paper, we report results of a metabolomic analysis of blood and urine samples from RSV patients and healthy controls, which complements a previous transcriptomics study performed on the same cohort [5]. The current study is the combined analysis of plasma and urine, which allows us to explore associations between metabolites extracted from the two biofluids in our multivariate analysis models.
We found significant metabolite differences in the plasma and urinary metabolites of RSV patients compared to non-RSV controls. Among urinary metabolites, citrate/isocitrate strongly and negatively correlated with RSV severity (Figure 3A). This result is consistent with previous studies comparing urinary metabolites between RSV patients and healthy controls [28,29]. The studies also found perturbations in the TCA cycle that included decreased succinate in urine samples of RSV patients; in agreement with this, we found lower succinate levels in RSV samples relative to controls (Appendix A; Figure A5).
Other significant perturbations to urinary metabolites appear to involve nucleotides, specifically increased monophosphate nucleotides (AMP, UMP, IMP) and decreased diphosphate nucleotides (GDP, CDP) in RSV samples. Monophosphate nucleotides are converted to diphosphate nucleotides by nucleoside–monophosphate kinases (NMPK) [30], while the opposite conversion from diphosphate to monophosphate nucleotides is mediated by nucleoside pyrophosphatase/phosphodiesterases (NPPs) [31]. Hence, accumulation of monophosphate nucleotides and decrease in diphosphate nucleotides may be due to either decreased NMPK activity, and/or increased NPP activity in urine-associated tissues such as the kidney or urinary tract.
We also observed apparent perturbations in the plasma metabolites of RSV patients compared to controls, particularly to the hexosamine and sialic acid biosynthesis pathways. N-acetylglucosamine-1-phosphate is an intermediate in hexosamine biosynthesis, which produces UDP-N-acetylglucosamine that in turn feeds into the sialic acid biosynthesis pathway. The product of sialic acid biosynthesis, N-acetylneuraminate, is further conjugated with CMP to form CMP-N-acetylneuraminate, the activated form that is used as substrate by sialyltransferases; sialyltransferase reactions attach the N-acetylneuraminate moiety onto glycan chains of glycoproteins and glycolipids destined for the cell surface. In addition to N-acetylglucosamine-1-phosphate, which had the largest negative coefficient in the RSV model (Figure 2A), both UDP-N-acetylglucosamine and CMP-N-acetylneuraminate were also decreased in RSV samples (Appendix A Figure A6). The decreased levels of these metabolites suggest a decrease in hexosamine and sialic acid pathway activities. Free N-acetylneuraminate positively correlated with RSV severity (Figure 2B). Since N-acetylneuraminate can be recycled from endocytosed cell surface proteins or taken up directly from the extracellular environment, it is possible that dysregulation in the pathway resulted in N-acetylneuraminate accumulation despite decreased biosynthetic activity. We previously reported an increase in N-acetylneuraminate and xanthosine levels in highly metastatic mouse mammary tumors relative to less metastatic tumors [17]. In the present study, we also identified a positive correlation between xanthosine and RSV severity, which echoes the results from the previous study. Free xanthosine can arise from degradation of xanthosine monophosphate (XMP), an intermediate in de novo guanosine nucleotide biosynthesis [32]. XMP is produced from IMP by inosine 5′-monophosphate dehydrogenase (IMPDH); this reaction is important for DNA synthesis particularly in rapidly proliferating cells including microbes, cancer cells and cells of the immune system [33]. Thus, an increase in xanthosine levels may reflect increased IMPDH activity following immune system activation in the RSV patients [34].
Respiratory viral infections are the most important risk factors for the onset of wheezing in infants and small children [22,35]. It remains unclear whether bronchiolitis causes chronic respiratory symptoms, or if it is a marker for children with a genetic/environmental predisposition for developing asthma. Here, we found that metabolite profiles from this patient cohort could be used to predict whether patients hospitalized for acute RSV bronchiolitis will become classified as asthmatic within a two-year follow-up period (Figure 5B). Although most patients who developed asthma within this cohort are RSV patients (Figure 5A), the metabolite coefficients in the RSV characterization and asthma prediction models do not show a high degree of correlation (Appendix A Figure A3), indicating that asthma predisposition is driven by a separate set of underlying metabolic factors independent of RSV severity. We additionally found that urinary metabolites have high predictive power for future asthma risk (Figure 5C). The initially surprising result that urinary metabolites have stronger predictive value than plasma could be due to the extensive contact of urine with tubule epithelial cells of the kidney, which share many of the cell properties of lung epithelial cells. As urine is one of the simplest to obtain materials from patients, these correlations could have high future clinical utility, especially if these metabolites are studied in a broader non-RSV infant cohort to establish if the RSV and metabolite correlations are dependent on each other.
Symptoms of both RSV and asthma are characterized by interferon responses within the respiratory system, with a complex interaction of immune cells with lung epithelial cells [36]. As nearly every infant is exposed to RSV, but only a small percentage are hospitalized, the priming of the innate immune system by other factors is likely responsible for the severe phenotype in many cases. This same priming is likely shared in asthma. Culture experiments have shown that Th1 or Th2 cytokine priming can overstimulate the immune system and independently increase risks for both severe RSV responses and asthma development [37]. Therefore, it has been suggested that asthma development is likely independent of RSV infections. The notion of a shared underlying cause of immune priming but distinct pathways of diseases development for severe RSV and asthma is reflected in our study results, which showed that patients hospitalized for severe RSV have a higher risk of asthma development (Figure 5A) but metabolic perturbations due to RSV and metabolic signatures predictive of asthma are distinct (Appendix A Figure A4).
One caveat of the current study is that we performed relative quantification of metabolites. We identified metabolites perturbed in RSV infection as well as metabolites predictive of asthma development; our data in the present state do not allow for establishment of concentration thresholds to be used as biomarkers for RSV diagnosis or asthma prediction. Absolute quantification of metabolite concentrations in biofluids require labeled internal standards for compounds of interest, which were not available for this study. As a future direction, we propose validating the results of this study in a separate cohort and then performing absolute quantification of the metabolites of interest to determine threshold concentrations for biomarkers.
In conclusion, we identified a perturbation in plasma metabolites in the sialic acid pathway that correlate with severity of RSV infection. Changes in urinary metabolites, especially nucleotides, also correlate with RSV severity. Additionally, we determined that urinary metabolites predict future development of asthma. These results will inform further mechanistic studies into RSV pathology as well as development of biomarkers for predicting future asthma risk to inform early detection and intervention.

4. Materials and Methods

4.1. Study Population, Site and Sample Collection

Samples were collected under the IRB protocol 2107-049-SH/HDVCH. Samples were collected at the Helen DeVos Children’s Hospital (HDVCH), a quaternary-care, urban, pediatric hospital in Western, Michigan. In brief, study patients were admitted to the PICU with bronchiolitis, were positive for RSV by nasopharyngeal test, and were less than six months of age (excluding premature births less than 34 weeks gestation). Consent was obtained, and then blood and urine samples were collected. At least 0.5 mL of blood was collected in EDTA tubes on ice and immediately centrifuged at 4 °C to obtain plasma. Plasma samples were stored at −20 °C overnight and then at −80 °C until batch processing. Urine samples were collected under clinically sterile conditions at the bedside using a substrate-free vacutainer (BD® #364992-1) according to manufacturer specifications and frozen at −20 °C for short-term storage, and −80 °C for long-term storage.
RSV patient samples were collected at two time points; within 24 h of PICU admission, and directly prior to intravenous (IV) catheter removal and discharge from the PICU/hospital (on average three-four days later), which reflected the sampling of the acute and stabilization/recovery phase of the illness. In total, 20 patients with RSV bronchiolitis consented. Patients were divided into either severe or moderate groupings based on need for invasive mechanical ventilation.
A total of 10 healthy, age-matched, sedation-control patients with normal airways and lungs, presenting for routine IV sedation had blood and urine samples collected at one time point. Patient cohort details are provided in Table 1. Patient data were collected, de-identified and managed using Research Electronic Data Capture (REDCap) tools.

4.2. Metabolite Extraction

Plasma and urine samples were thawed on ice, and 20 µL plasma or 50 µL urine were added to 500 µL ice-cold HPLC-grade methanol containing 1 µM piperazine-N,N′-bis (PIPES) and 1 µM camphorsulfonic acid as internal standards. Samples were vortexed and centrifuged at 16,100× g for 5 min at 4 °C to precipitate protein and other cell debris. The metabolite extracts were dried under nitrogen gas and stored at −80 °C until analysis. Protein precipitated during the extraction was dissolved in 0.2 mM KOH overnight and then quantified by BCA protein assay using Pierce BCA Protein Assay Kit (Thermo Fisher, Waltham, MA, USA).

4.3. Metabolite Analysis by LC-MS/MS

Dried metabolite extracts were suspended in 100 µL HPLC-grade water and centrifuged at 16,100× g for 5 min at 4 °C to remove any precipitate. The aqueous resuspensions were diluted 1:1 with 20 mM perfluoroheptanoic acid (PFHA) to obtain metabolite suspensions in 10 mM PFHA for amino acid analysis. Samples were analyzed by ion-pairing reverse phase chromatography using an Acquity HSS T3 column (2.1 × 100 mm, 1.8 µm, Waters) for separation and a Waters Quattro Micro triple quadrupole mass spectrometer operated in positive mode as mass analyzer. The LC parameters were as follows: autosampler temperature, 10 °C; injection volume, 10 µL; column temperature, 40 °C. The LC solvents were A: 1 mM PFHA in water and B: acetonitrile. Elution from the column was performed over 13 min with flow rate fixed at 0.3 mL/min and using the following gradient, defined in terms of acetonitrile %: t = 0, 0%; t = 1, 0%; t = 8, 65%; t = 8.01, 90%; t = 9, 90%; t = 9.01, 0%; t = 13, 0%. Mass spectra were acquired using positive-mode electrospray ionization operating in multiple reaction monitoring (MRM) mode. The capillary voltage was 1000 V, and cone voltage was 45 V. Nitrogen was used as cone gas and desolvation gas, with flow rates of 40 and 800 L/h, respectively. The source temperature was 120 °C, and desolvation temperature was 350 °C. Argon was used as collision gas at a manifold pressure of 2.31 × 10−3 mbar.
A second half of the aqueous resuspensions were diluted 1:1 with 2× tributylamine (TBA) solvent (see below) to obtain metabolite suspensions in 1× TBA solvent for analysis of other metabolites, including glycolytic intermediates, TCA cycle intermediates, nucleotides, and nucleotide sugars. These samples were analyzed by ion-pairing reverse phase chromatography using an Ascentis Express C18 column (5 cm × 2.1 mm, 2.7 µm, Sigma-Aldrich, St. Louis, MO, USA) for separation and a Waters Xevo TQ-S triple quadrupole mass spectrometer as mass analyzer. The LC parameters were as follows: autosampler temperature, 4 °C; injection volume, 5 µL; column temperature, 50 °C. The LC solvents were 10 mM TBA and 15 mM acetic acid in 97:3 water:methanol (“TBA solvent”) and methanol. Elution from the column was performed over 12 min with the following varying flow rate gradient, defined in terms of methanol %: t = 0, 0%, flow rate 0.4 mL/min; t = 1, 0%, flow rate 0.4 mL/min; t = 2, 20%, flow rate 0.3 mL/min; t = 3, 20%, flow rate 0.25 mL/min; t = 5, 55%, flow rate 0.15 mL/min; t = 8, 95%, flow rate 0.15 mL/min; t = 9.5, 95%, flow rate 0.15 mL/min; t = 10, 0%, flow rate 0.4 mL/min; t = 12, 0%, flow rate 0.4 mL/min. Mass spectra were acquired using negative-mode electrospray ionization operating in MRM mode. The capillary voltage was 3000 V, and cone voltage was 50 V. Nitrogen was used as cone gas and desolvation gas, with flow rates of 150 and 600 L/h, respectively. The source temperature was 150 °C, and desolvation temperature was 500 °C. Argon was used as collision gas at a manifold pressure of 4.3 × 10−3 mbar.
On both instruments, sample order was randomized to avoid systematic bias, and each sample was run twice as analytical replicates. MRM transitions for the two instruments are listed in Supplementary Table S7.

4.4. Data Processing and Analysis

Peak processing was performed using MAVEN [38]. For each plasma or urine sample, data from the Quattro Micro (amino acids) and TQ-S (other metabolites) runs were separately scaled to internal standard peak (camphorsulfonic acid) measured in the respective run, and then data from both platforms were combined. Samples were further normalized within the set (plasma or urine) using Probabilistic Quotient Normalization [39] to account for concentration differences in the biofluids. Analytical replicates were then averaged to obtain the final metabolite profile for each sample. Principal component analysis (PCA) and partial least squares (PLS) regression were performed in R using the ropls package [40] using default settings (unit variance scaling, no log-transformation, leave-one-out cross validation). Metabolite boxplots were plotted in R and statistical comparisons were performed using the Mann–Whitney U test.

4.5. Clinical Covariates

Urine and plasma protein levels were quantified by BCA assay following methanol precipitation and metabolite extraction as described above (Section 4.2).
Classification of patients as asthmatic was performed by pediatric pulmonologists at the hospital during outpatient visits in the two-year period following enrollment in this study. Asthma diagnosis was based on responsiveness to bronchodilators following episodic symptoms of airflow obstruction or airway hyperresponsiveness, as standardized pulmonary function testing such as spirometry could not be performed in patients under 5 years of age [41].
White blood cell count values, including total WBC, lymphocytes, eosinophils and neutrophils were obtained from Complete Blood Count (CBC) tests performed on a Sysmex XN 10 Automated Hematology Analyzer as part of the routine care for inpatients.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo12020178/s1, Supplementary Table S1: Patient demographics of RSV Study (N = 30), Supplementary Table S2: Diagnosis of RSV positive patients (N = 20), Supplementary Table S3: PLS coefficients of plasma metabolites in RSV model, Supplementary Table S4: PLS coefficients of urinary metabolites in RSV model, Supplementary Table S5: PCA loadings in the covariate model, Supplementary Table S6: PLS coefficients of metabolites in asthma prediction model, Supplementary Table S7: MRM transitions for target metabolites.

Author Contributions

Conceptualization, S.R. and S.Y.L.; methodology, S.R., A.S.B., S.Y.L., M.L.L.-L. and S.T.T.; validation, S.R., S.Y.L., S.T.T. and M.L.L.-L.; formal analysis, S.T.T., S.S.C. and J.W.P.; investigation, S.T.T. and J.W.W.; resources, S.R. and S.Y.L.; data curation, M.L.L.-L., S.T.T. and S.R.; writing—original draft preparation, S.T.T.; writing—review and editing, all authors; visualization, S.T.T.; supervision, S.R., M.L.L.-L. and S.Y.L.; project administration, S.R., A.S.B., S.Y.L. and M.L.L.-L.; funding acquisition, S.R., S.Y.L., A.S.B. and M.L.L.-L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Spectrum Health Office of Research Administration, grant number RI0013-18 (S.Y.L., S.R., A.S.B.), and the APC was funded by Spectrum Health. S.Y.L. and S.T.T. were supported by a METAvivor Early Career Investigator Grant to S.Y.L.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Institutional Review Board of Helen DeVos Children’s Hospital (IRB protocol 2017-049-SH/HDVCH).

Informed Consent Statement

Informed consent was obtained from at least one parent or guardian of all subjects involved in the study.

Data Availability Statement

The data presented in this study are available in article and Supplementary Materials.

Acknowledgments

We thank the MSU Mass Spectrometry and Metabolomics Core for supporting our metabolomics studies.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

AMP: adenosine monophosphate; AUC: area under the curve; CDP: cytidine diphosphate; ENT: ear nose throat; GDP: guanosine diphosphate; IMP: inosine monophosphate; IMPDH: inosine 5′-monophosphate dehydrogenase; IV: intravenous; LC-MS/MS: liquid chromatography coupled with tandem mass spectrometry; LOS: length of stay; MRM: multiple reaction monitoring; NMPK: nucleoside-monophosphate kinases; NPPs: nucleoside pyrophosphatase/phosphodiesterases; PCA: principal component analysis; PFHA: perfluoroheptanoic acid; PICU: pediatric intensive care unit; PIPES: piperazine-N,N′-bis; PLS: partial least squares; ROC: receiver operator curve; RSV: respiratory syncytial viral; TBA: tributylamine; TCA: tricarboxylic acid; UMP: uridine monophosphate; WBC: white blood cells; XMP: xanthosine monophosphate.

Appendix A

Table
Table A1. Data values used in principal components analysis of clinical data.
Table A1. Data values used in principal components analysis of clinical data.
ValueValue TypeCategorical to Numeric Conversion
Asthma classificationcategoricalno (0), yes (1)
Dietcategoricalbreast-fed (0), mixed diet (1), formula-fed (2)
ENT 1 problemscategoricalno (0), yes (1)
Gendercategoricalmale (1), female (2)
Prematuritycategoricalno (0), yes (1)
RacecategoricalWhite (1), Other (2)
RSV severitycategoricalcontrol (0), moderate (1), severe (2)
Sample timepointcategoricalbaseline (0), 72 h (1)
Agenumeric (weeks)
Birth weightnumeric (kg)
Eosinophilsnumeric (count)
Length of staynumeric (h)
Length on ventilatornumeric (h)
Lymphocytesnumeric (count)
Neutrophilsnumeric (count)
Plasma proteinnumeric (µg/mL)
Urine proteinnumeric (µg/mL)
WBC 2numeric (count)
Weight at hospitalizationnumeric (kg)
1 ENT: ear nose throat; 2 WBC: white blood cells.
Metabolites 12 00178 g0a1 550
Figure A1. The first two principal components are driven by strong outliers and do not reflect RSV severity. Score plot of principal components analysis (PCA) performed using plasma and urinary metabolite data all samples. Horizontal and vertical axes represent the 1st and 2nd principal component scores, respectively. PC 1 is driven by a single outlier (PDM06), while PC 2 is driven by two outliers (PDM17 and PDM40).
Figure A1. The first two principal components are driven by strong outliers and do not reflect RSV severity. Score plot of principal components analysis (PCA) performed using plasma and urinary metabolite data all samples. Horizontal and vertical axes represent the 1st and 2nd principal component scores, respectively. PC 1 is driven by a single outlier (PDM06), while PC 2 is driven by two outliers (PDM17 and PDM40).
Metabolites 12 00178 g0a1
Metabolites 12 00178 g0a2 550
Figure A2. PCA of clinical covariates reveals inverse relationship between prematurity and ENT problems. (A) Covariate loadings on the 2nd principal component of the covariates PCA model. Covariates with positive loadings have higher values in samples with positive 2nd principal component scores, and vice versa; hence, covariates on the same side of the loading plot are positively correlated with each other, and covariates on the opposite side are inversely correlated. (B–E) Score plots all showing the 1st and 2nd principal components in the covariates-only PCA model, with sample data points color-coded according to (B) prematurity, (C) ENT problems, (D) eosinophils count (empty squares indicate that eosinophils were not quantified for these patients), and (E) birth weight.
Figure A2. PCA of clinical covariates reveals inverse relationship between prematurity and ENT problems. (A) Covariate loadings on the 2nd principal component of the covariates PCA model. Covariates with positive loadings have higher values in samples with positive 2nd principal component scores, and vice versa; hence, covariates on the same side of the loading plot are positively correlated with each other, and covariates on the opposite side are inversely correlated. (B–E) Score plots all showing the 1st and 2nd principal components in the covariates-only PCA model, with sample data points color-coded according to (B) prematurity, (C) ENT problems, (D) eosinophils count (empty squares indicate that eosinophils were not quantified for these patients), and (E) birth weight.
Metabolites 12 00178 g0a2
Metabolites 12 00178 g0a3 550
Figure A3. Asthma and severe RSV have distinct metabolic signatures. Scatterplot of RSV model PLS coefficients (vertical axis) against asthma prediction model PLS coefficients (horizontal axis). Positive and negative values on each axis represent direction and strength of correlation of individual metabolites with the response variable represented by that axis. Although severe RSV and asthma classification are loosely correlated (Figure 4A and Figure 5A), the individual metabolites most strongly associated with RSV severity and asthma classification are distinct.
Figure A3. Asthma and severe RSV have distinct metabolic signatures. Scatterplot of RSV model PLS coefficients (vertical axis) against asthma prediction model PLS coefficients (horizontal axis). Positive and negative values on each axis represent direction and strength of correlation of individual metabolites with the response variable represented by that axis. Although severe RSV and asthma classification are loosely correlated (Figure 4A and Figure 5A), the individual metabolites most strongly associated with RSV severity and asthma classification are distinct.
Metabolites 12 00178 g0a3
Metabolites 12 00178 g0a4 550
Figure A4. Urinary metabolites are more predictive of future asthma occurrence. Receiver–operator characteristic (ROC) curves for asthma prediction, using either (A) the top positive and negative urinary metabolites (lysine and xanthosine) or (B) the top positive and negative plasma metabolites (uracil and citrate/isocitrate) in the asthma PLS model. The predictor for each ROC curve was calculated by subtracting the negative-coefficient metabolite value from the positive-coefficient metabolite value. The values in parentheses for the areas under the curve (AUC) represent 95% confidence intervals determined by bootstrap validation.
Figure A4. Urinary metabolites are more predictive of future asthma occurrence. Receiver–operator characteristic (ROC) curves for asthma prediction, using either (A) the top positive and negative urinary metabolites (lysine and xanthosine) or (B) the top positive and negative plasma metabolites (uracil and citrate/isocitrate) in the asthma PLS model. The predictor for each ROC curve was calculated by subtracting the negative-coefficient metabolite value from the positive-coefficient metabolite value. The values in parentheses for the areas under the curve (AUC) represent 95% confidence intervals determined by bootstrap validation.
Metabolites 12 00178 g0a4
Metabolites 12 00178 g0a5 550
Figure A5. Succinate levels are decreased in urine samples of RSV patients. Vertical axis represents normalized relative metabolite level. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. p values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Figure A5. Succinate levels are decreased in urine samples of RSV patients. Vertical axis represents normalized relative metabolite level. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. p values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Metabolites 12 00178 g0a5
Metabolites 12 00178 g0a6 550
Figure A6. Metabolites in the hexosamine and sialic acid biosynthesis pathways are decreased in plasma of RSV patients. Vertical axes represent normalized relative metabolite levels. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. p values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Figure A6. Metabolites in the hexosamine and sialic acid biosynthesis pathways are decreased in plasma of RSV patients. Vertical axes represent normalized relative metabolite levels. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. p values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Metabolites 12 00178 g0a6

References

  1. Scheltema, N.M.; Gentile, A.; Lucion, F.; Nokes, D.J.; Munywoki, P.K.; Madhi, S.A.; Groome, M.J.; Cohen, C.; Moyes, J.; Thorburn, K.; et al. Global Respiratory Syncytial Virus-Associated Mortality in Young Children (RSV GOLD): A Retrospective Case Series. Lancet Glob. Health 2017, 5, e984–e991. [Google Scholar] [CrossRef] [Green Version]
  2. Leader, S.; Kohlhase, K. Respiratory Syncytial Virus-Coded Pediatric Hospitalizations, 1997 to 1999. Pediatr. Infect. Dis. J. 2002, 21, 629–632. [Google Scholar] [CrossRef]
  3. Hall, C.B.; Weinberg, G.A.; Iwane, M.K.; Blumkin, A.K.; Edwards, K.M.; Staat, M.A.; Auinger, P.; Griffin, M.R.; Poehling, K.A.; Erdman, D.; et al. The Burden of Respiratory Syncytial Virus Infection in Young Children. N. Engl. J. Med. 2009, 360, 588–598. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Priante, E.; Cavicchiolo, M.E.; Baraldi, E. RSV Infection and Respiratory Sequelae. Minerva Pediatr. 2018, 70, 623–633. [Google Scholar] [CrossRef] [PubMed]
  5. Gupta, R.; Leimanis, M.L.; Adams, M.; Bachmann, A.S.; Uhl, K.L.; Bupp, C.P.; Hartog, N.L.; Kort, E.J.; Olivero, R.; Comstock, S.S.; et al. Balancing Precision versus Cohort Transcriptomic Analysis of Acute and Recovery Phase of Viral Bronchiolitis. Am. J. Physiol. Lung Cell. Mol. Physiol. 2021, 320, L1147–L1157. [Google Scholar] [CrossRef]
  6. Shommu, N.S.; Jenne, C.N.; Blackwood, J.; Martin, D.-A.; Joffe, A.R.; Eccles, R.; Brindle, M.; Khanafer, I.; Vogel, H.J.; Thompson, G.C. The Use of Metabolomics and Inflammatory Mediator Profiling Provides a Novel Approach to Identifying Pediatric Appendicitis in the Emergency Department. Sci. Rep. 2018, 8, 4083. [Google Scholar] [CrossRef]
  7. Daniluk, U.; Daniluk, J.; Kucharski, R.; Kowalczyk, T.; Pietrowska, K.; Samczuk, P.; Filimoniuk, A.; Kretowski, A.; Lebensztejn, D.; Ciborowski, M. Untargeted Metabolomics and Inflammatory Markers Profiling in Children With Crohn’s Disease and Ulcerative Colitis—A Preliminary Study. Inflamm. Bowel Dis. 2019, 25, 1120–1128. [Google Scholar] [CrossRef] [PubMed]
  8. Zhou, B.; Lou, B.; Liu, J.; She, J. Serum Metabolite Profiles as Potential Biochemical Markers in Young Adults with Community-Acquired Pneumonia Cured by Moxifloxacin Therapy. Sci. Rep. 2020, 10, 4436. [Google Scholar] [CrossRef] [Green Version]
  9. Ahmed, W.; White, I.R.; Wilkinson, M.; Johnson, C.F.; Rattray, N.; Kishore, A.K.; Goodacre, R.; Smith, C.J.; Fowler, S.J. Breath and Plasma Metabolomics to Assess Inflammation in Acute Stroke. Sci. Rep. 2021, 11, 21949. [Google Scholar] [CrossRef]
  10. Chen, R.; Mias, G.I.; Li-Pook-Than, J.; Jiang, L.; Lam, H.Y.K.; Chen, R.; Miriami, E.; Karczewski, K.J.; Hariharan, M.; Dewey, F.E.; et al. Personal Omics Profiling Reveals Dynamic Molecular and Medical Phenotypes. Cell 2012, 148, 1293–1307. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Rodríguez-Morató, J.; Pozo, Ó.J.; Marcos, J. Targeting Human Urinary Metabolome by LC–MS/MS: A Review. Bioanalysis 2018, 10, 489–516. [Google Scholar] [CrossRef]
  12. Stewart, C.J.; Mansbach, J.M.; Ajami, N.J.; Petrosino, J.F.; Zhu, Z.; Liang, L.; Camargo, C.A.; Hasegawa, K. Serum Metabolome Is Associated with the Nasopharyngeal Microbiota and Disease Severity among Infants with Bronchiolitis. J. Infect. Dis. 2019, 219, 2005–2014. [Google Scholar] [CrossRef]
  13. Zheng, J.; Zhang, L.; Johnson, M.; Mandal, R.; Wishart, D.S. Comprehensive Targeted Metabolomic Assay for Urine Analysis. Anal. Chem. 2020, 92, 10627–10634. [Google Scholar] [CrossRef]
  14. Roberts, L.D.; Souza, A.L.; Gerszten, R.E.; Clish, C.B. Targeted Metabolomics. Curr. Protoc. Mol. Biol. 2012, 98, 30.2.1–30.2.24. [Google Scholar] [CrossRef] [PubMed]
  15. Schwaiger-Haber, M.; Stancliffe, E.; Arends, V.; Thyagarajan, B.; Sindelar, M.; Patti, G.J. A Workflow to Perform Targeted Metabolomics at the Untargeted Scale on a Triple Quadrupole Mass Spectrometer. ACS Meas. Sci. Au 2021, 1, 35–45. [Google Scholar] [CrossRef] [PubMed]
  16. Lunt, S.Y.; Muralidhar, V.; Hosios, A.M.; Israelsen, W.J.; Gui, D.Y.; Newhouse, L.; Ogrodzinski, M.; Hecht, V.; Xu, K.; Acevedo, P.N.M.; et al. Pyruvate Kinase Isoform Expression Alters Nucleotide Synthesis to Impact Cell Proliferation. Mol. Cell 2015, 57, 95–107. [Google Scholar] [CrossRef] [Green Version]
  17. Teoh, S.T.; Ogrodzinski, M.P.; Ross, C.; Hunter, K.W.; Lunt, S.Y. Sialic Acid Metabolism: A Key Player in Breast Cancer Metastasis Revealed by Metabolomics. Front. Oncol. 2018, 8, 174. [Google Scholar] [CrossRef]
  18. Yu, L.; Teoh, S.T.; Ensink, E.; Ogrodzinski, M.P.; Yang, C.; Vazquez, A.I.; Lunt, S.Y. Cysteine Catabolism and the Serine Biosynthesis Pathway Support Pyruvate Production during Pyruvate Kinase Knockdown in Pancreatic Cancer Cells. Cancer Metab. 2019, 7, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Ogrodzinski, M.P.; Teoh, S.T.; Lunt, S.Y. Targeting Subtype-Specific Metabolic Preferences in Nucleotide Biosynthesis Inhibits Tumor Growth in a Breast Cancer Model. Cancer Res. 2021, 81, 303–314. [Google Scholar] [CrossRef]
  20. Evans, E.D.; Duvallet, C.; Chu, N.D.; Oberst, M.K.; Murphy, M.A.; Rockafellow, I.; Sontag, D.; Alm, E.J. Predicting Human Health from Biofluid-Based Metabolomics Using Machine Learning. Sci. Rep. 2020, 10, 17635. [Google Scholar] [CrossRef]
  21. Wittig, H.J.; Glaser, J. The Relationship between Bronchiolitis and Childhood Asthma; a Follow-up Study of 100 Cases of Bronchiolitis. J. Allergy 1959, 30, 19–23. [Google Scholar] [CrossRef]
  22. Stein, R.T.; Sherrill, D.; Morgan, W.J.; Holberg, C.J.; Halonen, M.; Taussig, L.M.; Wright, A.L.; Martinez, F.D. Respiratory Syncytial Virus in Early Life and Risk of Wheeze and Allergy by Age 13 Years. Lancet Lond. Engl. 1999, 354, 541–545. [Google Scholar] [CrossRef]
  23. Holt, P.G.; Sly, P.D. Viral Infections and Atopy in Asthma Pathogenesis: New Rationales for Asthma Prevention and Treatment. Nat. Med. 2012, 18, 726–735. [Google Scholar] [CrossRef] [PubMed]
  24. Di Cicco, M.; D’Elios, S.; Peroni, D.G.; Comberiati, P. The Role of Atopy in Asthma Development and Persistence. Curr. Opin. Allergy Clin. Immunol. 2020, 20, 131–137. [Google Scholar] [CrossRef]
  25. Leteurtre, S.; Duhamel, A.; Grandbastien, B.; Proulx, F.; Cotting, J.; Gottesman, R.; Joffe, A.; Wagner, B.; Hubert, P.; Martinot, A.; et al. Daily Estimation of the Severity of Multiple Organ Dysfunction Syndrome in Critically Ill Children. CMAJ Can. Med. Assoc. J. J. Assoc. Med. Can. 2010, 182, 1181–1187. [Google Scholar] [CrossRef] [Green Version]
  26. Slater, A.; Shann, F.; Pearson, G. Paediatric Index of Mortality (PIM) Study Group PIM2: A Revised Version of the Paediatric Index of Mortality. Intensive Care Med. 2003, 29, 278–285. [Google Scholar] [CrossRef] [PubMed]
  27. Pollack, M.M.; Patel, K.M.; Ruttimann, U.E. PRISM III: An Updated Pediatric Risk of Mortality Score. Crit. Care Med. 1996, 24, 743–752. [Google Scholar] [CrossRef]
  28. Adamko, D.J.; Saude, E.; Bear, M.; Regush, S.; Robinson, J.L. Urine Metabolomic Profiling of Children with Respiratory Tract Infections in the Emergency Department: A Pilot Study. BMC Infect. Dis. 2016, 16, 439. [Google Scholar] [CrossRef] [Green Version]
  29. Turi, K.N.; Romick-Rosendale, L.; Gebretsadik, T.; Watanabe, M.; Brunwasser, S.; Anderson, L.J.; Moore, M.L.; Larkin, E.K.; Peebles, R.S.; Hartert, T.V. Using Urine Metabolomics to Understand the Pathogenesis of Infant Respiratory Syncytial Virus (RSV) Infection and Its Role in Childhood Wheezing. Metabolomics 2018, 14, 135. [Google Scholar] [CrossRef] [PubMed]
  30. Biswas, A.; Shukla, A.; Chaudhary, S.K.; Santhosh, R.; Jeyakanthan, J.; Sekar, K. Structural Studies of a Hyperthermophilic Thymidylate Kinase Enzyme Reveal Conformational Substates along the Reaction Coordinate. FEBS J. 2017, 284, 2527–2544. [Google Scholar] [CrossRef] [PubMed]
  31. Mackenzie, N.C.W.; Huesa, C.; Rutsch, F.; MacRae, V.E. New Insights into NPP1 Function: Lessons from Clinical and Animal Studies. Bone 2012, 51, 961–968. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Torrecilla, A.; Marques, A.F.P.; Buscalioni, R.D.; Oliveira, J.M.A.; Teixeira, N.A.; Atencia, E.A.; Sillero, M.A.G.; Sillero, A. Metabolic Fate of AMP, IMP, GMP and XMP in the Cytosol of Rat Brain: An Experimental and Theoretical Analysis. J. Neurochem. 2001, 76, 1291–1307. [Google Scholar] [CrossRef]
  33. Zimmermann, A.G.; Gu, J.J.; Laliberté, J.; Mitchell, B.S. Inosine-5′-Monophosphate Dehydrogenase: Regulation of Expression and Role in Cellular Proliferation and T Lymphocyte Activation. Prog. Nucleic Acid Res. Mol. Biol. 1998, 61, 181–209. [Google Scholar] [CrossRef]
  34. Calise, S.J.; Abboud, G.; Kasahara, H.; Morel, L.; Chan, E.K.L. Immune Response-Dependent Assembly of IMP Dehydrogenase Filaments. Front. Immunol. 2018, 9, 2789. [Google Scholar] [CrossRef] [PubMed]
  35. Wu, P.; Hartert, T.V. Evidence for a Causal Relationship between Respiratory Syncytial Virus Infection and Asthma. Expert Rev. Anti-Infect. Ther. 2011, 9, 731–745. [Google Scholar] [CrossRef] [Green Version]
  36. Lemanske, R.F. The Childhood Origins of Asthma (COAST) Study. Pediatr. Allergy Immunol. 2002, 13, 38–43. [Google Scholar] [CrossRef] [PubMed]
  37. Yamada, Y.; Matsumoto, K.; Hashimoto, N.; Saikusa, M.; Homma, T.; Yoshihara, S.; Saito, H. Effect of Th1/Th2 Cytokine Pretreatment on RSV-Induced Gene Expression in Airway Epithelial Cells. Int. Arch. Allergy Immunol. 2011, 154, 185–194. [Google Scholar] [CrossRef] [PubMed]
  38. Clasquin, M.F.; Melamud, E.; Rabinowitz, J.D.; Clasquin, M.F.; Melamud, E.; Rabinowitz, J.D. LC-MS Data Processing with MAVEN: A Metabolomic Analysis and Visualization Engine. In Current Protocols in Bioinformatics, Current Protocols in Bioinformatics; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2012; pp. 14.11.1–14.11.23. ISBN 978-0-471-25095-1. [Google Scholar]
  39. 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]
  40. Thévenot, E.A.; Roux, A.; Xu, Y.; Ezan, E.; Junot, C. Analysis of the Human Adult Urinary Metabolome Variations with Age, Body Mass Index, and Gender by Implementing a Comprehensive Workflow for Univariate and OPLS Statistical Analyses. J. Proteome Res. 2015, 14, 3322–3335. [Google Scholar] [CrossRef] [PubMed]
  41. The National Assessment of Educational Progress (NAEP); Third Expert Panel on the Diagnosis and Management of Asthma. Section 3, Component 1: Measures of Asthma Assessment and Monitoring; National Heart, Lung, and Blood Institute (US): Bethesda, MD, USA, 2007.
Metabolites 12 00178 g001 550
Figure 1. Plasma and urinary metabolite profiles reflect RSV severity. (A) Score plot of principal components analysis (PCA) performed using plasma and urinary metabolite data for all samples. Horizontal and vertical axes represent 3rd and 4th principal component scores, respectively. (B) Score plot of partial least squares (PLS) regression performed using metabolite data from control and RSV baseline samples, with RSV severity used as response variable. Horizontal and vertical axes represent 1st and 2nd PLS component scores, respectively.
Figure 1. Plasma and urinary metabolite profiles reflect RSV severity. (A) Score plot of principal components analysis (PCA) performed using plasma and urinary metabolite data for all samples. Horizontal and vertical axes represent 3rd and 4th principal component scores, respectively. (B) Score plot of partial least squares (PLS) regression performed using metabolite data from control and RSV baseline samples, with RSV severity used as response variable. Horizontal and vertical axes represent 1st and 2nd PLS component scores, respectively.
Metabolites 12 00178 g001
Metabolites 12 00178 g002 550
Figure 2. Sialic acid pathway metabolites in plasma are correlated with RSV severity. (A) PLS coefficients of plasma metabolites in the RSV model. Metabolites with positive coefficients are positively correlated with RSV severity, i.e., enriched in RSV samples; metabolites with negative coefficients are negatively correlated with RSV severity, i.e., enriched in controls relative to RSV samples. (B,C) Box plots of top plasma metabolites with (B) positive or (C) negative coefficients. Vertical axes represent normalized relative metabolite levels. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. P values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Figure 2. Sialic acid pathway metabolites in plasma are correlated with RSV severity. (A) PLS coefficients of plasma metabolites in the RSV model. Metabolites with positive coefficients are positively correlated with RSV severity, i.e., enriched in RSV samples; metabolites with negative coefficients are negatively correlated with RSV severity, i.e., enriched in controls relative to RSV samples. (B,C) Box plots of top plasma metabolites with (B) positive or (C) negative coefficients. Vertical axes represent normalized relative metabolite levels. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. P values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Metabolites 12 00178 g002
Metabolites 12 00178 g003 550
Figure 3. Nucleotides in urine are correlated with RSV severity. (A) PLS coefficients of urinary metabolites in the RSV model. Metabolites with positive coefficients are positively correlated with RSV severity, i.e., enriched in RSV samples; metabolites with negative coefficients are negatively correlated with RSV severity, i.e., enriched in controls relative to RSV samples. (B,C). Box plots of top urinary metabolites with (B) positive or (C) negative coefficients. Vertical axes represent normalized relative metabolite levels. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. p values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Figure 3. Nucleotides in urine are correlated with RSV severity. (A) PLS coefficients of urinary metabolites in the RSV model. Metabolites with positive coefficients are positively correlated with RSV severity, i.e., enriched in RSV samples; metabolites with negative coefficients are negatively correlated with RSV severity, i.e., enriched in controls relative to RSV samples. (B,C). Box plots of top urinary metabolites with (B) positive or (C) negative coefficients. Vertical axes represent normalized relative metabolite levels. Lower and upper ranges of the boxes represent 1st and 3rd quartiles, respectively; the center lines represent median values, and whiskers indicate maximum and minimum values. p values displayed underneath the sample class labels are calculated by the Mann–Whitney U test; significant (p < 0.05) values are highlighted in red.
Metabolites 12 00178 g003
Metabolites 12 00178 g004 550
Figure 4. Patient weight and plasma protein are inversely correlated with RSV severity. (A) Covariate loadings on the 1st principal component of the covariates PCA model. Covariates with positive loadings have higher values in samples with positive 1st principal component scores, and vice versa; hence, covariates on the same side of the loading plot are positively correlated with each other, and covariates on the opposite side are inversely correlated. (B–E) Score plots all showing the 1st and 2nd principal components in the covariates-only PCA model, with sample data points color-coded according to (B) RSV severity, (C) patient weight, (D) plasma protein quantification, (E) classification of the patient as asthmatic within the two-year follow-up period after hospital visit.
Figure 4. Patient weight and plasma protein are inversely correlated with RSV severity. (A) Covariate loadings on the 1st principal component of the covariates PCA model. Covariates with positive loadings have higher values in samples with positive 1st principal component scores, and vice versa; hence, covariates on the same side of the loading plot are positively correlated with each other, and covariates on the opposite side are inversely correlated. (B–E) Score plots all showing the 1st and 2nd principal components in the covariates-only PCA model, with sample data points color-coded according to (B) RSV severity, (C) patient weight, (D) plasma protein quantification, (E) classification of the patient as asthmatic within the two-year follow-up period after hospital visit.
Metabolites 12 00178 g004
Metabolites 12 00178 g005 550
Figure 5. Future development of asthma is predictable from combined plasma and urinary metabolite profiles. (A) Stacked bar plot showing breakdown of numbers of control, moderate RSV and severe RSV patients who were diagnosed as asthmatic in the two-year follow-up period post hospital visit. (B) Score plot of PLS model using metabolite data from control and RSV 72 h timepoint samples, with asthma classification as response variable. Triangles represent patients classified as asthmatic while rectangles represent patients not classified as asthmatic. Horizontal and vertical axes represent 1st and 2nd PLS component scores, respectively. (C) PLS coefficients of metabolites in the RSV model. Metabolites with positive coefficients are positively correlated with asthma classification, i.e., enriched in patients classified as asthmatic; conversely, metabolites with negative coefficients are enriched patients not classified as asthmatic. Plasma and urinary metabolite names include “PL.” and “UR.”, respectively, and the bar plot is further color-coded (pink bars for plasma, yellow bars for urinary metabolites).
Figure 5. Future development of asthma is predictable from combined plasma and urinary metabolite profiles. (A) Stacked bar plot showing breakdown of numbers of control, moderate RSV and severe RSV patients who were diagnosed as asthmatic in the two-year follow-up period post hospital visit. (B) Score plot of PLS model using metabolite data from control and RSV 72 h timepoint samples, with asthma classification as response variable. Triangles represent patients classified as asthmatic while rectangles represent patients not classified as asthmatic. Horizontal and vertical axes represent 1st and 2nd PLS component scores, respectively. (C) PLS coefficients of metabolites in the RSV model. Metabolites with positive coefficients are positively correlated with asthma classification, i.e., enriched in patients classified as asthmatic; conversely, metabolites with negative coefficients are enriched patients not classified as asthmatic. Plasma and urinary metabolite names include “PL.” and “UR.”, respectively, and the bar plot is further color-coded (pink bars for plasma, yellow bars for urinary metabolites).
Metabolites 12 00178 g005
Table
Table 1. Patient cohort information for samples used in this study.
Table 1. Patient cohort information for samples used in this study.
Range (n)%MeanMedianSt. dev.
Age (months)0.5–7301002.242.01.60
  Controls0.5–610 3.203.01.83
  RSV0–720 2.002.01.52
Gender
  Female 1756.7
  Controls 413.3
  RSV 1343.3
Birth weight (kg)
  Controls 10 3.653.70.40
  RSV 20 3.103.20.84
Hospital LOS 1 (days) 20 14.9311.209.37
  Severe 15 15.7314.588.63
  Moderate 5 12.536.7212.13
PICU 2 LOS 1 (days) 20 10.6010.045.56
  Severe 15 11.2811.414.86
  Moderate 5 8.565.537.58
Respiratory Support (days)
  Mechanical ventilation 15 8.957.084.30
  Mechanical ventilation (non-invs.) 1 2.002.0N/A
  High flow oxygen 4 0.910.191.56
1 LOS: length of stay; 2 PICU: pediatric intensive care unit.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Teoh, S.T.; Leimanis-Laurens, M.L.; Comstock, S.S.; Winters, J.W.; Vandenbosch, N.L.; Prokop, J.W.; Bachmann, A.S.; Lunt, S.Y.; Rajasekaran, S. Combined Plasma and Urinary Metabolomics Uncover Metabolic Perturbations Associated with Severe Respiratory Syncytial Viral Infection and Future Development of Asthma in Infant Patients. Metabolites 2022, 12, 178. https://doi.org/10.3390/metabo12020178

AMA Style

Teoh ST, Leimanis-Laurens ML, Comstock SS, Winters JW, Vandenbosch NL, Prokop JW, Bachmann AS, Lunt SY, Rajasekaran S. Combined Plasma and Urinary Metabolomics Uncover Metabolic Perturbations Associated with Severe Respiratory Syncytial Viral Infection and Future Development of Asthma in Infant Patients. Metabolites. 2022; 12(2):178. https://doi.org/10.3390/metabo12020178

Chicago/Turabian Style

Teoh, Shao Thing, Mara L. Leimanis-Laurens, Sarah S. Comstock, John W. Winters, Nikita L. Vandenbosch, Jeremy W. Prokop, André S. Bachmann, Sophia Y. Lunt, and Surender Rajasekaran. 2022. "Combined Plasma and Urinary Metabolomics Uncover Metabolic Perturbations Associated with Severe Respiratory Syncytial Viral Infection and Future Development of Asthma in Infant Patients" Metabolites 12, no. 2: 178. https://doi.org/10.3390/metabo12020178

APA Style

Teoh, S. T., Leimanis-Laurens, M. L., Comstock, S. S., Winters, J. W., Vandenbosch, N. L., Prokop, J. W., Bachmann, A. S., Lunt, S. Y., & Rajasekaran, S. (2022). Combined Plasma and Urinary Metabolomics Uncover Metabolic Perturbations Associated with Severe Respiratory Syncytial Viral Infection and Future Development of Asthma in Infant Patients. Metabolites, 12(2), 178. https://doi.org/10.3390/metabo12020178

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