Metabolite T1 relaxation times decrease across the adult lifespan

Relaxation correction is an integral step in quantifying brain metabolite concentrations measured by in vivo magnetic resonance spectroscopy (MRS). While most quantification routines assume constant T1 relaxation across age, it is possible that aging alters T1 relaxation rates, as is seen for T2 relaxation. Here, we investigate the age dependence of metabolite T1 relaxation times at 3 T in both gray- and white-matter-rich voxels using publicly available metabolite and metabolite-nulled (single inversion recovery TI = 600 ms) spectra acquired at 3 T using Point RESolved Spectroscopy (PRESS) localization. Data were acquired from voxels in the posterior cingulate cortex (PCC) and centrum semiovale (CSO) in 102 healthy volunteers across 5 decades of life (aged 20-69 years). All spectra were analyzed in Osprey v.2.4.0. To estimate T1 relaxation times for total N-acetyl aspartate at 2.0 ppm (tNAA2.0) and total creatine at 3.0 ppm (tCr3.0), the ratio of modeled metabolite residual amplitudes in the metabolite-nulled spectrum to the full metabolite signal was calculated using the single-inversion-recovery signal equation. Correlations between T1 and subject age were evaluated. Spearman correlations revealed that estimated T1 relaxation times of tNAA2.0 (rs = -0.27; p < 0.006) and tCr3.0 (rs = -0.40; p < 0.001) decreased significantly with age in white-matter-rich CSO, and less steeply for tNAA2.0 (rs = -0.228; p = 0.005) and (not significantly for) tCr3.0 (rs = -0.13; p = 0.196) in graymatter-rich PCC. The analysis harnessed a large publicly available cross-sectional dataset to test an important hypothesis, that metabolite T1 relaxation times change with age. This preliminary study stresses the importance of further work to measure age-normed metabolite T1 relaxation times for accurate quantification of metabolite levels in studies of aging.

Impact of acquisition and modeling parameters on the test-retest reproducibility of edited GABA

Literature values vary widely for within-subject test-retest reproducibility of gamma-aminobutyric acid (GABA) measured with edited magnetic resonance spectroscopy (MRS). Reasons for this variation remain unclear. Here, we tested whether three acquisition parameters-(1) sequence complexity (two-experiment MEscher-GArwood Point RESolved Spectroscopy [MEGA-PRESS] vs. four-experiment Hadamard Encoding and Reconstruction of MEGA-Edited Spectroscopy [HERMES]); (2) editing pulse duration (14 vs. 20 ms); and (3) scanner frequency drift (interleaved water referencing [IWR] turned ON vs. OFF)-and two linear combination modeling variations-(1) three different coedited macromolecule models (called "1to1GABA", "1to1GABAsoft", and "3to2MM" in the Osprey software package); and (2) 0.55- versus 0.4-ppm spline baseline knot spacing-affected the within-subject coefficient of variation of GABA + macromolecules (GABA+). We collected edited MRS data from the dorsal anterior cingulate cortex from 20 participants (mean age: 30.8 ± 9.5 years; 10 males). Test and retest scans were separated by removing the participant from the scanner for 5-10 min. Each acquisition consisted of two MEGA-PRESS and two HERMES sequences with editing pulse durations of 14 and 20 ms (referred to here as MEGA-14, MEGA-20, HERMES-14, and HERMES-20; all TE = 80 ms, 224 averages). We identified the best test-retest reproducibility following postprocessing with a composite model of the 0.9- and 3-ppm macromolecules ("3to2MM"); this model performed particularly well for the HERMES data. Furthermore, sparser (0.55- compared with 0.4-ppm) spline baseline knot spacing yielded generally better test-retest reproducibility for GABA+. Replicating our prior results, linear combination modeling in Osprey compared with simple peak fitting in Gannet resulted in substantially better test-retest reproducibility. However, reproducibility did not consistently differ for MEGA-PRESS compared with HERMES, for 14- compared with 20-ms editing pulses, or for IWR-ON versus IWR-OFF. These results highlight the importance of model selection for edited MRS studies of GABA+, particularly for clinical studies that focus on individual patient differences in GABA+ or changes following an intervention.

Neurochemical predictors of generalised learning induced by brain stimulation and training

Methods of cognitive enhancement for humans are most impactful when they generalise across tasks. However, the extent to which such "transfer" is possible via interventions is widely debated. In addition, the contribution of excitatory and inhibitory processes to such transfer is unknown. Here, in a large-scale neuroimaging individual differences study with humans (both sexes), we paired multitasking training and non-invasive brain stimulation (transcranial direct current stimulation; tDCS) over multiple days and assessed performance across a range of paradigms. In addition, we varied tDCS dosage (1.0 mA and 2.0 mA), electrode montage (left or right prefrontal regions), and training task (multitasking versus a control task) and assessed GABA and glutamate concentrations via ultra-high field 7T magnetic resonance spectroscopy. Generalised benefits were observed in spatial attention, indexed by visual search performance, when multitasking training was combined with 1.0 mA stimulation targeting either the left or right prefrontal cortex. This transfer effect persisted for ∼30 days post-intervention. Critically, the transferred benefits associated with right prefrontal tDCS were predicted by pre-training concentrations of glutamate in the prefrontal cortex. Thus, the effects of this combined stimulation and training protocol appears to be linked predominantly to excitatory brain processes.Significance statement Despite the general public's fascination with cognitive training, performance benefits rarely extend beyond the trained task, i.e., 'transfer'. Our study examines the impact of combining executive function training and transcranial direct current stimulation (tDCS) on human cognitive performance and identifies a functional neural metabolite marker (glutamate concentrations in prefrontal cortex assessed via 7 T MR Spectroscopy) that predicts outcomes. In the largest study of its kind to date (178 individuals), we find generalised performance benefits induced by frontal tDCS for an untrained spatial attention task. Further, the degree of transfer correlated with concentrations of glutamate in the frontal cortex. Thus, excitatory neural processes in this region are implicated in the transfer of paired stimulation and training benefits.

Neurotransmitter levels in the basal ganglia are associated with intracortical circuit activity of the primary motor cortex in healthy humans

BACKGROUND

The basal ganglia are strongly connected to the primary motor cortex (M1) and play a crucial role in movement control. Interestingly, several disorders showing abnormal neurotransmitter levels in basal ganglia also present concomitant anomalies in intracortical function within M1.

OBJECTIVE/HYPOTHESIS

The main aim of this study was to clarify the relationship between neurotransmitter content in the basal ganglia and intracortical function at M1 in healthy individuals. We hypothesized that neurotransmitter content of the basal ganglia would be significant predictors of M1 intracortical function.

METHODS

We combined magnetic resonance spectroscopy (MRS) and transcranial magnetic stimulation (TMS) to test this hypothesis in 20 healthy adults. An extensive TMS battery probing common measures of intracortical, and corticospinal excitability was administered, and GABA and glutamate-glutamine levels were assessed from voxels placed over the basal ganglia and the occipital cortex (control region).

RESULTS

Regression models using metabolite concentration as predictor and TMS metrics as outcome measures showed that glutamate level in the basal ganglia significantly predicted short interval intracortical inhibition (SICI) and intracortical facilitation (ICF), while GABA content did not. No model using metabolite measures from the occipital control voxel was significant.

CONCLUSIONS

Taken together, these results converge with those obtained in clinical populations and suggest that intracortical circuits in human M1 are associated with the neurotransmitter content of connected but distal subcortical structures crucial for motor function.

Feasibility and implications of using subject-specific macromolecular spectra to model short echo time magnetic resonance spectroscopy data

