Fast refacing of MR images with a generative neural network lowers re-identification risk and preserves volumetric consistency.

With the rise of open data, identifiability of individuals based on 3D renderings obtained from routine structural magnetic resonance imaging (MRI) scans of the head has become a growing privacy concern. To protect subject privacy, several algorithms have been developed to de-identify imaging data using blurring, defacing or refacing. Completely removing facial structures provides the best re-identification protection but can significantly impact post-processing steps, like brain morphometry. As an alternative, refacing methods that replace individual facial structures with generic templates have a lower effect on the geometry and intensity distribution of original scans, and are able to provide more consistent post-processing results by the price of higher re-identification risk and computational complexity. In the current study, we propose a novel method for anonymized face generation for defaced 3D T1-weighted scans based on a 3D conditional generative adversarial network. To evaluate the performance of the proposed de-identification tool, a comparative study was conducted between several existing defacing and refacing tools, with two different segmentation algorithms (FAST and Morphobox). The aim was to evaluate (i) impact on brain morphometry reproducibility, (ii) re-identification risk, (iii) balance between (i) and (ii), and (iv) the processing time. The proposed method takes 9 s for face generation and is suitable for recovering consistent post-processing results after defacing.

White adipose tissue distribution and amount are associated with increased white matter connectivity.

Obesity represents a significant public health concern and is linked to various comorbidities and cognitive impairments. Previous research indicates that elevated body mass index (BMI) is associated with structural changes in white matter (WM). However, a deeper characterization of body composition is required, especially considering the links between abdominal obesity and metabolic dysfunction. This study aims to enhance our understanding of the relationship between obesity and WM connectivity by directly assessing the amount and distribution of fat tissue. Whole-body magnetic resonance imaging (MRI) was employed to evaluate total adipose tissue (TAT), visceral adipose tissue (VAT), and subcutaneous adipose tissue (SAT), while MR liver spectroscopy measured liver fat content in 63 normal-weight, overweight, and obese males. WM connectivity was quantified using microstructure-informed tractography. Connectome-based predictive modeling was used to predict body composition metrics based on WM connectomes. Our analysis revealed a positive dependency between BMI, TAT, SAT, and WM connectivity in brain regions involved in reward processing and appetite regulation, such as the insula, nucleus accumbens, and orbitofrontal cortex. Increased connectivity was also observed in cognitive control and inhibition networks, including the middle frontal gyrus and anterior cingulate cortex. No significant associations were found between WM connectivity and VAT or liver fat. Our findings suggest that altered neural communication between these brain regions may affect cognitive processes, emotional regulation, and reward perception in individuals with obesity, potentially contributing to weight gain. While our study did not identify a link between WM connectivity and VAT or liver fat, further investigation of the role of various fat depots and metabolic factors in brain networks is required to advance obesity prevention and treatment approaches.

Pore size estimation in axon-mimicking microfibers with diffusion-relaxation MRI.

PURPOSE: This study aims to evaluate two distinct approaches for fiber radius estimation using diffusion-relaxation MRI data acquired in biomimetic microfiber phantoms that mimic hollow axons. The methods considered are the spherical mean power-law approach and a T2-based pore size estimation technique. THEORY AND METHODS: A general diffusion-relaxation theoretical model for the spherical mean signal from water molecules within a distribution of cylinders with varying radii was introduced, encompassing the evaluated models as particular cases. Additionally, a new numerical approach was presented for estimating effective radii (i.e., MRI-visible mean radii) from the ground truth radii distributions, not reliant on previous theoretical approximations and adaptable to various acquisition sequences. The ground truth radii were obtained from scanning electron microscope images. RESULTS: Both methods show a linear relationship between effective radii estimated from MRI data and ground-truth radii distributions, although some discrepancies were observed. The spherical mean power-law method overestimated fiber radii. Conversely, the T2-based method exhibited higher sensitivity to smaller fiber radii, but faced limitations in accurately estimating the radius in one particular phantom, possibly because of material-specific relaxation changes. CONCLUSION: The study demonstrates the feasibility of both techniques to predict pore sizes of hollow microfibers. The T2-based technique, unlike the spherical mean power-law method, does not demand ultra-high diffusion gradients, but requires calibration with known radius distributions. This research contributes to the ongoing development and evaluation of neuroimaging techniques for fiber radius estimation, highlights the advantages and limitations of both methods, and provides datasets for reproducible research.

