Philipp
Nitschke
*a,
Samantha
Lodge
a,
Drew
Hall
a,
Hartmut
Schaefer
b,
Manfred
Spraul
b,
Nieves
Embade
c,
Oscar
Millet
c,
Elaine
Holmes
ad,
Julien
Wist
*ae and
Jeremy K.
Nicholson
*af
aAustralian National Phenome Centre and Computational and Systems Medicine, Health Futures Institute, Murdoch University, Harry Perkins Building, Perth, WA6150, Australia. E-mail: Philipp.Nitschke@murdoch.edu.au; Julien.Wist@murdoch.edu.au; Jeremy.Nicholson@murdoch.edu.au
bBruker Biospin GmbH, Rudolf-Plank Strasse 23, 76275 Ettlingen, Germany
cPrecision Medicine and Metabolism Laboratory, CIC bioGUNE, Parque Tecnológico de Bizkaia, Bld. 800, 48160, Derio, Spain
dDepartment of Metabolism, Digestion and Reproduction, Faculty of Medicine, Imperial College London, Sir Alexander Fleming Building, South Kensington, London SW7 2AZ, UK
eChemistry Department, Universidad del Valle, Cali 76001, Colombia
fInstitute of Global Health Innovation, Faculty of Medicine, Imperial College London, Level 1, Faculty Building, South Kensington Campus, London, SW7 2NA, UK
First published on 16th August 2022
A JEDI NMR pulse experiment incorporating relaxational, diffusional and J-modulation peak editing has been implemented for a low field (80 MHz proton resonance frequency) spectrometer system to measure quantitatively two recently discovered plasma markers of SARS-CoV-2 infection and general inflammation. JEDI spectra capture a unique signature of two biomarker signals from acetylated glycoproteins (Glyc) and the supramolecular phospholipid composite (SPC) signals that are relatively enhanced by the combination of relaxation, diffusion and J-editing properties of the JEDI experiment that strongly attenuate contributions from the other molecular species in plasma. The SPC/Glyc ratio data were essentially identical in the 600 MHz and 80 MHz spectra obtained (R2 = 0.97) and showed significantly different ratios for control (n = 28) versus SARS-CoV-2 positive patients (n = 29) (p = 5.2 × 10−8 and 3.7 × 10−8 respectively). Simplification of the sample preparation allows for data acquisition in a similar time frame to high field machines (∼4 min) and a high-throughput version with 1 min experiment time could be feasible. These data show that these newly discovered inflammatory biomarkers can be measured effectively on low field NMR instruments that do not not require housing in a complex laboratory environment, thus lowering the barrier to clinical translation of this diagnostic technology.
COVID-19 has presented major scientific, medical and socioeconomic challenges to the world that require multilevel solutions including new diagnostics and methods to evaluate the disease process20–23 and recovery from it.24,25 We previously described the complex metabolic phenoconversion process that characterises the multi-organ involvement of the SARS-CoV-2 infection process; e.g. SARS-CoV-2-(+) patients could be effectively distinguished from healthy controls modelling untargeted 1H NMR data (Fig. 1A). Phenoconversion from a healthy into an infected state is associated with changes in the concentrations of a wide range of metabolic entities (lipoproteins, glycoproteins, amino acids, lipids and other metabolites) that can be derived from NMR spectroscopic and mass spectrometric data,26–28 and which result in distinctive embedded biomarker features including some of those previously observed in diabetes,29,30 cardiovascular disease,31–35 liver dysfunction,36,37 neurological disruption38 and inflammation.39,40 Many of these pathological changes can be observed using a variety of NMR experiments on blood plasma and these changes are robust to sample handling with standardised protocols41,42 and were also observed by others.22,43–45
Fig. 1 Schematic showing the timeline and information summary for SPC and Glyc since their definition as candidate inflammatory markers of COVID-19 in 2020 (A and B), their subsequent quantitation in 2021 (C and D) and translation of detection methods for SPC and Glyc from 600 MHz to 80 MHz in 2022 (E and F). (A) OPLS of healthy controls (blue triangles) and SARS-CoV-2-positive patients using 1D 1H NMR input data effectively distinguishes between both groups (B) OPLS loadings show a high contribution and an inverse relationship between SPC and Glyc with SPC decreasing and Glyc increasing during COVID-19 due to acute inflammation. (A) and (B) was published in ref. 54 and is reproduced here under the terms of Creative Commons Attribution 4.0 licence. (C) 1D 1H with solvent suppression at 600 MHz (310 K) demonstrating the high peak convolution for SPC (red) and Glyc (blue) in an unedited experiment. (D) JEDI spectrum at 600 MHz (310 K) using relaxation, diffusion and J-editing to enable direct quantitation of SPC and Glyc by integration (E) 1D 1H NMR with solvent suppression at 80 MHz (298 K) showing the decreased signal dispersion and lower sensitivity compared to 600 MHz (C). At 80 MHz SPC is fully overlapped with Glucose peaks and Glyc fully overlapped with lipid resonances. (F) JEDI spectrum at 80 MHz (298 K) eliminates all overlapping peaks. More importantly, the lower dispersion at the benchtop field is not an issue, due to the extensive spectral editing of JEDI. |
Well resolved N-acetyl signals from glycosylated amino sugar residues from side-chains of acute phase reactive proteins such as α-1 N-acetyl-glycoprotein in NMR spectra of blood plasma were first reported by Bell et al.,46 and are elevated in multiple inflammatory states, such as obesity,47 diabetes,48 cardiovascular disease,49,50 rheumatoid arthritis,51 systemic immune-pathological conditions such as HIV infection,52 systemic lupus erythematosus53 and COVID-19.26,54 These signals are now referred to as GlycA and GlycB21 or just Glyc as the sum of the GlycA and GlycB. GlycA (δ 2.03 ppm) and GlycB (δ 2.07 ppm) are highly correlated and are postulated to be composites of N-acetyl signals from five main proteins: α-1-acid glycoprotein, α-1-antichymotrypsin, α-1-antitrypsin, haptoglobin and transferrin.53 We recently reported that Glyc signals are significantly elevated in COVID-19 patients and are strongly associated with infection positivity.26,27 In addition to Glyc, the supramolecular phospholipid composite (SPC) stemming from trimethylammonium residues in lipoprotein bound phospholipids was recently discovered as another inflammatory marker defining SARS-CoV-2 infection and SPC can also be used as a surrogate to gauge cardiovascular disease risk by analysis of its respective composite sub-regions.21,54 SPC was first quantified by DIffusion and Relaxation Editing (DIRE) NMR spectroscopy, which allows for integration of the peak by combining diffusion editing and T2-relaxation.21 Both candidate biomarkers SPC and Glyc were found to be distinguishing features under inflammatory conditions showing a distinct inverse relationship in the case of SARS-CoV-2 infection, with an increase of Glyc and a decrease of SPC (Fig. 1B), which can be expressed as the SPC/Glyc ratio for inflammation assessment.20
DIRE NMR spectroscopy was originally proposed as an experiment to enhance signals from molecules exhibiting slow translational diffusion, but with high segmental motional freedom.55 These requirements are satisfied by plasma glycoproteins and molecules constrained within certain lipoprotein sub-compartments.56,57 In a recent study we showed that measurement precision and sensitivity for both Glyc and SPC peaks could be improved using an additional J filter.20 This J-edited diffusion and relaxation (JEDI) NMR technique introduces destructive interference for signals with complex coupling interactions (Fig. 1C and D) thus minimising spectral interference and particularly enhances the quantitative integration accuracy of the Glyc peak(s). In addition, JEDI also yields increased signal to noise over DIRE due to shorter relaxation periods during the sequence. Our previous studies were performed on high field (600 MHz) IVDr NMR instruments that are widely used in clinical biochemistry studies to measure lipoprotein subfractions. Given the spectral simplification that JEDI experiments confer on complex mixtures, it makes them an ideal approach for low field NMR applications as they compensate both for the decreased dispersion and the more pronounced strong coupling effects at lower field. Adapted versions of JEDI experiments could provide comparable diagnostic capabilities, e.g. the SPC/Glyc ratio in reasonable NMR scanning times using more accessible low field NMR systems that are less costly and more readily field deployable.
Acquisition of the aliphatic region of a 1D 1H NMR spectrum (600 MHz, 310 K) of human serum allows for the detection of SPC (red) and Glyc (blue), but any quantitation of SPC and Glyc is hindered due to extensive convolution with small molecules, lipoproteins and protein background (Fig. 1C and S1A†). SPC and Glyc can be made spectroscopically accessible using JEDI spectroscopy, which combines diffusion, relaxation and J-editing to effectively remove all resonances interfering with SPC and Glyc. Fig. 1D shows the JEDI spectrum of 1C at 600 MHz, which enables direct quantitation of SPC and Glyc by simple integration. Apart from the SPC and Glyc signals, the JEDI spectrum is essentially void of other resonances. JEDI only contains residual lipoprotein peaks at δ 0.88, δ 1.28 and δ 2.77 ppm, the residual –CH– sugar residues belonging to the sugar moieties of N-acetylglucosamine (GlcNAc) and N-acetylneuraminic acid (Neu5Ac) of Glyc from δ 3.5–4.1 ppm, the choline moiety methylene groups at δ 3.69 and δ 4.32, and residual protein peaks in the aromatic region in addition to the main peaks SPC and Glyc (Fig. S1B†), whereas a 1D 1H NMR spectrum of serum usually contains over a hundred peaks.58,59 Considering that the JEDI experiment only refocuses uncoupled protons, i.e. singlets like SPC and Glyc, and given that the heavy pulse editing removes almost all resonances from the spectrum, JEDI presents an ideal candidate pulse sequence for application at lower fields as it intrinsically avoids two of the main issues encountered with benchtop NMR systems: namely, severely decreased spectral dispersion and more pronounced interference of strong coupling.60,61 This can be confirmed by acquisition of a regular 1D 1H NMR at 80 MHz (Fig. 1E) which also demonstrates the markedly lowered signal to noise (S/N) at lower field. Comparing the 80 MHz spectrum (Fig. 1E) to the 600 MHz equivalent (Fig. 1C) it is evident that the low field spectrum presents an amalgamation of broad features rather than many distinct peaks with a broad lipoprotein + protein background. SPC (red) and Glyc (blue) can still be gauged but are not just partially but indistinguishably overlapped with other resonances rendering attempts for relative quantification impractical. Other features like the main lipoprotein resonances, glucose and various other small molecules can still be detected, but also suffer from heavy convolution and/or low S/N (Fig. S1C†). Recording a spectrum with the JEDI method at 80 MHz gives a much clearer picture (Fig. 1F). Similar to the JEDI spectrum at 600 MHz (Fig. 1D), the spectrum shows SPC and Glyc to be the dominant features. As both the –+NCH3 residues of SPC and the N-acetyl residues of Glyc stem from uncoupled spin systems, they are unaffected by strong coupling effects at lower fields and are Fourier transformed into singlet peaks. In addition, due to the heavy editing of JEDI based on diffusion, relaxation and J-editing, almost all other resonances are eliminated compensating for the reduced dispersion at 80 MHz. This results in clean peaks for SPC and Glyc, which can be simply integrated for quantitation.
For the successful implementation and translation of the JEDI experiments from 600 to 80 MHz, various practical and theoretical considerations had to be taken into account. Currently the 80 MHz only allows for operation at a fixed temperature of 298 K, whereas serum or plasma measurements according to IVDr procedures are carried out at 310 K.62 Lower temperatures in serum lead to a change in chemical shifts,63 different coupling constants and a different lipoprotein profile (shape of the composite peaks), because the lipoprotein signal response depends on the mobility of the varying sub-components, which are further immobilised at reduced temperature.56 In addition, lower temperatures lead to smaller diffusion coefficients and relaxation times due to increased solvent viscosity,64 which are critical for JEDI editing. The effect of lower temperature on human serum can be clearly seen when comparing solvent suppressed 1H NMR spectra at 310 and 298 K (Fig. S2†). The signal responses from all macromolecular components i.e. protein background and lipoproteins are significantly reduced at 298 K. This includes the two biomarkers of SPC and Glyc, which stem from superimposed lipoprotein (SPC) and glycoprotein glycan (Glyc) peaks respectively.20,21 A more detailed picture of the effects of lowered temperature to SPC and Glyc can be taken from the JEDI spectra (Fig. S3†). Here, SPC is significantly reduced by ∼30% in peak intensity at 298 K compared to 310 K (Fig. S3A†). In addition, the different regions of SPC are affected differently. The high frequency part of SPC=(Fig. S3B†), which stems from LDL lipoproteins is almost completely eliminated at 298 K compared to the low frequency part of SPC, stemming from HDL lipoproteins.31 Similar to SPC, Glyc is also reduced by ∼20% in intensity at lower temperature (Fig. S3C†). Here, the intensity of the Glyc sub-regions, GlycA and GlycB are affected the same and only their frequencies are slightly altered with GlycB shifting to 2.4 Hz higher frequency with respect to GlycA. Notably, the residual lipoprotein peaks (–CH2CHCH–), which are overlapped with Glyc, are also strongly reduced at lower temperatures further improving the editing capabilities of JEDI (Fig. S3D†). The diffusional strength of JEDI was still found to be sufficient at 298 K and no small molecule resonances were detected.
For the translation to 80 MHz operating at 298 K, two possible JEDI pulse sequences were considered; the JEDI pulsed gradient spin echo (JEDI-PGSE; Fig. S4†) and the pulsed gradient perfect echo (JEDI-PGPE; Fig. S5†) (ESI:† basic pulse sequences). Both were introduced previously, with the JEDI-PGPE being deemed the go-to sequence at 600 MHz yielding a high S/N ratio combined with good lipoprotein suppression.20 The JEDI-PGSE yields a higher S/N, which is offset by reduced lipoprotein suppression and some baseline issues, but due to the lower sensitivity of benchtop NMR systems, S/N is one of the most crucial factors. Porting the methods to 80 MHz is straightforward and only requires adjustment of the time domain data points, due to the lower field (see Experimental section). Additionally, the previously proposed z-filter for both JEDI-PGSE and JEDI-PGPE had to be removed as the potency of the purge gradient was not strong enough on the benchtop system to remove a re-introduced water signal during the z-filter. The respective JEDI spectra of healthy control serum acquired with JEDI-PGSE and JEDI-PGPE are shown in Fig. S6.† Both yield clear peaks for SPC and Glyc void of obstructive overlap preventing analysis. The JEDI-PGSE demonstrated a higher S/N over the JEDI-PGPE (factor ∼1.6), which is in line with experiments at 600 MHz.20 In contrast, the JEDI-PGPE yielded a flatter baseline and omitted an antiphase residue from δ 3.4–3.8 ppm, which might make the JEDI-PGPE more robust to outliers and thus more suitable for high-throughput applications. The second biomarker, Glyc, can be detected as a lone but uneven singlet (Fig. S7†). This is due to the reduced frequency difference of the Glyc subregions GlycA and GlycB, which coalesce into a single peak at 80 MHz. Using a Gaussian apodization and longer experimental time (2k scans), both GlycA and GlycB can still be discerned at 80 MHz (Fig. S7†). In summary, both JEDI pulse sequences produce a spectrum with clear peaks for SPC and Glyc, that allow for straightforward quantification by integration.
The performance of the 80 MHz JEDI experiment at 298 K against its corresponding 600 MHz counterparts at 298 and 310 K were investigated using a cohort of SARS-CoV-2-(+) patients (n = 29) and healthy controls (n = 28) described elsewhere.54 The cohort is expected to show highly increased values for Glyc and decreased values for SPC in the COVID patients resulting in a lowered SPC/Glyc ratio compared to the healthy controls.
The JEDI spectra at 600 MHz were acquired with the JEDI-PGPE sequence at the two temperatures of 310 K, which is the standard temperature for IVDr acquisition and 298 K, which is the fixed operating temperature of the benchtop device. The 80 MHz spectra were acquired with the JEDI-PGPE and the JEDI-PGSE sequence. Comparison of the SPC/Glyc ratio at 600 MHz, 310 K and 80 MHz at 298 K (Fig. 2A) shows a clear distinction between healthy controls and SARS-CoV-2-(+) patients for both fields (p = 5.2 × 10−8 and 3.7 × 10−8 respectively) using the JEDI-PGSE experiment on the 80 MHz. A similar result is achieved comparing the distinction for 600 MHz, to 80 MHz both at 298 K (p = 3.7 × 10−8 and 3.7 × 10−8 respectively) (Fig. 2B). A linear fit of the data for 600 MHz, 310 K against 80 MHz indicated high correlation of the SPC/Glyc ratios (R2 of 0.93) at both fields (Fig. 2C) suggesting that the method can not only provide an overall distinction between COVID-(+) and healthy controls, but also the individual SPC/Glyc ratios per sample can be faithfully reproduced. A higher correlation (R2 of 0.97) was achieved when comparing the two fields at the same temperature (Fig. 2D). The higher correlation at 298 K is expected due to the slight spectral changes between 310 and 298 K like reduced/eliminated SPC-B contribution in SPC at lower temperatures, which leads to slight offsets when comparing 600 MHz at 310 K to 80 MHz at 298 K. Repeating the same experiments with the JEDI-PGPE experiment revealed a similar picture (Fig. S8†). For both temperatures, 310 and 298 K, a clear distinction was witnessed between SARS-Cov-2-(+) patients and healthy controls (Fig. S8A and B†). In addition, the linear fits of the SPC/Glyc ratios showed the same trend as for the JEDI-PGSE yielding a high correlation (R2 = 0.93) for 600 MHz, 310 K and (R2 = 0.97) for 600 MHz, 298 K against the 80 MHz SPC/Glyc ratios (Fig. S8C and D†). This confirms that the JEDI-PGPE is also a valuable option despite its lower S/N, while providing the option with better baseline stability and residual lipoprotein suppression for outliers. Although the linear fits for the SPC/Glyc ratio for 600 MHz against 80 MHz with JEDI-PGSE and JEDI-PGPE yielded high correlations, T1 relaxation changes according to the field strength increased with lower field. Moreover, the sequences for the low field application were slightly modified to fit the low field application. This explains the subtle variations in the mean SPC/Glyc ratio between high and lower fields observed in Fig. 2A and B.
After such promising results, we went on to investigate the impact of sample preparation and experimental parameters on 80 MHz spectra. Although, JEDI yields the same results at 600 and 80 MHz, samples prepared according to IVDr procedures (serum/plasma and buffer in a 1:1 ratio) lead to ∼15 min experimental times at 80 MHz as the lower S/N at low field has to be compensated by a higher number of scans (see section Materials and Methods). In order to reduce the measurement time to make high throughput applications feasible, we excluded the addition of a buffer to double the amount of plasma in the sample. Usually, a buffer is added for NMR serum analysis to ensure a high chemical shift reproducibility for accurate lipoprotein analysis and small molecule quantification; pH control might not be necessary because of the high chemical shift stability of SPC and Glyc. Basic NMR theory tells us that twice the amount of sample leads to a four times reduction of measurement time to achieve the same S/N ratio.65 This was confirmed by acquiring a JEDI spectrum for the same sample with (Fig. 3A) and without buffer dilution (Fig. 3B). As expected, both spectra yielded similar S/N ratios for SPC and Glyc with the advantage of the undiluted plasma reducing the possible measurement time for JEDI at 80 MHz to ∼4 min, while retaining sufficient signal for analyses. Furthermore, we shortened the measurement times until reaching the theoretical limit of quantification for JEDI experiments at 80 MHz, which is commonly accepted at S/N ≥ 10. Fig. 3C depicts a series of JEDI-PGSE spectra from a healthy control serum sample (i.e. SPC high, Glyc low) with increasing number of scans (no dummy scans) leading to a congruent experimental time. SPC as the more intensive peak reaches a SN ≥ 10 with as little as 2 scans or a measurement time of 7 s. A sufficient S/N for Glyc is achieved after 16 scans or 57 s resulting in a S/N of 28.71 for SPC under these conditions (Fig. 3C, green projection). Thus, 16 scans could suffice for the diagnostic test as Glyc usually increases and SPC (S/N of ∼30) is not expected to decrease by more than ∼30–40%, under inflammatory conditions. This suggests that an experiment time of under 1 minute might be feasible for adequate measurement of the SPC/Glyc ratio to determine inflammation in human serum. In addition, it should be noted that, as the 80 MHz system does not require individual sample shimming and due to the chemical shift stability of SPC and Glyc, the required experimental times will not be increased by additional “sample preparation” in the magnet and also extensive sample temperature equilibration might be unnecessary (see Fig. S9†). Hence, including the time for 90° pulse calibration (∼10 s) and sample exchange by robotics (∼75 s for combined injection and ejection) a real time sample throughput of ∼1 sample per 2.5 min could be envisioned. Repeating the procedure of increasing number of scans with the JEDI-PGPE experiment showed the expected result of a lower S/N experiment (Fig. S10†). Instead of reaching sufficient S/N ≥ 10 for SPC and Glyc after 16 scans in under one minute, the JEDI-PGPE requires 32 scans or two minutes for comparable results.
Footnote |
† Electronic supplementary information (ESI) available: Fig. S1: Comparison of solvent suppressed proton spectra and JEDI for a healthy control at 600 and 80 MHz. Fig. S2: Aliphatic 600 MHz 1H NMR region of serum at 310 and 298 K. Fig. S3: 600 MHz JEDI spectra of SPC and Glyc at 310 and 298 K. Fig. S4: Sequence of Pulsed Gradient Spin Echo (PGSE). Fig. S5: Sequence of Pulsed Gradient Perfect Echo (PGPE). Fig. S6: JEDI-PGSE and JEDI-PGPE of serum at 80 MHz. Fig. S7: Comparison of the Glyc region in a serum sample at 80 MHz and 600 MHz by JEDI. Fig. S8: Statistical analyses comparing SPC/Glyc JEDI-PGPE measurements at 600 MHz 310 K (IVDr) and 600 MHz 298 K with SPC/Glyc measurements at 80 MHz 298 K. Fig. S9: Chemical shift sensitivity of SPC and Glyc in a serum sample at 80 MHz. Fig. S10: Effects of varying plasma concentration and number of scans/experimental time on JEDI-PGPE at 80 MHz. See DOI: https://doi.org/10.1039/d2an01097f |
This journal is © The Royal Society of Chemistry 2022 |