Search

Abstract Associations between datasets can be discovered through multivariate methods like Canonical Correlation Analysis (CCA) or Partial Least Squares (PLS). A requisite property for interpretability and generalizability of CCA/PLS associations is stability of their feature patterns. However, stability of CCA/PLS in high-dimensional datasets is questionable, as found in empirical characterizations. To study these issues systematically, we developed a generative modeling framework to simulate synthetic datasets. We found that when sample size is relatively small, but comparable to typical studies, CCA/PLS associations are highly unstable and inaccurate; both in their magnitude and importantly in the feature pattern underlying the association. We confirmed these trends across two neuroimaging modalities and in independent datasets with n ≈ 1000 and n = 20,000, and found that only the latter comprised sufficient observations for stable mappings between imaging-derived and behavioral features. We further developed a power calculator to provide sample sizes required for stability and reliability of multivariate analyses. Collectively, we characterize how to limit detrimental effects of overfitting on CCA/PLS stability, and provide recommendations for future studies.

IntroductionCerebral microbleeds (CMBs) are associated with white matter damage, and various neurodegenerative and cerebrovascular diseases. CMBs occur as small, circular hypointense lesions on T2*-weighted gradient recalled echo (GRE) and susceptibility-weighted imaging (SWI) images, and hyperintense on quantitative susceptibility mapping (QSM) images due to their paramagnetic nature. Accurate automated detection of CMBs would help to determine quantitative imaging biomarkers (e.g., CMB count) on large datasets. In this work, we propose a fully automated, deep learning-based, 3-step algorithm, using structural and anatomical properties of CMBs from any single input image modality (e.g., GRE/SWI/QSM) for their accurate detections.MethodsIn our method, the first step consists of an initial candidate detection step that detects CMBs with high sensitivity. In the second step, candidate discrimination step is performed using a knowledge distillation framework, with a multi-tasking teacher network that guides the student network to classify CMB and non-CMB instances in an offline manner. Finally, a morphological clean-up step further reduces false positives using anatomical constraints. We used four datasets consisting of different modalities specified above, acquired using various protocols and with a variety of pathological and demographic characteristics.ResultsOn cross-validation within datasets, our method achieved a cluster-wise true positive rate (TPR) of over 90% with an average of 80 % with different modalities. The python implementation of the proposed method is openly available.

IntroductionNeuroimaging technology has experienced explosive growth and transformed the study of neural mechanisms across health and disease. However, given the diversity of sophisticated tools for handling neuroimaging data, the field faces challenges in method integration, particularly across multiple modalities and species. Specifically, researchers often have to rely on siloed approaches which limit reproducibility, with idiosyncratic data organization and limited software interoperability.MethodsTo address these challenges, we have developed Quantitative Neuroimaging Environment & Toolbox (QuNex), a platform for consistent end-to-end processing and analytics. QuNex provides several novel functionalities for neuroimaging analyses, including a “turnkey” command for the reproducible deployment of custom workflows, from onboarding raw data to generating analytic features.ResultsThe platform enables interoperable integration of multi-modal, community-developed neuroimaging software through an extension framework with a software development kit (SDK) for seamless integration of community tools. Critically, it supports high-throughput, parallel processing in high-performance compute environments, either locally or in the cloud. Notably, QuNex has successfully processed over 10,000 scans across neuroimaging consortia, including multiple clinical datasets. Moreover, QuNex enables integration of human and non-human workflows via a cohesive translational platform.DiscussionCollectively, this effort stands to significantly impact neuroimaging method integration across acquisition approaches, pipelines, datasets, computational environments, and species. Building on this platform will enable more rapid, scalable, and reproducible impact of neuroimaging technology across health and disease.

