Comparison of quantitative susceptibility mapping methods on evaluating radiation-induced cerebral microbleeds and basal ganglia at 3T and 7T.

Quantitative susceptibility mapping (QSM) has the potential for being a biomarker for various diseases because of its ability to measure tissue susceptibility related to iron deposition, myelin, and hemorrhage from the phase signal of a T2 *-weighted MRI. Despite its promise as a quantitative marker, QSM is faced with many challenges, including its dependence on preprocessing of the raw phase data, the relatively weak tissue signal, and the inherently ill posed relationship between the magnetic dipole and measured phase. The goal of this study was to evaluate the effects of background field removal and dipole inversion algorithms on noise characteristics, image uniformity, and structural contrast for cerebral microbleed (CMB) quantification at both 3T and 7T. We selected four widely used background phase removal and five dipole field inversion algorithms for QSM and applied them to volunteers and patients with CMBs, who were scanned at two different field strengths, with ground truth QSM reference calculated using multiple orientation scans. 7T MRI provided QSM images with lower noise than did 3T MRI. QSIP and VSHARP + iLSQR achieved the highest white matter homogeneity and vein contrast, with QSIP also providing the highest CMB contrast. Compared with ground truth COSMOS QSM images, overall good correlations between susceptibility values of dipole inversion algorithms and the COSMOS reference were observed in basal ganglia regions, with VSHARP + iLSQR achieving the susceptibility values most similar to COSMOS across all regions. This study can provide guidance for selecting the most appropriate QSM processing pipeline based on the application of interest and scanner field strength.

On the accuracy and reproducibility of a novel probabilistic atlas-based generation for calculation of head attenuation maps on integrated PET/MR scanners.

PURPOSE: To propose an MR-based method for generating continuous-valued head attenuation maps and to assess its accuracy and reproducibility. Demonstrating that novel MR-based photon attenuation correction methods are both accurate and reproducible is essential prior to using them routinely in research and clinical studies on integrated PET/MR scanners. METHODS: Continuous-valued linear attenuation coefficient maps ("µ-maps") were generated by combining atlases that provided the prior probability of voxel positions belonging to a certain tissue class (air, soft tissue, or bone) and an MR intensity-based likelihood classifier to produce posterior probability maps of tissue classes. These probabilities were used as weights to generate the µ-maps. The accuracy of this probabilistic atlas-based continuous-valued µ-map ("PAC-map") generation method was assessed by calculating the voxel-wise absolute relative change (RC) between the MR-based and scaled CT-based attenuation-corrected PET images. To assess reproducibility, we performed pair-wise comparisons of the RC values obtained from the PET images reconstructed using the µ-maps generated from the data acquired at three time points. RESULTS: The proposed method produced continuous-valued µ-maps that qualitatively reflected the variable anatomy in patients with brain tumor and agreed well with the scaled CT-based µ-maps. The absolute RC comparing the resulting PET volumes was 1.76 ± 2.33 %, quantitatively demonstrating that the method is accurate. Additionally, we also showed that the method is highly reproducible, the mean RC value for the PET images reconstructed using the µ-maps obtained at the three visits being 0.65 ± 0.95 %. CONCLUSION: Accurate and highly reproducible continuous-valued head µ-maps can be generated from MR data using a probabilistic atlas-based approach.

Intravoxel incoherent motion analysis of renal allograft diffusion with clinical and histopathological correlation in pediatric kidney transplant patients: A preliminary cross-sectional observational study.

The purpose of this study was to compare IVIM values in pediatric renal transplants with histopathology and clinical management change. Fifteen pediatric renal transplant recipients (mean 15.7±2.9 years) were prospectively scanned on a 3T MR scanner with multi-b DTI, prior to same-day transplant biopsy. IVIM maps from 14 subjects were analyzed (one excluded due to motion). Mean values were computed from cortical ROIs and medullary ROIs corresponding to the biopsy site. Subjects were also grouped according to whether or not the biopsy resulted in a change in clinical management. Cortico-medullary IVIM estimates and histopathologic Banff scores were correlated with KT. Cortico-medullary IVIM differences between the "change" and "no change" groups was compared with Mann-Whitney U test. Cortical Dp showed significant moderate negative correlation with Banff t and ci scores (KT=-0.497, P=.035 and KT=-0.46, P=.046) and moderate positive correlation with Banff i score (KT=0.527, P=.028). Cortical Pf showed significant moderate correlation with ci and ct scores (KT=0.489, P=.035 and KT=0.457, P=.043). Tissue diffusivity, Dt , estimated with IVIM was significantly different between the "change" and "no change" groups in medullary ROIs (U=6, P=.021). IVIM analysis has potential as a noninvasive biomarker in assessment of pediatric renal allograft pathology.

