Next Article in Journal
A Deep Mining Strategy for Peptide Rapid Identification in Lactobacillus reuteri Based on LC–MS/MS Integrated with FBMN and De Novo Sequencing
Previous Article in Journal
The Effects of Maternal Nutrient Restriction during Mid to Late Gestation with Realimentation on Fetal Metabolic Profiles in the Liver, Skeletal Muscle, and Blood in Sheep
Previous Article in Special Issue
Association of Indoleamine 2,3-Dioxygenase (IDO) Activity with Outcome after Cardiac Surgery in Adult Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Untargeted Blood Lipidomics Analysis in Critically Ill Pediatric Patients with Ventilator-Associated Pneumonia: A Pilot Study

1
Analytical Chemistry Laboratory, Department of Chemical Engineering, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
2
Biomic Auth, Bioanalysis and Omics Lab, Centre for Interdisciplinary Research of Aristotle University of Thessaloniki, Innovation Area of Thessaloniki, 57001 Thermi, Greece
3
Department of Chemistry, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
4
Infectious Diseases Unit, 3rd Department Pediatrics, School of Medicine, Faculty of Health Sciences, Aristotle University of Thessaloniki, Hippokration General Hospital, 54642 Thessaloniki, Greece
5
Pediatric Intensive Care Unit, Hippokration General Hospital, 54642 Thessaloniki, Greece
6
Laboratory of Forensic Medicine and Toxicology, School of Medicine, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
*
Author to whom correspondence should be addressed.
Metabolites 2024, 14(9), 466; https://doi.org/10.3390/metabo14090466
Submission received: 25 July 2024 / Revised: 18 August 2024 / Accepted: 19 August 2024 / Published: 23 August 2024

Abstract

:
This study aims to explore the diagnostic potential of blood lipid profiles in suspected ventilator-associated pneumonia (VAP). Early detection of VAP remains challenging for clinicians due to subjective clinical criteria and the limited effectiveness of current diagnostic tests. Blood samples from 20 patients, with ages between 6 months and 15 years, were collected at days 1, 3, 6, and 12, and an untargeted lipidomics analysis was performed using a Ultra high Pressure Liquid Chromatography hyphenated with High Resolution Mass Spectrometry UPLC-HRMS (TIMS-TOF/MS) platform. Patients were stratified based on modified pediatric clinical pulmonary index score (mCPIS) into high (mCPIS ≥ 6, n = 12) and low (mCPIS < 6, n = 8) VAP suspicion groups. With the untargeted lipid profiling, we were able to identify 144 lipid species from different lipid groups such as glycerophospholipids, glycerolipids, and sphingolipids, in the blood of children with VAP. Multivariate and univariate statistical analyses revealed a distinct distribution of blood lipid profiles between the studied groups, indicating the potential utility of lipid biomarkers in discriminating VAP presence. Additionally, specific lipids were associated with pharyngeal culture results, notably the presence of Klebsiella pneumoniae and Staphylococcus aureus, underscoring the importance of lipid profiling in identifying the microbial etiology of VAP.

1. Introduction

Ventilator-associated pneumonia (VAP) is one of the most common nosocomial infections in critically ill children in pediatric intensive care units (PICUs). VAP is a severe complication, leading to increased morbidity, mortality, and healthcare costs [1]. Timely and accurate diagnosis is essential for early initiation of appropriate antimicrobial treatment [2,3,4]. Clinicians face significant challenges in diagnosing VAP early on, primarily due to the reliance on subjective clinical criteria and the limited accuracy of existing diagnostic tests. While most VAP definitions are based on adult ICU research, pediatric patients require distinct classifications due to their unique physiologies and pathologies [1]. Few studies have focused on the risk factors for pediatric VAP, and there is still no gold standard for a VAP definition and diagnosis in children. Current diagnostics rely on clinical assessment, radiological imaging, and microbiological cultures of respiratory samples. These methods exhibit high inter-variability and only moderate sensitivity and specificity. Microbiological confirmation can take several days, resulting in overtreatment with antibiotics until the specific pathogen is identified [5]. New diagnostic tools—faster, more reliable and robust, and less invasive—are required to differentiate patients with different severities or characteristics of VAP, enabling potential therapeutic targets and personalized treatment strategies.
Omics approaches can offer new advantages in diagnostics as complex pathologies can be more accurately diagnosed via the readout of the whole metabolic phenotype than with the current one-dimensional, single-marker diagnostic procedures. Metabolomics represents a rapidly evolving discipline focused on the holistic analysis and study of small endogenous molecules (molecular weight, MW < 1500 Da), commonly referred to as “metabolites” [6]. Lipidomics, a subfield of metabolomics, involves the study of the lipidome or the total lipid content within a cell, organ, or biological system [7]; however, only a few studies have been conducted in pneumonia and, more particularly, in children cohorts. The majority of the studies have been focused on the analysis of small, volatile organic compounds (VOCs) in the exhaled breath of hospitalized ICU adult patients or in strain liquid cultures using Gas Chromatography coupled to different Mass Spectrometry detectors, GC-MS or GC-qTOFMS analytical platforms [8,9,10,11,12,13,14]. The VOC’s metabolic profile was found to be differentiated between various disease states, including VAP; thus, breath samples are found to be an ideal, non-invasive matrix for VAP diagnosis and/or prognosis. A few studies explored urine and serum metabolic profiles of the same population groups using Reversed Phase Liquid Chromatography RPLC-qTOFMS and Nuclear Magnetic Resonance 1D 1H NMR, respectively [15,16]. Only Laiakis et al. investigated possible alterations in the plasma and urine metabolic profiles of Gambian children with VAP compared to healthy controls using LC-qTOF MS [17].
The objective of this study was to investigate the diagnostic potential of blood lipidomic biomarker signatures for identifying VAP in hospitalized, pediatric, critically ill patients. By examining lipidomic profiles, this study aimed to provide a more reliable and objective method for early VAP diagnosis that could potentially improve patient management and ultimately prognosis.

2. Materials and Methods

2.1. Study Population