Expert consensus recommends linear-combination modeling (LCM) of 1 H MR spectra with sequence-specific simulated metabolite basis function and experimentally derived macromolecular (MM) basis functions. Measured MM basis functions are usually derived from metabolite-nulled spectra averaged across a small cohort. The use of subject-specific instead of cohort-averaged measured MM basis functions has not been studied widely. Furthermore, measured MM basis functions are not widely available to non-expert users, who commonly rely on parameterized MM signals internally simulated by LCM software. To investigate the impact of the choice of MM modeling, this study, therefore, compares metabolite level estimates between different MM modeling strategies (cohort-mean measured; subject-specific measured; parameterized) in a lifespan cohort and characterizes its impact on metabolite-age associations. 100 conventional (TE = 30 ms) and metabolite-nulled (TI = 650 ms) PRESS datasets, acquired from the medial parietal lobe in a lifespan cohort (20-70 years of age), were analyzed in Osprey. Short-TE spectra were modeled in Osprey using six different strategies to consider the MM baseline. Fully tissue- and relaxation-corrected metabolite levels were compared between MM strategies. Model performance was evaluated by model residuals, the Akaike information criterion (AIC), and the impact on metabolite-age associations. The choice of MM strategy had a significant impact on the mean metabolite level estimates and no major impact on variance. Correlation analysis revealed moderate-to-strong agreement between different MM strategies (r > 0.6). The lowest relative model residuals and AIC values were found for the cohort-mean measured MM. Metabolite-age associations were consistently found for two major singlet signals (total creatine (tCr])and total choline (tCho)) for all MM strategies; however, findings for metabolites that are less distinguishable from the background signals associations depended on the MM strategy. A variance partition analysis indicated that up to 44% of the total variance was related to the choice of MM strategy. Additionally, the variance partition analysis reproduced the metabolite-age association for tCr and tCho found in the simpler correlation analysis. In summary, the inclusion of a single high signal-to-noise ratio MM basis function (cohort-mean) in the short-TE LCM leads to more lower model residuals and AIC values compared with MM strategies with more degrees of freedom (Gaussian parametrization) or subject-specific MM information. Integration of multiple LCM analyses into a single statistical model potentially allows to identify the robustness in the detection of underlying effects (e.g., metabolite vs. age), reduces algorithm-based bias, and estimates algorithm-related variance.

Neurochemical and functional reorganization of the cognitive-ear link underlies cognitive impairment in presbycusis

Recent studies suggest that the interaction between presbycusis and cognitive impairment may be partially explained by the cognitive-ear link. However, the underlying neurophysiological mechanisms remain largely unknown. In this study, we combined magnetic resonance spectroscopy (MRS) and resting-state functional magnetic resonance imaging (fMRI) to investigate auditory gamma-aminobutyric acid (GABA) and glutamate (Glu) levels, intra- and inter-network functional connectivity, and their relationships with auditory and cognitive function in 51 presbycusis patients and 51 well-matched healthy controls. Our results confirmed reorganization of the cognitive-ear link in presbycusis, including decreased auditory GABA and Glu levels and aberrant functional connectivity involving auditory networks (AN) and cognitive-related networks, which were associated with reduced speech perception or cognitive impairment. Moreover, mediation analyses revealed that decreased auditory GABA levels and dysconnectivity between the AN and default mode network (DMN) mediated the association between hearing loss and impaired information processing speed in presbycusis. These findings highlight the importance of AN-DMN dysconnectivity in cognitive-ear link reorganization leading to cognitive impairment, and hearing loss may drive reorganization via decreased auditory GABA levels. Modulation of GABA neurotransmission may lead to new treatment strategies for cognitive impairment in presbycusis patients.

Impact of gradient scheme and non-linear shimming on out-of-voxel echo artifacts in edited MRS

Out-of-voxel (OOV) signals are common spurious echo artifacts in MRS. These signals often manifest in the spectrum as very strong "ripples," which interfere with spectral quantification by overlapping with targeted metabolite resonances. Dephasing optimization through coherence order pathway selection (DOTCOPS) gradient schemes are algorithmically optimized to suppress all potential alternative coherence transfer pathways (CTPs), and should suppress unwanted OOV echoes. In addition, second-order shimming uses non-linear gradient fields to maximize field homogeneity inside the voxel, which unfortunately increases the diversity of local gradient fields outside of the voxel. Given that strong local spatial B0 gradients can refocus unintended CTPs, it is possible that OOVs are less prevalent when only linear first-order shimming is applied. Here we compare the size of unwanted OOV signals in Hadamard-edited (HERMES) data acquired with either a local gradient scheme (which we refer to here as "Shared") or DOTCOPS, and with first- or second-order shimming. We collected data from 15 healthy volunteers in two brain regions (voxel size 30 × 26 × 26 mm3 ) from which it is challenging to acquire MRS data: medial prefrontal cortex and left temporal cortex. Characteristic OOV echoes were seen in both GABA- and GSH-edited spectra for both brain regions, gradient schemes, and shimming approaches. A linear mixed-effect model revealed a statistically significant difference in the average residual based on the gradient scheme in both GABA- (p < 0.001) and GSH-edited (p < 0.001) spectra: that is, the DOTCOPS gradient scheme resulted in smaller OOV artifacts compared with the Shared scheme. There were no significant differences in OOV artifacts associated with shimming method. Thus, these results suggest that the DOTCOPS gradient scheme for J-difference-edited PRESS acquisitions yields spectra with smaller OOV echo artifacts than the Shared gradient scheme implemented in a widely disseminated editing sequence.

Brain total creatine differs between primary progressive aphasia (PPA) subtypes and correlates with disease severity

Primary progressive aphasia (PPA) is comprised of three subtypes: logopenic (lvPPA), non-fluent (nfvPPA), and semantic (svPPA). We used magnetic resonance spectroscopy (MRS) to measure tissue-corrected metabolite levels in the left inferior frontal gyrus (IFG) and right sensorimotor cortex (SMC) from 61 PPA patients. We aimed to: (1) characterize subtype differences in metabolites; and (2) test for metabolite associations with symptom severity. tCr differed by subtype across the left IFG and right SMC. tCr levels were lowest in lvPPA and highest in svPPA. tCr levels predicted lvPPA versus svPPA diagnosis. Higher IFG tCr and lower Glx correlated with greater disease severity. As tCr is involved in brain energy metabolism, svPPA pathology might involve changes in specific cellular energy processes. Perturbations to cellular energy homeostasis in language areas may contribute to symptoms. Reduced cortical excitatory capacity (i.e. lower Glx) in language regions may also contribute to symptoms. Thus, tCr may be useful for differentiating between PPA subtypes, and both tCr and Glx might have utility in understanding PPA mechanisms and tracking progression.

Impact of acquisition and modeling parameters on test-retest reproducibility of edited GABA

Literature values for within-subject test-retest reproducibility of gamma-aminobutyric acid (GABA), measured with edited magnetic resonance spectroscopy (MRS), vary widely. Reasons for this variation remain unclear. Here we tested whether sequence complexity (two-experiment MEGA-PRESS versus four-experiment HERMES), editing pulse duration (14 versus 20 ms), scanner frequency drift (interleaved water referencing (IWR) turned ON versus OFF), and linear combination modeling variations (three different co-edited macromolecule models and 0.55 versus 0.4 ppm spline baseline knot spacing) affected the within-subject coefficient of variation of GABA + macromolecules (GABA+). We collected edited MRS data from the dorsal anterior cingulate cortex from 20 participants (30.8 ± 9.5 years; 10 males). Test and retest scans were separated by removing the participant from the scanner for 5-10 minutes. Each acquisition consisted of two MEGA-PRESS and two HERMES sequences with editing pulse durations of 14 and 20 ms (referred to here as: MEGA-14, MEGA-20, HERMES-14, and HERMES-20; all TE = 80 ms, 224 averages). Reproducibility did not consistently differ for MEGA-PRESS compared with HERMES or for 14 compared with 20 ms editing pulses. A composite model of the 0.9 and 3 ppm macromolecules (particularly for HERMES) and sparser (0.55 compared with 0.4 ppm) spline baseline knot spacing yielded generally better test-retest reproducibility for GABA+. Replicating our prior results, linear combination modeling in Osprey compared with simple peak fitting in Gannet resulted in substantially better test-retest reproducibility. These results highlight the importance of model selection for edited MRS studies of GABA+, particularly for clinical studies which focus on individual patient differences in GABA+ or changes following an intervention.

