Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Cerebral phosphoester signals measured by 31P magnetic resonance spectroscopy at 3 and 7 Tesla

  • Shizhe Li ,

    Roles Data curation, Formal analysis, Methodology, Writing – original draft

    steveli@mail.nih.gov

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • Jan Willem van der Veen,

    Roles Data curation, Formal analysis, Methodology, Writing – review & editing

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • Li An,

    Roles Data curation, Formal analysis, Methodology

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • JoEllyn Stolinski,

    Roles Data curation, Resources

    Affiliation National Institute of Neurological Disorders and Stroke, Bethesda, MD, United States of America

  • Christopher Johnson,

    Roles Resources

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • Maria Ferraris-Araneta,

    Roles Project administration, Resources

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • Milalynn Victorino,

    Roles Resources

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • Jyoti Singh Tomar,

    Roles Resources

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

  • Jun Shen

    Roles Conceptualization, Formal analysis, Project administration, Writing – review & editing

    Affiliation National Institute of Mental Health, Bethesda, MD, United States of America

Abstract

Abnormal cell membrane metabolism is associated with many neuropsychiatric disorders. Free phosphomonoesters and phosphodiesters, which can be detected by in vivo 31P magnetic resonance spectroscopy (MRS), are important cell membrane building blocks. However, the quantification of phosphoesters has been highly controversial even in healthy individuals due to overlapping signals from macromolecule membrane phospholipids (MP). In this study, high signal-to-noise ratio (SNR) cerebral 31P MRS spectra were acquired from healthy volunteers at both 3 and 7 Tesla. Our results indicated that, with minimal spectral interference from MP, the [phosphocreatine (PCr)]/[phosphocholine (PC) + glycerophosphocholine (GPC)] ratio measured at 7 Tesla agreed with its value expected from biochemical constraints. In contrast, the 3 Tesla [PCr]/[PC+GPC] ratio obtained using standard spectral fitting procedures was markedly smaller than the 7 Tesla ratio and than the expected value. The analysis suggests that the commonly used spectral model for MP may fail to capture its complex spectral features at 3 Tesla, and that additional prior knowledge is necessary to reliably quantify the phosphoester signals at low magnetic field strengths when spectral overlapping is significant.

Introduction

31P magnetic resonance spectroscopy (MRS) allows the evaluation of several important aspects of brain energetics and metabolism, from levels of 31P-containing metabolites to metabolic fluxes associated with adenosine triphosphate (ATP)-generating enzymes. Since the 1970s, many studies have also used 31P MRS and functional 31P MRS techniques to investigate a broad range of both normal brain functions and neuropsychiatric disorders [15]. These 31P MRS studies have provided considerable insight into the role of mitochondrial metabolism and cell membrane synthesis in normal physiology and disease development as well as therapeutic interventions.

Phospholipids are known to have a phosphodiester (PDE) structure whose major resonances overlap with the resonances of free PDE and phosphomonoester (PME). The in vivo phospholipid signals are very pronounced at low magnetic fields [6, 7]. These signals, often referred to as membrane phospholipids (MP), originate from 31P-containing macromolecules, including phospholipid bilayers. They decrease with increasing magnetic field strength [79]. At 7 Tesla and higher magnetic field strengths, the signals of glycerophosphocholine (GPC), glycerophosphoethanolamine (GPE), phosphocholine (PC), and phosphoethanolamine (PE) are well resolved, with little or no spectral interference from the much reduced MP signals [8, 9]. However, at the more commonly used 3 Tesla or lower field strengths, large MP signals have consistently been observed [10, 11]. Typically, spectral fitting techniques are necessary to extract the individual components of PDE, PME, and inorganic phosphates (Pi). With proton decoupling, the MP signals at 1.5 Tesla were found to be broad and highly asymmetrical [10]. High resolution 31P nuclear magnetic resonance (NMR) studies of lipid extracts of animal brain also revealed multiple 31P signals resonating at different frequencies [12]. Nevertheless, it has been customary to use a single symmetrical peak centered at 2.3 ppm to fit the in vivo MP resonances. This is of concern because many pathological conditions involving abnormal cell membrane metabolism may alter the composition of the MP signals and, by inference, their lineshape [13].

The detrimental effects of the strong macromolecule signals on accurate quantification of short echo time proton MRS spectra using the LCModel have recently been recognized [14]. Because the separation between small metabolite molecules and macromolecules by spectral fitting relies on differences in their linewidths, the broad linewidths of the composite PDE and PME signals at 1.5–3 Tesla may make it difficult to reliably remove spectral contamination from MP. As an example, the reported PDE/ATP ratio in brain ranged from 0.4 to 4.0, as measured in healthy adult volunteers using 1.5–3 Tesla [10, 1520]; a span of this order of magnitude cannot be reconciled by potential differences in MRS techniques, T1 saturation, nuclear Overhauser enhancement, or variations in spatial distribution of phosphoesters.

In this work we acquired high signal-to-noise ratio (SNR) cerebral 31P MRS spectra from healthy subjects at both 3 and 7 Tesla. We found that the [PCr]/[PC+GPC] ratio extracted from 7 Tesla 31P MRS spectra of healthy subjects is consistent with biochemical constraints from the literature. In comparison, the 3 Tesla [PCr]/[PC+GPC] ratio obtained using standard spectral fitting procedure is significantly smaller than the expected value. Our results suggest that the commonly used spectral model of MP may not be suitable for quantifying the overlapping phosphoester signals at low magnetic field strengths.

Materials and methods

The human study was performed under the protocol NCT01266577 that was approved by the National Institutes of Health (NIH) Institutional Review Board (IRB). Written consent forms were obtained from all participants.

