On the use of Cramér-Rao minimum variance bounds for the design of magnetic resonance spectroscopy experiments

Localized Magnetic Resonance Spectroscopy (MRS) is in widespread use for clinical brain research. Standard acquisition sequences to obtain one-dimensional spectra suffer from substantial overlap of spectral contributions from many metabolites. Therefore, specially tuned editing sequences or two-dimensional acquisition schemes are applied to extend the information content. Tuning specific acquisition parameters allows to make the sequences more efficient or more specific for certain target metabolites. Cramér-Rao bounds have been used in other fields for optimization of experiments and are now shown to be very useful as design criteria for localized MRS sequence optimization. The principle is illustrated for one- and two-dimensional MRS, in particular the 2D separation experiment, where the usual restriction to equidistant echo time spacings and equal acquisition times per echo time can be abolished. Particular emphasis is placed on optimizing experiments for quantification of GABA and glutamate. The basic principles are verified by Monte Carlo simulations and in vivo for repeated acquisitions of generalized two-dimensional separation brain spectra obtained from healthy subjects and expanded by bootstrapping for better definition of the quantification uncertainties.

Quantitative 1H-magnetic resonance spectroscopy of human brain: Influence of composition and parameterization of the basis set in linear combination model-fitting

Localized short-echo-time (1)H-MR spectra of human brain contain contributions of many low-molecular-weight metabolites and baseline contributions of macromolecules. Two approaches to model such spectra are compared and the data acquisition sequence, optimized for reproducibility, is presented. Modeling relies on prior knowledge constraints and linear combination of metabolite spectra. Investigated was what can be gained by basis parameterization, i.e., description of basis spectra as sums of parametric lineshapes. Effects of basis composition and addition of experimentally measured macromolecular baselines were investigated also. Both fitting methods yielded quantitatively similar values, model deviations, error estimates, and reproducibility in the evaluation of 64 spectra of human gray and white matter from 40 subjects. Major advantages of parameterized basis functions are the possibilities to evaluate fitting parameters separately, to treat subgroup spectra as independent moieties, and to incorporate deviations from straightforward metabolite models. It was found that most of the 22 basis metabolites used may provide meaningful data when comparing patient cohorts. In individual spectra, sums of closely related metabolites are often more meaningful. Inclusion of a macromolecular basis component leads to relatively small, but significantly different tissue content for most metabolites. It provides a means to quantitate baseline contributions that may contain crucial clinical information.

Brain metabolite composition during early human brain development as measured by quantitative in vivo 1H magnetic resonance spectroscopy

Biochemical maturation of the brain can be studied noninvasively by (1)H magnetic resonance spectroscopy (MRS) in human infants. Detailed time courses of cerebral tissue contents are known for the most abundant metabolites only, and whether or not premature birth affects biochemical maturation of the brain is disputed. Hence, the last trimester of gestation was observed in infants born prematurely, and their cerebral metabolite contents at birth and at expected term were compared with those of fullterm infants. Successful quantitative short-TE (1)H MRS was performed in three cerebral locations in 21 infants in 28 sessions (gestational age 32-43 weeks). The spectra were analyzed with linear combination model fitting, considerably extending the range of observable metabolites to include acetate, alanine, aspartate, cholines, creatines, gamma-aminobutyrate, glucose, glutamine, glutamate, glutathione, glycine, lactate, myo-inositol, macromolecular contributions, N-acetylaspartate, N-acetylaspartylglutamate, o-phosphoethanolamine, scyllo-inositol, taurine, and threonine. Significant effects of age and location were found for many metabolites, including the previously observed neuronal maturation reflected by an increase in N-acetylaspartate. Absolute brain metabolite content in premature infants at term was not considerably different from that in fullterm infants, indicating that prematurity did not affect biochemical brain maturation substantially in the studied population, which did not include infants of extremely low birthweight.

Magnetic resonance spectroscopy investigations of functionally defined language areas in schizophrenia patients with and without auditory hallucinations

