UNITY: A low-field magnetic resonance neuroimaging initiative to characterize neurodevelopment in low and middle-income settings.

Measures of physical growth, such as weight and height have long been the predominant outcomes for monitoring child health and evaluating interventional outcomes in public health studies, including those that may impact neurodevelopment. While physical growth generally reflects overall health and nutritional status, it lacks sensitivity and specificity to brain growth and developing cognitive skills and abilities. Psychometric tools, e.g., the Bayley Scales of Infant and Toddler Development, may afford more direct assessment of cognitive development but they require language translation, cultural adaptation, and population norming. Further, they are not always reliable predictors of future outcomes when assessed within the first 12-18 months of a child's life. Neuroimaging may provide more objective, sensitive, and predictive measures of neurodevelopment but tools such as magnetic resonance (MR) imaging are not readily available in many low and middle-income countries (LMICs). MRI systems that operate at lower magnetic fields (< 100mT) may offer increased accessibility, but their use for global health studies remains nascent. The UNITY project is envisaged as a global partnership to advance neuroimaging in global health studies. Here we describe the UNITY project, its goals, methods, operating procedures, and expected outcomes in characterizing neurodevelopment in sub-Saharan Africa and South Asia.

A novel MRI phantom to study interstitial fluid transport in the glymphatic system.

The glymphatic system is a recently discovered transport system, mediated by cerebral spinal fluid (CSF), that clears metabolic and cellular waste products in the brain. This system's function in the brain is analogous to that of the lymphatic system in the rest of the mammalian body. It is hypothesized that CSF clears harmful chemicals from the brain by flowing through interstitial spaces in the brain during sleep. While there is growing recognition of the critical role the glymphatic system plays in maintaining normal brain health and in explaining pathology, there are few noninvasive imaging methods that measure and characterize the efficacy of glymphatic transport in vivo. In this study we designed, constructed, and tested a glymphatic transport magnetic resonance imaging (MRI) flow phantom, which combines regions that mimic CSF-filled ventricles and brain interstitial space. We tested high- and low-q space diffusion MRI and diffusion tensor imaging (DTI) acquisitions to determine if they could detect, measure, and map interstitial glymphatic flows. The results suggest that, under certain flow conditions, diffusion-weighted MRI can detect the enhanced mixing that occurs during glymphatic clearance.

Anisotropic phantom to calibrate high-q diffusion MRI methods.

A silicon oil-filled glass capillary array is proposed as an anisotropic diffusion MRI phantom. Together with a computational/theoretical pipeline these provide a gold standard for calibrating and validating high-q diffusion MRI experiments. The phantom was used to test high angular resolution diffusion imaging (HARDI) and double pulsed-field gradient (d-PFG) MRI acquisition schemes. MRI-based predictions of microcapillary diameter using both acquisition schemes were compared with results from optical microscopy. This phantom design can be used for quality control and quality assurance purposes and for testing and validating proposed microstructure imaging experiments and the processing pipelines used to analyze them.

Role of myelin plasticity in oscillations and synchrony of neuronal activity.

Conduction time is typically ignored in computational models of neural network function. Here we consider the effects of conduction delays on the synchrony of neuronal activity and neural oscillators, and evaluate the consequences of allowing conduction velocity (CV) to be regulated adaptively. We propose that CV variation, mediated by myelin, could provide an important mechanism of activity-dependent nervous system plasticity. Even small changes in CV, resulting from small changes in myelin thickness or nodal structure, could have profound effects on neuronal network function in terms of spike-time arrival, oscillation frequency, oscillator coupling, and propagation of brain waves. For example, a conduction delay of 5ms could change interactions of two coupled oscillators at the upper end of the gamma frequency range (∼100Hz) from constructive to destructive interference; delays smaller than 1ms could change the phase by 30°, significantly affecting signal amplitude. Myelin plasticity, as another form of activity-dependent plasticity, is relevant not only to nervous system development but also to complex information processing tasks that involve coupling and synchrony among different brain rhythms. We use coupled oscillator models with time delays to explore the importance of adaptive time delays and adaptive synaptic strengths. The impairment of activity-dependent myelination and the loss of adaptive time delays may contribute to disorders where hyper- and hypo-synchrony of neuronal firing leads to dysfunction (e.g., dyslexia, schizophrenia, epilepsy).

Mapping average axon diameters in porcine spinal cord white matter and rat corpus callosum using d-PFG MRI.

