Integrated data acquisition and processing to determine metabolite contents, relaxation times, and macromolecule baseline in single examinations of individual subjects

Absolute quantitation of clinical (1)H-MR spectra is virtually always incomplete for single subjects because the separate determination of spectrum, baseline, and transverse and longitudinal relaxation times in single subjects is prohibitively long. Integrated Processing and Acquisition of Data (IPAD) based on a combined 2-dimensional experimental and fitting strategy is suggested to substantially improve the information content from a given measurement time. A series of localized saturation-recovery spectra was recorded and combined with 2-dimensional prior-knowledge fitting to simultaneously determine metabolite T(1) (from analysis of the saturation-recovery time course), metabolite T(2) (from lineshape analysis based on metabolite and water peak shapes), macromolecular baseline (based on T(1) differences and analysis of the saturation-recovery time course), and metabolite concentrations (using prior knowledge fitting and conventional procedures of absolute standardization). The procedure was tested on metabolite solutions and applied in 25 subjects (15-78 years old). Metabolite content was comparable to previously found values. Interindividual variation was larger than intraindividual variation in repeated spectra for metabolite content as well as for some relaxation times. Relaxation times were different for various metabolite groups. Parts of the interindividual variation could be explained by significant age dependence of relaxation times.

Visibility of vascular phenylalanine in dynamic uptake studies in humans using magnetic resonance spectroscopy

In general, vascular contributions to the in vivo magnetic resonance (MR) brain spectrum are too small to be relevant. In cerebral uptake studies, however, vascular contributions may constitute a major confounder. MR visibility of vascular Phe was investigated by recording localized spectra from fully oxygenated and well-mixed whole blood. Blood Phe levels determined by MR spectroscopy (MRS) and ion-exchange chromatography showed excellent correlation. In addition, effects of blood flow were shown to have a small effect on signal amplitude with the MRS methodology used. Hence, blood Phe is almost completely MR visible at 1.5 T, even though it is severely broadened at higher fields. Without appropriate correction, cerebral Phe influx in studies of brain Phe uptake in phenylketonuria patients or healthy subjects would appear to be faster and lead to higher levels. Similar effects are envisaged for studies of ethanol or glucose uptake across the blood-brain barrier.

Distribution of intramyocellular lipids in human calf muscles as determined by MR spectroscopic imaging

In this study the distribution of intramyocellular lipids (IMCL) in human calf muscles was determined by 1H-MR spectroscopic imaging (MRSI) measurements. An obstacle for MRSI measurements in the calf, including different muscles, is the inevitable inclusion of regions with high concentrations of extramyocellular lipids (EMCL). This can lead to signal bleeding and consequently to unpredictable overlaps of IMCL resonances with EMCL in voxels of interest. The results of this study show that signal bleeding from EMCL can be substantially reduced in voxels from calf muscles by the application of a lipid extrapolation (LE) procedure (Haupt et al., Magn Reson Med 1996;35:678). The spectra of all voxels located within muscle tissue were fitted, and the metabolite values were assigned to one of 10 different muscles based on image segmentation. Significant IMCL differences between some muscles were obtained, with high values in m. soleus and two to three times lower values in the tibialis anterior, tibialis posterior, and gastrocnemius muscles. In addition to gross differences between muscles, significant intersubject differences were observed in both IMCL content and distribution over different muscles. A significant correlation between fiber orientation (obtained from orientation-dependent dipolar coupling of creatine and taurine resonances) and IMCL content was found, indicating that IMCL content is directly correlated to biomechanical properties.

Reproducibility of cerebral phenylalanine levels in patients with phenylketonuria determined by 1H-MR spectroscopy

The reproducibility of metabolite content determined by MR spectroscopy (MRS) is usually at best a few percent for the prominent singlets. When studying low-concentration metabolites, like phenylalanine (Phe), where tissue content can be <100 micromol/kg, better reproducibility is paramount-particularly in view of using MRS results for potential individual treatment advice. An optimized, targeted spectroscopy method was established at 1.5T and reproducibility was established in 21 patients with phenylketonuria (PKU) where three spectra were recorded in each of three independent sessions, two of which were in immediate succession to minimize physiologic variation. Intersession variation was found to be only 7 micromol/kg Phe for back-to-back repetition of sessions, in close agreement with the variation of 16 micromol/kg observed for single spectra within a session. Analysis of variance proved the individuality of the blood/brain Phe ratio-though this ratio seems to be influenced by physiologic factors that are not stable in time. The excellent reproducibility was achieved through optimization of various factors, including signal-to-noise ratio, repositioning, and prescan calibrations, but also by enforcing as much prior information as possible (e.g., lineshape and phase from reference scans, constant prior-knowledge-locked baseline). While the application of maximum general prior knowledge is a general method to reduce fluctuations, one should remember that it may introduce systematic errors.

