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.

Multi-site genetic analysis of diffusion images and voxelwise heritability analysis: A pilot project of the ENIGMA-DTI working group

The ENIGMA (Enhancing NeuroImaging Genetics through Meta-Analysis) Consortium was set up to analyze brain measures and genotypes from multiple sites across the world to improve the power to detect genetic variants that influence the brain. Diffusion tensor imaging (DTI) yields quantitative measures sensitive to brain development and degeneration, and some common genetic variants may be associated with white matter integrity or connectivity. DTI measures, such as the fractional anisotropy (FA) of water diffusion, may be useful for identifying genetic variants that influence brain microstructure. However, genome-wide association studies (GWAS) require large populations to obtain sufficient power to detect and replicate significant effects, motivating a multi-site consortium effort. As part of an ENIGMA-DTI working group, we analyzed high-resolution FA images from multiple imaging sites across North America, Australia, and Europe, to address the challenge of harmonizing imaging data collected at multiple sites. Four hundred images of healthy adults aged 18-85 from four sites were used to create a template and corresponding skeletonized FA image as a common reference space. Using twin and pedigree samples of different ethnicities, we used our common template to evaluate the heritability of tract-derived FA measures. We show that our template is reliable for integrating multiple datasets by combining results through meta-analysis and unifying the data through exploratory mega-analyses. Our results may help prioritize regions of the FA map that are consistently influenced by additive genetic factors for future genetic discovery studies. Protocols and templates are publicly available at (http://enigma.loni.ucla.edu/ongoing/dti-working-group/).

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.

Automatic clustering of white matter fibers in brain diffusion mri with an application to genetics

To understand factors that affect brain connectivity and integrity, it is beneficial to automatically cluster white matter (WM) fibers into anatomically recognizable tracts. Whole brain tractography, based on diffusion-weighted MRI, generates vast sets of fibers throughout the brain; clustering them into consistent and recognizable bundles can be difficult as there are wide individual variations in the trajectory and shape of WM pathways. Here we introduce a novel automated tract clustering algorithm based on label fusion - a concept from traditional intensity-based segmentation. Streamline tractography generates many incorrect fibers, so our top-down approach extracts tracts consistent with known anatomy, by mapping multiple hand-labeled atlases into a new dataset. We fuse clustering results from different atlases, using a mean distance fusion scheme. We reliably extracted the major tracts from 105-gradient high angular resolution diffusion images (HARDI) of 198 young normal twins. To compute population statistics, we use a pointwise correspondence method to match, compare, and average WM tracts across subjects. We illustrate our method in a genetic study of white matter tract heritability in twins.

Comparison of fractional and geodesic anisotropy in diffusion tensor images of 90 monozygotic and dizygotic twins

We used diffusion tensor magnetic resonance imaging (DTI) to reveal the extent of genetic effects on brain fiber microstructure, based on tensor-derived measures, in 22 pairs of monozygotic (MZ) twins and 23 pairs of dizygotic (DZ) twins (90 scans). After Log-Euclidean denoising to remove rank-deficient tensors, DTI volumes were fluidly registered by high-dimensional mapping of co-registered MP-RAGE scans to a geometrically-centered mean neuroanatomical template. After tensor reorientation using the strain of the 3D fluid transformation, we computed two widely used scalar measures of fiber integrity: fractional anisotropy (FA), and geodesic anisotropy (GA), which measures the geodesic distance between tensors in the symmetric positive-definite tensor manifold. Spatial maps of intraclass correlations (r) between MZ and DZ twins were compared to compute maps of Falconer's heritability statistics, i.e. the proportion of population variance explainable by genetic differences among individuals. Cumulative distribution plots (CDF) of effect sizes showed that the manifold measure, GA, comparably the Euclidean measure, FA, in detecting genetic correlations. While maps were relatively noisy, the CDFs showed promise for detecting genetic influences on brain fiber integrity as the current sample expands.

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.

Strain- and host species-specific inflammasome activation, IL-1β release, and cell death in macrophages infected with uropathogenic Escherichia coli

Uropathogenic Escherichia coli (UPEC) is the main etiological agent of urinary tract infections (UTIs). Little is known about interactions between UPEC and the inflammasome, a key innate immune pathway. Here we show that UPEC strains CFT073 and UTI89 trigger inflammasome activation and lytic cell death in human macrophages. Several other UPEC strains, including two multidrug-resistant ST131 isolates, did not kill macrophages. In mouse macrophages, UTI89 triggered cell death only at a high multiplicity of infection, and CFT073-mediated inflammasome responses were completely NLRP3-dependent. Surprisingly, CFT073- and UTI89-mediated responses only partially depended on NLRP3 in human macrophages. In these cells, NLRP3 was required for interleukin-1β (IL-1β) maturation, but contributed only marginally to cell death. Similarly, caspase-1 inhibition did not block cell death in human macrophages. In keeping with such differences, the pore-forming toxin α-hemolysin mediated a substantial proportion of CFT073-triggered IL-1β secretion in mouse but not human macrophages. There was also a more substantial α-hemolysin-independent cell death response in human vs. mouse macrophages. Thus, in mouse macrophages, CFT073-triggered inflammasome responses are completely NLRP3-dependent, and largely α-hemolysin-dependent. In contrast, UPEC activates an NLRP3-independent cell death pathway and an α-hemolysin-independent IL-1β secretion pathway in human macrophages. This has important implications for understanding UTI in humans.

The importance of connexin hemichannels during chondroprogenitor cell differentiation in hydrogel versus microtissue culture models

Appropriate selection of scaffold architecture is a key challenge in cartilage tissue engineering. Gap junction-mediated intercellular contacts play important roles in precartilage condensation of mesenchymal cells. However, scaffold architecture could potentially restrict cell-cell communication and differentiation. This is particularly important when choosing the appropriate culture platform as well as scaffold-based strategy for clinical translation, that is, hydrogel or microtissues, for investigating differentiation of chondroprogenitor cells in cartilage tissue engineering. We, therefore, studied the influence of gap junction-mediated cell-cell communication on chondrogenesis of bone marrow-derived mesenchymal stromal cells (BM-MSCs) and articular chondrocytes. Expanded human chondrocytes and BM-MSCs were either (re-) differentiated in micromass cell pellets or encapsulated as isolated cells in alginate hydrogels. Samples were treated with and without the gap junction inhibitor 18-α glycyrrhetinic acid (18αGCA). DNA and glycosaminoglycan (GAG) content and gene expression levels (collagen I/II/X, aggrecan, and connexin 43) were quantified at various time points. Protein localization was determined using immunofluorescence, and adenosine-5'-triphosphate (ATP) was measured in conditioned media. While GAG/DNA was higher in alginate compared with pellets for chondrocytes, there were no differences in chondrogenic gene expression between culture models. Gap junction blocking reduced collagen II and extracellular ATP in all chondrocyte cultures and in BM-MSC hydrogels. However, differentiation capacity was not abolished completely by 18αGCA. Connexin 43 levels were high throughout chondrocyte cultures and peaked only later during BM-MSC differentiation, consistent with the delayed response of BM-MSCs to 18αGCA. Alginate hydrogels and microtissues are equally suited culture platforms for the chondrogenic (re-)differentiation of expanded human articular chondrocytes and BM-MSCs. Therefore, reducing direct cell-cell contacts does not affect in vitro chondrogenesis. However, blocking gap junctions compromises cell differentiation, pointing to a prominent role for hemichannel function in this process. Therefore, scaffold design strategies that promote an increasing distance between single chondroprogenitor cells do not restrict their differentiation potential in tissue-engineered constructs.

EVERSUN: A phase 2 trial of alternating sunitinib and everolimus as first-line therapy for advanced renal cell carcinoma

Background We hypothesised that alternating inhibitors of the vascular endothelial growth factor receptor (VEGFR) and mammalian target of rapamycin pathways would delay the development of resistance in advanced renal cell carcinoma (aRCC). Patients and methods A single-arm, two-stage, multicentre, phase 2 trial to determine the activity, feasibility, and safety of 12-week cycles of sunitinib 50 mg daily 4 weeks on / 2 weeks off, alternating with everolimus 10 mg daily for 5 weeks on / 1 week off, until disease progression or prohibitive toxicity in favourable or intermediate-risk aRCC. The primary end point was proportion alive and progression-free at 6 months (PFS6m). The secondary end points were feasibility, tumour response, overall survival (OS), and adverse events (AEs). The correlative objective was to assess biomarkers and correlate with clinical outcome. Results We recruited 55 eligible participants from September 2010 to August 2012. Demographics: mean age 61, 71% male, favourable risk 16%, intermediate risk 84%. Cycle 2 commenced within 14 weeks for 80% of participants; 64% received ≥22 weeks of alternating therapy; 78% received ≥22 weeks of any treatment. PFS6m was 29/55 (53%; 95% confidence interval [CI] 40% to 66%). Tumour response rate was 7/55 (13%; 95% CI 4% to 22%, all partial responses). After median follow-up of 20 months, 47 of 55 (86%) had progressed with a median progression-free survival of 8 months (95% CI 5–10), and 30 of 55 (55%) had died with a median OS of 17 months (95% CI 12–undefined). AEs were consistent with those expected for each single agent. No convincing prognostic biomarkers were identified. Conclusions The EVERSUN regimen was feasible and safe, but its activity did not meet pre-specified values to warrant further research. This supports the current approach of continuing anti-VEGF therapy until progression or prohibitive toxicity before changing treatment.

Innovation in Road Building: Removing Obstacles for Diffusion of Novel Building Products

How can obstacles to innovation be overcome in road construction? Using a focus group methodology, and based on two prior rounds of empirical work, the analysis in this chapter generates a set of four key solutions to two main construction innovation obstacles: (1) restrictive tender assessment and (2) disagreement over who carries the risk of new product failure. The four key solutions uncovered were: 1) pre-project product certification; 2) past innovation performance assessment; 3) earlier involvement of product suppliers and road asset operators; and 4) performance-based specifications. Additional research is suggested in order to illicit deeper insights into possible solutions to construction innovation obstacles, and should emphasise furthering the theoretical interpretation of empirical phenomena.