The study design and patient population enrolled in the study are mentioned in detail by Sdougka et al. [18]. Briefly, children aged from 1 month to 16 years with clinical suspicion of VAP and on mechanical ventilation for at least 48 h were eligible for the study, and informed consent was asked by parents/caregivers.More information about study population can be found in Table 1.
Enrollment criteria included different levels of clinical suspicion. Step one was a purulent respiratory or positive bronchial aspirate culture and initiation of antimicrobial agents; an increased oxygen requirement (>20%) with fever, hypothermia, leukocytosis, or leukopenia, and an initiation of antimicrobial agents; and new lung infiltrates on radiology with at least two additional criteria such as fever, hypothermia, increased oxygen requirement, purulent secretions, abnormal white blood cell count, and elevated C-reactive protein CRP. The exclusions criteria included a lack of consent, unlikely survival beyond 48 h, pregnancy in adolescent females, and a body weight < 3 kg.
On day 1 (step 2), all the patients enrolled were further assessed using the modified CPIS (mCPIS) tool and categorized into high (mCPIS ≥ 6) and low (mCPIS < 6) VAP suspicion groups. Twelve children were allocated to high VAP suspicion and eight to low VAP suspicion, according to the mCPIS.
Pharyngeal cultures were obtained on days 1, 3, 6, and 12. Pathogen identification was performed using the Vitek® 2 automated system (Biomerieux, Marcyl’Etoile, France).

2.2. Chemicals and Reagents

LC-MS grade acetonitrile (ACN), methanol (MeOH), and 2-propanol (IPA) were procured from Sigma-Aldrich (St. Louis, MO, USA), while methyl-tert-butyl-ether (MTBE; ≥99%), ammonium formate, and formic acid (FA, >99% LC-MS grade) were supplied by Chem-Lab NV (Zedelgem, Belgium). Pure water (H2O, 18.2 MΩ cm−1) was purified using a Milli-Q device (Millipore Purification System, Merck Darmstadt, Darmstadt, Germany).

2.3. Blood Samples

Blood samples of the children were collected at four specific time points: days 1, 3, 6, and 12 of hospitalization at the ICU following the onset of the event (VAP suspicion). They were then divided into respective aliquots for all analyses performed (hematological, biochemical, metabolomics/lipidomics), and stored at −70 °C until the LC-HRMS analysis.

2.4. Sample Preparation

Blood samples were left to thaw at 4 °C prior to the analysis. Fifty microliters of sample were diluted with 700 μL of ice-cold MeOH:MTBE, 1:3 v/v (−20 °C), vortexed for 5 min, and centrifuged at 25,200× g for 20 min at 4 °C. Sixty hundreds microliters of the organic phase were then transferred to an Eppendorf tube and evaporated to dryness under a vacuum (SpeedVac, Eppendorf Austria GmbH, Vienna, Austria). The residue was reconstituted with 150 μL of H2O:ACN:IPA (1:1:3, v/v), and the clear supernatant was transferred into a LC-MS/MS glass vial, placed in an autosampler at 6 °C, and subjected to an RPLC-HRMS analysis.
Quality control (QC) samples were prepared by mixing equal volumes of all analyzed clinical samples (20 μL). Phenotypic QC samples were also prepared for high and low VAP suspicion groups by mixing equal volumes of the respective samples. All QC samples were treated as mentioned above. From the initial QC samples, diluted QCs (1:2, 1:4 and 1:8, v/v) were also prepared for the evaluation of the dilution factor of the detected features in both positive and negative ionization mode.

2.5. RPLC-TOFMS Analysis

An analysis was performed in a UHPLC Elute chromatography system coupled to a tims (time of flight) TOF mass spectrometer (Bruker, Billerica, MA, USA), operating in positive and negative electrospray ionisation ESI ionization mode. For chromatographic separation of all analytes, an ACQUITY UPLC CSH C18 (2.1 × 100 mm, 1.7 μm) column was used, equipped with an ACQUITY UPLC Van-Guard pre-column of the same packing material, (Waters Ltd., Elstree, UK). The binary mobile phase system consisted of (a) ACN:H2O (60:40, v/v), 10 mM ammonium formate and 0.1% formic acid, and (b) IPA:ACN (90:10, v/v), 0.1% formic acid. The gradient elution program was set as follows: 60–57% A (0.0–2.0 min), 57–50% A (2.0–2.1 min), 50–46% A (2.1–12.0 min), 46–30% A (12.0–12.1 min), 30–1%A (12.1–18 min), 1–60% A (18.0–18.1 min), and 60% A (18.1–20.0 min). Flow rate and column temperature were 0.3 mL/min and 55 °C, respectively. Before and after each injection, the needle was conditioned with both a weak and strong wash cycle (1500 μL each), consisting of ACN:H2O (60/40, v/v) and IPA:ACN (90:10, v/v), respectively. Injection volume was 3 μL for both positive and negative ESI mode.
Data-dependent acquisition (DDA) experiments were performed for the analyses of all blood samples. The parameters in ESI included capillary voltage for both modes at ±4.5 kV, and dessolvation parameters were 200 °C, 10 L/min and 2 bar for dry temperature, dry gas, and nebulizer gas, respectively. Dynamic spectra acquisition was applied for auto MS/MS analysis with a frequency of minimum 6 and maximum 10 Hz. Regarding mass fragmentation, collision energy was set at 20 V for m/z below 100, 30 V for m/z 100–1000, and 40 V for m/z 1000–2000 m/z. Deflection delta and transfer time were set to 80 V and 54 μs, respectively, and collision RF (Radio Frequency) was set to 1100 Vpp. Sodium formate 10 mM was used as a calibrant mixture and was directly infused into the MS at a 30 µL/h and acquired in the first 0.5 min of each injection.
Quality control samples (QCs) were used throughout the analytical batch in order to evaluate the system’s repeatability and analytical performance. QC samples were analyzed seven times at the beginning of the analytical batch for system equilibration and then at every other 10 blood samples. Phenotype QCs and diluted QCs were also analyzed after the system’s equilibration.

2.6. Data Handling and Statistical Analysis

