A Monte Carlo-Based Fiber Tracking Algorithm using Diffusion Tensor MRI

Diffusion tensor magnetic resonance imaging, which measures directional information of water diffusion in the brain, has emerged as a powerful tool for human brain studies. In this paper, we introduce a new Monte Carlo-based fiber tracking approach to estimate brain connectivity. One of the main characteristics of this approach is that all parameters of the algorithm are automatically determined at each point using the entropy of the eigenvalues of the diffusion tensor. Experimental results show the good performance of the proposed approach

Tensor LRR and sparse coding-based subspace clustering

Subspace clustering groups a set of samples from a union of several linear subspaces into clusters, so that the samples in the same cluster are drawn from the same linear subspace. In the majority of the existing work on subspace clustering, clusters are built based on feature information, while sample correlations in their original spatial structure are simply ignored. Besides, original high-dimensional feature vector contains noisy/redundant information, and the time complexity grows exponentially with the number of dimensions. To address these issues, we propose a tensor low-rank representation (TLRR) and sparse coding-based (TLRRSC) subspace clustering method by simultaneously considering feature information and spatial structures. TLRR seeks the lowest rank representation over original spatial structures along all spatial directions. Sparse coding learns a dictionary along feature spaces, so that each sample can be represented by a few atoms of the learned dictionary. The affinity matrix used for spectral clustering is built from the joint similarities in both spatial and feature spaces. TLRRSC can well capture the global structure and inherent feature information of data, and provide a robust subspace segmentation from corrupted data. Experimental results on both synthetic and real-world data sets show that TLRRSC outperforms several established state-of-the-art methods.

Interpretation of the pressuremeter test using numerical models based on deformation tensor equations

The pressuremeter test in boreholes has proven itself as a useful tool in geotechnical explorations, especially comparing its results with those obtained from a mathematical model ruled by a soil representative constitutive equation. The numerical model shown in this paper is aimed to be the reference framework for the interpretation of this test. The model analyses variables such as: the type of response, the initial state, the drainage regime and the constitutive equations. It is a model of finite elements able to work with a mesh without deformation or one adapted to it.

Incremental tensor biased discriminant analysis:a new color-based visual tracking method

Most existing color-based tracking algorithms utilize the statistical color information of the object as the tracking clues, without maintaining the spatial structure within a single chromatic image. Recently, the researches on the multilinear algebra provide the possibility to hold the spatial structural relationship in a representation of the image ensembles. In this paper, a third-order color tensor is constructed to represent the object to be tracked. Considering the influence of the environment changing on the tracking, the biased discriminant analysis (BDA) is extended to the tensor biased discriminant analysis (TBDA) for distinguishing the object from the background. At the same time, an incremental scheme for the TBDA is developed for the tensor biased discriminant subspace online learning, which can be used to adapt to the appearance variant of both the object and background. The experimental results show that the proposed method can track objects precisely undergoing large pose, scale and lighting changes, as well as partial occlusion. © 2009 Elsevier B.V.

Gender identification of five genera of stingless bees (Apidae, Meliponini) based on wing morphology

Currently, the identification of pollinators is a critical necessity of conservation programs. After it was found that features extracted from patterns of wing venation are sufficient to discriminate among insect species, various studies have focused on this structure. We examined wing venation patterns of males and workers of five stingless bee species in order to determine if there are differences between sexes and if these differences are greater within than between species. Geometric morphometric analyses were made of the forewings of males and workers of Nannotrigona testaceicornis, Melipona quadrifasciata, Frieseomelitta varia, and Scaptotrigona aff. depilis and Plebeia remota. The patterns of males and workers from the same species were more similar than the patterns of individuals of the same sex from different species, and the patterns of both males and workers, when analyzed alone, were sufficiently different to distinguish among these five species. This demonstrates that we can use this kind of analysis for the identification of stingless bee species and that the sex of the individual does not impede identification. Computer-assisted morphometric analysis of bee wing images can be a useful tool for biodiversity studies and conservation programs.

Fast Structure-Based Assignment of 15N HSQC Spectra of Selectively 15N-Labeled Paramagnetic Proteins

A novel strategy for fast NMR resonance assignment of N-15 HSQC spectra of proteins is presented. It requires the structure coordinates of the protein, a paramagnetic center, and one or more residue-selectively N-15-labeled samples. Comparison of sensitive undecoupled N-15 HSQC spectra recorded of paramagnetic and diamagnetic samples yields data for every cross-peak on pseudocontact shift, paramagnetic relaxation enhancement, cross-correlation between Curie-spin and dipole-dipole relaxation, and residual dipolar coupling. Comparison of these four different paramagnetic quantities with predictions from the three-dimensional structure simultaneously yields the resonance assignment and the anisotropy of the susceptibility tensor of the paramagnetic center. The method is demonstrated with the 30 kDa complex between the N-terminal domain of the epsilon subunit and the theta subunit of Escherichia Coll DNA polymerase III. The program PLATYPUS was developed to perform the assignment, provide a measure of reliability of the assignment, and determine the susceptibility tensor anisotropy.

