A statistical method for image-mediated association studies discovers genes and pathways associated with four brain disorders.

Brain imaging and genomics are critical tools enabling characterization of the genetic basis of brain disorders. However, imaging large cohorts is expensive and may be unavailable for legacy datasets used for genome-wide association studies (GWASs). Using an integrated feature selection/aggregation model, we developed an image-mediated association study (IMAS), which utilizes borrowed imaging/genomics data to conduct association mapping in legacy GWAS cohorts. By leveraging the UK Biobank image-derived phenotypes (IDPs), the IMAS discovered genetic bases underlying four neuropsychiatric disorders and verified them by analyzing annotations, pathways, and expression quantitative trait loci (eQTLs). A cerebellar-mediated mechanism was identified to be common to the four disorders. Simulations show that, if the goal is identifying genetic risk, our IMAS is more powerful than a hypothetical protocol in which the imaging results were available in the GWAS dataset. This implies the feasibility of reanalyzing legacy GWAS datasets without conducting additional imaging, yielding cost savings for integrated analysis of genetics and imaging.

Image Translation for Estimating Two-Dimensional Axial Amyloid-Beta PET From Structural MRI.

BACKGROUND: Amyloid-beta and brain atrophy are hallmarks for Alzheimer's Disease that can be targeted with positron emission tomography (PET) and MRI, respectively. MRI is cheaper, less-invasive, and more available than PET. There is a known relationship between amyloid-beta and brain atrophy, meaning PET images could be inferred from MRI. PURPOSE: To build an image translation model using a Conditional Generative Adversarial Network able to synthesize Amyloid-beta PET images from structural MRI. STUDY TYPE: Retrospective. POPULATION: Eight hundred eighty-two adults (348 males/534 females) with different stages of cognitive decline (control, mild cognitive impairment, moderate cognitive impairment, and severe cognitive impairment). Five hundred fifty-two subjects for model training and 331 for testing (80%:20%). FIELD STRENGTH/SEQUENCE: 3 T, T1-weighted structural (T1w). ASSESSMENT: The testing cohort was used to evaluate the performance of the model using the Structural Similarity Index Measure (SSIM) and Peak Signal-to-Noise Ratio (PSNR), comparing the likeness of the overall synthetic PET images created from structural MRI with the overall true PET images. SSIM was computed in the overall image to include the luminance, contrast, and structural similarity components. Experienced observers reviewed the images for quality, performance and tried to determine if they could tell the difference between real and synthetic images. STATISTICAL TESTS: Pixel wise Pearson correlation was significant, and had an R2 greater than 0.96 in example images. From blinded readings, a Pearson Chi-squared test showed that there was no significant difference between the real and synthetic images by the observers (P = 0.68). RESULTS: A high degree of likeness across the evaluation set, which had a mean SSIM = 0.905 and PSNR = 2.685. The two observers were not able to determine the difference between the real and synthetic images, with accuracies of 54% and 46%, respectively. CONCLUSION: Amyloid-beta PET images can be synthesized from structural MRI with a high degree of similarity to the real PET images. EVIDENCE LEVEL: 3 TECHNICAL EFFICACY: Stage 1.

Direct machine learning reconstruction of respiratory variation waveforms from resting state fMRI data in a pediatric population.

In many functional magnetic resonance imaging (fMRI) studies, respiratory signals are unavailable or do not have acceptable quality due to issues with subject compliance, equipment failure or signal error. In large databases, such as the Human Connectome Projects, over half of the respiratory recordings may be unusable. As a result, the direct removal of low frequency respiratory variations from the blood oxygen level-dependent (BOLD) signal time series is not possible. This study proposes a deep learning-based method for reconstruction of respiratory variation (RV) waveforms directly from BOLD fMRI data in pediatric participants (aged 5 to 21 years old), and does not require any respiratory measurement device. To do this, the Lifespan Human Connectome Project in Development (HCP-D) dataset, which includes respiratory measurements, was used to both train a convolutional neural network (CNN) and evaluate its performance. Results show that a CNN can capture informative features from the BOLD signal time course and reconstruct accurate RV timeseries, especially when the subject has a prominent respiratory event. This work advances the use of direct estimation of physiological parameters from fMRI, which will eventually lead to reduced complexity and decrease the burden on participants because they may not be required to wear a respiratory bellows.

MRI of healthy brain aging: A review.