All obtained raw data underwent recalibration based on sodium formate clusters using a vendor’s software “Data Analysis” (version 5.3, Bruker, Bremen, Germany) and then was converted to a .mzML file format using MSConvert v. 3.0.20133 (ProteoWizard). Chromatographic peak detection, retention time alignment, and feature grouping were performed with XCMS (version 3.2.0) in the R programming environment, where extracted feature values were reported as peak areas. Normalization of the raw data was performed using quality control (QC) samples, support vector regression (QC-SVRC), and a radial basis function kernel to correct the intra-batch effect as proposed by Kuligowski et al. [19]. Features further assessed were present in more than 80% of real samples in each group, and their coefficient variation (CV) in the QC sample was less than 0.3%. Zero or missing values were replaced with half of the respective minimum feature value of the respective group, and all raw data were log-transformed. For univariate statistical analysis, the advanced software GraphPad Prism 8.0.1 for Windows (GraphPad Software, La Jolla, CA, USA) and SPSS Statistics for Windows, version 29.0.2.0 (IBM Corp., Armonk, NY, USA) were used. Data were checked for normal distribution based on the Shapiro–Wilk test (significance level, alpha = 0.05), and respective parametric or non-parametric tests (Mann-Whitney or Kruskal-Wallis) were applied, followed by false discovery rate (FDR) correction (Benjamini, Krieger, and Yekutieli). The statistical significance threshold value was defined as p-value < 0.05, and variation between lipids was illustrated as a box and whisker plot. The area under the curve–receiver operating characteristic (AUC–ROC) curves and the logarithm base two-fold change (log2FC) were calculated setting the initial/control value to low, or Klebsiella pneumoniae group, and the final value to high, or Staphylococcus aureus group.
For multivariate statistical analysis, unsupervised principal component analysis (PCA) and orthogonal partial least squares-discriminate analysis (OPLS-DA) in UV scaling were performed using SIMCA P+ 13.0 software (Umetrics, Malmö, Sweden). The validity of the constructed models was based on their R2X, R2Y, and Q2Y values, permutation plots, and p-value of CV-ANOVA (cross validation, R2X: fraction of the variation in X explained by the model, R2Y: total sum of variation in Y explained by the model, and Q2Y: goodness of prediction). Different visualization plots, as VIP and “S-plot” graphs, were also assessed for the evaluation of possible biomarkers, with lipids meeting the criteria of p > |0.05| and p (corr) > |0.5| were considered statistically significant.

2.7. Lipid Annotation

Lipid species were firstly identified using Lipostar2 (v.2.0.2, Molecular Discovery Ltd., Hertfordshire, UK) [20] with the LIPID MAPS mass spectra library (version September, 2021) [20] and MSDial software using MS-DIAL LipidBlast library (version 68) [21]. Raw or converted .abf files were imported and aligned with the default settings. Peak picking was executed automatically by the Savitzky–Golay algorithm, employing window size of seven, degree of two, one multi-pass iteration, and a minimum signal-to-noise (S/N) ratio of three. The alignment of lipid profiles was performed by configuring the mass accuracy to 10 ppm and the retention time (RT) tolerance to 0.2 min. Other parameters, such as “Retain lipids with isotopic pattern” and “Retain lipids with MS/MS”, were employed to retain only the compounds exhibiting isotopic patterns and MS/MS in order to be used for possible annotation. Only the lipids with a quality rating of 3–4 stars were accepted and used further for statistical analysis. MS Dial parameters were set as follows: mass accuracy of 0.01 Da and 0.05 Da in MS1 and MS2, respectively, for spectral centroiding, and the identification score cutoff was set to 70%. Identified lipids and possible annotations of both software integrated, and their accuracy, was verified by plotting the retention time of specific lipid species against their Kendrick mass defect relative to the hydrogen base using an in-house script written in R. To gain a thorough understanding of the retention time mapping process for different lipid (sub) classes, one can refer to the detailed description provided by Lange et al. [22].

3. Results

3.1. Lipid Profiling by LC-TOF-MS

Raw data processing by XCMS revealed a wealth of features; 2835 ion signals in positive and 1446 in negative ion mode were considered for further analysis after quality control filtering. An untargeted lipidomics workflow facilitated the identification of 144 lipid species, which were then subjected to univariate and multivariate statistical analyses. The identified lipid species in blood plasma comprised 2.7% fatty acyls (4 carnitines), 61.1% glycerophospholipids (88 species), 18% glycerolipids (26 species), 15% sphingolipids (21 species), and 3.4% sterols (5 species). Figure 1 illustrates the subclasses of lipids species quantified in the blood of children with VAP. A detailed table, with the identifications of detected lipids, the chemical formulas, the monoisotopic masses, and the retention time, is presented in Table S1.

3.2. Quality of the Analytical System

Initially, to assess the quality, repeatability, and robustness of the analytical system, QC samples were evaluated (Figure 2A) after QC-SVRC normalization, filtering, and exclusion of the ions with % RSD > 30. The score plot demonstrated the clustering of the QC samples, indicating the stability of the analytical system providing a map of how the samples relate to each other. The first component explains 36% of the variation, and the second component explains 44%.

3.3. Multivariate Analysis of Lipidomics and VAP Suspicion: High or Low

To investigate any possible lipidomic-based signature associated with the presence of VAP, a multivariate statistical analysis was performed based on the semi-quantified lipid species. The unsupervised PCA scores-plot model constructed could not clearly classify the studied groups. However, separation of the children according to the mCPIS to low and high VAP groups could be achieved by partial least squares-discriminant analysis (PLS-DA). OPLS-DA analysis between the studied groups revealed strong differentiation between children with different CPIS scores for VAP as illustrated in Figure 2B. A combined multivariate and univariate analysis revealed significant perturbations, particularly in response to high and low VAP suspicion groups. All the statistically significant lipids were observed at decreased levels in the blood of children with low VAP mCPIS score. Phosphatidylcholines (LPC 20:5, PC 32:1, PC 32:2, PC 34:1, PC 34:2, PC 34:2, PC 34:3, PC 35:2, PC 36:1, PC 37:3, PC 38:6) sphingomyelins (SM 32:1, SM 34:0, SM 38:1, SM 40:1, SM 40:2, SM 41:1, SM 42:2), and glycerolipids (DG 36:3, TG 46:1, TG 48:1, TG 48:2, TG 48:3, TG 50:4, TG 50:3) were among the main lipid species with their abundances significantly different between the studied groups.

3.4. Multivariate Analysis for Lipidomics and Isolation of Targeted Respiratory Bacterial Pathogens

