Automated delineation of white matter fiber tracts with a multiple region-of-interest approach

A Riemannian framework for incorporating white matter bundle prior in orientation distribution function based tractography algorithms

Diffusion magnetic resonance imaging (dMRI) tractography is a powerful approach to study brain structural connectivity. However, its reliability in a clinical context is still highly debated. Recent studies have shown that most classical algorithms achieve to recover the majority of existing true bundles. However, the generated tractograms contain many invalid bundles. This is due to the crossing fibers and bottleneck problems which increase the number of false positive fibers. In this work, we proposed to overpass this limitation with a novel method to guide the algorithms in those challenging regions with prior knowledge of the anatomy. We developed a method to create a combination of anatomical prior applicable to any orientation distribution function (ODF)-based tractography algorithms. The proposed method captures the tract orientation distribution (TOD) from an atlas of segmented fiber bundles and incorporates it during the tracking process, using a Riemannian framework. We tested the prior incorporation method on two ODF-based state-of-the-art algorithms, iFOD2 and Trekker PTT, on the diffusion-simulated connectivity (DiSCo) dataset and on the Human Connectome Project (HCP) data. We also compared our method with two bundles priors generated by the bundle specific tractography (BST) method. We showed that our method improves the overall spatial coverage and connectivity of a tractogram on the two datasets, especially in crossing fiber regions. Moreover, the fiber reconstruction may be improved on clinical data, informed by prior extracted on high quality data, and therefore could help in the study of brain anatomy and function.

Diffusion MRI with Machine Learning

Diffusion-weighted magnetic resonance imaging (dMRI) of the brain offers unique capabilities including noninvasive probing of tissue microstructure and structural connectivity. It is widely used for clinical assessment of disease and injury, and for neuroscience research. Analyzing the dMRI data to extract useful information for medical and scientific purposes can be challenging. The dMRI measurements may suffer from strong noise and artifacts, and may exhibit high inter-session and inter-scanner variability in the data, as well as inter-subject heterogeneity in brain structure. Moreover, the relationship between measurements and the phenomena of interest can be highly complex. Recent years have witnessed increasing use of machine learning methods for dMRI analysis. This manuscript aims to assess these efforts, with a focus on methods that have addressed data preprocessing and harmonization, microstructure mapping, tractography, and white matter tract analysis. We study the main findings, strengths, and weaknesses of the existing methods and suggest topics for future research. We find that machine learning may be exceptionally suited to tackle some of the difficult tasks in dMRI analysis. However, for this to happen, several shortcomings of existing methods and critical unresolved issues need to be addressed. There is a pressing need to improve evaluation practices, to increase the availability of rich training datasets and validation benchmarks, as well as model generalizability, reliability, and explainability concerns.

Comparison of T1-weighted landmark placement and ROI transfer onto diffusion-weighted EPI sequences for targeted tractography tasks in the optic nerve

Diffusion-based tractography in the optic nerve requires sampling strategies assisted by anatomical landmark information (regions of interest [ROIs]). We aimed to investigate the feasibility of expert-placed, high-resolution T1-weighted ROI-data transfer onto lower spatial resolution diffusion-weighted images. Slab volumes from 20 volunteers were acquired and preprocessed including distortion bias correction and artifact reduction. Constrained spherical deconvolution was used to generate a directional diffusion information grid (fibre orientation distribution-model [FOD]). Three neuroradiologists marked landmarks on both diffusion imaging variants and structural datasets. Structural ROI information (volumetric interpolated breath-hold sequence [VIBE]) was respectively registered (linear with 6/12 degrees of freedom [DOF]) onto single-shot EPI (ss-EPI) and readout-segmented EPI (rs-EPI) volumes, respectively. All eight ROI/FOD-combinations were compared in a targeted tractography task of the optic nerve pathway. Inter-rater reliability for placed ROIs among experts was highest in VIBE images (lower confidence interval 0.84 to 0.97, mean 0.91) and lower in both ss-EPI (0.61 to 0.95, mean 0.79) and rs-EPI (0.59 to 0.86, mean 0.70). Tractography success rate based on streamline selection performance was highest in VIBE-drawn ROIs registered (6-DOF) onto rs-EPI FOD (70.0% over 5%-threshold, capped to failed ratio 39/16) followed by both 12-DOF-registered (67.5%; 41/16) and nonregistered VIBE (67.5%; 40/23). On ss-EPI FOD, VIBE-ROI-datasets obtained fewer streamlines overall with each at 55.0% above 5%-threshold and with lower capped to failed ratio (6-DOF: 35/36; 12-DOF: 34/34, nonregistered 33/36). The combination of VIBE-placed ROIs (highest inter-rater reliability) with 6-DOF registration onto rs-EPI targets (best streamline selection performance) is most suitable for white matter template generation required in group studies.

Anatomically constrained tractography of the fetal brain

Diffusion-weighted Magnetic Resonance Imaging (dMRI) is increasingly used to study the fetal brain in utero. An important computation enabled by dMRI is streamline tractography, which has unique applications such as tract-specific analysis of the brain white matter and structural connectivity assessment. However, due to the low fetal dMRI data quality and the challenging nature of tractography, existing methods tend to produce highly inaccurate results. They generate many false streamlines while failing to reconstruct the streamlines that constitute the major white matter tracts. In this paper, we advocate for anatomically constrained tractography based on an accurate segmentation of the fetal brain tissue directly in the dMRI space. We develop a deep learning method to compute the segmentation automatically. Experiments on independent test data show that this method can accurately segment the fetal brain tissue and drastically improve the tractography results. It enables the reconstruction of highly curved tracts such as optic radiations. Importantly, our method infers the tissue segmentation and streamline propagation direction from a diffusion tensor fit to the dMRI data, making it applicable to routine fetal dMRI scans. The proposed method can facilitate the study of fetal brain white matter tracts with dMRI.

Population Atlas Analysis of Emerging Brain Structural Connections in the Human Fetus

BACKGROUND

A lack of in utero imaging data hampers our understanding of the connections in the human fetal brain. Generalizing observations from postmortem subjects and premature newborns is inaccurate due to technical and biological differences.

PURPOSE

To evaluate changes in fetal brain structural connectivity between 23 and 35 weeks postconceptional age using a spatiotemporal atlas of diffusion tensor imaging (DTI).

STUDY TYPE

Retrospective.

POPULATION

Publicly available diffusion atlases, based on 60 healthy women (age 18-45 years) with normal prenatal care, from 23 and 35 weeks of gestation.

FIELD STRENGTH/SEQUENCE

3.0 Tesla/DTI acquired with diffusion-weighted echo planar imaging (EPI).

ASSESSMENT

We performed whole-brain fiber tractography from DTI images. The cortical plate of each diffusion atlas was segmented and parcellated into 78 regions derived from the Edinburgh Neonatal Atlas (ENA33). Connectivity matrices were computed, representing normalized fiber connections between nodes. We examined the relationship between global efficiency (GE), local efficiency (LE), small-worldness (SW), nodal efficiency (NE), and betweenness centrality (BC) with gestational age (GA) and with laterality.

STATISTICAL TESTS

Linear regression was used to analyze changes in GE, LE, NE, and BC throughout gestation, and to assess changes in laterality. The t-tests were used to assess SW. P-values were corrected using Holm-Bonferroni method. A corrected P-value <0.05 was considered statistically significant.

