Highlights

• Spinal cord morphology estimated based on synthetic T₁-weighted MRI shows high repeatability.

• Spinal cord atrophy following spinal cord injury can reliably be measured on synthetic MR images.

• Synthetic MRI allows to reduce acquisition time in long imaging protocols and application in clinical studies.

Abstract

Short MRI acquisition time, high signal-to-noise ratio, and high reliability are crucial for image quality when scanning healthy volunteers and patients. Cross-sectional cervical cord area (CSA) has been suggested as a marker of neurodegeneration and potential outcome measure in clinical trials and is conventionally measured on T₁-weigthed 3D Magnetization Prepared Rapid Acquisition Gradient-Echo (MPRAGE) images. This study aims to reduce the acquisition time for the comprehensive assessment of the spinal cord, which is typically based on MPRAGE for morphometry and multi-parameter mapping (MPM) for microstructure. The MPRAGE is replaced by a synthetic T₁-w MRI (synT₁-w) estimated from the MPM, in order to measure CSA. SynT₁-w images were reconstructed using the MPRAGE signal equation based on quantitative maps of proton density (PD), longitudinal (R₁) and effective transverse (R₂*) relaxation rates. The reliability of CSA measurements from synT₁-w images was determined within a multi-center test-retest study format and validated against acquired MPRAGE scans by assessing the agreement between both methods. The response to pathological changes was tested by longitudinally measuring spinal cord atrophy following spinal cord injury (SCI) for synT₁-w and MPRAGE using linear mixed effect models. CSA measurements based on the synT₁-w MRI showed high intra-site (Coefficient of variation [CoV]: 1.43% to 2.71%) and inter-site repeatability (CoV: 2.90% to 5.76%), and only a minor deviation of -1.65 mm² compared to MPRAGE. Crucially, by assessing atrophy rates and by comparing SCI patients with healthy controls longitudinally, differences between synT₁-w and MPRAGE were negligible. These results demonstrate that reliable estimates of CSA can be obtained from synT₁-w images, thereby reducing scan time significantly.

Highlights

• We proposed a data augmentation approach CarveMix for brain lesion segmentation.

• CarveMix mixes pairs of annotated images to generate synthetic training images.

• The image mixing is performed according to the location and shape of the lesions.

• CarveMix was validated on multiple public and private datasets.

• The results show that CarveMix improves the quality of brain lesion segmentation.

Abstract

Brain lesion segmentation provides a valuable tool for clinical diagnosis and research, and convolutional neural networks (CNNs) have achieved unprecedented success in the segmentation task. Data augmentation is a widely used strategy to improve the training of CNNs. In particular, data augmentation approaches that mix pairs of annotated training images have been developed. These methods are easy to implement and have achieved promising results in various image processing tasks. However, existing data augmentation approaches based on image mixing are not designed for brain lesions and may not perform well for brain lesion segmentation. Thus, the design of this type of simple data augmentation method for brain lesion segmentation is still an open problem. In this work, we propose a simple yet effective data augmentation approach, dubbed as CarveMix, for CNN-based brain lesion segmentation. Like other mixing-based methods, CarveMix stochastically combines two existing annotated images (annotated for brain lesions only) to obtain new labeled samples. To make our method more suitable for brain lesion segmentation, CarveMix is lesion-aware, where the image combination is performed with a focus on the lesions and preserves the lesion information. Specifically, from one annotated image we carve a region of interest (ROI) according to the lesion location and geometry with a variable ROI size. The carved ROI then replaces the corresponding voxels in a second annotated image to synthesize new labeled images for network training, and additional harmonization steps are applied for heterogeneous data where the two annotated images can originate from different sources. Besides, we further propose to model the mass effect that is unique to whole brain tumor segmentation during image mixing. To evaluate the proposed method, experiments were performed on multiple publicly available or private datasets, and the results show that our method improves the accuracy of brain lesion segmentation. The code of the proposed method is available at https://github.com/ZhangxinruBIT/CarveMix.git.

Highlights

• 3D velocity selective inversion based arterial spin labeling (VSI-ASL) achieved high reliability during the test-retest CBF measurement for detecting between-subject variations among healthy subjects using absolute CBF.

• 3D VSI-ASL achieved high reliability during the test-retest CBF measurement for detecting between-region variations among healthy subjects using relative CBF.

• Physiology parameters partially explained the variability of CBF across subjects.

Abstract

