Brain fiber architecture, genetics, and intelligence: a high angular resolution diffusion imaging (HARDI) study

We developed an analysis pipeline enabling population studies of HARDI data, and applied it to map genetic influences on fiber architecture in 90 twin subjects. We applied tensor-driven 3D fluid registration to HARDI, resampling the spherical fiber orientation distribution functions (ODFs) in appropriate Riemannian manifolds, after ODF regularization and sharpening. Fitting structural equation models (SEM) from quantitative genetics, we evaluated genetic influences on the Jensen-Shannon divergence (JSD), a novel measure of fiber spatial coherence, and on the generalized fiber anisotropy (GFA) a measure of fiber integrity. With random-effects regression, we mapped regions where diffusion profiles were highly correlated with subjects' intelligence quotient (IQ). Fiber complexity was predominantly under genetic control, and higher in more highly anisotropic regions; the proportion of genetic versus environmental control varied spatially. Our methods show promise for discovering genes affecting fiber connectivity in the brain.

Mapping genetic influences on brain fiber architecture with high angular resolution diffusion imaging (HARDI)

We report the first 3D maps of genetic effects on brain fiber complexity. We analyzed HARDI brain imaging data from 90 young adult twins using an information-theoretic measure, the Jensen-Shannon divergence (JSD), to gauge the regional complexity of the white matter fiber orientation distribution functions (ODF). HARDI data were fluidly registered using Karcher means and ODF square-roots for interpol ation; each subject's JSD map was computed from the spatial coherence of the ODFs in each voxel's neighborhood. We evaluated the genetic influences on generalized fiber anisotropy (GFA) and complexity (JSD) using structural equation models (SEM). At each voxel, genetic and environmental components of data variation were estimated, and their goodness of fit tested by permutation. Color-coded maps revealed that the optimal models varied for different brain regions. Fiber complexity was predominantly under genetic control, and was higher in more highly anisotropic regions. These methods show promise for discovering factors affecting fiber connectivity in the brain.

Genetics of brain fiber architecture and intellectual performance

The study is the first to analyze genetic and environmental factors that affect brain fiber architecture and its genetic linkage with cognitive function. We assessed white matter integrity voxelwise using diffusion tensor imaging at high magnetic field (4 Tesla), in 92 identical and fraternal twins. White matter integrity, quantified using fractional anisotropy (FA), was used to fit structural equation models (SEM) at each point in the brain, generating three-dimensional maps of heritability. We visualized the anatomical profile of correlations between white matter integrity and full-scale, verbal, and performance intelligence quotients (FIQ, VIQ, and PIQ). White matter integrity (FA) was under strong genetic control and was highly heritable in bilateral frontal (a 2 = 0.55, p = 0.04, left; a 2 = 0.74, p = 0.006, right), bilateral parietal (a 2 = 0.85, p < 0.001, left; a 2 = 0.84, p < 0.001, right), and left occipital (a 2 = 0.76, p = 0.003) lobes, and was correlated with FIQ and PIQ in the cingulum, optic radiations, superior fronto- occipital fasciculus, internal capsule, callosal isthmus, and the corona radiata (p = 0.04 for FIQ and p = 0.01 for PIQ, corrected for multiple comparisons). In a cross-trait mapping approach, common genetic factors mediated the correlation between IQ and white matter integrity, suggesting a common physiological mechanism for both, and common genetic determination. These genetic brain maps reveal heritable aspects of white matter integrity and should expedite the discovery of single-nucleotide polymorphisms affecting fiber connectivity and cognition.

Information-theoretic analysis of brain white matter fiber orientation distribution functions

We propose a new information-theoretic metric, the symmetric Kullback-Leibler divergence (sKL-divergence), to measure the difference between two water diffusivity profiles in high angular resolution diffusion imaging (HARDI). Water diffusivity profiles are modeled as probability density functions on the unit sphere, and the sKL-divergence is computed from a spherical harmonic series, which greatly reduces computational complexity. Adjustment of the orientation of diffusivity functions is essential when the image is being warped, so we propose a fast algorithm to determine the principal direction of diffusivity functions using principal component analysis (PCA). We compare sKL-divergence with other inner-product based cost functions using synthetic samples and real HARDI data, and show that the sKL-divergence is highly sensitive in detecting small differences between two diffusivity profiles and therefore shows promise for applications in the nonlinear registration and multisubject statistical analysis of HARDI data.

Heritability of white matter fiber tract shapes: A HARDI study of 198 twins