Knowledge of microstructural features of nerve fascicles, such as their axon diameter, is crucial for understanding normal function in the central and peripheral nervous systems as well as assessing changes due to pathologies. In this study double-pulsed field gradient (d-PFG) filtered MRI was used to map the average axon diameter (AAD) in porcine spinal cord, which was then compared to AADs measured with optical microscopy of the same specimen, as a way to further validate this new MRI method. A novel 3D d-PFG acquisition scheme was used to obtain AADs in each voxel of a coronal slice of rat brain corpus callosum. AAD measurements were also acquired using optical microscopy performed on histological sections and validated using a glass capillary array phantom.

A spin echo sequence with a single-sided bipolar diffusion gradient pulse to obtain snapshot diffusion weighted images in moving media.

In vivo MRI data can be corrupted by motion. Motion artifacts are particularly troublesome in Diffusion Weighted MRI (DWI), since the MR signal attenuation due to Brownian motion can be much less than the signal loss due to dephasing from other types of complex tissue motion, which can significantly degrade the estimation of self-diffusion coefficients, diffusion tensors, etc. This paper describes a snapshot DWI sequence, which utilizes a novel single-sided bipolar diffusion sensitizing gradient pulse within a spin echo sequence. The proposed method shortens the diffusion time by applying a single refocused bipolar diffusion gradient on one side of a refocusing RF pulse, instead of a set of diffusion sensitizing gradients, separated by a refocusing RF pulse, while reducing the impact of magnetic field inhomogeneity by using a spin echo sequence. A novel MRI phantom that can exhibit a range of complex motions was designed to demonstrate the robustness of the proposed DWI sequence.

Determining which mechanisms lead to activation in the motor cortex: a modeling study of transcranial magnetic stimulation using realistic stimulus waveforms and sulcal geometry.

OBJECTIVE: To determine which mechanisms lead to activation of neurons in the motor cortex during transcranial magnetic stimulation (TMS) with different current directions and pulse waveforms. METHODS: The total electric field induced in a simplified model of a cortical sulcus by a figure-eight coil was calculated using the finite element method (FEM). This electric field was then used as the input to determine the response of compartmental models of several types of neurons. RESULTS: The modeled neurons were stimulated at different sites: fiber bends for pyramidal tract neurons, axonal terminations for cortical interneurons and axon collaterals, and a combination of both for pyramidal association fibers. All neurons were more easily stimulated by a PA - directed electric field, except association fibers. Additionally, the second phase of a biphasic pulse was found to be more efficient than the first phase of either monophasic or biphasic pulses. CONCLUSIONS: The stimulation threshold for different types of neurons depends on the pulse waveform and relative current direction. The reported results might account for the range of responses obtained in TMS of the motor cortex when using different stimulation parameters. SIGNIFICANCE: Modeling studies combining electric field calculations and neuronal models may lead to a deeper understanding of the effect of the TMS-induced electric field on cortical tissue, and may be used to improve TMS coil and waveform design.

Double pulsed field gradient (double-PFG) MR imaging (MRI) as a means to measure the size of plant cells.

Measurement of diffusion in porous materials and biological tissues with the pulsed field gradient (PFG) MR techniques has proven useful in characterizing the microstructure of such specimens noninvasively. A natural extension of the traditional PFG technique comprises multiple pairs of diffusion gradients. This approach has been shown to provide the ability to characterize anisotropy at different length scales without the need to employ very strong gradients. In this work, the double-PFG imaging technique was used on a specimen involving a series of glass capillary arrays with different diameters. The experiments on the phantom demonstrated the ability to create a quantitative and accurate map of pore sizes. The same technique was subsequently employed to image a celery stalk. A diffusion tensor image (DTI) of the same specimen was instrumental in accurately delineating the regions of vascular tissue and determining the local orientation of cells. This orientation information was incorporated into a theoretical double-PFG framework and the technique was employed to estimate the cell size in the vascular bundles of the celery stalk. The findings suggest that the double-PFG MRI framework could provide important new information regarding the microstructure of many plants and other food products.

Elucidating the mechanisms and loci of neuronal excitation by transcranial magnetic stimulation using a finite element model of a cortical sulcus.