RESULTS

Network analysis revealed a significant weekly increase in GE (5.83%/week, 95% CI 4.32-7.37), LE (5.43%/week, 95% CI 3.63-7.25), and presence of SW across GA. No significant hemisphere differences were found in GE (P = 0.971) or LE (P = 0.458). Increasing GA was significantly associated with increasing NE in 41 nodes, increasing BC in 3 nodes, and decreasing BC in 2 nodes.

DATA CONCLUSION

Extensive network development and refinement occur in the second and third trimesters, marked by a rapid increase in global integration and local segregation.

LEVEL OF EVIDENCE

3 TECHNICAL EFFICACY: Stage 2.

Serial neuroimaging of brain growth and development in very preterm infants receiving tailored neuropromotive support in the NICU. Protocol for a prospective cohort study

Introduction

Children born very preterm (VP) remain at risk for long-term neurodevelopmental impairment. Patterns of brain growth and injury, and how early neuropromotive therapies might mitigate developmental risk in VP infants remain insufficiently understood.

Methods

This is a prospective cohort study of VP infants born at/before 32 weeks gestation. The study will enroll n = 75 consecutively-born VP infants in a level-III NICU. Exposed infants will be categorized into two groups (group 1: low-risk, n = 25 or group 2: high-risk, n = 25) based on the degree of neurological injury on early brain magnetic resonance imaging (MRI) at enrollment. Infants in the low-risk group (i.e., without significant injury defined as intraventricular hemorrhage with dilation, moderate or severe white matter injury, or cerebellar hemorrhage) will receive neurodevelopmental support utilizing the Supporting and Enhancing NICU Sensory Experiences (SENSE) program, while infants in the high-risk group (with neurological injury) will receive more intensive neurorehabilitative support (SENSE-plus). Age-specific, tailored sensory experiences will be facilitated contingently, preferentially by the infant's family with coaching from NICU staff. VP infants in exposure groups will undergo a brain MRI approximately every 2 weeks from enrollment until term-equivalent to monitor brain growth and evolution of injury. Exposed infants will be compared with a reference group (group 3: n = 25), i.e. VP infants whose families decline initial enrollment in SENSE, and subsequently undergo a term-equivalent brain MRI for other purposes. The primary aim of this study is characterization of term-equivalent brain growth and development among VP infants receiving NICU-based neuropromotive interventions compared to VP infants receiving the standard of care. Secondary aims include defining the timing and factors associated with total and regional brain growth on serial brain MRI among VP infants, (Aim 2), and using early imaging to tailor developmental intervention in the NICU while exploring associations with outcomes in VP infants at discharge and at two years corrected age (Aim 3).

Discussion

This study will address gaps in understanding patterns of brain growth and injury drawing on serial MRI of hospitalized VP infants. These data will also explore the impact of intensive, tailored neuropromotive support delivered prior to term-equivalent on child and family outcomes.

TBSS++: A novel computational method for Tract-Based Spatial Statistics

Diffusion-weighted magnetic resonance imaging (dMRI) is widely used to assess the brain white matter. One of the most common computations in dMRI involves cross-subject tract-specific analysis, whereby dMRI-derived biomarkers are compared between cohorts of subjects. The accuracy and reliability of these studies hinges on the ability to compare precisely the same white matter tracts across subjects. This is an intricate and error-prone computation. Existing computational methods such as Tract-Based Spatial Statistics (TBSS) suffer from a host of shortcomings and limitations that can seriously undermine the validity of the results. We present a new computational framework that overcomes the limitations of existing methods via (i) accurate segmentation of the tracts, and (ii) precise registration of data from different subjects/scans. The registration is based on fiber orientation distributions. To further improve the alignment of cross-subject data, we create detailed atlases of white matter tracts. These atlases serve as an unbiased reference space where the data from all subjects is registered for comparison. Extensive evaluations show that, compared with TBSS, our proposed framework offers significantly higher reproducibility and robustness to data perturbations. Our method promises a drastic improvement in accuracy and reproducibility of cross-subject dMRI studies that are routinely used in neuroscience and medical research.

Direct segmentation of brain white matter tracts in diffusion MRI

The brain white matter consists of a set of tracts that connect distinct regions of the brain. Segmentation of these tracts is often needed for clinical and research studies. Diffusion-weighted MRI offers unique contrast to delineate these tracts. However, existing segmentation methods rely on intermediate computations such as tractography or estimation of fiber orientation density. These intermediate computations, in turn, entail complex computations that can result in unnecessary errors. Moreover, these intermediate computations often require dense multi-shell measurements that are unavailable in many clinical and research applications. As a result, current methods suffer from low accuracy and poor generalizability. Here, we propose a new deep learning method that segments these tracts directly from the diffusion MRI data, thereby sidestepping the intermediate computation errors. Our experiments show that this method can achieve segmentation accuracy that is on par with the state of the art methods (mean Dice Similarity Coefficient of 0.826). Compared with the state of the art, our method offers far superior generalizability to undersampled data that are typical of clinical studies and to data obtained with different acquisition protocols. Moreover, we propose a new method for detecting inaccurate segmentations and show that it is more accurate than standard methods that are based on estimation uncertainty quantification. The new methods can serve many critically important clinical and scientific applications that require accurate and reliable non-invasive segmentation of white matter tracts.

Associations of body composition with regional brain volumes and white matter microstructure in very preterm infants

OBJECTIVE

To determine associations between body composition and concurrent measures of brain development including (1) Tissue-specific brain volumes and (2) White matter microstructure, among very preterm infants at term equivalent age.

DESIGN

Prospective observational study.

SETTING

Single-centre academic level III neonatal intensive care unit.

PATIENTS

We studied 85 infants born <33 weeks' gestation.

METHODS

At term equivalent age, infants underwent air displacement plethysmography to determine body composition, and brain MRI from which we quantified tissue-specific brain volumes and fractional anisotropy (FA) of white matter tracts. We estimated associations of fat and lean mass Z-scores with each brain outcome, using linear mixed models adjusted for intrafamilial correlation among twins and potential confounding variables.

RESULTS

Median gestational age was 29 weeks (range 23.4-32.9). One unit greater lean mass Z-score was associated with larger total brain volume (10.5 cc, 95% CI 3.8 to 17.2); larger volumes of the cerebellum (1.2 cc, 95% CI 0.5 to 1.9) and white matter (4.5 cc, 95% CI 0.7 to 8.3); and greater FA in the left cingulum (0.3%, 95% CI 0.1% to 0.6%), right uncinate fasciculus (0.2%, 95% CI 0.0% to 0.5%), and right posterior limb of the internal capsule (0.3%, 95% CI 0.03% to 0.6%). Fat Z-scores were not associated with any outcome.

CONCLUSIONS

Lean mass-but not fat-at term was associated with larger brain volume and white matter microstructure differences that suggest improved maturation. Lean mass accrual may index brain growth and development.

Associations of Macronutrient Intake Determined by Point-of-Care Human Milk Analysis with Brain Development among very Preterm Infants