7 Tesla 31P MRS

In vivo 31P MRS experiments were performed on a Siemens Magnetom 7 Tesla scanner (Siemens Healthcare, Erlangen, Germany). A 1H and 31P coil assembly was built in-house (Fig 1A), comprising a circular 31P coil (diameter = 7.0 cm), a quadrature half-volume proton coil, and a slotted radio frequency (RF) shield. These RF devices were mounted on three vertically stacked, thin, semi-cylindrical plastic tubes. A 31P surface coil was mounted on the outer surface of the upper tube (outer diameter = 20.3 cm). A 1H frequency blocking L-C tank circuit was placed at the middle point of the 31P circular loop. The proton coil, which consisted of two overlapped octagon loops with nominal length/width of 12.3 cm, was mounted on the outer surface of the middle tube (outer diameter = 20.3 cm), and the slotted RF shield was mounted on the inner surface of the lower tube (inner diameter = 22.2 cm). Each proton loop had a single-tuned 1H cable trap constructed using an RG-316 cable. A 1H/31P dual-tuned cable trap, built inside an RF-shielded box, was connected to the 31P coil. The coil assembly was connected to the 7 Tesla scanner via an interface box (Quality ElectroDynamics, Mayfield Village, Ohio, USA) that included transmit-receive switches, pre-amplifiers, RF filters for both channels, and a quadrature combiner for the proton channel.

thumbnail
Fig 1. A 7 Tesla radio frequency (RF) coil assembly and its relative position to the human head.

(A) The circular 31P surface coil, the 1H coil (two overlapping octagonal loops), and the RF shield were mounted on the bottom surface of the upper, middle, and lower tubes, respectively. Each proton loop had a 1H cable trap constructed using an RG-316 cable. A 31P cable trap, built inside an RF shielded box, was connected to the 31P coil. (B) Gradient echo localizer image from a healthy volunteer showing the relative position of the 31P surface coil (red arc, 7.5 cm dia.) and the shim voxel (5 x 5 x 5 cm3 box in while color).

https://doi.org/10.1371/journal.pone.0248632.g001

31P MRS data were acquired from six healthy volunteers (mean age: 30 years; age range: 22–48 years). A voxel of 5 x 5 x 5 cm3 located right above the 31P coil (Fig 1B) was selected for performing B0 shimming in the occipital lobe using the Siemens 3D shim tool with correction of all first- and second-order shims and five third-order shims. Shimming results were evaluated by measuring the water linewidth (full width at half height) from the 5 x 5 x 5 cm3 voxel using proton MRS with a point resolved spectroscopy (PRESS) sequence. Water linewidth was found to be 13.9 ± 1.2 Hz (n = 6). 31P spectra were acquired using the Siemens excite-acquire sequence: hard pulse = 250 μs, number of data points = 1024, and spectral width (SW) = 5 kHz. The transmit frequency was set to the frequency of PCr at 0 ppm. Two repetition time (TR) were used to evaluate the effect of T1 saturation. At TR = 3 seconds, the number of averages (NA) was 64, and at TR = 30 seconds, the NA was 32 (n = 6 for both TRs). The transmit voltage for 31P was pre-optimized in vivo to obtain maximum signal-to-noise ratio for a given TR.

3 Tesla 31P MRS

In vivo 31P MRS scans were also performed on a Siemens Skyra 3 Tesla scanner (Siemens Healthcare, Erlangen, Germany). A 3 Tesla 1H and 31P coil assembly was built in-house that had the same geometrical layout as the 7 Tesla coil assembly described above except that the coils were tuned to 3 Tesla frequencies and there was no RF shield. The coil assembly was connected to the 3 Tesla scanner via an interface box (Stark Contrast MRI Coils Research, Erlangen, Germany). Eight heathy volunteers (mean age = 29 years, age range: 21–45 years) were recruited and scanned. The same shimming procedure used in the 7 Tesla MRS was applied to optimize field homogeneity of a 4 x 4 x 4 cm3 voxel prescribed above the 31P surface coil without the third order corrections. Water linewidth was found to be 15.5 ± 0.9 Hz (n = 8). 31P spectra were acquired using the Siemens excite-acquire sequence: hard pulse = 500 μs, number of data points = 1024, and SW = 5 kHz. The transmit frequency was set to the frequency of PCr at 0 ppm. At TR = 2 seconds, the NA was 128, and at TR = 25 seconds, the NA was 64 (n = 8 for both TRs). The 31P RF power at 3 Tesla was calibrated in the same way as in the 7 Tesla experiments.

31P MRS spectral fitting

Raw data were first preprocessed using an IDL (Harris Geospatial Solutions, Boulder, CO, USA) program to zero the first two (TR = 2, 3 seconds) or three (TR = 25, 30 seconds) complex data points in free induction decay for baseline suppression. The chemical shift of phosphocreatine (PCr) was set to 0 ppm. The data were subsequently phase-corrected and saved in the Siemens RDA raw data format. The RDA data were then processed using the jMRUI software package with the Amares routine [21] and prior spectral knowledge of 31P metabolites [2124]. The 31P basis set consisted of α-, β-, γ-adenosine triphosphate (α-, β-, γ-ATP), nicotinamide adenine dinucleotide (NAD), PCr, MP, GPC, GPE, intracellular inorganic phosphate (Pi), extracellular inorganic phosphate (Pi extra), PC, PE, and uridine diphosphate glucose (UDPG). Four different spectral models of MP (Lorentzian, Gaussian lineshapes, fixed and freely adjustable linewidths) were evaluated to assess the effects of MP on 31P MRS data quantification. For MP fitting with fixed linewidths, the fixed linewidth values were calculated from the median of the freely fitted MP linewidths.

