Improving 3D edge detection for visual inspection of MRI coregistration and alignment

BACKGROUND

Visualizing edges is critical for neuroimaging. For example, edge maps enable quality assurance for the automatic alignment of an image from one modality (or individual) to another.

NEW METHOD

We suggest that using the second derivative (difference of Gaussian, or DoG) provides robust edge detection. This method is tuned by size (which is typically known in neuroimaging) rather than intensity (which is relative).

RESULTS

We demonstrate that this method performs well across a broad range of imaging modalities. The edge contours produced consistently form closed surfaces, whereas alternative methods may generate disconnected lines, introducing potential ambiguity in contiguity.

COMPARISON WITH EXISTING METHODS

Current methods for computing edges are based on either the first derivative of the image (FSL), or a variation of the Canny Edge detection method (AFNI). These methods suffer from two primary limitations. First, the crucial tuning parameter for each of these methods relates to the image intensity. Unfortunately, image intensity is relative for most neuroimaging modalities making the performance of these methods unreliable. Second, these existing approaches do not necessarily generate a closed edge/surface, which can reduce the ability to determine the correspondence between a represented edge and another image.

CONCLUSION

The second derivative is well suited for neuroimaging edge detection. We include this method as part of both the AFNI and FSL software packages, standalone code and online.

Multimodal active subspace analysis for computing assessment oriented subspaces from neuroimaging data

BACKGROUND

For successful biomarker discovery, it is essential to develop computational frameworks that summarize high-dimensional neuroimaging data in terms of involved sub-systems of the brain, while also revealing underlying heterogeneous functional and structural changes covarying with specific cognitive and biological traits. However, unsupervised decompositions do not inculcate clinical assessment information, while supervised approaches extract only individual feature importance, thereby impeding qualitative interpretation at the level of subspaces.

NEW METHOD

We present a novel framework to extract robust multimodal brain subspaces associated with changes in a given cognitive or biological trait. Our approach involves active subspace learning on the gradients of a trained machine learning model followed by clustering to extract and summarize the most salient and consistent subspaces associated with the target variable.

RESULTS

Through a rigorous cross-validation procedure on an Alzheimer's disease (AD) dataset, our framework successfully extracts multimodal subspaces specific to a given clinical assessment (e.g., memory and other cognitive skills), and also retains predictive performance in standard machine learning algorithms. We also show that the salient active subspace directions occur consistently across randomly sub-sampled repetitions of the analysis.

COMPARISON WITH EXISTING METHOD(S)

Compared to existing unsupervised decompositions based on principle component analysis, the subspace components in our framework retain higher predictive information.

CONCLUSIONS

As an important step towards biomarker discovery, our framework not only uncovers AD-related brain regions in the associated brain subspaces, but also enables automated identification of multiple underlying structural and functional sub-systems of the brain that collectively characterize changes in memory and proficiency in cognitive skills related to brain disorders like AD.

Potential of neuroimaging as a biomarker in amyotrophic lateral sclerosis: from structure to metabolism

Amyotrophic lateral sclerosis (ALS) is a progressive neurodegenerative disease characterized by motor neuron degeneration. The development of ALS involves metabolite alterations leading to tissue lesions in the nervous system. Recent advances in neuroimaging have significantly improved our understanding of the underlying pathophysiology of ALS, with findings supporting the corticoefferent axonal disease progression theory. Current studies on neuroimaging in ALS have demonstrated inconsistencies, which may be due to small sample sizes, insufficient statistical power, overinterpretation of findings, and the inherent heterogeneity of ALS. Deriving meaningful conclusions solely from individual imaging metrics in ALS studies remains challenging, and integrating multimodal imaging techniques shows promise for detecting valuable ALS biomarkers. In addition to giving an overview of the principles and techniques of different neuroimaging modalities, this review describes the potential of neuroimaging biomarkers in the diagnosis and prognostication of ALS. We provide an insight into the underlying pathology, highlighting the need for standardized protocols and multicenter collaborations to advance ALS research.

Combining neuropsychological assessment and structural neuroimaging to identify early Alzheimer's disease in a memory clinic cohort

INTRODUCTION

The current study examined the contributions of comprehensive neuropsychological assessment and volumetric assessment of selected mesial temporal subregions on structural magnetic resonance imaging (MRI) to identify patients with amnestic mild cognitive impairment (aMCI) and mild probable Alzheimer's disease (AD) dementia in a memory clinic cohort.

METHODS

Comprehensive neuropsychological assessment and automated entorhinal, transentorhinal, and hippocampal volume measurements were conducted in 40 healthy controls, 38 patients with subjective memory symptoms, 16 patients with aMCI, 16 patients with mild probable AD dementia. Multinomial logistic regression was used to compare the neuropsychological and MRI measures.

RESULTS

Combining the neuropsychological and MRI measures improved group membership prediction over the MRI measures alone but did not improve group membership prediction over the neuropsychological measures alone.

CONCLUSION

Comprehensive neuropsychological assessment was an important tool to evaluate cognitive impairment. The mesial temporal volumetric MRI measures contributed no diagnostic value over and above the determinations made through neuropsychological assessment.

SAN: Mitigating spatial covariance heterogeneity in cortical thickness data collected from multiple scanners or sites

In neuroimaging studies, combining data collected from multiple study sites or scanners is becoming common to increase the reproducibility of scientific discoveries. At the same time, unwanted variations arise by using different scanners (inter-scanner biases), which need to be corrected before downstream analyses to facilitate replicable research and prevent spurious findings. While statistical harmonization methods such as ComBat have become popular in mitigating inter-scanner biases in neuroimaging, recent methodological advances have shown that harmonizing heterogeneous covariances results in higher data quality. In vertex-level cortical thickness data, heterogeneity in spatial autocorrelation is a critical factor that affects covariance heterogeneity. Our work proposes a new statistical harmonization method called spatial autocorrelation normalization (SAN) that preserves homogeneous covariance vertex-level cortical thickness data across different scanners. We use an explicit Gaussian process to characterize scanner-invariant and scanner-specific variations to reconstruct spatially homogeneous data across scanners. SAN is computationally feasible, and it easily allows the integration of existing harmonization methods. We demonstrate the utility of the proposed method using cortical thickness data from the Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) study. SAN is publicly available as an R package.

Searching Reproducible Brain Features using NeuroMark: Templates for Different Age Populations and Imaging Modalities

A primary challenge to the data-driven analysis is the balance between poor generalizability of population-based research and characterizing more subject-, study- and population-specific variability. We previously introduced a fully automated spatially constrained independent component analysis (ICA) framework called NeuroMark and its functional MRI (fMRI) template. NeuroMark has been successfully applied in numerous studies, identifying brain markers reproducible across datasets and disorders. The first NeuroMark template was constructed based on young adult cohorts. We recently expanded on this initiative by creating a standardized normative multi-spatial-scale functional template using over 100,000 subjects, aiming to improve generalizability and comparability across studies involving diverse cohorts. While a unified template across the lifespan is desirable, a comprehensive investigation of the similarities and differences between components from different age populations might help systematically transform our understanding of the human brain by revealing the most well-replicated and variable network features throughout the lifespan. In this work, we introduced two significant expansions of NeuroMark templates first by generating replicable fMRI templates for infants, adolescents, and aging cohorts, and second by incorporating structural MRI (sMRI) and diffusion MRI (dMRI) modalities. Specifically, we built spatiotemporal fMRI templates based on 6,000 resting-state scans from four datasets. This is the first attempt to create robust ICA templates covering dynamic brain development across the lifespan. For the sMRI and dMRI data, we used two large publicly available datasets including more than 30,000 scans to build reliable templates. We employed a spatial similarity analysis to identify replicable templates and investigate the degree to which unique and similar patterns are reflective in different age populations. Our results suggest remarkably high similarity of the resulting adapted components, even across extreme age differences. With the new templates, the NeuroMark framework allows us to perform age-specific adaptations and to capture features adaptable to each modality, therefore facilitating biomarker identification across brain disorders. In sum, the present work demonstrates the generalizability of NeuroMark templates and suggests the potential of new templates to boost accuracy in mental health research and advance our understanding of lifespan and cross-modal alterations.

Accurate brain-age models for routine clinical MRI examinations

