Combining meta- and mega-analytic approaches for multi-site diffusion imaging based genetic studies: From the enigma-DTI working group

Meta-analyses estimate a statistical effect size for a test or an analysis by combining results from multiple studies without necessarily having access to each individual study's raw data. Multi-site meta-analysis is crucial for imaging genetics, as single sites rarely have a sample size large enough to pick up effects of single genetic variants associated with brain measures. However, if raw data can be shared, combining data in a "mega-analysis" is thought to improve power and precision in estimating global effects. As part of an ENIGMA-DTI investigation, we use fractional anisotropy (FA) maps from 5 studies (total N=2, 203 subjects, aged 9-85) to estimate heritability. We combine the studies through meta-and mega-analyses as well as a mixture of the two - combining some cohorts with mega-analysis and meta-analyzing the results with those of the remaining sites. A combination of mega-and meta-approaches may boost power compared to meta-analysis alone.

Diffusion imaging protocol effects on genetic associations

Large multi-site image-analysis studies have successfully discovered genetic variants that affect brain structure in tens of thousands of subjects scanned worldwide. Candidate genes have also associated with brain integrity, measured using fractional anisotropy in diffusion tensor images (DTI). To evaluate the heritability and robustness of DTI measures as a target for genetic analysis, we compared 417 twins and siblings scanned on the same day on the same high field scanner (4-Tesla) with two protocols: (1) 94-directions; 2mm-thick slices, (2) 27-directions; 5mm-thickness. Using mean FA in white matter ROIs and FA skeletons derived using FSL, we (1) examined differences in voxelwise means, variances, and correlations among the measures; and (2) assessed heritability with structural equation models, using the classical twin design. FA measures from the genu of the corpus callosum were highly heritable, regardless of protocol. Genome-wide analysis of the genu mean FA revealed differences across protocols in the top associations.

Assessment of dynamic susceptibility contrast cerebral blood flow response to amphetamine challenge: A human pharmacological magnetic resonance imaging study at 1.5 and 4 T

Pharmacological MRI (phMRI) techniques can be used to monitor the neurophysiological effects of central nervous system (CNS) active drugs. In this study, we investigated whether dynamic susceptibility contrast (DSC) perfusion imaging employing the use of superparamagnetic iron oxide nanoparticles (Resovist) could be used to measure hemodynamic response to d-amphetamine challenge in human subjects at both 1.5 and 4 T. Significant changes in cerebral blood flow (CBF) were found in focal regions associated with the nigrostriatal circuit and mesolimbic and mesocortical dopaminergic pathways. More significant CBF responses were found at higher field strength, mainly within striatal structures. The results from this study indicate that DSC perfusion imaging using Resovist can be used to assess the efficacy of CNS-active drugs and may play a role in the development of novel psychiatric therapies at the preclinical level. © 2005 Wiley-Liss, Inc.

Diffusion indices on magnetic resonance imaging and neuropsychological performance in amnestic mild cognitive impairment

Background: Magnetic resonance diffusion tensor imaging (DTI) shows promise in the early detection of microstructural pathophysiological changes in the brain. Objectives: To measure microstructural differences in the brains of participants with amnestic mild cognitive impairment (MCI) compared with an age-matched control group using an optimised DTI technique with fully automated image analysis tools and to investigate the correlation between diffusivity measurements and neuropsychological performance scores across groups. Methods: 34 participants (17 participants with MCI, 17 healthy elderly adults) underwent magnetic resonance imaging (MRI)-based DTI. To control for the effects of anatomical variation, diffusion images of all participants were registered to standard anatomical space. Significant statistical differences in diffusivity measurements between the two groups were determined on a pixel-by-pixel basis using gaussian random field theory. Results: Significantly raised mean diffusivity measurements (p<0.001) were observed in the left and right entorhinal cortices (BA28), posterior occipital-parietal cortex (BA18 and BA19), right parietal supramarginal gyrus (BA40) and right frontal precentral gyri (BA4 and BA6) in participants with MCI. With respect to fractional anisotropy, participants with MCI had significantly reduced measurements (p<0.001) in the limbic parahippocampal subgyral white matter, right thalamus and left posterior cingulate. Pearson's correlation coefficients calculated across all participants showed significant correlations between neuropsychological assessment scores and regional measurements of mean diffusivity and fractional anisotropy. Conclusions: DTI-based diffusivity measures may offer a sensitive method of detecting subtle microstructural brain changes associated with preclinical Alzheimer's disease.