Anatomy-specific Progression Classification in Chest Radiographs via Weakly Supervised Learning.

Purpose To develop a machine learning approach for classifying disease progression in chest radiographs using weak labels automatically derived from radiology reports. Materials and Methods In this retrospective study, a twin neural network was developed to classify anatomy-specific disease progression into four categories: improved, unchanged, worsened, and new. A two-step weakly supervised learning approach was employed, pretraining the model on 243 008 frontal chest radiographs from 63 877 patients (mean age, 51.7 years ± 17.0 [SD]; 34 813 [55%] female) included in the MIMIC-CXR database and fine-tuning it on the subset with progression labels derived from consecutive studies. Model performance was evaluated for six pathologic observations on test datasets of unseen patients from the MIMIC-CXR database. Area under the receiver operating characteristic (AUC) analysis was used to evaluate classification performance. The algorithm is also capable of generating bounding-box predictions to localize areas of new progression. Recall, precision, and mean average precision were used to evaluate the new progression localization. One-tailed paired t tests were used to assess statistical significance. Results The model outperformed most baselines in progression classification, achieving macro AUC scores of 0.72 ± 0.004 for atelectasis, 0.75 ± 0.007 for consolidation, 0.76 ± 0.017 for edema, 0.81 ± 0.006 for effusion, 0.7 ± 0.032 for pneumonia, and 0.69 ± 0.01 for pneumothorax. For new observation localization, the model achieved mean average precision scores of 0.25 ± 0.03 for atelectasis, 0.34 ± 0.03 for consolidation, 0.33 ± 0.03 for edema, and 0.31 ± 0.03 for pneumothorax. Conclusion Disease progression classification models were developed on a large chest radiograph dataset, which can be used to monitor interval changes and detect new pathologic conditions on chest radiographs. Keywords: Prognosis, Unsupervised Learning, Transfer Learning, Convolutional Neural Network (CNN), Emergency Radiology, Named Entity Recognition Supplemental material is available for this article. © RSNA, 2024 See also commentary by Alves and Venkadesh in this issue.

Quantitative Interstitial Abnormality Progression and Outcomes in the Genetic Epidemiology of COPD and Pittsburgh Lung Screening Study Cohorts.

BACKGROUND: The risk factors and clinical outcomes of quantitative interstitial abnormality progression over time have not been characterized. RESEARCH QUESTIONS: What are the associations of quantitative interstitial abnormality progression with lung function, exercise capacity, and mortality? What are the demographic and genetic risk factors for quantitative interstitial abnormality progression? STUDY DESIGN AND METHODS: Quantitative interstitial abnormality progression between visits 1 and 2 was assessed from 4,635 participants in the Genetic Epidemiology of COPD (COPDGene) cohort and 1,307 participants in the Pittsburgh Lung Screening Study (PLuSS) cohort. We used multivariable linear regression to determine the risk factors for progression and the longitudinal associations between progression and FVC and 6-min walk distance, and Cox regression models for the association with mortality. RESULTS: Age at enrollment, female sex, current smoking status, and the MUC5B minor allele were associated with quantitative interstitial abnormality progression. Each percent annual increase in quantitative interstitial abnormalities was associated with annual declines in FVC (COPDGene: 8.5 mL/y; 95% CI, 4.7-12.4 mL/y; P < .001; PLuSS: 9.5 mL/y; 95% CI, 3.7-15.4 mL/y; P = .001) and 6-min walk distance, and increased mortality (COPDGene: hazard ratio, 1.69; 95% CI, 1.34-2.12; P < .001; PLuSS: hazard ratio, 1.28; 95% CI, 1.10-1.49; P = .001). INTERPRETATION: The objective, longitudinal measurement of quantitative interstitial abnormalities may help identify people at greatest risk for adverse events and most likely to benefit from early intervention.

Quantitative susceptibility mapping by inversion of a perturbation field model: correlation with brain iron in normal aging.