Cortical alterations after very preterm birth and the association with socio-emotional abilities from childhood to early adolescence.

Very preterm birth (VPT; <32 weeks' gestation) leads to a situation where crucial steps of brain development occur in an abnormal ex utero environment, translating to vulnerable cortical and subcortical development. Associated with this atypical brain development, children and adolescents born VPT are at a high risk of socio-emotional difficulties. In the current study, we unravel developmental changes in cortical gray matter (GM) concentration in VPT and term-born controls aged 6-14 years, together with their associations with socio-emotional abilities. T1-weighted images were used to estimate signal intensities of brain tissue types in a single voxel (GM, white matter, and cortico-spinal fluid) and extract GM concentration disentangled from the presence of partial volume effects (PVEs). General linear model analysis was used to compare groups. Socio-emotional abilities were assessed and associations with GM concentration were explored using univariate and multivariate analyses. The effects of prematurity were far-reaching, with intricated patterns of increases and decreases of GM concentration mainly in frontal, temporal, parietal, and cingular regions. Better socio-emotional abilities were associated with increased GM concentration in regions known to be involved in such process for both groups. Our findings suggest that the trajectory of brain development following VPT birth may be fundamentally distinctive and impact socio-emotional abilities.

Acceptability of artificial intelligence for cervical cancer screening in Dschang, Cameroon: a qualitative study on patient perspectives.

BACKGROUND: Cervical cancer is the fourth most frequent cancer among women, with 90% of cervical cancer-related deaths occurring in low- and middle-income countries like Cameroon. Visual inspection with acetic acid is often used in low-resource settings to screen for cervical cancer; however, its accuracy can be limited. To address this issue, the Swiss Federal Institute of Technology Lausanne and the University Hospitals of Geneva are collaborating to develop an automated smartphone-based image classifier that serves as a computer aided diagnosis tool for cancerous lesions. The primary objective of this study is to explore the acceptability and perspectives of women in Dschang regarding the usage of a screening tool for cervical cancer relying on artificial intelligence. A secondary objective is to understand the preferred form and type of information women would like to receive regarding this artificial intelligence-based screening tool. METHODS: A qualitative methodology was employed to gain better insight into the women's perspectives. Participants, aged between 30 and 49 were invited from both rural and urban regions and semi-structured interviews using a pre-tested interview guide were conducted. The focus groups were divided on the basis of level of education, as well as HPV status. The interviews were audio-recorded, transcribed, and coded using the ATLAS.ti software. RESULTS: A total of 32 participants took part in the six focus groups, and 38% of participants had a primary level of education. The perspectives identified were classified using an adapted version of the Technology Acceptance Model. Key factors influencing the acceptability of artificial intelligence include privacy concerns, perceived usefulness, and trust in the competence of providers, accuracy of the tool as well as the potential negative impact of smartphones. CONCLUSION: The results suggest that an artificial intelligence-based screening tool for cervical cancer is mostly acceptable to the women in Dschang. By ensuring patient confidentiality and by providing clear explanations, acceptance can be fostered in the community and uptake of cervical cancer screening can be improved. TRIAL REGISTRATION: Ethical Cantonal Board of Geneva, Switzerland (CCER, N°2017-0110 and CER-amendment n°4) and Cameroonian National Ethics Committee for Human Health Research (N°2022/12/1518/CE/CNERSH/SP). NCT: 03757299.

Globally, cervical cancer is the fourth most frequent cancer among women. However, 90% of all deaths caused by cervical cancer occur in low-and middle-income countries. Methods traditionally used in settings like Cameroon to detect cervical cancer unfortunately lack accuracy. Therefore, researchers at the Swiss Federal Institute of Technology Lausanne and the University Hospitals of Geneva are developing an artificial intelligence-based computer aided diagnosis tool to detect pre-cancerous lesions using a smartphone application. The aim of this study was to explore the acceptability and perspectives regarding an AI-based tool for cervical cancer screening for women in Dschang, a city in the west of Cameroon. A qualitative methodology was conducted with six focus groups and a total of 32 participants. The main concerns highlighted by the study are related to privacy, trust in the ability of the healthcare providers, accuracy of the tool as well as the potential negative impact of smartphones. In conclusion, our results show that a computer aided diagnosis tool using artificial intelligence is mostly acceptable to women in Dschang, as long as their confidentiality is preserved, and they are provided with clear explanations beforehand.

Computer-Assisted Diagnosis of Lymph Node Metastases in Colorectal Cancers Using Transfer Learning With an Ensemble Model.

Screening of lymph node metastases in colorectal cancer (CRC) can be a cumbersome task, but it is amenable to artificial intelligence (AI)-assisted diagnostic solution. Here, we propose a deep learning-based workflow for the evaluation of CRC lymph node metastases from digitized hematoxylin and eosin-stained sections. A segmentation model was trained on 100 whole-slide images (WSIs). It achieved a Matthews correlation coefficient of 0.86 (±0.154) and an acceptable Hausdorff distance of 135.59 µm (±72.14 µm), indicating a high congruence with the ground truth. For metastasis detection, 2 models (Xception and Vision Transformer) were independently trained first on a patch-based breast cancer lymph node data set and were then fine-tuned using the CRC data set. After fine-tuning, the ensemble model showed significant improvements in the F1 score (0.797-0.949; P <.00001) and the area under the receiver operating characteristic curve (0.959-0.978; P <.00001). Four independent cohorts (3 internal and 1 external) of CRC lymph nodes were used for validation in cascading segmentation and metastasis detection models. Our approach showed excellent performance, with high sensitivity (0.995, 1.0) and specificity (0.967, 1.0) in 2 validation cohorts of adenocarcinoma cases (n = 3836 slides) when comparing slide-level labels with the ground truth (pathologist reports). Similarly, an acceptable performance was achieved in a validation cohort (n = 172 slides) with mucinous and signet-ring cell histology (sensitivity, 0.872; specificity, 0.936). The patch-based classification confidence was aggregated to overlay the potential metastatic regions within each lymph node slide for visualization. We also applied our method to a consecutive case series of lymph nodes obtained over the past 6 months at our institution (n = 217 slides). The overlays of prediction within lymph node regions matched 100% when compared with a microscope evaluation by an expert pathologist. Our results provide the basis for a computer-assisted diagnostic tool for easy and efficient lymph node screening in patients with CRC.

Cellular Exchange Imaging (CEXI): Evaluation of a diffusion model including water exchange in cells using numerical phantoms of permeable spheres.

PURPOSE: Biophysical models of diffusion MRI have been developed to characterize microstructure in various tissues, but existing models are not suitable for tissue composed of permeable spherical cells. In this study we introduce Cellular Exchange Imaging (CEXI), a model tailored for permeable spherical cells, and compares its performance to a related Ball & Sphere (BS) model that neglects permeability. METHODS: We generated DW-MRI signals using Monte-Carlo simulations with a PGSE sequence in numerical substrates made of spherical cells and their extracellular space for a range of membrane permeability. From these signals, the properties of the substrates were inferred using both BS and CEXI models. RESULTS: CEXI outperformed the impermeable model by providing more stable estimates cell size and intracellular volume fraction that were diffusion time-independent. Notably, CEXI accurately estimated the exchange time for low to moderate permeability levels previously reported in other studies ( κ < 25 µ m / s $$ \kappa <25\kern0.3em \mu \mathrm{m}/\mathrm{s} $$ ). However, in highly permeable substrates ( κ = 50 µ m / s $$ \kappa =50\kern0.3em \mu \mathrm{m}/\mathrm{s} $$ ), the estimated parameters were less stable, particularly the diffusion coefficients. CONCLUSION: This study highlights the importance of modeling the exchange time to accurately quantify microstructure properties in permeable cellular substrates. Future studies should evaluate CEXI in clinical applications such as lymph nodes, investigate exchange time as a potential biomarker of tumor severity, and develop more appropriate tissue models that account for anisotropic diffusion and highly permeable membranes.

Submillimeter T1 atlas for subject-specific abnormality detection at 7T.

PURPOSE: Studies at 3T have shown that T1 relaxometry enables characterization of brain tissues at the single-subject level by comparing individual physical properties to a normative atlas. In this work, an atlas of normative T1 values at 7T is introduced with 0.6 mm isotropic resolution and its clinical potential is explored in comparison to 3T. METHODS: T1 maps were acquired in two separate healthy cohorts scanned at 3T and 7T. Using transfer learning, a template-based brain segmentation algorithm was adapted to ultra-high field imaging data. After segmenting brain tissues, volumes were normalized into a common space, and an atlas of normative T1 values was established by modeling the T1 inter-subject variability. A method for single-subject comparisons restricted to white matter and subcortical structures was developed by computing Z-scores. The comparison was applied to eight patients scanned at both field strengths for proof of concept. RESULTS: The proposed method for morphometry delivered segmentation masks without statistically significant differences from those derived with the original pipeline at 3T and achieved accurate segmentation at 7T. The established normative atlas allowed characterizing tissue alterations in single-subject comparisons at 7T, and showed greater anatomical details compared with 3T results. CONCLUSION: A high-resolution quantitative atlas with an adapted pipeline was introduced and validated. Several case studies on different clinical conditions showed the feasibility, potential and limitations of high-resolution single-subject comparisons based on quantitative MRI atlases. This method in conjunction with 7T higher resolution broadens the range of potential applications of quantitative MRI in clinical practice.

Examining the synergistic effects of a cognitive control video game and a home-based, self-administered non-invasive brain stimulation on alleviating depression: the DiSCoVeR trial protocol.

Enhanced behavioral interventions are gaining increasing interest as innovative treatment strategies for major depressive disorder (MDD). In this study protocol, we propose to examine the synergistic effects of a self-administered home-treatment, encompassing transcranial direct current stimulation (tDCS) along with a video game based training of attentional control. The study is designed as a two-arm, double-blind, randomized and placebo-controlled multi-center trial (ClinicalTrials.gov: NCT04953208). At three study sites (Israel, Latvia, and Germany), 114 patients with a primary diagnosis of MDD undergo 6 weeks of intervention (30 × 30 min sessions). Patients assigned to the intervention group receive active tDCS (anode F3 and cathode F4; 2 mA intensity) and an action-like video game, while those assigned to the control group receive sham tDCS along with a control video game. An electrode-positioning algorithm is used to standardize tDCS electrode positioning. Participants perform their designated treatment at the clinical center (sessions 1-5) and continue treatment at home under remote supervision (sessions 6-30). The endpoints are feasibility (primary) and safety, treatment efficacy (secondary, i.e., change of Montgomery-Åsberg Depression Rating Scale (MADRS) scores at week six from baseline, clinical response and remission, measures of social, occupational, and psychological functioning, quality of life, and cognitive control (tertiary). Demonstrating the feasibility, safety, and efficacy of this novel combined intervention could expand the range of available treatments for MDD to neuromodulation enhanced interventions providing cost-effective, easily accessible, and low-risk treatment options.ClinicalTrials.gov: NCT04953208.

Axon morphology is modulated by the local environment and impacts the noninvasive investigation of its structure-function relationship.

Axonal conduction velocity, which ensures efficient function of the brain network, is related to axon diameter. Noninvasive, in vivo axon diameter estimates can be made with diffusion magnetic resonance imaging, but the technique requires three-dimensional (3D) validation. Here, high-resolution, 3D synchrotron X-ray nano-holotomography images of white matter samples from the corpus callosum of a monkey brain reveal that blood vessels, cells, and vacuoles affect axonal diameter and trajectory. Within single axons, we find that the variation in diameter and conduction velocity correlates with the mean diameter, contesting the value of precise diameter determination in larger axons. These complex 3D axon morphologies drive previously reported 2D trends in axon diameter and g-ratio. Furthermore, we find that these morphologies bias the estimates of axon diameter with diffusion magnetic resonance imaging and, ultimately, impact the investigation and formulation of the axon structure-function relationship.

Evaluating reproducibility and subject-specificity of microstructure-informed connectivity.

Tractography enables identifying and evaluating the healthy and diseased brain's white matter pathways from diffusion-weighted magnetic resonance imaging data. As previous evaluation studies have reported significant false-positive estimation biases, recent microstructure-informed tractography algorithms have been introduced to improve the trade-off between specificity and sensitivity. However, a major limitation for characterizing the performance of these techniques is the lack of ground truth brain data. In this study, we compared the performance of two relevant microstructure-informed tractography methods, SIFT2 and COMMIT, by assessing the subject specificity and reproducibility of their derived white matter pathways. Specifically, twenty healthy young subjects were scanned at eight different time points at two different sites. Subject specificity and reproducibility were evaluated using the whole-brain connectomes and a subset of 29 white matter bundles. Our results indicate that although the raw tractograms are more vulnerable to the presence of false-positive connections, they are highly reproducible, suggesting that the estimation bias is subject-specific. This high reproducibility was preserved when microstructure-informed tractography algorithms were used to filter the raw tractograms. Moreover, the resulting track-density images depicted a more uniform coverage of streamlines throughout the white matter, suggesting that these techniques could increase the biological meaning of the estimated fascicles. Notably, we observed an increased subject specificity by employing connectivity pre-processing techniques to reduce the underlaying noise and the data dimensionality (using principal component analysis), highlighting the importance of these tools for future studies. Finally, no strong bias from the scanner site or time between measurements was found. The largest intraindividual variance originated from the sole repetition of data measurements (inter-run).

Insights from the IronTract challenge: Optimal methods for mapping brain pathways from multi-shell diffusion MRI.

Limitations in the accuracy of brain pathways reconstructed by diffusion MRI (dMRI) tractography have received considerable attention. While the technical advances spearheaded by the Human Connectome Project (HCP) led to significant improvements in dMRI data quality, it remains unclear how these data should be analyzed to maximize tractography accuracy. Over a period of two years, we have engaged the dMRI community in the IronTract Challenge, which aims to answer this question by leveraging a unique dataset. Macaque brains that have received both tracer injections and ex vivo dMRI at high spatial and angular resolution allow a comprehensive, quantitative assessment of tractography accuracy on state-of-the-art dMRI acquisition schemes. We find that, when analysis methods are carefully optimized, the HCP scheme can achieve similar accuracy as a more time-consuming, Cartesian-grid scheme. Importantly, we show that simple pre- and post-processing strategies can improve the accuracy and robustness of many tractography methods. Finally, we find that fiber configurations that go beyond crossing (e.g., fanning, branching) are the most challenging for tractography. The IronTract Challenge remains open and we hope that it can serve as a valuable validation tool for both users and developers of dMRI analysis methods.

Data-driven myelin water imaging based on T1 and T2 relaxometry.

Long acquisition times preclude the application of multiecho spin echo (MESE) sequences for myelin water fraction (MWF) mapping in daily clinical practice. In search of alternative methods, previous studies of interest explored the biophysical modeling of MWF from measurements of different tissue properties that can be obtained in scan times shorter than those required for the MESE. In this work, a novel data-driven estimation of MWF maps from fast relaxometry measurements is proposed and investigated. T1 and T2 relaxometry maps were acquired in a cohort of 20 healthy subjects along with a conventional MESE sequence. Whole-brain quantitative mapping was achieved with a fast protocol in 6 min 24 s. Reference MWF maps were derived from the MESE sequence (TA = 11 min 17 s) and their data-driven estimation from relaxometry measurements was investigated using three different modeling strategies: two general linear models (GLMs) with linear and quadratic regressors, respectively; a random forest regression model; and two deep neural network architectures, a U-Net and a conditional generative adversarial network (cGAN). Models were validated using a 10-fold crossvalidation. The resulting maps were visually and quantitatively compared by computing the root mean squared error (RMSE) between the estimated and reference MWF maps, the intraclass correlation coefficients (ICCs) between corresponding MWF values in different brain regions, and by performing Bland-Altman analysis. Qualitatively, the estimated maps appear to generally provide a similar, yet more blurred MWF contrast in comparison with the reference, with the cGAN model best capturing MWF variabilities in small structures. By estimating the average adjusted coefficient of determination of the GLM with quadratic regressors, we showed that 87% of the variability in the MWF values can be explained by relaxation times alone. Further quantitative analysis showed an average RMSE smaller than 0.1% for all methods. The ICC was greater than 0.81 for all methods, and the bias smaller than 2.19%. It was concluded that this work confirms the notion that relaxometry parameters contain a large part of the information on myelin water and that MWF maps can be generated from T1 /T2 data with minimal error. Among the investigated modeling approaches, the cGAN provided maps with the best trade-off between accuracy and blurriness. Fast relaxometry, like the 6 min 24 s whole-brain protocol used in this work in conjunction with machine learning, may thus have the potential to replace time-consuming MESE acquisitions.

Multiple sclerosis cortical lesion detection with deep learning at ultra-high-field MRI.

Manually segmenting multiple sclerosis (MS) cortical lesions (CLs) is extremely time consuming, and past studies have shown only moderate inter-rater reliability. To accelerate this task, we developed a deep-learning-based framework (CLAIMS: Cortical Lesion AI-Based Assessment in Multiple Sclerosis) for the automated detection and classification of MS CLs with 7 T MRI. Two 7 T datasets, acquired at different sites, were considered. The first consisted of 60 scans that include 0.5 mm isotropic MP2RAGE acquired four times (MP2RAGE×4), 0.7 mm MP2RAGE, 0.5 mm T2 *-weighted GRE, and 0.5 mm T2 *-weighted EPI. The second dataset consisted of 20 scans including only 0.75 × 0.75 × 0.9 mm3 MP2RAGE. CLAIMS was first evaluated using sixfold cross-validation with single and multi-contrast 0.5 mm MRI input. Second, the performance of the model was tested on 0.7 mm MP2RAGE images after training with either 0.5 mm MP2RAGE×4, 0.7 mm MP2RAGE, or alternating the two. Third, its generalizability was evaluated on the second external dataset and compared with a state-of-the-art technique based on partial volume estimation and topological constraints (MSLAST). CLAIMS trained only with MP2RAGE×4 achieved results comparable to those of the multi-contrast model, reaching a CL true positive rate of 74% with a false positive rate of 30%. Detection rate was excellent for leukocortical and subpial lesions (83%, and 70%, respectively), whereas it reached 53% for intracortical lesions. The correlation between disability measures and CL count was similar for manual and CLAIMS lesion counts. Applying a domain-scanner adaptation approach and testing CLAIMS on the second dataset, the performance was superior to MSLAST when considering a minimum lesion volume of 6 µL (lesion-wise detection rate of 71% versus 48%). The proposed framework outperforms previous state-of-the-art methods for automated CL detection across scanners and protocols. In the future, CLAIMS may be useful to support clinical decisions at 7 T MRI, especially in the field of diagnosis and differential diagnosis of MS patients.

The structural connectome and motor recovery after stroke: predicting natural recovery.

Stroke patients vary considerably in terms of outcomes: some patients present 'natural' recovery proportional to their initial impairment (fitters), while others do not (non-fitters). Thus, a key challenge in stroke rehabilitation is to identify individual recovery potential to make personalized decisions for neuro-rehabilitation, obviating the 'one-size-fits-all' approach. This goal requires (i) the prediction of individual courses of recovery in the acute stage; and (ii) an understanding of underlying neuronal network mechanisms. 'Natural' recovery is especially variable in severely impaired patients, underscoring the special clinical importance of prediction for this subgroup. Fractional anisotropy connectomes based on individual tractography of 92 patients were analysed 2 weeks after stroke (TA) and their changes to 3 months after stroke (TC - TA). Motor impairment was assessed using the Fugl-Meyer Upper Extremity (FMUE) scale. Support vector machine classifiers were trained to separate patients with natural recovery from patients without natural recovery based on their whole-brain structural connectomes and to define their respective underlying network patterns, focusing on severely impaired patients (FMUE < 20). Prediction accuracies were cross-validated internally, in one independent dataset and generalized in two independent datasets. The initial connectome 2 weeks after stroke was capable of segregating fitters from non-fitters, most importantly among severely impaired patients (TA: accuracy = 0.92, precision = 0.93). Secondary analyses studying recovery-relevant network characteristics based on the selected features revealed (i) relevant differences between networks contributing to recovery at 2 weeks and network changes over time (TC - TA); and (ii) network properties specific to severely impaired patients. Important features included the parietofrontal motor network including the intraparietal sulcus, premotor and primary motor cortices and beyond them also attentional, somatosensory or multimodal areas (e.g. the insula), strongly underscoring the importance of whole-brain connectome analyses for better predicting and understanding recovery from stroke. Computational approaches based on structural connectomes allowed the individual prediction of natural recovery 2 weeks after stroke onset, especially in the difficult to predict group of severely impaired patients, and identified the relevant underlying neuronal networks. This information will permit patients to be stratified into different recovery groups in clinical settings and will pave the way towards personalized precision neurorehabilitative treatment.

Revisiting the T2 spectrum imaging inverse problem: Bayesian regularized non-negative least squares.

Multi-echo T2 magnetic resonance images contain information about the distribution of T2 relaxation times of compartmentalized water, from which we can estimate relevant brain tissue properties such as the myelin water fraction (MWF). Regularized non-negative least squares (NNLS) is the tool of choice for estimating non-parametric T2 spectra. However, the estimation is ill-conditioned, sensitive to noise, and highly affected by the employed regularization weight. The purpose of this study is threefold: first, we want to underline that the apparently innocuous use of two alternative parameterizations for solving the inverse problem, which we called the standard and alternative regularization forms, leads to different solutions; second, to assess the performance of both parameterizations; and third, to propose a new Bayesian regularized NNLS method (BayesReg). The performance of BayesReg was compared with that of two conventional approaches (L-curve and Chi-square (X2) fitting) using both regularization forms. We generated a large dataset of synthetic data, acquired in vivo human brain data in healthy participants for conducting a scan-rescan analysis, and correlated the myelin content derived from histology with the MWF estimated from ex vivo data. Results from synthetic data indicate that BayesReg provides accurate MWF estimates, comparable to those from L-curve and X2, and with better overall stability across a wider signal-to-noise range. Notably, we obtained superior results by using the alternative regularization form. The correlations reported in this study are higher than those reported in previous studies employing the same ex vivo and histological data. In human brain data, the estimated maps from L-curve and BayesReg were more reproducible. However, the T2 spectra produced by BayesReg were less affected by over-smoothing than those from L-curve. These findings suggest that BayesReg is a good alternative for estimating T2 distributions and MWF maps.

Resolving bundle-specific intra-axonal T2 values within a voxel using diffusion-relaxation tract-based estimation.

At the typical spatial resolution of MRI in the human brain, approximately 60-90% of voxels contain multiple fiber populations. Quantifying microstructural properties of distinct fiber populations within a voxel is therefore challenging but necessary. While progress has been made for diffusion and T1-relaxation properties, how to resolve intra-voxel T2 heterogeneity remains an open question. Here a novel framework, named COMMIT-T2, is proposed that uses tractography-based spatial regularization with diffusion-relaxometry data to estimate multiple intra-axonal T2 values within a voxel. Unlike previously-proposed voxel-based T2 estimation methods, which (when applied in white matter) implicitly assume just one fiber bundle in the voxel or the same T2 for all bundles in the voxel, COMMIT-T2 can recover specific T2 values for each unique fiber population passing through the voxel. In this approach, the number of recovered unique T2 values is not determined by a number of model parameters set a priori, but rather by the number of tractography-reconstructed streamlines passing through the voxel. Proof-of-concept is provided in silico and in vivo, including a demonstration that distinct tract-specific T2 profiles can be recovered even in the three-way crossing of the corpus callosum, arcuate fasciculus, and corticospinal tract. We demonstrate the favourable performance of COMMIT-T2 compared to that of voxelwise approaches for mapping intra-axonal T2 exploiting diffusion, including a direction-averaged method and AMICO-T2, a new extension to the previously-proposed Accelerated Microstructure Imaging via Convex Optimization (AMICO) framework.

Tractography dissection variability: What happens when 42 groups dissect 14 white matter bundles on the same dataset?

White matter bundle segmentation using diffusion MRI fiber tractography has become the method of choice to identify white matter fiber pathways in vivo in human brains. However, like other analyses of complex data, there is considerable variability in segmentation protocols and techniques. This can result in different reconstructions of the same intended white matter pathways, which directly affects tractography results, quantification, and interpretation. In this study, we aim to evaluate and quantify the variability that arises from different protocols for bundle segmentation. Through an open call to users of fiber tractography, including anatomists, clinicians, and algorithm developers, 42 independent teams were given processed sets of human whole-brain streamlines and asked to segment 14 white matter fascicles on six subjects. In total, we received 57 different bundle segmentation protocols, which enabled detailed volume-based and streamline-based analyses of agreement and disagreement among protocols for each fiber pathway. Results show that even when given the exact same sets of underlying streamlines, the variability across protocols for bundle segmentation is greater than all other sources of variability in the virtual dissection process, including variability within protocols and variability across subjects. In order to foster the use of tractography bundle dissection in routine clinical settings, and as a fundamental analytical tool, future endeavors must aim to resolve and reduce this heterogeneity. Although external validation is needed to verify the anatomical accuracy of bundle dissections, reducing heterogeneity is a step towards reproducible research and may be achieved through the use of standard nomenclature and definitions of white matter bundles and well-chosen constraints and decisions in the dissection process.

Probing myelin content of the human brain with MRI: A review.

Rapid and efficient transmission of electric signals among neurons of vertebrates is ensured by myelin-insulating sheaths surrounding axons. Human cognition, sensation, and motor functions rely on the integrity of these layers, and demyelinating diseases often entail serious cognitive and physical impairments. Magnetic resonance imaging radically transformed the way these disorders are monitored, offering an irreplaceable tool to noninvasively examine the brain structure. Several advanced techniques based on MRI have been developed to provide myelin-specific contrasts and a quantitative estimation of myelin density in vivo. Here, the vast offer of acquisition strategies developed to date for this task is reviewed. Advantages and pitfalls of the different approaches are compared and discussed.

Fast and high-resolution myelin water imaging: Accelerating multi-echo GRASE with CAIPIRINHA.

PURPOSE: Although several MRI methods have been explored to achieve in vivo myelin quantification, imaging the whole brain in clinically acceptable times and sufficiently high resolution remains challenging. To address this problem, this work investigates the acceleration of multi-echo T2 acquisitions based on the multi-echo gradient and spin echo (GRASE) sequence using CAIPIRINHA undersampling and adapted k-space reordering patterns. METHODS: A prototype multi-echo GRASE sequence supporting CAIPIRINHA parallel imaging was implemented. Multi-echo T2 data were acquired from 12 volunteers using the implemented sequence (1.6 × 1.6 × 1.6 mm3 , 84 slices, acquisition time [TA] = 10:30 min) and a multi-echo spin echo (MESE) sequence as reference (1.6 × 1.6 × 3.2 mm3 , single-slice, TA = 5:41 min). Myelin water fraction (MWF) maps derived from both acquisitions were compared via correlation and Bland-Altman analyses. In addition, scan-rescan datasets were acquired to evaluate the repeatability of the derived maps. RESULTS: Resulting maps from the MESE and multi-echo GRASE sequences were found to be correlated (r = 0.83). The Bland-Altman analysis revealed a mean bias of -0.2% (P = .24) with the limits of agreement ranging from -3.7% to 3.3%. The Pearson's correlation coefficient among MWF values obtained from the scan-rescan datasets was found to be 0.95 and the mean bias equal to 0.11% (P = .32), indicating good repeatability of the retrieved maps. CONCLUSION: By combining a 3D multi-echo GRASE sequence with CAIPIRINHA sampling, whole-brain MWF maps were obtained in 10:30 min with 1.6 mm isotropic resolution. The good correlation with conventional MESE-based maps demonstrates that the implemented sequence may be a promising alternative to time-consuming MESE acquisitions.