[PCr]/[GPC+PC]

The choline singlet at 3.2 ppm observed in the proton MRS spectrum of healthy volunteers originates almost entirely from the trimethylamine group of PC and GPC [10, 25], given that the concentrations of choline, acetylcholine, and phosphatidylcholine in their free forms are negligibly low at the sensitivity level of proton MRS [2628]. That is, (1) where tCho is total choline. The ratio of total creatine (tCr) to total choline ([tCr]/[tCho]) and [PCr]/[GPC+PC] can be obtained from proton and 31P MRS data, respectively, without absolute quantification. Because of Eq 1, (2) The creatine kinase equilibrium is tightly regulated, with [PCr] ≈ [Cr] predicted by the creatine kinase equilibrium constant and in vitro assays [10, 29]. The methylene proton signals of PCr and creatine (Cr) resonate at 3.93 and 3.91 ppm, respectively, in the proton MRS spectra. MRS studies of rodent brain performed at very high magnetic field strengths (≥ 9.4 Tesla) provided visual confirmation that the spectrally resolved PCr and Cr methylene protons were of similar intensity. Previous studies noted that the [PCr]/[tCr] ratio in the cerebral cortex was 0.46–0.59 from the spectrally resolved PCr and Cr in proton MRS spectra acquired at 9.4–16.4 Tesla [3032], with an average value of 0.53. From the above analysis, the expected [PCr]/[GPC+PC] ratio is given by (3) Because the reported phosphoester levels of healthy volunteers in the 31P MRS literature vary over an order of magnitude (0.4–4.0) [10, 1520], the relationship given by Eq 3 may be used as an approximate internal standard to check 31P spectral fitting of data acquired from healthy volunteers.

Results

7 Tesla 31P MRS

Typical 7 Tesla in vivo 31P spectra from the occipital lobe of human brain are shown in Fig 2, where spectrum (A) was acquired with TR = 3 seconds and NA = 64, and (B) with TR = 30 seconds and NA = 64. Resonances of PE (6.8 ppm), PC (6.2 ppm), Pi (4.8 ppm), GPE (3.5 ppm), GPC (2.9 ppm), MP (2.3 ppm), PCr (0 ppm), γ-ATP (-2.5 ppm), α-ATP (-7.6 ppm), NAD, UDPG Pβ (-9.8 ppm), and β-ATP (-16.2 ppm) were detected. The signal of oxidized NAD (NAD+, -8.2 and -8.4 ppm) and reduced form of NAD (NADH, -8.1 ppm) were collectively labeled as NAD. Notably, the MP signals at 2.3 ppm were quite weak and spectrally resolved from GPC at 2.9 ppm. The intensity variations between the two TR settings were expected from the known T1 values of 31P-containing metabolites [33]. An example of jMRUI fitting of a 7 Tesla 31P spectrum acquired with TR = 30 seconds is shown in Fig 3, where the MP signals at 2.3 ppm were modeled using a Gaussian curve with freely adjustable linewidth. Table 1 lists the corresponding concentrations of 31P-containing metabolites (normalized to [γ-ATP]) extracted from jMRUI fitting for both TR = 3 seconds and 30 seconds (see S1S4 Figs in Supplementary Materials for additional examples of jMRUI fitting of 7 Tesla spectra and a comparison of spectral analysis results using the four different spectral models of MP described in Materials and Methods). From the jMRUI analysis results given in Table 1, [PCr]/[PC+GPC] was found to be 2.4 ± 0.3 (TR = 3 seconds, n = 6) and 2.6 ± 0.3 (TR = 30 seconds, n = 6). The [MP]/[γ-ATP] ratio was found to be 0.13 ± 0.09 (TR = 3 seconds, n = 6) and 0.03 ± 0.04 (TR = 30 seconds, n = 6), consistent with the small MP signals observed at 2.3 ppm in Fig 2.

thumbnail
Fig 2. Typical 7 Tesla 31P spectra of a healthy volunteer.

A 7 Tesla spectrum at (A) TR = 3 seconds, NA = 64 and (B) TR = 30 s, NA = 32 x 2. PE: phosphoethanolamine; PC: phosphocholine; Pi: inorganic phosphate; GPE: glycerophosphoethanolamine; GPC: glycerophosphocholine; MP: membrane phospholipids; PCr: phosphocreatine; α-, β-, γ-ATP, α-, β-, γ-adenosine triphosphate; NAD: nicotinamide adenine dinucleotide; UDPG: uridine diphosphate glucose.

https://doi.org/10.1371/journal.pone.0248632.g002

thumbnail
Fig 3. jMRUI analysis of a typical 31P spectrum acquired at 7 Tesla.

The original spectrum was from Fig 2B with TR = 30 seconds, NA = 64. The individual spectral components, the overall fit and fit residual were also shown. A Gaussian curve with freely adjustable linewidth was used to fit the MP signal at 2.3 ppm. PCr: phosphocreatine; α-, β-, γ-ATP: α-, β-, γ-adenosine triphosphate; NAD+: oxidized nicotinamide adenine dinucleotide; UDPG: uridine diphosphate glucose; NADH: reduced form of nicotinamide adenine dinucleotide; Pi: inorganic phosphate; Pi extra: extracellular inorganic phosphate; GPE: glycerophosphosethanolamine; GPC: glycerophosphocholine; PE: phosphoethanolamine; PC: phosphocholine; MP: membrane phospholipids.

https://doi.org/10.1371/journal.pone.0248632.g003

thumbnail
Table 1. 31P metabolite concentration ratios determined at 7 Tesla (n = 6).