In addition, we assessed correlations between the blood lipid content of children with VAP suspicion and the bacterial species isolated in pharyngeal swab cultures. The pharyngeal cultures of children showed the presence of pathogens: K. pneumoniae in four children, S. aureus in four children, and Acinetobacter spp., Stenotrophomonas maltophilia, and Candida spp. in one child each. The rest of the cohort had a negative culture, suggesting that VAP suspicion was clinically, rather than microbiologically, driven. Multivariate statistical analysis was also performed to assess if the bacterial species correlated with the lipid profile in children with VAP suspicion. The PCA scores plot model did not reveal any clear classification of the samples according to the detected bacteria; however, OPLS-DA (Figure 2C) displayed distinct clustering of the studied groups, indicating the correlation of the lipid profile with the presence of different bacterial strains. Table 2 presents the cross-validation parameters of the aforementioned scores plot models. Based on univariate and multivariate analyses, nine lipids were identified to contribute to this classification. Three diglyceride species, namely, DG 16:1_18:0, DG 18:0_18:1, DG 18:1_18:1, eight phosphatidylcholines PC 33:2, PC 34:3, PC 36:4, PC 36:4, PC 37:3, PC 38:4, PC O-34:0 and PC O-38:4, four sphingomyelins SM 34:1, SM 36:2, SM 38:2, SM 40:1, five triacylglycerols TG 48:0, TG 50:1, TG 52:5, TG 54:5, TG 54:6, and the cholesterol ester CE 20:4, demonstrated a significant impact in the discrimination of the two groups, as all of them, with the exception of SM 34:1, were found to be elevated in the blood plasma of children with cultures positive for S. aureus.
Statistically differentiated lipid species, derived from all comparisons, are provided in Table 3 and Table 4, along with the respective p-values < 0.05, VIP scores, Log2FC, and their lower and upper bounds of the 95% confidence intervals (CI). Figure 3A,B illustrate box plots of all altered metabolites found in the respective comparisons.

4. Discussion

Lipids were found to participate in a variety of processes intimately connected with the pathogenesis of pneumonia and related infectious diseases. The lungs are essential respiratory organs, and lipids play crucial roles in the development and progression of lung diseases. Inflammation and infections induce a plethora of alterations in lipid metabolism that could initially reduce inflammation or fight infection. To the best of our knowledge, this is the first time that the contribution of blood lipids in VAP diagnosis is assessed in critically ill patients. Only a few studies have assessed the relationship between lipid metabolism and pneumonia, and in most of the cases, these studies have focused on community-acquired pneumonia. In these cases, distinct lipid profiles consistently appeared in numerous patients, accompanied by varying airway inflammatory statuses and disease severities. Similarly, in our study, there were differences in lipid metabolism in patients with high or low VAP probability.
In this study, we characterized the blood lipidome of 20 children with VAP suspicion, including those with high- and low-grades of suspicion according to predefined clinical, laboratory, and imaging criteria (mCPIS score). In addition, given the significant number of children whom tested positive for Klebsiella and Staphylococcus aureus, these groups were further analyzed to identify potential lipid markers associated with these bacteria. Blood untargeted lipidomic analysis revealed extensive variations in phospholipids, sphingolipids, and glycerolipids levels in all these patient groups. In particular, this study identified certain phospholipids, sphingolipids, and glycerolipids at significantly increased abundances in children with low or high VAP suspicion according to mCPIS classification score.
Cholesterol, sphingomyelins (SM), and phosphatidylcholines (PC) contribute to the creation of the immunological synapse. Additionally, lipids are crucial in activating macrophages, supporting NK (natural killer) cell functionality, and guiding the differentiation and activity of T and B effector cells. Lipid-based signaling is known to regulate various processes linked to inflammation and damage related to infections, such as apoptosis, autophagy, and pulmonary fibrosis [23].
Sphingolipids are essential components found in the mucous secreted by the alveolar epithelium, which helps protect lung tissue from invading pathogens [24]. Studies have indicated that Mycoplasma pneumoniae infection in the lungs can trigger the production of autoantibodies against glycosphingolipids, suggesting a potential role of sphingolipids in promoting lung inflammation [25,26]. This raises the possibility that sphingolipids could contribute to pulmonary inflammation during infection. Sphingomyelins are essential cellular components that participate in numerous processes such as cell division, cell proliferation, and autophagy. They also help balance pro- and anti-inflammatory lipids and regulate immune responses in lung tissues [27,28]. Inflammatory lung diseases like pneumonia and chronic obstructive pulmonary disease involve altered sphingomyelin metabolism. This alteration triggers the activation of inflammatory cells and the release of mediators, worsening tissue damage and inflammation [29]. Moreover, abnormalities in sphingomyelin metabolism are associated with pulmonary fibrosis, a serious lung condition characterized by an excessive growth of connective tissue that impairs respiratory function. In pulmonary fibrosis tissues, levels of sphingomyelin are heightened, potentially promoting fibroblast proliferation and collagen synthesis, contributing to the progression of pulmonary fibrosis [30]. In our study, all the SM species, with the exception of SM 42:2, presented AUC values > 0.70, indicating that SMs might be a good indicator for the discrimination of children according to VAP severity.
In order to identify significant lipid markers that could aid in VAP diagnosis, the entire set of blood samples was utilized for comprehensive multivariate and univariate statistical analyses to evaluate potential biomarkers. Further longitudinal comparisons among the different time points, along with dynamic monitoring of the lipid abundance, were conducted to provide additional insights and enhance the predictive accuracy for disease severity. Lipid markers with an AUC greater than 0.8 were subjected to further evaluation to identify any trends related to the timing of the sample collection (Figure 4). We found that all triglycerides (TGs) displayed a statistically significant upward trend over time (comparing day 1 to day 12) in the high VAP suspicion group (p-value < 0.05). This trend was not observed in the low suspicion group, indicating differences between these two groups. This finding could suggest that monitoring TG levels over time could be a valuable tool in assessing VAP diagnosis in pediatric patients. Such dynamic changes in lipid abundance could potentially serve as robust biomarkers for VAP suspicion and diagnosis, ultimately contributing to better clinical outcomes.
The relationship between phospholipid metabolism and inflammatory diseases has attracted increasing attention in recent years. Our study found that blood phosphocholines and lyso-phosphocholine concentrations were significantly increased in children with high VAP suspicion and significantly decreased in children with positive K. pneumoniae cultures. Fluctuations in the abundance of these lipids lead to the change of glycerophospholipid metabolism, which is closely related to the occurrence and development of many diseases including inflammation [31]. The results of the present study suggest that blood phospholipids not only have a strong VAP biomarker potential with AUC ranging from 0.69 for PC 36:1 and PC 32:1 to 0.81 for PC 32:2, but fluctuations in its levels may also be closely related to bacterial infections. Different pathogens may induce different phospholipid metabolism alterations [32].
Inflammation and infections can cause shifts in lipid metabolism. While initially beneficial for combating infection or reducing inflammation, prolonged changes may increase the risk of developing atherosclerosis. Common effects include decreased levels of HDL (high-density lipoprotein) cholesterol and increased triglycerides. Inflammation leads to an increase in lipoprotein (a) levels due to an enhanced apolipoprotein (a) synthesis. LDL levels often decrease, but there is an increase in low density LDL particles due to the exchange of triglycerides from triglyceride-rich lipoproteins to LDL, followed by triglyceride hydrolysis. In addition to altering blood lipid levels, inflammation also impairs lipoprotein function. LDL becomes more susceptible to oxidation because HDL’s ability to prevent LDL oxidation diminishes. Furthermore, inflammation adversely affects several steps in the reverse cholesterol transport pathway. These lipids and lipoprotein abnormalities are consistently observed with greater severity in inflammatory diseases. Treatment that reduces inflammation typically restores the lipid profile towards normal. These changes in lipids and lipoproteins are integral to the innate immune response and likely play a crucial role in host protection during infection and inflammation [33]. TGs and DGs were observed as significant lipids in both studied groups, showing the ability to serve as indicators in both VAP severity and differentiation of K. pneumoniae from S. aureus bacteria.
In conclusion, there were significant differences in lipid metabolism in critically ill children with low and high VAP suspicion. The small sample size and the lack of an external validation cohort could be a limitation of the study; thus, the predictive performance of these lipids and their usefulness as biomarkers requires further verification. The study underscores the potential of lipidomics-based approaches to improve VAP diagnosis in PICU, necessitating further exploration of omics data integration with microbiological, molecular, and immune response studies. Such integrated investigations are crucial for elucidating the underlying biochemical pathways associated with VAP pathogenesis and diagnosis, offering insights into potential therapeutic targets and personalized treatment strategies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo14090466/s1, Table S1: Peak areas of annotated lipids.