Neurometabolic timecourse of healthy aging

PURPOSE

The neurometabolic timecourse of healthy aging is not well-established, in part due to diversity of quantification methodology. In this study, a large structured cross-sectional cohort of male and female subjects throughout adulthood was recruited to investigate neurometabolic changes as a function of age, using consensus-recommended magnetic resonance spectroscopy quantification methods.

METHODS

102 healthy volunteers, with approximately equal numbers of male and female participants in each decade of age from the 20s, 30s, 40s, 50s, and 60s, were recruited with IRB approval. MR spectroscopic data were acquired on a 3T MRI scanner. Metabolite spectra were acquired using PRESS localization (TE=30 ms; 96 transients) in the centrum semiovale (CSO) and posterior cingulate cortex (PCC). Water-suppressed spectra were modeled using the Osprey algorithm, employing a basis set of 18 simulated metabolite basis functions and a cohort-mean measured macromolecular spectrum. Pearson correlations were conducted to assess relationships between metabolite concentrations and age for each voxel; Spearman correlations were conducted where metabolite distributions were non-normal. Paired t-tests were run to determine whether metabolite concentrations differed between the PCC and CSO. Finally, robust linear regressions were conducted to assess both age and sex as predictors of metabolite concentrations in the PCC and CSO and separately, to assess age, signal-noise ratio, and full width half maximum (FWHM) linewidth as predictors of metabolite concentrations.

RESULTS

Data from four voxels were excluded (2 ethanol; 2 unacceptably large lipid signal). Statistically-significant age*metabolite Pearson correlations were observed for tCho (r(98)=0.33, p<0.001), tCr (r(98)=0.60, p<0.001), and mI (r(98)=0.32, p=0.001) in the CSO and for NAAG (r(98)=0.26, p=0.008), tCho(r(98)=0.33, p<0.001), tCr (r(98)=0.39, p<0.001), and Gln (r(98)=0.21, p=0.034) in the PCC. Spearman correlations for non-normal variables revealed a statistically significant correlation between sI and age in the CSO (r(86)=0.26, p=0.013). No significant correlations were seen between age and tNAA, NAA, Glx, Glu, GSH, PE, Lac, or Asp in either region (all p>0.20). Age associations for tCho, tCr, mI and sI in the CSO and for NAAG, tCho, and tCr in the PCC remained when controlling for sex in robust regressions. CSO NAAG and Asp, as well as PCC tNAA, sI, and Lac were higher in women; PCC Gln was higher in men. When including an age*sex interaction term in robust regression models, a significant age*sex interaction was seen for tCho (F(1,96)=11.53, p=0.001) and GSH (F(1,96)=7.15, p=0.009) in the CSO and tCho (F(1,96)=9.17, p=0.003), tCr (F(1,96)=9.59, p=0.003), mI (F(1,96)=6.48, p=0.012), and Lac (F(1,78)=6.50, p=0.016) in the PCC. In all significant interactions, metabolite levels increased with age in females, but not males. There was a significant positive correlation between linewidth and age. Age relationships with tCho, tCr, and mI in the CSO and tCho, tCr, mI, and sI in the PCC were significant after controlling for linewidth and FWHM in robust regressions.

CONCLUSION

The primary (correlation) results indicated age relationships for tCho, tCr, mI, and sI in the CSO and for NAAG, tCho, tCr, and Gln in the PCC, while no age correlations were found for tNAA, NAA, Glx, Glu, GSH, PE, Lac, or Asp in either region. Our results provide a normative foundation for future work investigating the neurometabolic time course of healthy aging using MRS.

Edited magnetic resonance spectroscopy in the neonatal brain

J-difference-edited spectroscopy is a valuable approach for the detection of low-concentration metabolites with magnetic resonance spectroscopy (MRS). Currently, few edited MRS studies are performed in neonates due to suboptimal signal-to-noise ratio, relatively long acquisition times, and vulnerability to motion artifacts. Nonetheless, the technique presents an exciting opportunity in pediatric imaging research to study rapid maturational changes of neurotransmitter systems and other metabolic systems in early postnatal life. Studying these metabolic processes is vital to understanding the widespread and rapid structural and functional changes that occur in the first years of life. The overarching goal of this review is to provide an introduction to edited MRS for neonates, including the current state-of-the-art in editing methods and editable metabolites, as well as to review the current literature applying edited MRS to the neonatal brain. Existing challenges and future opportunities, including the lack of age-specific reference data, are also discussed.

Importance of Linear Combination Modeling for Quantification of Glutathione and γ-Aminobutyric Acid Levels Using Hadamard-Edited Magnetic Resonance Spectroscopy

BACKGROUND

J-difference-edited ¹H-MR spectra require modeling to quantify signals of low-concentration metabolites. Two main approaches are used for this spectral modeling: simple peak fitting and linear combination modeling (LCM) with a simulated basis set. Recent consensus recommended LCM as the method of choice for the spectral analysis of edited data.

PURPOSE

The aim of this study is to compare the performance of simple peak fitting and LCM in a test-retest dataset, hypothesizing that the more sophisticated LCM approach would improve quantification of Hadamard-edited data compared with simple peak fitting.

METHODS

A test-retest dataset was re-analyzed using Gannet (simple peak fitting) and Osprey (LCM). These data were obtained from the dorsal anterior cingulate cortex of twelve healthy volunteers, with TE = 80 ms for HERMES and TE = 120 ms for MEGA-PRESS of glutathione (GSH). Within-subject coefficients of variation (CVs) were calculated to quantify between-scan reproducibility of each metabolite estimate.

RESULTS

The reproducibility of HERMES GSH estimates was substantially improved using LCM compared to simple peak fitting, from a CV of 19.0-9.9%. For MEGA-PRESS GSH data, reproducibility was similar using LCM and simple peak fitting, with CVs of 7.3 and 8.8%. GABA + CVs from HERMES were 16.7 and 15.2%, respectively for the two models.

CONCLUSION

LCM with simulated basis functions substantially improved the reproducibility of GSH quantification for HERMES data.

Comparison of linear combination modeling strategies for edited magnetic resonance spectroscopy at 3 T

J-difference-edited spectroscopy is a valuable approach for the in vivo detection of γ-aminobutyric-acid (GABA) with magnetic resonance spectroscopy (MRS). A recent expert consensus article recommends linear combination modeling (LCM) of edited MRS but does not give specific details regarding implementation. This study explores different modeling strategies to adapt LCM for GABA-edited MRS. Sixty-one medial parietal lobe GABA-edited MEGA-PRESS spectra from a recent 3-T multisite study were modeled using 102 different strategies combining six different approaches to account for co-edited macromolecules (MMs), three modeling ranges, three baseline knot spacings, and the use of basis sets with or without homocarnosine. The resulting GABA and GABA+ estimates (quantified relative to total creatine), the residuals at different ranges, standard deviations and coefficients of variation (CVs), and Akaike information criteria, were used to evaluate the models' performance. Significantly different GABA+ and GABA estimates were found when a well-parameterized MM_3co basis function was included in the model. The mean GABA estimates were significantly lower when modeling MM_3co , while the CVs were similar. A sparser spline knot spacing led to lower variation in the GABA and GABA+ estimates, and a narrower modeling range-only including the signals of interest-did not substantially improve or degrade modeling performance. Additionally, the results suggest that LCM can separate GABA and the underlying co-edited MM_3co . Incorporating homocarnosine into the modeling did not significantly improve variance in GABA+ estimates. In conclusion, GABA-edited MRS is most appropriately quantified by LCM with a well-parameterized co-edited MM_3co basis function with a constraint to the nonoverlapped MM_0.93 , in combination with a sparse spline knot spacing (0.55 ppm) and a modeling range of 0.5-4 ppm.

The macromolecular MR spectrum does not change with healthy aging

PURPOSE

To acquire the mobile macromolecule (MM) spectrum from healthy participants, and to investigate changes in the signals with age and sex.

METHODS

