A Stereotactic Probabilistic Atlas for the Major Cerebral Arteries

Segmentation of stroke lesions using transformers-augmented MRI analysis

Accurate segmentation of chronic stroke lesions from mono-spectral magnetic resonance imaging scans (e.g., T1-weighted images) is a difficult task due to the arbitrary shape, complex texture, variable size and intensities, and varied locations of the lesions. Due to this inherent spatial heterogeneity, existing machine learning methods have shown moderate performance for chronic lesion delineation. In this study, we introduced: (1) a method that integrates transformers' deformable feature attention mechanism with convolutional deep learning architecture to improve the accuracy and generalizability of stroke lesion segmentation, and (2) an ecological data augmentation technique based on inserting real lesions into intact brain regions. Our combination of these two approaches resulted in a significant increase in segmentation performance, with a Dice index of 0.82 (±0.39), outperforming the existing methods trained and tested on the same Anatomical Tracings of Lesions After Stroke (ATLAS) 2022 dataset. Our method performed relatively well even for cases with small stroke lesions. We validated the robustness of our method through an ablation study and by testing it on new unseen brain scans from the Ischemic Stroke Lesion Segmentation (ISLES) 2015 dataset. Overall, our proposed approach of transformers with ecological data augmentation offers a robust way to delineate chronic stroke lesions with clinically relevant accuracy. Our method can be extended to other challenging tasks that require automated detection and segmentation of diverse brain abnormalities from clinical scans.

Robust data-driven segmentation of pulsatile cerebral vessels using functional magnetic resonance imaging

Functional magnetic resonance imaging (fMRI) captures rich physiological and neuronal information that can offer insights into neurofluid dynamics, vascular health, and waste clearance function. The availability of cerebral vessel segmentation could facilitate fluid dynamics research in fMRI. However, without magnetic resonance angiography scans, cerebral vessel segmentation is challenging and time-consuming. This study leverages cardiac-induced pulsatile fMRI signal to develop a data-driven, automatic segmentation of large cerebral arteries and the superior sagittal sinus (SSS). The method was validated in a local dataset by comparing it to ground truth cerebral artery and SSS segmentations. Using the Human Connectome Project (HCP) aging dataset, the method's reproducibility was tested on 422 participants aged 36 to 100 years, each with four repeated fMRI scans. The method demonstrated high reproducibility, with an intraclass correlation coefficient > 0.7 in both cerebral artery and SSS segmentation volumes. This study demonstrates that the large cerebral arteries and SSS can be reproducibly and automatically segmented in fMRI datasets, facilitating the investigation of fluid dynamics in these regions.

Static and dynamic analysis of cerebral blood flow in fifty-six large arterial vessel networks

Objective.The cerebral vasculature is formed of an intricate network of blood vessels over many different length scales. Changes in their structure and connection are implicated in multiple cerebrovascular and neurological disorders. In this study, we present a novel approach to the quantitative analysis of the cerebral macrovasculature using computational and mathematical tools in a large dataset.Approach.We analysed a publicly available vessel dataset from a cohort of 56 (32/24F/M) healthy subjects. This dataset includes digital reconstructions of human brain macrovasculatures. We then propose a new mathematical model to compute blood flow dynamics and pressure distributions within these 56-representative cerebral macrovasculatures and quantify the results across this cohort.Main results.Statistical analysis showed that the steady state level of cerebrovascular resistance (CVR) gradually increases with age in both men and women. These age-related changes in CVR are in good agreement with previously reported values. All subjects were found to have only small phase angles (<6°) between blood pressure and blood flow at the cardiac frequency.Significance.These results showed that the dynamic component of blood flow adds very little phase shift at the cardiac frequency, which implies that the cerebral macrocirculation can be regarded as close to steady state in its behaviour, at least in healthy populations, irrespective of age or sex. This implies that the phase shift observed in measurements of blood flow in cerebral vessels is caused by behaviour further down the vascular bed. This behaviour is important for future statistical models of the dynamic maintenance of oxygen and nutrient supply to the brain.

Reliability of the TMS-evoked potential in dorsolateral prefrontal cortex