https://doi.org/10.1371/journal.pone.0248632.t001

3 Tesla 31P MRS

Typical 3 Tesla 31P spectra acquired from a healthy volunteer are shown in Fig 4 at TR = 2 seconds (A) and 25 seconds (B); each spectrum was averaged with NA = 128. Resonances of PE, PC, Pi, GPE, GPC, PCr, NAD, and α-, β-, and γ-ATP were detected although the individual components of PME and PDE were not spectrally resolved. In contrast to the nearly flat baseline at 7 Tesla seen in Figs 2 and 3, a prominent baseline was observed in the 1–7 ppm region in the 3 Tesla 31P MRS spectra (Fig 4). From TR = 2 seconds to TR = 25 seconds, the relative increase in signal intensities of 31P metabolites measured at 3 Tesla was similar to the 7 Tesla results. An example of jMRUI fitting of a 3 Tesla 31P spectrum acquired with TR = 25 seconds, where the MP signal at 2.3 ppm was also modeled using a Gaussian curve, is shown in Fig 5. The corresponding concentrations of 31P-containing metabolites (normalized to [γ-ATP]) extracted from jMRUI fitting of the 3 Tesla spectra are listed in Table 2 for both TR = 2 seconds and 25 seconds. The [PCr]/[PC+GPC] ratio extracted from our 3 Tesla data was 1.1 ± 0.1 (TR = 2 seconds, n = 8) and 1.3 ± 0.3 (TR = 25 seconds, n = 8). The fitted [MP]/[γ-ATP] ratio was found to be 0.77 ± 0.27 (TR = 2 seconds, n = 8) and 0.43 ± 0.22 (TR = 25 seconds, n = 8) at 3 Tesla, which is substantially greater than the values detected at 7 Tesla. Notably, the 3 Tesla [PCr]/[PC+GPC] ratio extracted from spectral analysis using LCModel [24] was similar to our jMRUI results obtained at the same field strength. For additional examples of 3 Tesla spectral fitting and a comparison of spectral analysis results using the four different spectral models of MP at 3 Tesla see S5S8 Figs in Supplementary Materials.

thumbnail
Fig 4. Typical 3 Tesla 31P spectra of a healthy volunteer.

A 3 Tesla spectrum at (A) TR = 2.5 seconds, NA = 128 and (B) TR = 25 seconds, NA = 128. PE: phosphoethanolamine; PC: phosphocholine; Pi: inorganic phosphate; GPE: glycerophosphoethanolamine; GPC: glycerophosphocholine; MP: membrane phospholipids; PCr: phosphocreatine; α-, β-, γ-ATP: α-, β-, γ-adenosine triphosphate; NAD: nicotinamide adenine dinucleotide.

https://doi.org/10.1371/journal.pone.0248632.g004

thumbnail
Fig 5. jMRUI analysis of a typical 31P spectrum acquired at 3 Tesla.

The original spectrum was from Fig 4B with TR = 25 seconds, NA = 64. The individual spectral components, the overall fit and fit residual were also shown. A Gaussian curve with freely adjustable linewidth was used to fit the MP signal at 2.3 ppm. PCr: phosphocreatine; α-, β-, γ-ATP: α-, β-, γ—adenosine triphosphate; NAD+: oxidized nicotinamide adenine dinucleotide; UDPG: uridine diphosphate glucose; Pi: inorganic phosphate; Pi extra: extracellular inorganic phosphate; GPE: glycerophosphosethanolamine; GPC: glycerophosphocholine; PE: phosphoethanolamine; PC: phosphocholine; MP: membrane phospholipids.

https://doi.org/10.1371/journal.pone.0248632.g005

thumbnail
Table 2. 31P metabolite concentration ratios determined at 3 Tesla (n = 8).

https://doi.org/10.1371/journal.pone.0248632.t002

Fig 6 compares [31P metabolite]/[γ-ATP] ratios measured at 3 and 7 Tesla with long TRs and freely adjustable linewidths for MP. Student t-tests were performed for the values averaged over the four spectral models of MP (see Table 3). Although there was a general agreement among the upfield signals (PCr, α-, β-, γ-ATP, and NAD) across field strengths and MP lineshape models, significant differences were seen between 3 and 7 Tesla downfield [31P metabolite]/[γ-ATP] ratios. The differences between the fitted 3 and 7 Tesla downfield [31P metabolite]/[γ-ATP] ratios remain highly significant after Bonferroni corrections of the p values are made.

thumbnail
Fig 6. Comparison of [31P metabolite]/[γ-ATP] ratios measured at 3 Tesla and 7 Tesla.

3 Tesla spectra (TR = 25 seconds) and 7 Tesla spectra (TR = 30 seconds) were fitted with freely adjustable linewidths for MP. Blue: 3 Tesla, fitted using Gaussian lineshape; orange: 3 Tesla, fitted using Lorentzian lineshape; grey: 7 Tesla, fitted using Gaussian lineshape; yellow: 7 Tesla, fitted using Lorentzian lineshape.

https://doi.org/10.1371/journal.pone.0248632.g006

thumbnail
Table 3. Student t-test results for 3 Tesla (n = 8) and 7 Tesla (n = 6) 31P metabolite ratios.

https://doi.org/10.1371/journal.pone.0248632.t003

Comparison with Eq 3

The total creatine to total choline ratio [tCr]/[tCho] has been well documented in a large number of in vivo proton MRS and chemical shift imaging studies, and results from various studies broadly agree [34, 35]. Because the volume of brain tissue measured by 31P MRS is usually large, we averaged tCr and tCho values from the existing literature [34, 35] and arrived at [tCr]/[tCho] ≈ 4.8 for cerebral cortex of healthy adult volunteers (see Table 4). The [PCr]/[GPC+PC] ratio expected from Eq 3 is therefore approximately 0.53*4.8 = 2.5.

