Disease Markers
Volume 2018 (2018), Article ID 6964529, 16 pages
https://doi.org/10.1155/2018/6964529

Research Article

Novel Biomarker Candidates for Febrile Neutropenia in Hematological Patients Using Nontargeted Metabolomics

Marika Lappalainen,1,2 Jenna Jokkala,³ Auni Juutilainen,² Sari Hämäläinen,⁴ Irma Koivula,⁴ Esa Jantunen,2,4,5 Kati Hanhineva,³ and Kari Pulkki6,7,8

¹Department of Internal Medicine, Central Hospital of Central Finland, Jyväskylä, Finland
²Institute of Clinical Medicine/Internal Medicine, University of Eastern Finland, Kuopio, Finland
³Institute of Public Health and Clinical Nutrition, University of Eastern Finland, Kuopio, Finland
⁴Department of Medicine, Kuopio University Hospital, Kuopio, Finland
⁵Siun Sote-Hospital District of North Carelia, Joensuu, Finland
⁶Eastern Finland Laboratory Centre, Kuopio, Finland
⁷Laboratory Division, Turku University Hospital, Turku, Finland
⁸Clinical Chemistry, Faculty of Medicine, University of Turku, Turku, Finland

Correspondence should be addressed to Kari Pulkki; kari.pulkki@utu.fi

Received 3 November 2017; Revised 4 February 2018; Accepted 22 February 2018; Published 12 April 2018

Academic Editor: Hyundoo Hwang

Copyright © 2018 Marika Lappalainen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Background. Novel potential small molecular biomarkers for sepsis were analyzed with nontargeted metabolite profiling to find biomarkers for febrile neutropenia after intensive chemotherapy for hematological malignancies. Methods. Altogether, 85 patients were included into this prospective study at the start of febrile neutropenia after intensive chemotherapy for acute myeloid leukemia or after autologous stem cell transplantation. The plasma samples for the nontargeted metabolite profiling analysis by liquid chromatography-mass spectrometry were taken when fever rose over 38° and on the next morning. Results. Altogether, 90 differential molecular features were shown to explain the differences between patients with complicated (bacteremia, severe sepsis, or fatal outcome) and noncomplicated courses of febrile neutropenia. The most differential compounds were an androgen hormone, citrulline, and phosphatidylethanolamine PE(18:0/20:4). The clinical relevance of the findings was evaluated by comparing them with conventional biomarkers like C-reactive protein and procalcitonin. Conclusion. These results hold promise to find out novel biomarkers for febrile neutropenia, including citrulline. Furthermore, androgen metabolism merits further studies.

1. Introduction

Febrile neutropenia is a common complication in hematological patients receiving intensive chemotherapy. Although a minority of these patients develop septic shock [1, 2], sepsis is still a major cause of morbidity and mortality during the neutropenic phase [3, 4]. In these patients, life-threatening complications can develop in hours depending on the pathogen. Not only C-reactive protein (CRP) and procalcitonin (PCT) but also several other biomarkers have been explored to identify patients at risk for complicated course of febrile neutropenia [5]. CRP and PCT both have some limitations such as nonspecificity and delayed response. PCT is superior to CRP for predictive purposes and is slightly more pathogen-dependent, especially in gram-negative bacteremia [3, 6]. However, there still is a need for more rapid and accurate markers, which also could be used for de-escalation strategies of broad-spectrum antibiotics.

Nontargeted metabolite profiling or metabolomics is a hypothesis-free study approach that focuses on finding differences in metabolite profiles between study subjects contributing to the discovery of novel small-sized molecular biomarkers for disease progression or prevention [7, 8]. In previous studies, metabolomics has been used to differentiate sepsis from systemic inflammatory response syndrome [9–12] and to detect prognostic biomarkers for septic shock [11, 13–15]. There is only one previous study including metabolomics in patients with febrile neutropenia. Richter et al. [16] found twenty-one biomarker candidates including three peptides, six proteins, and six phosphatidylcholines (PC), which identified febrile neutropenic patients with proven infection from those without it.

The use of metabolomics as a tool for discovery of diagnostic markers in febrile neutropenia is a novel approach. The purpose of our study was to evaluate the differences found in a metabolic profile and to identify potentially useful biomarker candidates for further validation to recognize early increased risk of adverse outcome in hematological patients with febrile neutropenia after intensive chemotherapy.

2. Patients and Methods

2.1. Patients

Between December 2009 and November 2012, altogether 85 hematological patients treated at the adult hematology ward of the Department of Medicine, Kuopio University Hospital, and who gave written permission were included into this prospective study. The inclusion criteria were fulfilled if the patient was ≤70 years old, if the patient received intensive chemotherapy for acute myeloid leukemia (AML), or if the patient was an autologous stem cell transplant (ASCT) recipient and had febrile neutropenia (see definitions later). There were 54 males and 31 females with a median age of 61 years (18–70 years). Twenty-three patients had AML, and 62 were ASCT recipients (40 non-Hodgkin lymphoma, 19 multiple myeloma, and 3 Hodgkin lymphoma). Only the first induction course of AML patients was included in the study.

All patients were carefully followed up at the hematology ward from the start of febrile neutropenia until recovery of neutropenia. Blood pressure, oxygen saturation, respiratory frequency, heart rate, skin temperature, urine output, and fluid intake were closely followed up. Each patient was examined daily thoroughly for clinical signs and sources of infections. Serum and plasma samples for laboratory analyses were taken at the onset of neutropenic fever on day 0 (d0) and further samples on the next morning (d1). Broad-spectrum antibiotics were started as soon as the samples for blood cultures had been taken. Fifty-four patients (64%) received granulocyte colony-stimulating factor to shorten the length of the neutropenic period.

The study setting was to investigate differences in the responses of patients with and without complicated course of febrile neutropenia on two consecutive days. The noncomplicated patient group was regarded as a control group.

2.2. Data Collection

Clinical data, including the hour and date of the start of fever, possible sites of infections, and hemodynamic parameters suggesting development of septic shock were recorded on a structured data collection form. Laboratory findings, including microbiological blood culture results, were registered.

2.3. Definitions

Febrile neutropenia was defined using the criteria from IDSA (Infectious Diseases Society of America) [17]. Neutropenia was defined as a neutrophil count < 0.5 × 10⁹/L or with a predicted decrease to <0.5 × 10⁹/L during the next 48 h. Fever was defined as a single oral temperature of ≥38.3°C or a temperature of ≥38.0°C sustained over a 1-hour period.

Sepsis was defined as a syndrome in which systemic inflammatory response was present with infection. Severe sepsis was defined as sepsis with organ dysfunction. Septic shock was defined if hypoperfusion (systolic arterial , a mean arterial , or a reduction in systolic blood pressure of >40 mmHg from baseline) was present despite adequate volume resuscitation, in the absence of other causes of hypotension [18, 19].

Complicated course of febrile neutropenia was defined as a positive blood culture finding and/or development of severe sepsis or septic shock during the period from the onset of febrile neutropenia until the recovery of neutropenia.

2.4. Laboratory Measurements

The concentration of serum CRP was measured with a Konelab 60i Clinical Chemistry Analyzer (Lab systems CLD, Konelab, Helsinki, Finland) or Cobas 6000 analyzer (Hitachi, Tokyo, Japan). The intra- and interassay CV% were 2.3–4.3%. The upper reference limit of serum or plasma CRP of a healthy reference population is 3 mg/L. Plasma PCT was measured from EDTA plasma using a Cobas 6000 analyzer (Hitachi, Tokyo, Japan) with a sensitivity of 0.06 μg/L. The CVs (intra- and interassay) were 1.4% and 3.0% for 0.46 μg/L and 1.1% and 2.6% for 9.4 μg/L PCT, respectively. Blood cultures (2-3 sets including two bottles/set) were drawn immediately at the beginning of neutropenic fever (day 0), and an additional sampling was done if fever persisted for 3–5 days. They were processed using the automated blood culture system Bactec 9240 (Becton Dickinson, Sparks, MD, USA). The incubation episode for aerobic and anaerobic bottles was 7 days and for MYCO F/Lytic bottles 42 days. The plasma samples for metabolomics assays were stored frozen at −80°C until analyzed.