We currently lack a reliable method to probe cortical excitability noninvasively from the human dorsolateral prefrontal cortex (dlPFC). We recently found that the strength of early and local dlPFC transcranial magnetic stimulation (TMS)-evoked potentials (EL-TEPs) varied widely across dlPFC subregions. Despite these differences in response amplitude, reliability at each target is unknown. Here we quantified within-session reliability of dlPFC EL-TEPs after TMS to six left dlPFC subregions in 15 healthy subjects. We evaluated reliability (concordance correlation coefficient [CCC]) across targets, time windows, quantification methods, regions of interest, sensor- vs. source-space, and number of trials. On average, the medial target was most reliable (CCC = 0.78) and the most anterior target was least reliable (CCC = 0.24). However, all targets except the most anterior were reliable (CCC > 0.7) using at least one combination of the analytical parameters tested. Longer (20 to 60 ms) and later (30 to 60 ms) windows increased reliability compared to earlier and shorter windows. Reliable EL-TEPs (CCC up to 0.86) were observed using only 25 TMS trials at a medial dlPFC target. Overall, medial dlPFC targeting, wider windows, and peak-to-peak quantification improved reliability. With careful selection of target and analytic parameters, highly reliable EL-TEPs can be extracted from the dlPFC after only a small number of trials.

Magnetoencephalography Atlas Viewer for Dipole Localization and Viewing

Magnetoencephalography (MEG) is a noninvasive neuroimaging technique widely recognized for epilepsy and tumor mapping. MEG clinical reporting requires a multidisciplinary team, including expert input regarding each dipole's anatomic localization. Here, we introduce a novel tool, the "Magnetoencephalography Atlas Viewer" (MAV), which streamlines this anatomical analysis. The MAV normalizes the patient's Magnetic Resonance Imaging (MRI) to the Montreal Neurological Institute (MNI) space, reverse-normalizes MNI atlases to the native MRI, identifies MEG dipole files, and matches dipoles' coordinates to their spatial location in atlas files. It offers a user-friendly and interactive graphical user interface (GUI) for displaying individual dipoles, groups, coordinates, anatomical labels, and a tri-planar MRI view of the patient with dipole overlays. It evaluated over 273 dipoles obtained in clinical epilepsy subjects. Consensus-based ground truth was established by three neuroradiologists, with a minimum agreement threshold of two. The concordance between the ground truth and MAV labeling ranged from 79% to 84%, depending on the normalization method. Higher concordance rates were observed in subjects with minimal or no structural abnormalities on the MRI, ranging from 80% to 90%. The MAV provides a straightforward MEG dipole anatomic localization method, allowing a nonspecialist to prepopulate a report, thereby facilitating and reducing the time of clinical reporting.

Automated anatomical labeling of the intracranial arteries via deep learning in computed tomography angiography

