A multivariate approach for processing magnetization effects in triggered event-related functional magnetic resonance imaging time series

Triggered event-related functional magnetic resonance imaging requires sparse intervals of temporally resolved functional data acquisitions, whose initiation corresponds to the occurrence of an event, typically an epileptic spike in the electroencephalographic trace. However, conventional fMRI time series are greatly affected by non-steady-state magnetization effects, which obscure initial blood oxygen level-dependent (BOLD) signals. Here, conventional echo-planar imaging and a post-processing solution based on principal component analysis were employed to remove the dominant eigenimages of the time series, to filter out the global signal changes induced by magnetization decay and to recover BOLD signals starting with the first functional volume. This approach was compared with a physical solution using radiofrequency preparation, which nullifies magnetization effects. As an application of the method, the detectability of the initial transient BOLD response in the auditory cortex, which is elicited by the onset of acoustic scanner noise, was used to demonstrate that post-processing-based removal of magnetization effects allows to detect brain activity patterns identical with those obtained using the radiofrequency preparation. Using the auditory responses as an ideal experimental model of triggered brain activity, our results suggest that reducing the initial magnetization effects by removing a few principal components from fMRI data may be potentially useful in the analysis of triggered event-related echo-planar time series. The implications of this study are discussed with special caution to remaining technical limitations and the additional neurophysiological issues of the triggered acquisition.

Statistical evaluation of time-dependent metabolite concentrations: estimation of post-mortem intervals based on in situ 1H-MRS of the brain

Knowledge of the time interval from death (post-mortem interval, PMI) has an enormous legal, criminological and psychological impact. Aiming to find an objective method for the determination of PMIs in forensic medicine, 1H-MR spectroscopy (1H-MRS) was used in a sheep head model to follow changes in brain metabolite concentrations after death. Following the characterization of newly observed metabolites (Ith et al., Magn. Reson. Med. 2002; 5: 915-920), the full set of acquired spectra was analyzed statistically to provide a quantitative estimation of PMIs with their respective confidence limits. In a first step, analytical mathematical functions are proposed to describe the time courses of 10 metabolites in the decomposing brain up to 3 weeks post-mortem. Subsequently, the inverted functions are used to predict PMIs based on the measured metabolite concentrations. Individual PMIs calculated from five different metabolites are then pooled, being weighted by their inverse variances. The predicted PMIs from all individual examinations in the sheep model are compared with known true times. In addition, four human cases with forensically estimated PMIs are compared with predictions based on single in situ MRS measurements. Interpretation of the individual sheep examinations gave a good correlation up to 250 h post-mortem, demonstrating that the predicted PMIs are consistent with the data used to generate the model. Comparison of the estimated PMIs with the forensically determined PMIs in the four human cases shows an adequate correlation. Current PMI estimations based on forensic methods typically suffer from uncertainties in the order of days to weeks without mathematically defined confidence information. In turn, a single 1H-MRS measurement of brain tissue in situ results in PMIs with defined and favorable confidence intervals in the range of hours, thus offering a quantitative and objective method for the determination of PMIs.

Observation and identification of metabolites emerging during postmortem decomposition of brain tissue by means of in situ 1H-magnetic resonance spectroscopy

Postmortem decomposition of brain tissue was investigated by (1)H-magnetic resonance spectroscopy (MRS) in a sheep head model and selected human cases. Aiming at the eventual estimation of postmortem intervals in forensic medicine, this study focuses on the characterization and identification of newly observed metabolites. In situ single-voxel (1)H-MRS at 1.5 T was complemented by multidimensional homo- and heteronuclear high-resolution NMR spectroscopy of an extract of sheep brain tissue. The inclusion of spectra of model solutions in the program LC Model confirmed the assignments in situ. The first postmortem phase was characterized mainly by changes in the concentrations of metabolites usually observed in vivo and by the appearance of previously reported decay products. About 3 days postmortem, new metabolites, including free trimethylammonium, propionate, butyrate, and iso-butyrate, started to appear in situ. Since the observed metabolites and the time course is comparable in sheep and human brain tissue, the model system seems to be appropriate.