Evaluation of texture features in hepatic tissue characterization from non-enhanced CT images

Aim of this paper is to evaluate the diagnostic contribution of various types of texture features in discrimination of hepatic tissue in abdominal non-enhanced Computed Tomography (CT) images. Regions of Interest (ROIs) corresponding to the classes: normal liver, cyst, hemangioma, and hepatocellular carcinoma were drawn by an experienced radiologist. For each ROI, five distinct sets of texture features are extracted using First Order Statistics (FOS), Spatial Gray Level Dependence Matrix (SGLDM), Gray Level Difference Method (GLDM), Laws' Texture Energy Measures (TEM), and Fractal Dimension Measurements (FDM). In order to evaluate the ability of the texture features to discriminate the various types of hepatic tissue, each set of texture features, or its reduced version after genetic algorithm based feature selection, was fed to a feed-forward Neural Network (NN) classifier. For each NN, the area under Receiver Operating Characteristic (ROC) curves (Az) was calculated for all one-vs-all discriminations of hepatic tissue. Additionally, the total Az for the multi-class discrimination task was estimated. The results show that features derived from FOS perform better than other texture features (total Az: 0.802+/-0.083) in the discrimination of hepatic tissue.

Accounting for baseline differences and measurement error in the analysis of change over time

If change over time is compared in several groups, it is important to take into account baseline values so that the comparison is carried out under the same preconditions. As the observed baseline measurements are distorted by measurement error, it may not be sufficient to include them as covariate. By fitting a longitudinal mixed-effects model to all data including the baseline observations and subsequently calculating the expected change conditional on the underlying baseline value, a solution to this problem has been provided recently so that groups with the same baseline characteristics can be compared. In this article, we present an extended approach where a broader set of models can be used. Specifically, it is possible to include any desired set of interactions between the time variable and the other covariates, and also, time-dependent covariates can be included. Additionally, we extend the method to adjust for baseline measurement error of other time-varying covariates. We apply the methodology to data from the Swiss HIV Cohort Study to address the question if a joint infection with HIV-1 and hepatitis C virus leads to a slower increase of CD4 lymphocyte counts over time after the start of antiretroviral therapy.

Transition from nerve stimulator to sonographically guided axillary brachial plexus anesthesia in hand surgery: block quality and patient satisfaction during the transition period

OBJECTIVES Sonographic guidance for peripheral nerve anesthesia has proven increasingly successful in clinical practice; however, fears that a change to sonographically guided regional anesthesia may impair the block quality and operating room work flow persist in certain units. In this retrospective cohort study, block quality and patient satisfaction during the transition period from nerve stimulator to sonographic guidance for axillary brachial plexus anesthesia in a tertiary referral center were investigated. METHODS Anesthesia records of all patients who had elective surgery of the wrist or hand during the transition time (September 1, 2006-August 25, 2007) were reviewed for block success, placement time, anesthesiologist training level, local anesthetic volume, and requirement of additional analgesics. Postoperative records were reviewed, and patient satisfaction was assessed by telephone interviews in matched subgroups. RESULTS Of 415 blocks, 341 were sonographically guided, and 74 were nerve stimulator guided. Sonographically guided blocks were mostly performed by novices, whereas nerve stimulator-guided blocks were performed by advanced users (72.3% versus 14%; P < .001). Block performance times and success rates were similar in both groups. In sonographically guided blocks, significantly less local anesthetics were applied compared to nerve stimulator-guided blocks (mean ± SD, 36.1 ± 7.1 versus 43.9 ± 6.1 mL; P< .001), and less opioids were required (fentanyl, 66.1 ± 30 versus 90 ± 62 μg; P< .001). Interviewed patients reported significantly less procedure-related discomfort, pain, and prolonged procedure time when block placement was sonographically guided (2% versus 20%; P = .002). CONCLUSIONS Transition from nerve stimulator to sonographic guidance for axillary brachial plexus blocks did not change block performance times or success rates. Patient satisfaction was improved even during the early institutional transition period.

Proton diffusion tensor spectroscopy of metabolites in human muscle in vivo