Point-of-care human milk analysis is now feasible in the neonatal intensive care unit (NICU) and allows accurate measurement of macronutrient delivery. Higher macronutrient intakes over this period may promote brain growth and development. In a prospective, observational study of 55 infants born at <32 weeks’ gestation, we used a mid-infrared spectroscopy-based human milk analyzer to measure the macronutrient content in repeated samples of human milk over the NICU hospitalization. We calculated daily nutrient intakes from unfortified milk and assigned infants to quintiles based on median intakes over the hospitalization. Infants underwent brain magnetic resonance imaging at term equivalent age to quantify total and regional brain volumes and fractional anisotropy of white matter tracts. Infants in the highest quintile of energy intake from milk, as compared with the lower four quintiles, had larger total brain volume (31 cc, 95% confidence interval [CI]: 5, 56), cortical gray matter (15 cc, 95%CI: 1, 30), and white matter volume (23 cc, 95%CI: 12, 33). Higher protein intake was associated with larger total brain (36 cc, 95%CI: 7, 65), cortical gray matter (22 cc, 95%CI: 6, 38) and deep gray matter (1 cc, 95%CI: 0.1, 3) volumes. These findings suggest innovative strategies to close nutrient delivery gaps in the NICU may promote brain growth for preterm infants.

A Cross-Sectional Study on the Impact of Arterial Stiffness on the Corpus Callosum, a Key White Matter Tract Implicated in Alzheimer's Disease

BACKGROUND

Vascular risk factors such as arterial stiffness play an important role in the etiology of Alzheimer's disease (AD), presumably due to the emergence of white matter lesions. However, the impact of arterial stiffness to white matter structure involved in the etiology of AD, including the corpus callosum remains poorly understood.

OBJECTIVE

The aims of the study are to better understand the relationship between arterial stiffness, white matter microstructure, and perfusion of the corpus callosum in older adults.

METHODS

Arterial stiffness was estimated using the gold standard measure of carotid-femoral pulse wave velocity (cfPWV). Cognitive performance was evaluated with the Trail Making Test part B-A. Neurite orientation dispersion and density imaging was used to obtain microstructural information such as neurite density and extracellular water diffusion. The cerebral blood flow was estimated using arterial spin labelling.

RESULTS

cfPWV better predicts the microstructural integrity of the corpus callosum when compared with other index of vascular aging (the augmentation index, the systolic blood pressure, and the pulse pressure). In particular, significant associations were found between the cfPWV, an alteration of the extracellular water diffusion, and a neuronal density increase in the body of the corpus callosum which was also correlated with the performance in cognitive flexibility.

CONCLUSION

Our results suggest that arterial stiffness is associated with an alteration of brain integrity which impacts cognitive function in older adults.

Cortically constrained shape recognition: Automated white matter tract segmentation validated in the pediatric brain

BACKGROUND AND PURPOSE

Manual segmentation of white matter (WM) bundles requires extensive training and is prohibitively labor-intensive for large-scale studies. Automated segmentation methods are necessary in order to eliminate operator dependency and to enable reproducible studies. Significant changes in the WM landscape throughout childhood require flexible methods to capture the variance across the span of brain development.

METHODS

Here, we describe a novel automated segmentation tool called Cortically Constrained Shape Recognition (CCSR), which combines two complementary approaches: (1) anatomical connectivity priors based on FreeSurfer-derived regions of interest and (2) shape priors based on 3-dimensional streamline bundle atlases applied using RecoBundles. We tested the performance and repeatability of this approach by comparing volume and diffusion metrics of the main language WM tracts that were both manually and automatically segmented in a pediatric cohort acquired at the UCSF Dyslexia Center (n = 59; 25 females; average age: 11 ± 2; range: 7-14).

RESULTS

The CCSR approach showed high agreement with the expert manual segmentations: across all tracts, the spatial overlap between tract volumes showed an average Dice Similarity Coefficient (DSC) of 0.76, and the fractional anisotropy (FA) on average had a Lin's Concordance Correlation Coefficient (CCC) of 0.81. The CCSR's repeatability in a subset of this cohort achieved a DSC of 0.92 on average across all tracts.

CONCLUSION

This novel automated segmentation approach is a promising tool for reproducible large-scale tractography analyses in pediatric populations and particularly for the quantitative assessment of structural connections underlying various clinical presentations in neurodevelopmental disorders.

DTI Atlases Evaluations

The cerebral atlas of diffusion tensor magnetic resonance image (DT-MRI, shorted as DTI) is one of the key issues in neuroimaging research. It is crucial for comparisons of neuronal structural integrity and connectivity across populations. Usually, the atlas is constructed by iteratively averaging the registered individual image. In tradition, the fuzzy group average image is easily generated in the initial stage, which is harmful to providing clear guidance for subsequent registration, to the performance of the final atlas. To solve this problem, an improved unbiased DTI atlas construction algorithm based on adaptive weights is proposed in this paper. The adaptive weighted strategy based on diffeomorphic deformable tensor registration is introduced. At the same time, the distance measure for tensors is used as a constraint condition, which ensures the unbiasedness of the atlas. Then, using 77 DTIs from the dataset in http://www.brain-development.org , three study-specific atlases, i.e. the constructed atlases of the proposed algorithm and two open-sourced algorithms (DTIAtlasBuilder and DTI-TK), are compared with two standardized atlases (IIT v. 4.1 and NTU-DSI-122-DTI). The performances of the atlases were evaluated in spatial normalization way with six region-based criteria (including Euclidean distances between diffusion tensors, Euclidean distances of the deviatoric tensors, standard deviation, overlaps of eigenvalue-eigenvector, cross-correlations and three sets angles of eigenvalue-eigenvector pairs between diffusion tensors) and three fiber-based criteria (including distances between fiber bundles, angles between fiber bundles and fiber property profile-based criteria). The experimental results showed that the overall performances of the study-specific atlases are better than those of the standardized atlases for specific datasets, and the comprehensive performance of the improved algorithm proposed in this paper is the best.

Interpolation and Averaging of Diffusion MRI Multi-Compartment Models

Multi-compartment models (MCM) are increasingly used to characterize the brain white matter microstructure from diffusion-weighted imaging (DWI). Their use in clinical studies is however limited by the inability to resample an MCM image towards a common reference frame, or to construct atlases from such brain microstructure models. We propose to solve this problem by first identifying that these two tasks amount to the same problem. We propose to tackle it by viewing it as a simplification problem, solved thanks to spectral clustering and the definition of semi-metrics between several usual compartments encountered in the MCM literature. This generic framework is evaluated for two models: the multi-tensor model where individual fibers are modeled as individual tensors and the diffusion direction imaging (DDI) model that differentiates intra- and extra-axonal components of each fiber. Results on simulated data, simulated transformations and real data show the ability of our method to well interpolate MCM images of these types. We finally present as an application an MCM template of normal controls constructed using our approach.

Multiple sclerosis lesions in motor tracts from brain to cervical cord: spatial distribution and correlation with disability