OBJECTIVE: This work aims to elucidate by what physical mechanisms and where stimulation occurs in the brain during transcranial magnetic stimulation (TMS), taking into account cortical geometry and tissue heterogeneity. METHODS: An idealized computer model of TMS was developed, comprising a stimulation coil, a cortical sulcus, and surrounding tissues. The distribution of the induced electric field was computed, and estimates of the relevant parameters were generated to predict the locus and type of neurons stimulated during TMS, assuming three different stimulation mechanisms. RESULTS: Tissue heterogeneity strongly affects the spatial distribution of the induced electric field and hence which stimulation mechanism is dominant and where it acts. Stimulation of neurons may occur in the gyrus, in the lip of the gyrus, and in the walls of the sulcus. The stimulated cells can be either pyramidal cells having medium to large caliber axons, or intracortical fibers of medium caliber. CONCLUSIONS: The results highlight the influence of cortical folding on the action of magnetic and electric fields on cortical tissue. SIGNIFICANCE: Tissue geometry and heterogeneity in electrical conductivity both must be taken into account to predict accurately stimulation loci and mechanism in TMS.

Observation of microscopic diffusion anisotropy in the spinal cord using double-pulsed gradient spin echo MRI.

A double-pulsed gradient spin echo (d-PGSE) filtered MRI sequence is proposed to detect microscopic diffusion anisotropy in heterogeneous specimen. The technique was developed, in particular, to characterize local microscopic anisotropy in specimens that are macroscopically isotropic, such as gray matter. In such samples, diffusion tensor MRI (DTI) produces an isotropic or nearly isotropic diffusion tensor despite the fact that the medium may be anisotropic at a microscopic length scale. Using d-PGSE filtered MRI, microscopic anisotropy was observed in a "gray matter" phantom consisting of randomly oriented tubes filled with water, as well as in fixed pig spinal cord, within a range of b-values that can be readily achieved on clinical and small animal MR scanners. These findings suggest a potential use for this new contrast mechanism in clinical studies and biological research applications.

Detection of microscopic anisotropy in gray matter and in a novel tissue phantom using double Pulsed Gradient Spin Echo MR.

A double Pulsed Gradient Spin Echo (d-PGSE) MR experiment was used to measure and assess the degree of local diffusion anisotropy in brain gray matter, and in a novel "gray matter" phantom that consists of randomly oriented tubes filled with water. In both samples, isotropic diffusion was observed at a macroscopic scale while anisotropic diffusion was observed at a microscopic scale, however, the nature of the resulting echo attenuation profiles were qualitatively different. Gray matter, which contains multiple cell types and fibers, exhibits a more complicated echo attenuation profile than the phantom. Since microscopic anisotropy was observed in both samples in the low q regime comparable to that achievable in clinical scanner, it may offer a new potential contrast mechanism for characterizing gray matter microstructure in medical and biological applications.

Tissue heterogeneity as a mechanism for localized neural stimulation by applied electric fields.

We investigate the heterogeneity of electrical conductivity as a new mechanism to stimulate excitable tissues via applied electric fields. In particular, we show that stimulation of axons crossing internal boundaries can occur at boundaries where the electric conductivity of the volume conductor changes abruptly. The effectiveness of this and other stimulation mechanisms was compared by means of models and computer simulations in the context of transcranial magnetic stimulation. While, for a given stimulation intensity, the largest membrane depolarization occurred where an axon terminates or bends sharply in a high electric field region, a slightly smaller membrane depolarization, still sufficient to generate action potentials, also occurred at an internal boundary where the conductivity jumped from 0.143 S m(-1) to 0.333 S m(-1), simulating a white-matter-grey-matter interface. Tissue heterogeneity can also give rise to local electric field gradients that are considerably stronger and more focal than those impressed by the stimulation coil and that can affect the membrane potential, albeit to a lesser extent than the two mechanisms mentioned above. Tissue heterogeneity may play an important role in electric and magnetic 'far-field' stimulation.

Comprehensive approach for correction of motion and distortion in diffusion-weighted MRI.

Patient motion and image distortion induced by eddy currents cause artifacts in maps of diffusion parameters computed from diffusion-weighted (DW) images. A novel and comprehensive approach to correct for spatial misalignment of DW imaging (DWI) volumes acquired with different strengths and orientations of the diffusion sensitizing gradients is presented. This approach uses a mutual information-based registration technique and a spatial transformation model containing parameters that correct for eddy current-induced image distortion and rigid body motion in three dimensions. All parameters are optimized simultaneously for an accurate and fast solution to the registration problem. The images can also be registered to a normalized template with a single interpolation step without additional computational cost. Following registration, the signal amplitude of each DWI volume is corrected to account for size variations of the object produced by the distortion correction, and the b-matrices are properly recalculated to account for any rotation applied during registration. Both qualitative and quantitative results show that this approach produces a significant improvement of diffusion tensor imaging (DTI) data acquired in the human brain.