Author Contributions

Conceptualization, C.V., E.R. and E.I.; Methodology, C.V., O.B. and M.S.; Software, O.B.; Formal analysis, C.V.; Investigation, C.V., G.T. and E.I.; Resources, M.S. (Maria Sdougka), G.T. and H.G.; Data curation, C.V., A.F. and M.S. (Maria Simitsopoulouand); Writing—original draft, C.V. and O.B.; Writing—review & editing, C.V., E.R., G.T., H.G. and E.I.; Visualization, O.B.; Supervision, E.R. and H.G.; Project administration, E.I. All authors have read and agreed to the published version of the manuscript.

Funding

This research project was supported by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the “2nd Call for H.F.R.I. Research Projects to support Post-Doctoral Researchers” (Project Number: 98270).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Hippokration General Hospital of Thessaloniki, protocol code: 58291/7.12.2020, date of approval: 14 April 2021.

Informed Consent Statement

Informed consent was obtained from all subjects enrolled in the study.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Acknowledgments

We thank the patients enrolled in the study, their parents/legal guardians, and the personnel of the PICU caring for these critically ill children. This work was financially supported by the Hellenic Foundation for Research and Innovation (H.F.R.I.) under the “2nd Call for H.F.R.I. Research Projects to support Post-Doctoral Researchers” (Project Number: 98270). HG is supported in the framework of the project “Biomic_AUTh, Center of Excellence in Metabolomics research (BiACEM)” (Project number: 101079370), which is implemented under the Action “Horizon Europe Twinning grant”.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kalamitsou, S.; Violaki, A.; Iosifidis, E.; Avramidou, V.; Mantzafleri, P.-E.; Karaiskou, E.; Karantaglis, N.; Sdougka, M. Ventilator associated pneumonia (VAP) in children: A diagnostic challenge. Signa Vitae 2023, 19, 6–19. [Google Scholar] [CrossRef]
  2. da Silva, P.S.L.; de Aguiar, V.E.; de Carvalho, W.B.; Fonseca, M.C.M. Value of clinical pulmonary infection score in critically ill children as a surrogate for diagnosis of ventilator-associated pneumonia. J. Crit. Care 2014, 29, 545–550. [Google Scholar] [CrossRef] [PubMed]
  3. Ahmed, W.; Burke, A.; Gotts, N.; Van Oort, P.; Bannard-Smith, J.; Morton, B.; Welters, I.; Mcmullan, R.; Goodacre, R.; Dark, P.; et al. Metabolomics of blood plasma from critically ill patients with suspected ventilator-associated pneumonia. Eur. Respir. J. 2023, 62, PA2962. [Google Scholar] [CrossRef]
  4. Wildman, E.; Mickiewicz, B.; Vogel, H.J.; Thompson, G.C. Metabolomics in pediatric lower respiratory tract infections and sepsis: A literature review. Pediatr. Res. 2022, 93, 492–502. [Google Scholar] [CrossRef] [PubMed]
  5. Bajo-Fernández, M.; Souza-Silva, É.A.; Barbas, C.; Rey-Stolle, M.F.; García, A. GC-MS-based metabolomics of volatile organic compounds in exhaled breath: Applications in health and disease. A review. Front. Mol. Biosci. 2024, 10, 1295955. [Google Scholar] [CrossRef]
  6. Fiehn, O. Metabolomics—The link between genotypes and phenotypes. Plant Mol. Biol. 2022, 48, 155–171. [Google Scholar] [CrossRef]
  7. Belhaj, M.R.; Lawler, N.G.; Hoffman, N.J. Metabolomics and Lipidomics: Expanding the Molecular Landscape of Exercise Biology. Metabolites 2021, 11, 151. [Google Scholar] [CrossRef] [PubMed]
  8. van Oort, P.M.P.; Nijsen, T.; Weda, H.; Knobel, H.; Dark, P.; Felton, T.; Rattray, N.J.W.; Lawal, O.; Ahmed, W.; Portsmouth, C.; et al. BreathDx—Molecular analysis of exhaled breath as a diagnostic test for ventilator–associated pneumonia: Protocol for a European multicentre observational study. BMC Pulm. Med. 2017, 17, 1. [Google Scholar] [CrossRef]
  9. Filipiak, W.; Żuchowska, K.; Marszałek, M.; Depka, D.; Bogiel, T.; Warmuzińska, N.; Bojko, B. GC-MS profiling of volatile metabolites produced by Klebsiella pneumoniae. Front. Mol. Biosci. 2022, 9, 1019290. [Google Scholar] [CrossRef]
  10. Xu, W.; Zou, X.; Ding, Y.; Zhang, J.; Zheng, L.; Zuo, H.; Yang, M.; Zhou, Q.; Liu, Z.; Ge, D.; et al. Rapid screen for ventilator associated pneumonia using exhaled volatile organic compounds. Talanta 2023, 253, 124069. [Google Scholar] [CrossRef]
  11. Schnabel, R.; Fijten, R.; Smolinska, A.; Dallinga, J.; Boumans, M.-L.; Stobberingh, E.; Boots, A.; Roekaerts, P.; Bergmans, D.; van Schooten, F.J. Analysis of volatile organic compounds in exhaled breath to diagnose ventilator-associated pneumonia. Sci. Rep. 2015, 5, 17179. [Google Scholar] [CrossRef] [PubMed]
  12. van Oort, P.M.P.; de Bruin, S.; Weda, H.; Knobel, H.H.; Schultz, M.J.; Bos, L.D. Exhaled Breath Metabolomics for the Diagnosis of Pneumonia in Intubated and Mechanically-Ventilated Intensive Care Unit (ICU)-Patients. Int. J. Mol. Sci. 2017, 18, 449. [Google Scholar] [CrossRef] [PubMed]
  13. Lawal, O.; Muhamadali, H.; Ahmed, W.M.; White, I.R.; Nijsen, T.M.E.; Goodacre, R.; Fowler, S.J. Headspace volatile organic compounds from bacteria implicated in ventilator-associated pneumonia analysed by TD-GC/MS. J. Breath Res. 2018, 12, 026002. [Google Scholar] [CrossRef]
  14. Mourani, P.M.; Sontag, M.K. Ventilator-Associated Pneumonia in Critically Ill Children. Pediatr. Clin. 2017, 64, 1039–1056. [Google Scholar] [CrossRef]
  15. Jongers, B.; Hotterbeekx, A.; Bielen, K.; Vervliet, P.; Boddaert, J.; Lammens, C.; Fransen, E.; Baggerman, G.; Covaci, A.; Goossens, H.; et al. Identification of Potential Urinary Metabolite Biomarkers of Pseudomonas aeruginosa Ventilator-Associated Pneumonia. Biomark. Insights 2022, 17. [Google Scholar] [CrossRef]
  16. Antcliffe, D.; Jiménez, B.; Veselkov, K.; Holmes, E.; Gordon, A.C. Metabolic Profiling in Patients with Pneumonia on Intensive Care. EBioMedicine 2017, 18, 244–253. [Google Scholar] [CrossRef] [PubMed]
  17. Laiakis, E.C.; Morris, G.A.J.; Fornace, A.J.; Howie, S.R.C. Metabolomic Analysis in Severe Childhood Pneumonia in The Gambia, West Africa: Findings from a Pilot Study. PLoS ONE 2010, 5, e12655. [Google Scholar] [CrossRef]
  18. Sdougka, M.; Simitsopoulou, M.; Volakli, E.; Violaki, A.; Georgopoulou, V.; Ftergioti, A.; Roilides, E.; Iosifidis, E. Evaluation of Five Host Inflammatory Biomarkers in Early Diagnosis of Ventilator-Associated Pneumonia in Critically Ill Children: A Prospective Single Center Cohort Study. Antibiotics 2023, 12, 921. [Google Scholar] [CrossRef]
  19. Kuligowski, J.; Sánchez-Illana, Á.; Sanjuán-Herráez, D.; Vento, M.; Quintás, G. Intra-batch effect correction in liquid chromatography-mass spectrometry using quality control samples and support vector regression (QC-SVRC). Analyst 2015, 140, 7810–7817. [Google Scholar] [CrossRef]
  20. Goracci, L.; Tortorella, S.; Tiberi, P.; Pellegrino, R.M.; Di Veroli, A.; Valeri, A.; Cruciani, G. Lipostar, a Comprehensive Platform-Neutral Cheminformatics Tool for Lipidomics. Anal. Chem. 2017, 89, 6257–6264. [Google Scholar] [CrossRef]
  21. Tsugawa, H.; Cajka, T.; Kind, T.; Ma, Y.; Higgins, B.; Ikeda, K.; Kanazawa, M.; VanderGheynst, J.; Fiehn, O.; Arita, M. MS-DIAL: Data-independent MS/MS deconvolution for comprehensive metabolome analysis. Nat. Methods 2015, 12, 523–526. [Google Scholar] [CrossRef] [PubMed]
  22. Lange, M.; Fedorova, M. Evaluation of lipid quantification accuracy using HILIC and RPLC MS on the example of NIST® SRM® 1950 metabolites in human plasma. Anal. Bioanal. Chem. 2020, 412, 3573–3584. [Google Scholar] [CrossRef] [PubMed]
  23. Spadaro, F.; Cecchetti, S.; Fantuzzi, L. Macrophages and Phospholipases at the Intersection between Inflammation and the Pathogenesis of HIV-1 Infection. Int. J. Mol. Sci. 2017, 18, 1390. [Google Scholar] [CrossRef] [PubMed]
  24. Wheelock, C.E.; Goss, V.M.; Balgoma, D.; Nicholas, B.; Brandsma, J.; Skipp, P.J.; Snowden, S.; Burg, D.; D’Amico, A.; Horvath, I.; et al. Application of ’omics technologies to biomarker discovery in inflammatory lung diseases. Eur. Respir. J. 2013, 42, 802–825. [Google Scholar] [CrossRef] [PubMed]
  25. Sharma, L.; Prakash, H. Sphingolipids Are Dual Specific Drug Targets for the Management of Pulmonary Infections: Perspective. Front. Immunol. 2017, 8, 378. [Google Scholar] [CrossRef]
  26. Nishimura, M.; Saida, T.; Kuroki, S.; Kawabata, T.; Obayashi, H.; Saida, K.; Uchiyama, T. Post-infectious encephalitis with anti-galactocerebroside antibody subsequent to Mycoplasma pneumoniae infection. J. Neurol. Sci. 1996, 140, 91–95. [Google Scholar] [CrossRef]
  27. Harvald, E.B.; Olsen, A.S.B.; Færgeman, N.J. Autophagy in the light of sphingolipid metabolism. Apoptosis 2015, 20, 658–670. [Google Scholar] [CrossRef]
  28. Ghidoni, R.; Caretti, A.; Signorelli, P. Role of Sphingolipids in the Pathobiology of Lung Inflammation. Mediat. Inflamm. 2015, 2015, 487508. [Google Scholar] [CrossRef]
  29. Wu, W.; Li, Z.; Wang, Y.; Huang, C.; Zhang, T.; Zhao, H. Advances in metabolomics of chronic obstructive pulmonary disease. Chin. Med. J. Pulm. Crit. Care Med. 2023, 1, 223–230. [Google Scholar] [CrossRef]
  30. Shi, X.; Chen, Y.; Shi, M.; Gao, F.; Huang, L.; Wang, W.; Wei, D.; Shi, C.; Yu, Y.; Xia, X.; et al. The novel molecular mechanism of pulmonary fibrosis: Insight into lipid metabolism from reanalysis of single-cell RNA-seq databases. Lipids Health Dis. 2024, 23, 98. [Google Scholar] [CrossRef]
  31. Boggara, M.B.; Krishnamoorti, R. Partitioning of Nonsteroidal Antiinflammatory Drugs in Lipid Membranes: A Molecular Dynamics Simulation Study. Biophys. J. 2010, 98, 586–595. [Google Scholar] [CrossRef] [PubMed]
  32. Zhang, C.; Wang, K.; Yang, L.; Liu, R.; Chu, Y.; Qin, X.; Yang, P.; Yu, H. Lipid metabolism in inflammation-related diseases. Analyst 2018, 143, 4526–4536. [Google Scholar] [CrossRef] [PubMed]
  33. Feingold, K.R.; Grunfeld, C. The Effect of Inflammation and Infection on Lipids and Lipoproteins. In Endotext; Feingold, K.R., Anawalt, B., Blackman, M.R., Boyce, A., Chrousos, G., Corpas, E., de Herder, W.W., Dhatariya, K., Dungan, K., Hofland, J., et al., Eds.; MDText.com, Inc.: South Dartmouth, MA, USA, 2000. [Google Scholar]