Despite important efforts to solve the clinico-radiological paradox, correlation between lesion load and physical disability in patients with multiple sclerosis remains modest. One hypothesis could be that lesion location in corticospinal tracts plays a key role in explaining motor impairment. In this study, we describe the distribution of lesions along the corticospinal tracts from the cortex to the cervical spinal cord in patients with various disease phenotypes and disability status. We also assess the link between lesion load and location within corticospinal tracts, and disability at baseline and 2-year follow-up. We retrospectively included 290 patients (22 clinically isolated syndrome, 198 relapsing remitting, 39 secondary progressive, 31 primary progressive multiple sclerosis) from eight sites. Lesions were segmented on both brain (T2-FLAIR or T2-weighted) and cervical (axial T2- or T2*-weighted) MRI scans. Data were processed using an automated and publicly available pipeline. Brain, brainstem and spinal cord portions of the corticospinal tracts were identified using probabilistic atlases to measure the lesion volume fraction. Lesion frequency maps were produced for each phenotype and disability scores assessed with Expanded Disability Status Scale score and pyramidal functional system score. Results show that lesions were not homogeneously distributed along the corticospinal tracts, with the highest lesion frequency in the corona radiata and between C2 and C4 vertebral levels. The lesion volume fraction in the corticospinal tracts was higher in secondary and primary progressive patients (mean = 3.6 ± 2.7% and 2.9 ± 2.4%), compared to relapsing-remitting patients (1.6 ± 2.1%, both P < 0.0001). Voxel-wise analyses confirmed that lesion frequency was higher in progressive compared to relapsing-remitting patients, with significant bilateral clusters in the spinal cord corticospinal tracts (P < 0.01). The baseline Expanded Disability Status Scale score was associated with lesion volume fraction within the brain (r = 0.31, P < 0.0001), brainstem (r = 0.45, P < 0.0001) and spinal cord (r = 0.57, P < 0.0001) corticospinal tracts. The spinal cord corticospinal tracts lesion volume fraction remained the strongest factor in the multiple linear regression model, independently from cord atrophy. Baseline spinal cord corticospinal tracts lesion volume fraction was also associated with disability progression at 2-year follow-up (P = 0.003). Our results suggest a cumulative effect of lesions within the corticospinal tracts along the brain, brainstem and spinal cord portions to explain physical disability in multiple sclerosis patients, with a predominant impact of intramedullary lesions.

Diffusion Tensor Imaging Abnormalities in the Uncinate Fasciculus and Inferior Longitudinal Fasciculus in Phelan-McDermid Syndrome

BACKGROUND

This cohort study utilized diffusion tensor imaging tractography to compare the uncinate fasciculus and inferior longitudinal fasciculus in children with Phelan-McDermid syndrome with age-matched controls and investigated trends between autism spectrum diagnosis and the integrity of the uncinate fasciculus and inferior longitudinal fasciculus white matter tracts.

METHODS