2.5. Nontargeted LC-MS Metabolite Profiling Analyses

The methods including sample preparation and analysis were similar as described previously [20]. In detail, the samples were prepared in 96-well plates by mixing an aliquot of 100 μL of EDTA plasma samples with 400 μL of acetonitrile (VWR International), mixed on a vortex at maximum speed for 15 s, incubated on an ice bath for 15 min to precipitate proteins, and centrifuged at 16,000 ×g for 10 min to filter the samples (0.2 μm polytetrafluoroethylene filters in a 96-well plate) and collect the supernatant. Aliquots of 2 μL were taken from at least half of the plasma samples, mixed together in a tube, and used as the quality control (QC) sample in the analysis; a solvent blank was prepared in the same manner.

The samples were analyzed by nontargeted liquid chromatography quadrupole time-of-flight mass spectrometry (LC-QTOF-MS) using the UHPLC-QTOF-MS system (Agilent Technologies, Karlsruhe, Germany), which consisted of a 1290 LC system, Jet Stream electrospray ionization (ESI), and 6540 UHD accurate-mass QTOF spectrometry. The samples were analyzed by using two different chromatographic techniques: reversed phase (RP) and hydrophilic interaction (HILIC) liquid chromatography. The sample tray was kept at 4°C during the analysis. The data acquisition software was the MassHunter Acquisition B.04.00 (Agilent Technologies). The QC samples were injected after every 12 samples and 10 samples at the beginning of the analysis. The order of the analysis of the samples was random.

In the RP technique, 2 μL of the sample solution was injected onto the column (Zorbax Eclipse XDB-C18 Rapid-Resolution HD 1.8 μm, 2.1 × 100 mm; Agilent Technologies, Palo Alto, CA, USA) and maintained at 50°C. The mobile phases, delivered at 0.4 mL/min, consisted of water (eluent A, Milli-Q purified; Millipore) and methanol (eluent B; Sigma-Aldrich), both containing 0.1% (v/v) of formic acid (Sigma-Aldrich). The following gradient profile was used: 2% → 100% B (0–10 min), 100% B (10–14.50 min), 100% → 2%B (14.50–14.51 min), and 2% B (14.51–16.50 min).

In the HILIC technique, 2 μL of the sample solution was injected onto the column (Acquity UPLC BEH Amide 1.7 μm, 2.1 × 100 mm, Waters Corporation, Milford, MA, USA) and maintained at 45°C. The mobile phases, delivered at 0.6 mL/min, consisted of 50% acetonitrile in water (v/v; eluent A) and 90% acetonitrile in water (v/v; eluent B), both containing 20 mmol/L ammonium formate, pH 3 (Sigma-Aldrich). The following gradient profile was used: 100% B (0–2.5 min), 100% → 0% B (2.5–10 min), 0% → 100% B (10.0–10.01 min), and 100% B (10.01–12.50 min).

The MS conditions after both chromatographic analyses were as follows: Jet Stream ESI source, operated in positive and negative ionization mode, a drying gas temperature 325°C, gas flow 10 L/min, a sheath gas temperature 350°C, sheath gas flow 11 L/min, nebulizer pressure 45 pounds per square inch, capillary voltage 3500 V, nozzle voltage 1000 V, fragmentor voltage 100 V, and a skimmer 45 V. For data acquisition, a 2 GHz extended dynamic range mode was used, and the instrument was set to acquire over the m/z 20–1600. Data were collected in the centroid mode at the acquisition rate of 1.67 spectra/s (i.e., 600 ms/spectrum) with an abundance threshold of 150. For the automatic data-dependent MS/MS analyses performed on the QC samples, the 4 most abundant ions were selected for fragmentation from every precursor scan cycle with a scan rate 3.33 spectra/s. These ions were excluded after two product ion spectra and released again for fragmentation after a 0.25 min hold. The precursor scan time was based on ion intensity, ending at 20,000 counts or after 200 ms in HILIC and 25,000 counts or 200 ms in RP, respectively. The product ion scan time was 200 ms. Collision energies used were 10, 20, and 40 V in subsequent assays. The continuous mass axis calibration was performed monitoring two reference ions from an infusion solution throughout the assays: m/z 121.05087300 and 922.00979800 in positive mode and m/z 112.98558700 and 966.00072500 in negative mode.

2.6. Metabolomics Data Analysis

Liquid chromatography-mass spectrometry (LC-MS) data were collected first with the “Find by Molecular Feature” algorithm (MassHunter Qualitative Analysis B.06.00, Agilent Technologies, USA). The peak collection threshold was 1000–4000 counts depending on chromatography mode, and the allowed ion species were restricted to [M+H]+ in ESI(+) and [M−H]− in ESI(−). Data files (.cef format) were exported to Mass Profiler Professional version 14.0 (Agilent Technologies) for peak alignment to create a list of potential molecular features. The molecular features were restricted to those present at least in 50% of samples within one study group, and the resulting entity list was used for feature-specific data reanalyzation back from raw data with the “Find by Formula” algorithm (MassHunter Qualitative Analysis B.06.00). For this recursive analysis, compound mass tolerance was ±15.00 ppm, retention time ± 0.1 min, and symmetric expansion value for chromatograms ± 35.0 ppm. The resulting peak data was again aligned with Mass Profiler Professional software and cleaned by filtering (metabolite features that were present in at least 80% of samples in any of the four study groups) resulting in 417, 420, 2385, and 1276 molecular features from HILIC ESI(+), HILIC ESI(−), RP ESI(+), and RP ESI(−), respectively.

The metabolite features were further subjected for statistical analysis by a pair-wise comparison of the case group consisting of patients with complicated course of febrile neutropenia and the control group consisting of noncomplicated febrile neutropenia either at day 0 or day 1 by Student’s t-test. The resulting values were adjusted for multiple comparisons by Benjamini-Hochberg false discovery rate (FDR) correction within each of the four analytical approaches [21]. Finally, the four datasets were exported into Microsoft Excel.

The data from each of the four analytical approaches were subjected to unsupervised principal component analysis (PCA) and supervised classification algorithm partial least-squares discriminant analysis (PLS-DA; SIMCA 14, Umetrics, Sweden). The data were log10-transformed and Pareto-scaled, and the model was validated by SIMCA 13 internal cross validation [22, 23]. PLS-DA illustrates the differences between the two study groups at either day 0 or day 1, separately, and gives variable importance projection (VIP) values: the larger the VIP value is, the more significant contributor the metabolite is in the model. The resulting VIP values for each metabolite were integrated in the data and used for classifying out the most important metabolite features. Due to small group size and high variability in the group, the cut-off VIP value was set to 1.5.

The metabolite features were further filtered according to an average peak and molecular to exclude small and insignificant features from the analysis resulting in a set of 1935 molecular features (144 HILIC ESI(−), 152 HILIC ESI(+), 517 RP ESI(−), and 1122 RP ESI(+)). The molecular features with a VIP value on day 0 and day and the corrected value < 0.05 on either day 0 or day 1 were considered the most significant markers. Molecular features with a VIP on day 0 and day 1 but with a noncorrected value < 0.05 also on both days were considered as the second most important molecular features. Also, molecular features with a VIP either on day 0 or on day 1 and a value < 0.05 on day 0 and day 1 were considered important and subjected for identification.

The identification of metabolites was based on MS/MS fragmentation spectra acquired in the automatic MS/MS analysis during the initial data acquisition or later on via reinjection of the samples. The spectra were compared with the in-house standard-compound library, METLIN Metabolomics Database [24], the Human Metabolome Database [25], and LIPID MAPS [26] or fragmentation patterns characteristic for certain metabolite types including phospholipids [27, 28]. Shortly, the identification of PCs was based on the presence of a protonated head group (m/z 184.073) in positive ESI mode and presence of formic acid adduct [M+COO]−, neutral loss of M+COO−60.0182, and the size of fatty acyl side-chain fragments in the negative ESI mode spectra. Phosphatidylethanolamines (PE) were identified based on the characteristic neutral loss of 141.027 Da in the positive ESI mode and fatty acyl side-chain fragments in the negative ESI mode spectra. The identification of androsterone/5α-dihydrotestosterone sulfate (ADTS/DHTS) was based on sulfate fragment SO₄H₂ of 96.9591 Da and exact molecular mass 370.1815 Da. However, as the sulfate fragment is highly dominant and the fragmentation pattern shows no other clear fragments, the metabolite was annotated as ADTS/DHTS. MS/MS fragmentation data for all the identified metabolites is presented in Table 1.