We present a review of the characterization of healthy brain aging using MRI with an emphasis on morphology, lesions, and quantitative MR parameters. A scope review found 6612 articles encompassing the keywords "Brain Aging" and "Magnetic Resonance"; papers involving functional MRI or not involving imaging of healthy human brain aging were discarded, leaving 2246 articles. We first consider some of the biogerontological mechanisms of aging, and the consequences of aging in terms of cognition and onset of disease. Morphological changes with aging are reviewed for the whole brain, cerebral cortex, white matter, subcortical gray matter, and other individual structures. In general, volume and cortical thickness decline with age, beginning in mid-life. Prevalent silent lesions such as white matter hyperintensities, microbleeds, and lacunar infarcts are also observed with increasing frequency. The literature regarding quantitative MR parameter changes includes T1 , T2 , T2 *, magnetic susceptibility, spectroscopy, magnetization transfer, diffusion, and blood flow. We summarize the findings on how each of these parameters varies with aging. Finally, we examine how the aforementioned techniques have been used for age prediction. While relatively large in scope, we present a comprehensive review that should provide the reader with sound understanding of what MRI has been able to tell us about how the healthy brain ages.

Interdatabase Variability in Cortical Thickness Measurements.

The phenomenon of cortical thinning with age has been well established; however, the measured rate of change varies between studies. The source of this variation could be image acquisition techniques including hardware and vendor specific differences. Databases are often consolidated to increase the number of subjects but underlying differences between these datasets could have undesired effects. We explore differences in cerebral cortex thinning between 4 databases, totaling 1382 subjects. We investigate several aspects of these databases, including: 1) differences between databases of cortical thinning rates versus age, 2) correlation of cortical thinning rates between regions for each database, and 3) regression bootstrapping to determine the effect of the number of subjects included. We also examined the effect of different databases on age prediction modeling. Cortical thinning rates were significantly different between databases in all 68 parcellated regions (ANCOVA, P < 0.001). Subtle differences were observed in correlation matrices and bootstrapping convergence. Age prediction modeling using a leave-one-out cross-validation approach showed varying prediction performance (0.64 < R2 < 0.82) between databases. When a database was used to calibrate the model and then applied to another database, prediction performance consistently decreased. We conclude that there are indeed differences in the measured cortical thinning rates between these large-scale databases.

Whole head quantitative susceptibility mapping using a least-norm direct dipole inversion method.

A new dipole field inversion method for whole head quantitative susceptibility mapping (QSM) is proposed. Instead of performing background field removal and local field inversion sequentially, the proposed method performs dipole field inversion directly on the total field map in a single step. To aid this under-determined and ill-posed inversion process and obtain robust QSM images, Tikhonov regularization is implemented to seek the local susceptibility solution with the least-norm (LN) using the L-curve criterion. The proposed LN-QSM does not require brain edge erosion, thereby preserving the cerebral cortex in the final images. This should improve its applicability for QSM-based cortical grey matter measurement, functional imaging and venography of full brain. Furthermore, LN-QSM also enables susceptibility mapping of the entire head without the need for brain extraction, which makes QSM reconstruction more automated and less dependent on intermediate pre-processing methods and their associated parameters. It is shown that the proposed LN-QSM method reduced errors in a numerical phantom simulation, improved accuracy in a gadolinium phantom experiment, and suppressed artefacts in nine subjects, as compared to two-step and other single-step QSM methods. Measurements of deep grey matter and skull susceptibilities from LN-QSM are consistent with established reconstruction methods.

Gas-free calibrated fMRI with a correction for vessel-size sensitivity.

Calibrated functional magnetic resonance imaging (fMRI) is a method to independently measure the metabolic and hemodynamic contributions to the blood oxygenation level dependent (BOLD) signal. This technique typically requires the use of a respiratory challenge, such as hypercapnia or hyperoxia, to estimate the calibration constant, M. There has been a recent push to eliminate the gas challenge from the calibration procedure using asymmetric spin echo (ASE) based techniques. This study uses simulations to better understand spin echo (SE) and ASE signals, analytical modelling to characterize the signal evolution, and in vivo imaging to validate the modelling. Using simulations, it is shown how ASE imaging generally underestimates M and how this depends on several parameters of the acquisition, including echo time and ASE offset, as well as the vessel size. This underestimation is the result of imperfect SE refocusing due to diffusion of water through the extravascular environment surrounding the microvasculature. By empirically characterizing this SE attenuation as an exponential decay that increases with echo time, we have proposed a quadratic ASE biophysical signal model. This model allows for the characterization and compensation of the SE attenuation if SE and ASE signals are acquired at multiple echo times. This was tested in healthy subjects and was found to significantly increase the estimates of M across grey matter. These findings show promise for improved gas-free calibration and can be extended to other relaxation-based imaging studies of brain physiology.