This research was conducted under a longitudinal study that aims to map the genotype, phenotype, and natural history of Phelan-McDermid syndrome and identify biomarkers using neuroimaging (ClinicalTrial NCT02461420). Patients were aged three to 21 years and underwent longitudinal neuropsychologic assessment over 24 months. MRI processing and analyses were completed using previously validated image analysis software distributed as the Computational Radiology Kit (http://crl.med.harvard.edu/). Whole-brain connectivity was generated for each subject using a stochastic streamline tractography algorithm, and automatically defined regions of interest were used to map the uncinate fasciculus and inferior longitudinal fasciculus.

RESULTS

There were 10 participants (50% male; mean age 11.17 years) with Phelan-McDermid syndrome (n = 8 with autism). Age-matched controls, enrolled in a separate longitudinal study (NIH R01 NS079788), underwent the same neuroimaging protocol. There was a statistically significant decrease in the uncinate fasciculus fractional anisotropy measure and a statistically significant increase in uncinate fasciculus mean diffusivity measure, in the patient group versus controls in both right and left tracts (P ≤ 0.024).

CONCLUSION

Because the uncinate fasciculus plays a critical role in social and emotional interaction, this tract may underlie some deficits seen in the Phelan-McDermid syndrome population. These findings need to be replicated in a larger cohort.

Early white matter development is abnormal in tuberous sclerosis complex patients who develop autism spectrum disorder

Background

Autism spectrum disorder (ASD) is prevalent in tuberous sclerosis complex (TSC), occurring in approximately 50% of patients, and is hypothesized to be caused by disruption of neural circuits early in life. Tubers, or benign hamartomas distributed stochastically throughout the brain, are the most conspicuous of TSC neuropathology, but have not been consistently associated with ASD. Widespread neuropathology of the white matter, including deficits in myelination, neuronal migration, and axon formation, exist and may underlie ASD in TSC. We sought to identify the neural circuits associated with ASD in TSC by identifying white matter microstructural deficits in a prospectively recruited, longitudinally studied cohort of TSC infants.

Methods

TSC infants were recruited within their first year of life and longitudinally imaged at time of recruitment, 12 months of age, and at 24 months of age. Autism was diagnosed at 24 months of age with the ADOS-2. There were 108 subjects (62 TSC-ASD, 55% male; 46 TSC+ASD, 52% male) with at least one MRI and a 24-month ADOS, for a total of 187 MRI scans analyzed (109 TSC-ASD; 78 TSC+ASD). Diffusion tensor imaging properties of multiple white matter fiber bundles were sampled using a region of interest approach. Linear mixed effects modeling was performed to test the hypothesis that infants who develop ASD exhibit poor white matter microstructural integrity over the first 2 years of life compared to those who do not develop ASD.

Results

Subjects with TSC and ASD exhibited reduced fractional anisotropy in 9 of 17 white matter regions, sampled from the arcuate fasciculus, cingulum, corpus callosum, anterior limbs of the internal capsule, and the sagittal stratum, over the first 2 years of life compared to TSC subjects without ASD. Mean diffusivity trajectories did not differ between groups.

Conclusions

Underconnectivity across multiple white matter fiber bundles develops over the first 2 years of life in subjects with TSC and ASD. Future studies examining brain-behavior relationships are needed to determine how variation in the brain structure is associated with ASD symptoms.

Reproducibility of Structural and Diffusion Tensor Imaging in the TACERN Multi-Center Study

BACKGROUND

Multi-site MRI studies are often necessary for recruiting sufficiently sized samples when studying rare conditions. However, they require pooling data from multiple scanners into a single data set, and therefore it is critical to evaluate the variability of quantitative MRI measures within and across scanners used in multi-site studies. The aim of this study was to evaluate the reproducibility of structural and diffusion weighted (DW) MRI measurements acquired on seven scanners at five medical centers as part of the Tuberous Sclerosis Complex Autism Center of Excellence Research Network (TACERN) multisite study.

METHODS

The American College of Radiology (ACR) phantom was imaged monthly to measure reproducibility of signal intensity and uniformity within and across seven 3T scanners from General Electric, Philips, and Siemens vendors. One healthy adult male volunteer was imaged repeatedly on all seven scanners under the TACERN structural and DW protocol (5 b = 0 s/mm2 and 30 b = 1000 s/mm2) over a period of 5 years (age 22-27 years). Reproducibility of inter- and intra-scanner brain segmentation volumes and diffusion tensor imaging metrics fractional anisotropy (FA) and mean diffusivity (MD) within white matter regions was quantified with coefficient of variation.

RESULTS

The American College of Radiology Phantom signal intensity and uniformity were similar across scanners and changed little over time, with a mean intra-scanner coefficient of variation of 3.6 and 1.8%, respectively. The mean inter- and intra-scanner coefficients of variation of brain structure volumes derived from T1-weighted (T1w) images of the human phantom were 3.3 and 1.1%, respectively. The mean inter- and intra-scanner coefficients of variation of FA in white matter regions were 4.5 and 2.5%, while the mean inter- and intra-scanner coefficients of variation of MD in white matter regions were 5.4 and 1.5%.

CONCLUSION

Our results suggest that volumetric and diffusion tensor imaging (DTI) measurements are highly reproducible between and within scanners and provide typical variation amplitudes that can be used as references to interpret future findings in the TACERN network.

An anatomically curated fiber clustering white matter atlas for consistent white matter tract parcellation across the lifespan

This work presents an anatomically curated white matter atlas to enable consistent white matter tract parcellation across different populations. Leveraging a well-established computational pipeline for fiber clustering, we create a tract-based white matter atlas including information from 100 subjects. A novel anatomical annotation method is proposed that leverages population-based brain anatomical information and expert neuroanatomical knowledge to annotate and categorize the fiber clusters. A total of 256 white matter structures are annotated in the proposed atlas, which provides one of the most comprehensive tract-based white matter atlases covering the entire brain to date. These structures are composed of 58 deep white matter tracts including major long range association and projection tracts, commissural tracts, and tracts related to the brainstem and cerebellar connections, plus 198 short and medium range superficial fiber clusters organized into 16 categories according to the brain lobes they connect. Potential false positive connections are annotated in the atlas to enable their exclusion from analysis or visualization. In addition, the proposed atlas allows for a whole brain white matter parcellation into 800 fiber clusters to enable whole brain connectivity analyses. The atlas and related computational tools are open-source and publicly available. We evaluate the proposed atlas using a testing dataset of 584 diffusion MRI scans from multiple independently acquired populations, across genders, the lifespan (1 day-82 years), and different health conditions (healthy control, neuropsychiatric disorders, and brain tumor patients). Experimental results show successful white matter parcellation across subjects from different populations acquired on multiple scanners, irrespective of age, gender or disease indications. Over 99% of the fiber tracts annotated in the atlas were detected in all subjects on average. One advantage in terms of robustness is that the tract-based pipeline does not require any cortical or subcortical segmentations, which can have limited success in young children and patients with brain tumors or other structural lesions. We believe this is the first demonstration of consistent automated white matter tract parcellation across the full lifespan from birth to advanced age.

Corpus Callosum White Matter Diffusivity Reflects Cumulative Neurological Comorbidity in Tuberous Sclerosis Complex

INTRODUCTION

Neurological manifestations in Tuberous Sclerosis Complex (TSC) are highly variable. Diffusion tensor imaging (DTI) may reflect the neurological disease burden. We analyzed the association of autism spectrum disorder (ASD), intellectual disability (ID) and epilepsy with callosal DTI metrics in subjects with and without TSC.

METHODS

186 children underwent 3T MRI DTI: 51 with TSC (19 with concurrent ASD), 46 with non-syndromic ASD and 89 healthy controls (HC). Subgroups were based on presence of TSC, ASD, ID, and epilepsy. Density-weighted DTI metrics obtained from tractography of the corpus callosum were fitted using a 2-parameter growth model. We estimated distributions using bootstrapping and calculated half-life and asymptote of the fitted curves.

RESULTS

TSC was associated with a lower callosal fractional anisotropy (FA) than ASD, and ASD with a lower FA than HC. ID, epilepsy and ASD diagnosis were each associated with lower FA values, demonstrating additive effects. In TSC, the largest change in FA was related to a comorbid diagnosis of ASD. Mean diffusivity (MD) showed an inverse relationship to FA. Some subgroups were too small for reliable data fitting.

CONCLUSIONS

Using a cross-disorder approach, this study demonstrates cumulative abnormality of callosal white matter diffusion with increasing neurological comorbidity.

Groupwise Multichannel Image Registration

Multichannel image registration is an important challenge in medical image analysis. Multichannel images result from modalities such as dual-energy CT or multispectral microscopy. Besides, multichannel feature images can be derived from acquired images, for instance, by applying multiscale feature banks to the original images to register. Multichannel registration techniques have been proposed, but most of them are applicable to only two multichannel images at a time. In the present study, we propose to formulate multichannel registration as a groupwise image registration problem. In this way, we derive a method that allows the registration of two or more multichannel images in a fully symmetric manner (i.e., all images play the same role in the registration procedure), and therefore, has transitive consistency by definition. The method that we introduce is applicable to any number of multichannel images, any number of channels per image, and it allows to take into account correlation between any pair of images and not just corresponding channels. In addition, it is fully modular in terms of dissimilarity measure, transformation model, regularisation method, and optimisation strategy. For two multimodal datasets, we computed feature images from the initially acquired images, and applied the proposed registration technique to the newly created sets of multichannel images. MIND descriptors were used as feature images, and we chose total correlation as groupwise dissimilarity measure. Results show that groupwise multichannel image registration is a competitive alternative to the pairwise multichannel scheme, in terms of registration accuracy and insensitivity towards registration reference spaces.

A comparison of three fiber tract delineation methods and their impact on white matter analysis

Diffusion magnetic resonance imaging (dMRI) is an important method for studying white matter connectivity in the brain in vivo in both healthy and clinical populations. Improvements in dMRI tractography algorithms, which reconstruct macroscopic three-dimensional white matter fiber pathways, have allowed for methodological advances in the study of white matter; however, insufficient attention has been paid to comparing post-tractography methods that extract white matter fiber tracts of interest from whole-brain tractography. Here we conduct a comparison of three representative and conceptually distinct approaches to fiber tract delineation: 1) a manual multiple region of interest-based approach, 2) an atlas-based approach, and 3) a groupwise fiber clustering approach, by employing methods that exemplify these approaches to delineate the arcuate fasciculus, the middle longitudinal fasciculus, and the uncinate fasciculus in 10 healthy male subjects. We enable qualitative comparisons across methods, conduct quantitative evaluations of tract volume, tract length, mean fractional anisotropy, and true positive and true negative rates, and report measures of intra-method and inter-method agreement. We discuss methodological similarities and differences between the three approaches and the major advantages and drawbacks of each, and review research and clinical contexts for which each method may be most apposite. Emphasis is given to the means by which different white matter fiber tract delineation approaches may systematically produce variable results, despite utilizing the same input tractography and reliance on similar anatomical knowledge.

Neighborhood resolved fiber orientation distributions (NRFOD) in automatic labeling of white matter fiber pathways

Accurate digital representation of major white matter bundles in the brain is an important goal in neuroscience image computing since the representations can be used for surgical planning, intra-patient longitudinal analysis and inter-subject population connectivity studies. Reconstructing desired fiber bundles generally involves manual selection of regions of interest by an expert, which is subject to user bias and fatigue, hence an automation is desirable. To that end, we first present a novel anatomical representation based on Neighborhood Resolved Fiber Orientation Distributions (NRFOD) along the fibers. The resolved fiber orientations are obtained by generalized q-sampling imaging (GQI) and a subsequent diffusion decomposition method. A fiber-to-fiber distance measure between the proposed fiber representations is then used in a density-based clustering framework to select the clusters corresponding to the major pathways of interest. In addition, neuroanatomical priors are utilized to constrain the set of candidate fibers before density-based clustering. The proposed fiber clustering approach is exemplified on automation of the reconstruction of the major fiber pathways in the brainstem: corticospinal tract (CST); medial lemniscus (ML); middle cerebellar peduncle (MCP); inferior cerebellar peduncle (ICP); superior cerebellar peduncle (SCP). Experimental results on Human Connectome Project (HCP)'s publicly available "WU-Minn 500 Subjects + MEG2 dataset" and expert evaluations demonstrate the potential of the proposed fiber clustering method in brainstem white matter structure analysis.

