G.
Diserens
a,
M.
Vermathen
b,
C.
Precht
c,
N. T.
Broskey
d,
C.
Boesch
a,
F.
Amati
ad,
J.-F.
Dufour
e and
P.
Vermathen
*a
aDepts. Clinical Research and Radiology, University of Bern, Bern, Switzerland. E-mail: peter.vermathen@insel.ch
bDept. of Chemistry & Biochemistry, University of Bern, Bern, Switzerland
cDept. of Clinical Veterinary Medicine, Clinical Radiology, Vetsuisse-Faculty, University of Bern, Bern, Switzerland
dDept. of Physiology, University of Lausanne, Lausanne, Switzerland
eHepatology, Dept. Clinical Research and University Clinic for Visceral Surgery and Medicine, Inselspital, University of Bern, Bern, Switzerland
First published on 6th October 2014
High Resolution Magic Angle Spinning (HR-MAS) NMR allows metabolic characterization of biopsies. HR-MAS spectra from tissues of most organs show strong lipid contributions that are overlapping metabolite regions, which hamper metabolite estimation. Metabolite quantification and analysis would benefit from a separation of lipids and small metabolites. Generally, a relaxation filter is used to reduce lipid contributions. However, the strong relaxation filter required to eliminate most of the lipids also reduces the signals for small metabolites. The aim of our study was therefore to investigate different diffusion editing techniques in order to employ diffusion differences for separating lipid and small metabolite contributions in the spectra from different organs for unbiased metabonomic analysis. Thus, 1D and 2D diffusion measurements were performed, and pure lipid spectra that were obtained at strong diffusion weighting (DW) were subtracted from those obtained at low DW, which include both small metabolites and lipids. This subtraction yielded almost lipid free small metabolite spectra from muscle tissue. Further improved separation was obtained by combining a 1D diffusion sequence with a T2-filter, with the subtraction method eliminating residual lipids from the spectra. Similar results obtained for biopsies of different organs suggest that this method is applicable in various tissue types. The elimination of lipids from HR-MAS spectra and the resulting less biased assessment of small metabolites have potential to remove ambiguities in the interpretation of metabonomic results. This is demonstrated in a reproducibility study on biopsies from human muscle.
However, HR-MAS spectra from tissues of most organs show strong contributions from broad lipid signals. Although these may convey important information themselves,2 they are overlapping metabolite regions. This could hamper metabolite estimation and may lead to misinterpretation. Metabolite quantification and metabonomic analysis would therefore benefit from a separation of lipid resonances from those of small metabolites.
Spectral fitting may be applied for separating lipids from small metabolites.3 However, solely fitting the lipid contribution in spectra that contain also small metabolites is prone to artificial inclusion of metabolite contributions into the lipid fit. Fitting small metabolites and lipids simultaneously has many advantages, but it also introduces bias as to which metabolites to include and also on the included metabolites themselves, which is avoided in metabonomic analysis.
The Carr–Purcell–Meiboom–Gill (CPMG) sequence is commonly used as a T2-filter to reduce broad components and therefore also fat contributions from spectra, taking advantage of the short T2 relaxation times of lipids.2,4 However, long echo times would be required to eliminate most of the lipids,5 with consequently (A) reduced signal to noise also for small metabolites, (B) a T2 dependency of peak intensities impeding metabolite quantification, (C) heating of the sample, (D) J-evolution of spins, and (E) losing all information on lipids.
Besides shorter T2, lipids differ also from small metabolites in terms of much slower diffusion. Fast diffusing small metabolites may therefore be separated from slow diffusing lipids by NMR diffusion measurements. It has previously been shown that subtracting a strong diffusion weighted (DW) spectrum from a low DW spectrum can lead to a strong reduction or even elimination of spectral contributions from lipids. This was applied in high resolution liquid state NMR measurements of proteins,6 carbohydrates,7 as well as in rats8 and human blood plasma,9 and in in vivo measurements of the brain.10 HR-MAS NMR diffusion editing was also applied in heterogeneous mixtures containing liquid-, gel- and solid-like phases,11 as well as in different seeds.12
Moreover, an almost complete separation of small molecules and lipids may be obtained by a combination of both physical differences, i.e. combining diffusion gradient editing and the T2-filter.13 The two effects have been combined in oil–water samples14 for separate measurements of individual T2 attenuations of the two components, and conversely in vegetables for improved separated determination of diffusion constants.15 In liver biopsies, the two effects have been combined for distinguishing lipid and low molecular weight resonances for both improved relaxation and diffusion constant estimations, as this can provide information of metabolite changes inside the tissue.5 Liu et al. showed the possibility to separate the high molecular weight molecules by diffusion editing, the low molecular weight molecules by relaxation editing and applied a combined editing method (combined diffusion and T2-filter sequence) in human blood plasma for better lipoprotein analysis.9
However, to the best of our knowledge, diffusion gradient editing alone or its combination with a T2-filter has not yet been applied in HR-MAS spectra of biopsies for obtaining lipid-free small metabolite spectra.
The aim of our study was therefore to investigate different diffusion editing techniques in order to employ diffusion differences for separating lipid from small metabolite contributions in the spectra of various tissues for unbiased metabonomic analysis. Thus NMR-diffusion measurements applying different gradient strengths were performed on biopsies of human muscle. Spectra obtained at strong DW solely consisting of lipid signals were subtracted from those obtained at low DW including both small metabolites and lipids. In addition, the study aimed at simultaneously utilizing diffusion and T2 differences15 between small metabolites and lipids for an efficient separation. Therefore additional measurements were performed on biopsies from different organs and animals combining both effects.
In order to probe the impact of this editing method, it was applied in a reproducibility study of muscle spectra including chemometric analyses. The metabolite variability of two specimens of the same muscle was investigated before and after removing the lipid peaks. The secondary aim was to determine to what extent the metabolite profile of a single muscle specimen can be considered representative for the tissue.
Additional CPMG measurements were performed on pork biopsies for parameter optimization, modifying the TE, the number of loops and the interpulse delay time for each loop. The four following combinations were used: (1) a TE of 76.8 ms with 192 loops of 0.4 ms; (2) a TE of 152 ms with 380 loops of 0.4 ms; (3) a TE of 152 ms with 38 loops of 4 ms; (4) a TE of 760 ms with 380 individual loops of 2 ms each.
The 2D-diffusion measurements were processed using Bruker Dynamics Center 2.0.8. The peaks were fitted using a bi-exponential model, taking all 64 points in F2 into account. In addition the individual spectra of the 64 different DWs were investigated. Spectra at stronger gradients, which exhibited purely lipid signals, were added, yielding a pure lipid spectrum with a good SNR. Depending on the number of added spectra, the weightings of different lipid components vary, if the diffusion constants are different. It turned out that for Δ = 50 ms ideally the last 20 spectra could be added to yield a pure lipid spectrum [see Results section 3.2.2]. This summed spectrum was then subtracted from the summed spectrum at low gradient strength (the first 20 individual spectra of the series were selected for the best SNR).
For chemometric analysis of the muscle biopsies for the reproducibility study, a total of 63 buckets (between 0.95 and 8.62 ppm) were selected, with a variable size according to the peak width. Spectral regions comprising only noise and pure lipid regions were excluded from all analyses to facilitate a comparison of the small metabolite content from the original and lipid-subtracted CPMG spectra (see ESI†, Fig. S1). In addition, some spectra were contaminated by lidocaine, the anaesthetic agent. Five spectral regions with potential contributions from lidocaine resonances (peaks at 1.37, 2.2, 3.37, 4.45, 7.21, 7.26 ppm) were also excluded from the analysis. The buckets were normalized by probabilistic quotient normalization (PQN)25 and scaled with Pareto scaling, which keeps the influence of noise low, while decreasing contributions from intense peaks and increasing the weight of intermediate peaks.26 MatLab (R2011b, The MathsWorks Inc.) was used for statistical analysis.
The whole muscle analysis (bucketing, normalization, scaling and statistical analysis) was first done on the pure CPMG spectra, and then repeated with the lipid-subtracted CPMG spectra [see 2.2.2.2], to evaluate the influence of the lipid removal on the reproducibility metabonomic study.
Different numbers of spectra with strong and weak DW, respectively, were added from the pseudo-2D data, (A) to obtain a sufficient SNR, (B) to determine and potentially adjust for slightly different spectral contributions of lipid moieties due to different diffusion constants. It turned out that the number of added spectra had only a minor influence on the spectral lipid contributions.
An example for the subtraction is shown in Fig. 2 (see Table 1 for lipid assignments). The lipid spectrum subtraction removes the shoulder from the creatine peak at 3.0 ppm, and effectively suppresses the lipid broad component underneath the choline, taurine and glycerophosphocholine (GPC) peaks (Fig. 2, inset). The strong overlap of the lactate peak at 1.3 ppm or leucine, isoleucine and valine at ∼1.0 ppm with lipids is similarly strongly reduced. The lipid peaks at 0.9 and 2.0 ppm are well suppressed, but small lipid residues remain, which are – besides residual artefacts – probably partly due to mobile lipids.22 Nevertheless, the residues still impede metabolite quantification or inclusion in chemometric methods. Therefore, for a chemometric analysis of small metabolites, the pure lipid regions (0.9 and 1.2 ppm) need to be omitted. However, the procedure yields clearly reduced lipid contributions for several metabolites like lactate at 1.33 ppm, choline, GPC, taurine, or creatine and may thus lead to more reliable metabolite estimations, while at the same time maintaining the maximum signal strength without relaxation damping, which also hampers metabolite quantification. Furthermore, the pure lipid spectrum from the 2D diffusion data can be analyzed separately, besides the possibility to perform a regular diffusion analysis.
![]() | ||
Fig. 2 Example of human muscle diffusion spectra between 0.5 and 4.5 ppm. Blue: the summed spectrum with low DW showing small metabolites and strong lipid resonances; red: the summed spectrum with strong DW showing lipids only; green: the difference spectrum. Inset: the zoomed region. Lipid resonances are assigned according to Table 1. |
# | Lipid part | ppm |
---|---|---|
1 | ω-CH3 | 0.91 |
2 | –(CH2)n | 1.3 |
3 | β-CH2 | 1.59 |
4 | –CH2–CH![]() |
2.03 |
5 | α-CH2 | 2.26 |
6 | ![]() ![]() |
2.78 |
7 | Tentatively assigned to CH2 lysyl protein | 3.03 |
8 | (–N+(CH3)3)phosphatidylcholine | 3.25 |
9 | CH2-glyceryl | 4.15 |
10 | Cho–a–CH2; CH2-glyceryl | 4.35 |
11 | –CH![]() |
5.32 |
T2 differences between lipid resonances5 result in different signal intensities in CPMG spectra compared to those in standard 2D diffusion measurements. Therefore, a subtraction of the entire lipid spectrum after scaling with a single factor would leave remaining lipid signals, depending on their respective T2 differences. Instead, a specific intensity scaling for each resonance in the high DW spectra was chosen to minimize the residue (but retaining frequency and line width). This procedure allows for an effective lipid removal (Fig. 3) and is similar to fitting the lipid resonances, however, without the problem to include resonances from small metabolites into the fit. The figure demonstrates that the lipid contributions underlying some of the small metabolites are removed to a large extent. However, as for the subtraction of the high from the low DW spectra, residual artefacts remain at 0.9 and 1.2 ppm.
Spectrum in Fig. 4 | TE [ms] | Number of loops | Interpulse delay [ms] | Comment on spectra |
---|---|---|---|---|
1 – red | 76.8 | 192 | 0.4 | No visual J-modulation, good SNR |
2 – green | 152 | 380 | 0.4 | No visual J-modulation, slightly decreased intensities, stronger decreased lipids |
3 – cyan | 152 | 38 | 4 | Strong J-modulation |
4 – blue | 760 | 380 | 2 | Lipids entirely suppressed, but J-modulation and strong reduction of metabolite intensities |
![]() | ||
Fig. 4 Example of pork muscle CPMG spectra acquired with different parameters (see Table 2), all spectra are scaled to their creatine peak. Top: spectra between 0.8 and 3.1 ppm. Bottom: spectra between 2.6 and 4.7 ppm. |
Spectrum 1 (TE = 76.8 ms) and spectrum 2 (TE = 152 ms) appear quite similar despite the doubled echo time, with only a small decrease of the peak intensity for a longer TE. Both spectra were acquired with 0.4 ms interpulse delay time and do not demonstrate substantial J-modulation.
Spectrum 3 was acquired at the same TE as spectrum 2, but with a 10-time longer interpulse delay time (4 ms). The strong J-evolution leads to complicated spectral patterns and reduced signal intensities (visually e.g. at 1.5, 2.7 or 3.15 ppm) hampering metabolite quantification and inclusion in metabonomic analysis. The much greater difference of spectrum 3 compared to spectra 1 and 2 points out the strong impact of the interpulse delay time for each loop (see also ref. 27).
Increasing the echo time to 760 ms reduced considerably the lipid content and yielded an almost fat free spectrum (spectrum 4). However, metabolite resonances were also clearly reduced at this long TE, as seen e.g. at the cysteine multiplet at 3.1 ppm which disappears entirely in the noise, as well as clear intensity decreases of resonances between 3.5 and 3.8 ppm or at 2.7 ppm. Thus, increasing only the TE appears not well suited for complete removal of lipid contributions, because of a low SNR.
In conclusion, the results of the interpulse delay time and TE variation in the original CPMG sequence suggest that the combination of a CPMG sequence at moderate TE but short interpulse spacing with a DW module and subsequent lipid subtraction appears most promising for complete lipid removal while maintaining signals from small metabolites for analysis.
![]() | ||
Fig. 5 1D diffusion-CPMG spectra of a pork muscle (blue: low DW; red: strong DW; green: difference). For the lipid assignment, see Table 1. |
Similar results were obtained for spectra from sheep brain biopsies, though with much less prominent lipids (Fig. 6) and in the spectra of human liver with very strong lipid peaks (Fig. 7), thus, suggesting that this method is applicable in various tissue types.
![]() | ||
Fig. 6 1D diffusion–CPMG spectra of a sheep brain (blue: low DW; red: strong DW; green: difference). |
![]() | ||
Fig. 7 1D diffusion–CPMG spectra of a human liver sample (blue: low DW; red: strong DW; green: difference). For lipid assignment, see Table 1. Inset: zoom on the 1.15–1.45 ppm region – the subtracted spectrum is increased 6 times for a better visualization of the lactate peak. |
For the liver spectra, a residual lipid artefact can be seen at 1.2 ppm after subtraction (Fig. 7), probably due to the very high amount of lipid with the –(CH2)n lipid peak being 20 times higher than the creatine peak at 3.05 ppm. However, the lactate peak is clearly detectable in the subtracted spectrum and appears unaffected from lipids (see the inset). Lactate could thus be included in the metabonomic analysis, while it was almost hidden by lipids in the original spectrum (Fig. 7 – inset).
The 2D diffusion experiments permitted a flexible adjustment of strong and low diffusion weighted spectra for the best balance between metabolite attenuation and SNR. However, after establishing the required DW for optimal lipid removal, 1D-diffusion measurements may be preferred in the future for a better SNR and shorter scan times. The combined diffusion-CPMG sequence yielded high quality small metabolite and lipid spectra with no obvious spectral impairment compared to a regular CPMG sequence and benefits from amplification of simultaneous T2 and diffusion filters.
Principal Component Analysis (PCA) of the entire spectrum (Fig. 8) without lipid removal demonstrated clustering for most of the same needle biopsies, both on PC1/PC2 and PC2/PC3, indicating spectral similarities between them. However, because of considerable differences among the same needle biopsies in few subjects there are no significant group differences for PC1 (p = 0.15), PC2 (p = 0.23) and PC3 (p = 0.05).
![]() | ||
Fig. 8 Same-needle biopsy analysis without lipid subtraction (left) and with lipid subtraction (right). Up: PCA, down: loading plots. Red ellipse: 95% confidence interval. |
The PCA results reveal one biopsy (S4) outside of the 95% confidence interval (Fig. 8a), and the loading plot (Fig. 8b) indicates the bucket at 1.3 ppm to be responsible for this separation, which is most likely due to the –(CH2)n lipid group besides lactate. The loading plot from the PCA using non-subtracted spectra is dominated by lipid peaks (despite the exclusion of pure lipid regions from the analysis) and therefore reveals the possible influence of the lipid content on chemometric analyses for metabolites that are biased by the underlying lipid.
The analysis was repeated after lipid subtraction of the CPMG spectra in order to reduce the lipid contribution in the principal component analysis separation. The PCA of the lipid-removed spectra (Fig. 8c) shows better clustering of same-needle biopsies (significant for PC2, p = 0.03), as well as a more homogeneous loading plot (Fig. 8d). Thus, the PCA separation is due to a combination of different metabolites and not dominated by lipids as in the non-subtracted spectra. The two biopsies from sample S4 are much closer to each other after the lipid removal than before, indicating that the two biopsies have a similar metabolite but different lipid contents. The results suggest that lipid removal in this reproducibility study leads to clustering and group separation that is based predominantly on small metabolite differences instead of lipid contributions. If on the other hand lipid components are the focus of interest, the pure lipid spectra used for subtraction can be analysed accordingly in the absence of small metabolites.
A secondary aim of this reproducibility study was to investigate how far the metabolite profile of a single biopsy can be considered representative for the tissue under investigation. Same needle biopsies harvested from the same muscle demonstrated not only similarities but also quite surprisingly strong differences in the metabolite content as obtained from HR-MAS. This variability needs to be taken into account not only in longitudinal HR-MAS studies on muscle tissue but also for any other study that relies on the assumption that a biopsy sample represents a specific muscle type. It should be noted that the quite high variability also among the same needle biopsies was not due to the spectral quality, but may be due to the muscle fiber type composition.28 It has been suggested to use large biopsy samples to reduce intra-individual variation in the percentages of fast and slow fibers,29 although other studies have found only low variation in fiber type distribution within single biopsies.30 As discussed, a high variability in the lipid region may occur and is probably due to varying inclusions of small residual contributions from extramyocellular lipids despite careful prevention during biopsy harvesting.
A reproducibility study of the same-needle muscle biopsy measurements demonstrated the applicability of the method and its potential value, by shifting the separation between groups from lipid contributions to different other metabolites. The elimination of lipids from the HR-MAS spectra and the resulting less biased assessment of small metabolites thus have potential to remove ambiguities in the interpretation of metabonomic results. Moreover, the diffusion scans permit a separate lipid analysis and a regular diffusion analysis without expanding the measurement time.
Footnote |
† Electronic supplementary information (ESI) available: Table of spectral regions that were excluded from the analysis of the same-needle biopsy; the corresponding muscle spectrum indicating the excluded regions. See DOI: 10.1039/c4an01663g |
This journal is © The Royal Society of Chemistry 2015 |