Genetic analysis of diffusion tensor images (DTI) shows great promise in revealing specific genetic variants that affect brain integrity and connectivity. Most genetic studies of DTI analyze voxel-based diffusivity indices in the image space (such as 3D maps of fractional anisotropy) and overlook tract geometry. Here we propose an automated workflow to cluster fibers using a white matter probabilistic atlas and perform genetic analysis on the shape characteristics of fiber tracts. We apply our approach to large study of 4-Tesla high angular resolution diffusion imaging (HARDI) data from 198 healthy, young adult twins (age: 20-30). Illustrative results show heritability for the shapes of several major tracts, as color-coded maps.

The multivariate A/C/E model and the genetics of fiber architecture

We present a new algorithm to compute the voxel-wise genetic contribution to brain fiber microstructure using diffusion tensor imaging (DTI) in a dataset of 25 monozygotic (MZ) twins and 25 dizygotic (DZ) twin pairs (100 subjects total). First, the structural and DT scans were linearly co-registered. Structural MR scans were nonlinearly mapped via a 3D fluid transformation to a geometrically centered mean template, and the deformation fields were applied to the DTI volumes. After tensor re-orientation to realign them to the anatomy, we computed several scalar and multivariate DT-derived measures including the geodesic anisotropy (GA), the tensor eigenvalues and the full diffusion tensors. A covariance-weighted distance was measured between twins in the Log-Euclidean framework [2], and used as input to a maximum-likelihood based algorithm to compute the contributions from genetics (A), common environmental factors (C) and unique environmental ones (E) to fiber architecture. Quanititative genetic studies can take advantage of the full information in the diffusion tensor, using covariance weighted distances and statistics on the tensor manifold.

White matter integrity measured by fractional anisotropy correlates poorly with actual individual fiber anisotropy

Fractional anisotropy (FA), a very widely used measure of fiber integrity based on diffusion tensor imaging (DTI), is a problematic concept as it is influenced by several quantities including the number of dominant fiber directions within each voxel, each fiber's anisotropy, and partial volume effects from neighboring gray matter. High-angular resolution diffusion imaging (HARDI) can resolve more complex diffusion geometries than standard DTI, including fibers crossing or mixing. The tensor distribution function (TDF) can be used to reconstruct multiple underlying fibers per voxel, representing the diffusion profile as a probabilistic mixture of tensors. Here we found that DTIderived mean diffusivity (MD) correlates well with actual individual fiber MD, but DTI-derived FA correlates poorly with actual individual fiber anisotropy, and may be suboptimal when used to detect disease processes that affect myelination. Analysis of the TDFs revealed that almost 40% of voxels in the white matter had more than one dominant fiber present. To more accurately assess fiber integrity in these cases, we here propose the differential diffusivity (DD), which measures the average anisotropy based on all dominant directions in each voxel.

Brain network efficiency and topology depend on the fiber tracking method: 11 tractography algorithms compared in 536 subjects

As connectivity analyses become more popular, claims are often made about how the brain's anatomical networks depend on age, sex, or disease. It is unclear how results depend on tractography methods used to compute fiber networks. We applied 11 tractography methods to high angular resolution diffusion images of the brain (4-Tesla 105-gradient HARDI) from 536 healthy young adults. We parcellated 70 cortical regions, yielding 70×70 connectivity matrices, encoding fiber density. We computed popular graph theory metrics, including network efficiency, and characteristic path lengths. Both metrics were robust to the number of spherical harmonics used to model diffusion (4th-8th order). Age effects were detected only for networks computed with the probabilistic Hough transform method, which excludes smaller fibers. Sex and total brain volume affected networks measured with deterministic, tensor-based fiber tracking but not with the Hough method. Each tractography method includes different fibers, which affects inferences made about the reconstructed networks.

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.

Eye-safe UV stand-off Raman spectroscopy for the ranged detection of explosives in the field

Increasing worldwide terrorist attacks involving explosives presents a growing need for a rapid and ranged explosive detection method that can safely be deployed in the field. Stand-off Raman spectroscopy shows great promise; however, the radiant exposures of lasers required for adequate signal generation are often much greater than what is safe for the eye or the skin, restricting use of the technique to un-populated areas. Here, by determining the safe exposure levels for lasers typically used in Raman spectroscopy, optimal parameter values are identified, which produce the largest possible detection range using power densities that do not exceed the eye-safe limit. It is shown that safe ultraviolet pulse energies can be more than three orders of magnitude greater than equivalent safe visible pulse energies. Coupling this to the 16-fold increase in Raman signal obtained in the ultraviolet at 266 nm over that at 532 nm results in a 131 times larger detection range for the eye-safe 266-nm system over an equivalent eye-safe 532-nm laser system. For the Raman system described here, this translates to a maximum range of 42 m for detecting Teflon with a 266-nm laser emitting a 100-mm diameter beam of 23.5-mJ nanosecond pulses.