Velocity-selective inversion (VSI) based velocity-selective arterial spin labeling (VSASL) has been developed to measure cerebral blood flow (CBF) with low susceptibility to the prolonged arterial transit time and high sensitivity to brain perfusion signal. The purpose of this magnetic resonance imaging study is to evaluate the test-retest reliability of a VSI-prepared 3D VSASL protocol with whole-brain coverage to detect baseline CBF variations among cognitively normal participants in different brain regions. Coefficients of variation (CoV) of both absolute and relative CBF across scans or sessions, subjects, and gray matter regions were calculated, and corresponding intraclass correlation coefficients (ICC) were computed. The higher between-subject CoV of absolute CBF (13.4 ± 2.0%) over within-subject CoV (within-session: 3.8 ± 1.1%; between-session: 4.9 ± 0.9%) yielded moderate to excellent ICC (within-session: 0.88±0.08; between-session: 0.77±0.14) to detect normal variations of individual CBF. The higher between-region CoV of relative CBF (11.4 ± 3.0%) over within-region CoV (within-session: 2.3 ± 0.9%; between-session: 3.3 ± 1.0%) yielded excellent ICC (within-session: 0.92±0.06; between-session: 0.85±0.12) to detect normal variations of regional CBF. Age, blood pressure, end-tidal CO₂, and hematocrit partially explained the variability of CBF across subjects. Together these results show excellent test-retest reliability of VSASL to detect both between-subject and between-region variations supporting its clinical utility.

HIGHLIGHTS

• ModelArray is an R package for statistical analysis of fixel-wise data.

• ModelArray supports linear and nonlinear modeling and is extensible to more models.

• ModelArray is scalable for large-scale datasets.

• ModelArray facilitates easy statistical analysis of large-scale fixel-wise data.

ABSTRACT

Diffusion MRI is the dominant non-invasive imaging method used to characterize white matter organization in health and disease. Increasingly, fiber-specific properties within a voxel are analyzed using fixels. While tools for conducting statistical analyses of fixel-wise data exist, currently available tools support only a limited number of statistical models. Here we introduce ModelArray, an R package for mass-univariate statistical analysis of fixel-wise data. At present, ModelArray supports linear models as well as generalized additive models (GAMs), which are particularly useful for studying nonlinear effects in lifespan data. In addition, ModelArray also aims for scalable analysis. With only several lines of code, even large fixel-wise datasets can be analyzed using a standard personal computer. Detailed memory profiling revealed that ModelArray required only limited memory even for large datasets. As an example, we applied ModelArray to fixel-wise data derived from diffusion images acquired as part of the Philadelphia Neurodevelopmental Cohort (n = 938). ModelArray revealed anticipated nonlinear developmental effects in white matter. Moving forward, ModelArray is supported by an open-source software development model that can incorporate additional statistical models and other imaging data types. Taken together, ModelArray provides a flexible and efficient platform for statistical analysis of fixel-wise data.

Graphical abstract

Highlights

• Assessment of inter-organ metabolic connectivity.

• Automated and manual organ delineation.

• Validation of metabolic connectivity approach.

• Liver and kidney strongest connectivity with brain.

Abstract

The nervous and circulatory system interconnects the various organs of the human body, building hierarchically organized subsystems, enabling fine-tuned, metabolically expensive brain-body and inter-organ crosstalk to appropriately adapt to internal and external demands. A deviation or failure in the function of a single organ or subsystem could trigger unforeseen biases or dysfunctions of the entire network, leading to maladaptive physiological or psychological responses. Therefore, quantifying these networks in healthy individuals and patients may help further our understanding of complex disorders involving body-brain crosstalk.

Here we present a generalized framework to automatically estimate metabolic inter-organ connectivity utilizing whole-body functional positron emission tomography (fPET). The developed framework was applied to 16 healthy subjects (mean age ± SD, 25 ± 6 years; 13 female) that underwent one dynamic ¹⁸F-FDG PET/CT scan. Multiple procedures of organ segmentation (manual, automatic, circular volumes) and connectivity estimation (polynomial fitting, spatiotemporal filtering, covariance matrices) were compared to provide an optimized thorough overview of the workflow.

The proposed approach was able to estimate the metabolic connectivity patterns within brain regions and organs as well as their interactions. Automated organ delineation, but not simplified circular volumes, showed high agreement with manual delineation. Polynomial fitting yielded similar connectivity as spatiotemporal filtering at the individual subject level. Furthermore, connectivity measures and group-level covariance matrices did not match. The strongest brain-body connectivity was observed for the liver and kidneys.

The proposed framework offers novel opportunities towards analyzing metabolic function from a systemic, hierarchical perspective in a multitude of physiological pathological states.