Figure 1. Lipid subclasses quantified in the blood of children with VAP based on lipid analyses. Abbreviations: CAR: Carnitines, SM: Sphingomyelins, LPC: Monoacylglycerophosphocholines, PC: Diacylglycerophosphocholines, PE: Diacylglycerophosphoethanolamines, CE: Cholesterol Esters, DG: Diglycerides, MG Monoglycerides, and TG: Triglycerides.
Figure 1. Lipid subclasses quantified in the blood of children with VAP based on lipid analyses. Abbreviations: CAR: Carnitines, SM: Sphingomyelins, LPC: Monoacylglycerophosphocholines, PC: Diacylglycerophosphocholines, PE: Diacylglycerophosphoethanolamines, CE: Cholesterol Esters, DG: Diglycerides, MG Monoglycerides, and TG: Triglycerides.
Metabolites 14 00466 g001
Figure 2. (A) PCA-X score plots of all blood samples collected at different time points and QC samples. QC samples are presented with a red color and are clustered together, and OPLS−DA score plot showing the classification of (B) high and low VAP suspicion and (C) different bacteria strain isolation in pharyngeal cultures based on the blood lipidome.
Figure 2. (A) PCA-X score plots of all blood samples collected at different time points and QC samples. QC samples are presented with a red color and are clustered together, and OPLS−DA score plot showing the classification of (B) high and low VAP suspicion and (C) different bacteria strain isolation in pharyngeal cultures based on the blood lipidome.
Metabolites 14 00466 g002
Figure 3. (A) Box plots of all significantly altered plasma lipids in children with high (red) and low (green) suspicion in VAP. Data are plotted as median with range. (B) Box plots of all significantly altered plasma lipids in children with Klebsiella (blue) and S. aureus (orange) bacteria according to pharyngeal culture results. Data are plotted as median with range. (phoscholines PCs, sphyngomyelins SMs, triacylgleycerolsin TGs, Diacylglycerols (DGs, and cholester ester CE).
Figure 3. (A) Box plots of all significantly altered plasma lipids in children with high (red) and low (green) suspicion in VAP. Data are plotted as median with range. (B) Box plots of all significantly altered plasma lipids in children with Klebsiella (blue) and S. aureus (orange) bacteria according to pharyngeal culture results. Data are plotted as median with range. (phoscholines PCs, sphyngomyelins SMs, triacylgleycerolsin TGs, Diacylglycerols (DGs, and cholester ester CE).
Metabolites 14 00466 g003aMetabolites 14 00466 g003b
Figure 4. Box plots of significantly altered PCs (phosphocholines) and TGs (triacylglycerols) lipids that presented an AUC value > 0.8 at the two studied groups at the different time points. Asterisks indicate a p-value < 0.05 compared to time point 1.
Figure 4. Box plots of significantly altered PCs (phosphocholines) and TGs (triacylglycerols) lipids that presented an AUC value > 0.8 at the two studied groups at the different time points. Asterisks indicate a p-value < 0.05 compared to time point 1.
Metabolites 14 00466 g004
Table 1. Summarizes patients’ characteristics enrolled in the study. Patients’ characteristics adapted from [18].
Table 1. Summarizes patients’ characteristics enrolled in the study. Patients’ characteristics adapted from [18].
TotalHigh VAP Suspicion
(mCPIS ≥ 6)
Low VAP Suspicion
(mCPIS < 6)
n = 20n = 12n = 8
Age, months93 (6–184)24.5 (6–141)129 (28–184)
male, n (%)13 (65)9 (75)5 (63)
PRISM score12.1 (5.45)13.1 (6.49)11 (4.14)
time to event, days9.37 (7.65)8.73 (7.43)8.73 (7.43)
duration of mechanical ventilation, days26.3 (14.1)24.4 (16.2)29 (11)
Hospitalization, days67.3 (47.4)61.7 (46.0)75.8 (51.2)
mortality, n (%)2 (10)2 (17)0
Table 2. Characteristics of the constructed unsupervised and supervised models; UV scale was applied in all models.
Table 2. Characteristics of the constructed unsupervised and supervised models; UV scale was applied in all models.
ModelTypeNR2XR2YQ2CV ANOVA
Samples_QCPCA-X1310.580 0.262
LOW vs HIGH VAP SuspicionOPLS-DA740.2670.9650.4365.2 × 10−6
K. pneumoniae vs. S. aureusOPLS-DA300.2860.9870.5284.8 × 10−3
Table 3. Lipids with statistical significance according to their VAP score (low and high suspicion).
Table 3. Lipids with statistical significance according to their VAP score (low and high suspicion).
Blood LOW HIGH Suspecion
Lipidsp Valuep Value Adj. WaistVIPAUCLog2FC95% Lower CI95% Upper CI
CE 18:36.99 × 10−33.32 × 10−21.070.65−0.580.520.77
DG 36:32.34 × 10−31.62 × 10−21.410.73−0.690.600.85
LPC 20:51.10 × 10−24.79 × 10−21.180.71−0.620.570.81
PC 32:12.02 × 10−31.62 × 10−20.980.69−0.520.560.80
PC 32:22.17 × 10−51.46 × 10−31.800.81−0.690.700.91
PC 34:13.76 × 10−32.20 × 10−21.350.75−0.270.630.86
PC 34:21.61 × 10−31.46 × 10−21.450.75−0.330.620.85
PC 34:21.04 x 10−24.70 × 10−21.650.77−0.330.660.87
PC 34:31.29 × 10−24.96 × 10−21.320.68−0.960.570.80
PC 35:23.77 × 10−45.38 × 10−31.510.74−0.490.620.84
PC 36:11.22 × 10−24.96 × 10−20.930.69−0.310.560.79
PC 37:33.14 × 10−31.96 × 10−21.500.76−0.370.650.87
PC 38:61.38 × 10−31.46 × 10−21.730.80−0.870.690.90
SM 32:11.25 × 10−24.96 × 10−21.460.74−0.600.610.84
SM 34:05.33 × 10−32.66 × 10−21.440.74−0.430.630.85
SM 38:12.80 × 10−44.65 × 10−31.680.77−0.560.660.86
SM 40:11.49 × 10−31.46 × 10−21.710.79−0.620.680.89
SM 40:21.43 × 10−31.46 × 10−21.350.71−0.410.590.82
SM 41:11.32 × 10−43.29 × 10−31.730.79−0.590.680.89
SM 42:25.33 × 10−32.66 × 10−21.320.69−0.440.570.81
TG 46:16.01 × 10−52.00 × 10−31.150.81−1.090.690.90
TG 48:11.80 × 10−43.60 × 10−31.550.78−0.850.660.87
TG 48:22.92 × 10−51.46 × 10−31.520.82−1.390.720.91
TG 48:32.18 × 10−31.62 × 10−21.220.82−2.080.720.91
TG 50:42.43 × 10−31.62 × 10−21.580.77−1.320.650.87
TG 54:34.80 × 10−32.66 × 10−21.640.75−0.880.640.86
Table 4. Lipids with statistical significance according to the bacteria strains K. pneumoniae and S. aureus.
Table 4. Lipids with statistical significance according to the bacteria strains K. pneumoniae and S. aureus.
Bacteria Strain
Lipidsp ValueVIPAUCLog2FC95% Lower CI95% Upper CI
CE 20:43.45 × 10−21.10.68−0.380.470.89
PC 33:21.27 × 10−21.70.75−0.660.580.91
PC 34:31.47 × 10−21.20.71−0.570.500.88
PC 36:41.76 × 10−21.40.73−0.470.550.89
PC 36:42.46 × 10−21.20.73−0.410.530.91
PC 37:36.23 × 10−31.90.76−0.730.600.90
PC 38:42.97 × 10−21.40.72−0.380.530.89
PC O-34:01.55 × 10−21.80.76−0.400.590.91
PC O-38:43.29 × 10−21.70.72−0.590.530.89
SM 34:11.04 × 10−21.70.675.760.480.87
SM 36:21.39 × 10−21.70.76−0.350.570.91
SM 38:21.61 × 10−21.70.75−0.320.560.91
SM 40:13.15 × 10−21.50.72−0.440.530.88
TG 48:08.22 × 10−31.90.77−0.270.610.92
TG 50:17.98 × 10−32.00.76−0.600.580.91
TG 52:54.95 × 10−21.50.72−0.760.510.88
TG 54:52.31 × 10−31.90.81−0.760.640.95
TG 54:62.98 × 10−32.00.80−0.360.630.93
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Virgiliou, C.; Begou, O.; Ftergioti, A.; Simitsopoulou, M.; Sdougka, M.; Roilides, E.; Theodoridis, G.; Gika, H.; Iosifidis, E. Untargeted Blood Lipidomics Analysis in Critically Ill Pediatric Patients with Ventilator-Associated Pneumonia: A Pilot Study. Metabolites 2024, 14, 466. https://doi.org/10.3390/metabo14090466

AMA Style

Virgiliou C, Begou O, Ftergioti A, Simitsopoulou M, Sdougka M, Roilides E, Theodoridis G, Gika H, Iosifidis E. Untargeted Blood Lipidomics Analysis in Critically Ill Pediatric Patients with Ventilator-Associated Pneumonia: A Pilot Study. Metabolites. 2024; 14(9):466. https://doi.org/10.3390/metabo14090466

Chicago/Turabian Style

Virgiliou, Christina, Olga Begou, Argyro Ftergioti, Maria Simitsopoulou, Maria Sdougka, Emmanuel Roilides, Georgios Theodoridis, Helen Gika, and Elias Iosifidis. 2024. "Untargeted Blood Lipidomics Analysis in Critically Ill Pediatric Patients with Ventilator-Associated Pneumonia: A Pilot Study" Metabolites 14, no. 9: 466. https://doi.org/10.3390/metabo14090466

APA Style

Virgiliou, C., Begou, O., Ftergioti, A., Simitsopoulou, M., Sdougka, M., Roilides, E., Theodoridis, G., Gika, H., & Iosifidis, E. (2024). Untargeted Blood Lipidomics Analysis in Critically Ill Pediatric Patients with Ventilator-Associated Pneumonia: A Pilot Study. Metabolites, 14(9), 466. https://doi.org/10.3390/metabo14090466

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