There is increasing evidence that iron deposition occurs in specific regions of the brain in normal aging and neurodegenerative disorders such as Parkinson's, Huntington's, and Alzheimer's disease. Iron deposition changes the magnetic susceptibility of tissue, which alters the MR signal phase, and allows estimation of susceptibility differences using quantitative susceptibility mapping (QSM). We present a method for quantifying susceptibility by inversion of a perturbation model, or "QSIP." The perturbation model relates phase to susceptibility using a kernel calculated in the spatial domain, in contrast to previous Fourier-based techniques. A tissue/air susceptibility atlas is used to estimate B0 inhomogeneity. QSIP estimates in young and elderly subjects are compared to postmortem iron estimates, maps of the Field-Dependent Relaxation Rate Increase, and the L1-QSM method. Results for both groups showed excellent agreement with published postmortem data and in vivo FDRI: statistically significant Spearman correlations ranging from Rho=0.905 to Rho=1.00 were obtained. QSIP also showed improvement over FDRI and L1-QSM: reduced variance in susceptibility estimates and statistically significant group differences were detected in striatal and brainstem nuclei, consistent with age-dependent iron accumulation in these regions.

Probabilistic atlas-based segmentation of combined T1-weighted and DUTE MRI for calculation of head attenuation maps in integrated PET/MRI scanners.

We present a new MRI-based attenuation correction (AC) approach for integrated PET/MRI systems that combines both segmentation- and atlas-based methods by incorporating dual-echo ultra-short echo-time (DUTE) and T1-weighted (T1w) MRI data and a probabilistic atlas. Segmented atlases were constructed from CT training data using a leave-one-out framework and combined with T1w, DUTE, and CT data to train a classifier that computes the probability of air/soft tissue/bone at each voxel. This classifier was applied to segment the MRI of the subject of interest and attenuation maps (µ-maps) were generated by assigning specific linear attenuation coefficients (LACs) to each tissue class. The µ-maps generated with this "Atlas-T1w-DUTE" approach were compared to those obtained from DUTE data using a previously proposed method. For validation of the segmentation results, segmented CT µ-maps were considered to the "silver standard"; the segmentation accuracy was assessed qualitatively and quantitatively through calculation of the Dice similarity coefficient (DSC). Relative change (RC) maps between the CT and MRI-based attenuation corrected PET volumes were also calculated for a global voxel-wise assessment of the reconstruction results. The µ-maps obtained using the Atlas-T1w-DUTE classifier agreed well with those derived from CT; the mean DSCs for the Atlas-T1w-DUTE-based µ-maps across all subjects were higher than those for DUTE-based µ-maps; the atlas-based µ-maps also showed a lower percentage of misclassified voxels across all subjects. RC maps from the atlas-based technique also demonstrated improvement in the PET data compared to the DUTE method, both globally as well as regionally.

Morphometric differences in planum temporale in schizophrenia and bipolar disorder revealed by statistical analysis of labeled cortical depth maps.