102 volunteers (49 M/53 F) between 20 and 69 years were recruited for in vivo data acquisition in the centrum semiovale (CSO) and posterior cingulate cortex (PCC). Spectral data were acquired at 3T using PRESS localization with a voxel size of 30 × 26 × 26 mm³ , pre-inversion (TR/TI 2000/600 ms) and CHESS water suppression. Metabolite-nulled spectra were modeled to eliminate residual metabolite signals, which were then subtracted out to yield a "clean" MM spectrum using the Osprey software. Pearson's correlation coefficient was calculated between integrals and age for the 14 MM signals. One-way ANOVA was performed to determine differences between age groups. An independent t-test was carried out to determine differences between sexes.

RESULTS

MM spectra were successfully acquired in 99 (CSO) and 96 (PCC) of 102 subjects. No significant correlations were seen between age and MM signals. One-way ANOVA also suggested no age-group differences for any MM peak (all p > .004). No differences were observed between sex groups. WM and GM voxel fractions showed a significant (p < .05) negative linear association with age in the WM-predominant CSO (R = -0.29) and GM-predominant PCC regions (R = -0.57) respectively while CSF increased significantly with age in both regions.

CONCLUSION

Our findings suggest that a pre-defined MM basis function can be used for linear combination modeling of metabolite data from different age and sex groups.

Frequency drift in MR spectroscopy at 3T

PURPOSE

Heating of gradient coils and passive shim components is a common cause of instability in the B0 field, especially when gradient intensive sequences are used. The aim of the study was to set a benchmark for typical drift encountered during MR spectroscopy (MRS) to assess the need for real-time field-frequency locking on MRI scanners by comparing field drift data from a large number of sites.

METHOD

A standardized protocol was developed for 80 participating sites using 99 3T MR scanners from 3 major vendors. Phantom water signals were acquired before and after an EPI sequence. The protocol consisted of: minimal preparatory imaging; a short pre-fMRI PRESS; a ten-minute fMRI acquisition; and a long post-fMRI PRESS acquisition. Both pre- and post-fMRI PRESS were non-water suppressed. Real-time frequency stabilization/adjustment was switched off when appropriate. Sixty scanners repeated the protocol for a second dataset. In addition, a three-hour post-fMRI MRS acquisition was performed at one site to observe change of gradient temperature and drift rate. Spectral analysis was performed using MATLAB. Frequency drift in pre-fMRI PRESS data were compared with the first 5:20 minutes and the full 30:00 minutes of data after fMRI. Median (interquartile range) drifts were measured and showed in violin plot. Paired t-tests were performed to compare frequency drift pre- and post-fMRI. A simulated in vivo spectrum was generated using FID-A to visualize the effect of the observed frequency drifts. The simulated spectrum was convolved with the frequency trace for the most extreme cases. Impacts of frequency drifts on NAA and GABA were also simulated as a function of linear drift. Data from the repeated protocol were compared with the corresponding first dataset using Pearson's and intraclass correlation coefficients (ICC).

RESULTS

Of the data collected from 99 scanners, 4 were excluded due to various reasons. Thus, data from 95 scanners were ultimately analyzed. For the first 5:20 min (64 transients), median (interquartile range) drift was 0.44 (1.29) Hz before fMRI and 0.83 (1.29) Hz after. This increased to 3.15 (4.02) Hz for the full 30 min (360 transients) run. Average drift rates were 0.29 Hz/min before fMRI and 0.43 Hz/min after. Paired t-tests indicated that drift increased after fMRI, as expected (p < 0.05). Simulated spectra convolved with the frequency drift showed that the intensity of the NAA singlet was reduced by up to 26%, 44 % and 18% for GE, Philips and Siemens scanners after fMRI, respectively. ICCs indicated good agreement between datasets acquired on separate days. The single site long acquisition showed drift rate was reduced to 0.03 Hz/min approximately three hours after fMRI.

DISCUSSION

This study analyzed frequency drift data from 95 3T MRI scanners. Median levels of drift were relatively low (5-min average under 1 Hz), but the most extreme cases suffered from higher levels of drift. The extent of drift varied across scanners which both linear and nonlinear drifts were observed.

Implications of Abdominal Adipose Tissue Distribution on Nonalcoholic Fatty Liver Disease and Metabolic Syndrome: A Chinese General Population Study

INTRODUCTION

Visceral adipose tissue (VAT) has been found to play a critical role in the development of metabolic syndrome and nonalcoholic fatty liver disease (NAFLD) independent of generalized obesity.

METHODS

In this secondary study of prospectively acquired data, 625 participants underwent magnetic resonance spectroscopy and chemical shift fat-water separation MRI (2-point Dixon) of the liver and whole abdomen, respectively, in a 3 Tesla magnet. Whole abdominal VAT and subcutaneous adipose tissue (SAT) were extracted from the 2-point Dixon image series using an automated method. Clinical/anthropometric/blood biochemistry parameters were measured. Using region-specific body mass index, participants were classified into 3 paired subgroups (lean, overweight, and obese) and presence of NAFLD (liver fat content ≥ 5.5%).

RESULTS

All relevant clinical/anthropometric/blood biochemistry characteristics and liver enzymes were statistically significant between groups (P < 0.001). NAFLD was found in 12.1%, 43.8%, and 68.3% and metabolic syndrome in 51.1%, 61.9%, and 65% of the lean, overweight, and obese, respectively. Odds ratio for metabolic syndrome and NAFLD was increased by 2.73 (95% confidence interval [CI] 2.18-3.40) and 2.53 (95% CI 2.04-3.12), respectively, for 1SD increase in VAT volume while prevalence of metabolic syndrome was increased by 2.26 (95% CI 1.83-2.79) for 1SD increase in liver fat content (%). VAT/SAT ratio in the lean with fatty liver showed the highest ratio (0.54) among all the subgroups, without a significant difference between the lean and obese with NAFLD (P = 0.127).

DISCUSSION

Increased VAT volume/disproportional distribution of VAT/SAT may be vital drivers to the development of metabolic syndrome and NAFLD irrespective of body mass index category.

Frequency and phase correction of J-difference edited MR spectra using deep learning

PURPOSE

To investigate whether a deep learning-based (DL) approach can be used for frequency-and-phase correction (FPC) of MEGA-edited MRS data.

METHODS

Two neural networks (1 for frequency, 1 for phase) consisting of fully connected layers were trained and validated using simulated MEGA-edited MRS data. This DL-FPC was subsequently tested and compared to a conventional approach (spectral registration [SR]) and to a model-based SR implementation (mSR) using in vivo MEGA-edited MRS datasets. Additional artificial offsets were added to these datasets to further investigate performance.

RESULTS

The validation showed that DL-based FPC was capable of correcting within 0.03 Hz of frequency and 0.4°of phase offset for unseen simulated data. DL-based FPC performed similarly to SR for the unmanipulated in vivo test datasets. When additional offsets were added to these datasets, the networks still performed well. However, although SR accurately corrected for smaller offsets, it often failed for larger offsets. The mSR algorithm performed well for larger offsets, which was because the model was generated from the in vivo datasets. In addition, the computation times were much shorter using DL-based FPC or mSR compared to SR for heavily distorted spectra.

CONCLUSION

These results represent a proof of principle for the use of DL for preprocessing MRS data.

Osprey: Open-source processing, reconstruction & estimation of magnetic resonance spectroscopy data

BACKGROUND

Processing and quantitative analysis of magnetic resonance spectroscopy (MRS) data are far from standardized and require interfacing with third-party software. Here, we present Osprey, a fully integrated open-source data analysis pipeline for MRS data, with seamless integration of pre-processing, linear-combination modelling, quantification, and data visualization.

NEW METHOD

Osprey loads multiple common MRS data formats, performs phased-array coil combination, frequency-and phase-correction of individual transients, signal averaging and Fourier transformation. Linear combination modelling of the processed spectrum is carried out using simulated basis sets and a spline baseline. The MRS voxel is coregistered to an anatomical image, which is segmented for tissue correction and quantification is performed based upon modelling parameters and tissue segmentation. The results of each analysis step are visualized in the Osprey GUI. The analysis pipeline is demonstrated in 12 PRESS, 11 MEGA-PRESS, and 8 HERMES datasets acquired in healthy subjects.