Lateral ventricle differences between first-episode schizophrenia and first-episode psychotic bipolar disorder: A population-based morphometric MRI study

Objectives. The extent to which psychotic disorders fall into distinct diagnostic categories or can be regarded as lying on a single continuum is controversial. We compared lateral ventricle volumes between a large sample of patients with first-episode schizophrenia or bipolar disorder and a healthy control group from the same neighbourhood. Methods. Population-based MRI study with 88 first-episode psychosis (FEP) patients, grouped into those with schizophrenia/schizophreniform disorder (N = 62), bipolar disorder (N = 26) and 94 controls. Results. Right and left lateral ventricular and right temporal horn volumes were larger in FEP subjects than controls. Within the FEP sample, post-hoc tests revealed larger left lateral ventricles and larger right and left temporal horns in schizophrenia subjects relative to controls, while there was no difference between patients with bipolar disorder and controls. None of the findings was attributable to effects of antipsychotics. Conclusions. This large-sample population-based MRI study showed that neuroanatomical abnormalities in subjects with schizophrenia relative to controls from the same neighbourhood are evident at the first episode of illness, but are not detectable in bipolar disorder patients. These data are consistent with a model of psychosis in which early brain insults of neurodevelopmental origin are more relevant to schizophrenia than to bipolar disorder.

Structural connectivity based on diffusion Kurtosis imaging

Structural connectivity models based on Diffusion Tensor Imaging (DTI) are strongly affected by the technique’s inability to resolve crossing fibres, either intra- or inter-hemispherical connections. Several models have been proposed to address this issue, including an algorithm aiming to resolve crossing fibres which is based on Diffusion Kurtosis Imaging (DKI). This technique is clinically feasible, even when multi-band acquisitions are not available, and compatible with multi-shell acquisition schemes. DKI is an extension of DTI enabling the estimation of diffusion tensor and diffusion kurtosis metrics. In this study we compare the performance of DKI and DTI in performing structural brain connectivity. Six healthy subjects were recruited, aged between 25 and 35 (three females). The MRI experiments were performed using a 3T Siemens Trio with a 32-channel head coil. The scans included a T1-weighted sequence (1mm3), and a DWI with b-values 0, 1000 and 2000 s:mm

Multispectral brain morphometry in Tourette syndrome persisting into adulthood.

Tourette syndrome is a childhood-onset neuropsychiatric disorder with a high prevalence of attention deficit hyperactivity and obsessive-compulsive disorder co-morbidities. Structural changes have been found in frontal cortex and striatum in children and adolescents. A limited number of morphometric studies in Tourette syndrome persisting into adulthood suggest ongoing structural alterations affecting frontostriatal circuits. Using cortical thickness estimation and voxel-based analysis of T1- and diffusion-weighted structural magnetic resonance images, we examined 40 adults with Tourette syndrome in comparison with 40 age- and gender-matched healthy controls. Patients with Tourette syndrome showed relative grey matter volume reduction in orbitofrontal, anterior cingulate and ventrolateral prefrontal cortices bilaterally. Cortical thinning extended into the limbic mesial temporal lobe. The grey matter changes were modulated additionally by the presence of co-morbidities and symptom severity. Prefrontal cortical thickness reduction correlated negatively with tic severity, while volume increase in primary somatosensory cortex depended on the intensity of premonitory sensations. Orbitofrontal cortex volume changes were further associated with abnormal water diffusivity within grey matter. White matter analysis revealed changes in fibre coherence in patients with Tourette syndrome within anterior parts of the corpus callosum. The severity of motor tics and premonitory urges had an impact on the integrity of tracts corresponding to cortico-cortical and cortico-subcortical connections. Our results provide empirical support for a patho-aetiological model of Tourette syndrome based on developmental abnormalities, with perturbation of compensatory systems marking persistence of symptoms into adulthood. We interpret the symptom severity related grey matter volume increase in distinct functional brain areas as evidence of ongoing structural plasticity. The convergence of evidence from volume and water diffusivity imaging strengthens the validity of our findings and attests to the value of a novel multimodal combination of volume and cortical thickness estimations that provides unique and complementary information by exploiting their differential sensitivity to structural change.

Correlation of diffusion tensor tractography and intraoperative macrostimulation during deep brain stimulation for Parkinson disease.