thumbnail
Table 4. Literature values used to calculate the expected cerebral [PCr]/[GPC+PC] ratio using Eq 2.

https://doi.org/10.1371/journal.pone.0248632.t004

In short echo proton MRS, spectral fitting of the macromolecule baseline using LCModel was found to be strongly influenced by SNR and linewidth [14]. Despite the large number of 31P MRS studies reported in the literature, our comparison was limited to studies of healthy volunteers that measured PCr, GPC, and PC intensities using a large voxel or localization by surface coil for maximum SNR [10, 33, 36, 37]. General agreement was observed between the experimental [PCr]/[PC+GPC] ratios measured at 7 Tesla and that expected from Eq 3 (Table 5). Thus, it appears that one of the key advantages of 7 Tesla is that PE, PC, GPE, GPC, and MP can be spectrally resolved due to the dramatically reduced MP intensity and increased chemical shift dispersion at the much higher magnetic field strength.

For lower magnetic field strengths, proton decoupling is needed to improve spectral resolution of 31P MRS to resolve the phosphoesters acquired from large tissue volumes [10]. Due to safety concerns and nonstandard hardware requirements, few previous studies attempted proton decoupling. In particular, the 1.5 Tesla studies listed in Table 5 employed proton decoupling, nuclear Overhauser enhancement, and MP baseline removal with PE, PC, GPE, and GPC peaks fully resolved. The [PCr]/[PC+GPC] ratios obtained from those 1.5 Tesla studies agreed reasonably well with our 7 Tesla and Eq 3 results.

Although spectral fitting can extract individual components of MRS signals, sufficient prior knowledge is crucial when significant spectral overlapping is present. Because the ground truth lineshape of MP at 3 Tesla is unknown, a single symmetrical MP basis peak centered at 2.3 ppm is unlikely to fully capture the complex MP signals. The 3 Tesla [PCr]/[PC+GPC] ratios extracted using both LCModel [24] and jMRUI ([24], and Tables 2 and 5 of the current study) were ~50% of the value expected from Eq 3; thus, both are markedly smaller than the [PCr]/[PC+GPC] ratios determined from studies with spectrally resolved phosphoesters (Tables 1 and 5).

Discussion

This study, which found that the [PCr]/[PC + GPC] ratio measured at 7 Tesla, but not 3 Tesla, agreed with values expected from biochemical constraints. The results suggest that a single symmetrical peak may fail to capture the complex spectral features of MP at 3 Tesla, and that additional prior knowledge is necessary to reliably quantify the downfield 31P signals at low magnetic field strengths when spectral overlapping is significant.

The first two or three data points of the FIDs were discarded in this study, a practice commonly used in processing 31P MRS data in order to suppress the broad baseline (e.g., [33]). Since the spectral width of our study is 5 kHz, this is equivalent to adding a time delay of 400–600 μs before recording the FIDs. This procedure has no significant effect on quantification of MP as the time delay corresponds to signals with kHz linewidths.

Abnormal levels of PME and PDE—both of which are cell membrane building blocks—have been associated with many neuropsychiatric disorders [38, 39]. In particular, the PDE/PME ratio was found to be elevated in cancer patients and to fall with effective treatment [5]; the PDE/PME ratio has also been considered a marker of MP turnover equilibrium, which is altered in many diseases, including cancer [10, 38]. Furthermore, cell line and tumor model studies have found that PC/GPC and PE/GPE ratios, which appear to be significantly higher in primitive neuroectodermal tumors [4], may serve as indicators of malignancy.

Despite the importance of these markers, 31P MRS studies in many neuropsychiatric disorders that involve phosphoesters have obtained mixed results. For example, many 31P MRS studies reported disturbances in phospholipid metabolism in Alzheimer’s disease that correlated with pathological biomarkers [40]; in contrast, other studies [1] found no changes in metabolite concentrations or ratios in patients compared with controls in whole axial sections of the brain and no correlation between 31P MRS indexes and severity of dementia. Notably, it is likely that contamination by MP has contributed to many of the mixed findings.

In the present study, the low [PCr]/[PC+GPC] ratio measured at 3 Tesla is consistent with significant spectral contamination by MP that inflates its denominator while PCr is affected very little due to its narrow linewidth and to the large resonant frequency separation between PCr and MP. Therefore, when the [PCr]/[PC+GPC] ratio extracted from the 31P MRS spectra of healthy volunteers is conspicuously smaller than the value expected from Eq 3, MP contamination should be suspected. Caution is needed if the same spectral fitting procedure is to be applied to data acquired from patients since their MP lineshape may have been altered under pathological conditions. It also should be noted that Eq 3 is derived for healthy subjects. When data from healthy controls are unavailable additional experiments including complementary 1H MRS studies of patients are needed to validate spectral fitting procedures at 3 Tesla or lower magnetic field strengths.

The advent of 7 Tesla clinical scanners provides an opportunity to revisit many of the early controversies because not only is MP signal intensity greatly reduced at 7 Tesla, it is also spectrally resolved from GPE and GPC. Furthermore, PE, PC, GPC, and GPE are all resolved at 7 Tesla with reasonable B0 shimming without the need to perform proton decoupling. Despite these advantages, most 31P clinical studies are expected to continue to be conducted on the more prevalent 3 Tesla scanners and without proton decoupling for the foreseeable future. Importantly, our analysis of spectral fitting with symmetrical MP models suggests that additional prior knowledge is required to avoid significant contamination of downfield signals, as evidenced by the lower [PCr]/[PC+GPC] ratios obtained from both LCModel [24] and jMRUI (see Tables 2 and 5) analyses of 3 Tesla 31P data.