Convolutional neural networks (CNN) can accurately predict chronological age in healthy individuals from structural MRI brain scans. Potentially, these models could be applied during routine clinical examinations to detect deviations from healthy ageing, including early-stage neurodegeneration. This could have important implications for patient care, drug development, and optimising MRI data collection. However, existing brain-age models are typically optimised for scans which are not part of routine examinations (e.g., volumetric T1-weighted scans), generalise poorly (e.g., to data from different scanner vendors and hospitals etc.), or rely on computationally expensive pre-processing steps which limit real-time clinical utility. Here, we sought to develop a brain-age framework suitable for use during routine clinical head MRI examinations. Using a deep learning-based neuroradiology report classifier, we generated a dataset of 23,302 'radiologically normal for age' head MRI examinations from two large UK hospitals for model training and testing (age range = 18-95 years), and demonstrate fast (< 5 s), accurate (mean absolute error [MAE] < 4 years) age prediction from clinical-grade, minimally processed axial T2-weighted and axial diffusion-weighted scans, with generalisability between hospitals and scanner vendors (Δ MAE < 1 year). The clinical relevance of these brain-age predictions was tested using 228 patients whose MRIs were reported independently by neuroradiologists as showing atrophy 'excessive for age'. These patients had systematically higher brain-predicted age than chronological age (mean predicted age difference = +5.89 years, 'radiologically normal for age' mean predicted age difference = +0.05 years, p < 0.0001). Our brain-age framework demonstrates feasibility for use as a screening tool during routine hospital examinations to automatically detect older-appearing brains in real-time, with relevance for clinical decision-making and optimising patient pathways.

Mental health in the UK Biobank: A roadmap to self-report measures and neuroimaging correlates

The UK Biobank (UKB) is a highly promising dataset for brain biomarker research into population mental health due to its unprecedented sample size and extensive phenotypic, imaging, and biological measurements. In this study, we aimed to provide a shared foundation for UKB neuroimaging research into mental health with a focus on anxiety and depression. We compared UKB self-report measures and revealed important timing effects between scan acquisition and separate online acquisition of some mental health measures. To overcome these timing effects, we introduced and validated the Recent Depressive Symptoms (RDS-4) score which we recommend for state-dependent and longitudinal research in the UKB. We furthermore tested univariate and multivariate associations between brain imaging-derived phenotypes (IDPs) and mental health. Our results showed a significant multivariate relationship between IDPs and mental health, which was replicable. Conversely, effect sizes for individual IDPs were small. Test-retest reliability of IDPs was stronger for measures of brain structure than for measures of brain function. Taken together, these results provide benchmarks and guidelines for future UKB research into brain biomarkers of mental health.

Reliability of structural MRI measurements: The effects of scan session, head tilt, inter-scan interval, acquisition sequence, FreeSurfer version and processing stream

BACKGROUND

Large-scale longitudinal and multi-centre studies are used to explore neuroimaging markers of normal ageing, and neurodegenerative and mental health disorders. Longitudinal changes in brain structure are typically small, therefore the reliability of automated techniques is crucial. Determining the effects of different factors on reliability allows investigators to control those adversely affecting reliability, calculate statistical power, or even avoid particular brain measures with low reliability. This study examined the impact of several image acquisition and processing factors and documented the test-retest reliability of structural MRI measurements.

METHODS

In Phase I, 20 healthy adults (11 females; aged 20-30 years) were scanned on two occasions three weeks apart on the same scanner using the ADNI-3 protocol. On each occasion, individuals were scanned twice (repetition), after re-entering the scanner (reposition) and after tilting their head forward. At one year follow-up, nine returning individuals and 11 new volunteers were recruited for Phase II (11 females; aged 22-31 years). Scans were acquired on two different scanners using the ADNI-2 and ADNI-3 protocols. Structural images were processed using FreeSurfer (v5.3.0, 6.0.0 and 7.1.0) to provide subcortical and cortical volume, cortical surface area and thickness measurements. Intra-class correlation coefficients (ICC) were calculated to estimate test-retest reliability. We examined the effect of repetition, reposition, head tilt, time between scans, MRI sequence and scanner on reliability of structural brain measurements. Mean percentage differences were also calculated in supplementary analyses.

RESULTS

Using the FreeSurfer v7.1.0 longitudinal pipeline, we observed high reliability for subcortical and cortical volumes, and cortical surface areas at repetition, reposition, three weeks and one year (mean ICCs>0.97). Cortical thickness reliability was lower (mean ICCs>0.82). Head tilt had the greatest adverse impact on ICC estimates, for example reducing mean right cortical thickness to ICC=0.74. In contrast, changes in ADNI sequence or MRI scanner had a minimal effect. We observed an increase in reliability for updated FreeSurfer versions, with the longitudinal pipeline consistently having a higher reliability than the cross-sectional pipeline.

DISCUSSION

Longitudinal studies should monitor or control head tilt to maximise reliability. We provided the ICC estimates and mean percentage differences for all FreeSurfer brain regions, which may inform power analyses for clinical studies and have implications for the design of future longitudinal studies.

Testing a convolutional neural network-based hippocampal segmentation method in a stroke population

As stroke mortality rates decrease, there has been a surge of effort to study poststroke dementia (PSD) to improve long-term quality of life for stroke survivors. Hippocampal volume may be an important neuroimaging biomarker in poststroke dementia, as it has been associated with many other forms of dementia. However, studying hippocampal volume using MRI requires hippocampal segmentation. Advances in automated segmentation methods have allowed for studying the hippocampus on a large scale, which is important for robust results in the heterogeneous stroke population. However, most of these automated methods use a single atlas-based approach and may fail in the presence of severe structural abnormalities common in stroke. Hippodeep, a new convolutional neural network-based hippocampal segmentation method, does not rely solely on a single atlas-based approach and thus may be better suited for stroke populations. Here, we compared quality control and the accuracy of segmentations generated by Hippodeep and two well-accepted hippocampal segmentation methods on stroke MRIs (FreeSurfer 6.0 whole hippocampus and FreeSurfer 6.0 sum of hippocampal subfields). Quality control was performed using a stringent protocol for visual inspection of the segmentations, and accuracy was measured as volumetric correlation with manual segmentations. Hippodeep performed significantly better than both FreeSurfer methods in terms of quality control. All three automated segmentation methods had good correlation with manual segmentations and no one method was significantly more correlated than the others. Overall, this study suggests that both Hippodeep and FreeSurfer may be useful for hippocampal segmentation in stroke rehabilitation research, but Hippodeep may be more robust to stroke lesion anatomy.

Reproducibility in the absence of selective reporting: An illustration from large-scale brain asymmetry research

The problem of poor reproducibility of scientific findings has received much attention over recent years, in a variety of fields including psychology and neuroscience. The problem has been partly attributed to publication bias and unwanted practices such as p-hacking. Low statistical power in individual studies is also understood to be an important factor. In a recent multisite collaborative study, we mapped brain anatomical left-right asymmetries for regional measures of surface area and cortical thickness, in 99 MRI datasets from around the world, for a total of over 17,000 participants. In the present study, we revisited these hemispheric effects from the perspective of reproducibility. Within each dataset, we considered that an effect had been reproduced when it matched the meta-analytic effect from the 98 other datasets, in terms of effect direction and significance threshold. In this sense, the results within each dataset were viewed as coming from separate studies in an "ideal publishing environment," that is, free from selective reporting and p hacking. We found an average reproducibility rate of 63.2% (SD = 22.9%, min = 22.2%, max = 97.0%). As expected, reproducibility was higher for larger effects and in larger datasets. Reproducibility was not obviously related to the age of participants, scanner field strength, FreeSurfer software version, cortical regional measurement reliability, or regional size. These findings constitute an empirical illustration of reproducibility in the absence of publication bias or p hacking, when assessing realistic biological effects in heterogeneous neuroscience data, and given typically-used sample sizes.

Mega-analysis methods in ENIGMA: The experience of the generalized anxiety disorder working group

The ENIGMA group on Generalized Anxiety Disorder (ENIGMA-Anxiety/GAD) is part of a broader effort to investigate anxiety disorders using imaging and genetic data across multiple sites worldwide. The group is actively conducting a mega-analysis of a large number of brain structural scans. In this process, the group was confronted with many methodological challenges related to study planning and implementation, between-country transfer of subject-level data, quality control of a considerable amount of imaging data, and choices related to statistical methods and efficient use of resources. This report summarizes the background information and rationale for the various methodological decisions, as well as the approach taken to implement them. The goal is to document the approach and help guide other research groups working with large brain imaging data sets as they develop their own analytic pipelines for mega-analyses.

Ironsmith: An Automated Pipeline for QSM-based Data Analyses