Cartilage calcification studied by proton nuclear magnetic resonance microscopy.

A three-dimensional (3D) mineralizing culture system using hollow fiber bioreactors has been developed to study the early stages of endochondral ossification by proton nuclear magnetic resonance (NMR) microscopy. Chondrocytes harvested from the cephalic half of the sterna from 17-day-old chick embryos were terminally differentiated with 33 nM of retinoic acid for 1 week and mineralization was initiated by the addition of 1% beta-glycerophosphate to the culture medium. Histological sections taken after 6 weeks of development in culture confirmed calcification of the cartilage matrix formed in bioreactors. Calcium to phosphorus ratios (1.62-1.68) from X-ray microanalysis supported electron diffraction of thin tissue sections showing the presence of a poorly crystalline hydroxyapatite mineral phase in the cultures. After 4 weeks of culture, quantitative proton NMR images showed water proton magnetization transfer rate constants (km) were higher in premineralized cartilage compared with uncalcified cartilage, a result suggesting collagen enrichment of the matrix. Notably after 5 weeks mineral deposits formed in bioreactors principally in the collagen-enriched zones of the cartilage with increased km values. This caused marked reductions in water proton longitudinal (T1) and transverse (T2) relaxation times and water diffusion coefficients (D). These results support the hypothesis that mineralization proceeds in association with a collagen template. After 6 weeks of culture development, the water proton T2 values decreased by 13% and D increased by 7% in uncalcified areas, compared with the same regions of tissue examined 1 week earlier. These changes could be attributed to the formation of small mineral inclusions in the cartilage, possibly mediated by matrix vesicles, which may play an important role in cartilage calcification. In summary, NMR images acquired before and after the onset of mineralization of the same tissue provide unique insights into the matrix events leading to endochondral mineral formation.

Effect of monovalent-divalent cation exchange on the swelling of polyacrylate hydrogels in physiological salt solutions.

The volume transition induced by monovalent-divalent cation exchange of fully neutralized polyacrylate hydrogels was investigated in aqueous NaCl solutions. The variation of the osmotic swelling pressure, shear modulus, and mixing pressure was measured when Na(+) ions were substituted by divalent or trivalent cations. Alkali metal salts move freely throughout the entirely network, and alkaline earth metal salts (CaCl(2), SrCl(2)) promote aggregation of polyacrylate chains, but these aggregates are relatively weak. Transition metal salts (CoCl(2), NiCl(2)) form stronger interchain associates. Rare earth cations (La(3+) and Ce(3+)) bind practically irreversibly to the polymer. Experimental data indicate that transition metal cations modify both the elastic and mixing components of the free energy, while alkaline earth metal cations affect primarily the mixing term. The behavior of freely swollen gels was compared with similar gels subjected to uniaxial compression. In uniaxially compressed gels, volume transition occurs at lower cation concentrations than in the corresponding undeformed gels. The shift of the transition point increases with the deformation ratio and is larger for Co(2+) than for Ca(2+).

Spatial transformations of diffusion tensor magnetic resonance images.

We address the problem of applying spatial transformations (or "image warps") to diffusion tensor magnetic resonance images. The orientational information that these images contain must be handled appropriately when they are transformed spatially during image registration. We present solutions for global transformations of three-dimensional images up to 12-parameter affine complexity and indicate how our methods can be extended for higher order transformations. Several approaches are presented and tested using synthetic data. One method, the preservation of principal direction algorithm, which takes into account shearing, stretching and rigid rotation, is shown to be the most effective. Additional registration experiments are performed on human brain data obtained from a single subject, whose head was imaged in three different orientations within the scanner. All of our methods improve the consistency between registered and target images over naïve warping algorithms.

Analytical computation of the eigenvalues and eigenvectors in DT-MRI.

In this paper a noniterative algorithm to be used for the analytical determination of the sorted eigenvalues and corresponding orthonormalized eigenvectors obtained by diffusion tensor magnetic resonance imaging (DT-MRI) is described. The algorithm uses the three invariants of the raw water spin self-diffusion tensor represented by a 3 x 3 positive definite matrix and certain math functions that do not require iteration. The implementation requires a positive definite mask to preserve the physical meaning of the eigenvalues. This algorithm can increase the speed of eigenvalue/eigenvector calculations by a factor of 5-40 over standard iterative Jacobi or singular-value decomposition techniques. This approach may accelerate the computation of eigenvalues, eigenvalue-dependent metrics, and eigenvectors especially when having high-resolution measurements with large numbers of slices and large fields of view.