Summary: Background: Idiopathic pulmonary fibrosis (IPF) is a serious disease of the lung parenchyma. It has a known polygenetic risk, with at least seventeen regions of the genome implicated to date. Growing evidence suggests linked multimorbidity of IPF with neurodegenerative or affective disorders. However, no study so far has explicitly explored links between IPF, associated genetic risk profiles, and specific brain features. Methods: We exploited imaging and genetic data from more than 32,000 participants available through the UK Biobank population-level resource to explore links between IPF genetic risk and imaging-derived brain endophenotypes. We performed a brain-wide imaging-genetics association study between the presence of 17 known IPF risk variants and 1248 multi-modal imaging-derived features, which characterise brain structure and function. Findings: We identified strong associations between cortical morphological features, white matter microstructure and IPF risk loci in chromosomes 17 (17q21.31) and 8 (DEPTOR). Through co-localisation analysis, we confirmed that cortical thickness in the anterior cingulate and more widespread white matter microstructure changes share a single causal variant with IPF at the chromosome 8 locus. Post-hoc preliminary analysis suggested that forced vital capacity may partially mediate the association between the DEPTOR variant and white matter microstructure, but not between the DEPTOR risk variant and cortical thickness. Interpretation: Our results reveal the associations between IPF genetic risk and differences in brain structure, for both cortex and white matter. Differences in tissue-specific imaging signatures suggest distinct underlying mechanisms with focal cortical thinning in regions with known high DEPTOR expression, unrelated to lung function, and more widespread microstructural white matter changes consistent with hypoxia or neuroinflammation with potential mediation by lung function. Funding: This study was supported by the NIHR Nottingham Biomedical Research Centre and the UK Medical Research Council.

Cerebral microbleeds (CMBs) appear as small, circular, well defined hypointense lesions of a few mm in size on T2*-weighted gradient recalled echo (T2*-GRE) images and appear enhanced on susceptibility weighted images (SWI). Due to their small size, contrast variations and other mimics (e.g., blood vessels), CMBs are highly challenging to detect automatically. In large datasets (e.g., the UK Biobank dataset), exhaustively labelling CMBs manually is difficult and time consuming. Hence it would be useful to preselect candidate CMB subjects in order to focus on those for manual labelling, which is essential for training and testing automated CMB detection tools on these datasets. In this work, we aim to detect CMB candidate subjects from a larger dataset, UK Biobank, using a machine learning-based, computationally light pipeline. For our evaluation, we used 3 different datasets, with different intensity characteristics, acquired with different scanners. They include the UK Biobank dataset and two clinical datasets with different pathological conditions. We developed and evaluated our pipelines on different types of images, consisting of SWI or GRE images. We also used the UK Biobank dataset to compare our approach with alternative CMB preselection methods using non-imaging factors and/or imaging data. Finally, we evaluated the pipeline's generalisability across datasets. Our method provided subject-level detection accuracy > 80% on all the datasets (within-dataset results), and showed good generalisability across datasets, providing a consistent accuracy of over 80%, even when evaluated across different modalities.

The Human Connectome Project (HCP) was launched in 2010 as an ambitious effort to accelerate advances in human neuroimaging, particularly for measures of brain connectivity; apply these advances to study a large number of healthy young adults; and freely share the data and tools with the scientific community. NIH awarded grants to two consortia; this retrospective focuses on the “WU-Minn-Ox” HCP consortium centered at Washington University, the University of Minnesota, and University of Oxford. In just over 6 years, the WU-Minn-Ox consortium succeeded in its core objectives by: 1) improving MR scanner hardware, pulse sequence design, and image reconstruction methods, 2) acquiring and analyzing multimodal MRI and MEG data of unprecedented quality together with behavioral measures from more than 1100 HCP participants, and 3) freely sharing the data (via the ConnectomeDB database) and associated analysis and visualization tools. To date, more than 27 Petabytes of data have been shared, and 1538 papers acknowledging HCP data use have been published. The “HCP-style” neuroimaging paradigm has emerged as a set of best-practice strategies for optimizing data acquisition and analysis. This article reviews the history of the HCP, including comments on key events and decisions associated with major project components. We discuss several scientific advances using HCP data, including improved cortical parcellations, analyses of connectivity based on functional and diffusion MRI, and analyses of brain-behavior relationships. We also touch upon our efforts to develop and share a variety of associated data processing and analysis tools along with detailed documentation, tutorials, and an educational course to train the next generation of neuroimagers. We conclude with a look forward at opportunities and challenges facing the human neuroimaging field from the perspective of the HCP consortium.