Differences in cortical thickness in the lateral temporal lobe, including the planum temporale (PT), have been reported in MRI studies of schizophrenia (SCZ) and bipolar disorder (BPD) patients. Most of these studies have used a single-valued global or local measure for thickness. However, additional and complementary information can be obtained by generating labeled cortical distance maps (LCDMs), which are distances of labeled gray matter (GM) voxels from the nearest point on the GM/white matter (WM) (inner) cortical surface. Statistical analyses of pooled and censored LCDM distances reveal subtle differences in PT between SCZ and BPD groups from data generated by Ratnanather et al. (Schizophrenia Research, http://dx.doi.org/10.1016/j.schres.2013.08.014). These results confirm that the left planum temporale (LPT) is more sensitive than the right PT in distinguishing between SCZ, BPD, and healthy controls. Also confirmed is a strong gender effect, with a thicker PT seen in males than in females. The differences between groups at smaller distances in the LPT revealed by pooled and censored LCDM analysis suggest that SCZ and BPD have different effects on the cortical mantle close to the GM/WM surface. This is consistent with reported subtle changes in the cortical mantle observed in post-mortem studies.

Morphometry of superior temporal gyrus and planum temporale in schizophrenia and psychotic bipolar disorder.

Structural abnormalities in temporal lobe, including the superior temporal gyrus (STG) and planum temporale (PT), have been reported in schizophrenia (SCZ) and bipolar disorder (BPD) patients. While most MRI studies have suggested gray matter volume and surface area reduction in temporal lobe regions, few have explored changes in laminar thickness in PT and STG in SCZ and BPD. ROI subvolumes of the STG from 94 subjects were used to yield gray matter volume, gray/white surface area and laminar thickness for STG and PT cortical regions. Morphometric analysis suggests that there may be gender and laterality effects on the size and shape of the PT in BPD (n=36) and SCZ (n=31) with reduced laterality in PT in subjects with SCZ but not in BPD. In addition, PT surface area was seen to be larger in males, and asymmetry in PT surface area was larger in BPD. Subjects with SCZ had reduced thickness and smaller asymmetry in PT volume. Thus, the PT probably plays a more sensitive role than the STG in structural abnormalities seen in SCZ.

Cortico-cortical, cortico-striatal, and cortico-thalamic white matter fiber tracts generated in the macaque brain via dynamic programming.

Probabilistic methods have the potential to generate multiple and complex white matter fiber tracts in diffusion tensor imaging (DTI). Here, a method based on dynamic programming (DP) is introduced to reconstruct fibers pathways whose complex anatomical structures cannot be resolved beyond the resolution of standard DTI data. DP is based on optimizing a sequentially additive cost function derived from a Gaussian diffusion model whose covariance is defined by the diffusion tensor. DP is used to determine the optimal path between initial and terminal nodes by efficiently searching over all paths, connecting the nodes, and choosing the path in which the total probability is maximized. An ex vivo high-resolution scan of a macaque hemi-brain is used to demonstrate the advantages and limitations of DP. DP can generate fiber bundles between distant cortical areas (superior longitudinal fasciculi, arcuate fasciculus, uncinate fasciculus, and fronto-occipital fasciculus), neighboring cortical areas (dorsal and ventral banks of the principal sulcus), as well as cortical projections to the hippocampal formation (cingulum bundle), neostriatum (motor cortical projections to the putamen), thalamus (subcortical bundle), and hippocampal formation projections to the mammillary bodies via the fornix. Validation is established either by comparison with in vivo intracellular transport of horseradish peroxidase in another macaque monkey or by comparison with atlases. DP is able to generate known pathways, including crossing and kissing tracts. Thus, DP has the potential to enhance neuroimaging studies of cortical connectivity.

Atlas-based improved prediction of magnetic field inhomogeneity for distortion correction of EPI data.

We describe a method for atlas-based segmentation of structural MRI for calculation of magnetic fieldmaps. CT data sets are used to construct a probabilistic atlas of the head and corresponding MR is used to train a classifier that segments soft tissue, air, and bone. Subject-specific fieldmaps are computed from the segmentations using a perturbation field model. Previous work has shown that distortion in echo-planar images can be corrected using predicted fieldmaps. We obtain results that agree well with acquired fieldmaps: 90% of voxel shifts from predicted fieldmaps show subvoxel disagreement with those computed from acquired fieldmaps. In addition, our fieldmap predictions show statistically significant improvement following inclusion of the atlas.

Fieldmap-free retrospective registration and distortion correction for EPI-based functional imaging.

We describe a method for correcting the distortions present in echo planar images (EPI) and registering the EPI to structural MRI. A fieldmap is predicted from an air / tissue segmentation of the MRI using a perturbation method and subsequently used to unwarp the EPI data. Shim and other missing parameters are estimated by registration. We obtain results that are similar to those obtained using fieldmaps, however neither fieldmaps, nor knowledge of shim coefficients is required.

tGolgin-1 (p230, golgin-245) modulates Shiga-toxin transport to the Golgi and Golgi motility towards the microtubule-organizing centre.

tGolgin-1 (trans-Golgi p230, golgin-245) is a member of a family of large peripheral membrane proteins that associate with the trans-Golgi network (TGN) via a C-terminal GRIP domain. Some GRIP-domain proteins have been implicated in endosome-to-TGN transport but no function for tGolgin-1 has been described. Here, we show that tGolgin-1 production is required for efficient retrograde distribution of Shiga toxin from endosomes to the Golgi. Surprisingly, we also found an indirect requirement for tGolgin-1 in Golgi positioning. In HeLa cells depleted of tGolgin-1, the normally centralized Golgi and TGN membranes were displaced to the periphery, forming 'mini stacks'. These stacks resembled those in cells with disrupted microtubules or dynein-dynactin motor, in that they localized to endoplasmic-reticulum exit sites, maintained their secretory capacity and cis-trans polarity, and were relatively immobile by video microscopy. The mini stacks formed concomitant with a failure of pre-Golgi elements to migrate along microtubules towards the microtubule-organizing centre. The requirement for tGolgin-1 in Golgi positioning did not appear to reflect direct binding of tGolgin-1 to motile pre-Golgi membranes, because distinct Golgi and tGolgin-1-containing TGN elements that formed after recovery of HeLa cells from brefeldin-A treatment moved independently toward the microtubule-organizing centre. These data demonstrate that tGolgin-1 functions in Golgi positioning indirectly, probably by regulating retrograde movement of cargo required for recruitment or activation of dynein-dynactin complexes on newly formed Golgi elements.