We present a multivariate metabotyping approach to assess the functional recovery of nonhospitalized COVID-19 patients and the possible biochemical sequelae of “Post-Acute COVID-19 Syndrome”, colloquially known as long-COVID. Blood samples were taken from patients ca. 3 months after acute COVID-19 infection with further assessment of symptoms at 6 months. Some 57% of the patients had one or more persistent symptoms including respiratory-related symptoms like cough, dyspnea, and rhinorrhea or other nonrespiratory symptoms including chronic fatigue, anosmia, myalgia, or joint pain. Plasma samples were quantitatively analyzed for lipoproteins, glycoproteins, amino acids, biogenic amines, and tryptophan pathway intermediates using Nuclear Magnetic Resonance (NMR) spectroscopy and mass spectrometry. Metabolic data for the follow-up patients (n = 27) were compared with controls (n = 41) and hospitalized severe acute respiratory syndrome SARS-CoV-2 positive patients (n = 18, with multiple time-points). Univariate and multivariate statistics revealed variable patterns of functional recovery with many patients exhibiting residual COVID-19 biomarker signatures. Several parameters were persistently perturbed, e.g., elevated taurine (p = 3.6 × 10^–3 versus controls) and reduced glutamine/glutamate ratio (p = 6.95 × 10^–8 versus controls), indicative of possible liver and muscle damage and a high energy demand linked to more generalized tissue repair or immune function. Some parameters showed near-complete normalization, e.g., the plasma apolipoprotein B100/A1 ratio was similar to that of healthy controls but significantly lower (p = 4.2 × 10^–3) than post-acute COVID-19 patients, reflecting partial reversion of the metabolic phenotype (phenoreversion) toward the healthy metabolic state. Plasma neopterin was normalized in all follow-up patients, indicative of a reduction in the adaptive immune activity that has been previously detected in active SARS-CoV-2 infection. Other systemic inflammatory biomarkers such as GlycA and the kynurenine/tryptophan ratio remained elevated in some, but not all, patients. Correlation analysis, principal component analysis (PCA), and orthogonal-partial least-squares discriminant analysis (O-PLS-DA) showed that the follow-up patients were, as a group, metabolically distinct from controls and partially comapped with the acute-phase patients. Significant systematic metabolic differences between asymptomatic and symptomatic follow-up patients were also observed for multiple metabolites. The overall metabolic variance of the symptomatic patients was significantly greater than that of nonsymptomatic patients for multiple parameters (χ²p = 0.014). Thus, asymptomatic follow-up patients including those with post-acute COVID-19 Syndrome displayed a spectrum of multiple persistent biochemical pathophysiology, suggesting that the metabolic phenotyping approach may be deployed for multisystem functional assessment of individual post-acute COVID-19 patients.

We have applied nuclear magnetic resonance spectroscopy based plasma phenotyping to reveal diagnostic molecular signatures of SARS-CoV-2 infection via combined diffusional and relaxation editing (DIRE). We compared plasma from healthy age-matched controls (n = 26) with SARS-CoV-2 negative non-hospitalized respiratory patients and hospitalized respiratory patients (n = 23 and 11 respectively) with SARS-CoV-2 rRT-PCR positive respiratory patients (n = 17, with longitudinal sampling time-points). DIRE data were modelled using principal component analysis and orthogonal projections to latent structures discriminant analysis (O-PLS-DA), with statistical cross-validation indices indicating excellent model generalization for the classification of SARS-CoV-2 positivity for all comparator groups (area under the receiver operator characteristic curve = 1). DIRE spectra show biomarker signal combinations conferred by differential concentrations of metabolites with selected molecular mobility properties. These comprise the following: (a) composite N-acetyl signals from α-1-acid glycoprotein and other glycoproteins (designated GlycA and GlycB) that were elevated in SARS-CoV-2 positive patients [p = 2.52 × 10^–10 (GlycA) and 1.25 × 10^–9 (GlycB) vs controls], (b) two diagnostic supramolecular phospholipid composite signals that were identified (SPC-A and SPC-B) from the –⁺N–(CH₃)₃ choline headgroups of lysophosphatidylcholines carried on plasma glycoproteins and from phospholipids in high-density lipoprotein subfractions (SPC-A) together with a phospholipid component of low-density lipoprotein (SPC–B). The integrals of the summed SPC signals (SPC_total) were reduced in SARS-CoV-2 positive patients relative to both controls (p = 1.40 × 10^–7) and SARS-CoV-2 negative patients (p = 4.52 × 10^–8) but were not significantly different between controls and SARS-CoV-2 negative patients. The identity of the SPC signal components was determined using one and two dimensional diffusional, relaxation, and statistical spectroscopic experiments. The SPC_total/GlycA ratios were also significantly different for control versus SARS-CoV-2 positive patients (p = 1.23 × 10^–10) and for SARS-CoV-2 negatives versus positives (p = 1.60 × 10^–9). Thus, plasma SPC_total and SPC_total/GlycA are proposed as sensitive molecular markers for SARS-CoV-2 positivity that could effectively augment current COVID-19 diagnostics and may have value in functional assessment of the disease recovery process in patients with long-term symptoms.

Quantitative plasma lipoprotein and metabolite profiles were measured on an autonomous community of the Basque Country (Spain) cohort consisting of hospitalized COVID-19 patients (n = 72) and a matched control group (n = 75) and a Western Australian (WA) cohort consisting of (n = 17) SARS-CoV-2 positives and (n = 20) healthy controls using 600 MHz ¹H nuclear magnetic resonance (NMR) spectroscopy. Spanish samples were measured in two laboratories using one-dimensional (1D) solvent-suppressed and T₂-filtered methods with in vitro diagnostic quantification of lipoproteins and metabolites. SARS-CoV-2 positive patients and healthy controls from both populations were modeled and cross-projected to estimate the biological similarities and validate biomarkers. Using the top 15 most discriminatory variables enabled construction of a cross-predictive model with 100% sensitivity and specificity (within populations) and 100% sensitivity and 82% specificity (between populations). Minor differences were observed between the control metabolic variables in the two cohorts, but the lipoproteins were virtually indistinguishable. We observed highly significant infection-related reductions in high-density lipoprotein (HDL) subfraction 4 phospholipids, apolipoproteins A1 and A2,that have previously been associated with negative regulation of blood coagulation and fibrinolysis. The Spanish and Australian diagnostic SARS-CoV-2 biomarkers were mathematically and biologically equivalent, demonstrating that NMR-based technologies are suitable for the study of the comparative pathology of COVID-19 via plasma phenotyping.

Proton nuclear magnetic resonance (NMR) N-acetyl signals (Glyc) from glycoproteins and supramolecular phospholipids composite peak (SPC) from phospholipid quaternary nitrogen methyls in subcompartments of lipoprotein particles) can give important systemic metabolic information, but their absolute quantification is compromised by overlap with interfering resonances from lipoprotein lipids themselves. We present a J-Edited DIffusional (JEDI) proton NMR spectroscopic approach to selectively augment signals from the inflammatory marker peaks Glyc and SPCs in blood serum NMR spectra, which enables direct integration of peaks associated with molecules found in specific compartments. We explore a range of pulse sequences that allow editing based on peak J-modulation, translational diffusion, and T₂ relaxation time and validate them for untreated blood serum samples from SARS-CoV-2 infected patients (n = 116) as well as samples from healthy controls and pregnant women with physiological inflammation and hyperlipidemia (n = 631). The data show that JEDI is an improved approach to selectively investigate inflammatory signals in serum and may have widespread diagnostic applicability to disease states associated with systemic inflammation.