Many brain imaging studies aim to measure structural connectivity with diffusion tractography. However, biases in tractography data, particularly near the boundary between white matter and cortical grey matter can limit the accuracy of such studies. When seeding from the white matter, streamlines tend to travel parallel to the convoluted cortical surface, largely avoiding sulcal fundi and terminating preferentially on gyral crowns. When seeding from the cortical grey matter, streamlines generally run near the cortical surface until reaching deep white matter. These so-called “gyral biases” limit the accuracy and effective resolution of cortical structural connectivity profiles estimated by tractography algorithms, and they do not reflect the expected distributions of axonal densities seen in invasive tracer studies or stains of myelinated fibres. We propose an algorithm that concurrently models fibre density and orientation using a divergence-free vector field within gyral blades to encourage an anatomically-justified streamline density distribution along the cortical white/grey-matter boundary while maintaining alignment with the diffusion MRI estimated fibre orientations. Using in vivo data from the Human Connectome Project, we show that this algorithm reduces tractography biases. We compare the structural connectomes to functional connectomes from resting-state fMRI, showing that our model improves cross-modal agreement. Finally, we find that after parcellation the changes in the structural connectome are very minor with slightly improved interhemispheric connections (i.e, more homotopic connectivity) and slightly worse intrahemispheric connections when compared to tracers.

Deep learning (DL) has shown great potential in medical image enhancement problems, such as super-resolution or image synthesis. However, to date, most existing approaches are based on deterministic models, neglecting the presence of different sources of uncertainty in such problems. Here we introduce methods to characterise different components of uncertainty, and demonstrate the ideas using diffusion MRI super-resolution. Specifically, we propose to account for intrinsic uncertainty through a heteroscedastic noise model and for parameter uncertainty through approximate Bayesian inference, and integrate the two to quantify predictive uncertainty over the output image. Moreover, we introduce a method to propagate the predictive uncertainty on a multi-channelled image to derived scalar parameters, and separately quantify the effects of intrinsic and parameter uncertainty therein. The methods are evaluated for super-resolution of two different signal representations of diffusion MR images—Diffusion Tensor images and Mean Apparent Propagator MRI—and their derived quantities such as mean diffusivity and fractional anisotropy, on multiple datasets of both healthy and pathological human brains. Results highlight three key potential benefits of modelling uncertainty for improving the safety of DL-based image enhancement systems. Firstly, modelling uncertainty improves the predictive performance even when test data departs from training data (“out-of-distribution” datasets). Secondly, the predictive uncertainty highly correlates with reconstruction errors, and is therefore capable of detecting predictive “failures”. Results on both healthy subjects and patients with brain glioma or multiple sclerosis demonstrate that such an uncertainty measure enables subject-specific and voxel-wise risk assessment of the super-resolved images that can be accounted for in subsequent analysis. Thirdly, we show that the method for decomposing predictive uncertainty into its independent sources provides high-level “explanations” for the model performance by separately quantifying how much uncertainty arises from the inherent difficulty of the task or the limited training examples. The introduced concepts of uncertainty modelling extend naturally to many other imaging modalities and data enhancement applications.

We present a new software package with a library of standardised tractography protocols devised for the robust automated extraction of white matter tracts both in the human and the macaque brain. Using in vivo data from the Human Connectome Project (HCP) and the UK Biobank and ex vivo data for the macaque brain datasets, we obtain white matter atlases, as well as atlases for tract endpoints on the white-grey matter boundary, for both species. We illustrate that our protocols are robust against data quality, generalisable across two species and reflect the known anatomy. We further demonstrate that they capture inter-subject variability by preserving tract lateralisation in humans and tract similarities stemming from twinship in the HCP cohort. Our results demonstrate that the presented toolbox will be useful for generating imaging-derived features in large cohorts, and in facilitating comparative neuroanatomy studies. The software, tractography protocols, and atlases are publicly released through FSL, allowing users to define their own tractography protocols in a standardised manner, further contributing to open science.