Object The purpose of this study was to investigate whether diffusion tensor imaging (DTI) of the corticospinal tract (CST) is a reliable surrogate for intraoperative macrostimulation through the deep brain stimulation (DBS) leads. The authors hypothesized that the distance on MRI from the DBS lead to the CST as determined by DTI would correlate with intraoperative motor thresholds from macrostimulations through the same DBS lead. Methods The authors retrospectively reviewed pre- and postoperative MRI studies and intraoperative macrostimulation recordings in 17 patients with Parkinson disease (PD) treated by DBS stimulation. Preoperative DTI tractography of the CST was coregistered with postoperative MRI studies showing the position of the DBS leads. The shortest distance and the angle from each contact of each DBS lead to the CST was automatically calculated using software-based analysis. The distance measurements calculated for each contact were evaluated with respect to the intraoperative voltage thresholds that elicited a motor response at each contact. Results There was a nonsignificant trend for voltage thresholds to increase when the distances between the DBS leads and the CST increased. There was a significant correlation between the angle and the voltage, but the correlation was weak (coefficient of correlation [R] = 0.36). Conclusions Caution needs to be exercised when using DTI tractography information to guide DBS lead placement in patients with PD. Further studies are needed to compare DTI tractography measurements with other approaches such as microelectrode recordings and conventional intraoperative MRI-guided placement of DBS leads.

Regional specificity of MRI contrast parameter changes in normal ageing revealed by voxel-based quantification (VBQ).

Normal ageing is associated with characteristic changes in brain microstructure. Although in vivo neuroimaging captures spatial and temporal patterns of age-related changes of anatomy at the macroscopic scale, our knowledge of the underlying (patho)physiological processes at cellular and molecular levels is still limited. The aim of this study is to explore brain tissue properties in normal ageing using quantitative magnetic resonance imaging (MRI) alongside conventional morphological assessment. Using a whole-brain approach in a cohort of 26 adults, aged 18-85years, we performed voxel-based morphometric (VBM) analysis and voxel-based quantification (VBQ) of diffusion tensor, magnetization transfer (MT), R1, and R2* relaxation parameters. We found age-related reductions in cortical and subcortical grey matter volume paralleled by changes in fractional anisotropy (FA), mean diffusivity (MD), MT and R2*. The latter were regionally specific depending on their differential sensitivity to microscopic tissue properties. VBQ of white matter revealed distinct anatomical patterns of age-related change in microstructure. Widespread and profound reduction in MT contrasted with local FA decreases paralleled by MD increases. R1 reductions and R2* increases were observed to a smaller extent in overlapping occipito-parietal white matter regions. We interpret our findings, based on current biophysical models, as a fingerprint of age-dependent brain atrophy and underlying microstructural changes in myelin, iron deposits and water. The VBQ approach we present allows for systematic unbiased exploration of the interaction between imaging parameters and extends current methods for detection of neurodegenerative processes in the brain. The demonstrated parameter-specific distribution patterns offer insights into age-related brain structure changes in vivo and provide essential baseline data for studying disease against a background of healthy ageing.

Diffusion tensor imaging shows structural remodeling of stroke mirror region: results from a pilot study.

Background: The role of the non-injured hemisphere in stroke recovery is poorly understood. In this pilot study, we sought to explore the presence of structural changes detectable by diffusion tensor imaging (DTI) in the contralesional hemispheres of patients who recovered well from ischemic stroke. Methods: We analyzed serial DTI data from 16 stroke patients who had moderate initial neurological deficits (NIHSS scores 3-12) and good functional outcome at 3-6 months (NIHSS score 0 or modified Rankin Score ≤1). We segmented the brain tissue in gray and white matter (GM and WM) and measured the apparent diffusion coefficient (ADC) and fractional anisotropy in the infarct, in the contralesional infarct mirror region as well as in concentrically expanding regions around them. Results: We found that GM and WM ADC significantly increased in the infarct region (p < 0.01) from acute to chronic time points, whereas in the infarct mirror region, GM and WM ADC increased (p < 0.01) and WM fractional anisotropy decreased (p < 0.05). No significant changes were detected in other regions. Conclusion: DTI-based metrics are sensitive to regional structural changes in the contralesional hemisphere during stroke recovery. Prospective studies in larger cohorts with varying levels of recovery are needed to confirm our findings.

Proposed gravitational wave observatory based on solid elastic spheres

Spherical gravitational wave (GW) detectors offer a wealth of so far unexplored possibilities to detect gravitational radiation. We find that a sphere can be used as a powerful testbed for any metric theory of gravity, not only general relativity as considered so far, by making use of a deconvolution procedure for all the electric components of the Riemann tensor. We also find that the spheres cross section is large at two frequencies, and advantageous at higher frequencies in the sense that a single antenna constitutes a real xylophone in its own. Proposed GW networks will greatly benefit from this. The main features of a two large sphere observatory are reported.

A connectome-based comparison of diffusion MRI schemes.