Background and purpose: Anatomical labeling of the cerebral vasculature is a crucial topic in determining the morphological nature and characterizing the vital variations of vessels, yet precise labeling of the intracranial arteries is time-consuming and challenging, given anatomical structural variability and surging imaging data. We present a U-Net-based deep learning (DL) model to automatically label detailed anatomical segments in computed tomography angiography (CTA) for the first time. The trained DL algorithm was further tested on a clinically relevant set for the localization of intracranial aneurysms (IAs). Methods: 457 examinations with varying degrees of arterial stenosis were used to train, validate, and test the model, aiming to automatically label 42 segments of the intracranial arteries [e.g., 7 segments of the internal carotid artery (ICA)]. Evaluation metrics included Dice similarity coefficient (DSC), mean surface distance (MSD), and Hausdorff distance (HD). Additionally, 96 examinations containing at least one IA were enrolled to assess the model's potential in enhancing clinicians' precision in IA localization. A total of 5 clinicians with different experience levels participated as readers in the clinical experiment and identified the precise location of IA without and with algorithm assistance, where there was a washout period of 14 days between two interpretations. The diagnostic accuracy, time, and mean interrater agreement (Fleiss' Kappa) were calculated to assess the differences in clinical performance of clinicians. Results: The proposed model exhibited notable labeling performance on 42 segments that included 7 anatomical segments of ICA, with the mean DSC of 0.88, MSD of 0.82 mm and HD of 6.59 mm. Furthermore, the model demonstrated superior labeling performance in healthy subjects compared to patients with stenosis (DSC: 0.91 vs. 0.89, p < 0.05; HD: 4.75 vs. 6.19, p < 0.05). Concurrently, clinicians with model predictions achieved significant improvements when interpreting the precise location of IA. The clinicians' mean accuracy increased by 0.04 (p = 0.003), mean time to diagnosis reduced by 9.76 s (p < 0.001), and mean interrater agreement (Fleiss' Kappa) increased by 0.07 (p = 0.029). Conclusion: Our model stands proficient for labeling intracranial arteries using the largest CTA dataset. Crucially, it demonstrates clinical utility, helping prioritize the patients with high risks and ease clinical workload.

Using deep learning for an automatic detection and classification of the vascular bifurcations along the Circle of Willis

Most of the intracranial aneurysms (ICA) occur on a specific portion of the cerebral vascular tree named the Circle of Willis (CoW). More particularly, they mainly arise onto fifteen of the major arterial bifurcations constituting this circular structure. Hence, for an efficient and timely diagnosis it is critical to develop some methods being able to accurately recognize each Bifurcation of Interest (BoI). Indeed, an automatic extraction of the bifurcations presenting the higher risk of developing an ICA would offer the neuroradiologists a quick glance at the most alarming areas. Due to the recent efforts on Artificial Intelligence, Deep Learning turned out to be the best performing technology for many pattern recognition tasks. Moreover, various methods have been particularly designed for medical image analysis purposes. This study intends to assist the neuroradiologists to promptly locate any bifurcation presenting a high risk of ICA occurrence. It can be seen as a Computer Aided Diagnosis scheme, where the Artificial Intelligence facilitates the access to the regions of interest within the MRI. In this work, we propose a method for a fully automatic detection and recognition of the bifurcations of interest forming the Circle of Willis. Several neural networks architectures have been tested, and we thoroughly evaluate the bifurcation recognition rate.

Reliability of the TMS-evoked potential in dorsolateral prefrontal cortex

Background

We currently lack a robust and reliable method to probe cortical excitability noninvasively from the human dorsolateral prefrontal cortex (dlPFC), a region heavily implicated in psychiatric disorders. We recently found that the strength of early and local dlPFC single pulse transcranial magnetic stimulation (TMS)-evoked potentials (EL-TEPs) varied widely depending on the anatomical subregion probed, with more medial regions eliciting stronger responses than anterolateral sites. Despite these differences in amplitude of response, the reliability at each target is not known.

Objective

To evaluate the reliability of EL-TEPs across the dlPFC.

Methods

In 15 healthy subjects, we quantified within-session reliability of dlPFC EL-TEPs after single pulse TMS to six dlPFC subregions. We evaluated the concordance correlation coefficient (CCC) across targets and analytical parameters including time window, quantification method, region of interest, sensor-vs. source-space, and number of trials.

Results

At least one target in the anterior and posterior dlPFC produced reliable EL-TEPs (CCC>0.7). The medial target was most reliable (CCC = 0.78) and the most anterior target was least reliable (CCC = 0.24). ROI size and type (sensor vs. source space) did not affect reliability. Longer (20-60 ms, CCC = 0.62) and later (30-60 ms, CCC = 0.61) time windows resulted in higher reliability compared to earlier and shorter (20-40 ms, CCC 0.43; 20-50 ms, CCC = 0.55) time windows. Peak-to-peak quantification resulted in higher reliability than the mean of the absolute amplitude. Reliable EL-TEPs (CCC up to 0.86) were observed using only 25 TMS trials for a medial dlPFC target.

Conclusions

Medial TMS location, wider time window (20-60ms), and peak-to-peak quantification improved reliability. Highly reliable EL-TEPs can be extracted from dlPFC after only a small number of trials.

Highlights

Medial dlPFC target improved EL-TEP reliability compared to anterior targets.After optimizing analytical parameters, at least one anterior and one posterior target was reliable (CCC>0.7).Longer (20-60 ms) and later (30-60 ms) time windows were more reliable than earlier and shorter (20-40 ms or 20-50 ms) latencies.Peak-to-peak quantification resulted in higher reliability compared to the mean of the absolute amplitude.As low as 25 trials can yield reliable EL-TEPs from the dlPFC.

Automated in-depth cerebral arterial labelling using cerebrovascular vasculature reframing and deep neural networks

Identifying the cerebral arterial branches is essential for undertaking a computational approach to cerebrovascular imaging. However, the complexity and inter-individual differences involved in this process have not been thoroughly studied. We used machine learning to examine the anatomical profile of the cerebral arterial tree. The method is less sensitive to inter-subject and cohort-wise anatomical variations and exhibits robust performance with an unprecedented in-depth vessel range. We applied machine learning algorithms to disease-free healthy control subjects (n = 42), patients with stroke with intracranial atherosclerosis (ICAS) (n = 46), and patients with stroke mixed with the existing controls (n = 69). We trained and tested 70% and 30% of each study cohort, respectively, incorporating spatial coordinates and geometric vessel feature vectors. Cerebral arterial images were analyzed based on the 'segmentation-stacking' method using magnetic resonance angiography. We precisely classified the cerebral arteries across the exhaustive scope of vessel components using advanced geometric characterization, redefinition of vessel unit conception, and post-processing algorithms. We verified that the neural network ensemble, with multiple joint models as the combined predictor, classified all vessel component types independent of inter-subject variations in cerebral arterial anatomy. The validity of the categorization performance of the model was tested, considering the control, ICAS, and control-blended stroke cohorts, using the area under the receiver operating characteristic (ROC) curve and precision-recall curve. The classification accuracy rarely fell outside each image's 90-99% scope, independent of cohort-dependent cerebrovascular structural variations. The classification ensemble was calibrated with high overall area rates under the ROC curve of 0.99-1.00 [0.97-1.00] in the test set across various study cohorts. Identifying an all-inclusive range of vessel components across controls, ICAS, and stroke patients, the accuracy rates of the prediction were: internal carotid arteries, 91-100%; middle cerebral arteries, 82-98%; anterior cerebral arteries, 88-100%; posterior cerebral arteries, 87-100%; and collections of superior, anterior inferior, and posterior inferior cerebellar arteries, 90-99% in the chunk-level classification. Using a voting algorithm on the queued classified vessel factors and anatomically post-processing the automatically classified results intensified quantitative prediction performance. We employed stochastic clustering and deep neural network ensembles. Ma-chine intelligence-assisted prediction of vessel structure allowed us to personalize quantitative predictions of various types of cerebral arterial structures, contributing to precise and efficient decisions regarding the cerebrovascular disease.

4D flow MRI hemodynamic biomarkers for cerebrovascular diseases

Alterations in cerebral blood flow are common in several neurological diseases among the elderly including stroke, cerebral small vessel disease, vascular dementia, and Alzheimer's disease. 4D flow magnetic resonance imaging (MRI) is a relatively new technique to investigate cerebrovascular disease, and makes it possible to obtain time-resolved blood flow measurements of the entire cerebral arterial venous vasculature and can be used to derive a repertoire of hemodynamic biomarkers indicative of cerebrovascular health. The information that can be obtained from one single 4D flow MRI scan allows both the investigation of aberrant flow patterns at a focal location in the vasculature as well as estimations of brain-wide disturbances in blood flow. Such focal and global hemodynamic biomarkers show the potential of being sensitive to impending cerebrovascular disease and disease progression and can also become useful during planning and follow-up of interventions aiming to restore a normal cerebral circulation. Here, we describe 4D flow MRI approaches for analyzing the cerebral vasculature. We then survey key hemodynamic biomarkers that can be reliably assessed using the technique. Finally, we highlight cerebrovascular diseases where one or multiple hemodynamic biomarkers are of central interest.

4D flow MRI for non-invasive measurement of blood flow in the brain: A systematic review

The brain's vasculature is essential for brain health and its dysfunction contributes to the onset and development of many dementias and neurological disorders. While numerous in vivo imaging techniques exist to investigate cerebral haemodynamics in humans, phase-contrast magnetic resonance imaging (MRI) has emerged as a reliable, non-invasive method of quantifying blood flow within intracranial vessels. In recent years, an advanced form of this method, known as 4D flow, has been developed and utilised in patient studies, where its ability to capture complex blood flow dynamics within any major vessel across the acquired volume has proved effective in collecting large amounts of information in a single scan. While extremely promising as a method of examining the vascular system's role in brain-related diseases, the collection of 4D data can be time-consuming, meaning data quality has to be traded off against the acquisition time. Here, we review the available literature to examine 4D flow's capabilities in assessing physiological and pathological features of the cerebrovascular system. Emerging techniques such as dynamic velocity-encoding and advanced undersampling methods, combined with increasingly high-field MRI scanners, are likely to bring 4D flow to the forefront of cerebrovascular imaging studies in the years to come.

Evolution of Human Brain Atlases in Terms of Content, Applications, Functionality, and Availability

Human brain atlases have been evolving tremendously, propelled recently by brain big projects, and driven by sophisticated imaging techniques, advanced brain mapping methods, vast data, analytical strategies, and powerful computing. We overview here this evolution in four categories: content, applications, functionality, and availability, in contrast to other works limited mostly to content. Four atlas generations are distinguished: early cortical maps, print stereotactic atlases, early digital atlases, and advanced brain atlas platforms, and 5 avenues in electronic atlases spanning the last two generations. Content-wise, new electronic atlases are categorized into eight groups considering their scope, parcellation, modality, plurality, scale, ethnicity, abnormality, and a mixture of them. Atlas content developments in these groups are heading in 23 various directions. Application-wise, we overview atlases in neuroeducation, research, and clinics, including stereotactic and functional neurosurgery, neuroradiology, neurology, and stroke. Functionality-wise, tools and functionalities are addressed for atlas creation, navigation, individualization, enabling operations, and application-specific. Availability is discussed in media and platforms, ranging from mobile solutions to leading-edge supercomputers, with three accessibility levels. The major application-wise shift has been from research to clinical practice, particularly in stereotactic and functional neurosurgery, although clinical applications are still lagging behind the atlas content progress. Atlas functionality also has been relatively neglected until recently, as the management of brain data explosion requires powerful tools. We suggest that the future human brain atlas-related research and development activities shall be founded on and benefit from a standard framework containing the core virtual brain model cum the brain atlas platform general architecture.

Computer Modeling of Clot Retrieval-Circle of Willis

Endovascular clot retrieval, often referred to as mechanical thrombectomy, has transformed the treatment of patients with ischemic stroke based on an underlying large cerebral vessel occlusion, ranging from the extracranial internal carotid artery (ICA) to the M1 (proximal) segment of the middle cerebral artery (MCA). The aim of this study was to evaluate the effect of a progressive occlusion of the extracranial portion of the ICA on the cerebral blood flow either with a conventional guiding catheter or a balloon-guiding catheter, which enables the operator to completely occlude the parent artery by inflating the balloon around the tip of this type of guiding catheter. We evaluated the impact of flow reduction in the ICA in the setting of ipsilateral MCA occlusion given the different configurations of the circle of Willis (CoW). The computer model of cerebral arteries was based on anatomical works by Rhoton (1) and van der Eecken (2). The interactive experimental results are available on the web at https://gntem3.shinyapps.io/ecrsim. In the setting of left MCA occlusion, compensation from the anterior and posterior communicating artery preserved the flow in the left anterior cerebral artery (ACA) but not the left MCA branches. Under selected CoW configurations, such as classic, missing Acom, or missing A1 segment of the ACA and concurrent right ICA occlusion, there was a progressive decrease of flow in the left ACA to a minimum of 78% when the simulated catheter fully occluded the left ICA. Flow collapsed (<10%) in the left ACA and MCA branches under CoW configurations, such as bilateral fetal PCA. In summary, compensatory flow collapsed under certain clot retrieval scenarios and unusual configurations of CoW.

Direct comparison of activation maps during galvanic vestibular stimulation: A hybrid H2[15 O] PET-BOLD MRI activation study

Previous unimodal PET and fMRI studies in humans revealed a reproducible vestibular brain activation pattern, but with variations in its weighting and expansiveness. Hybrid studies minimizing methodological variations at baseline conditions are rare and still lacking for task-based designs. Thus, we applied for the first time hybrid 3T PET-MRI scanning (Siemens mMR) in healthy volunteers using galvanic vestibular stimulation (GVS) in healthy volunteers in order to directly compare H₂¹⁵O-PET and BOLD MRI responses. List mode PET acquisition started with the injection of 750 MBq H₂¹⁵O simultaneously to MRI EPI sequences. Group-level statistical parametric maps were generated for GVS vs. rest contrasts of PET, MR-onset (event-related), and MR-block. All contrasts showed a similar bilateral vestibular activation pattern with remarkable proximity of activation foci. Both BOLD contrasts gave more bilateral wide-spread activation clusters than PET; no area showed contradictory signal responses. PET still confirmed the right-hemispheric lateralization of the vestibular system, whereas BOLD-onset revealed only a tendency. The reciprocal inhibitory visual-vestibular interaction concept was confirmed by PET signal decreases in primary and secondary visual cortices, and BOLD-block decreases in secondary visual areas. In conclusion, MRI activation maps contained a mixture of CBF measured using H₂¹⁵O-PET and additional non-CBF effects, and the activation-deactivation pattern of the BOLD-block appears to be more similar to the H₂¹⁵O-PET than the BOLD-onset.

Boundary and vulnerability estimation of the internal borderzone using ischemic stroke lesion mapping

Distinction between deep and superficial middle cerebral artery (MCA) territories and their junctional vascular area (the internal borderzone or IBZ) constitutes a predictor of stroke patient outcome. However, the IBZ boundaries are not well-defined because of substantial anatomical variance. Here, we built a statistical estimate of the IBZ and tested its vulnerability to ischemia using an independent sample. First, we used delineated lesions of 122 patients suffering of chronic ischemic stroke grouped in deep, superficial and territorial topographies and statistical comparisons to generate a probabilistic estimate of the IBZ. The IBZ extended from the insular cortex to the internal capsule and the anterior part of the caudate nucleus head. The IBZ showed the highest lesion frequencies (~30% on average across IBZ voxels) in our chronic stroke patients but also in an independent sample of 87 acute patients. Additionally, the most important apparent diffusion coefficient reductions (−6%), which reflect stroke severity, were situated within our IBZ estimate. The IBZ was most severely injured in case of a territorial infarction. Then, our results are in favour of an increased IBZ vulnerability to ischemia. Moreover, our probabilistic estimates of deep, superficial and IBZ regions can help the everyday spatial classification of lesions.

Cerebral arterial pulsatility is associated with features of small vessel disease in patients with acute stroke and TIA: a 4D flow MRI study

Cerebral small vessel disease (SVD) is a major cause of stroke and cognitive impairment. However, the underlying mechanisms behind SVD are still poorly understood. High cerebral arterial pulsatility has been suggested as a possible cause of SVD. In population studies, arterial pulsatility has been linked to white matter hyperintensities (WMH), cerebral atrophy, and cognitive impairment, all features of SVD. In stroke, pulsatility data are scarce and contradictory. The aim of this study was to investigate the relationship between arterial pulsatility and SVD in stroke patients. With a cross-sectional design, 89 patients with acute ischemic stroke or TIA were examined with MRI. A neuropsychological assessment was performed 1 year later. Using 4D flow MRI, pulsatile indices (PI) were calculated for the internal carotid artery (ICA) and middle cerebral artery (M1, M3). Flow volume pulsatility (FVP), a measure corresponding to the cyclic expansion of the arterial tree, was calculated for the same locations. These parameters were assessed for associations with WMH volume, brain volume and cognitive function. ICA-FVP was associated with WMH volume (β = 1.67, 95% CI: [0.1, 3.24], p = 0.037). M1-PI and M1-FVP were associated with decreasing cognitive function (β = − 4.4, 95% CI: [− 7.7, − 1.1], p = 0.009 and β = − 13.15, 95% CI: [− 24.26, − 2.04], p = 0.02 respectively). In summary, this supports an association between arterial pulsatility and SVD in stroke patients, and provides a potential target for further research and preventative treatment. FVP may become a useful biomarker for assessing pulsatile stress with PCMRI and 4D flow MRI.

Assessment of Cerebral Blood Flow Pulsatility and Cerebral Arterial Compliance With 4D Flow MRI

BACKGROUND

Four-dimensional flow magnetic resonance imaging (4D flow MRI) enables efficient investigation of cerebral blood flow pulsatility in the cerebral arteries. This is important for exploring hemodynamic mechanisms behind vascular diseases associated with arterial pulsations.

PURPOSE

To investigate the feasibility of pulsatility assessments with 4D flow MRI, its agreement with reference two-dimensional phase-contrast MRI (2D PC-MRI) measurements, and to demonstrate how 4D flow MRI can be used to assess cerebral arterial compliance and cerebrovascular resistance in major cerebral arteries.

STUDY TYPE

Prospective.

SUBJECTS

Thirty-five subjects (20 women, 79 ± 5 years, range 70-91 years).

FIELD STRENGTH/SEQUENCE

4D flow MRI (PC-VIPR) and 2D PC-MRI acquired with a 3T scanner.

ASSESSMENT

Time-resolved flow was assessed in nine cerebral arteries. From the pulsatile flow waveform in each artery, amplitude (ΔQ), volume load (ΔV), and pulsatility index (PI) were calculated. To reduce high-frequency noise in the 4D flow MRI data, the flow waveforms were low-pass filtered. From the total cerebral blood flow, total PI (PI_tot ), total volume load (ΔV_tot ), cerebral arterial compliance (C), and cerebrovascular resistance (R) were calculated.

STATISTICAL TESTS

Two-tailed paired t-test, intraclass correlation (ICC).

RESULTS

There was no difference in ΔQ between 4D flow MRI and the reference (0.00 ± 0.022 ml/s, mean ± SEM, P = 0.97, ICC = 0.95, n = 310) with a cutoff frequency of 1.9 Hz and 15 cut plane long arterial segments. For ΔV, the difference was -0.006 ± 0.003 ml (mean ± SEM, P = 0.07, ICC = 0.93, n = 310) without filtering. Total R was 11.4 ± 2.41 mmHg/(ml/s) (mean ± SD) and C was 0.021 ± 0.009 ml/mmHg (mean ± SD). ΔV_tot was 1.21 ± 0.29 ml (mean ± SD) with an ICC of 0.82 compared with the reference. PI_tot was 1.08 ± 0.21 (mean ± SD).

DATA CONCLUSION

We successfully assessed 4D flow MRI cerebral arterial pulsatility, cerebral arterial compliance, and cerebrovascular resistance. Averaging of multiple cut planes and low-pass filtering was necessary to assess accurate peak-to-peak features in the flow rate waveforms.

LEVEL OF EVIDENCE

2 Technical Efficacy Stage: 2 J. Magn. Reson. Imaging 2020;51:1516-1525.

A statistical atlas of cerebral arteries generated using multi-center MRA datasets from healthy subjects

Magnetic resonance angiography (MRA) can capture the variation of cerebral arteries with high spatial resolution. These measurements include valuable information about the morphology, geometry, and density of brain arteries, which may be useful to identify risk factors for cerebrovascular and neurological diseases at an early time point. However, this requires knowledge about the distribution and morphology of vessels in healthy subjects. The statistical arterial brain atlas described in this work is a free and public neuroimaging resource that can be used to identify vascular morphological changes. The atlas was generated based on 544 freely available multi-center MRA and T1-weighted MRI datasets. The arteries were automatically segmented in each MRA dataset and used for vessel radius quantification. The binary segmentation and vessel size information were non-linearly registered to the MNI brain atlas using the T1-weighted MRI datasets to construct atlases of artery occurrence probability, mean artery radius, and artery radius standard deviation. This public neuroimaging resource improves the understanding of the distribution and size of arteries in the healthy human brain.

Accuracy of blood flow assessment in cerebral arteries with 4D flow MRI: Evaluation with three segmentation methods

Background

Accelerated 4D flow MRI allows for high‐resolution velocity measurements with whole‐brain coverage. Such scans are increasingly used to calculate flow rates of individual arteries in the vascular tree, but detailed information about the accuracy and precision in relation to different postprocessing options is lacking.

Purpose

To evaluate and optimize three proposed segmentation methods and determine the accuracy of in vivo 4D flow MRI blood flow rate assessments in major cerebral arteries, with high‐resolution 2D PCMRI as a reference.

Study Type

Prospective.

Subjects

Thirty‐five subjects (20 women, 79 ± 5 years, range 70–91 years).

Field Strength/Sequence

4D flow MRI with PC‐VIPR and 2D PCMRI acquired with a 3 T scanner.

Assessment

We compared blood flow rates measured with 4D flow MRI, to the reference, in nine main cerebral arteries. Lumen segmentation in the 4D flow MRI was performed with k‐means clustering using four different input datasets, and with two types of thresholding methods. The threshold was defined as a percentage of the maximum intensity value in the complex difference image. Local and global thresholding approaches were used, with evaluated thresholds from 6–26%.

Statistical Tests

Paired t‐test, F‐test, linear correlation (P < 0.05 was considered significant) along with intraclass correlation (ICC).

Results

With the thresholding methods, the lowest average flow difference was obtained for 20% local (0.02 ± 15.0 ml/min, ICC = 0.97, n = 310) or 10% global (0.08 ± 17.3 ml/min, ICC = 0.97, n = 310) thresholding with a significant lower standard deviation for local (F‐test, P = 0.01). For all clustering methods, we found a large systematic underestimation of flow compared with 2D PCMRI (16.1–22.3 ml/min).

Data Conclusion

A locally adapted threshold value gives a more stable result compared with a globally fixed threshold. 4D flow with the proposed segmentation method has the potential to become a useful reliable clinical tool for assessment of blood flow in the major cerebral arteries.

Level of Evidence: 2

Technical Efficacy: Stage 2

J. Magn. Reson. Imaging 2019;50:511–518.

The morphology of the human cerebrovascular system

While several methodologies exist for quantifying gray and white matter properties in humans, relatively little is known regarding the spatial organization and the intersubject variability of cerebral vessels. To resolve this, we developed a fast, open-source processing algorithm using advanced vessel segmentation schemes and iterative nonlinear registration to isolate, extract, and quantify cerebral vessels in susceptibility weighting imaging (SWI) and time-of-flight angiography (TOF-MRA) datasets acquired in a large cohort (n = 42) of healthy individuals. From this, whole-brain venous and arterial probabilistic maps were generated along with the computation of regional densities and diameters within regions based on popular anatomical and functional atlases. The results show that cerebral vasculature is highly heterogeneous, displaying disproportionally large vessel densities in brain areas such as the anterior and posterior cingulate, cuneus, precuneus, parahippocampus, insula, and temporal gyri. On average, venous densities were slightly higher and less variable across subjects than arterial. Moreover, regional variations in both venous and arterial density were significantly correlated to cortical thickness (R = 0.42). This publicly available new atlas of the human cerebrovascular system provides a first step toward quantifying morphological changes in the diseased brain and serving as a potential regression tool in fMRI analysis.

Human brain atlasing: past, present and future

We have recently witnessed an explosion of large-scale initiatives and projects addressing mapping, modeling, simulation and atlasing of the human brain, including the BRAIN Initiative, the Human Brain Project, the Human Connectome Project (HCP), the Big Brain, the Blue Brain Project, the Allen Brain Atlas, the Brainnetome, among others. Besides these large and international initiatives, there are numerous mid-size and small brain atlas-related projects. My contribution to these global efforts has been to create adult human brain atlases in health and disease, and to develop atlas-based applications. For over two decades with my R&D lab I developed 35 brain atlases, licensed to 67 companies and made available in about 100 countries. This paper has two objectives. First, it provides an overview of the state of the art in brain atlasing. Second, as it is already 20 years from the release of our first brain atlas, I summarise my past and present efforts, share my experience in atlas creation, validation and commercialisation, compare with the state of the art, and propose future directions.