Measuring fibre dispersion in white matter with diffusion magnetic resonance imaging (MRI) is limited by an inherent degeneracy between fibre dispersion and microscopic diffusion anisotropy (i.e., the diffusion anisotropy expected for a single fibre orientation). This means that estimates of fibre dispersion rely on strong assumptions, such as constant microscopic anisotropy throughout the white matter or specific biophysical models. Here we present a simple approach for resolving this degeneracy using measurements that combine linear (conventional) and spherical tensor diffusion encoding.To test the accuracy of the fibre dispersion when our microstructural model is only an approximation of the true tissue structure, we simulate multi-compartment data and fit this with a single-compartment model. For such overly simplistic tissue assumptions, we show that the bias in fibre dispersion is greatly reduced (~5x) for single-shell linear and spherical tensor encoding data compared with single-shell or multi-shell conventional data. In in-vivo data we find a consistent estimate of fibre dispersion as we reduce the b-value from 3 to 1.5 ms/μm2, increase the repetition time, increase the echo time, or increase the diffusion time. We conclude that the addition of spherical tensor encoded data to conventional linear tensor encoding data greatly reduces the sensitivity of the estimated fibre dispersion to the model assumptions of the tissue microstructure.

Macaque monkeys are an important animal model where invasive investigations can lead to a better understanding of the cortical organization of primates including humans. However, the tools and methods for noninvasive image acquisition (e.g. MRI RF coils and pulse sequence protocols) and image data preprocessing have lagged behind those developed for humans. To resolve the structural and functional characteristics of the smaller macaque brain, high spatial, temporal, and angular resolutions combined with high signal-to-noise ratio are required to ensure good image quality. To address these challenges, we developed a macaque 24-channel receive coil for 3-T MRI with parallel imaging capabilities. This coil enables adaptation of the Human Connectome Project (HCP) image acquisition protocols to the in-vivo macaque brain. In addition, we adapted HCP preprocessing methods to the macaque brain, including spatial minimal preprocessing of structural, functional MRI (fMRI), and diffusion MRI (dMRI). The coil provides the necessary high signal-to-noise ratio and high efficiency in data acquisition, allowing four- and five-fold accelerations for dMRI and fMRI. Automated FreeSurfer segmentation of cortex, reconstruction of cortical surface, removal of artefacts and nuisance signals in fMRI, and distortion correction of dMRI all performed well, and the overall quality of basic neurobiological measures was comparable with those for the HCP. Analyses of functional connectivity in fMRI revealed high sensitivity as compared with those from publicly shared datasets. Tractography-based connectivity estimates correlated with tracer connectivity similarly to that achieved using ex-vivo dMRI. The resulting HCP-style in vivo macaque MRI data show considerable promise for analyzing cortical architecture and functional and structural connectivity using advanced methods that have previously only been available in studies of the human brain.

Diffusion MRI is an exquisitely sensitive probe of tissue microstructure, and is currently the only non-invasive measure of the brain's fibre architecture. As this technique becomes more sophisticated and microstructurally informative, there is increasing value in comparing diffusion MRI with microscopic imaging in the same tissue samples. This study compared estimates of fibre orientation dispersion in white matter derived from diffusion MRI to reference measures of dispersion obtained from polarized light imaging and histology.Three post-mortem brain specimens were scanned with diffusion MRI and analyzed with a two-compartment dispersion model. The specimens were then sectioned for microscopy, including polarized light imaging estimates of fibre orientation and histological quantitative estimates of myelin and astrocytes. Dispersion estimates were correlated on region – and voxel-wise levels in the corpus callosum, the centrum semiovale and the corticospinal tract.The region-wise analysis yielded correlation coefficients of r = 0.79 for the diffusion MRI and histology comparison, while r = 0.60 was reported for the comparison with polarized light imaging. In the corpus callosum, we observed a pattern of higher dispersion at the midline compared to its lateral aspects. This pattern was present in all modalities and the dispersion profiles from microscopy and diffusion MRI were highly correlated. The astrocytes appeared to have minor contribution to dispersion observed with diffusion MRI.These results demonstrate that fibre orientation dispersion estimates from diffusion MRI represents the tissue architecture well. Dispersion models might be improved by more faithfully incorporating an informed mapping based on microscopy data.