Development of a switch mode assisted linear amplifier for use as a high fidelity voltage source

Frequency Domain Spectroscopy (FDS) is one of the major techniques used for determining the condition of the cellulose based paper and pressboard components in large oil/paper insulated power transformers. This technique typically makes use of a sinusoidal voltage source swept from 0.1 mHz to 1 kHz. The excitation test voltage source used must meet certain characteristics, such as high output voltage, high fidelity, low noise and low harmonic content. The amplifier used; in the test voltage source; must be able to drive highly capacitive loads. This paper proposes that a switch-mode assisted linear amplifier (SMALA) can be used in the test voltage source to meet these criteria. A three level SMALA prototype amplifier was built to experimentally demonstrate the effectiveness of this proposal. The developed SMALA prototype shows no discernable harmonic distortion in the output voltage waveform, or the need for output filters, and is therefore seen as a preferable option to pulse width modulated digital amplifiers. The lack of harmonic distortion and high frequency switching noise in the output voltage of this SMALA prototype demonstrates its feasibility for applications in FDS, particularly on highly capacitive test objects such as transformer insulation systems.

Formation of highly conductive composite coatings and their applications to broadband antennas and mechanical transducers

Tight networks of interwoven carbon nanotube bundles are formed in our highly conductive composite. The composite possesses propertiessuggesting a two-dimensional percolative network rather than other reported dispersions displaying three-dimensional networks. Binding nanotubes into large but tight bundles dramatically alters the morphology and electronic transport dynamics of the composite. This enables itto carry higher levels of charge in the macroscale leading to conductivities as high as 1600 S/cm. We now discuss in further detail, the electronic and physical properties of the nanotube composites through Raman spectroscopy and transmission electron microscopy analysis. When controlled and usedappropriately, the interesting properties of these composites reveal their potential for practical device applications. For instance, we used this composite to fabricate coatings, whic improve the properties of an electromagnetic antenna/amplifier transducer. The resulting transducer possesses a broadband range up to GHz frequencies. A strain gauge transducer was also fabricated using changes in conductivity to monitor structural deformations in the composite coatings.

Fiber Bragg grating accelerometer based on a transversely rotating stick

Fiber Bragg Grating (FBG) accelerometers using transverse forces with an inertial object placed at the middle of the FBG have a high sensitivity but low resonant frequency. The resonant frequency 26 Hz and sensitivity at 6 Hz 1.29 nm/g were reported based on a 50mm-long FBG accelerometer. We demonstrate that the first FBG accelerometer based on a transversely rotating stick, which can, at the same or even larger size, keep the high sensitivity and significantly increase the low resonant frequency. In our experiments, a 77.5mm-long FBG accelerometer has achieved a similar sensitivity but 65% higher resonant frequency. This novel structure not only significantly widens the potential applications of FBG accelerometers by increasing their resonant frequencies but also provides a new route to design other accelerometers, e.g. micro accelerometers.

Trace vapour detection at room temperature using Raman spectroscopy

A miniaturized flow-through system consisting of a gold coated silicon substrate based on enhanced Raman spectroscopy has been used to study the detection of vapour from model explosive compounds. The measurements show that the detectability of the vapour molecules at room temperature depends sensitively on the interaction between the molecule and the substrate. The results highlight the capability of a flow system combined with Raman spectroscopy for detecting low vapour pressure compounds with a limit of detection of 0.2 ppb as demonstrated by the detection of bis(2-ethylhexyl)phthalate, a common polymer additive emitted from a commercial polyvinyl chloride (PVC) tubing at room temperature.

Raman spectroscopy techniques for detection of energetic materials at high and low concentration

This thesis presents the development of a rapid, sensitive and reproducible spectroscopic method for the detection of TNT in forensic and environmental applications. Simple nano sensors prepared by cost effective methods were utilized as sensitive platforms for the detection of TNT by surface enhanced Raman spectroscopy. The optimization of the substrate and the careful selection of a suitable recognition molecule contributed to the significant improvements of sensitive and selective targeting over current detection methods. The work presented in this thesis paves the way for effective detection and monitoring of explosives residues in law enforcement and environmental health applications.