Quantitative (1)H magnetic resonance spectroscopy of myoglobin de- and reoxygenation in skeletal muscle: reproducibility and effects of location and disease

1H-magnetic resonance spectroscopy ((1)H-MRS) of deoxymyoglobin (DMb) provides a means to noninvasively monitor the oxygenation state of human skeletal muscle in work and disease. As shown in this work, it also offers the opportunity to measure the absolute tissue content of DMb, the basic oxygen consumption of resting muscle, and the reperfusion characteristics after release of a pressure cuff. The methodology to determine these tissue properties simultaneously at two positions along the calf is presented. The obtained values are in agreement with invasive determinations. The reproducibility of the (1)H-MRS measurements is established for healthy controls and patients with peripheral arterial disease (PAD). A location dependence in axial direction, as well as differences between controls and patients are demonstrated for all parameters. The reoxygenation time in particular is expected to provide a means to quantitatively monitor therapies aimed at improving muscular perfusion in these patients.

Prevalence and time course of acute mountain sickness in older children and adolescents after rapid ascent to 3450 meters

OBJECTIVE: Acute mountain sickness is a frequent and debilitating complication of high-altitude exposure, but there is little information on the prevalence and time course of acute mountain sickness in children and adolescents after rapid ascent by mechanical transportation to 3500 m, an altitude at which major tourist destinations are located throughout the world. METHODS: We performed serial assessments of acute mountain sickness (Lake Louise scores) in 48 healthy nonacclimatized children and adolescents (mean +/- SD age: 13.7 +/- 0.3 years; 20 girls and 28 boys), with no previous high-altitude experience, 6, 18, and 42 hours after arrival at the Jungfraujoch high-altitude research station (3450 m), which was reached through a 2.5-hour train ascent. RESULTS: We found that the overall prevalence of acute mountain sickness during the first 3 days at high altitude was 37.5%. Rates were similar for the 2 genders and decreased progressively during the stay (25% at 6 hours, 21% at 18 hours, and 8% at 42 hours). None of the subjects needed to be evacuated to lower altitude. Five subjects needed symptomatic treatment and responded well. CONCLUSION: After rapid ascent to high altitude, the prevalence of acute mountain sickness in children and adolescents was relatively low; the clinical manifestations were benign and resolved rapidly. These findings suggest that, for the majority of healthy nonacclimatized children and adolescents, travel to 3500 m is safe and pharmacologic prophylaxis for acute mountain sickness is not needed.

Leg ischemia: assessment with MR angiography and spectroscopy

PURPOSE: To prospectively determine reproducibility of magnetic resonance (MR) angiography and MR spectroscopy of deoxymyoglobin in assessment of collateral vessels and tissue perfusion in patients with critical limb ischemia (CLI) and to follow changes in patients undergoing intramuscular vascular endothelial growth factor (pVEGF)-C gene therapy, percutaneous transluminal angioplasty, supervised exercise training, or no therapy. MATERIALS AND METHODS: Study and gene therapy protocols were approved, and all patients gave written informed consent. To determine repeatability and reproducibility, seven patients underwent MR angiography and five underwent MR spectroscopy. The techniques were used to judge disease progress in 12 other patients with or without therapy: MR angiography to help determine change in visualization of collateral vessels and MR spectroscopy to help assess change in perfusion at proximal and distal calf levels. MR angiographic results were subjectively analyzed by three blinded readers. Intraobserver variability was expressed as 95% confidence interval (CI) (n=7); interobserver variability, as kappa statistic (n=15). Reexamination variability of MR spectroscopy was given as 95% CI for subsequent recovery times, and correlation with disease extent was calculated with Kendall taub rank correlation. Fisher-Yates test was used to correlate changes with pressure measurements and clinical course. RESULTS: Intraobserver and interobserver concordance was sensitive for detection of collateral vessels. Intraobserver agreement was 85.7% (95% CI: 42.1%, 99.6%). Interobserver agreement was high for small collateral vessels (kappa=0.74, P <.001) and fair for large collateral vessels (kappa=0.36, P=.002). MR spectroscopy was reproducible (95% CI: +/-26 seconds for proximal, +/-21 seconds for distal) and showed a correlation with disease extent (proximal calf, taub=0.84, P <.001; distal calf, taub=0.68, P=.04). Small collateral vessels increased over time (P=.04) but did not correlate with pressure measurements and clinical course. Recovery time correlated with clinical course (proximal calf, P=.03; distal calf, P=.005). CONCLUSION: MR angiography and MR spectroscopy of deoxymyoglobin can help document changes in visualization of collateral vessels and tissue perfusion in patients with CLI.

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.