Table 1: Characteristics for the tentatively identified differential compounds in liquid chromatography-mass spectrometry analysis.

Column
ESI
MW
RT (min)
Tentative identification
ID
Identification based on MS/MS
ANOVA
Groupwise t-test Corrected value
Groupwise t-test noncorrected value

D0
D1
Fold change

Corrected
value
D0
D1
D0
D1
D0
D1
VIP D0
VIP D1

RP
−
370.1815
8.957
ADTS or DHTS
MID3559
20ev: 369.1746 [M−H]−; 369.1747 (100), 96.9591 (27), 370.1794 (7), 209.6079 (7)
1.08E − 07
8.49E − 11
6.80E − 04
0.030864
5.33E − 07
1.94E − 04
−2.37
−1.93
2.65
1.74

HILIC
+
175.0954
6.357
Citrulline
MID16
Verified with standard 10ev: 176.1043 [M+H]+; 159.0762 (100), 70.0653 (68), 113.0692 (46), 176.1043 (15)
1.17E − 04
1.05E − 06
0.77889
0.042103
0.025201
2.02E − 04
1.327
1.572
1.67
1.69

RP
−
767.5465
12.93
PE(18:0/20:4)
LMGP0201
20ev: 766.5635 [M−H]−; 766.5635 (100), 303.2348 (39), 767.5412 (34), 283.2612 (13), 59.0151 (7), 10ev: 768.5594 [M+H]+; 768.5635 (100), 627.5359 (32)
1.43E − 06
2.25E − 09
0.24644
0.016158
0.002157
4.34E − 05
−1.42
−1.63
1.71
1.98

RP
−
837.5525
11.98
PC(15:0/22:6) [M+FA]
LMGP0101
40ev: 836.546 [M+FA−H]−; 327.2326 (100), 241.2187 (42), 283.242 (39), 776.5189 (29), 168.0414 (27), 20ev: 792.5457 [M+H]+; 184.0729 (100), 792.5537 (81), 86.0693 (3)
0.007421
2.95E − 04
0.325565
0.199228
0.003948
0.016788
1.225
1.209
1.56
1.56

RP
−
827.5697
12.3
PC(16:0/20:4) [M+FA]
LMGP0101
20ev: 826.5624 [M+FA−H]−; 766.5390 (100), 767.5404 (34), 826.5603 (25)
0.007421
3.00E − 04
0.366474
0.180537
0.010339
0.007452
−1.09
−1.09
1.58
1.58

HILIC
+
763.5173
0.706
PE(16:0/22:6)
LMGP0201
20ev: 762.5062 [M−H]−; 327.2371 (100), 762.5071 (72), 344.5837 (55), 255.2369 (38), 764.5264 [M+H]+; 764.5264 (100), 623.4989 (50)
0.024721
0.002518
0.663723
0.142026
0.003572
0.007017
−1.49
−1.44
2.18
1.42

RP
+
805.5637
12.22
PC(16:0/22:6)
LMGP0101
20ev: 806.5696 [M+H]+; 184.0729 (100), 806.5701 (86), 524.3132 (5), 20ev: 850.5625 [M+FA−H]−; 93 (100), 850.5600 (39), 791.5447 (39), 327.2323 (18), 44.9986 (14), 851.5603 (13), 255.2322 (6)
0.027808
0.001026
0.599027
0.187622
0.116353
0.004307
1.099
1.163
1.37
1.71

HILIC
+
809.5935
0.632
PC(18:0/20:4)
LMGP0101
20ev: 854.588 [M+FA−H]−; 794.5697 (100), 854.5934 (29), 795.5785 (21), 303.2318 (11), 283.26 (7); 10ev: 810.6012 [M+H]+; 810.6045 (100), 184.073 (17)
0.09676
0.028541
0.784922
0.504751
0.047911
0.210616
−1.1
−1.06
1.59
0.99

HILIC
+
523.3635
1.08
LysoPC(0:0/18:0)
LMGP0105
40ev: 568.3607 [M+FA−H]−; 508.3418 (100), 283.2632 (30), 224.0657 (10), 44.9993 (9), 20ev: 524.3706 [M+H]+; 184.0731 (100), 524.3689 (10)
0.047961
0.008856
0.84261
0.146463
0.279208
0.014867
1.137
1.354
1.48
1.60

HILIC
+
521.3482
1.091
LysoPC(0:0/18:1)
LMGP0105
20ev: 522.3570 [M+H]+; 184.0738 (100), 522.3606 (7), 86.0964 (3)
0.040335
0.005835
0.84261
0.15821
0.230192
0.017452
1.157
1.355
1.47
1.66

HILIC
+
495.334
1.15
LysoPC(0:0/16:0)
LMGP0105
20ev: 540.3303 [M+FA−H]−; 255.2323 (100), 480.3042 (62), 242.0762 (7); 20ev: 496.3492 [M + H]+; 184.0764 (100), 496.3477 (5), 86.0974 (2)
0.047031
0.007895
0.887095
0.142026
0.436041
0.01111
1.085
1.318
1.47
1.52

HILIC
+
523.3653
1.2
LysoPC(18:0/0:0)
LMGP01050026
0.037508
0.005217
0.84261
0.142026
0.289454
0.008929
1.129
1.365
1.42
1.54

HILIC
+
521.3502
1.215
LysoPC(18:1/0:0)
LMGP0105
40ev: 566.345 [M+FA−H]−; 281.2481 (100), 224.0735 (14), 44.999 (9)
0.024721
0.002743
0.84261
0.146463
0.24612
0.014391
1.159
1.392
1.34
1.54

HILIC
+
519.3347
1.227
LysoPC(18:2/0:0)
LMGP0105
20ev: 520.3400 [M+H]+; 184.0733 (100), 520.339 (10), 104.1067 (2), 86.095 (2), 20ev: 564.3326 [M+FA−H]−; 504.3113 (100), 279.2337 (68), 44.9995 (10)
0.009186
4.33E − 04
0.887095
0.142026
0.45312
0.006973
1.093
1.45
1.29
1.50

HILIC
+
495.3345
1.284
LysoPC(16:0/0:0)
LMGP01050018
0.103087
0.031116
0.909764
0.185013
0.609742
0.02529
1.055
1.276
1.48
1.59

HILIC
−
222.02
1.302
Unknown

10ev: 221.0126; 221.012 (100), 148.9898 (35), 130.9797 (12), 94.9242 (5), 92.9242 (3)
0.016106
9.59E − 04
0.585329
0.33396
0.013859
0.057424
−1.16
−1.13
1.64
1.14

RP
+
231.1466
2.065
Butyryl-L-carnitine
MID964
Verified with standard. 10ev: 232.1545 [M+H]+; 232.1549 (100), 85.0276 (95), 173.0738 (16), 60.0784 (3)
0.49639
0.218745
0.46469
0.782517
0.047231
0.482634
−1.39
−1.14
1.57
0.81

RP
−
173.9987
2.408
Unknown

20ev: 172.9929 [M−H]−?; 93.0348 (100), 34.9705 (6)
0.404675
0.116603
0.861434
0.273229
0.594768
0.030835
−1.19
−2.01
1.36
1.82

RP
+
250.1785
5.666
Unknown

10ev: 251.1860; 59.05 (100), 117.0906 (19), 41.0397 (7), 251.1881 (6)
0.009147
1.34E − 04
0.824211
0.216978
0.380635
0.008023
−1.02
1.27
0.86
1.55