Background: Cerebral dysfunction occurring in mental disorders can show metabolic disturbances which are limited to circumscribed brain areas. Auditory hallucinations have been shown to be related to defined cortical areas linked to specific language functions. Here, we investigated if the study of metabolic changes in auditory hallucinations requires a functional rather than an anatomical definition of their location and size to allow a reliable investigation by magnetic resonance spectroscopy (MRS). Methods: Schizophrenia patients with (AH; n = 12) and without hallucinations (NH; n = 8) and healthy controls (HC; n = 11) underwent a verbal fluency task in functional MRI (fMRI) to functionally define Broca's and Wernicke's areas. Left and right Heschl's gyri were defined anatomically. Results: The mean distances in native space between the fMRI-defined regions and a corresponding anatomically defined area were 12.4 ± 6.1 mm (range: 2.7–36.1 mm) for Broca's area and 16.8 ± 6.2 mm (range: 4.5–26.4 mm) for Wernicke's area, respectively. Hence, the spatial variance was of similar extent as the size of the investigated regions. Splitting the investigations into a single voxel examination in the frontal brain and a spectroscopic imaging part for the more homogeneous field areas led to good spectral quality for almost all spectra. In Broca's area, there was a significant group effect (p = 0.03) with lower levels of N-acetyl-aspartate (NAA) in NH compared to HC (p = 0.02). There were positive associations of NAA levels in the left Heschl's gyrus with total (p = 0.03) and negative (p = 0.006) PANSS scores. In Broca's area, there was a negative association of myo-inositol levels with total PANSS scores (p = 0.008). Conclusion: This study supports the neurodegenerative hypothesis of schizophrenia only in a frontal region whereas the results obtained from temporal regions are in contrast to the majority of previous studies. Future research should test the hypothesis raised by this study that a functional definition of language regions is needed if neurochemical imbalances are expected to be restricted to functional foci.

Proton magnetic resonance spectroscopy in multiple sclerosis.

Proton magnetic resonance spectroscopy ((1)H-MRS) provides tissue metabolic information in vivo. This article reviews the role of MRS-determined metabolic alterations in lesions, normal-appearing white matter, gray matter, and spinal cord in advancing our knowledge of pathologic changes in multiple sclerosis (MS). In addition, the role of MRS in objectively evaluating therapeutic efficacy is reviewed. This potential metabolic information makes MRS a unique tool to follow MS disease evolution, understand its pathogenesis, evaluate the disease severity, establish a prognosis, and objectively evaluate the efficacy of therapeutic interventions.

Neuronal and axonal degeneration in experimental spinal cord injury: in vivo proton magnetic resonance spectroscopy and histology.

Longitudinal in vivo proton magnetic resonance spectroscopy (1H-MRS) and immunohistochemistry were performed to investigate the tissue degeneration in traumatically injured rat spinal cord rostral and caudal to the lesion epicenter. On 1H-MRS significant decreases in N-acetyl aspartate (NAA) and total creatine (Cr) levels in the rostral, epicenter, and caudal segments were observed by 14 days, and levels remained depressed up to 56 days post-injury (PI). In contrast, the total choline (Cho) levels increased significantly in all three segments by 14 days PI, but recovered in the epicenter and caudal, but not the rostral region, at 56 days PI. Immunohistochemistry demonstrated neuronal cell death in the gray matter, and reactive astrocytes and axonal degeneration in the dorsal, lateral, and ventral white-matter columns. These results suggest delayed tissue degeneration in regions both rostrally and caudally from the epicenter in the injured spinal cord tissue. A rostral-caudal asymmetry in tissue recovery was seen both on MRI-observed hyperintense lesion volume and the Cho, but not NAA and Cr, levels at 56 days PI. These studies suggest that dynamic metabolic changes take place in regions away from the epicenter in injured spinal cord.

Increased levels of choline metabolites are an early marker of docetaxel treatment response in BRCA1-mutated mouse mammary tumors: an assessment by ex vivo proton magnetic resonance spectroscopy.