High-resolution summer temperature reconstruction from Lake Silvaplana based on in-situ reflectance spectroscopy

Annually laminated (varved) sediments of proglacial Lake Silvaplana (46 ̊27’N, 9 ̊48’E, 1791 m a.s.l., Engadine, eastern Swiss Alps) provide an excellent archive for quantitative high-resolution (seasonal – annual) reconstruction of high- and lowfrequency climate signals back to AD 1580. The chronology of the core is based on varve counting, Cs-137, Pb-210 and event stratigraphy. In this study we present a reconstruction based on in-situ reflectance spectroscopy. In situ reflectance spectroscopy is known as a cost- and time-effective non destructtive method for semi-quantitative analysis of pigments (e.g., chlorines and carotenoids) and of lithoclastic sediment fractions. Reflectance-dependent absorption (RDA) was measured with a Gretac Macbeth spectrolino at 2 mm resolution. The spectral coverage ranges from 380 nm to 730 nm at 10 nm band resolution. In proglacial Lake Silvaplana, 99% of the sediment is lithoclastic prior to AD 1950. Therefore, we concentrate on absorption features that are characteristic for lithoclastic sediment fractions. In Lake Silvaplana, two significant correlations that are stable in time were found between RDA typical for lithoclastics and meteorological data: (1) the time series R 570 /R 630 (ratio between RDA at 570 nm and 630 nm) of varves in Lake Silvaplana and May to October temperatures at nearby station of Sils correlate highly significantly (calibration period AD 1864 – 1951, r = 0.74, p < 0.01 for 5ptsmoothed series; RMSE is 0.28 ̊C, RE = 0.41 and CE = 0.38), and (2) the minimum reflectance within the 690nm band (min690) data correlate with May to October (calibration period AD 1864 – 1951, r = 0.68, p < 0.01 for 5pt-smoothed series; RMSE = 0.22 ̊C, RE = 0.5, CE = 0.31). Both proxy series (min690nm and R 570 /R 630 values) are internally highly consistent (r = 0.8, p < 0.001). In proglacial Lake Silvaplana the largest amount of sediment is transported by glacial meltwater. The melting season spans approximately from May to October, which gives us a good understanding of the geophysical processes explaining the correlations between lithoclastic proxies and the meteorological data. The reconstructions were extended back to AD 1580 and show a broad corresponddence with fully independent reconstructions from tree rings and documentary data.

Electro-oxidation of Au(1 1 1) in contact with aqueous electrolytes: New insight from in situ vibration spectroscopy

We carried out a comprehensive study of Au(1 1 1) oxidation–reduction in the presence of (hydrogen-) sulfate ions on ideally smooth and stepped Au(S)[n(1 1 1)-(1 1 1)] single crystal electrodes using cyclic voltammetry, in situ scanning tunneling microscopy (STM) and vibration spectroscopy, such as surface-enhanced infrared absorption spectroscopy (SEIRAS) and shell-isolated nanoparticle-enhanced Raman spectroscopy (SHINERS). Surface structure changes and the role of surface defects in the potential regions of double layer charging and gold oxidation/reduction are discussed based on cyclic voltammetry and in situ STM data. SEIRAS and SHINERS provide complementary information on the chemical nature of adsorbates. In particular, the potential-dependent formation and stability ranges of adsorbed sulfate, hydroxide-species and of gold surface oxide could be resolved in detail. Based on our experimental observations, we proposed new and extended mechanisms of gold surface oxidation and reduction in 1.0 M H2SO4 and 1.0 M Na2SO4.

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.