Evaluation of standardized and study-specific diffusion tensor imaging templates of the adult human brain: Template characteristics, spatial normalization accuracy, and detection of small inter-group FA differences

Digital diffusion tensor imaging (DTI) templates of the adult human brain are commonly used in neuroimaging research, and their characteristics influence the accuracy of the application. However, a systematic evaluation of the characteristics and performance of standardized and study-specific DTI templates has not been conducted. The purpose of this work was to compare eight available standardized DTI templates to each other (ICBM81, ENIGMA, FMRIB58, SRI24, IIT2, NTU-DSI-122-DTI, IIT v.3.0, Eve), as well as to study-specific templates, in terms of template characteristics (image sharpness, ability to identify small brain structures, artifacts, mean values, noise properties) and performance in spatial normalization and detection of small inter-group FA differences. The IIT v.3.0 template was shown to combine a number of desirable characteristics: includes full-tensor information, is population-based, has high image sharpness, shows no visible artifacts, has low noise levels, has diffusion tensor properties and spatial features representative of data from the average individual adult brain. Furthermore, the IIT v.3.0 template was shown to allow higher inter-subject DTI spatial normalization accuracy, and detection of smaller inter-group FA differences, compared to all other templates, including study-specific templates. These findings were consistent when evaluating the templates in younger as well as older adult cohorts.

A tract-specific approach to assessing white matter in preterm infants

Diffusion-weighted imaging (DWI) is becoming an increasingly important tool for studying brain development. DWI analyses relying on manually-drawn regions of interest and tractography using manually-placed waypoints are considered to provide the most accurate characterisation of the underlying brain structure. However, these methods are labour-intensive and become impractical for studies with large cohorts and numerous white matter (WM) tracts. Tract-specific analysis (TSA) is an alternative WM analysis method applicable to large-scale studies that offers potential benefits. TSA produces a skeleton representation of WM tracts and projects the group's diffusion data onto the skeleton for statistical analysis. In this work we evaluate the performance of TSA in analysing preterm infant data against results obtained from native space tractography and tract-based spatial statistics. We evaluate TSA's registration accuracy of WM tracts and assess the agreement between native space data and template space data projected onto WM skeletons, in 12 tracts across 48 preterm neonates. We show that TSA registration provides better WM tract alignment than a previous protocol optimised for neonatal spatial normalisation, and that TSA projects FA values that match well with values derived from native space tractography. We apply TSA for the first time to a preterm neonatal population to study the effects of age at scan on WM tracts around term equivalent age. We demonstrate the effects of age at scan on DTI metrics in commissural, projection and association fibres. We demonstrate the potential of TSA for WM analysis and its suitability for infant studies involving multiple tracts.

•

Evaluation of tract-specific analysis (TSA) for white matter studies in infants.

•

TSA improves white matter tract alignment over scalar-based registration.

•

TSA closely approximates native space tractography DTI values.

•

The first application of TSA to a neonatal population.

Segmentation of the Cerebellar Peduncles Using a Random Forest Classifier and a Multi-object Geometric Deformable Model: Application to Spinocerebellar Ataxia Type 6

The cerebellar peduncles, comprising the superior cerebellar peduncles (SCPs), the middle cerebellar peduncle (MCP), and the inferior cerebellar peduncles (ICPs), are white matter tracts that connect the cerebellum to other parts of the central nervous system. Methods for automatic segmentation and quantification of the cerebellar peduncles are needed for objectively and efficiently studying their structure and function. Diffusion tensor imaging (DTI) provides key information to support this goal, but it remains challenging because the tensors change dramatically in the decussation of the SCPs (dSCP), the region where the SCPs cross. This paper presents an automatic method for segmenting the cerebellar peduncles, including the dSCP. The method uses volumetric segmentation concepts based on extracted DTI features. The dSCP and noncrossing portions of the peduncles are modeled as separate objects, and are initially classified using a random forest classifier together with the DTI features. To obtain geometrically correct results, a multi-object geometric deformable model is used to refine the random forest classification. The method was evaluated using a leave-one-out cross-validation on five control subjects and four patients with spinocerebellar ataxia type 6 (SCA6). It was then used to evaluate group differences in the peduncles in a population of 32 controls and 11 SCA6 patients. In the SCA6 group, we have observed significant decreases in the volumes of the dSCP and the ICPs and significant increases in the mean diffusivity in the noncrossing SCPs, the MCP, and the ICPs. These results are consistent with a degeneration of the cerebellar peduncles in SCA6 patients.

Automatic whole brain tract-based analysis using predefined tracts in a diffusion spectrum imaging template and an accurate registration strategy

Automated tract-based analysis of diffusion MRI is an important tool for investigating tract integrity of the cerebral white matter. Current template-based automatic analyses still lack a comprehensive list of tract atlas and an accurate registration method. In this study, tract-based automatic analysis (TBAA) was developed to meet the demands. Seventy-six major white matter tracts were reconstructed on a high-quality diffusion spectrum imaging (DSI) template, and an advanced two-step registration strategy was proposed by incorporating anatomical information of the gray matter from T1-weighted images in addition to microstructural information of the white matter from diffusion-weighted images. The automatic analysis was achieved by establishing a transformation between the DSI template and DSI dataset of the subject derived from the registration strategy. The tract coordinates in the template were transformed to native space in the individual's DSI dataset, and the microstructural properties of major tract bundles were sampled stepwise along the tract coordinates of the subject's DSI dataset. In a validation study of eight well-known tracts, our results showed that TBAA had high geometric agreement with manual tracts in both deep and superficial parts but significantly smaller measurement variability than manual method in functional difference. Additionally, the feasibility of the method was demonstrated by showing tracts with altered microstructural properties in patients with schizophrenia. Fifteen major tract bundles were found to have significant differences after controlling the family-wise error rate. In conclusion, the proposed TBAA method is potentially useful in brain-wise investigations of white matter tracts, particularly for a large cohort study.

A Model of Population and Subject (MOPS) Intensities With Application to Multiple Sclerosis Lesion Segmentation

White matter (WM) lesions are thought to play an important role in multiple sclerosis (MS) disease burden. Recent work in the automated segmentation of white matter lesions from magnetic resonance imaging has utilized a model in which lesions are outliers in the distribution of tissue signal intensities across the entire brain of each patient. However, the sensitivity and specificity of lesion detection and segmentation with these approaches have been inadequate. In our analysis, we determined this is due to the substantial overlap between the whole brain signal intensity distribution of lesions and normal tissue. Inspired by the ability of experts to detect lesions based on their local signal intensity characteristics, we propose a new algorithm that achieves lesion and brain tissue segmentation through simultaneous estimation of a spatially global within-the-subject intensity distribution and a spatially local intensity distribution derived from a healthy reference population. We demonstrate that MS lesions can be segmented as outliers from this intensity model of population and subject. We carried out extensive experiments with both synthetic and clinical data, and compared the performance of our new algorithm to those of state-of-the art techniques. We found this new approach leads to a substantial improvement in the sensitivity and specificity of lesion detection and segmentation.

CranialCloud: a cloud-based architecture to support trans-institutional collaborative efforts in neurodegenerative disorders

PURPOSE

Neurological diseases have a devastating impact on millions of individuals and their families. These diseases will continue to constitute a significant research focus for this century. The search for effective treatments and cures requires multiple teams of experts in clinical neurosciences, neuroradiology, engineering, and industry. Hence, the need to communicate a large amount of information with accuracy and precision is more necessary than ever for this specialty.