Modeling hyperoxia-induced BOLD signal dynamics to estimate cerebral blood flow, volume and mean transit time.

A new method is proposed for obtaining cerebral perfusion measurements whereby blood oxygen level dependent (BOLD) MRI is used to dynamically monitor hyperoxia-induced changes in the concentration of deoxygenated hemoglobin in the cerebral vasculature. The data is processed using kinetic modeling to yield perfusion metrics, namely: cerebral blood flow (CBF), cerebral blood volume (CBV), and mean transit time (MTT). Ten healthy human subjects were continuously imaged with BOLD sequence while a hyperoxic (70% O2) state was interspersed with baseline periods of normoxia. The BOLD time courses were fit with exponential uptake and decay curves and a biophysical model of the BOLD signal was used to estimate oxygen concentration functions. The arterial input function was derived from end-tidal oxygen measurements, and a deconvolution operation between the tissue and arterial concentration functions was used to yield CBF. The venous component of the CBV was calculated from the ratio of the integrals of the estimated tissue and arterial concentration functions. Mean gray and white matter measurements were found to be: 61.6 ±â€¯13.7 and 24.9 ±â€¯4.0 ml 100 g-1 min-1 for CBF; 1.83 ±â€¯0.32 and 1.10 ±â€¯0.19 ml 100 g-1 for venous CBV; and 2.94 ±â€¯0.52 and 3.73 ±â€¯0.60 s for MTT, respectively. We conclude that it is possible to derive CBF, CBV and MTT metrics within expected physiological ranges via analysis of dynamic BOLD fMRI acquired during a period of hyperoxia.

Closed-loop vasculature network design for bioprinting large, solid tissue scaffolds.

Vascularization is an indispensable requirement for fabricating large solid tissues and organs. The natural vasculature derived from medical imaging modalities for large tissues and organs are highly complex and convoluted. However, the present bioprinting capabilities limit the fabrication of such complex natural vascular networks. Simplified bioprinted vascular networks, on the other hand, lack the capability to sustain large solid tissues. This work proposes a generalized and adaptable numerical model to design the vasculature by utilizing the tissue/organ anatomy. Starting with processing the patient's medical images, organ structure, tissue-specific cues, and key vasculature tethers are determined. An open-source abdomen magnetic resonance image dataset was used in this work. The extracted properties and cues are then used in a mathematical model for guiding the vascular network formation comprising arterial and venous networks. Next, the generated three-dimensional networks are used to simulate the nutrient transport and consumption within the organ over time and the regions deprived of the nutrients are identified. These regions provide cues to evolve and optimize the vasculature in an iterative manner to ensure the availability of the nutrient transport throughout the bioprinted scaffolds. The mass transport of six components of cell culture media-glucose, glycine, glutamine, riboflavin, human serum albumin, and oxygen was studied within the organ with designed vasculature. As the vascular structure underwent iterations, the organ regions deprived of these key components decreased significantly highlighting the increase in structural complexity and efficacy of the designed vasculature. The numerical method presented in this work offers a valuable tool for designing vascular scaffolds to guide the cell growth and maturation of the bioprinted tissues for faster regeneration post bioprinting.

Correspondence between BOLD fMRI task response and cerebrovascular reactivity across the cerebral cortex.

BOLD sensitivity to baseline perfusion and blood volume is a well-acknowledged fMRI confound. Vascular correction techniques based on cerebrovascular reactivity (CVR) might reduce variance due to baseline cerebral blood volume, however this is predicated on an invariant linear relationship between CVR and BOLD signal magnitude. Cognitive paradigms have relatively low signal, high variance and involve spatially heterogenous cortical regions; it is therefore unclear whether the BOLD response magnitude to complex paradigms can be predicted by CVR. The feasibility of predicting BOLD signal magnitude from CVR was explored in the present work across two experiments using different CVR approaches. The first utilized a large database containing breath-hold BOLD responses and 3 different cognitive tasks. The second experiment, in an independent sample, calculated CVR using the delivery of a fixed concentration of carbon dioxide and a different cognitive task. An atlas-based regression approach was implemented for both experiments to evaluate the shared variance between task-invoked BOLD responses and CVR across the cerebral cortex. Both experiments found significant relationships between CVR and task-based BOLD magnitude, with activation in the right cuneus (R 2 = 0.64) and paracentral gyrus (R 2 = 0.71), and the left pars opercularis (R 2 = 0.67), superior frontal gyrus (R 2 = 0.62) and inferior parietal cortex (R 2 = 0.63) strongly predicted by CVR. The parietal regions bilaterally were highly consistent, with linear regressions significant in these regions for all four tasks. Group analyses showed that CVR correction increased BOLD sensitivity. Overall, this work suggests that BOLD signal response magnitudes to cognitive tasks are predicted by CVR across different regions of the cerebral cortex, providing support for the use of correction based on baseline vascular physiology.

