Continuous evaluation of denoising strategies in resting-state fMRI connectivity using fMRIPrep and Nilearn.

Reducing contributions from non-neuronal sources is a crucial step in functional magnetic resonance imaging (fMRI) connectivity analyses. Many viable strategies for denoising fMRI are used in the literature, and practitioners rely on denoising benchmarks for guidance in the selection of an appropriate choice for their study. However, fMRI denoising software is an ever-evolving field, and the benchmarks can quickly become obsolete as the techniques or implementations change. In this work, we present a denoising benchmark featuring a range of denoising strategies, datasets and evaluation metrics for connectivity analyses, based on the popular fMRIprep software. The benchmark prototypes an implementation of a reproducible framework, where the provided Jupyter Book enables readers to reproduce or modify the figures on the Neurolibre reproducible preprint server (https://neurolibre.org/). We demonstrate how such a reproducible benchmark can be used for continuous evaluation of research software, by comparing two versions of the fMRIprep. Most of the benchmark results were consistent with prior literature. Scrubbing, a technique which excludes time points with excessive motion, combined with global signal regression, is generally effective at noise removal. Scrubbing was generally effective, but is incompatible with statistical analyses requiring the continuous sampling of brain signal, for which a simpler strategy, using motion parameters, average activity in select brain compartments, and global signal regression, is preferred. Importantly, we found that certain denoising strategies behave inconsistently across datasets and/or versions of fMRIPrep, or had a different behavior than in previously published benchmarks. This work will hopefully provide useful guidelines for the fMRIprep users community, and highlight the importance of continuous evaluation of research methods.

Impact of B 0 $$ {\mathrm{B}}_0 $$ field imperfections correction on BOLD sensitivity in 3D-SPARKLING fMRI data.

PURPOSE: Static and dynamic B 0 $$ {\mathrm{B}}_0 $$ field imperfections are detrimental to functional MRI (fMRI) applications, especially at ultra-high magnetic fields (UHF). In this work, a field camera is used to assess the benefits of retrospectively correcting B 0 $$ {\mathrm{B}}_0 $$ field perturbations on Blood Oxygen Level Dependent (BOLD) sensitivity in non-Cartesian three-dimensional (3D)-SPARKLING fMRI acquisitions. METHODS: fMRI data were acquired at 1 mm 3 $$ {}^3 $$ and for a 2.4s-TR while concurrently monitoring in real-time field perturbations using a Skope Clip-on field camera in a novel experimental setting involving a shorter TR than the required minimal TR of the field probes. Measurements of the dynamic field deviations were used along with a static Δ B 0 $$ \Delta {\mathrm{B}}_0 $$ map to retrospectively correct static and dynamic field imperfections, respectively. In order to evaluate the impact of such a correction on fMRI volumes, a comparative study was conducted on healthy volunteers. RESULTS: Correction of B 0 $$ {\mathrm{B}}_0 $$ deviations improved image quality and yielded between 20% and 30% increase in median temporal signal-to-noise ratio (tSNR).Using fMRI data collected during a retinotopic mapping experiment, we demonstrated a significant increase in sensitivity to the BOLD contrast and improved accuracy of the BOLD phase maps: 44% (resp., 159%) more activated voxels were retrieved when using a significance control level based on a p-value of 0.001 without correcting for multiple comparisons (resp., 0.05 with a false discovery rate correction). CONCLUSION: 3D-SPARKLING fMRI hugely benefits from static and dynamic B 0 $$ {\mathrm{B}}_0 $$ imperfections correction. However, the proposed experimental protocol is flexible enough to be deployed on a large spectrum of encoding schemes, including arbitrary non-Cartesian readouts.

Verbal fluency predicts work resumption after awake surgery in low-grade glioma patients.

BACKGROUND: Resuming professional activity after awake surgery for diffuse low-grade glioma (DLGG) is an important goal, which is not reached in every patient. Cognitive deficits can occur and persist after surgery. In this study, we analyzed the impact of mild cognitive impairments on the work resumption. METHODS: Fifty-four surgeries (including five redo surgeries) performed between 2012 and 2020 for grade 2 (45) and 3 (nine) DLGG in 49 professionally active patients (mean age 40 [range 23-58.) were included. We retrospectively extracted the results of semantic and phonemic verbal fluency tests from preoperative and 4-month postoperative cognitive assessments. Patients were interviewed about their working life after surgery, between April and June 2021. RESULTS: Patients (85%) returned to work, most within 3 to 6 months. Patients (76%) reported subjective complaints (primarily fatigue). Self-reported symptoms and individual and clinical variables had no impact on the work resumption. Late-postoperative average Z-scores in verbal fluency tasks were significantly lower than preoperative for the entire cohort (Wilcoxon test, p < 0.001 for semantic and p = 0.008 for phonemic fluency). The decrease in Z-scores was significantly greater (Mann Whitney U-test, semantic, p = 0.018; phonemic, p = 0.004) in the group of patients who did not return to work than in the group of patients who did. CONCLUSION: The proportion of patients returning to work was comparable to similar studies. A decrease in verbal fluency tasks could predict the inability to return to work.

Notip: Non-parametric true discovery proportion control for brain imaging.

Cluster-level inference procedures are widely used for brain mapping. These methods compare the size of clusters obtained by thresholding brain maps to an upper bound under the global null hypothesis, computed using Random Field Theory or permutations. However, the guarantees obtained by this type of inference - i.e. at least one voxel is truly activated in the cluster - are not informative with regards to the strength of the signal therein. There is thus a need for methods to assess the amount of signal within clusters; yet such methods have to take into account that clusters are defined based on the data, which creates circularity in the inference scheme. This has motivated the use of post hoc estimates that allow statistically valid estimation of the proportion of activated voxels in clusters. In the context of fMRI data, the All-Resolutions Inference framework introduced in Rosenblatt et al. (2018) provides post hoc estimates of the proportion of activated voxels. However, this method relies on parametric threshold families, which results in conservative inference. In this paper, we leverage randomization methods to adapt to data characteristics and obtain tighter false discovery control. We obtain Notip, for Non-parametric True Discovery Proportion control: a powerful, non-parametric method that yields statistically valid guarantees on the proportion of activated voxels in data-derived clusters. Numerical experiments demonstrate substantial gains in number of detections compared with state-of-the-art methods on 36 fMRI datasets. The conditions under which the proposed method brings benefits are also discussed.

Synthetic FLAIR as a Substitute for FLAIR Sequence in Acute Ischemic Stroke.

Background In acute ischemic stroke (AIS), fluid-attenuated inversion recovery (FLAIR) is used for treatment decisions when onset time is unknown. Synthetic FLAIR could be generated with deep learning from information embedded in diffusion-weighted imaging (DWI) and could replace acquired FLAIR sequence (real FLAIR) and shorten MRI duration. Purpose To compare performance of synthetic and real FLAIR for DWI-FLAIR mismatch estimation and identification of patients presenting within 4.5 hours from symptom onset. Materials and Methods In this retrospective study, all pretreatment and early follow-up (<48 hours after symptom onset) MRI data sets including DWI (b = 0-1000 sec/mm2) and FLAIR sequences obtained in consecutive patients with AIS referred for reperfusion therapies between January 2002 and May 2019 were included. On the training set (80%), a generative adversarial network was trained to produce synthetic FLAIR with DWI as input. On the test set (20%), synthetic FLAIR was computed without real FLAIR knowledge. The DWI-FLAIR mismatch was evaluated on both FLAIR data sets by four independent readers. Interobserver reproducibility and DWI-FLAIR mismatch concordance between synthetic and real FLAIR were evaluated with κ statistics. Sensitivity and specificity for identification of AIS within 4.5 hours were compared in patients with known onset time by using McNemar test. Results The study included 1416 MRI scans (861 patients; median age, 71 years [interquartile range, 57-81 years]; 375 men), yielding 1134 and 282 scans for training and test sets, respectively. Regarding DWI-FLAIR mismatch, interobserver reproducibility was substantial for real and synthetic FLAIR (κ = 0.80 [95% CI: 0.74, 0.87] and 0.80 [95% CI: 0.74, 0.87], respectively). After consensus, concordance between real and synthetic FLAIR was almost perfect (κ = 0.88; 95% CI: 0.82, 0.93). Diagnostic value for identifying AIS within 4.5 hours did not differ between real and synthetic FLAIR (sensitivity: 107 of 131 [82%] vs 111 of 131 [85%], P = .2; specificity: 96 of 104 [92%] vs 96 of 104 [92%], respectively, P > .99). Conclusion Synthetic fluid-attenuated inversion recovery (FLAIR) had diagnostic performances similar to real FLAIR in depicting diffusion-weighted imaging-FLAIR mismatch and in helping to identify early acute ischemic stroke, and it may accelerate MRI protocols. © RSNA, 2022 Online supplemental material is available for this article. See also the editorial by Carroll and Hurley in this issue.

Extracting representations of cognition across neuroimaging studies improves brain decoding.

Cognitive brain imaging is accumulating datasets about the neural substrate of many different mental processes. Yet, most studies are based on few subjects and have low statistical power. Analyzing data across studies could bring more statistical power; yet the current brain-imaging analytic framework cannot be used at scale as it requires casting all cognitive tasks in a unified theoretical framework. We introduce a new methodology to analyze brain responses across tasks without a joint model of the psychological processes. The method boosts statistical power in small studies with specific cognitive focus by analyzing them jointly with large studies that probe less focal mental processes. Our approach improves decoding performance for 80% of 35 widely-different functional-imaging studies. It finds commonalities across tasks in a data-driven way, via common brain representations that predict mental processes. These are brain networks tuned to psychological manipulations. They outline interpretable and plausible brain structures. The extracted networks have been made available; they can be readily reused in new neuro-imaging studies. We provide a multi-study decoding tool to adapt to new data.

Tackling the Complexity of Lesion-Symptoms Mapping: How to Bridge the Gap Between Data Scientists and Clinicians?

Accurate and predictive lesion-symptoms mapping is a major goal in the field of clinical neurosciences. Recent studies have called for a reappraisal of the results given by the standard univariate voxel-based lesion-symptom mapping technique, emphasizing the need of developing multivariate methods. While the organization of large datasets and their analysis with machine learning (ML) approaches represents an opportunity to increase prediction accuracy, the complexity and dimensionality of the problem remain a major obstacle. Acknowledging the difficulty of inferring individual outcomes from the observation of spatial patterns of lesions, we propose here to base prediction on new individuals on models of brain connectivity, whereby the disruption of a given network predicts the occurrence of selective deficits. Well-suited ML tools are necessary to capture the relevant information from limited datasets and perform reliable inference.

Functional annotation of human cognitive states using deep graph convolution.

A key goal in neuroscience is to understand brain mechanisms of cognitive functions. An emerging approach is "brain decoding", which consists of inferring a set of experimental conditions performed by a participant, using pattern classification of brain activity. Few works so far have attempted to train a brain decoding model that would generalize across many different cognitive tasks drawn from multiple cognitive domains. To tackle this problem, we proposed a multidomain brain decoder that automatically learns the spatiotemporal dynamics of brain response within a short time window using a deep learning approach. We evaluated the decoding model on a large population of 1200 participants, under 21 different experimental conditions spanning six different cognitive domains, acquired from the Human Connectome Project task-fMRI database. Using a 10s window of fMRI response, the 21 cognitive states were identified with a test accuracy of 90% (chance level 4.8%). Performance remained good when using a 6s window (82%). It was even feasible to decode cognitive states from a single fMRI volume (720ms), with the performance following the shape of the hemodynamic response. Moreover, a saliency map analysis demonstrated that the high decoding performance was driven by the response of biologically meaningful brain regions. Together, we provide an automated tool to annotate human brain activity with fine temporal resolution and fine cognitive granularity. Our model shows potential applications as a reference model for domain adaptation, possibly making contributions in a variety of domains, including neurological and psychiatric disorders.

Decoding with confidence: Statistical control on decoder maps.

In brain imaging, decoding is widely used to infer relationships between brain and cognition, or to craft brain-imaging biomarkers of pathologies. Yet, standard decoding procedures do not come with statistical guarantees, and thus do not give confidence bounds to interpret the pattern maps that they produce. Indeed, in whole-brain decoding settings, the number of explanatory variables is much greater than the number of samples, hence classical statistical inference methodology cannot be applied. Specifically, the standard practice that consists in thresholding decoding maps is not a correct inference procedure. We contribute a new statistical-testing framework for this type of inference. To overcome the statistical inefficiency of voxel-level control, we generalize the Family Wise Error Rate (FWER) to account for a spatial tolerance δ, introducing the δ-Family Wise Error Rate (δ-FWER). Then, we present a decoding procedure that can control the δ-FWER: the Ensemble of Clustered Desparsified Lasso (EnCluDL), a procedure for multivariate statistical inference on high-dimensional structured data. We evaluate the statistical properties of EnCluDL with a thorough empirical study, along with three alternative procedures including decoder map thresholding. We show that EnCluDL exhibits the best recovery properties while ensuring the expected statistical control.

An empirical evaluation of functional alignment using inter-subject decoding.

Inter-individual variability in the functional organization of the brain presents a major obstacle to identifying generalizable neural coding principles. Functional alignment-a class of methods that matches subjects' neural signals based on their functional similarity-is a promising strategy for addressing this variability. To date, however, a range of functional alignment methods have been proposed and their relative performance is still unclear. In this work, we benchmark five functional alignment methods for inter-subject decoding on four publicly available datasets. Specifically, we consider three existing methods: piecewise Procrustes, searchlight Procrustes, and piecewise Optimal Transport. We also introduce and benchmark two new extensions of functional alignment methods: piecewise Shared Response Modelling (SRM), and intra-subject alignment. We find that functional alignment generally improves inter-subject decoding accuracy though the best performing method depends on the research context. Specifically, SRM and Optimal Transport perform well at both the region-of-interest level of analysis as well as at the whole-brain scale when aggregated through a piecewise scheme. We also benchmark the computational efficiency of each of the surveyed methods, providing insight into their usability and scalability. Taking inter-subject decoding accuracy as a quantification of inter-subject similarity, our results support the use of functional alignment to improve inter-subject comparisons in the face of variable structure-function organization. We provide open implementations of all methods used.

Brain topography beyond parcellations: Local gradients of functional maps.

Functional neuroimaging provides the unique opportunity to characterize brain regions based on their response to tasks or ongoing activity. As such, it holds the premise to capture brain spatial organization. Yet, the conceptual framework to describe this organization has remained elusive: on the one hand, parcellations build implicitly on a piecewise constant organization, i.e. flat regions separated by sharp boundaries; on the other hand, the recently popularized concept of functional gradient hints instead at a smooth structure. Noting that both views converge to a topographic scheme that pieces together local variations of functional features, we perform a quantitative assessment of local gradient-based models. Using as a driving case the prediction of functional Magnetic Resonance Imaging (fMRI) data -concretely, the prediction of task-fMRI from rest-fMRI maps across subjects- we develop a parcel-wise linear regression model based on a dictionary of reference topographies. Our method uses multiple random parcellations -as opposed to a single fixed parcellation- and aggregates estimates across these parcellations to predict functional features in left-out subjects. Our experiments demonstrate the existence of an optimal cardinality of the parcellation to capture local gradients of functional maps.

Subject-specific segregation of functional territories based on deep phenotyping.

Functional magnetic resonance imaging (fMRI) has opened the possibility to investigate how brain activity is modulated by behavior. Most studies so far are bound to one single task, in which functional responses to a handful of contrasts are analyzed and reported as a group average brain map. Contrariwise, recent data-collection efforts have started to target a systematic spatial representation of multiple mental functions. In this paper, we leverage the Individual Brain Charting (IBC) dataset-a high-resolution task-fMRI dataset acquired in a fixed environment-in order to study the feasibility of individual mapping. First, we verify that the IBC brain maps reproduce those obtained from previous, large-scale datasets using the same tasks. Second, we confirm that the elementary spatial components, inferred across all tasks, are consistently mapped within and, to a lesser extent, across participants. Third, we demonstrate the relevance of the topographic information of the individual contrast maps, showing that contrasts from one task can be predicted by contrasts from other tasks. At last, we showcase the benefit of contrast accumulation for the fine functional characterization of brain regions within a prespecified network. To this end, we analyze the cognitive profile of functional territories pertaining to the language network and prove that these profiles generalize across participants.

Subspecialization within default mode nodes characterized in 10,000 UK Biobank participants.

The human default mode network (DMN) is implicated in several unique mental capacities. In this study, we tested whether brain-wide interregional communication in the DMN can be derived from population variability in intrinsic activity fluctuations, gray-matter morphology, and fiber tract anatomy. In a sample of 10,000 UK Biobank participants, pattern-learning algorithms revealed functional coupling states in the DMN that are linked to connectivity profiles between other macroscopical brain networks. In addition, DMN gray matter volume was covaried with white matter microstructure of the fornix. Collectively, functional and structural patterns unmasked a possible division of labor within major DMN nodes: Subregions most critical for cortical network interplay were adjacent to subregions most predictive of fornix fibers from the hippocampus that processes memories and places.

Multi-subject MEG/EEG source imaging with sparse multi-task regression.

Magnetoencephalography and electroencephalography (M/EEG) are non-invasive modalities that measure the weak electromagnetic fields generated by neural activity. Estimating the location and magnitude of the current sources that generated these electromagnetic fields is an inverse problem. Although it can be cast as a linear regression, this problem is severely ill-posed as the number of observations, which equals the number of sensors, is small. When considering a group study, a common approach consists in carrying out the regression tasks independently for each subject using techniques such as MNE or sLORETA. An alternative is to jointly localize sources for all subjects taken together, while enforcing some similarity between them. By pooling S subjects in a single joint regression, the number of observations is S times larger, potentially making the problem better posed and offering the ability to identify more sources with greater precision. Here we show how the coupling of the different regression problems can be done through a multi-task regularization that promotes focal source estimates. To take into account intersubject variabilities, we propose the Minimum Wasserstein Estimates (MWE). Thanks to a new joint regression method based on optimal transport (OT) metrics, MWE does not enforce perfect overlap of activation foci for all subjects but rather promotes spatial proximity on the cortical mantle. Besides, by estimating the noise level of each subject, MWE copes with the subject-specific signal-to-noise ratios with only one regularization parameter. On realistic simulations, MWE decreases the localization error by up to 4 â€‹mm per source compared to individual solutions. Experiments on the Cam-CAN dataset show improvements in spatial specificity in population imaging compared to individual models such as dSPM as well as a state-of-the-art Bayesian group level model. Our analysis of a multimodal dataset shows how multi-subject source localization reduces the gap between MEG and fMRI for brain mapping.

Fine-grain atlases of functional modes for fMRI analysis.

Population imaging markedly increased the size of functional-imaging datasets, shedding new light on the neural basis of inter-individual differences. Analyzing these large data entails new scalability challenges, computational and statistical. For this reason, brain images are typically summarized in a few signals, for instance reducing voxel-level measures with brain atlases or functional modes. A good choice of the corresponding brain networks is important, as most data analyses start from these reduced signals. We contribute finely-resolved atlases of functional modes, comprising from 64 to 1024 networks. These dictionaries of functional modes (DiFuMo) are trained on millions of fMRI functional brain volumes of total size 2.4 â€‹TB, spanned over 27 studies and many research groups. We demonstrate the benefits of extracting reduced signals on our fine-grain atlases for many classic functional data analysis pipelines: stimuli decoding from 12,334 brain responses, standard GLM analysis of fMRI across sessions and individuals, extraction of resting-state functional-connectomes biomarkers for 2500 individuals, data compression and meta-analysis over more than 15,000 statistical maps. In each of these analysis scenarii, we compare the performance of our functional atlases with that of other popular references, and to a simple voxel-level analysis. Results highlight the importance of using high-dimensional "soft" functional atlases, to represent and analyze brain activity while capturing its functional gradients. Analyses on high-dimensional modes achieve similar statistical performance as at the voxel level, but with much reduced computational cost and higher interpretability. In addition to making them available, we provide meaningful names for these modes, based on their anatomical location. It will facilitate reporting of results.

Multi-scale network regression for brain-phenotype associations.

Brain networks are increasingly characterized at different scales, including summary statistics, community connectivity, and individual edges. While research relating brain networks to behavioral measurements has yielded many insights into brain-phenotype relationships, common analytical approaches only consider network information at a single scale. Here, we designed, implemented, and deployed Multi-Scale Network Regression (MSNR), a penalized multivariate approach for modeling brain networks that explicitly respects both edge- and community-level information by assuming a low rank and sparse structure, both encouraging less complex and more interpretable modeling. Capitalizing on a large neuroimaging cohort (n = 1, 051), we demonstrate that MSNR recapitulates interpretable and statistically significant connectivity patterns associated with brain development, sex differences, and motion-related artifacts. Compared to single-scale methods, MSNR achieves a balance between prediction performance and model complexity, with improved interpretability. Together, by jointly exploiting both edge- and community-level information, MSNR has the potential to yield novel insights into brain-behavior relationships.

Benchmarking functional connectome-based predictive models for resting-state fMRI.

Functional connectomes reveal biomarkers of individual psychological or clinical traits. However, there is great variability in the analytic pipelines typically used to derive them from rest-fMRI cohorts. Here, we consider a specific type of studies, using predictive models on the edge weights of functional connectomes, for which we highlight the best modeling choices. We systematically study the prediction performances of models in 6 different cohorts and a total of 2000 individuals, encompassing neuro-degenerative (Alzheimer's, Post-traumatic stress disorder), neuro-psychiatric (Schizophrenia, Autism), drug impact (Cannabis use) clinical settings and psychological trait (fluid intelligence). The typical prediction procedure from rest-fMRI consists of three main steps: defining brain regions, representing the interactions, and supervised learning. For each step we benchmark typical choices: 8 different ways of defining regions -either pre-defined or generated from the rest-fMRI data- 3 measures to build functional connectomes from the extracted time-series, and 10 classification models to compare functional interactions across subjects. Our benchmarks summarize more than 240 different pipelines and outline modeling choices that show consistent prediction performances in spite of variations in the populations and sites. We find that regions defined from functional data work best; that it is beneficial to capture between-region interactions with tangent-based parametrization of covariances, a midway between correlations and partial correlation; and that simple linear predictors such as a logistic regression give the best predictions. Our work is a step forward to establishing reproducible imaging-based biomarkers for clinical settings.

The functional database of the ARCHI project: Potential and perspectives.

More than two decades of functional magnetic resonance imaging (fMRI) of the human brain have succeeded to identify, with a growing level of precision, the neural basis of multiple cognitive skills within various domains (perception, sensorimotor processes, language, emotion and social cognition …). Progress has been made in the comprehension of the functional organization of localized brain areas. However, the long time required for fMRI acquisition limits the number of experimental conditions performed in a single individual. As a consequence, distinct brain localizations have mostly been studied in separate groups of participants, and their functional relationships at the individual level remain poorly understood. To address this issue, we report here preliminary results on a database of fMRI data acquired on 78 individuals who each performed a total of 29 experimental conditions, grouped in 4 cross-domains functional localizers. This protocol has been designed to efficiently isolate, in a single session, the brain activity associated with language, numerical representation, social perception and reasoning, premotor and visuomotor representations. Analyses are reported at the group and at the individual level, to establish the ability of our protocol to selectively capture distinct regions of interest in a very short time. Test-retest reliability was assessed in a subset of participants. The activity evoked by the different contrasts of the protocol is located in distinct brain networks that, individually, largely replicate previous findings and, taken together, cover a large proportion of the cortical surface. We provide detailed analyses of a subset of regions of relevance: the left frontal, left temporal and middle frontal cortices. These preliminary analyses highlight how combining such a large set of functional contrasts may contribute to establish a finer-grained brain atlas of cognitive functions, especially in regions of high functional overlap. Detailed structural images (structural connectivity, micro-structures, axonal diameter) acquired in the same individuals in the context of the ARCHI database provide a promising situation to explore functional/structural interdependence. Additionally, this protocol might also be used as a way to establish individual neurofunctional signatures in large cohorts.

Brain-based ranking of cognitive domains to predict schizophrenia.

Schizophrenia is a devastating brain disorder that disturbs sensory perception, motor action, and abstract thought. Its clinical phenotype implies dysfunction of various mental domains, which has motivated a series of theories regarding the underlying pathophysiology. Aiming at a predictive benchmark of a catalog of cognitive functions, we developed a data-driven machine-learning strategy and provide a proof of principle in a multisite clinical dataset (n = 324). Existing neuroscientific knowledge on diverse cognitive domains was first condensed into neurotopographical maps. We then examined how the ensuing meta-analytic cognitive priors can distinguish patients and controls using brain morphology and intrinsic functional connectivity. Some affected cognitive domains supported well-studied directions of research on auditory evaluation and social cognition. However, rarely suspected cognitive domains also emerged as disease relevant, including self-oriented processing of bodily sensations in gustation and pain. Such algorithmic charting of the cognitive landscape can be used to make targeted recommendations for future mental health research.

SmartPulse, a machine learning approach for calibration-free dynamic RF shimming: Preliminary study in a clinical environment.

PURPOSE: A calibration-free pulse design method is introduced to alleviate B1+ artifacts in clinical routine with parallel transmission at high field, dealing with significant inter-subject variability, found for instance in the abdomen. THEORY AND METHODS: From a dual-transmit 3T scanner, a database of B1+ and off-resonance abdominal maps from 50 subjects was first divided into 3 clusters based on mutual affinity between their respective tailored kT -points pulses. For each cluster, a kT -points pulse was computed, minimizing normalized root-mean-square flip angle deviations simultaneously for all subjects comprised in it. Using 30 additional subjects' field distributions, a machine learning classifier was trained on this 80-labeled-subject database to recognize the best pulse from the 3 ones available, relying only on patient features accessible from the preliminary localizer sequence present in all protocols. This so-called SmartPulse process was experimentally tested on an additional 53-subject set and compared with other pulse types: vendor's hard calibration-free dual excitation, tailored static radiofrequency shimming, universal and tailored kT -points pulses. RESULTS: SmartPulse outperformed both calibration-free approaches. Tailored static radiofrequency shimming yielded similar flip angle homogeneity for most patients but broke down for some while SmartPulse remained robust. Although flip angle homogeneity was systematically better with tailored kT -points, the difference was barely noticeable on in vivo images. CONCLUSION: The proposed method paves the way toward an efficient trade-off between tailored and universal pulse design approaches for large inter-subject variability. With no need for on-line field mapping or pulse design, it can fit seamlessly into a clinical protocol.