Visualization tools for high angular resolution diffusion imaging

There is a major effort in medical imaging to develop algorithms to extract information from DTI and HARDI, which provide detailed information on brain integrity and connectivity. As the images have recently advanced to provide extraordinarily high angular resolution and spatial detail, including an entire manifold of information at each point in the 3D images, there has been no readily available means to view the results. This impedes developments in HARDI research, which need some method to check the plausibility and validity of image processing operations on HARDI data or to appreciate data features or invariants that might serve as a basis for new directions in image segmentation, registration, and statistics. We present a set of tools to provide interactive display of HARDI data, including both a local rendering application and an off-screen renderer that works with a web-based viewer. Visualizations are presented after registration and averaging of HARDI data from 90 human subjects, revealing important details for which there would be no direct way to appreciate using conventional display of scalar images.

Investigating brain connectivity heritability in a twin study using diffusion imaging data

Heritability of brain anatomical connectivity has been studied with diffusion-weighted imaging (DWI) mainly by modeling each voxel's diffusion pattern as a tensor (e.g., to compute fractional anisotropy), but this method cannot accurately represent the many crossing connections present in the brain. We hypothesized that different brain networks (i.e., their component fibers) might have different heritability and we investigated brain connectivity using High Angular Resolution Diffusion Imaging (HARDI) in a cohort of twins comprising 328 subjects that included 70 pairs of monozygotic and 91 pairs of dizygotic twins. Water diffusion was modeled in each voxel with a Fiber Orientation Distribution (FOD) function to study heritability for multiple fiber orientations in each voxel. Precision was estimated in a test-retest experiment on a sub-cohort of 39 subjects. This was taken into account when computing heritability of FOD peaks using an ACE model on the monozygotic and dizygotic twins. Our results confirmed the overall heritability of the major white matter tracts but also identified differences in heritability between connectivity networks. Inter-hemispheric connections tended to be more heritable than intra-hemispheric and cortico-spinal connections. The highly heritable tracts were found to connect particular cortical regions, such as medial frontal cortices, postcentral, paracentral gyri, and the right hippocampus.

How does angular resolution affect diffusion imaging measures?

A key question in diffusion imaging is how many diffusion-weighted images suffice to provide adequate signal-to-noise ratio (SNR) for studies of fiber integrity. Motion, physiological effects, and scan duration all affect the achievable SNR in real brain images, making theoretical studies and simulations only partially useful. We therefore scanned 50 healthy adults with 105-gradient high-angular resolution diffusion imaging (HARDI) at 4T. From gradient image subsets of varying size (6 ≤ N ≤ 94) that optimized a spherical angular distribution energy, we created SNR plots (versus gradient numbers) for seven common diffusion anisotropy indices: fractional and relative anisotropy (FA, RA), mean diffusivity (MD), volume ratio (VR), geodesic anisotropy (GA), its hyperbolic tangent (tGA), and generalized fractional anisotropy (GFA). SNR, defined in a region of interest in the corpus callosum, was near-maximal with 58, 66, and 62 gradients for MD, FA, and RA, respectively, and with about 55 gradients for GA and tGA. For VR and GFA, SNR increased rapidly with more gradients. SNR was optimized when the ratio of diffusion-sensitized to non-sensitized images was 9.13 for GA and tGA, 10.57 for FA, 9.17 for RA, and 26 for MD and VR. In orientation density functions modeling the HARDI signal as a continuous mixture of tensors, the diffusion profile reconstruction accuracy rose rapidly with additional gradients. These plots may help in making trade-off decisions when designing diffusion imaging protocols.

Analyzing multi-fiber reconstruction in high angular resolution diffusion imaging using the tensor distribution function