In addition to potential contamination by MP, other factors may also affect metabolite quantification by 31P MRS. For example, nuclear Overhauser enhancement is not uniform across the 31P signals [41]; specifically, in contrast to the nearly uniform T1 of most metabolites detected by proton MRS [42], there is a large dispersion in T1s of 31P MRS signals [33, 43], making their intensities sensitive to the commonly used TR values. Such variations may also have contributed to the different quantification results reported in the 31P MRS literature. Furthermore, although it is relatively easy to use the same MRS parameter settings for control and patient groups to reduce variability, MP contamination is more worrisome as many diseases may also alter MP signals, as may specific changes in cell membrane metabolism [13]. Therefore, it is important to minimize or eliminate spectral contamination from MP signals when free phosphoesters are of clinical interest.

Many metabolite ratios reported in Tables 1 and 2 may not be generalized because the current study used surface transceiver coils to acquire 31P MRS spectra from a small number of human subjects and no gradient-based spatial localization techniques were employed. Because of the severe B1 inhomogeneity associated with surface coils the 31P signals are expected to be affected by many factors specific to our experiments such as head size, size of the surface coil, RF calibration, and placement of transmitter frequency. For example, many ratios reported in Table 1 are significantly different from the corresponding values reported by a similar 7 Tesla study that used a larger surface coil [33]. Nevertheless, it should be noted that the [PCr]/[PC+GPC] ratio is insensitive to spatial location of 31P signals for the same tissue type due to the intrinsic biological constraints contained in Eq [3]. As such it is expected that this internal standard is relatively immune to many variations originated from the use of surface transceiver coils.

Taken together, the results of the present study suggest that contamination by MP signals is minimal at 7 Tesla because of the diminished intensity of MP and increased spectral resolution at high magnetic fields. In contrast, the free phosphoester signals in 31P MRS spectra acquired at low magnetic fields (≤ 3 Tesla) overlap with the large MP signals. Although MP lineshape has been revealed to be highly asymmetrical at 1.5 Tesla with proton decoupling [35], the true lineshape and spectral pattern of the complex MP signal at 3 Tesla remain undetermined. The biochemical constraints described here suggest that fitting the MP signals at 3 Tesla using a symmetrical peak centered at 2.3 ppm leads both to its underestimation and to a significant overestimation of phosphoesters. In contrast, the [PCr]/[PC+GPC] ratio expected from Eq 3 agrees reasonably well with 1.5 Tesla data acquired with MP removal and 7 Tesla data, lending confidence to characterization of membrane metabolism in brain disorders using 7 Tesla 31P MRS without proton decoupling.

Conclusion

When spectral interference from MP is minimized, the experimentally measured [PCr]/[PC+GPC] ratios were found to be in reasonable agreement with the value expected from biochemical constraints. The common practice of fitting the complex MP signal by a single symmetrical peak can lead to overestimation of phosphoester signals at low magnetic field strengths including 3 Tesla.

Supporting information

S1 Fig. 31P spectral fitting of the 7 Tesla spectrum shown in Fig 2B with TR = 30 seconds.

A Gaussian curve with fixed linewidth of 160 Hz was used to fit the MP signal.

https://doi.org/10.1371/journal.pone.0248632.s001

(TIF)

S2 Fig. 31P spectral fitting of the 7 Tesla spectrum shown in Fig 2B with TR = 30 seconds.

A Lorentzian curve with a freely adjustable linewidth was used to fit the MP signal.

https://doi.org/10.1371/journal.pone.0248632.s002

(TIF)

S3 Fig. 31P spectral fitting of the 7 Tesla spectrum shown in Fig 2B with TR = 30 seconds.

A Lorentzian curve with a fixed linewidth of 320 Hz was used to fit the MP signal.

https://doi.org/10.1371/journal.pone.0248632.s003

(TIF)

S4 Fig. 31P metabolite ratios extracted from 7 Tesla spectra at TR = 30 seconds.

The MP signals were fitted using the four spectral models: Blue: Gaussian with a fixed linewidth of 160 Hz; orange: Gaussian with a freely adjustable linewidth; grey: Lorentzian with a fixed linewidth of 320 Hz; yellow: Lorentzian with a freely adjustable linewidth.

https://doi.org/10.1371/journal.pone.0248632.s004

(TIF)

S5 Fig. 31P spectral fitting of the 3 Tesla spectrum shown in Fig 4B with TR = 25 seconds.

A Gaussian curve with a fixed linewidth of 100 Hz was used to fit the MP signals.

https://doi.org/10.1371/journal.pone.0248632.s005

(TIF)

S6 Fig. 31P spectral fitting of the 3 Tesla spectrum shown in Fig 4B with TR = 25 seconds.

A Lorentzian curve with a freely adjustable linewidth was used to fit the MP signals.

https://doi.org/10.1371/journal.pone.0248632.s006

(TIF)

S7 Fig. 31P spectral fitting of the 3 Tesla spectrum shown in Fig 4B with TR = 25 seconds.

A Lorentzian curve with a fixed linewidth of 200 Hz was used to fit the MP signals.

https://doi.org/10.1371/journal.pone.0248632.s007

(TIF)

S8 Fig. 31P metabolite ratios extracted from 3 Tesla spectra TR = 25 seconds.

The MP signals were fitted using the four spectral models: Blue: Gaussian with a fixed linewidth of 100 Hz; orange: Gaussian with a freely adjustable linewidth; grey: Lorentzian with a fixed linewidth of 200 Hz; yellow: Lorentzian with a freely adjustable linewidth.