HILIC
+
155.0695
6.772
L-Histidine
MID21
Verified with standard 10ev: 156.0765 [M+H]+; 110.0712 (100), 156.0758 (30), 95.0609 (11), 83.0602 (5)
0.016822
0.001302
0.670033
0.826355
0.008034
0.614635
1.207
1.035
1.76
0.86

RP
−
396.1973
8.611
Pregnenolone sulfate
MID5740
20ev: 395.1903 [M−H]−; 395.1903 (100), 96.9607 (23), 248.8949 (7), 305.7866 (6)
0.082724
0.010827
0.384093
0.398782
0.015051
0.079233
−1.5
−1.38
1.71
1.12

RP
+
567.3333
10.12
LysoPC(22:6/0:0) (RP a)
LMGP0105
20ev: 612.3301 [M+FA−H]−; 552.3112 (100), 327.2352 (25), 283.243 (8), 44.9985 (4), 20ev: 568.3392 [M+H]+; 104.1063 (100), 184.0726 (98), 568.3392 (86), 86.096 (14)
0.008933
1.03E − 04
0.457269
0.105666
0.041413
5.08E − 04
1.291
1.498
1.23
1.51

RP
+
505.4133
11.69
Unknown

10ev: 506.4203; 60.0816 (100), 61.083 (4)
0.159382
0.029154
0.314458
0.799741
0.01188
0.505812
1.307
1.071
1.79
0.87

RP
−
849.5526
11.87
PC(16:1/22:6) [M+FA]
LMGP0101
40ev: 848.5490 [M+FA−H]−; 327.2322 (100), 253.2164 (52), 283.2444 (49), 478.2944 (31), 44.9987 (30), 20ev: 804.5531 [M+H]+; 184.0718 (100), 804.5512 (67)
0.03311
0.002648
0.366474
0.309525
0.008646
0.041723
1.309
1.235
1.56
1.32

RP
−
801.5533
12.02
PC(16:1/18:2) [M+FA]
LMGP0101
20ev: 756.5562 [M+H]+; 184.0737 (100), 756.5533 (78), 20ev: 800.5511 [M+FA−H]−; 740.5234 (100), 800.545 (53), 279.2312 (15), 253.2161 (11), 44.9988 (6)
0.035314
0.002909
0.372759
0.312787
0.01227
0.043172
−1.5
−1.44
1.76
1.19

RP
+
415.3572
12.15
Unknown

10ev: 416.3633; 416.3617 (100), 164.1129 (4), 191.1039 (3)
0.067283
0.006009
0.513659
0.333893
0.070642
0.049803
−1.2
−1.2
1.33
1.60

RP
+
831.5781
12.33
PC(18:1/22:6)
LMGP0101
20ev: 832.5857 [M+H]+; 184.0738 (100), 832.585 (97), 40ev: 876.5774 [M+FA−H]−; 327.2333 (100), 281.2478 (67), 283.243 (32), 816.5537 (18), 44.9985 (18)
0.009147
1.40E − 04
0.534883
0.105666
0.07969
3.23E − 04
1.163
1.333
1.44
1.66

RP
+
771.5775
12.58
Unknown PC

20ev: 772.5868 [M+H]+; 772.587 (100), 184.0736 (94), 86.0975 (5), double peak, no MS/MS from the first one
0.356662
0.11337
0.314458
0.848946
0.011944
0.591385
−1.24
1.046
1.98
0.98

RP
+
759.5784
12.8
PC(16:0/18:1)
LMGP0101
20ev: 804.5795 [M+FA−H]−; 745.5581 (28), 804.5754 (26), 281.2484 (14), 255.2315 (6), 44.9982 (3), 20ev: 760.5859 [M+H]+; 184.0732 (100), 760.5862 (82)
0.015602
4.25E − 04
0.298009
0.596521
0.008491
0.213038
−1.15
−1.06
1.58
1.24

RP
+
833.5943
12.81
PC(18:0/22:6)
LMGP0101
20ev:834.6028 [M+H]+; 834.5985 (100), 184.0727 (91), 86.0979 (2), 40ev: 878.5934 [M+FA−H]; 327.2332 (100), 283.2625 (48), 283.2436 (36), 818.5716 (20), 508.3402 (20), 44.9983 (10)
0.045589
0.002876
0.653825
0.216978
0.165991
0.007318
1.125
1.251
1.23
1.51

RP
+
837.6245
13.57
Unknown PC

20ev: 838.6312 [M+H]+; 838.6324 (100), 184.0728 (56), 185.0744 (6)
0.077038
0.007688
0.285943
0.660969
0.007543
0.291104
−1.42
−1.13
1.57
1.11

RP
+
787.6087
13.61
PC(18:0/18:1)
LMGP0101
20ev: 788.6156 [M+H]+; 184.0728 (100), 788.6134 (100), 84.0787 (2); 40ev:832.6118 [M+FA−H]−; 281.2489 (100), 283.2635 (46), 772.5866 (14), 44.9979 (11)
0.091692
0.010657
0.301815
0.839344
0.008949
0.568381
−1.32
−1.06
1.60
0.84