High-angular resolution diffusion imaging (HARDI) can reconstruct fiber pathways in the brain with extraordinary detail, identifying anatomical features and connections not seen with conventional MRI. HARDI overcomes several limitations of standard diffusion tensor imaging, which fails to model diffusion correctly in regions where fibers cross or mix. As HARDI can accurately resolve sharp signal peaks in angular space where fibers cross, we studied how many gradients are required in practice to compute accurate orientation density functions, to better understand the tradeoff between longer scanning times and more angular precision. We computed orientation density functions analytically from tensor distribution functions (TDFs) which model the HARDI signal at each point as a unit-mass probability density on the 6D manifold of symmetric positive definite tensors. In simulated two-fiber systems with varying Rician noise, we assessed how many diffusionsensitized gradients were sufficient to (1) accurately resolve the diffusion profile, and (2) measure the exponential isotropy (EI), a TDF-derived measure of fiber integrity that exploits the full multidirectional HARDI signal. At lower SNR, the reconstruction accuracy, measured using the Kullback-Leibler divergence, rapidly increased with additional gradients, and EI estimation accuracy plateaued at around 70 gradients.

Lateralization of temporal lobe epilepsy based on resting-state functional magnetic resonance imaging and machine learning

Lateralization of temporal lobe epilepsy (TLE) is critical for successful outcome of surgery to relieve seizures. TLE affects brain regions beyond the temporal lobes and has been associated with aberrant brain networks, based on evidence from functional magnetic resonance imaging. We present here a machine learning-based method for determining the laterality of TLE, using features extracted from resting-state functional connectivity of the brain. A comprehensive feature space was constructed to include network properties within local brain regions, between brain regions, and across the whole network. Feature selection was performed based on random forest and a support vector machine was employed to train a linear model to predict the laterality of TLE on unseen patients. A leave-one-patient-out cross validation was carried out on 12 patients and a prediction accuracy of 83% was achieved. The importance of selected features was analyzed to demonstrate the contribution of resting-state connectivity attributes at voxel, region, and network levels to TLE lateralization.

Meta-analysis of structural imaging findings in Attention-Deficit/Hyperactivity Disorder

Background Although there are many structural neuroimaging studies of attention-deficit/hyperactivity disorder (ADHD) in children, there are inconsistencies across studies and no consensus regarding which brain regions show the most robust area or volumetric reductions relative to control subjects. Our goal was to statistically analyze structural imaging data via a meta-analysis to help resolve these issues. Methods We searched the MEDLINE and PsycINFO databases through January 2005. Studies must have been written in English, used magnetic resonance imaging, and presented the means and standard deviations of regions assessed. Data were extracted by one of the authors and verified independently by another author. Results Analyses were performed using STATA with metan, metabias, and metainf programs. A meta-analysis including all regions across all studies indicated global reductions for ADHD subjects compared with control subjects, standardized mean difference equal to .408, p less than .001. Regions most frequently assessed and showing the largest differences included cerebellar regions, the splenium of the corpus callosum, total and right cerebral volume, and right caudate. Several frontal regions assessed in only two studies also showed large significant differences. Conclusions This meta-analysis provides a quantitative analysis of neuroanatomical abnormalities in ADHD and information that can be used to guide future studies.

Multimodal imaging of the collagen and elastic fibre networks in the bovine intervertebral disc

INTRODUCTION. The intervertebral disc is the largest avascular structure in the human body, withstanding transient loads of up to nine times body weight during rigorous physical activity. The key structural elements of the disc are a gel-like nucleus pulposus surrounded by concentric lamellar rings containing criss-crossed collagen fibres. The disc also contains an elastic fiber network which has been suggested to play a structural role, but to date the relationship between the collagen and elastic fiber networks is unclear. CONCLUSION. The multimodal transmitted and reflected polarized light microscopy technique developed here allows clear differentiation between the collagen and elastic fiber networks of the intervertebral disc. The ability to image unstained specimens avoids concerns with uneven stain penetration or specificity of staining. In bovine tail discs, the elastic fiber network is intimately associated with the collagen network.

Spinal inflammation in the absence of sacroiliac joint inflammation on magnetic resonance imaging in patients with active nonradiographic axial spondyloarthritis