BACKGROUND Docetaxel is one of the most frequently used drugs to treat breast cancer. However, resistance or incomplete response to docetaxel is a major challenge. The aim of this study was to utilize MR metabolomics to identify potential biomarkers of docetaxel resistance in a mouse model for BRCA1-mutated breast cancer. METHODOLOGY High resolution magic angle spinning (HRMAS) (1)H MR spectroscopy was performed on tissue samples obtained from docetaxel-sensitive or -resistant BRCA1-mutated mammary tumors in mice. Measurements were performed on samples obtained before treatment and at 1-2, 3-5 and 6-7 days after a 25 mg/kg dose of docetaxel. The MR spectra were analyzed by multivariate analysis, followed by analysis of the signals of individual compounds by peak fitting and integration with normalization to the integral of the creatine signal and of all signals between 2.9 and 3.6 ppm. RESULTS The HRMAS spectra revealed significant metabolic differences between sensitive and resistant tissue samples. In particular choline metabolites were higher in resistant tumors by more than 50% with respect to creatine and by more than 30% with respect to all signals between 2.9 and 3.6 ppm. Shortly after treatment (1-2 days) the normalized choline metabolite levels were significantly increased by more than 30% in the sensitive group coinciding with the time of highest apoptotic activity induced by docetaxel. Thereafter, choline metabolites in these tumors returned towards pre-treatment levels. No change in choline compounds was observed in the resistant tumors over the whole time of investigation. CONCLUSIONS Relative tissue concentrations of choline compounds are higher in docetaxel resistant than in sensitive BRCA1-mutated mouse mammary tumors, but in the first days after docetaxel treatment only in the sensitive tumors an increase of these compounds is observed. Thus both pre- and post-treatment tissue levels of choline compounds have potential to predict response to docetaxel treatment.

A dual resonance, image -guided, volume localization technique for magnetic resonance spectroscopy

An interleaved, dual resonance, volume localization technique for $\sp1$H/$\sp{31}$P magnetic resonance spectroscopy has been designed, implemented on a 2 T imager/spectrometer, and verified with phantom studies.^ Localization techniques, including several single voxel techniques and spectroscopic imaging, were implemented, and studies were performed to compare the efficiency of each sequence of $\sp1$H/$\sp{31}$P spectral acquisitions. The sequence chosen was a hybrid of the stimulated echo single voxel technique and the spectroscopic imaging technique.^ Water suppression during the $\sp1$H spectral acquisitions was accomplished by the use of three narrow bandwidth RF saturation pulses in combination with three spoiler gradients. The spoiler gradient amplitudes were selected on the basis of a numerical solution of the Bloch equations. A post-acquisition water suppression algorithm was used to minimize any residual water signal.^ For interleaved $\sp1$H/$\sp{31}$P acquisitions, a dual resonance RF coil was constructed and interfaced to the existing RF detection system via a custom-designed dual resonance transcoupler and switching system. Programmable attenuators were incorporated to allow for changes in receiver and transmitter attenuation "on the fly".^ To provide the rapidly switched gradient fields required for the $\sp1$H/$\sp{31}$P acquisitions, an actively screened gradient coil system was designed and implemented. With this system, gradient field rise times on the order of 100 $\mu$s were obtained. These rapid switching times were necessary for minimizing intrasequence delays and for improving localization quality and water suppression efficiency.^ The interleaved $\sp1$H/$\sp{31}$P volume localization technique was tested using a two-compartment phantom. Analysis of the data showed that the spectral contamination was less than three percent. One-to-one spatial correspondence of the $\sp1$H and $\sp{31}$P spectra was verified and allowed for direct correlation of the spectral data with a standard magnetic resonance image. ^

Biochemical changes in the frontal lobe of HIV-infected individuals detected by magnetic resonance spectroscopy

We have developed a proton magnetic resonance spectroscopy method that selectively can sample cortical gray matter and adjacent white matter in the frontal lobe. We have used this approach to study a group of patients (n = 7) infected with HIV and clinical manifestations of the AIDS dementia complex (ADC), a group of patients (n = 8) infected with HIV without any indications of ADC, and seven controls. The patients without ADC had a statistically significant increase in the ratio of myo-inositol to creatine in white matter compared with normal controls. In contrast, the group of patients with ADC had almost normal levels of myo-inositol to creatine in both gray matter and white matter and showed a statistically significant decrease in the N-acetylaspartate to creatine ratio in gray matter compared with either the normal controls or the patients without ADC. Patterns of spectral abnormalities correlated with neuropsychological measures of frontal lobe dysfunction, suggesting that the evaluation of frontal lobe metabolism by magnetic resonance spectroscopy can play a role in the early detection of ADC, in determining its progression, and in assessing responses to therapeutic interventions.

A Cross-Polarization, Magic-Angle-Spinning, 13C-Nuclear-Magnetic-Resonance Study of Polysaccharides in Sugar Beet Cell Walls1