Real-time observation of the charge-transfer-to-solvent dynamics

Intermolecular electron-transfer reactions have a crucial role in biology, solution chemistry and electrochemistry. The first step of such reactions is the expulsion of the electron to the solvent, whose mechanism is determined by the structure and dynamical response of the latter. Here we visualize the electron transfer to water using ultrafast fluorescence spectroscopy with polychromatic detection from the ultraviolet to the visible region, upon photo-excitation of the so-called charge transfer to solvent states of aqueous iodide. The initial emission is short lived (~60 fs) and it relaxes to a broad distribution of lower-energy charge transfer to solvent states upon rearrangement of the solvent cage. This distribution reflects the inhomogeneous character of the solvent cage around iodide. Electron ejection occurs from the relaxed charge transfer to solvent states with lifetimes of 100–400 fs that increase with decreasing emission energy.

Electronic band dispersion of graphene nanoribbons via Fourier-transformed scanning tunneling spectroscopy

The electronic structure of atomically precise armchair graphene nanoribbons of width N=7 (7-AGNRs) are investigated by scanning tunneling spectroscopy (STS) on Au(111). We record the standing waves in the local density of states of finite ribbons as a function of sample bias and extract the dispersion relation of frontier electronic states by Fourier transformation. The wave-vector-dependent contributions from these states agree with density functional theory calculations, thus enabling the unambiguous assignment of the states to the valence band, the conduction band, and the next empty band with effective masses of 0.41±0.08me,0.40±0.18me, and 0.20±0.03me, respectively. By comparing the extracted dispersion relation for the conduction band to corresponding height-dependent tunneling spectra, we find that the conduction band edge can be resolved only at small tip-sample separations and has not been observed before. As a result, we report a band gap of 2.37±0.06 eV for 7-AGNRs adsorbed on Au(111).

Excitonic Splitting, Delocalization, and Vibronic Quenching in the Benzonitrile Dimer

The excitonic S1/S2 state splitting and the localization/delocalization of the S1 and S2 electronic states are investigated in the benzonitrile dimer (BN)2 and its 13C and d5 isotopomers by mass-resolved two-color resonant two-photon ionization spectroscopy in a supersonic jet, complemented by calculations. The doubly hydrogen-bonded (BN-h5)2 and (BN-d5)2 dimers are C2h symmetric with equivalent BN moieties. Only the S0 → S2 electronic origin is observed, while the S0 → S1 excitonic component is electric-dipole forbidden. A single 12C/13C or 5-fold h5/d5 isotopic substitution reduce the dimer symmetry to Cs, so that the heteroisotopic dimers (BN)2-(h5 – h513C), (BN)2-(h5 – d5), and (BN)2-(h5 – h513C) exhibit both S0 → S1 and S0 → S2 origins. Isotope-dependent contributions Δiso to the excitonic splittings arise from the changes of the BN monomer zero-point vibrational energies; these range from Δiso(12C/13C) = 3.3 cm–1 to Δiso(h5/d5) = 155.6 cm–1. The analysis of the experimental S1/S2 splittings of six different isotopomeric dimers yields the S1/S2 exciton splitting Δexc = 2.1 ± 0.1 cm–1. Since Δiso(h5/d5) ≫ Δexc and Δiso(12C/13C) > Δexc, complete and near-complete exciton localization occurs upon 12C/13C and h5/d5 substitutions, respectively, as diagnosed by the relative S0 → S1 and S0 → S2 origin band intensities. The S1/S2 electronic energy gap of (BN)2 calculated by the spin-component scaled approximate second-order coupled-cluster (SCS-CC2) method is Δelcalc = 10 cm–1. This electronic splitting is reduced by the vibronic quenching factor Γ. The vibronically quenched exciton splitting Δelcalc·Γ = Δvibroncalc = 2.13 cm–1 is in excellent agreement with the observed splitting Δexc = 2.1 cm–1. The excitonic splittings can be converted to semiclassical exciton hopping times; the shortest hopping time is 8 ps for the homodimer (BN-h5)2, the longest is 600 ps for the (BN)2(h5 – d5) heterodimer.