Quantitative susceptibility mapping (QSM) is an MRI-based, computational method for anatomically localizing and measuring concentrations of specific biomarkers in tissue such as iron. Growing research suggests QSM is a viable method for evaluating the impact of iron overload in neurological disorders and on cognitive performance in aging. Several software toolboxes are currently available to reconstruct QSM maps from 3D GRE MR Images. However, few if any software packages currently exist that offer fully automated pipelines for QSM-based data analyses: from DICOM images to region-of-interest (ROI) based QSM values. Even less QSM-based software exist that offer quality control measures for evaluating the QSM output. Here, we address these gaps in the field by introducing and demonstrating the reliability and external validity of Ironsmith; an open-source, fully automated pipeline for creating and processing QSM maps, extracting QSM values from subcortical and cortical brain regions (89 ROIs) and evaluating the quality of QSM data using SNR measures and assessment of outlier regions on phase images. Ironsmith also features automatic filtering of QSM outlier values and precise CSF-only QSM reference masks that minimize partial volume effects. Testing of Ironsmith revealed excellent intra- and inter-rater reliability. Finally, external validity of Ironsmith was demonstrated via an anatomically selective relationship between motor performance and Ironsmith-derived QSM values in motor cortex. In sum, Ironsmith provides a freely-available, reliable, turn-key pipeline for QSM-based data analyses to support research on the impact of brain iron in aging and neurodegenerative disease.

Beyond the average patient: how neuroimaging models can address heterogeneity in dementia

Dementia is a highly heterogeneous condition, with pronounced individual differences in age of onset, clinical presentation, progression rates and neuropathological hallmarks, even within a specific diagnostic group. However, the most common statistical designs used in dementia research studies and clinical trials overlook this heterogeneity, instead relying on comparisons of group average differences (e.g. patient versus control or treatment versus placebo), implicitly assuming within-group homogeneity. This one-size-fits-all approach potentially limits our understanding of dementia aetiology, hindering the identification of effective treatments. Neuroimaging has enabled the characterization of the average neuroanatomical substrates of dementias; however, the increasing availability of large open neuroimaging datasets provides the opportunity to examine patterns of neuroanatomical variability in individual patients. In this update, we outline the causes and consequences of heterogeneity in dementia and discuss recent research that aims to tackle heterogeneity directly, rather than assuming that dementia affects everyone in the same way. We introduce spatial normative modelling as an emerging data-driven technique, which can be applied to dementia data to model neuroanatomical variation, capturing individualized neurobiological 'fingerprints'. Such methods have the potential to detect clinically relevant subtypes, track an individual's disease progression or evaluate treatment responses, with the goal of moving towards precision medicine for dementia.

Comparison of B-Scan Ultrasound and MRI-Based Optic Nerve Sheath Diameter (ONSD) Measurements in Children

BACKGROUND

Qualitative, noninvasive assessment of intracranial pressure is of eminent importance in pediatric patients in many clinical situations and can reliably be performed using transorbital ultrasonographic measurement of the optic nerve sheath diameter (ONSD). MRI-based determination of ONSD can serve as an alternative if ultrasound (US) is not possible or available for various reasons, for example, in small, incompliant children. This study investigates repeatability and observer reliability of US ONSD and correlation and bias of US- versus MRI-based ONSD assessment in pediatric patients.

METHODS

One hundred fifty children diagnosed with tumor (n = 40), hydrocephalus (n = 42), and other cranial pathologies (n = 68) were included. Bilateral ONSD was quantified by US using a 12-MHz linear array transducer. This was compared with ONSD measured in simultaneously acquired (≤24 h) T2-weighted MRI scans of the orbit.

RESULTS