Objective: To evaluate the presence of spinal inflammation with and without sacroiliac (SI) joint inflammation on magnetic resonance imaging (MRI) in patients with active nonradiographic axial spondyloarthritis (SpA), and to compare the disease characteristics of these subgroups. Methods: ABILITY-1 is a multicenter, randomized, controlled trial of adalimumab versus placebo in patients with nonradiographic axial SpA classified using the Assessment of SpondyloArthritis international Society axial SpA criteria. Baseline MRIs were centrally scored independently by 2 readers using the Spondyloarthritis Research Consortium of Canada (SPARCC) method for the SI joints and the SPARCC 6-discovertebral unit method for the spine. Positive evidence of inflammation on MRI was defined as a SPARCC score of >2 for either the SI joints or the spine. Results: Among patients with baseline SPARCC scores, 40% had an SI joint score of >2 and 52% had a spine score of >2. Forty-nine percent of patients with baseline SI joint scores of <2, and 58% of those with baseline SI joint scores of >2, had a spine score of >2. Comparison of baseline disease characteristics by baseline SI joint and spine scores showed that a greater proportion of patients in the subgroup with a baseline SPARCC score of >2 for both SI joints and spine were male, and patients with spine and SI joint scores of <2 were younger and had shorter symptom duration. SPARCC spine scores correlated with baseline symptom duration, and SI joint scores correlated negatively with the baseline Bath Ankylosing Spondylitis Disease Activity Index, but neither correlated with the baseline Ankylosing Spondylitis Disease Activity Score, total back pain, the patient's global assessment of disease activity, the Bath Ankylosing Spondylitis Functional Index, morning stiffness, nocturnal pain, or C-reactive protein level. Conclusion: Assessment by experienced readers showed that spinal inflammation on MRI might be observed in half of patients with nonradiographic axial SpA without SI joint inflammation.

SERS-barcoded colloidal gold NP assemblies as imaging agents for use in biodiagnostics

There is a growing need for new biodiagnostics that combine high throughput with enhanced spatial resolution and sensitivity. Gold nanoparticle (NP) assemblies with sub-10 nm particle spacing have the benefits of improving detection sensitivity via Surface enhanced Raman scattering (SERS) and being of potential use in biomedicine due to their colloidal stability. A promising and versatile approach to form solution-stable NP assemblies involves the use of multi-branched molecular linkers which allows tailoring of the assembly size, hot-spot density and interparticle distance. We have shown that linkers with multiple anchoring end-groups can be successfully employed as a linker to assemble gold NPs into dimers, linear NP chains and clustered NP assemblies. These NP assemblies with diameters of 30-120 nm are stable in solution and perform better as SERS substrates compared with single gold NPs, due to an increased hot-spot density. Thus, tailored gold NP assemblies are potential candidates for use as biomedical imaging agents. We observed that the hot-spot density and in-turn the SERS enhancement is a function of the linker polymer concentration and polymer architecture. New deep Raman techniques like Spatially Offset Raman Spectroscopy (SORS) have emerged that allow detection from beneath diffusely scattering opaque materials, including biological media such as animal tissue. We have been able to demonstrate that the gold NP assemblies could be detected from within both proteinaceous and high lipid containing animal tissue by employing a SORS technique with a backscattered geometry.

Analysis of focal adhesion turnover: A quantitative live-cell imaging example

Recent advances in optical and fluorescent protein technology have rapidly raised expectations in cell biology, allowing quantitative insights into dynamic intracellular processes like never before. However, quantitative live-cell imaging comes with many challenges including how best to translate dynamic microscopy data into numerical outputs that can be used to make meaningful comparisons rather than relying on representative data sets. Here, we use analysis of focal adhesion turnover dynamics as a straightforward specific example on how to image, measure, and analyze intracellular protein dynamics, but we believe this outlines a thought process and can provide guidance on how to understand dynamic microcopy data of other intracellular structures.

Advanced CUBIC protocols for whole-brain and whole-body clearing and imaging

Here we describe a protocol for advanced CUBIC (Clear, Unobstructed Brain/Body Imaging Cocktails and Computational analysis). The CUBIC protocol enables simple and efficient organ clearing, rapid imaging by light-sheet microscopy and quantitative imaging analysis of multiple samples. The organ or body is cleared by immersion for 1–14 d, with the exact time required dependent on the sample type and the experimental purposes. A single imaging set can be completed in 30–60 min. Image processing and analysis can take <1 d, but it is dependent on the number of samples in the data set. The CUBIC clearing protocol can process multiple samples simultaneously. We previously used CUBIC to image whole-brain neural activities at single-cell resolution using Arc-dVenus transgenic (Tg) mice. CUBIC informatics calculated the Venus signal subtraction, comparing different brains at a whole-organ scale. These protocols provide a platform for organism-level systems biology by comprehensively detecting cells in a whole organ or body.