RESULTS

Osprey successfully loads, processes, models, and quantifies MRS data acquired with a variety of conventional and spectral editing techniques.

COMPARISON WITH EXISTING METHOD(S)

Osprey is the first MRS software to combine uniform pre-processing, linear-combination modelling, tissue correction and quantification into a coherent ecosystem. Compared to existing compiled, often closed-source modelling software, Osprey's open-source code philosophy allows researchers to integrate state-of-the-art data processing and modelling routines, and potentially converge towards standardization of analysis.

CONCLUSIONS

Osprey combines robust, peer-reviewed data processing methods into a modular workflow that is easily augmented by community developers, allowing the rapid implementation of new methods.

The relationship between pancreas steatosis and the risk of metabolic syndrome and insulin resistance in Chinese adolescents with concurrent obesity and non-alcoholic fatty liver disease

BACKGROUND

The incidence of childhood obesity and associated comorbidities are on an increasing trend worldwide. More than 340 million children and adolescents aged between 5 and 19 years old were overweight or had obesity in 2016, from which over 124 million children and adolescents (6% of girls and 8% of boys) had obesity.

OBJECTIVE

To describe the relationship between pancreas steatosis, body fat and the risk of metabolic syndrome, insulin resistance in Hong Kong Chinese adolescents with both obesity and non-alcoholic fatty liver disease (NAFLD).

METHODS

Fifty two adolescents with obesity and NAFLD were analysed (14-18 years), stratified into fatty and non-fatty pancreas groups using chemical shift encoded MRI-pancreas proton density fat fraction ≥5%. Pancreatic, abdominal subcutaneous adipose tissue (SAT)/visceral adipose tissue (VAT) volumes, biochemical and anthropometric parameters were measured. Mann-Whitney U test, multiple linear/binary logistic regression analyses and odds ratios were used.

RESULTS

Fifty percent had fatty pancreas, 38% had metabolic syndrome and 81% had insulin resistance. Liver proton density fat fraction (PDFF) and VAT were independent predictors of insulin resistance (P = .006, .016). Pancreas and liver PDFF were both independent predictors of beta cells dysfunction (P = .015, .050) and metabolic syndrome (P = .021, .041). Presence of fatty pancreas in obesity was associated with insulin resistance (OR = 1.58, 95% CI = 0.39-6.4) and metabolic syndrome (OR = 1.70, 95% CI = 0.53-5.5).

CONCLUSION

A significant causal relationship exists between fatty pancreas, fatty liver, body fat and the risk of developing metabolic syndrome and insulin resistance.

KEY POINTS

Fatty pancreas is a common finding in adolescents with obesity, with a prevalence rate of 50% in this study cohort. Liver PDFF and VAT are independent predictors of insulin resistance while pancreas PDFF and liver PDFF are independent predictors of both beta cells dysfunction and metabolic syndrome. Presence of fatty pancreas at imaging should not be considered as a benign finding but rather as an imaging biomarker of emerging pancreatic metabolic and endocrine dysfunction.

Effect of Age on GABA+ and Glutathione in a Pediatric Sample

BACKGROUND AND PURPOSE

Gamma-aminobutyric acid (GABA) is the primary inhibitory neurotransmitter in the human brain and is implicated in several neuropathologies. Glutathione is a major antioxidant in the brain and is considered a marker of oxidative stress. Several studies have reported age-related declines in GABA levels in adulthood, but the trajectory of both GABA and glutathione during childhood has not been well explored. The aim of this study is to establish how GABA and glutathione vary with age during early development.

MATERIALS AND METHODS

Twenty-three healthy children (5.6-13.9 years of age) were recruited for this study. MR imaging/MR spectroscopy experiments were conducted on a 3T MR scanner. A 27-mL MR spectroscopy voxel was positioned in the frontal lobe. J-difference edited MR spectroscopy was used to spectrally edit GABA and glutathione. Data were analyzed using the Gannet software, and GABA+ (GABA + macromolecules/homocarnosine) and glutathione were quantified using water (GABA+_H2O and Glutathione_H2O) and Cr (GABA+/Cr and glutathione/Cr) as concentration references. Also, the relative gray matter contribution to the voxel volume (GM_ratio) was estimated from structural images. Pearson correlation coefficients were used to examine the association between age and GABA+_H2O (and glutathione_H2O), between age and GABA+/Cr (and glutathione/Cr), and between age and GM_ratio.

RESULTS

Both GABA+_H2O (r = 0.63, P = .002) and GABA+/Cr (r = 0.48, P = .026) significantly correlated with age, whereas glutathione measurements and GM_ratio did not.

CONCLUSIONS

We demonstrate increases in GABA and no differences in glutathione with age in a healthy pediatric sample. This study provides insight into neuronal maturation in children and may facilitate better understanding of normative behavioral development and the pathophysiology of developmental disorders.

Anterior-posterior length discrepancy of the spinal column in adolescent idiopathic scoliosis-a 3D CT study

BACKGROUND CONTEXT

One of the characteristics of reported observations in adolescent idiopathic scoliosis (AIS) is that the thoracic spine is longer anteriorly than posteriorly, more pronounced around the apex than the transitional zones. This reversal of the normal kyphotic anatomy of the thoracic spine is related to questions of etiopathogenesis of AIS. The changes in the anatomy of the anterior column have been described rather in detail; however, the role of the posterior spinal column and the laminae has so far not been elucidated. If the posterior column exhibits a longitudinal growth disturbance, it could act as a tether, leading to a more or less normal anterior column with a deformed and shorter posterior aspect of the spine. So far, it has remained unclear whether this anterior-posterior length discrepancy is the result of relative anterior lengthening or relative posterior shortening, and which tissues (bone, disc, intervertebral soft tissue) are involved.

PURPOSE

The present study aimed to compare the discrepancy of the anterior-posterior length of the spinal column in the "true" midsagittal plane of each vertebra in patients with idiopathic scoliosis versus controls, using three-dimensional computed tomography (CT) scans.

STUDY DESIGN/SETTING

This is a cross-sectional study.

PATIENT SAMPLE

The sample consisted of computed tomography scans of 80 patients with moderate to severe AIS (Cobb angle: 46°-109°) before scoliosis navigation surgery and 30 non-scoliotic age-matched controls.

OUTCOME MEASURES

The height of the osseous and non-osseous structures from anterior to posterior in the "true" midsagittal plane has been determined: the anterior side of the vertebral body and disc, the posterior side of the vertebral body and disc, the lamina and interlaminar space and the spinous process and interspinous space, as well as the height ratios between the anterior column and posterior structures of the primary thoracic and lumbar AIS curves and corresponding levels in non-scoliotic controls.

METHODS

Semiautomatic software was used to reconstruct and measure the parameters in the true midsagittal plane of each vertebra and intervertebral structure that are rotated and tilted in a different way.

RESULTS

In AIS, the anterior height of the thoracic curve was 3.6±2.8% longer than the posterior height, 2.0±6.1% longer than the length along the laminae, and 8.7±7.1% longer than the length along the spinous processes, and this differed significantly from controls (-2.7±2.4%, -7.4±5.2%, and +0.7±7.8%; p<.001). The absolute height of the osseous parts did not differ significantly between AIS and controls in the midsagittal plane. In contrast, the intervertebral structures contributed significantly to the observed length discrepancies. In absolute lengths, the anterior side of the disc of the thoracic curve was higher in AIS (5.4±0.8 mm) than controls (4.8±1.0 mm; p<.001), whereas the interspinous space was smaller in AIS (12.3±1.4 mm vs. 14.0±1.6 mm; p<.001).

CONCLUSIONS

Based on this in vivo analysis, the true three-dimensional anterior-posterior length discrepancy of AIS curves was found to occur through both anterior column lengthening and posterior column shortening, with the facet joints functioning as the fulcrum. The vertebrae contribute partly to the anterior-posterior length discrepancy accompanied by more significant and possibly secondary increased anterior intervertebral discs height.