Diffusion MRI has evolved towards an important clinical diagnostic and research tool. Though clinical routine is using mainly diffusion weighted and tensor imaging approaches, Q-ball imaging and diffusion spectrum imaging techniques have become more widely available. They are frequently used in research-oriented investigations in particular those aiming at measuring brain network connectivity. In this work, we aim at assessing the dependency of connectivity measurements on various diffusion encoding schemes in combination with appropriate data modeling. We process and compare the structural connection matrices computed from several diffusion encoding schemes, including diffusion tensor imaging, q-ball imaging and high angular resolution schemes, such as diffusion spectrum imaging with a publically available processing pipeline for data reconstruction, tracking and visualization of diffusion MR imaging. The results indicate that the high angular resolution schemes maximize the number of obtained connections when applying identical processing strategies to the different diffusion schemes. Compared to the conventional diffusion tensor imaging, the added connectivity is mainly found for pathways in the 50-100mm range, corresponding to neighboring association fibers and long-range associative, striatal and commissural fiber pathways. The analysis of the major associative fiber tracts of the brain reveals striking differences between the applied diffusion schemes. More complex data modeling techniques (beyond tensor model) are recommended 1) if the tracts of interest run through large fiber crossings such as the centrum semi-ovale, or 2) if non-dominant fiber populations, e.g. the neighboring association fibers are the subject of investigation. An important finding of the study is that since the ground truth sensitivity and specificity is not known, the comparability between results arising from different strategies in data reconstruction and/or tracking becomes implausible to understand.

Quantitative morphometry analysis of the fetal brain using clinical MR imaging

In vivo fetal magnetic resonance imaging provides aunique approach for the study of early human braindevelopment [1]. In utero cerebral morphometry couldpotentially be used as a marker of the cerebralmaturation and help to distinguish between normal andabnormal development in ambiguous situations. However,this quantitative approach is a major challenge becauseof the movement of the fetus inside the amniotic cavity,the poor spatial resolution provided by very fast MRIsequences and the partial volume effect. Extensiveefforts are made to deal with the reconstruction ofhigh-resolution 3D fetal volumes based on severalacquisitions with lower resolution [2,3,4]. Frameworkswere developed for the segmentation of specific regionsof the fetal brain such as posterior fossa, brainstem orgerminal matrix [5,6], or for the entire brain tissue[7,8], applying the Expectation-Maximization MarkovRandom Field (EM-MRF) framework. However, many of theseprevious works focused on the young fetus (i.e. before 24weeks) and use anatomical atlas priors to segment thedifferent tissue or regions. As most of the gyraldevelopment takes place after the 24th week, acomprehensive and clinically meaningful study of thefetal brain should not dismiss the third trimester ofgestation. To cope with the rapidly changing appearanceof the developing brain, some authors proposed a dynamicatlas [8]. To our opinion, this approach however faces arisk of circularity: each brain will be analyzed /deformed using the template of its biological age,potentially biasing the effective developmental delay.Here, we expand our previous work [9] to proposepost-processing pipeline without prior that allow acomprehensive set of morphometric measurement devoted toclinical application. Data set & Methods: Prenatal MRimaging was performed with a 1-T system (GE MedicalSystems, Milwaukee) using single shot fast spin echo(ssFSE) sequences (TR 7000 ms, TE 180 ms, FOV 40 x 40 cm,slice thickness 5.4mm, in plane spatial resolution1.09mm). For each fetus, 6 axial volumes shifted by 1 mmwere acquired under motherâeuro?s sedation (about 1min pervolume). First, each volume is segmentedsemi-automatically using region-growing algorithms toextract fetal brain from surrounding maternal tissues.Inhomogeneity intensity correction [10] and linearintensity normalization are then performed. Brain tissues(CSF, GM and WM) are then segmented based on thelow-resolution volumes as presented in [9]. Ahigh-resolution image with isotropic voxel size of 1.09mm is created as proposed in [2] and using B-splines forthe scattered data interpolation [11]. Basal gangliasegmentation is performed using a levet setimplementation on the high-resolution volume [12]. Theresulting white matter image is then binarized and givenas an input in FreeSurfer software(http://surfer.nmr.mgh.harvard.edu) to providetopologically accurate three-dimensional reconstructionsof the fetal brain according to the local intensitygradient. References: [1] Guibaud, Prenatal Diagnosis29(4) (2009). [2] Rousseau, Acad. Rad. 13(9), 2006. [3]Jiang, IEEE TMI 2007. [4] Warfield IADB, MICCAI 2009. [5]Claude, IEEE Trans. Bio. Eng. 51(4) 2004. [6] Habas,MICCAI 2008. [7] Bertelsen, ISMRM 2009. [8] Habas,Neuroimage 53(2) 2010. [9] Bach Cuadra, IADB, MICCAI2009. [10] Styner, IEEE TMI 19(39 (2000). [11] Lee, IEEETrans. Visual. And Comp. Graph. 3(3), 1997. [12] BachCuadra, ISMRM 2010.