Lesion-symptom mapping with NIHSS sub-scores in ischemic stroke patients.

BACKGROUND: Lesion-symptom mapping (LSM) is a statistical technique to investigate the population-specific relationship between structural integrity and post-stroke clinical outcome. In clinical practice, patients are commonly evaluated using the National Institutes of Health Stroke Scale (NIHSS), an 11-domain clinical score to quantitate neurological deficits due to stroke. So far, LSM studies have mostly used the total NIHSS score for analysis, which might not uncover subtle structure-function relationships associated with the specific sub-domains of the NIHSS evaluation. Thus, the aim of this work was to investigate the feasibility to perform LSM analyses with sub-score information to reveal category-specific structure-function relationships that a total score may not reveal. METHODS: Employing a multivariate technique, LSM analyses were conducted using a sample of 180 patients with NIHSS assessment at 48-hour post-stroke from the ESCAPE trial. The NIHSS domains were grouped into six categories using two schemes. LSM was conducted for each category of the two groupings and the total NIHSS score. RESULTS: Sub-score LSMs not only identify most of the brain regions that are identified as critical by the total NIHSS score but also reveal additional brain regions critical to each function category of the NIHSS assessment without requiring extensive, specialised assessments. CONCLUSION: These findings show that widely available sub-scores of clinical outcome assessments can be used to investigate more specific structure-function relationships, which may improve predictive modelling of stroke outcomes in the context of modern clinical stroke assessments and neuroimaging. TRIAL REGISTRATION NUMBER: NCT01778335.

Invertible Modeling of Bidirectional Relationships in Neuroimaging With Normalizing Flows: Application to Brain Aging.

Many machine learning tasks in neuroimaging aim at modeling complex relationships between a brain's morphology as seen in structural MR images and clinical scores and variables of interest. A frequently modeled process is healthy brain aging for which many image-based brain age estimation or age-conditioned brain morphology template generation approaches exist. While age estimation is a regression task, template generation is related to generative modeling. Both tasks can be seen as inverse directions of the same relationship between brain morphology and age. However, this view is rarely exploited and most existing approaches train separate models for each direction. In this paper, we propose a novel bidirectional approach that unifies score regression and generative morphology modeling and we use it to build a bidirectional brain aging model. We achieve this by defining an invertible normalizing flow architecture that learns a probability distribution of 3D brain morphology conditioned on age. The use of full 3D brain data is achieved by deriving a manifold-constrained formulation that models morphology variations within a low-dimensional subspace of diffeomorphic transformations. This modeling idea is evaluated on a database of MR scans of more than 5000 subjects. The evaluation results show that our bidirectional brain aging model (1) accurately estimates brain age, (2) is able to visually explain its decisions through attribution maps and counterfactuals, (3) generates realistic age-specific brain morphology templates, (4) supports the analysis of morphological variations, and (5) can be utilized for subject-specific brain aging simulation.

High-resolution T2-FLAIR and non-contrast CT brain atlas of the elderly.

Normative brain atlases are a standard tool for neuroscience research and are, for example, used for spatial normalization of image datasets prior to voxel-based analyses of brain morphology and function. Although many different atlases are publicly available, they are usually biased with respect to an imaging modality and the age distribution. Both effects are well known to negatively impact the accuracy and reliability of the spatial normalization process using non-linear image registration methods. An important and very active neuroscience area that lacks appropriate atlases is lesion-related research in elderly populations (e.g. stroke, multiple sclerosis) for which FLAIR MRI and non-contrast CT are often the clinical imaging modalities of choice. To overcome the lack of atlases for these tasks and modalities, this paper presents high-resolution, age-specific FLAIR and non-contrast CT atlases of the elderly generated using clinical images.

Age-related differences in cerebral blood flow and cortical thickness with an application to age prediction.