Validation of water-fat MRI and proton MRS in assessment of hepatic fat and the heterogeneous distribution of hepatic fat and iron in subjects with non-alcoholic fatty liver disease

BACKGROUND

Research studies demonstrated pathologic lesions were unevenly distributed in patients with non-alcoholic fatty liver disease (NAFLD) and non-alcoholic steatohepatitis. As hepatic steatosis occurs prior to steatohepatitis and other late stage liver conditions, the distribution pattern of hepatic fat and iron concentration should be investigated to prevent sampling variability. The first purpose of this study was to perform comparison and validation of in-house hepatic fat measurements using water-fat MRI and MRS. The second objective was to quantify hepatic fat-fraction and T2* values in left and right liver lobes using water-fat MRI.

METHOD

Fifty-four non-alcoholic adults (27 NAFLD, age: 42.8 ± 11.8), 27 non-NAFLD, age: 45.5 ± 11.2) and 46 non-alcoholic teenagers (23 NAFLD (age: 15.4 ± 2.6), 23 non-NAFLD (age: 13.9 ± 2.3) were recruited. All participants underwent chemical shift water-fat MRI and 1H MRS at 3 T. Hepatic steatosis was defined by intrahepatic triglyceride more than the threshold of 5.56% using MRS (clinical reference) and non-alcoholic was defined by alcohol ingestion of no more than 30 g and 20 g per day for male and female respectively. Hepatic fat-fractions in left and right liver lobes were measured using regions-of-interest (ROIs) approach. Three ROIs were drawn on the fat-fraction images and duplicated on to the co-registered T2* images at the inferior right, superior right and superior left liver lobes. Comparison and validation of water-fat MRI and MRS were performed using intraclass correlation coefficient (ICC) and Bland-Altman plot. Hepatic fat-fraction and T2* measured from the ROIs were compared using repeated measures ANOVA. Independent t-test was used for between groups analysis.

RESULTS

Statistical analysis indicated good correlation (R = 0.987) and agreement (ICC = 0.982) between MRS and water-fat MRI in hepatic fat measurements. Results indicated that hepatic fat was significantly higher in the right lobe compared to the left in NAFLD adults (p < 0.001) and NAFLD teenagers (p < 0.001). For T2*, significant difference between left and right lobes was observed in NAFLD adults (p < 0.001) and non-NAFLD adults (p < 0.001) but not in teenagers.

CONCLUSION

Hepatic fat measurements using MRS and water-fat MRI are statistically equivalent. In subjects with NAFLD regardless of their age, hepatic fat is stored preferentially in the right live lobe probably due to the streamline of blood flow to the right liver. T2* value is significantly higher in the right liver lobe in adults but not in the teenagers regardless of their hepatic fat contents probably due to the longer time span of hepatic iron accumulation.

Altered White Matter Microstructure in the Corpus Callosum and Its Cerebral Interhemispheric Tracts in Adolescent Idiopathic Scoliosis: Diffusion Tensor Imaging Analysis

BACKGROUND AND PURPOSE

Neural system was one of the important contributors to the etiopathogenesis of adolescent idiopathic scoliosis; additionally, the morphology of corpus callosum interconnecting both hemispheres of the brain was found to be altered morphologically. Our aim was to evaluate and compare the microstructural changes of the corpus callosum and its interhemispheric white matter fiber tracts interconnecting both cerebral hemispheres in patients with adolescent idiopathic scoliosis and matched controls using diffusion tensor imaging.

MATERIALS AND METHODS

Brain DTI was performed in 69 patients with adolescent idiopathic scoliosis (female, right thoracic/thoracolumbar curve) and 40 age-matched controls without adolescent idiopathic scoliosis (female). 2D and 3D segmentation of the corpus callosum were performed using a region-growing method, and the corpus callosum was further divided into 6 regions, including the rostrum, genu, anterior and posterior midbodies, isthmus, and splenium. The laterality index was calculated to quantify the asymmetry of the corpus callosum. Interhemispheric fiber tractography were performed using the Brodmann atlas.

RESULTS

2D ROI analysis revealed reduced fractional anisotropy in the genu and splenium (P = .075 and P = .024, respectively). Consistently reduced fractional anisotropy on the left sides of the genu and splenium was also found in 3D ROI analysis (P = .03 and P = .012, respectively). The laterality index analysis revealed a pseudo-right lateralization of the corpus callosum in adolescent idiopathic scoliosis. Interhemispheric fibers via the splenium interconnecting Brodmann 3, 1, and 2; Brodmann 17; and Brodmann 18 (corresponding to the primary somatosensory cortex and primary and secondary visual cortices) were also found to have reduced fractional anisotropy (P ≤ .05).

CONCLUSIONS

Reduced fractional anisotropy was found in the genu and splenium of the corpus callosum and corresponding interhemispheric fiber tracts interconnecting the somatosensory and visual cortices via the splenium. Our results are suggestive of altered white matter microstructure within the brain of those with adolescent idiopathic scoliosis, which could be related to abnormal brain maturation during adolescence in adolescent idiopathic scoliosis and could possibly explain the previously documented somatosensory function impairment and visuo-oculomotor dysfunction in this condition.

Supplementary Addendum to "Radiation dose of digital radiography (DR) versus micro-dose x-ray (EOS) on patients with adolescent idiopathic scoliosis: 2016 SOSORT- IRSSD "John Sevastic Award" Winner in Imaging Research"

Regarding the publication entitle "Radiation dose of digital radiography (DR) versus micro-dose x-ray (EOS) on patients with adolescent idiopathic scoliosis: 2016 SOSORT- IRSSD "John Sevastic Award" Winner in Imaging Research" in Scoliosis and Spinal Disorders, we would like to provide more details to readers about the dose calculated by simulation using PCXMC 2.0. In this study, the data and results are correct based on the given parameters and calculation provided in the manuscript. In the simulation of EOS micro-dose, only a 0.1 mm copper filter was applied. We agree with the suggestion from Dr. Pedersen and colleagues that the inclusion of a 1.5 mm aluminum filter together with the 0.1 mm copper reflects more realistic representation of X-ray filtration which would improve the accuracy of the simulation. We believe that this supplementary addendum would be beneficial to other researchers who are planning to conduct a similar study.

Automated segmentation of abdominal subcutaneous adipose tissue and visceral adipose tissue in obese adolescent in MRI

PURPOSE

To develop a reliable and reproducible automatic technique to segment and measure SAT and VAT based on MRI.

MATERIALS AND METHODS

Chemical-shift water-fat MRI were taken on twelve obese adolescents (mean age: 16.1±0.6, BMI: 31.3±2.3) recruited under the health monitoring program. The segmentation applied a spoke template created using Midpoint Circle algorithm followed by Bresenham's Line algorithm to detect narrow connecting regions between subcutaneous and visceral adipose tissues. Upon satisfaction of given constrains, a cut was performed to separate SAT and VAT. Bone marrow was consisted in pelvis and femur. By using the intensity difference in T2*, a mask was created to extract bone marrow adipose tissue (MAT) from VAT. Validation was performed using a semi-automatic method. Pearson coefficient, Bland-Altman plot and intra-class coefficient (ICC) were applied to measure accuracy and reproducibility.

RESULTS

Pearson coefficient indicated that results from the proposed method achieved high correlation with the semi-automatic method. Bland-Altman plot and ICC showed good agreement between the two methods. Lowest ICC was obtained in VAT segmentation at lower regions of the abdomen while the rests were all above 0.80. ICC (0.98-0.99) also indicated the proposed method performed good reproducibility.

CONCLUSION

No user interaction was required during execution of the algorithm and the segmented images and volume results were given as output. This technique utilized the feature in the regions connecting subcutaneous and visceral fat and T2* intensity difference in bone marrow to achieve volumetric measurement of various types of adipose tissue in abdominal site.

Upright, prone, and supine spinal morphology and alignment in adolescent idiopathic scoliosis

BACKGROUND

Patients with adolescent idiopathic scoliosis (AIS) are usually investigated by serial imaging studies during the course of treatment, some imaging involves ionizing radiation, and the radiation doses are cumulative. Few studies have addressed the correlation of spinal deformity captured by these different imaging modalities, for which patient positioning are different. To the best of our knowledge, this is the first study to compare the coronal, axial, and sagittal morphology of the scoliotic spine in three different body positions (upright, prone, and supine) and between three different imaging modalities (X-ray, CT, and MRI).

METHODS

Sixty-two AIS patients scheduled for scoliosis surgery, and having undergone standard pre-operative work-up, were included. This work-up included upright full-spine radiographs, supine bending radiographs, supine MRI, and prone CT as is the routine in one of our institutions. In all three positions, Cobb angles, thoracic kyphosis (TK), lumbar lordosis (LL), and vertebral rotation were determined. The relationship among three positions (upright X-ray, prone CT, and supine MRI) was investigated according to the Bland-Altman test, whereas the correlation was described by the intraclass correlation coefficient (ICC).

RESULTS

Thoracic and lumbar Cobb angles correlated significantly between conventional radiographs (68° ± 15° and 44° ± 17°), prone CT (54° ± 15° and 33° ± 15°), and supine MRI (57° ± 14° and 35° ± 16°; ICC ≥0.96; P < 0.001). The thoracic and lumbar apical vertebral rotation showed a good correlation among three positions (upright, 22° ± 12° and 11° ± 13°; prone, 20° ± 9° and 8° ± 11°; supine, 16° ± 11° and 6° ± 14°; ICC ≥0.82; P < 0.001). The TK and LL correlated well among three different positions (TK 26° ± 11°, 22° ± 12°, and 17° ± 10°; P ≤ 0.004; LL 49° ± 12°, 45° ± 11°, and 44° ± 12°; P < 0.006; ICC 0.87 and 0.85).

CONCLUSIONS

Although there is a generalized underestimation of morphological parameters of the scoliotic deformity in the supine and prone positions as compared to the upright position, a significant correlation of these parameters is still evident among different body positions by different imaging modalities. Findings of this study suggest that severity of scoliotic deformity in AIS patients can be largely represented by different imaging modalities despite the difference in body positioning.

Asymmetry of the Vertebral Body and Pedicles in the True Transverse Plane in Adolescent Idiopathic Scoliosis: A CT-Based Study

STUDY DESIGN

Cross-sectional.

OBJECTIVES

To quantify the asymmetry of the vertebral bodies and pedicles in the true transverse plane in adolescent idiopathic scoliosis (AIS) and to compare this with normal anatomy.

SUMMARY OF BACKGROUND DATA

There is an ongoing debate about the existence and magnitude of the vertebral body and pedicle asymmetry in AIS and whether this is an expression of a primary growth disturbance, or secondary to asymmetrical loading.

METHODS

Vertebral body asymmetry, defined as left-right overlap of the vertebral endplates (ie, 100%: perfect symmetry, 0%: complete asymmetry) was evaluated in the true transverse plane on CT scans of 77 AIS patients and 32 non-scoliotic controls. Additionally, the pedicle width, length, and angle and the length of the ideal screw trajectory were calculated.

RESULTS

Scoliotic vertebrae were on average more asymmetric than controls (thoracic: AIS 96.0% vs. controls 96.4%; p = .005, lumbar: 95.8% vs. 97.2%; p < .001) and more pronounced around the thoracic apex (95.8%) than at the end vertebrae (96.3%; p = .031). In the thoracic apex; the concave pedicle was thinner (4.5 vs. 5.4 mm; p < .001) and longer (20.9 vs. 17.9 mm; p < .001), the length of the ideal screw trajectory was longer (43.0 vs. 37.3 mm; p < .001), and the transverse pedicle angle was greater (12.3° vs. 5.7°; p < .001) than the convex one. The axial rotation showed no clear correlation with the asymmetry.

CONCLUSIONS

Even in non-scoliotic controls is a degree of vertebral body and pedicle asymmetry, but scoliotic vertebrae showed slightly more asymmetry, mostly around the thoracic apex. In contrast to the existing literature, there is no major asymmetry in the true transverse plane in AIS and no uniform relation between the axial rotation and vertebral asymmetry could be observed in these moderate to severe patients, suggesting that asymmetrical vertebral growth does not initiate rotation, but rather follows it as a secondary phenomenon.

LEVEL OF EVIDENCE

Level 4.

Radiation dose of digital radiography (DR) versus micro-dose x-ray (EOS) on patients with adolescent idiopathic scoliosis: 2016 SOSORT- IRSSD "John Sevastic Award" Winner in Imaging Research

BACKGROUND

Patients with adolescent idiopathic scoliosis (AIS) frequently receive x-ray imaging at diagnosis and subsequent follow monitoring. The ionizing radiation exposure has accumulated through their development stage and the effect of radiation to this young vulnerable group of patients is uncertain. To achieve the ALARA (as low as reasonably achievable) concept of radiation dose in medical imaging, a slot-scanning x-ray technique by the EOS system has been adopted and the radiation dose using micro-dose protocol was compared with the standard digital radiography on patients with AIS.

METHODS

Ninety-nine participants with AIS underwent micro-dose EOS and 33 underwent standard digital radiography (DR) for imaging of the whole spine. Entrance-skin dose was measured using thermoluminescent dosimeters (TLD) at three regions (i.e. dorsal sites at the level of sternal notch, nipple line, symphysis pubis). Effective dose and organ dose were calculated by simulation using PCXMC 2.0. Data from two x-ray systems were compared using independent-samples t-test and significance level at 0.05. All TLD measurements were conducted on PA projection only. Image quality was also assessed by two raters using Cobb angle measurement and a set of imaging parameters for optimization purposes.

RESULTS

Entrance-skin dose from micro-dose EOS system was 5.9-27.0 times lower at various regions compared with standard DR. The calculated effective dose was 2.6 ± 0.5 (μSv) and 67.5 ± 23.3 (μSv) from micro-dose and standard DR, respectively. The reduction in the micro-dose was approximately 26 times. Organ doses at thyroid, lung and gonad regions were significantly lower in micro-dose (p < 0.001). Data were further compared within the different gender groups. Females received significantly higher (p < 0.001) organ dose at ovaries compared to the testes in males. Patients with AIS received approximately 16-34 times lesser organ dose from micro-dose x-ray as compared with the standard DR. There was no significant difference in overall rating of imaging quality between EOS and DR. Micro-dose protocol provided enough quality to perform consistent measurement on Cobb angle.

CONCLUSIONS

Entrance-skin dose, effective dose and organ dose were significantly reduced in micro-dose x-ray. The effective dose of a single micro-dose x-ray (2.6 μSv) was less than a day of background radiation. As AIS patients require periodic x-ray follow up for surveillance of curve progression, clinical use of micro-dose x-ray system is beneficial for these young patients to reduce the intake of ionizing radiation.

MRI-based morphological evidence of spinal cord tethering predicts curve progression in adolescent idiopathic scoliosis

BACKGROUND CONTEXT

Existing prognostic factors for adolescent idiopathic scoliosis (AIS) patients have focused mainly on curve, maturity, and bone-related factors. Previous studies have shown significant associations between curve severity and morphological evidences of relative shorter spinal cord tethering in AIS, and increased prevalence of abnormal somatosensory cortical-evoked potentials and low-lying cerebellar tonsil in severe AIS. Earlier evidence suggests that there might be neural morphological predictors for curve progression.

PURPOSE

The purpose of this study was to identify any morphological predictors associated with cord tethering, as measured by magnetic resonance imaging (MRI), for curve progression in AIS patients.

STUDY DESIGN/SETTING

This is a prospective cohort study.

PATIENT SAMPLE

A total of 81 female AIS subjects between 10 and 14 years were included, without surgical intervention during the follow-up period.

OUTCOME MEASURES

Magnetic resonance imaging scans of hindbrain and whole spine and areal bone mineral density (BMD) at bilateral femoral necks were performed.

METHODS

All AIS patients were longitudinally followed up starting from initiation of bracing beyond skeletal maturity in 6-month intervals. Clinical and radiographic data were recorded at each clinic visit. Bone mineral density and MRI measurements including ratio of spinal cord to vertebral column length, ratio of anteroposterior (AP) and transverse (TS) diameter of cord, lateral cord space (LCS) ratio, cerebellar tonsil level, and conus medullaris position were obtained at baseline. Only compliant patients with a minimum 2-year follow-up were analyzed. Adolescent idiopathic scoliosis girls were assigned into three groups according to bracing outcome: Group A, nonprogression (curvature increase of less than or equal to 5°); Group B, progression (curvature increase of greater than or equal to 6°); Group C, progression with surgery indication (Cobb angle of greater than or equal to 50° after skeletal maturity despite bracing). The predictors for curve progression were evaluated using univariate analysis and multivariate ordinal regression model.

RESULTS

The average duration of follow-up was 3.4 (range, 2.0-5.6) years. There were 46 girls (57%) in Group A, 19 (23%) in Group B, and 16 (20%) in Group C. No significant intergroup differences were found in spinal cord length, tonsil level, and conus position. Group C had significantly longer vertebral column length, smaller cord-vertebral length ratio, and higher AP/TS cord ratio compared with Group A, whereas LCS ratio in Group C was significantly increased compared with both Group A and Group B. In regression model, five significant independent predictors including cord-vertebral length ratio (odds ratio [OR]: 1.993 [95% confidence interval {CI}: 1.053-3.771], p=.034), LCS ratio (OR: 2.639 [95% CI: 1.128-6.174], p=.025), initial Cobb angle (OR: 1.156 [95% CI: 1.043-1.281], p=.006), menarche age (OR: 1.688 [95% CI: 1.010-2.823], p=.046), and BMD (OR: 2.960 [95% CI: 1.301-6.731], p=.010) and a marginally significant predictor namely AP/TS cord ratio (OR: 1.463 [95% CI: 0.791-2.706], p=.096) were obtained.

CONCLUSIONS

On baseline MRI measurement, cord-vertebral length ratio and LCS ratio are identified as new significant independent predictors for curve progression in AIS, whereas AP/TS cord ratio is suggested as a potential predictor requiring further validations. The earlier MRI parameters can be taken into accounts for prognostication of bracing outcome.

Variation in anisotropy and diffusivity along the medulla oblongata and the whole spinal cord in adolescent idiopathic scoliosis: a pilot study using diffusion tensor imaging

BACKGROUND AND PURPOSE

Disturbed somatosensory evoked potentials have been demonstrated in patients with adolescent idiopathic scoliosis (but this functional delay was found to originate above the C5-6 level, while the lower cord level was unaffected). Together with MR imaging observation of tonsillar ectopia and a relatively tethered cord, we hypothesized that there is disturbed mean diffusivity integrity along the spinal cord. In this study, advanced DTI was used to evaluate whether there was underlying decreased WM integrity within the brain stem and spinal cord in adolescent idiopathic scoliosis and any relationship to cerebellar tonsillar ectopia. Clinical impact on balance testing was also correlated.

MATERIALS AND METHODS

Thirteen girls with adolescent idiopathic scoliosis with right thoracic curves were compared with 13 age-matched healthy girls. DTI of the brain and whole spinal cord was performed. ROIs were manually defined for the medulla oblongata and along each intervertebral segment of the cord. Mean values of fractional anisotropy and mean diffusivity were computed at the defined regions. Between-group comparisons were performed by 1-way ANOVA.

RESULTS

Significantly decreased fractional anisotropy values and increased mean diffusivity values were found at the medulla oblongata and C1-2, C2-3, C3-4, and C4-5 segments in patients with adolescent idiopathic scoliosis compared with healthy subjects. No significant difference was found in the lower cord levels. Significant correlation was found between the tonsillar level and fractional anisotropy value at the C4-5 level in patients with adolescent idiopathic scoliosis only.

CONCLUSIONS

The findings from this study are in agreement with previous findings showing abnormal somatosensory evoked potential readings occurring only above the C5-6 level in patients with adolescent idiopathic scoliosis; these findings might partially explain the pathophysiology of the neural pathway involved.

Application of multimodal MR imaging on studying Alzheimer's disease: a survey

The aim of this review is to summarize the current magnetic resonance imaging (MRI) variants applied on the studies of Alzheimer's disease (AD). Experimental findings, advantages and disadvantages, and the prospect of every individual technique will be presented. MRI can be used to investigate the change of the brain in terms of volume, function, white matter track orientation and even mechanical properties and metabolic concentration.Results from volumetric and morphological analysis indicated that the hippocampus reduced by approximately 18%-24% or 3%-4% annually. Functional MRI (fMRI) detected the blood oxygen level dependent (BOLD) signaland suggested thatmedial temporal lobe (MTL) such as hippocampus and entorhinal cortex reduced at the early stage of AD, and frontal, temporal, and parietal cortices at the later stage. In diffusion tensor imaging (DTI), the white matter fiber track integrity is measured but the results are inconsistent. Besides,magnetic resonance spectroscopy (MRS), chemical shift imaging (CSI), arterial spin labeling (ASL), magnetic resonance elastography (MRE) and susceptibility weighted imaging (SWI) will also be included in this review.

Adaptive distance metric learning for diffusion tensor image segmentation

High quality segmentation of diffusion tensor images (DTI) is of key interest in biomedical research and clinical application. In previous studies, most efforts have been made to construct predefined metrics for different DTI segmentation tasks. These methods require adequate prior knowledge and tuning parameters. To overcome these disadvantages, we proposed to automatically learn an adaptive distance metric by a graph based semi-supervised learning model for DTI segmentation. An original discriminative distance vector was first formulated by combining both geometry and orientation distances derived from diffusion tensors. The kernel metric over the original distance and labels of all voxels were then simultaneously optimized in a graph based semi-supervised learning approach. Finally, the optimization task was efficiently solved with an iterative gradient descent method to achieve the optimal solution. With our approach, an adaptive distance metric could be available for each specific segmentation task. Experiments on synthetic and real brain DTI datasets were performed to demonstrate the effectiveness and robustness of the proposed distance metric learning approach. The performance of our approach was compared with three classical metrics in the graph based semi-supervised learning framework.

Altered topological organization of cortical network in adolescent girls with idiopathic scoliosis

Adolescent idiopathic scoliosis (AIS) is a multifactorial disease affecting approximately 1-4% of teenagers especially girls at the age of 10-16, but its etiopathogenesis remains uncertain. Previous study has revealed that the cortical thickness in AIS patients is different from that in normal controls. Cortical thickness measurements are known to be strongly correlated between regions that are axonally connected. Hence, a hypothesis is proposed to study the possibility to demonstrate abnormal structural network revealed by cortical thickness in AIS patients. The aim of the study is to investigate abnormalities in the organization of the brain cortical network in AIS patients. This study included 42 girls with severe idiopathic scoliosis (14.7±1.3 years old) and 41 age-matched normal controls (NC, 14.6±1.4 years old). The brain cortex was partitioned into 154 cortical regions based on gyral and sulcal structure. The interregional connectivity was measured as the statistical correlations between the regional mean thicknesses across the subjects. We employed the graph theoretic analysis to examine the alteration in interregional correlation, small-world efficiency, hub distribution, and regional nodal characteristics in AIS patients. We demonstrated that the cortical network of AIS patients fully preserved the small-world architecture and organization, and further verified the hemispheric asymmetry of AIS brain. Our results indicated increased central role of temporal and occipital cortex and decreased central role of limbic cortex in AIS patients compared with controls. Furthermore, decreased structural connectivity between hemispheres and increased connectivity in several cortical regions were observed. The findings of the study reveal the pattern of structural network alteration in AIS brain, and would help in understanding the mechanism and etiopathogenesis of AIS.