PURPOSE To study the apparent diffusivity and its directionality for metabolites of skeletal muscle in humans in vivo by (1) H magnetic resonance spectroscopy. METHODS The diffusion tensors were determined on a 3 Tesla MR system using optimized acquisition and processing methods including an adapted STEAM sequence with orientation-dependent diffusion weighting, pulse-triggering with individually adapted delays, eddy-current correction schemes, median filtering, and simultaneous prior-knowledge fitting of all related spectra. RESULTS The average apparent diffusivities, as well as the fractional anisotropies of taurine (ADCav  = 0.74 × 10(-3) s/mm(2) , FA = 0.46), creatine (ADCav  = 0.41 × 10(-3)  s/mm(2) , FA = 0.33), trimethylammonium compounds (ADCav  = 0.48 × 10(-3)  s/mm(2) , FA = 0.34), carnosine (ADCav  = 0.46 × 10(-3)  s/mm(2) , FA = 0.47), and water (ADCav  = 1.5 × 10(-3)  s/mm(2) , FA = 0.36) were estimated. The diffusivities of most metabolites and water were significantly different from each other. Diffusion was found to be anisotropic and the diffusion tensors showed tensor correlation coefficients close to 1 and were hence found to be essentially coaligned. The magnitudes of apparent metabolite diffusivities were largely ordered according to molecular weight, with taurine as the smallest molecule diffusing fastest, both along and across the fiber direction. CONCLUSION Diffusivities, directional dependence of diffusion and fractional anisotropies of (1) H MRS-visible muscle metabolites were presented. It was shown that metabolites share diffusion directionality with water and have similar fractional anisotropies, hinting at similar diffusion barriers. Magn Reson Med, 2014. © 2014 Wiley Periodicals, Inc.

Characterization of the macromolecule baseline in localized (1)H-MR spectra of human brain

Short-echo-time magnetic resonance spectra of human brain contain broad contributions from macromolecules. As they are a priori of unknown shape and intensity, they pose a problem if one wants to quantitate the overlying spectral features from low-molecular-weight metabolites. On the other hand, the macromolecular contributions may provide relevant clinical information themselves, if properly evaluated. Several methods, based on T(1), T(2), or spectral shape, have previously been suggested to suppress or edit the macromolecule contributions. Here, a method is presented based on a series of saturation recovery scans and that allows for simultaneous recording of the macromolecular baseline and the fully relaxed metabolite spectrum. In comparison to an inversion recovery technique aimed at nulling signals from long-T(1) components, the saturation recovery method is less susceptible to T(1) differences inherent in signals from different metabolites or introduced by pathology. The saturation recovery method was used to quantitate the macromolecular baseline in white and/or gray matter locations of the human brain in 40 subjects. It was found that the content and composition of MR visible macromolecules depends on cerebral location, as well as the age of the investigated subject, while no gender dependence could be found.

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.

Transconvolution and the virtual positron emission tomograph--a new method for cross calibration in quantitative PET∕CT imaging

PURPOSE Positron emission tomography (PET)∕computed tomography (CT) measurements on small lesions are impaired by the partial volume effect, which is intrinsically tied to the point spread function of the actual imaging system, including the reconstruction algorithms. The variability resulting from different point spread functions hinders the assessment of quantitative measurements in clinical routine and especially degrades comparability within multicenter trials. To improve quantitative comparability there is a need for methods to match different PET∕CT systems through elimination of this systemic variability. Consequently, a new method was developed and tested that transforms the image of an object as produced by one tomograph to another image of the same object as it would have been seen by a different tomograph. The proposed new method, termed Transconvolution, compensates for differing imaging properties of different tomographs and particularly aims at quantitative comparability of PET∕CT in the context of multicenter trials. METHODS To solve the problem of image normalization, the theory of Transconvolution was mathematically established together with new methods to handle point spread functions of different PET∕CT systems. Knowing the point spread functions of two different imaging systems allows determining a Transconvolution function to convert one image into the other. This function is calculated by convolving one point spread function with the inverse of the other point spread function which, when adhering to certain boundary conditions such as the use of linear acquisition and image reconstruction methods, is a numerically accessible operation. For reliable measurement of such point spread functions characterizing different PET∕CT systems, a dedicated solid-state phantom incorporating (68)Ge∕(68)Ga filled spheres was developed. To iteratively determine and represent such point spread functions, exponential density functions in combination with a Gaussian distribution were introduced. Furthermore, simulation of a virtual PET system provided a standard imaging system with clearly defined properties to which the real PET systems were to be matched. A Hann window served as the modulation transfer function for the virtual PET. The Hann's apodization properties suppressed high spatial frequencies above a certain critical frequency, thereby fulfilling the above-mentioned boundary conditions. The determined point spread functions were subsequently used by the novel Transconvolution algorithm to match different PET∕CT systems onto the virtual PET system. Finally, the theoretically elaborated Transconvolution method was validated transforming phantom images acquired on two different PET systems to nearly identical data sets, as they would be imaged by the virtual PET system. RESULTS The proposed Transconvolution method matched different PET∕CT-systems for an improved and reproducible determination of a normalized activity concentration. The highest difference in measured activity concentration between the two different PET systems of 18.2% was found in spheres of 2 ml volume. Transconvolution reduced this difference down to 1.6%. In addition to reestablishing comparability the new method with its parameterization of point spread functions allowed a full characterization of imaging properties of the examined tomographs. CONCLUSIONS By matching different tomographs to a virtual standardized imaging system, Transconvolution opens a new comprehensive method for cross calibration in quantitative PET imaging. The use of a virtual PET system restores comparability between data sets from different PET systems by exerting a common, reproducible, and defined partial volume effect.