Solid-state nuclear magnetic resonance relaxation experiments were used to study the rigidity and spatial proximity of polymers in sugar beet (Beta vulgaris) cell walls. Proton T1ρ decay and cross-polarization patterns were consistent with the presence of rigid, crystalline cellulose microfibrils with a diameter of approximately 3 nm, mobile pectic galacturonans, and highly mobile arabinans. A direct-polarization, magic-angle-spinning spectrum recorded under conditions adapted to mobile polymers showed only the arabinans, which had a conformation similar to that of beet arabinans in solution. These cell walls contained very small amounts of hemicellulosic polymers such as xyloglucan, xylan, and mannan, and no arabinan or galacturonan fraction closely associated with cellulose microfibrils, as would be expected of hemicelluloses. Cellulose microfibrils in the beet cell walls were stable in the absence of any polysaccharide coating.

Transport, Compartmentation, and Metabolism of Homoserine in Higher Plant Cells : Carbon-13- and Phosphorus-31-Nuclear Magnetic Resonance Studies

The transport, compartmentation, and metabolism of homoserine was characterized in two strains of meristematic higher plant cells, the dicotyledonous sycamore (Acer pseudoplatanus) and the monocotyledonous weed Echinochloa colonum. Homoserine is an intermediate in the synthesis of the aspartate-derived amino acids methionine, threonine (Thr), and isoleucine. Using 13C-nuclear magnetic resonance, we showed that homoserine actively entered the cells via a high-affinity proton-symport carrier (Km approximately 50–60 μm) at the maximum rate of 8 ± 0.5 μmol h−1 g−1 cell wet weight, and in competition with serine or Thr. We could visualize the compartmentation of homoserine, and observed that it accumulated at a concentration 4 to 5 times higher in the cytoplasm than in the large vacuolar compartment. 31P-nuclear magnetic resonance permitted us to analyze the phosphorylation of homoserine. When sycamore cells were incubated with 100 μm homoserine, phosphohomoserine steadily accumulated in the cytoplasmic compartment over 24 h at the constant rate of 0.7 μmol h−1 g−1 cell wet weight, indicating that homoserine kinase was not inhibited in vivo by its product, phosphohomoserine. The rate of metabolism of phosphohomoserine was much lower (0.06 μmol h−1 g−1 cell wet weight) and essentially sustained Thr accumulation. Similarly, homoserine was actively incorporated by E. colonum cells. However, in contrast to what was seen in sycamore cells, large accumulations of Thr were observed, whereas the intracellular concentration of homoserine remained low, and phosphohomoserine did not accumulate. These differences with sycamore cells were attributed to the presence of a higher Thr synthase activity in this strain of monocot cells.

31P magnetic resonance spectroscopy of the Sherpa heart: a phosphocreatine/adenosine triphosphate signature of metabolic defense against hypobaric hypoxia.

Of all humans thus far studied, Sherpas are considered by many high-altitude biomedical scientists as most exquisitely adapted for life under continuous hypobaric hypoxia. However, little is known about how the heart is protected in hypoxia. Hypoxia defense mechanisms in the Sherpa heart were explored by in vivo, noninvasive 31P magnetic resonance spectroscopy. Six Sherpas were examined under two experimental conditions [normoxic (21% FiO2) and hypoxic (11% FiO2) and in two adaptational states--the acclimated state (on arrival at low-altitude study sites) and the deacclimating state (4 weeks of ongoing exposure to low altitude). Four lowland subjects were used for comparison. We found that the concentration ratios of phosphocreatine (PCr)/adenosine triphosphate (ATP) were maintained at steady-state normoxic values (0.96, SEM = 0.22) that were about half those found in normoxic lowlanders (1.76, SEM = 0.03) monitored the same way at the same time. These differences in heart energetic status between Sherpas and lowlanders compared under normoxic conditions remained highly significant (P < 0.02) even after 4 weeks of deacclimation at low altitudes. In Sherpas under acute hypoxia, the heart rate increased by 20 beats per min from resting values of about 70 beats per min, and the percent saturation of hemoglobin decreased to about 75%. However, these perturbations did not alter the PCr/ATP concentration ratios, which remained at about 50% of the values expected in healthy lowlanders. Because the creatine phosphokinase reaction functions close to equilibrium, these steady-state PCr/ATP ratios presumably coincided with about 3-fold higher free adenosine diphosphate (ADP) concentrations. Higher ADP concentrations (i.e., lower [PCr]/[ATP] ratios) were interpreted to correlate with the Km values for ADP-requiring kinases of glycolysis and to reflect elevated carbohydrate contributions to heart energy needs. This metabolic organization is postulated as advantageous in hypobaria because the ATP yield per O2 molecule is 25-60% higher with glucose than with free fatty acids (the usual fuels utilized in the human heart in postfasting conditions).