https://doi.org/10.1371/journal.pone.0248632.s008

(TIF)

Acknowledgments

The authors thank Ioline Henter (NIMH) for invaluable editorial assistance.

References

  1. 1. Bottomley PA, Drayer BP, Smith LS. Chronic adult cerebral infarction studied by phosphorus NMR spectroscopy. Radiology 1986;160:763–766. pmid:3737915
  2. 2. Sappey-Marinier D, Calabrese G, Fein G, Hugg JW, Biggins C, Weiner MW. Effect of photic stimulation on human visual cortex lactate and phosphates using 1H and 31P magnetic resonance spectroscopy. J Cereb Blood Flow Metab. 1992;12:584–592. pmid:1618937
  3. 3. Levine SR, Helpern JA, Welch KM, Vande Linde AM, Sawaya KL, Brown EE, et al. Human focal cerebral ischemia: evaluation of brain pH and energy metabolism with P-31 NMR spectroscopy. Radiology 1992;185:537–44. pmid:1410369
  4. 4. Albers MJ, Krieger MD, Gonzalez-Gomez I, Gilles FH, McComb JG, Nelson MD Jr, et al. Proton-decoupled 31P MRS in untreated pediatric brain tumors. Magn Reson Med. 2005;53:22–29. pmid:15690498
  5. 5. Krikken E, van der Kemp WJM, van Diest PJ, van Dalen T, van Laarhoven HWM, Luijten PR, et al. Early detection of changes in phospholipid metabolism during neoadjuvant chemotherapy in breast cancer patients using phosphorus magnetic resonance spectroscopy at 7T. NMR Biomed. 2019;32:e4086. pmid:30924571
  6. 6. Gyulai L, Bolinger L, Leigh JS Jr, Barlow C, Chance B. Phosphorylethanolamine—the major constituent of the phosphomonoester peak observed by 31P-NMR on developing dog brain. FEBS Lett. 1984;1781:137–42. pmid:6500057
  7. 7. Kilby PM, Bolas NM, Radda GK. 31P-NMR study of brain phospholipid structures in vivo. Biochim Biophys Acta. 1991;1085:257–264. pmid:1654113
  8. 8. Qiao H, Zhang X, Zhu XH, Du F, Chen W. In vivo 31P MRS of human brain at high/ultrahigh fields: a quantitative comparison of NMR detection sensitivity and spectral resolution between 4 T and 7 T. Magn Reson Imaging. 2006;24:1281–1286. pmid:17145398
  9. 9. Lu M, Chen W, Zhu XH. Field dependence study of in vivo brain 31P MRS up to 16.4 T. NMR Biomed. 2014;27:1135–1141. pmid:25070004
  10. 10. Bluml S, Seymour KJ, Ross BD. Developmental changes in choline- and ethanolamine-containing compounds measured with proton-decoupled 31P MRS in in vivo human brain. Magn Reson Med. 1999;42:643–654. pmid:10502752
  11. 11. Novak J, Wilson M, MacPherson L, Arvanitis TN, Davies NP, Peet AC. Eur J Radiol. Clinical protocols for 31P MRS of the brain and their use in evaluating optic pathway gliomas in children. 2014;83:e106–e112.
  12. 12. Culeddu N, Bosco M, Toffanin R, Pollesello P. 31P NMR analysis of phospholipids in crude extracts from different sources: improved efficiency of the solvent system. Magn Reson Chem. 1998;36:907–912.
  13. 13. Muller CP, Reichel M, Muhle C, Rhein C, Gulbins E, Kornhuber J. Brain membrane lipids in major depression and anxiety disorders. Biochim Biophys Acta. 2015;1851:1052–1065. pmid:25542508
  14. 14. Zhang Y, Shen J. Effects of noise and linewidth on in vivo analysis of glutamate at 3 T. J Magn Reson. 2020;314:106732. pmid:32361510
  15. 15. Cadoux-Hudson TA, Blackledge MJ, Rajagopalan B, Taylor DJ, Radda GK. Human primary brain tumour metabolism in vivo: a phosphorus magnetic resonance spectroscopy study. Br J Cancer. 1989;60:430–436. pmid:2551360
  16. 16. Mathur R, Cox IJ, Oatridge A, Shephard DT, Shaw RJ, Taylor-Robinson SD. Cerebral bioenergetics in stable chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 1999;160:1994–1999. pmid:10588619
  17. 17. Buchli R, Martin E, Boesiger P, Rumpel H. Developmental changes of phosphorus metabolite concentrations in the human brain: a 31P magnetic resonance spectroscopy study in vivo. Pediatr Res. 1994;35:431–435. pmid:8047379
  18. 18. Wang X, Lee JH. IRIS-HSVD algorithm for automatic quantitation of in vivo 31P MRS. J Magn Reson. 2009;196:23–32. pmid:18926748
  19. 19. Andrade CS, Otaduy MC, Valente KD, Maia DF, Park EJ, Valerio RM, et al. Phosphorus magnetic resonance spectroscopy in malformations of cortical development. Epilepsia. 2011;52:2276–2284. pmid:21973076
  20. 20. Ha DH, Choi S, Oh JY, Yoon SK, Kang MJ, Kim KU. Application of 31P MR spectroscopy to the brain tumors. Korean J Radiol. 2013;14:477–486. pmid:23690717
  21. 21. Vanhamme L, van den Boogaart A, Van Huffel S. Improved method for accurate and efficient quantification of MRS data with use of prior knowledge. J Magn Reson. 1997;129:35–43. pmid:9405214
  22. 22. Lu M, Zhu XH, Zhang Y, Chen W. Intracellular redox state revealed by in vivo 31P MRS measurement of NAD+ and NADH contents in brains. Magn Reson Med. 2014;71:1959–1972. pmid:23843330
  23. 23. Smith SA, Levante TO, Meier BH, Ernst RR. Computer simulations in magnetic resonance. An object oriented programming approach. J Magn Reson. 1994;106A:75–105.
  24. 24. Deelchand DK, Nguyen TM, Zhu XH, Mochel F, Henry PG. Quantification of in vivo 31P NMR brain spectra using LCModel. NMR Biomed. 2015;28:633–641 pmid:25871439
  25. 25. Jope RS, Jenden DJ. Choline and phospholipid metabolism and the synthesis of acetylcholine in rat brain. J Neurosci Res. 1979;4:69–82. pmid:423317
  26. 26. Nirogi R, Mudigonda K, Kandikere V, Ponnamaneni R. Quantification of acetylcholine, an essential neurotransmitter, in brain microdialysis samples by liquid chromatography mass spectrometry. Biomed Chromatogr. 2010;24:39–48. pmid:19877295
  27. 27. Klein J, Gonzalez R, Koppen A, Loffelholz K. Free choline and choline metabolites in rat brain and body fluids: sensitive determination and implications for choline supply to the brain. Neurochem Int. 1993;22:293–300. pmid:8443570
  28. 28. Lakher M, Wurtman RJ, Blusztajn J, Holbrook P, Maire JC, Mauron C, et al. Brain phosphatidylcholine pools as possible sources of free choline for acetylcholine synthesis. In: Borchardt RT, Creveling CR, Ueland PMeditors. Biological methylation and drug design: experimental and clinical role of sadenosylmethionine. Humana Press Clifton; 1986. pmid:3025910
  29. 29. Veech RL, Lawson JWR, Cornell NW, Krebs HA. Cytosolic phosphorylation potential. J Biol Chem. 1979;254:6538–6547. pmid:36399
  30. 30. Pfeuffer J, Tkac I, Provencher SW, Gruetter R. Toward an in vivo neurochemical profile: quantification of 18 metabolites in short-echo-time 1H NMR spectra of the rat brain. J Magn Reson. 1999;141:104–120. pmid:10527748
  31. 31. Xu S, Yang J, Li CQ, Zhu W, Shen J. Metabolic alterations in focally activated primary somatosensory cortex of alpha-chloralose-anesthetized rats measured by 1H MRS at 11.7 T. Neuroimage 2005;28:401–409. pmid:16182571
  32. 32. Hong ST, Balla DZ, Pohmann R. Determination of regional variations and reproducibility in in vivo 1H NMR spectroscopy of the rat brain at 16.4 T. Magn Reson Med. 2011;66:11–17. pmid:21604299
  33. 33. Ren J, Sherry AD, Malloy CR. 31P-MRS of healthy human brain: ATP synthesis, metabolite concentrations, pH, and T1 relaxation times. NMR Biomed. 2015;28:1455–1462. pmid:26404723
  34. 34. Govindaraju V, Young K, Maudsley AA. Proton NMR chemical shifts and coupling constants for brain metabolites. NMR Biomed. 2000;13:129–153. pmid:10861994
  35. 35. Kreis R. Quantitative localized 1H MR spectroscopy for clinical use. J of Prog Nucl Magn Reson Spectroscopy. 1997;31:155–195.
  36. 36. Hoang TQ, Bluml S, Dubowitz DJ, Moats R, Kopyov O, Jacques D, et al. Quantitative proton-decoupled 31P MRS and 1H MRS in the evaluation of Huntington’s and Parkinson’s diseases. Neurology 1998;50:1033–1040. pmid:9566391
  37. 37. Bluml S, Tan J, Harris K, Adatia N, Karme A, Sproull T, et al. Brain quantitative proton-decoupled 31P MRS of the schizophrenic brain in vivo. J Comp Assist Tomo. 1999;23:272–275. pmid:10096336
  38. 38. Keshavan MS, Stanley JA, Montrose DM, Minshew NJ, Pettegrew JW. Prefrontal membrane phospholipid metabolism of child and adolescent offspring at risk for schizophrenia or schizoaffective disorder: an in vivo 31P MRS study. Mol Psychiatry. 2003;8:316–323. pmid:12660804
  39. 39. Shi XF, Carlson PJ, Sung YH, Fiedler KK, Forrest LN, Hellem TL, et al. Decreased brain PME/PDE ratio in bipolar disorder: a preliminary 31P magnetic resonance spectroscopy study. Bipolar Disord. 2015;17:743–752. pmid:26477793
  40. 40. Pettegrew JW, Panchalingam K, Klunk WE, McClure RJ, Muenz LR. Alterations of cerebral metabolism in probable Alzheimer’s disease: a preliminary study. Neurobiol Aging. 1994;15:117–132. pmid:8159258
  41. 41. Peeters TH, van Uden MJ, Rijpma A, Scheenen TWJ, Heerschap A. 3D 31P MR spectroscopic imaging of the human brain at 3 T with a 31P receive array: An assessment of 1H decoupling, T1 relaxation times, 1H-31P nuclear Overhauser effects and NAD. NMR Biomed. 2019;13:e4169.
  42. 42. An L, Li S, Shen J. Simultaneous determination of metabolite concentrations, T1 and T2 relaxation times. Magn Reson Med. 2017;78:2072–2081. pmid:28164364
  43. 43. Klomp DWJ, Wijnen JP, Scheenen TWJ, Heerschap A. Efficient 1H to 31P polarization transfer on a clinical 3T MR system. Magn Reson Med. 2008;60:1298–1305. pmid:19030163