Repeatability of individual US values and intraobserver ONSD was outstanding (Cronbach's α = 0.984 and 0.996, respectively). Overall mean values for ONSD were 5.8 ± 0.88 mm and 5.7 ± 0.89 mm for US and MRI, respectively. Correlation between US and MRI-based ONSD was strong (r = 0.976, P < 0.01). Bland and Altman analysis showed a mean bias of 0.078 mm. A repeated-measures correlation (r_rm) in 9 patients showed an excellent value (r_rm = 0.94, P < 0.01).

CONCLUSIONS

Repeatability and reliability of US ONSD determination is excellent. In case US ONSD assessment is not possible or available, MRI scans can serve as an excellent alternative. The difference of US and MRI ONSD is minimal and insignificant, and thus, both techniques can complement each other.

A simple permutation-based test of intermodal correspondence

Many key findings in neuroimaging studies involve similarities between brain maps, but statistical methods used to measure these findings have varied. Current state-of-the-art methods involve comparing observed group-level brain maps (after averaging intensities at each image location across multiple subjects) against spatial null models of these group-level maps. However, these methods typically make strong and potentially unrealistic statistical assumptions, such as covariance stationarity. To address these issues, in this article we propose using subject-level data and a classical permutation testing framework to test and assess similarities between brain maps. Our method is comparable to traditional permutation tests in that it involves randomly permuting subjects to generate a null distribution of intermodal correspondence statistics, which we compare to an observed statistic to estimate a p-value. We apply and compare our method in simulated and real neuroimaging data from the Philadelphia Neurodevelopmental Cohort. We show that our method performs well for detecting relationships between modalities known to be strongly related (cortical thickness and sulcal depth), and it is conservative when an association would not be expected (cortical thickness and activation on the n-back working memory task). Notably, our method is the most flexible and reliable for localizing intermodal relationships within subregions of the brain and allows for generalizable statistical inference.

Evaluation of movement and brain activity

Clinical neurophysiology studies can contribute important information about the physiology of human movement and the pathophysiology and diagnosis of different movement disorders. Some techniques can be accomplished in a routine clinical neurophysiology laboratory and others require some special equipment. This review, initiating a series of articles on this topic, focuses on the methods and techniques. The methods reviewed include EMG, EEG, MEG, evoked potentials, coherence, accelerometry, posturography (balance), gait, and sleep studies. Functional MRI (fMRI) is also reviewed as a physiological method that can be used independently or together with other methods. A few applications to patients with movement disorders are discussed as examples, but the detailed applications will be the subject of other articles.

Voxel-wise and spatial modelling of binary lesion masks: Comparison of methods with a realistic simulation framework

OBJECTIVES

White matter lesions are a very common finding on MRI in older adults and their presence increases the risk of stroke and dementia. Accurate and computationally efficient modelling methods are necessary to map the association of lesion incidence with risk factors, such as hypertension. However, there is no consensus in the brain mapping literature whether a voxel-wise modelling approach is better for binary lesion data than a more computationally intensive spatial modelling approach that accounts for voxel dependence.

METHODS

We review three regression approaches for modelling binary lesion masks including mass-univariate probit regression modelling with either maximum likelihood estimates, or mean bias-reduced estimates, and spatial Bayesian modelling, where the regression coefficients have a conditional autoregressive model prior to account for local spatial dependence. We design a novel simulation framework of artificial lesion maps to compare the three alternative lesion mapping methods. The age effect on lesion probability estimated from a reference data set (13,680 individuals from the UK Biobank) is used to simulate a realistic voxel-wise distribution of lesions across age. To mimic the real features of lesion masks, we propose matching brain lesion summaries (total lesion volume, average lesion size and lesion count) across the reference data set and the simulated data sets. Thus, we allow for a fair comparison between the modelling approaches, under a realistic simulation setting.

RESULTS

Our findings suggest that bias-reduced estimates for voxel-wise binary-response generalized linear models (GLMs) overcome the drawbacks of infinite and biased maximum likelihood estimates and scale well for large data sets because voxel-wise estimation can be performed in parallel across voxels. Contrary to the assumption of spatial dependence being key in lesion mapping, our results show that voxel-wise bias-reduction and spatial modelling result in largely similar estimates.

CONCLUSIONS

Bias-reduced estimates for voxel-wise GLMs are not only accurate but also computationally efficient, which will become increasingly important as more biobank-scale neuroimaging data sets become available.

Minimal specifications for non-human primate MRI: Challenges in standardizing and harmonizing data collection

Recent methodological advances in MRI have enabled substantial growth in neuroimaging studies of non-human primates (NHPs), while open data-sharing through the PRIME-DE initiative has increased the availability of NHP MRI data and the need for robust multi-subject multi-center analyses. Streamlined acquisition and analysis protocols would accelerate and improve these efforts. However, consensus on minimal standards for data acquisition protocols and analysis pipelines for NHP imaging remains to be established, particularly for multi-center studies. Here, we draw parallels between NHP and human neuroimaging and provide minimal guidelines for harmonizing and standardizing data acquisition. We advocate robust translation of widely used open-access toolkits that are well established for analyzing human data. We also encourage the use of validated, automated pre-processing tools for analyzing NHP data sets. These guidelines aim to refine methodological and analytical strategies for small and large-scale NHP neuroimaging data. This will improve reproducibility of results, and accelerate the convergence between NHP and human neuroimaging strategies which will ultimately benefit fundamental and translational brain science.

Strengths and challenges of longitudinal non-human primate neuroimaging

Longitudinal non-human primate neuroimaging has the potential to greatly enhance our understanding of primate brain structure and function. Here we describe its specific strengths, compared to both cross-sectional non-human primate neuroimaging and longitudinal human neuroimaging, but also its associated challenges. We elaborate on factors guiding the use of different analytical tools, subject-specific versus age-specific templates for analyses, and issues related to statistical power.

Beware of white matter hyperintensities causing systematic errors in FreeSurfer gray matter segmentations!

Volumetric estimates of subcortical and cortical structures, extracted from T1-weighted MRIs, are widely used in many clinical and research applications. Here, we investigate the impact of the presence of white matter hyperintensities (WMHs) on FreeSurfer gray matter (GM) structure volumes and its possible bias on functional relationships. T1-weighted images from 1,077 participants (4,321 timepoints) from the Alzheimer's Disease Neuroimaging Initiative were processed with FreeSurfer version 6.0.0. WMHs were segmented using a previously validated algorithm on either T2-weighted or Fluid-attenuated inversion recovery images. Mixed-effects models were used to assess the relationships between overlapping WMHs and GM structure volumes and overall WMH burden, as well as to investigate whether such overlaps impact associations with age, diagnosis, and cognitive performance. Participants with higher WMH volumes had higher overlaps with GM volumes of bilateral caudate, cerebral cortex, putamen, thalamus, pallidum, and accumbens areas (p < .0001). When not corrected for WMHs, caudate volumes increased with age (p < .0001) and were not different between cognitively healthy individuals and age-matched probable Alzheimer's disease patients. After correcting for WMHs, caudate volumes decreased with age (p < .0001), and Alzheimer's disease patients had lower caudate volumes than cognitively healthy individuals (p < .01). Uncorrected caudate volume was not associated with ADAS13 scores, whereas corrected lower caudate volumes were significantly associated with poorer cognitive performance (p < .0001). Presence of WMHs leads to systematic inaccuracies in GM segmentations, particularly for the caudate, which can also change clinical associations. While specifically measured for the Freesurfer toolkit, this problem likely affects other algorithms.

Replicating extensive brain structural heterogeneity in individuals with schizophrenia and bipolar disorder

Identifying brain processes involved in the risk and development of mental disorders is a major aim. We recently reported substantial interindividual heterogeneity in brain structural aberrations among patients with schizophrenia and bipolar disorder. Estimating the normative range of voxel-based morphometry (VBM) data among healthy individuals using a Gaussian process regression (GPR) enables us to map individual deviations from the healthy range in unseen datasets. Here, we aim to replicate our previous results in two independent samples of patients with schizophrenia (n1 = 94; n2 = 105), bipolar disorder (n1 = 116; n2 = 61), and healthy individuals (n1 = 400; n2 = 312). In line with previous findings with exception of the cerebellum our results revealed robust group level differences between patients and healthy individuals, yet only a small proportion of patients with schizophrenia or bipolar disorder exhibited extreme negative deviations from normality in the same brain regions. These direct replications support that group level-differences in brain structure disguise considerable individual differences in brain aberrations, with important implications for the interpretation and generalization of group-level brain imaging findings to the individual with a mental disorder.

Beyond linearity in neuroimaging: Capturing nonlinear relationships with application to longitudinal studies

The ubiquitous adoption of linearity for quantitative predictors in statistical modeling is likely attributable to its advantages of straightforward interpretation and computational feasibility. The linearity assumption may be a reasonable approximation especially when the variable is confined within a narrow range, but it can be problematic when the variable's effect is non-monotonic or complex. Furthermore, visualization and model assessment of a linear fit are usually omitted because of challenges at the whole brain level in neuroimaging. By adopting a principle of learning from the data in the presence of uncertainty to resolve the problematic aspects of conventional polynomial fitting, we introduce a flexible and adaptive approach of multilevel smoothing splines (MSS) to capture any nonlinearity of a quantitative predictor for population-level neuroimaging data analysis. With no prior knowledge regarding the underlying relationship other than a parsimonious assumption about the extent of smoothness (e.g., no sharp corners), we express the unknown relationship with a sufficient number of smoothing splines and use the data to adaptively determine the specifics of the nonlinearity. In addition to introducing the theoretical framework of MSS as an efficient approach with a counterbalance between flexibility and stability, we strive to (a) lay out the specific schemes for population-level nonlinear analyses that may involve task (e.g., contrasting conditions) and subject-grouping (e.g., patients vs controls) factors; (b) provide modeling accommodations to adaptively reveal, estimate and compare any nonlinear effects of a predictor across the brain, or to more accurately account for the effects (including nonlinear effects) of a quantitative confound; (c) offer the associated program 3dMSS to the neuroimaging community for whole-brain voxel-wise analysis as part of the AFNI suite; and (d) demonstrate the modeling approach and visualization processes with a longitudinal dataset of structural MRI scans.

The variability of MR axon radii estimates in the human white matter

The noninvasive quantification of axonal morphology is an exciting avenue for gaining understanding of the function and structure of the central nervous system. Accurate non-invasive mapping of micron-sized axon radii using commonly applied neuroimaging techniques, that is, diffusion-weighted MRI, has been bolstered by recent hardware developments, specifically MR gradient design. Here the whole brain characterization of the effective MR axon radius is presented and the inter- and intra-scanner test-retest repeatability and reproducibility are evaluated to promote the further development of the effective MR axon radius as a neuroimaging biomarker. A coefficient-of-variability of approximately 10% in the voxelwise estimation of the effective MR radius is observed in the test-retest analysis, but it is shown that the performance can be improved fourfold using a customized along-tract analysis.

The Use, Standardization, and Interpretation of Brain Imaging Data in Clinical Trials of Neurodegenerative Disorders

Imaging biomarkers play a wide-ranging role in clinical trials for neurological disorders. This includes selecting the appropriate trial participants, establishing target engagement and mechanism-related pharmacodynamic effect, monitoring safety, and providing evidence of disease modification. In the early stages of clinical drug development, evidence of target engagement and/or downstream pharmacodynamic effect-especially with a clear relationship to dose-can provide confidence that the therapeutic candidate should be advanced to larger and more expensive trials, and can inform the selection of the dose(s) to be further tested, i.e., to "de-risk" the drug development program. In these later-phase trials, evidence that the therapeutic candidate is altering disease-related biomarkers can provide important evidence that the clinical benefit of the compound (if observed) is grounded in meaningful biological changes. The interpretation of disease-related imaging markers, and comparability across different trials and imaging tools, is greatly improved when standardized outcome measures are defined. This standardization should not impinge on scientific advances in the imaging tools per se but provides a common language in which the results generated by these tools are expressed. PET markers of pathological protein aggregates and structural imaging of brain atrophy are common disease-related elements across many neurological disorders. However, PET tracers for pathologies beyond amyloid β and tau are needed, and the interpretability of structural imaging can be enhanced by some simple considerations to guard against the possible confound of pseudo-atrophy. Learnings from much-studied conditions such as Alzheimer's disease and multiple sclerosis will be beneficial as the field embraces rarer diseases.

MarkVCID cerebral small vessel consortium: II. Neuroimaging protocols

The MarkVCID consortium was formed under cooperative agreements with the National Institute of Neurologic Diseases and Stroke (NINDS) and National Institute on Aging (NIA) in 2016 with the goals of developing and validating biomarkers for the cerebral small vessel diseases associated with the vascular contributions to cognitive impairment and dementia (VCID). Rigorously validated biomarkers have consistently been identified as crucial for multicenter studies to identify effective strategies to prevent and treat VCID, specifically to detect increased VCID risk, diagnose the presence of small vessel disease and its subtypes, assess prognosis for disease progression or response to treatment, demonstrate target engagement or mechanism of action for candidate interventions, and monitor disease progression during treatment. The seven project sites and central coordinating center comprising MarkVCID, working with NINDS and NIA, identified a panel of 11 candidate fluid- and neuroimaging-based biomarker kits and established harmonized multicenter study protocols (see companion paper "MarkVCID cerebral small vessel consortium: I. Enrollment, clinical, fluid protocols" for full details). Here we describe the MarkVCID neuroimaging protocols with specific focus on validating their application to future multicenter trials. MarkVCID procedures for participant enrollment; clinical and cognitive evaluation; and collection, handling, and instrumental validation of fluid samples are described in detail in a companion paper. Magnetic resonance imaging (MRI) has long served as the neuroimaging modality of choice for cerebral small vessel disease and VCID because of its sensitivity to a wide range of brain properties, including small structural lesions, connectivity, and cerebrovascular physiology. Despite MRI's widespread use in the VCID field, there have been relatively scant data validating the repeatability and reproducibility of MRI-based biomarkers across raters, scanner types, and time intervals (collectively defined as instrumental validity). The MRI protocols described here address the core MRI sequences for assessing cerebral small vessel disease in future research studies, specific sequence parameters for use across various research scanner types, and rigorous procedures for determining instrumental validity. Another candidate neuroimaging modality considered by MarkVCID is optical coherence tomography angiography (OCTA), a non-invasive technique for directly visualizing retinal capillaries as a marker of the cerebral capillaries. OCTA has theoretical promise as a unique opportunity to visualize small vessels derived from the cerebral circulation, but at a considerably earlier stage of development than MRI. The additional OCTA protocols described here address procedures for determining OCTA instrumental validity, evaluating sources of variability such as pupil dilation, and handling data to maintain participant privacy. MRI protocol and instrumental validation The core sequences selected for the MarkVCID MRI protocol are three-dimensional T1-weighted multi-echo magnetization-prepared rapid-acquisition-of-gradient-echo (ME-MPRAGE), three-dimensional T2-weighted fast spin echo fluid-attenuated-inversion-recovery (FLAIR), two-dimensional diffusion-weighted spin-echo echo-planar imaging (DWI), three-dimensional T2*-weighted multi-echo gradient echo (3D-GRE), three-dimensional T2 -weighted fast spin-echo imaging (T2w), and two-dimensional T2*-weighted gradient echo echo-planar blood-oxygenation-level-dependent imaging with brief periods of CO2 inhalation (BOLD-CVR). Harmonized parameters for each of these core sequences were developed for four 3 Tesla MRI scanner models in widespread use at academic medical centers. MarkVCID project sites are trained and certified for their instantiation of the consortium MRI protocols. Sites are required to perform image quality checks every 2 months using the Alzheimer's Disease Neuroimaging Initiative phantom. Instrumental validation for MarkVCID MRI-based biomarkers is operationally defined as inter-rater reliability, test-retest repeatability, and inter-scanner reproducibility. Assessments of these instrumental properties are performed on individuals representing a range of cerebral small vessel disease from mild to severe. Inter-rater reliability is determined by distribution of an independent dataset of MRI scans to each analysis site. Test-retest repeatability is determined by repeat MRI scans performed on individual participants on a single MRI scanner after a short (1- to 14-day) interval. Inter-scanner reproducibility is determined by repeat MRI scans performed on individuals performed across four MRI scanner models. OCTA protocol and instrumental validation The MarkVCID OCTA protocol uses a commercially available, Food and Drug Administration-approved OCTA apparatus. Imaging is performed on one dilated and one undilated eye to assess the need for dilation. Scans are performed in quadruplicate. MarkVCID project sites participating in OCTA validation are trained and certified by this biomarker's lead investigator. Inter-rater reliability for OCTA is assessed by distribution of OCTA datasets to each analysis site. Test-retest repeatability is assessed by repeat OCTA imaging on individuals on the same day as their baseline OCTA and a different-day repeat session after a short (1- to 14-day) interval. Methods were developed to allow the OCTA data to be de-identified by the sites before transmission to the central data management system. The MarkVCID neuroimaging protocols, like the other MarkVCID procedures, are designed to allow translation to multicenter trials and as a template for outside groups to generate directly comparable neuroimaging data. The MarkVCID neuroimaging protocols are available to the biomedical community and intended to be shared. In addition to the instrumental validation procedures described here, each of the neuroimaging MarkVCID kits will undergo biological validation to determine its ability to measure important aspects of VCID such as cognitive function. The analytic methods for the neuroimaging-based kits and the results of these validation studies will be published separately. The results will ultimately determine the neuroimaging kits' potential usefulness for multicenter interventional trials in small vessel disease-related VCID.

DeepMapi: a Fully Automatic Registration Method for Mesoscopic Optical Brain Images Using Convolutional Neural Networks

The extreme complexity of mammalian brains requires a comprehensive deconstruction of neuroanatomical structures. Scientists normally use a brain stereotactic atlas to determine the locations of neurons and neuronal circuits. However, different brain images are normally not naturally aligned even when they are imaged with the same setup, let alone under the differing resolutions and dataset sizes used in mesoscopic imaging. As a result, it is difficult to achieve high-throughput automatic registration without manual intervention. Here, we propose a deep learning-based registration method called DeepMapi to predict a deformation field used to register mesoscopic optical images to an atlas. We use a self-feedback strategy to address the problem of imbalanced training sets (sampling at a fixed step size in nonuniform brains of structures and deformations) and use a dual-hierarchical network to capture the large and small deformations. By comparing DeepMapi with other registration methods, we demonstrate its superiority over a set of ground truth images, including both optical and MRI images. DeepMapi achieves fully automatic registration of mesoscopic micro-optical images, even macroscopic MRI datasets, in minutes, with an accuracy comparable to those of manual annotations by anatomists.

Flexible annotation atlas of the mouse brain: combining and dividing brain structures of the Allen Brain Atlas while maintaining anatomical hierarchy

A brain atlas is necessary for analyzing structure and function in neuroimaging research. Although various annotation volumes (AVs) for the mouse brain have been proposed, it is common in magnetic resonance imaging (MRI) of the mouse brain that regions-of-interest (ROIs) for brain structures (nodes) are created arbitrarily according to each researcher's necessity, leading to inconsistent ROIs among studies. One reason for such a situation is the fact that earlier AVs were fixed, i.e. combination and division of nodes were not implemented. This report presents a pipeline for constructing a flexible annotation atlas (FAA) of the mouse brain by leveraging public resources of the Allen Institute for Brain Science on brain structure, gene expression, and axonal projection. A mere two-step procedure with user-specified, text-based information and Python codes constructs FAA with nodes which can be combined or divided objectively while maintaining anatomical hierarchy of brain structures. Four FAAs with total node count of 4, 101, 866, and 1381 were demonstrated. Unique characteristics of FAA realized analysis of resting-state functional connectivity (FC) across the anatomical hierarchy and among cortical layers, which were thin but large brain structures. FAA can improve the consistency of whole brain ROI definition among laboratories by fulfilling various requests from researchers with its flexibility and reproducibility.

Neuroimaging in Acute Stroke

PURPOSE OF REVIEW

This article describes how imaging can be used by physicians in diagnosing, determining prognosis, and making appropriate treatment decisions in a timely manner in patients with acute stroke.

RECENT FINDINGS

Advances in acute stroke treatment, including the use of endovascular thrombectomy in patients with large vessel occlusion and, more recently, of IV thrombolysis in an extended time window, have resulted in a paradigm shift in how imaging is used in patients with acute stroke. This paradigm shift, combined with the understanding that "time is brain," means that imaging must be fast, reliable, and available around the clock for physicians to make appropriate clinical decisions. CT has therefore become the primary imaging modality of choice. Recognition of a large vessel occlusion using CT angiography has become essential in identifying patients for endovascular thrombectomy, and techniques such as imaging collaterals on CT angiography or measuring blood flow to predict tissue fate using CT perfusion have become useful tools in selecting patients for acute stroke therapy. Understanding the use of these imaging modalities and techniques in dealing with an emergency such as acute stroke has therefore become more important than ever for physicians treating patients with acute stroke.

SUMMARY

Imaging the brain and the blood vessels supplying it using modern tools and techniques is a key step in understanding the pathophysiology of acute stroke and making appropriate and timely clinical decisions.

Diagnostic accuracy of MR planimetry in clinically unclassifiable parkinsonism

INTRODUCTION

Quantitative MR planimetric measurements were reported to discriminate between progressive supranuclear palsy (PSP) and non-PSP parkinsonism, yet few data exist on the usefulness of these markers in early disease stages.

METHODS

The pons-to-midbrain area ratio (P/M) and the Magnetic Resonance Parkinsonism Index (MRPI) as well as new indices, termed P/M2.0 and MRPI2.0, multiplying the former by a ratio of the third ventricle (3rdV) width/frontal horns (FH) width, were calculated on T1-weighted images in 84 patients with clinically unclassifiable neurodegenerative parkinsonism (CUP) at the time of imaging. Areas under the curve (AUCs) of these markers for predicting future PSP was determined. The final clinical diagnosis was made after at least 24 months of follow-up.

RESULTS

Final diagnosis was Parkinson's disease in 55 patients, multiple system atrophy in 12 cases, and PSP in 17. At baseline imaging, patients with a final PSP diagnosis had significantly higher MRPI, P/M, MRPI2.0 and P/M2.0 values compared to the other groups. AUCs in discriminating between future PSP and non-PSP parkinsonism were 0.91 for both the P/M and the MRPI and 0.98 for the P/M2.0 and the MRPI2.0.

CONCLUSIONS

Brainstem-derived MR planimetric measures yield high diagnostic accuracy for separating PSP from non-PSP parkinsonism in early disease stages when clinical criteria are not yet fully met. Consistent with the underlying pathology in PSP, our study suggests that inclusion of 3^rdV width makes P/M2.0 and MRPI2.0 more accurate in diagnosing early stage PSP patients than the P/M and MRPI.

Harmonization of diffusion MRI data sets with adaptive dictionary learning

Diffusion magnetic resonance imaging can indirectly infer the microstructure of tissues and provide metrics subject to normal variability in a population. Potentially abnormal values may yield essential information to support analysis of controls and patients cohorts, but subtle confounds could be mistaken for purely biologically driven variations amongst subjects. In this work, we propose a new harmonization algorithm based on adaptive dictionary learning to mitigate the unwanted variability caused by different scanner hardware while preserving the natural biological variability of the data. Our harmonization algorithm does not require paired training data sets, nor spatial registration or matching spatial resolution. Overcomplete dictionaries are learned iteratively from all data sets at the same time with an adaptive regularization criterion, removing variability attributable to the scanners in the process. The obtained mapping is applied directly in the native space of each subject toward a scanner-space. The method is evaluated with a public database which consists of two different protocols acquired on three different scanners. Results show that the effect size of the four studied diffusion metrics is preserved while removing variability attributable to the scanner. Experiments with alterations using a free water compartment, which is not simulated in the training data, shows that the modifications applied to the diffusion weighted images are preserved in the diffusion metrics after harmonization, while still reducing global variability at the same time. The algorithm could help multicenter studies pooling their data by removing scanner specific confounds, and increase statistical power in the process.

The Preoperative Evaluation of Drug-Resistant Epilepsy

Antiepileptic drugs afford good seizure control for approximately 70% of individuals with epilepsy. Epilepsy surgery is extremely helpful for appropriate individuals with drug resistance. Since antiquity, trephination was a crude and invasive technique to manage epilepsy. The late 1800s saw the advent of a more evidence-based approach with attempts to define seizure foci and determine areas of function. Seizure localization initially required direct brain stimulation during surgery before resection. Fortunately, improved knowledge of seizure semiology and advancements in preoperative investigations have enabled epilepsy specialists to better analyze the benefit of seizure reduction versus risk of functional harm. This preoperative phase and the investigative techniques used to analyze surgical candidacy will be discussed in this article.

Finding specificity in structural brain alterations through Bayesian reverse inference

In the field of neuroimaging reverse inferences can lead us to suppose the involvement of cognitive processes from certain patterns of brain activity. However, the same reasoning holds if we substitute "brain activity" with "brain alteration" and "cognitive process" with "brain disorder." The fact that different brain disorders exhibit a high degree of overlap in their patterns of structural alterations makes forward inference-based analyses less suitable for identifying brain areas whose alteration is specific to a certain pathology. In the forward inference-based analyses, in fact, it is impossible to distinguish between areas that are altered by the majority of brain disorders and areas that are specifically affected by certain diseases. To address this issue and allow the identification of highly pathology-specific altered areas we used the Bayes' factor technique, which was employed, as a proof of concept, on voxel-based morphometry data of schizophrenia and Alzheimer's disease. This technique allows to calculate the ratio between the likelihoods of two alternative hypotheses (in our case, that the alteration of the voxel is specific for the brain disorder under scrutiny or that the alteration is not specific). We then performed temporal simulations of the alterations' spread associated with different pathologies. The Bayes' factor values calculated on these simulated data were able to reveal that the areas, which are more specific to a certain disease, are also the ones to be early altered. This study puts forward a new analytical instrument capable of innovating the methodological approach to the investigation of brain pathology.

Automated classification of depression from structural brain measures across two independent community-based cohorts

Major depressive disorder (MDD) has been the subject of many neuroimaging case-control classification studies. Although some studies report accuracies ≥80%, most have investigated relatively small samples of clinically-ascertained, currently symptomatic cases, and did not attempt replication in larger samples. We here first aimed to replicate previously reported classification accuracies in a small, well-phenotyped community-based group of current MDD cases with clinical interview-based diagnoses (from STratifying Resilience and Depression Longitudinally cohort, 'STRADL'). We performed a set of exploratory predictive classification analyses with measures related to brain morphometry and white matter integrity. We applied three classifier types-SVM, penalised logistic regression or decision tree-either with or without optimisation, and with or without feature selection. We then determined whether similar accuracies could be replicated in a larger independent population-based sample with self-reported current depression (UK Biobank cohort). Additional analyses extended to lifetime MDD diagnoses-remitted MDD in STRADL, and lifetime-experienced MDD in UK Biobank. The highest cross-validation accuracy (75%) was achieved in the initial current MDD sample with a decision tree classifier and cortical surface area features. The most frequently selected decision tree split variables included surface areas of bilateral caudal anterior cingulate, left lingual gyrus, left superior frontal, right precentral and paracentral regions. High accuracy was not achieved in the larger samples with self-reported current depression (53.73%), with remitted MDD (57.48%), or with lifetime-experienced MDD (52.68-60.29%). Our results indicate that high predictive classification accuracies may not immediately translate to larger samples with broader criteria for depression, and may not be robust across different classification approaches.

Deep learning-guided joint attenuation and scatter correction in multitracer neuroimaging studies

PET attenuation correction (AC) on systems lacking CT/transmission scanning, such as dedicated brain PET scanners and hybrid PET/MRI, is challenging. Direct AC in image-space, wherein PET images corrected for attenuation and scatter are synthesized from nonattenuation corrected PET (PET-nonAC) images in an end-to-end fashion using deep learning approaches (DLAC) is evaluated for various radiotracers used in molecular neuroimaging studies. One hundred eighty brain PET scans acquired using 18 F-FDG, 18 F-DOPA, 18 F-Flortaucipir (targeting tau pathology), and 18 F-Flutemetamol (targeting amyloid pathology) radiotracers (40 + 5, training/validation + external test, subjects for each radiotracer) were included. The PET data were reconstructed using CT-based AC (CTAC) to generate reference PET-CTAC and without AC to produce PET-nonAC images. A deep convolutional neural network was trained to generate PET attenuation corrected images (PET-DLAC) from PET-nonAC. The quantitative accuracy of this approach was investigated separately for each radiotracer considering the values obtained from PET-CTAC images as reference. A segmented AC map (PET-SegAC) containing soft-tissue and background air was also included in the evaluation. Quantitative analysis of PET images demonstrated superior performance of the DLAC approach compared to SegAC technique for all tracers. Despite the relatively low quantitative bias observed when using the DLAC approach, this approach appears vulnerable to outliers, resulting in noticeable local pseudo uptake and false cold regions. Direct AC in image-space using deep learning demonstrated quantitatively acceptable performance with less than 9% absolute SUV bias for the four different investigated neuroimaging radiotracers. However, this approach is vulnerable to outliers which result in large local quantitative bias.

Response assessment in diffuse intrinsic pontine glioma: recommendations from the Response Assessment in Pediatric Neuro-Oncology (RAPNO) working group

Optimising the conduct of clinical trials for diffuse intrinsic pontine glioma involves use of consistent, objective disease assessments and standardised response criteria. The Response Assessment in Pediatric Neuro-Oncology working group, consisting of an international panel of paediatric and adult neuro-oncologists, clinicians, radiologists, radiation oncologists, and neurosurgeons, was established to address issues and unique challenges in assessing response in children with CNS tumours. A working group was formed specifically to address response assessment in children and young adults with diffuse intrinsic pontine glioma and to develop a consensus on recommendations for response assessment. Response should be assessed using MRI of brain and spine, neurological examination, and anti-inflammatory or antiangiogenic drugs. Clinical imaging standards are defined. As with previous consensus recommendations, these recommendations will need to be validated in prospective clinical trials.

Longitudinal validity of PET-based staging of regional amyloid deposition

Positron emission tomography (PET)-based staging of regional amyloid deposition has recently emerged as a promising tool for sensitive detection and stratification of pathology progression in Alzheimer's Disease (AD). Here we present an updated methodological framework for PET-based amyloid staging using region-specific amyloid-positivity thresholds and assess its longitudinal validity using serial PET acquisitions. We defined region-specific thresholds of amyloid-positivity based on Florbetapir-PET data of 13 young healthy individuals (age ≤ 45y), applied these thresholds to Florbetapir-PET data of 179 cognitively normal older individuals to estimate a regional amyloid staging model, and tested this model in a larger sample of patients with mild cognitive impairment (N = 403) and AD dementia (N = 85). 2-year follow-up Florbetapir-PET scans from a subset of this sample (N = 436) were used to assess the longitudinal validity of the cross-sectional model based on individual stage transitions and data-driven longitudinal trajectory modeling. Results show a remarkable congruence between cross-sectionally estimated and longitudinally modeled trajectories of amyloid accumulation, beginning in anterior temporal areas, followed by frontal and medial parietal areas, the remaining associative neocortex, and finally primary sensory-motor areas and subcortical regions. Over 98% of individual amyloid deposition profiles and longitudinal stage transitions adhered to this staging scheme of regional pathology progression, which was further supported by corresponding changes in cerebrospinal fluid biomarkers. In conclusion, we provide a methodological refinement and longitudinal validation of PET-based staging of regional amyloid accumulation, which may help improving early detection and in-vivo stratification of pathologic disease progression in AD.

Fidelity imposed network edit (FINE) for solving ill-posed image reconstruction

Deep learning (DL) is increasingly used to solve ill-posed inverse problems in medical imaging, such as reconstruction from noisy and/or incomplete data, as DL offers advantages over conventional methods that rely on explicit image features and hand engineered priors. However, supervised DL-based methods may achieve poor performance when the test data deviates from the training data, for example, when it has pathologies not encountered in the training data. Furthermore, DL-based image reconstructions do not always incorporate the underlying forward physical model, which may improve performance. Therefore, in this work we introduce a novel approach, called fidelity imposed network edit (FINE), which modifies the weights of a pre-trained reconstruction network for each case in the testing dataset. This is achieved by minimizing an unsupervised fidelity loss function that is based on the forward physical model. FINE is applied to two important inverse problems in neuroimaging: quantitative susceptibility mapping (QSM) and under-sampled image reconstruction in MRI. Our experiments demonstrate that FINE can improve reconstruction accuracy.

A field-monitoring-based approach for correcting eddy-current-induced artifacts of up to the 2nd spatial order in human-connectome-project-style multiband diffusion MRI experiment at 7T: A pilot study

Over the recent years, significant advances in Spin-Echo (SE) Echo-Planar (EP) Diffusion MRI (dMRI) have enabled improved fiber tracking conspicuity in the human brain. At the same time, pushing the spatial resolution and using higher b-values inherently expose the acquired images to further eddy-current-induced distortion and blurring. Recently developed data-driven correction techniques, capable of significantly mitigating these defects, are included in the reconstruction pipelines developed for the Human Connectome Project (HCP) driven by the NIH BRAIN initiative. In this case, however, corrections are derived from the original diffusion-weighted (DW) magnitude images affected by distortion and blurring. Considering the complexity of k-space deviations in the presence of time varying high spatial order eddy currents, distortion and blurring may not be fully reversed when relying on magnitude DW images only. An alternative approach, consisting of iteratively reconstructing DW images based on the actual magnetic field spatiotemporal evolution measured with a magnetic field monitoring camera, has been successfully implemented at 3T in single band dMRI (Wilm et al., 2017, 2015). In this study, we aim to demonstrate the efficacy of this eddy current correction method in the challenging context of HCP-style multiband (MB ​= ​2) dMRI protocol. The magnetic field evolution was measured during the EP-dMRI readout echo train with a field monitoring camera equipped with 16 ¹⁹F NMR probes. The time variation of 0^th, 1^st and 2^nd order spherical field harmonics were used to reconstruct DW images. Individual DW images reconstructed with and without field correction were compared. The impact of eddy current correction was evaluated by comparing the corresponding direction-averaged DW images and fractional anisotropy (FA) maps. ¹⁹F field monitoring data confirmed the existence of significant field deviations induced by the diffusion-encoding gradients, with variations depending on diffusion gradient amplitude and direction. In DW images reconstructed with the field correction, residual aliasing artifacts were reduced or eliminated, and when high b-values were applied, better gray/white matter delineation and sharper gyri contours were observed, indicating reduced signal blurring. The improvement in image quality further contributed to sharper contours and better gray/white matter delineation in mean DW images and FA maps. In conclusion, we demonstrate that up-to-2^nd-order-eddy-current-induced field perturbation in multiband, in-plane accelerated HCP-style dMRI acquisition at 7T can be corrected by integrating the measured field evolution in image reconstruction.

Selective inversion recovery quantitative magnetization transfer imaging: Toward a 3 T clinical application in multiple sclerosis

BACKGROUND

Assessing the degree of myelin injury in patients with multiple sclerosis (MS) is challenging due to the lack of magnetic resonance imaging (MRI) methods specific to myelin quantity. By measuring distinct tissue parameters from a two-pool model of the magnetization transfer (MT) effect, quantitative magnetization transfer (qMT) may yield these indices. However, due to long scan times, qMT has not been translated clinically.

OBJECTIVES

We aim to assess the clinical feasibility of a recently optimized selective inversion recovery (SIR) qMT and to test the hypothesis that SIR-qMT-derived metrics are informative of radiological and clinical disease-related changes in MS.

METHODS

A total of 18 MS patients and 9 age- and sex-matched healthy controls (HCs) underwent a 3.0 Tesla (3 T) brain MRI, including clinical scans and an optimized SIR-qMT protocol. Four subjects were re-scanned at a 2-week interval to determine inter-scan variability.

RESULTS

SIR-qMT measures differed between lesional and non-lesional tissue (p < 0.0001) and between normal-appearing white matter (NAWM) of patients with more advanced disability and normal white matter (WM) of HCs (p < 0.05). SIR-qMT measures were associated with lesion volumes, disease duration, and disability scores (p ⩽ 0.002).

CONCLUSION

SIR-qMT at 3 T is clinically feasible and predicts both radiological and clinical disease severity in MS.

Towards Precision Medicine in Psychosis: Benefits and Challenges of Multimodal Multicenter Studies-PSYSCAN: Translating Neuroimaging Findings From Research into Clinical Practice

In the last 2 decades, several neuroimaging studies investigated brain abnormalities associated with the early stages of psychosis in the hope that these could aid the prediction of onset and clinical outcome. Despite advancements in the field, neuroimaging has yet to deliver. This is in part explained by the use of univariate analytical techniques, small samples and lack of statistical power, lack of external validation of potential biomarkers, and lack of integration of nonimaging measures (eg, genetic, clinical, cognitive data). PSYSCAN is an international, longitudinal, multicenter study on the early stages of psychosis which uses machine learning techniques to analyze imaging, clinical, cognitive, and biological data with the aim of facilitating the prediction of psychosis onset and outcome. In this article, we provide an overview of the PSYSCAN protocol and we discuss benefits and methodological challenges of large multicenter studies that employ neuroimaging measures.

Value of thrombus imaging in predicting the outcomes of patients with large-vessel occlusive strokes after endovascular therapy

BACKGROUND

Acute ischemic stroke leads to serious long-term disability and high mortality, especially in patients with large-vessel occlusive strokes. Nowadays, endovascular therapy is considered as an alternative treatment for these patients. Several studies have used thrombus characteristics based on non-contrast computed tomography (NCCT) and computed tomography angiography (CTA) to predict prognosis in ischemic stroke. We conducted a systematic review to identify potential imaging predictive factors for successful recanalization and improved clinical outcome after endovascular therapy in patients with large-vessel occlusion (LVO) in anterior arterial circulation.

METHODS

The PubMed databases were searched for related studies reported between September 18, 2009, and September 18, 2019.

RESULTS

We selected 11 studies on revascularization and 12 studies on clinical outcome. Patients with thrombus of higher Hounsfield unit (HU), shorter length, higher clot burden score, and increased thrombus permeability may achieve higher recanalization and improved clinical outcome, but the matter is still under debate.

CONCLUSION

Imaging of thrombus can be used as an aseessment tool to predict the outcomes and it needs further studies in the future.

Adaptive phase correction of diffusion-weighted images

Phase correction (PC) is a preprocessing technique that exploits the phase of images acquired in Magnetic Resonance Imaging (MRI) to obtain real-valued images containing tissue contrast with additive Gaussian noise, as opposed to magnitude images which follow a non-Gaussian distribution, e.g. Rician. PC finds its natural application to diffusion-weighted images (DWIs) due to their inherent low signal-to-noise ratio and consequent non-Gaussianity that induces a signal overestimation bias that propagates to the calculated diffusion indices. PC effectiveness depends upon the quality of the phase estimation, which is often performed via a regularization procedure. We show that a suboptimal regularization can produce alterations of the true image contrast in the real-valued phase-corrected images. We propose adaptive phase correction (APC), a method where the phase is estimated by using MRI noise information to perform a complex-valued image regularization that accounts for the local variance of the noise. We show, on synthetic and acquired data, that APC leads to phase-corrected real-valued DWIs that present a reduced number of alterations and a reduced bias. The substantial absence of parameters for which human input is required favors a straightforward integration of APC in MRI processing pipelines.

Hippocampal segmentation for brains with extensive atrophy using three-dimensional convolutional neural networks

Hippocampal volumetry is a critical biomarker of aging and dementia, and it is widely used as a predictor of cognitive performance; however, automated hippocampal segmentation methods are limited because the algorithms are (a) not publicly available, (b) subject to error with significant brain atrophy, cerebrovascular disease and lesions, and/or (c) computationally expensive or require parameter tuning. In this study, we trained a 3D convolutional neural network using 259 bilateral manually delineated segmentations collected from three studies, acquired at multiple sites on different scanners with variable protocols. Our training dataset consisted of elderly cases difficult to segment due to extensive atrophy, vascular disease, and lesions. Our algorithm, (HippMapp3r), was validated against four other publicly available state-of-the-art techniques (HippoDeep, FreeSurfer, SBHV, volBrain, and FIRST). HippMapp3r outperformed the other techniques on all three metrics, generating an average Dice of 0.89 and a correlation coefficient of 0.95. It was two orders of magnitude faster than some of the tested techniques. Further validation was performed on 200 subjects from two other disease populations (frontotemporal dementia and vascular cognitive impairment), highlighting our method's low outlier rate. We finally tested the methods on real and simulated "clinical adversarial" cases to study their robustness to corrupt, low-quality scans. The pipeline and models are available at: https://hippmapp3r.readthedocs.ioto facilitate the study of the hippocampus in large multisite studies.

Traumatic Parafalcine Subdural Hematoma: A Clinically Benign Finding

BACKGROUND

Guidelines for management of intracranial hemorrhage do not account for bleed location. We hypothesize that parafalcine subdural hematoma (SDH), as compared to convexity SDH, is a distinct clinical entity and these patients do not benefit from critical care monitoring or repeat imaging.

METHODS

We identified patients presenting to a single level I trauma center with isolated head injuries from February 2016 to August 2017. We identified 88 patients with isolated blunt traumatic parafalcine SDH and 228 with convexity SDH.

RESULTS

Demographics, comorbidities, and use of antiplatelet and anticoagulant agents were similar between the groups. As compared to patients with convexity SDH, patients with parafalcine SDH had a significantly lower incidence of radiographic progression, and had no cases of neurologic deterioration, neurosurgical intervention, or mortality (all P < 0.005). Compared to patients admitted to the intensive care unit, patients with parafalcine SDH admitted to the floor had a shorter length of stay (2.0 ± 1.6 versus 3.8 ± 2.9 d, P < 0.005) with no difference in outcomes.

CONCLUSIONS

Patients presenting with a parafalcine SDH are a distinct and relatively benign clinical entity as compared to convexity SDH and do not benefit from repeat imaging or intensive care unit admission.

Using Machine Learning and Structural Neuroimaging to Detect First Episode Psychosis: Reconsidering the Evidence

Despite the high level of interest in the use of machine learning (ML) and neuroimaging to detect psychosis at the individual level, the reliability of the findings is unclear due to potential methodological issues that may have inflated the existing literature. This study aimed to elucidate the extent to which the application of ML to neuroanatomical data allows detection of first episode psychosis (FEP), while putting in place methodological precautions to avoid overoptimistic results. We tested both traditional ML and an emerging approach known as deep learning (DL) using 3 feature sets of interest: (1) surface-based regional volumes and cortical thickness, (2) voxel-based gray matter volume (GMV) and (3) voxel-based cortical thickness (VBCT). To assess the reliability of the findings, we repeated all analyses in 5 independent datasets, totaling 956 participants (514 FEP and 444 within-site matched controls). The performance was assessed via nested cross-validation (CV) and cross-site CV. Accuracies ranged from 50% to 70% for surfaced-based features; from 50% to 63% for GMV; and from 51% to 68% for VBCT. The best accuracies (70%) were achieved when DL was applied to surface-based features; however, these models generalized poorly to other sites. Findings from this study suggest that, when methodological precautions are adopted to avoid overoptimistic results, detection of individuals in the early stages of psychosis is more challenging than originally thought. In light of this, we argue that the current evidence for the diagnostic value of ML and structural neuroimaging should be reconsidered toward a more cautious interpretation.

Assessing the Reliability of Reported Medical History in Older Adults

BACKGROUND

Age-associated increases in medical complexity, frailty, and cognitive impairment may compromise reliable reporting of medical history.

OBJECTIVE

To evaluate the influence of increasing age and cognitive impairment on concordance between reported history of stroke and cerebral infarction, and reported history of diabetes and elevated hemoglobinA1c in community-dwelling older adults.

METHODS

The association between participant-specific factors and accurate reporting of stroke or diabetes was evaluated using multivariable logistic regression in 1,401 participants enrolled in longitudinal studies of memory and aging, including 425 participants with dementia (30.3%). Stroke and diabetes were selected as index variables as gold standard measures of both were obtained in all participants: magnetic resonance neuroimaging for cerebral infarcts and hemoglobinA1c (≥6.5%) for diabetes.

RESULTS

Concordance between reported history of stroke and imaging-confirmed cerebral infarction was low (sensitivity: 17.4%, 8/46; specificity: 97.9%, 799/816). Small infarcts were strongly associated with inaccurate reporting (OR = 265.8; 95% CI: 86.2, 819.4), suggesting that occult/silent infarcts contributed to discordant reporting. Reporting accuracy was higher concerning diabetes (sensitivity: 83.5%, 147/176; specificity: 96.2%, 1100/1143). A history of hypertension (OR = 2.3; 95% CI: 1.3, 4.2), higher hemoglobinA1c (OR = 1.9; 95% CI: 1.5, 2.4), and hemoglobinA1c compatible with impaired glucose tolerance (OR = 3.1; 95% CI 1.8, 5.3) associated with increased odds of discordant reporting. Cognitive impairment and increased age were not independently associated with reliable reporting.

CONCLUSION

Factors beyond advancing age and cognitive impairment appear to drive discordance in reported medical history in older participants. Objective testing for cerebral infarcts or diabetes should be performed when relevant to diagnostic or therapeutic decisions in clinical and research settings.

Magnetic resonance imaging texture predicts progression to dementia due to Alzheimer disease earlier than hippocampal volume

BACKGROUND

Early identification of people at risk of imminent progression to dementia due to Alzheimer disease is crucial for timely intervention and treatment. We investigated whether the texture of MRI brain scans could predict the progression of mild cognitive impairment (MCI) to Alzheimer disease earlier than volume.

METHODS

We constructed a development data set (121 people who were cognitively normal and 145 who had mild Alzheimer disease) and a validation data set (113 patients with stable MCI who did not progress to Alzheimer disease for 3 years; 40 with early MCI who progressed to Alzheimer disease after 12–36 months; and 41 with late MCI who progressed to Alzheimer disease within 12 months) from the Alzheimer’s Disease Neuroimaging Initiative. We analyzed the texture of the hippocampus, precuneus and posterior cingulate cortex using a grey-level co-occurrence matrix. We constructed texture and volume indices from the development data set using logistic regression. Using area under the curve (AUC) of receiver operator characteristics, we compared the accuracy of hippocampal volume, hippocampal texture and the composite texture of the hippocampus, precuneus and posterior cingulate cortex in predicting conversion from MCI to Alzheimer disease in the validation data set.

RESULTS

Compared with hippocampal volume, hippocampal texture (0.790 v. 0.739, p = 0.047) and composite texture (0.811 v. 0.739, p = 0.007) showed larger AUCs for conversion to Alzheimer disease from both early and late MCI. Hippocampal texture showed a marginally larger AUC than hippocampal volume in early MCI (0.795 v. 0.726, p = 0.060). Composite texture showed a larger AUC for conversion to Alzheimer disease than hippocampal volume in both early (0.817 v. 0.726, p = 0.027) and late MCI (0.805 v. 0.753, p = 0.019).

LIMITATIONS

This study was limited by the absence of histological data, and the pathology reflected by the texture measures remains to be validated.

CONCLUSION

Textures of the hippocampus, precuneus and posterior cingulate cortex predicted conversion from MCI to Alzheimer disease at an earlier time point and with higher accuracy than hippocampal volume.

Unification of optimal targeting methods in transcranial electrical stimulation

One of the major questions in high-density transcranial electrical stimulation (TES) is: given a region of interest (ROI) and electric current limits for safety, how much current should be delivered by each electrode for optimal targeting of the ROI? Several solutions, apparently unrelated, have been independently proposed depending on how "optimality" is defined and on how this optimization problem is stated mathematically. The least squares (LS), weighted LS (WLS), or reciprocity-based approaches are the simplest ones and have closed-form solutions. An extended optimization problem can be stated as follows: maximize the directional intensity at the ROI, limit the electric fields at the non-ROI, and constrain total injected current and current per electrode for safety. This problem requires iterative convex or linear optimization solvers. We theoretically prove in this work that the LS, WLS and reciprocity-based closed-form solutions are specific solutions to the extended directional maximization optimization problem. Moreover, the LS/WLS and reciprocity-based solutions are the two extreme cases of the intensity-focality trade-off, emerging under variation of a unique parameter of the extended directional maximization problem, the imposed constraint to the electric fields at the non-ROI. We validate and illustrate these findings with simulations on an atlas head model. The unified approach we present here allows a better understanding of the nature of the TES optimization problem and helps in the development of advanced and more effective targeting strategies.