Cerebral cortex thinning and cerebral blood flow (CBF) reduction are typically observed during normal healthy aging. However, imaging-based age prediction models have primarily used morphological features of the brain. Complementary physiological CBF information might result in an improvement in age estimation. In this study, T1-weighted structural magnetic resonance imaging and arterial spin labeling CBF images were acquired in 146 healthy participants across the adult life span. Sixty-eight cerebral cortex regions were segmented, and the cortical thickness and mean CBF were computed for each region. Linear regression with age was computed for each region and data type, and laterality and correlation matrices were computed. Sixteen predictive models were trained with the cortical thickness and CBF data alone as well as a combination of both data types. The age explained more variance in the cortical thickness data (average R2 of 0.21) than in the CBF data (average R2 of 0.09). All 16 models performed significantly better when combining both measurement types and using feature selection, and thus, we conclude that the inclusion of CBF data marginally improves age estimation.

Supervised machine learning tools: a tutorial for clinicians.

In an increasingly data-driven world, artificial intelligence is expected to be a key tool for converting big data into tangible benefits and the healthcare domain is no exception to this. Machine learning aims to identify complex patterns in multi-dimensional data and use these uncovered patterns to classify new unseen cases or make data-driven predictions. In recent years, deep neural networks have shown to be capable of producing results that considerably exceed those of conventional machine learning methods for various classification and regression tasks. In this paper, we provide an accessible tutorial of the most important supervised machine learning concepts and methods, including deep learning, which are potentially the most relevant for the medical domain. We aim to take some of the mystery out of machine learning and depict how machine learning models can be useful for medical applications. Finally, this tutorial provides a few practical suggestions for how to properly design a machine learning model for a generic medical problem.

Phase Error Correction in Time-Averaged 3D Phase Contrast Magnetic Resonance Imaging of the Cerebral Vasculature.

PURPOSE: Volume flow rate (VFR) measurements based on phase contrast (PC)-magnetic resonance (MR) imaging datasets have spatially varying bias due to eddy current induced phase errors. The purpose of this study was to assess the impact of phase errors in time averaged PC-MR imaging of the cerebral vasculature and explore the effects of three common correction schemes (local bias correction (LBC), local polynomial correction (LPC), and whole brain polynomial correction (WBPC)). METHODS: Measurements of the eddy current induced phase error from a static phantom were first obtained. In thirty healthy human subjects, the methods were then assessed in background tissue to determine if local phase offsets could be removed. Finally, the techniques were used to correct VFR measurements in cerebral vessels and compared statistically. RESULTS: In the phantom, phase error was measured to be <2.1 ml/s per pixel and the bias was reduced with the correction schemes. In background tissue, the bias was significantly reduced, by 65.6% (LBC), 58.4% (LPC) and 47.7% (WBPC) (p < 0.001 across all schemes). Correction did not lead to significantly different VFR measurements in the vessels (p = 0.997). In the vessel measurements, the three correction schemes led to flow measurement differences of -0.04 ± 0.05 ml/s, 0.09 ± 0.16 ml/s, and -0.02 ± 0.06 ml/s. Although there was an improvement in background measurements with correction, there was no statistical difference between the three correction schemes (p = 0.242 in background and p = 0.738 in vessels). CONCLUSIONS: While eddy current induced phase errors can vary between hardware and sequence configurations, our results showed that the impact is small in a typical brain PC-MR protocol and does not have a significant effect on VFR measurements in cerebral vessels.

Flow and pressure measurements in aneurysms and arteriovenous malformations with phase contrast MR imaging.

PURPOSE: To explore phase contrast (PC) magnetic resonance imaging of aneurysms and arteriovenous malformations (AVM). PC imaging obtains a vector field of the velocity and can yield additional hemodynamic information, including: volume flow rate (VFR) and intravascular pressure. We expect to find lower VFR distal to aneurysms and higher VFR in vessels of the AVM. MATERIALS AND METHODS: Five cerebral aneurysm and three AVM patients were imaged with PC techniques and compared to VFR of a healthy cohort. VFR was calculated in vessel segments in each patient and compared statistically to the healthy cohort by computing the z-score. Intravascular pressure was calculated in the aneurysms and in the nidus of each AVM. RESULTS: We found that patients with aneurysm had z<-0.48 in vessels distal to the aneurysm (reduced flow), while AVM patients had z>6 in some vessels supplying and draining the nidus (increased flow). Pressures in aneurysms were highly variable between subjects and location, while in the nidus of the AVM patients; pressure trended higher in larger AVMs. CONCLUSION: The study findings confirm the expectation of lower distal flow in aneurysm and higher arterial and venous flow in AVM patients.