The characteristics include both uncorrected and FDR- (Benjamini-Hochberg false discovery rate-) corrected values on d0 and d1, fold changes, variable influence on projection (VIP) values, and identification references, together with parameters for the LC-MS analysis, including chromatography (column), ionization mode in mass spectrometry (Ioni), molecular weight (MW), retention time (RT), and fragment ions in tandem mass spectrometry (MS/MS fragments). (20 patients with complicated course of febrile neutropenia (complicated) and 65 patients without complicated course of febrile neutropenia (noncomplicated). ADTS: androsterone sulfate; DHTS: dihydrotestosterone sulfate; PC: phosphatidylcholine; PE: phosphatidylethanolamine; LysoPC: lysophosphatidylcholine. ANOVA (analysis of variance) and groupwise t-test (unpaired t-test) comparing the fold changes between complicated versus noncomplicated patients. values < 0.05. Fold changes = average fold changes when comparing complicated group against noncomplicated group, with values. Fold changes ≥ ± 1.2 were considered significant. Positive values indicate increased plasma levels in complicated versus noncomplicated patients, whereas negative values indicate decreased plasma levels in complicated versus noncomplicated patients. ID identification of metabolites is based on references: MID refers to METLIN database https://metlin.scripps.edu/index.php, LMGP refers to LIPID MAPS database http://lipidmaps.org/data/structure/index.html, and HMDB refers to Human Metabolome Database http://www.hmdb.ca/.

2.7. The Statistical Analysis of the Conventional Biomarkers

CRP and PCT values were expressed as medians and interquartiles due to the skewed distribution. The Mann–Whitney U test was used to detect differences between the groups in continuous variables. The association between categorical variables was studied by χ² test with Spearman’s correlation. A value less than 0.05 was considered significant. Data analyses were conducted with SPSS 21.0 Software (SPSS Inc., Chicago, Illinois, USA).

3. Results

3.1. Course of Febrile Neutropenia and Blood Culture Findings

The characteristics of patients with complicated course of febrile neutropenia are presented in Table 2. Altogether, the group included twenty patients (24%). Twelve patients had bacteremia without any other signs of complication. Eight patients fulfilled the criteria for complicated sepsis, and three of them developed septic shock. Altogether, six patients needed intensive care unit treatment and three of them died (mortality 3% in the whole series).

Table 2: Characteristics of the patients with complicated febrile neutropenia.

Gender/age (years)
Blood culture
Time
Characteristics of severe sepsis
qSOFA
ICU
Outcome

M/65
Pseudomonas aeruginosa
1
Pneumonia, respiratory failure
3
Yes
Died

F/70
Negative, RSV-pneumonia
6
Pneumonia, septic shock
2
Yes
Died

M/47
Enterococcus faecium
4
Pneumonia, respiratory failure
1
Yes
Died

M/67
Negative
7
Pneumonia, respiratory failure
2
Yes
Recovered

M/58
Staphylococcus epidermidis
2
Septic shock, respiratory failure
2
Yes
Recovered

M/63
Enterococcus faecium
5
Respiratory failure, ileus
2
Yes
Recovered

M/68
Streptococcus mitis
11
Respiratory failure (need for NIV)
2
No (CCU)
Recovered

M/51
Enterococcus faecium
1
Low blood pressure, need for plasma expander
1
No
Recovered

F/65
Streptococcus mitis
—
—
—
No
Recovered

M/41
Streptococcus salivarius
—
—
—

M/31
Gemella morbillorum
—
—
—
No
Recovered

M/62
Staphylococcus haemolyticus
—
—
—
No
Recovered

M/36
Staphylococcus haemolyticus
—
—
—
No
Recovered

F/62
Staphylococcus epidermidis
—
—
—
No
Recovered

M/65
Staphylococcus epidermidis
—
—
—
No
Recovered

M/69
Enterococcus faecium
—
—
—
No
Recovered

M/58
Enterococcus faecium
—
—
—
No
Recovered

F/49
Klebsiella oxytoca
—
—
—
No
Recovered

F/56
Escherichia coli
—
—
—
No
Recovered

M/61
Candida krusei
—
—
—
No
Recovered

M: male; F: female; Time: days from the onset of fever to severe sepsis; qSOFA: quick sepsis-related organ failure assessment; ICU: intensive care unit; NIV: noninvasive ventilation; CCU: cardiac care unit.

The blood cultures were positive in 18 out of 85 patients (21%) with fourteen gram-positive and three gram-negative bacteremias, respectively. The gram-positive findings included Enterococcus faecium (), Staphylococcus epidermidis (), Streptococcus mitis (), Staphylococcus haemolyticus (), Streptococcus salivarius (), and Gemella morbillorum (). The gram-negative findings were Klebsiella oxytoca (), Escherichia coli (), and Pseudomonas aeruginosa (). One case of fungemia was found (Candida krusei).

3.2. Conventional Biomarkers

The medians (IQ) of CRP for patients with complicated course of febrile neutropenia were 51 (26–100) on d0 and 100 (57–214) on d1. In patients without complications, CRP values were 36 (19–67) on d0 and 67 (35–95) on d1. There was a significant statistical difference on d1 (). The medians (IQ) of PCT for patients with complicated course of febrile neutropenia were 0.164 (0.120–0.279) on d0 and 0.350 (0.190–1.855) on d1. In patients without complications, PCT values were 0.122 (0.077–0.193) on d0 and 0.163 (0.104–0.320) on d1. There was a significant statistical difference between the groups on both days. values were 0.027 and 0.020 on day d0 and d1, respectively.

3.3. Nontargeted Metabolite Profiling

Principal component analysis of the molecular features collected at four analytical modes of LC-MS analysis is presented in the Supplementary Figure (available here). Altogether, 90 molecular features fulfilled the filtering criteria and were considered to be important for explaining the metabolic differences between patients with noncomplicated versus complicated course of febrile neutropenia. Of these 90 molecular features, 52 were tentatively identified corresponding to 25 different metabolites (Table 1) (Figure 1).

(a)

(b)

(c)

(d)

Figure 1: Panels (a) and (c) show the first two components of principal component analysis (PCA) for all the molecular features (log-transformed) from day 0 and day 1, respectively. Panels (b) and (d) show the first two components of a partial least-squares discriminant analysis (PLS-DA) for all the molecular features (log-transformed) from day 0 and day 1, respectively. In both PLS-DA models ((b) and (d)), six components cumulatively explain 97% of the variance in the data between the groups on day 0 (b) and on day 1 (d). These six components explain 49% and 58% of the variance in the data on day 0 and day 1, respectively. White dots are patients without complicated course of febrile neutropenia. Black dots are patients with complicated course of febrile neutropenia.

Six molecular features fulfilled the strictest filtering criteria and were considered as the most significant markers to differentiate between noncomplicated and complicated course of febrile neutropenia (Figure 2). Of these markers, one metabolite was unambiguously identified to be citrulline and two others tentatively identified to be ADTS/DHTS and PE(18:0/20:4), whereas three remained unidentified (Table 1). While ADTS/DHTS and PE(18:0/20:4) were significantly increased, citrulline demonstrated low levels in patients with complicated course of febrile neutropenia (Figure 2).

(a)

(b)

(c)

(d)

(e)

(f)

Figure 2: The most differential metabolites and their changes on day 0 (d0) and day 1 (d1) from beginning of the febrile neutropenia. (a) Androsterone sulfate (ADTS)/dihydrotestosterone sulfate (DHTS), (b) citrulline, (c) phosphatidylethanolamine PE(18:0/20:4), (d) phosphatidylcholine PC(15:0/22:6), (e) phosphatidylcholine PC(16:0/20:4), and (f) phosphatidylcholine PC(16:0/22:6). “No” indicates patients without complicated course of febrile neutropenia. “Yes” indicates patients with complicated course of febrile neutropenia. Day 0 indicates the onset of febrile neutropenia. Day 1 indicates the next day. Horizontal line represents median, box encompasses 25th and 75th percentiles, and error bar encompasses 10th and 90th percentiles.

Twelve molecular features fulfilled the second most strict criteria and were considered important in explaining the differences. These features were tentatively identified to correspond to four different lipid metabolites, namely, PC(15:0/22:6), PC(16:0/20:4), PC(16:0/22:6), and PE(16:0/22:6) (Table 1). While PC with 20 : 4 fatty-acyl side chain increased, the PCs with 22 : 6 side chain decreased in patients with complicated course of febrile neutropenia (Figure 2). PE(16:0/22:6) showed a similar increase in patients with complicated course of febrile neutropenia as PE(18:0/20:4), but it did not meet the criteria set for the most significant markers due to a large relative standard deviation (Table 1).

Additionally, 19 metabolites were identified which fulfilled the lowest criteria. These tentatively identified metabolites containing various lysophosphatidylcholines and phosphatidylcholines, butyryl-L-carnitine, L-histidine, and putative pregnenolone sulfate (Table 1). All of the LysoPCs showed decreased levels in patients with complicated course of febrile neutropenia. Whereas in the case of PCs, a similar behavior was seen as with more important metabolite markers as the PCs containing a 22 : 6 fatty-acyl side chain, as they showed decreased levels in patients with complicated course of febrile neutropenia, while other PCs were increased in the complicated group. Butyryl-L-carnitine and tentative pregnenolone sulfate showed increased levels and L-histidine decreased levels in patients with complicated course of febrile neutropenia (Table 1).

3.4. Correlation of the Metabolomics Findings with Conventional Biomarkers

ADTS/DHTS had significant correlations with CRP and PCT on both days. With CRP, the values were 0.002 and 0.003 on d0 and d1, respectively, and with PCT they were <0.001 on both days. Citrulline had a significant correlation on d1 with CRP () and PCT (). PE(18:0/22:4) had a significant correlation with PCT on both time points ( and <0.001 on d0 and d1, resp.), but with CRP only on d1 (). Naturally, ADTS/DHTS had a significant difference according to gender (Table 3). However, when looking at the differences only among men, the significant difference remained between groups of complicated versus noncomplicated course of febrile neutropenia ( values < 0.001 and 0.001 on days d0 and d1, resp.). Also, women with complicated course of febrile neutropenia had higher ADTS/DHTS values than women without complications on d0 had () (Table 4). The effects of andropause and menopause were evaluated by dividing patients to those under fifty years of age and those over fifty years (Table 5). Citrulline or PE(18:0/22:4) levels did not have any difference due to gender or age (data not shown).

Table 3: Influence of gender to androsterone sulfate (ADTS)/dihydrotestosterone sulfate (DHTS).

Women
Men
values

ADTS/DHTS d0
3.0 × 10⁴ (2.0–5.5 × 10⁴)
6.7 × 10⁴ (4.2–10.7 × 10⁴)
<0.001

ADTS/DHTS d1
3.0 × 10⁴ (2.5–4.2 × 10⁴)
7.0 × 10⁴ (4.1–12 × 10⁴)
<0.001

Medians (IQ). Data are expressed as medians (interquartile ranges).

Table 4: Influence of complicated course of fever to androsterone sulfate (ADTS)/dihydrotestosterone sulfate (DHTS).

Men with complication
Men without complication
values

ADTS/DHTS d0
12.6 × 10⁴ (7.9–14 × 10⁴)
4.9 × 10⁴ (3.6-7.8 × 10⁴)
<0.001

ADTS/DHTS d1
12.6 × 10⁴ (8.3–14 × 10⁴)
5.7 × 10⁴ (3.6–8.8 × 10⁴)
0.001

Women with complication
Women without complication

ADTS/DHTS d0
6.8 × 10⁴ (4.2–6.8 × 10⁴)
2.5 × 10⁴ (2.0–4.3 × 10⁴)
0.045

ADTS/DHTS d1
7.0 × 10⁴ (1.9–13 × 10⁴)
2.9 × 10⁴ (2.5–4.0 × 10⁴)
0.322

Data are expressed as medians (interquartile ranges).

Table 5: Influence of age to androsterone sulfate (ADTS)/dihydrotestosterone sulfate (DHTS).

Men under 50 years
Men over 50 years
values

ADTS/DHTS d0
8.5 × 10⁴ (5.6–13 × 10⁴)
5.9 × 10⁴ (3.7–9.3 × 10⁴)
0.085

ADTS/DHTS d1
10 × 10⁴ (6.2–12 × 10⁴)
6.4 × 10⁴ (3.6–10.2 × 10⁴)
0.121

Women under 50 years
Women over 50 years

ADTS/DHTS d0
5.9 × 10⁴ (3.5–13.7 × 10⁴)
2.5 × 10⁴ (2.0–4.0 × 10⁴)
0.068

ADTS/DHTS d1
8.5 × 10⁴ (2.8–14.5 × 10⁴)
2.9 × 10⁴ (2.4–4.1 × 10⁴)
0.144

Data are expressed as medians (interquartile ranges).

4. Discussion

In this study, we describe how plasma samples from patients with hematological malignancy were analyzed with nontargeted metabolite profiling to find out potential novel small molecular biomarker candidates for complicated course of febrile neutropenia. This study is until now the largest study evaluating metabolite profiles in hematological patients with febrile neutropenia after intensive chemotherapy. ADTS/DHTS, citrulline, and PE(18:0/22:4) were all found to be important metabolic features differentiating patients with complicated versus noncomplicated course of febrile neutropenia at the start of fever (d0 and d1). These metabolites need further studies and validation in this patient cohort, where early and predictive biomarkers are urgently needed.

The most significant molecular feature was tentatively identified to be either ADTS or DHTS, which increased significantly in patients with complicated course of febrile neutropenia. ADT and its sulfonated form androsterone sulfate (ADTS) are intermediates in the metabolic route of testosterone and its metabolite, dihydrotestosterone (DHT) [29]. The conversion of dehydroepiandrosterone DHEA to potent androgens and/or estrogens takes place in peripheral tissues by the 17β-HSD enzyme family [30]. In the testis and in the muscle, the conversion occurs via testosterone to DHT through the conventional front door pathway where DHEA is transformed to 4-androstenedione (4-dione) and further converted into testosterone. In tissues that express 5α-reductase, such as the prostate, liver, and skin, DHT is produced mostly from 4-dione without testosterone by the backdoor pathway (Figure 3) [31, 32]. Also, monocyte-derived macrophages have the ability to convert DHEAS to androgens [33].

Figure 3: Schematic representation for the metabolite route of testosterone and dihydrotestosterone. Front door pathway (∆4 and ∆5 via testosterone) and the backdoor pathway (5α and 5α, 3α routes not via testosterone). Modified from Fukami et al. [32] and Luu-The [31].

To the best of our knowledge, there are no previous studies of ADT or ADTS in sepsis patients. Androstanedione is the intermediate between ADT and DHT in the backdoor pathway, and it has been studied as a part of adrenal and testicular steroidogenesis in patients with burn trauma and in intensive care unit patients, but findings have been contradictory [34–36]. In animal and in vitro studies, immunomodulatory effects of testosterone on macrophage function are mediated via the 5α-reductase-dependent conversion of testosterone to DHT [37, 38].

In patients with sepsis, testosterone levels were below the normal range for men and estradiol levels were increased in both postmenopausal women and men [39]. Male sex steroids appear to be immunosuppressive, whereas female sex steroids increase the activity of humoral immune responses. The major source of enhanced estradiol production has been suggested to be the aromatization of testosterone to estradiol [39–41]. The potential role of ADTS/DHTS in sepsis remains to be clarified in future studies.

Another important finding in the present study was that circulating citrulline levels decreased significantly in patients with complicated course of febrile neutropenia. L-Citrulline is used as a biomarker of enterocyte functional mass [42]. Studies in patients with septic shock or those with multiple organ failure have shown that patients with the lowest citrulline nadirs had the highest risk of gastrointestinal tract translocation of bacteria, and plasma citrulline concentration at was an independent prognostic factor for mortality [43, 44]. Herbers et al. [45] studied citrulline concentrations in neutropenic patients after high-dose melphalan and ASCT. The patients with bacteremia had significantly lower citrulline concentrations on the first day of fever than did those without bacteremia. Our results are in line with these findings: citrulline levels were significantly lower in patients with complicated course of febrile neutropenia. Of note, citrulline serves also as a substrate for nitric oxide production and low levels are associated with acute respiratory distress syndrome (ARDS) [46].

In our study, phospholipid PE(18:0/20:4) increased significantly in patients with complicated course of febrile neutropenia. Also, 12 other molecular features, which fulfilled the second most strict criteria, were all phospholipids. PC is the most abundant and PE the second most abundant phospholipid in mammalian cells [47]. PE and phosphatidylserine (PS) are aminophospholipids, and because of their negative charge, they are localized in the inner leaflet on the plasma membrane. Neutrally charged PC and sphingomyelin localize on the outer leaflet [48, 49]. As a response to stimulus, cytoskeleton is reorganized and the asymmetric distribution of the phospholipid membrane is modified with exposure of PS and PE at the cell surface. Cellular blebbing then occurs, ultimately leading to the release of microparticles (MPs) [50, 51]. MPs originate from several cell types such as leukocytes, erythrocytes, platelets, or endothelial cells. They have been implicated in the pathophysiology of a range of diseases, including infectious diseases [49]. Also, a major role has been suggested for MPs in propagating proinflammatory and proanticoagulant activity in sepsis [50].

PE is abundantly present on the outer surface of mammalian MPs of a variety of human cells, and PE provides a sensitive molecular marker for detecting and characterizing MPs [50]. Clark et al. [48] studied more specifically activated human platelets and found that they externalize two PSs and five PEs: PS(18:0a/20:4), PS(18:0a/18:1), PE(16:0p/20:4), PE(18:0a/18:0), PE(18:1p/20:4), PE(18:0p/20:4), and PE(18:0a/20:4). These are also present in MPs from activated platelets. PS is required for coagulation to take place, and PE can further enhance these reactions. As the acyl-chain length of PE is an essential determinant of this enhancing effect to thrombin generation with PS, the most effective PE was found to be PE(18:0/20:4).

One of the major findings in our study was the positive association of PE(18:0/20:4) with the development of complicated course of neutropenic sepsis. Previously, PE(18:0/20:4) was found to be one of the relevant phospholipids to detect MPs, but the context with clinical data of sepsis has been unsettled. The studies of Larson et al. [51] and Clark et al. [48] were performed in vitro, and there are no studies on PE as predictors on sepsis or neutropenic sepsis in humans.

Changes in the plasma phospholipid profile have been found in sepsis. Previously, Rival et al. [52] observed lower total concentrations of phospholipids in septic patients compared to the reference values. This decrease concerned mainly n-6 and n-3 polyunsaturated fatty acids, especially those with a carbon , which was associated with mortality. Also in patients with ARDS, a significant decrease in the levels of various fatty acids especially in the polysaturated fatty acids (PUFA) including docosahexaenoic acid (DHA, 22:6n3) has been found. DHA levels were decreased also in patients who were at risk of ARDS [53]. In the experimental study, DHA and eicosapentaenoic acid (EPA, 20:5n3) decreased lipopolysaccharide-induced nuclear factor-kappaB activation and monocyte chemoattractant protein-1 production and increased peroxisome proliferator-activated receptor mRNA expression [54]. The decrease in PUFAs has been assumed to be due to their degeneration by reactive oxygen species or a higher synthesis of inflammatory lipid mediators because these PUFAs are the precursors of eicosanoids and docosanoids, which are involved in inflammation, vasomotricity, and capillary permeability [52].

5. Conclusion

Although nontargeted metabolite profiling is a hypothesis-generating method, these results already provide a new insight into metabolite markers to differentiate between noncomplicated and complicated courses of febrile neutropenia. Our study represents the largest study in this patient population, even though the relatively low number of patients reduces the statistical power of the results.

Our findings are consistent with previously published data about citrulline and phospholipids. The current findings are clinically relevant when compared with CRP and PCT, which are widely used biomarkers in febrile neutropenia. Most of the previous studies are performed in animals or in vitro in contrast to our prospective clinical study, which highlights the importance of current findings.

The pathways of androgen metabolism merit further studies in patients with febrile neutropenia. The results are preliminary and suggest metabolic changes in patients with complicated febrile neutropenia. Further targeted tests to identify and confirm these changes are needed. The whole pathway of metabolites of steroid synthesis in the adrenal gland should be quantified. If specific, consistent metabolite changes can be identified, a routine laboratory-compatible method may be developed.

Ethical Approval

This study is conducted according to the principles expressed in the Declaration of Helsinki and was approved by the local ethics committee (100/2006).

Consent

All patients gave a written informed consent at the enrollment.

Conflicts of Interest

The authors declare that they have no conflict of interest in relation to the work described.

Acknowledgments

The authors acknowledge Ms. Miia Reponen for skillful technical assistance in LC-MS analyses, and the technical assistance of Ms. Raija Isomäki and Anu Holopainen, MSc, is gratefully acknowledged. The study was financially supported by a study grant from the North Savo Hospital District (VTR). Financial support from Biocenter Finland and Academy of Finland is acknowledged.

Supplementary Materials

Supplementary Figure: principal component analysis of the molecular features collected at four analytical modes of liquid chromatography-mass spectrometry analysis. (Supplementary Materials)

References

O. Penack, C. Becker, D. Buchheidt et al., “Management of sepsis in neutropenic patients: 2014 updated guidelines from the Infectious Diseases Working Party of the German Society of Hematology and Medical Oncology (AGIHO),” Annals of Hematology, vol. 93, no. 7, pp. 1083–1095, 2014.

K. A. Thursky and L. J. Worth, “Can mortality of cancer patients with fever and neutropenia be improved?” Current Opinion in Infectious Diseases, vol. 28, no. 6, pp. 505–513, 2015.

S. Hämäläinen, T. Kuittinen, I. Matinlauri, T. Nousiainen, I. Koivula, and E. Jantunen, “Neutropenic fever and severe sepsis in adult acute myeloid leukemia (AML) patients receiving intensive chemotherapy: causes and consequences,” Leukemia & Lymphoma, vol. 49, no. 3, pp. 495–501, 2008.

S. Hämäläinen, A. Juutilainen, I. Matinlauri et al., “Serum vascular endothelial growth factor in adult haematological patients with neutropenic fever: a comparison with C-reactive protein,” European Journal of Haematology, vol. 83, no. 3, pp. 251–257, 2009.

C. Pierrakos and J. L. Vincent, “Sepsis biomarkers: a review,” Critical Care, vol. 14, no. 1, article R15, 2010.

I. Koivula, S. Hämäläinen, E. Jantunen et al., “Elevated procalcitonin predicts gram-negative sepsis in haematological patients with febrile neutropenia,” Scandinavian Journal of Infectious Diseases, vol. 43, no. 6-7, pp. 471–478, 2011.

C. H. Johnson, J. Ivanisevic, and G. Siuzdak, “Metabolomics: beyond biomarkers and towards mechanisms,” Nature Reviews Molecular Cell Biology, vol. 17, no. 7, pp. 451–459, 2016.

V. D. de Mello, J. Paananen, J. Lindström et al., “Indolepropionic acid and novel lipid metabolites are associated with a lower risk of type 2 diabetes in the Finnish Diabetes Prevention Study,” Scientific Reports, vol. 7, article 46337, 2017.

M. Kiehntopf, D. Schmerler, F. M. Brunkhorst et al., “Mass spectometry-based protein patterns in the diagnosis of sepsis/systemic inflammatory response syndrome,” Shock, vol. 36, no. 6, pp. 560–569, 2011.

D. Schmerler, S. Neugebauer, K. Ludewig, S. Bremer-Streck, F. M. Brunkhorst, and M. Kiehntopf, “Targeted metabolomics for discrimination of systemic inflammatory disorders in critically ill patients,” Journal of Lipid Research, vol. 53, no. 7, pp. 1369–1375, 2012.

L. Su, Y. Huang, Y. Zhu et al., “Discrimination of sepsis stage metabolic profiles with an LC/MS-MS-based metabolomics approach,” BMJ Open Respiratory Research, vol. 1, no. 1, article e000056, 2014.

S. Neugebauer, E. J. Giamarellos-Bourboulis, A. Pelekanou et al., “Metabolite profiles in sepsis: developing prognostic tools based on the type of infection,” Critical Care Medicine, vol. 44, no. 9, pp. 1649–1662, 2016.

M. Garcia-Simon, J. M. Morales, V. Modesto-Alapont et al., “Prognosis biomarkers of severe sepsis and septic shock by ¹H NMR urine metabolomics in the intensive care unit,” PLoS One, vol. 10, no. 11, article e0140993, 2015.

M. Ferrario, A. Cambiaghi, L. Brunelli et al., “Mortality prediction in patients with severe septic shock: a pilot study using a target metabolomics approach,” Scientific Reports, vol. 6, no. 1, article 20391, 2016.

Z. Liu, P. Yin, R. Amathieu, P. Savarin, and G. Xu, “Application of LC-MS-based metabolomics method in differentiating septic survivors from non-survivors,” Analytical and Bioanalytical Chemistry, vol. 408, no. 27, pp. 7641–7649, 2016.

M. E. Richter, S. Neugebauer, F. Engelmann et al., “Biomarker candidates for the detection of an infectious etiology of febrile neutropenia,” Infection, vol. 44, no. 2, pp. 175–186, 2016.

W. T. Hughes, D. Armstrong, G. P. Bodey et al., “2002 guidelines for the use of antimicrobial agents in neutropenic patients with cancer,” Clinical Infectious Diseases, vol. 34, no. 6, pp. 730–751, 2002.

R. C. Bone, R. A. Balk, F. B. Cerra et al., “Definitions for sepsis and organ failure and guidelines for the use of innovative therapies in sepsis,” Chest, vol. 101, no. 6, pp. 1644–1655, 1992.

M. Singer, C. S. Deutschman, C. W. Seymour et al., “The third international consensus definitions for sepsis and septic shock (Sepsis-3),” JAMA, vol. 315, no. 8, pp. 801–810, 2016.

K. Hanhineva, M. A. Lankinen, A. Pedret et al., “Nontargeted metabolite profiling discriminates diet-specific biomarkers for consumption of whole grains, fatty fish, and bilberries in a randomized controlled trial,” The Journal of Nutrition, vol. 145, no. 1, pp. 7–17, 2015.

Y. Benjamini and Y. Hochberg, “Controlling the false discovery rate: a practical and powerful approach to multiple testing,” Journal of the Royal Statistical Society Series B Methodological, vol. 57, no. 1, pp. 289–300, 1995.

M. Sugimoto, M. Kawakami, M. Robert, T. Soga, and M. Tomita, “Bioinformatics tools for mass spectroscopy-based metabolomic data processing and analysis,” Current Bioinformatics, vol. 7, no. 1, pp. 96–108, 2012.

R. G. Brereton and G. R. Lloyd, “Partial least squares discriminant analysis: taking the magic away,” Journal of Chemometrics, vol. 28, no. 4, pp. 213–225, 2014.

C. A. Smith, G. O’Maille, E. J. Want et al., “METLIN: a metabolite mass spectral database,” Therapeutic Drug Monitoring, vol. 27, no. 6, pp. 747–751, 2005.

D. S. Wishart, T. Jewison, A. C. Guo et al., “HMDB 3.0–the human metabolome database in 2013,” Nucleic Acids Research, vol. 41, no. D1, pp. D801–D807, 2013.

E. Fahy, M. Sud, D. Cotter, and S. Subramaniam, “LIPID MAPS online tools for lipid research,” Nucleic Acids Research, vol. 35, Supplement 2, pp. W606–W612, 2007.

Y. Q. Xia and M. Jemal, “Phospholipids in liquid chromatography/mass spectrometry bioanalysis: comparison of three tandem mass spectrometric techniques for monitoring plasma phospholipids, the effect of mobile phase composition on phospholipids elution and the association of phospholipids with matrix effects,” Rapid Communications in Mass Spectrometry, vol. 23, no. 14, pp. 2125–2138, 2009.

R. C. Murphy and P. H. Axelsen, “Mass spectrometric analysis of long-chain lipids,” Mass Spectrometry Reviews, vol. 30, no. 4, pp. 579–599, 2011.

K. Homma, T. Hasegawa, T. Nagai et al., “Urine steroid hormone profile analysis in cytochrome P450 oxidoreductase deficiency: implication for the backdoor pathway to dihydrotestosterone,” The Journal of Clinical Endocrinology & Metabolism, vol. 91, no. 7, pp. 2643–2649, 2006.

F. Labrie, A. Bélanger, J. Simard, and L. C. Van Luu-The, “DHEA and peripheral androgen and estrogen formation: intracrinology,” Annals of the New York Academy of Sciences, vol. 774, no. 1, pp. 16–28, 1995.

V. Luu-The, “Assessment of steroidogenesis and steroidogenic enzyme functions,” The Journal of Steroid Biochemistry and Molecular Biology, vol. 137, pp. 176–182, 2013.

M. Fukami, K. Homma, T. Hasegawa, and T. Ogata, “Backdoor pathway for dihydrotestosterone biosynthesis: implications for normal and abnormal human sex development,” Developmental Dynamics, vol. 242, no. 4, pp. 320–329, 2013.

M. Schmidt, M. Kreutz, G. Löffler, J. Schölmerich, and R. H. Straub, “Conversion of dehydroepiandrosterone to downstream steroid hormones in macrophages,” Journal of Endocrinology, vol. 164, no. 2, pp. 161–169, 2000.

E. D. Lephart, C. R. Baxter, and C. R. Parker Jr, “Effect of burn trauma on adrenal and testicular steroid hormone production,” The Journal of Clinical Endocrinology & Metabolism, vol. 64, no. 4, pp. 842–848, 1987.

C. E. Wade, J. S. Lindberg, J. L. Cockrell et al., “Upon-admission adrenal steroidogenesis is adapted to the degree of illness in intensive care unit patients,” The Journal of Clinical Endocrinology & Metabolism, vol. 67, no. 2, pp. 223–227, 1988.

P. Luppa, R. Munker, D. Nagel, M. Weber, and D. Engelhardt, “Serum androgens in intensive-care patients: correlations with clinical findings,” Clinical Endocrinology, vol. 34, no. 4, pp. 305–310, 1991.

C. P. Schneider, M. G. Schwacha, T. S. A. Samy, K. I. Bland, and I. H. Chaudry, “Androgen-mediated modulation of macrophage function after trauma-hemorrhage: central role of 5α-dihydrotestosterone,” Journal of Applied Physiology, vol. 95, no. 1, pp. 104–112, 2003.

M. Frink, Y. C. Hsieh, S. Hu et al., “Mechanism of salutary effects of finasteride on post-traumatic immune/inflammatory response: upregulation of estradiol synthesis,” Annals of Surgery, vol. 246, no. 5, pp. 836–843, 2007.

M. W. A. Angstwurm, R. Gaertner, and J. Schopohl, “Outcome in elderly patients with severe infection is influenced by sex hormones but not gender,” Critical Care Medicine, vol. 33, no. 12, pp. 2786–2793, 2005.

J. Schröder, V. Kahlke, K. H. Staubach, P. Zabel, and F. Stüber, “Gender differences in human sepsis,” Archives of Surgery, vol. 133, no. 11, pp. 1200–1205, 1998.

M. K. Angele, S. Pratschke, W. J. Hubbard, and I. H. Chaudry, “Gender differences in sepsis: cardiovascular and immunological aspects,” Virulence, vol. 5, no. 1, pp. 12–19, 2014.

G. Piton and G. Capellier, “Biomarkers of gut barrier failure in the ICU,” Current Opinion in Critical Care, vol. 22, no. 2, pp. 152–160, 2016.

G. Piton, C. Manzon, E. Monnet et al., “Plasma citrulline kinetics and prognostic value in critically ill patients,” Intensive Care Medicine, vol. 36, no. 4, pp. 702–706, 2010.

P. Crenn, N. Neveux, S. Chevret et al., “Plasma L-citrulline concentrations and its relationship with inflammation at the onset of septic shock: a pilot study,” Journal of Critical Care, vol. 29, no. 2, pp. 315.e1–315.e6, 2014.

A. H. E. Herbers, N. M. A. Blijlevens, J. P. Donnelly, and T. J. M. de Witte, “Bacteraemia coincides with low citrulline concentrations after high-dose melphalan in autologous HSCT recipients,” Bone Marrow Transplantation, vol. 42, no. 5, pp. 345–349, 2008.

L. B. Ware, J. A. Magarik, N. Wickersham et al., “Low plasma citrulline levels are associated with acute respiratory distress syndrome in patients with severe sepsis,” Critical Care, vol. 17, no. 1, article R10, 2013.

J. E. Vance and G. Tasseva, “Formation and function of phosphatidylserine and phosphatidylethanolamine in mammalian cells,” Biochimica et Biophysica Acta (BBA) – Molecular and Cell Biology of Lipids, vol. 1831, no. 3, pp. 543–554, 2013.

S. R. Clark, C. P. Thomas, V. J. Hammond et al., “Characterization of platelet aminophospholipid externalization reveals fatty acids as molecular determinants that regulate coagulation,” Proceedings of the National Academy of Sciences of the United States of America, vol. 110, no. 15, pp. 5875–5880, 2013.

F. Cognasse, H. Hamzeh-Cognasse, S. Laradi et al., “The role of microparticles in inflammation and transfusion: a concise review,” Transfusion and Apheresis Science, vol. 53, no. 2, pp. 159–167, 2015.

F. Meziani, X. Delabranche, P. Asfar, and F. Toti, “Bench-to-bedside review: circulating microparticles - a new player in sepsis?” Critical Care, vol. 14, no. 5, p. 236, 2010.

M. C. Larson, J. E. Woodliff, C. A. Hillery, T. J. Kearl, and M. Zhao, “Phosphatidylethanolamine is externalized at the surface of microparticles,” Biochimica et Biophysica Acta (BBA) – Molecular and Cell Biology of Lipids, vol. 1821, no. 12, pp. 1501–1507, 2012.

T. Rival, C. Cinq-Frais, S. Silva-Sifontes et al., “Alteration of plasma phospholipid fatty acid profile in patients with septic shock,” Biochimie, vol. 95, no. 11, pp. 2177–2181, 2013.

K. V. Kumar, S. M. Rao, R. Gayani, I. K. Mohan, and M. U. R. Naidu, “Oxidant stress and essential fatty acids in patients with risk and established ARDS,” Clinica Chimica Acta, vol. 298, no. 1-2, pp. 111–120, 2000.

H. Li, X. Z. Ruan, S. H. Powis et al., “EPA and DHA reduce LPS-induced inflammation responses in HK-2 cells: evidence for a PPAR-γ–dependent mechanism,” Kidney International, vol. 67, no. 3, pp. 867–874, 2005.