METHODS

In this paper, we present a distributed system that supports this vision, which we call the CranialVault Cloud (CranialCloud). It consists in a network of nodes, each with the capability to store and process data, that share the same spatial normalization processes, thus guaranteeing a common reference space. We detail and justify design choices, the architecture and functionality of individual nodes, the way these nodes interact, and how the distributed system can be used to support inter-institutional research.

RESULTS

We discuss the current state of the system that gathers data for more than 1,600 patients and how we envision it to grow.

CONCLUSION

We contend that the fastest way to find and develop promising treatments and cures is to permit teams of researchers to aggregate data, spatially normalize these data, and share them. The CranialVault system is a system that supports this vision.

Diffusion MRI abnormalities detection with orientation distribution functions: a multiple sclerosis longitudinal study

We propose a new algorithm for the voxelwise analysis of orientation distribution functions between one image and a group of reference images. It relies on a generic framework for the comparison of diffusion probabilities on the sphere, sampled from the underlying models. We demonstrate that this method, combined to dimensionality reduction through a principal component analysis, allows for more robust detection of lesions on simulated data when compared to classical tensor-based analysis. We then demonstrate the efficiency of this pipeline on the longitudinal comparison of multiple sclerosis patients at an early stage of the disease: right after their first clinically isolated syndrome (CIS) and three months later. We demonstrate the predictive value of ODF-based scores for the early detection of lesions that will appear or heal.

Multispectral image registration based on local canonical correlation analysis

Medical scans are today routinely acquired using multiple sequences or contrast settings, resulting in multispectral data. For the automatic analysis of this data, the evaluation of multispectral similarity is essential. So far, few concepts have been proposed to deal in a principled way with images containing multiple channels. Here, we present a new approach based on a well known statistical technique: canonical correlation analysis (CCA). CCA finds a mapping of two multidimensional variables into two new bases, which best represent the true underlying relations of the signals. In contrast to previously used metrics, it is therefore able to find new correlations based on linear combinations of multiple channels. We extend this concept to efficiently model local canonical correlation (LCCA) between image patches. This novel, more general similarity metric can be applied to images with an arbitrary number of channels. The most important property of LCCA is its invariance to affine transformations of variables. When used on local histograms, LCCA can also deal with multimodal similarity. We demonstrate the performance of our concept on challenging clinical multispectral datasets.

A logarithmic opinion pool based STAPLE algorithm for the fusion of segmentations with associated reliability weights

Pelvic floor dysfunction is common in women after childbirth and precise segmentation of magnetic resonance images (MRI) of the pelvic floor may facilitate diagnosis and treatment of patients. However, because of the complexity of its structures, manual segmentation of the pelvic floor is challenging and suffers from high inter and intra-rater variability of expert raters. Multiple template fusion algorithms are promising segmentation techniques for these types of applications, but they have been limited by imperfections in the alignment of templates to the target, and by template segmentation errors. A number of algorithms sought to improve segmentation performance by combining image intensities and template labels as two independent sources of information, carrying out fusion through local intensity weighted voting schemes. This class of approach is a form of linear opinion pooling, and achieves unsatisfactory performance for this application. We hypothesized that better decision fusion could be achieved by assessing the contribution of each template in comparison to a reference standard segmentation of the target image and developed a novel segmentation algorithm to enable automatic segmentation of MRI of the female pelvic floor. The algorithm achieves high performance by estimating and compensating for both imperfect registration of the templates to the target image and template segmentation inaccuracies. A local image similarity measure is used to infer a local reliability weight, which contributes to the fusion through a novel logarithmic opinion pooling. We evaluated our new algorithm in comparison to nine state-of-the-art segmentation methods and demonstrated our algorithm achieves the highest performance.

Three-dimensional hip cartilage quality assessment of morphology and dGEMRIC by planar maps and automated segmentation

OBJECTIVE

The quantitative interpretation of hip cartilage magnetic resonance imaging (MRI) has been limited by the difficulty of identifying and delineating the cartilage in a three-dimensional (3D) dataset, thereby reducing its routine usage. In this paper a solution is suggested by unfolding the cartilage to planar two-dimensional (2D) maps on which both morphology and biochemical degeneration patterns can be investigated across the entire hip joint.

DESIGN

Morphological TrueFISP and biochemical delayed gadolinium enhanced MRI of cartilage (dGEMRIC) hip images were acquired isotropically for 15 symptomatic subjects with mild or no radiographic osteoarthritis (OA). A multi-template based label fusion technique was used to automatically segment the cartilage tissue, followed by a geometric projection algorithm to generate the planar maps. The segmentation performance was investigated through a leave-one-out study, for two different fusion methods and as a function of the number of utilized templates.

RESULTS

For each of the generated planar maps, various patterns could be seen, indicating areas of healthy and degenerated cartilage. Dice coefficients for cartilage segmentation varied from 0.76 with four templates to 0.82 with 14 templates. Regional analysis suggests even higher segmentation performance in the superior half of the cartilage.

CONCLUSIONS

The proposed technique is the first of its kind to provide planar maps that enable straightforward quantitative assessment of hip cartilage morphology and dGEMRIC values. This technique may have important clinical applications for patient selection for hip preservation surgery, as well as for epidemiological studies of cartilage degeneration patterns. It is also shown that 10-15 templates are sufficient for accurate segmentation in this application.

Multi-parametric evaluation of the white matter maturation

In vivo evaluation of the brain white matter maturation is still a challenging task with no existing gold standards. In this article we propose an original approach to evaluate the early maturation of the white matter bundles, which is based on comparison of infant and adult groups using the Mahalanobis distance computed from four complementary MRI parameters: quantitative qT1 and qT2 relaxation times, longitudinal λ_║ and transverse λ_⊥ diffusivities from diffusion tensor imaging. Such multi-parametric approach is expected to better describe maturational asynchrony than conventional univariate approaches because it takes into account complementary dependencies of the parameters on different maturational processes, notably the decrease in water content and the myelination. Our approach was tested on 17 healthy infants (aged 3- to 21-week old) for 18 different bundles. It finely confirmed maturational asynchrony across the bundles: the spino-thalamic tract, the optic radiations, the cortico-spinal tract and the fornix have the most advanced maturation, while the superior longitudinal and arcuate fasciculi, the anterior limb of the internal capsule and the external capsule have the most delayed maturation. Furthermore, this approach was more reliable than univariate approaches as it revealed more maturational relationships between the bundles and did not violate a priori assumptions on the temporal order of the bundle maturation. Mahalanobis distances decreased exponentially with age in all bundles, with the only difference between them explained by different onsets of maturation. Estimation of these relative delays confirmed that the most dramatic changes occur during the first post-natal year.

Automated tract extraction via atlas based Adaptive Clustering