In vivo fiber tractography using DT-MRI data.

Fiber tract trajectories in coherently organized brain white matter pathways were computed from in vivo diffusion tensor magnetic resonance imaging (DT-MRI) data. First, a continuous diffusion tensor field is constructed from this discrete, noisy, measured DT-MRI data. Then a Frenet equation, describing the evolution of a fiber tract, was solved. This approach was validated using synthesized, noisy DT-MRI data. Corpus callosum and pyramidal tract trajectories were constructed and found to be consistent with known anatomy. The method's reliability, however, degrades where the distribution of fiber tract directions is nonuniform. Moreover, background noise in diffusion-weighted MRIs can cause a computed trajectory to hop from tract to tract. Still, this method can provide quantitative information with which to visualize and study connectivity and continuity of neural pathways in the central and peripheral nervous systems in vivo, and holds promise for elucidating architectural features in other fibrous tissues and ordered media.

Statistical artifacts in diffusion tensor MRI (DT-MRI) caused by background noise.

This work helps elucidate how background noise introduces statistical artifacts in the distribution of the sorted eigenvalues and eigenvectors in diffusion tensor MRI (DT-MRI) data. Although it was known that sorting eigenvalues (principal diffusivities) by magnitude introduces a bias in their sample mean within a homogeneous region of interest (ROI), here it is shown that magnitude sorting also introduces a significant bias in the variance of the sample mean eigenvalues. New methods are presented to calculate the mean and variance of the eigenvectors of the diffusion tensor, based on a dyadic tensor representation of eigenvalue-eigenvector pairs. Based on their use it is shown that sorting eigenvalues by magnitude also introduces a bias in the mean and the variance of the sample eigenvectors (principal directions). This required the development of new methods to calculate the mean and variance of the eigenvectors of the diffusion tensor, based on a dyadic tensor representation of eigenvalue-eigenvector pairs. Moreover, a new approach is proposed to order these pairs within an ROI. To do this, a correspondence between each principal axis of the diffusion ellipsoid, an eigenvalue-eigenvector pair, and a dyadic tensor constructed from it is exploited. A measure of overlap between principal axes of diffusion ellipsoids in different voxels is defined that employs projections between these dyadic tensors. The optimal eigenvalue assignment within an ROI maximizes this overlap. Bias in the estimate of the mean and of the variance of the eigenvalues and of their corresponding eigenvectors is reduced in DT-MRI experiments and in Monte Carlo simulations of such experiments. Improvement is most significant in isotropic regions, but some is also observed in anisotropic regions. This statistical framework should enhance our ability to characterize microstructure and architecture of healthy tissue, and help to assess its changes in development, disease, and degeneration. Mitigating these artifacts should also improve the characterization of diffusion anisotropy and the elucidation of fiber-tract trajectories in the brain and in other fibrous tissues. Magn Reson Med 44:41-50, 2000. Published 2000 Wiley-Liss, Inc.

Osmotic swelling of polyacrylate hydrogels in physiological salt solutions.

The swelling behavior of fully neutralized sodium polyacrylate gels was investigated in aqueous solutions of alkali metal (LiCl, NaCl, KCl, CsCl) and alkaline earth metal salts (CaCl2, SrCl2, BaCl2). The total salt concentration and the ratio of monovalent to divalent cations were varied in the biologically significant range. It is found that the concentrations of both monovalent and divalent cations vary continuously and smoothly in the gel despite the abrupt change in the gel volume. The individual elastic, mixing, and ionic contributions to the free energy of the gel were separately determined as a function of the degree of network swelling to elucidate the thermodynamics of swelling. Shear modulus measurements performed at different Ca2+ concentrations suggest that Ca2+ does not form stable cross-links between the polymer chains. At low and moderate swelling ratios the concentration dependence of the shear modulus follows a power law behavior, G variation of phi n, with n = 0.34 +/- 0.03. At high swelling degrees, however, the shear modulus increases with increasing swelling. The value of the Flory-Huggins interaction parameter, chi, determined from osmotic swelling pressure and shear modulus measurements, strongly depends on the ionic composition of the equilibrium solution and increases with increasing Ca2+ concentration.