Advancements in imaging protocols such as the high angular resolution diffusion-weighted imaging (HARDI) and in tractography techniques are expected to cause an increase in the tract-based analyses. Statistical analyses over white matter tracts can contribute greatly towards understanding structural mechanisms of the brain since tracts are representative of connectivity pathways. The main challenge with tract-based studies is the extraction of the tracts of interest in a consistent and comparable manner over a large group of individuals without drawing the inclusion and exclusion regions of interest. In this work, we design a framework for automated extraction of white matter tracts. The framework introduces three main components, namely a connectivity based fiber representation, a fiber bundle atlas, and a clustering approach called Adaptive Clustering. The fiber representation relies on the connectivity signatures of fibers to establish an easy correspondence between different subjects. A group-wise clustering of these fibers that are represented by the connectivity signatures is then used to generate a fiber bundle atlas. Finally, Adaptive Clustering incorporates the previously generated clustering atlas as a prior, to cluster the fibers of a new subject automatically. Experiments on the HARDI scans of healthy individuals acquired repeatedly, demonstrate the applicability, reliability and the repeatability of our approach in extracting white matter tracts. By alleviating the seed region selection and the inclusion/exclusion ROI drawing requirements that are usually handled by trained radiologists, the proposed framework expands the range of possible clinical applications and establishes the ability to perform tract-based analyses with large samples.

Imaging connectivity: MRI and the structural networks of the brain

Magnetic resonance imaging (MRI) is a flexible and widely available neuroimaging technique. Structural MRI and diffusion MRI, in particular, provide information about connectivity between brain regions which may be combined to obtain a picture of entire neural networks, or the so-called connectome. In this review we outline the principles of MR-based connectivity analysis, discuss what relevant information it can provide for clinical and non-clinical neuroscience research, and outline some of the outstanding needs which future work will aim to meet.

A mathematical framework for the registration and analysis of multi-fascicle models for population studies of the brain microstructure

Diffusion tensor imaging (DTI) is unable to represent the diffusion signal arising from multiple crossing fascicles and freely diffusing water molecules. Generative models of the diffusion signal, such as multi-fascicle models, overcome this limitation by providing a parametric representation for the signal contribution of each population of water molecules. These models are of great interest in population studies to characterize and compare the brain microstructural properties. Central to population studies is the construction of an atlas and the registration of all subjects to it. However, the appropriate definition of registration and atlasing methods for multi-fascicle models have proven challenging. This paper proposes a mathematical framework to register and analyze multi-fascicle models. Specifically, we define novel operators to achieve interpolation, smoothing and averaging of multi-fascicle models. We also define a novel similarity metric to spatially align multi-fascicle models. Our framework enables simultaneous comparisons of different microstructural properties that are confounded in conventional DTI. The framework is validated on multi-fascicle models from 24 healthy subjects and 38 patients with tuberous sclerosis complex, 10 of whom have autism. We demonstrate the use of the multi-fascicle models registration and analysis framework in a population study of autism spectrum disorder.

Imaging connectivity: MRI and the structural networks of the brain

Magnetic resonance imaging (MRI) is a flexible and widely available neuroimaging technique. Structural MRI and diffusion MRI, in particular, provide information about connectivity between brain regions which may be combined to obtain a picture of entire neural networks, or the so-called connectome. In this review we outline the principles of MR-based connectivity analysis, discuss what relevant information it can provide for clinical and non-clinical neuroscience research, and outline some of the outstanding needs which future work will aim to meet.

Geometrical constraints for robust tractography selection

Tract-based analysis from DTI has become a widely employed procedure to study the white matter of the brain and its alterations in neurological and neurosurgical pathologies. Automatic tractography selection methods, where a subset of detected tracts corresponding to a specific white matter structure are selected, are a key component of the DTI processing pipeline. Using automatic tractography selection, repeatable results free of intra and inter-expert variability can be obtained rapidly, without the need for cumbersome manual segmentation. Many of the current approaches for automatic tractography selection rely on a previous registration procedure using an atlas; hence, these methods are likely very sensitive to the accuracy of the registration. In this paper we show that the performance of the registration step is critical to the overall result. This effect can in turn affect the calculation of scalar parameters derived subsequently from the selected tracts and often used in clinical practice; we show that such errors may be comparable in magnitude to the subtle differences found in clinical studies to differentiate between healthy and pathological. As an alternative, we propose a tractography selection method based on the use of geometrical constraints specific for each fiber bundle. Our experimental results show that the approach proposed performs with increased robustness and accuracy with respect to other approaches in the literature, particularly in the presence of imperfect registration.

Optimization of tractography of the optic radiations

Imaging and delineation of the optic radiations (OpRs) remains challenging, despite repeated attempts to achieve reliable validated tractography of this complex structure. Previous studies have used varying methods to generate representations of the OpR which differ markedly from one another and, frequently, from the OpR's known structure. We systematically examined the influence of a key variable that has differed across previous studies, the tractography seed region, in 13 adult participants (nine male; mean age 31 years; SD 8.7 years; range 16-47). First, we compared six seed regions at the lateral geniculate nucleus (LGN) and sagittal stratum based on the literature and known OpR anatomy. Three of the LGN regions seeded streamlines consistent with the OpR's three "bundles," whereas a fourth seeded streamlines consistent with each of the three bundles. The remaining two generated OpR streamlines unreliably and inconsistently. Two stratum regions seeded the radiations. This analysis identified a set of optimal regions of interest (ROI) for seeding OpR tractography and important inclusion and exclusion ROI. An optimized approach was then used to seed LGN regions to the stratum. The radiations, including streamlines consistent with Meyer's Loop, were streamlined in all cases. Streamlines extended 0.2 ± 2.4 mm anterior to the tip of the anterior horn of the lateral ventricle. These data suggest some existing approaches likely seed representations of the OpR that are visually plausible but do not capture all OpR components, and that using an optimized combination of regions seeded previously allows optimal mapping of this complex structure.

Impaired language pathways in tuberous sclerosis complex patients with autism spectrum disorders

The purpose of this study was to examine the relationship between language pathways and autism spectrum disorders (ASDs) in patients with tuberous sclerosis complex (TSC). An advanced diffusion-weighted magnetic resonance imaging (MRI) was performed on 42 patients with TSC and 42 age-matched controls. Using a validated automatic method, white matter language pathways were identified and microstructural characteristics were extracted, including fractional anisotropy (FA) and mean diffusivity (MD). Among 42 patients with TSC, 12 had ASD (29%). After controlling for age, TSC patients without ASD had a lower FA than controls in the arcuate fasciculus (AF); TSC patients with ASD had even a smaller FA, lower than the FA for those without ASD. Similarly, TSC patients without ASD had a greater MD than controls in the AF; TSC patients with ASD had even a higher MD, greater than the MD in those without ASD. It remains unclear why some patients with TSC develop ASD, while others have better language and socio-behavioral outcomes. Our results suggest that language pathway microstructure may serve as a marker of the risk of ASD in TSC patients. Impaired microstructure in language pathways of TSC patients may indicate the development of ASD, although prospective studies of language pathway development and ASD diagnosis in TSC remain essential.

Geodesic shape-based averaging

A new method for the geometrical averaging of labels or landmarks is presented. This method expands the shape-based averaging framework from an Euclidean to a geodesic based distance, incorporating a spatially varying similarity term as time cost. This framework has unique geometrical properties, making it ideal for propagating very small structures following rigorous labelling protocols. The method is used to automate the seeding and way-pointing of optic radiation tractography in DTI imaging. The propagated seeds and waypoints follow a strict clinical protocol by being geometrically constrained to one single slice and by guaranteeing spatial contiguity. The proposed method not only reduces the fragmentation of the propagated areas but also significantly increases the seed positioning accuracy and subsequent tractography results when compared to state-of-the-art label fusion techniques.