Connectomic consistency: a systematic stability analysis of structural and functional connectivity

OBJECTIVE

Connectomics, the study of brain connectivity, has become an indispensable tool in neuroscientific research as it provides insights into brain organization. Connectomes are generated using different modalities such as diffusion MRI to capture structural organization of the brain or functional MRI to elaborate brain's functional organization. Understanding links between structural and functional organizations is crucial in explaining how observed behavior emerges from the underlying neurobiological mechanisms. Many studies have investigated how these two organizations relate to each other; however, we still lack a comparative understanding on how much variation should be expected in the two modalities, both between people and within a single person across scans.

APPROACH

In this study, we systematically analyzed the consistency of connectomes, that is the similarity between connectomes in terms of individual connections between brain regions and in terms of overall network topology. We present a comprehensive study of consistency in connectomes for a single subject examined longitudinally and across a large cohort of subjects cross-sectionally, in structure and function separately. Within structural connectomes, we compared connectomes generated by different tracking algorithms, parcellations, edge weighting schemes, and edge pruning techniques. In functional connectomes, we compared full, positive, and negative connectivity separately along with thresholding of weak edges. We evaluated consistency using correlation (incorporating information at the level of individual edges) and graph matching accuracy (evaluating connectivity at the level of network topology). We also examined the consistency of connectomes that are generated using different communication schemes.

MAIN RESULTS

Our results demonstrate varying degrees of consistency for the two modalities, with structural connectomes showing higher consistency than functional connectomes. Moreover, we observed a wide variation in consistency depending on how connectomes are generated.

SIGNIFICANCE

Our study sets a reference point for consistency of connectome types, which is especially important for structure-function coupling studies in evaluating mismatches between modalities.

Freewater estimatoR using iNtErpolated iniTialization (FERNET): Characterizing peritumoral edema using clinically feasible diffusion MRI data

Characterization of healthy versus pathological tissue in the peritumoral area is confounded by the presence of edema, making free water estimation the key concern in modeling tissue microstructure. Most methods that model tissue microstructure are either based on advanced acquisition schemes not readily available in the clinic or are not designed to address the challenge of edema. This underscores the need for a robust free water elimination (FWE) method that estimates free water in pathological tissue but can be used with clinically prevalent single-shell diffusion tensor imaging data. FWE in single-shell data requires the fitting of a bi-compartment model, which is an ill-posed problem. Its solution requires optimization, which relies on an initialization step. We propose a novel initialization approach for FWE, FERNET, which improves the estimation of free water in edematous and infiltrated peritumoral regions, using single-shell diffusion MRI data. The method has been extensively investigated on simulated data and healthy dataset. Additionally, it has been applied to clinically acquired data from brain tumor patients to characterize the peritumoral region and improve tractography in it.

Analysis of Consistency in Structural and Functional Connectivity of Human Brain

Analysis of structural and functional connectivity of brain has become a fundamental approach in neuroscientific research. Despite several studies reporting consistent similarities as well as differences for structural and resting state (rs) functional connectomes, a comparative investigation of connectomic consistency between the two modalities is still lacking. Nonetheless, connectomic analysis comprising both connectivity types necessitate extra attention as consistency of connectivity differs across modalities, possibly affecting the interpretation of the results. In this study, we present a comprehensive analysis of consistency in structural and rs-functional connectomes obtained from longitudinal diffusion MRI and rs-fMRI data of a single healthy subject. We contrast consistency of deterministic and probabilistic tracking with that of full, positive, and negative functional connectivities across various connectome generation schemes, using correlation as a measure of consistency.

Distance-based analysis of variance for brain connectivity

The field of neuroimaging dedicated to mapping connections in the brain is increasingly being recognized as key for understanding neurodevelopment and pathology. Networks of these connections are quantitatively represented using complex structures, including matrices, functions, and graphs, which require specialized statistical techniques for estimation and inference about developmental and disorder-related changes. Unfortunately, classical statistical testing procedures are not well suited to high-dimensional testing problems. In the context of global or regional tests for differences in neuroimaging data, traditional analysis of variance (ANOVA) is not directly applicable without first summarizing the data into univariate or low-dimensional features, a process that might mask the salient features of high-dimensional distributions. In this work, we consider a general framework for two-sample testing of complex structures by studying generalized within-group and between-group variances based on distances between complex and potentially high-dimensional observations. We derive an asymptotic approximation to the null distribution of the ANOVA test statistic, and conduct simulation studies with scalar and graph outcomes to study finite sample properties of the test. Finally, we apply our test to our motivating study of structural connectivity in autism spectrum disorder.

QC-Automator: Deep Learning-Based Automated Quality Control for Diffusion MR Images

Quality assessment of diffusion MRI (dMRI) data is essential prior to any analysis, so that appropriate pre-processing can be used to improve data quality and ensure that the presence of MRI artifacts do not affect the results of subsequent image analysis. Manual quality assessment of the data is subjective, possibly error-prone, and infeasible, especially considering the growing number of consortium-like studies, underlining the need for automation of the process. In this paper, we have developed a deep-learning-based automated quality control (QC) tool, QC-Automator, for dMRI data, that can handle a variety of artifacts such as motion, multiband interleaving, ghosting, susceptibility, herringbone, and chemical shifts. QC-Automator uses convolutional neural networks along with transfer learning to train the automated artifact detection on a labeled dataset of ∼332,000 slices of dMRI data, from 155 unique subjects and 5 scanners with different dMRI acquisitions, achieving a 98% accuracy in detecting artifacts. The method is fast and paves the way for efficient and effective artifact detection in large datasets. It is also demonstrated to be replicable on other datasets with different acquisition parameters.

Deviation from normative brain development is associated with symptom severity in autism spectrum disorder

Background

Autism spectrum disorder (ASD) is a heterogeneous neurodevelopmental condition. The degree to which the brain development in ASD deviates from typical brain development, and how this deviation relates to observed behavioral outcomes at the individual level are not well-studied. We hypothesize that the degree of deviation from typical brain development of an individual with ASD would relate to observed symptom severity.

Methods

The developmental changes in anatomical (cortical thickness, surface area, and volume) and diffusion metrics (fractional anisotropy and apparent diffusion coefficient) were compared between a sample of ASD (n = 247) and typically developing children (TDC) (n = 220) aged 6–25. Machine learning was used to predict age (brain age) from these metrics in the TDC sample, to define a normative model of brain development. This model was then used to compute brain age in the ASD sample. The difference between chronological age and brain age was considered a developmental deviation index (DDI), which was then correlated with ASD symptom severity.

Results

Machine learning model trained on all five metrics accurately predicted age in the TDC (r = 0.88) and the ASD (r = 0.85) samples, with dominant contributions to the model from the diffusion metrics. Within the ASD group, the DDI derived from fractional anisotropy was correlated with ASD symptom severity (r = − 0.2), such that individuals with the most advanced brain age showing the lowest severity, and individuals with the most delayed brain age showing the highest severity.

Limitations

This work investigated only linear relationships between five specific brain metrics and only one measure of ASD symptom severity in a limited age range. Reported effect sizes are moderate. Further work is needed to investigate developmental differences in other age ranges, other aspects of behavior, other neurobiological measures, and in an independent sample before results can be clinically applicable.

Conclusions

Findings demonstrate that the degree of deviation from typical brain development relates to ASD symptom severity, partially accounting for the observed heterogeneity in ASD. Our approach enables characterization of each individual with reference to normative brain development and identification of distinct developmental subtypes, facilitating a better understanding of developmental heterogeneity in ASD.

Neuroimaging Findings in US Government Personnel With Possible Exposure to Directional Phenomena in Havana, Cuba

IMPORTANCE

United States government personnel experienced potential exposures to uncharacterized directional phenomena while serving in Havana, Cuba, from late 2016 through May 2018. The underlying neuroanatomical findings have not been described.

OBJECTIVE

To examine potential differences in brain tissue volume, microstructure, and functional connectivity in government personnel compared with individuals not exposed to directional phenomena.

DESIGN, SETTING, AND PARTICIPANTS

Forty government personnel (patients) who were potentially exposed and experienced neurological symptoms underwent evaluation at a US academic medical center from August 21, 2017, to June 8, 2018, including advanced structural and functional magnetic resonance imaging analytics. Findings were compared with imaging findings of 48 demographically similar healthy controls.

EXPOSURES

Potential exposure to uncharacterized directional phenomena of unknown etiology, manifesting as pressure, vibration, or sound.

MAIN OUTCOMES AND MEASURES

Potential imaging-based differences between patients and controls with regard to (1) white matter and gray matter total and regional brain volumes, (2) cerebellar tissue microstructure metrics (eg, mean diffusivity), and (3) functional connectivity in the visuospatial, auditory, and executive control subnetworks.

RESULTS

Imaging studies were completed for 40 patients (mean age, 40.4 years; 23 [57.5%] men; imaging performed a median of 188 [range, 4-403] days after initial exposure) and 48 controls (mean age, 37.6 years; 33 [68.8%] men). Mean whole brain white matter volume was significantly smaller in patients compared with controls (patients: 542.22 cm3; controls: 569.61 cm3; difference, -27.39 [95% CI, -37.93 to -16.84] cm3; P < .001), with no significant difference in the whole brain gray matter volume (patients: 698.55 cm3; controls: 691.83 cm3; difference, 6.72 [95% CI, -4.83 to 18.27] cm3; P = .25). Among patients compared with controls, there were significantly greater ventral diencephalon and cerebellar gray matter volumes and significantly smaller frontal, occipital, and parietal lobe white matter volumes; significantly lower mean diffusivity in the inferior vermis of the cerebellum (patients: 7.71 × 10-4 mm2/s; controls: 8.98 × 10-4 mm2/s; difference, -1.27 × 10-4 [95% CI, -1.93 × 10-4 to -6.17 × 10-5] mm2/s; P < .001); and significantly lower mean functional connectivity in the auditory subnetwork (patients: 0.45; controls: 0.61; difference, -0.16 [95% CI, -0.26 to -0.05]; P = .003) and visuospatial subnetwork (patients: 0.30; controls: 0.40; difference, -0.10 [95% CI, -0.16 to -0.04]; P = .002) but not in the executive control subnetwork (patients: 0.24; controls: 0.25; difference: -0.016 [95% CI, -0.04 to 0.01]; P = .23).

CONCLUSIONS AND RELEVANCE

Among US government personnel in Havana, Cuba, with potential exposure to directional phenomena, compared with healthy controls, advanced brain magnetic resonance imaging revealed significant differences in whole brain white matter volume, regional gray and white matter volumes, cerebellar tissue microstructural integrity, and functional connectivity in the auditory and visuospatial subnetworks but not in the executive control subnetwork. The clinical importance of these differences is uncertain and may require further study.

Estimating degradation-related settlement in two landfill-reclaimed soils by sand-salt analogues

Landfill reclaimed soil here refers to largely degraded materials excavated from old landfill sites, which after processing can be reinstated as more competent fill, thereby restoring the former landfill space. The success of the process depends on the presence of remaining degradable particles and their influence on settlement. Tests on salt-sand mixtures, from which the salt is removed, have been used to quantify the impact of particle loss on settlement. Where the amount of particle loss is small, say 10% by mass or less, settlements are small and apparently independent of lost particle size. A conceptual model is presented to explain this behaviour in terms of nestling particles and strong force chains. At higher percentages of lost particles, greater rates of settlement together with some sensitivity to particle size were observed. The conceptual model was then applied to two landfill reclaimed soils, the long-term settlements of which were found to be consistent with the conceptual model suggesting that knowledge of particle content and relative size are sufficient to estimate the influence of degradable particles in landfill reclaimed soils.

Ca 125 and Survival in Ovarian Cancer: Preliminary Communication

CA 125 is an epithelial membrane marker which can be detected in the serum of patients with ovarian cancer and which may reflect tumour burden. In a group of 42 patients, the level of this marker has also been shown to be significantly related to survival.

Cancer imaging phenomics toolkit: quantitative imaging analytics for precision diagnostics and predictive modeling of clinical outcome

The growth of multiparametric imaging protocols has paved the way for quantitative imaging phenotypes that predict treatment response and clinical outcome, reflect underlying cancer molecular characteristics and spatiotemporal heterogeneity, and can guide personalized treatment planning. This growth has underlined the need for efficient quantitative analytics to derive high-dimensional imaging signatures of diagnostic and predictive value in this emerging era of integrated precision diagnostics. This paper presents cancer imaging phenomics toolkit (CaPTk), a new and dynamically growing software platform for analysis of radiographic images of cancer, currently focusing on brain, breast, and lung cancer. CaPTk leverages the value of quantitative imaging analytics along with machine learning to derive phenotypic imaging signatures, based on two-level functionality. First, image analysis algorithms are used to extract comprehensive panels of diverse and complementary features, such as multiparametric intensity histogram distributions, texture, shape, kinetics, connectomics, and spatial patterns. At the second level, these quantitative imaging signatures are fed into multivariate machine learning models to produce diagnostic, prognostic, and predictive biomarkers. Results from clinical studies in three areas are shown: (i) computational neuro-oncology of brain gliomas for precision diagnostics, prediction of outcome, and treatment planning; (ii) prediction of treatment response for breast and lung cancer, and (iii) risk assessment for breast cancer.

Three-Dimensional Printed Modeling of Diffuse Low-Grade Gliomas and Associated White Matter Tract Anatomy

BACKGROUND

Diffuse low-grade gliomas (DLGGs) represent several pathological entities that infiltrate and invade cortical and subcortical structures in the brain.

OBJECTIVE

To describe methods for rapid prototyping of DLGGs and surgically relevant anatomy.

METHODS

Using high-definition imaging data and rapid prototyping technologies, we were able to generate 3 patient DLGGs to scale and represent the associated white matter tracts in 3 dimensions using advanced diffusion tensor imaging techniques.

RESULTS

This report represents a novel application of 3-dimensional (3-D) printing in neurosurgery and a means to model individualized tumors in 3-D space with respect to subcortical white matter tract anatomy. Faculty and resident evaluations of this technology were favorable at our institution.

CONCLUSION

Developing an understanding of the anatomic relationships existing within individuals is fundamental to successful neurosurgical therapy. Imaging-based rapid prototyping may improve on our ability to plan for and treat complex neuro-oncologic pathology.

Harmonization of multi-site diffusion tensor imaging data

Diffusion tensor imaging (DTI) is a well-established magnetic resonance imaging (MRI) technique used for studying microstructural changes in the white matter. As with many other imaging modalities, DTI images suffer from technical between-scanner variation that hinders comparisons of images across imaging sites, scanners and over time. Using fractional anisotropy (FA) and mean diffusivity (MD) maps of 205 healthy participants acquired on two different scanners, we show that the DTI measurements are highly site-specific, highlighting the need of correcting for site effects before performing downstream statistical analyses. We first show evidence that combining DTI data from multiple sites, without harmonization, may be counter-productive and negatively impacts the inference. Then, we propose and compare several harmonization approaches for DTI data, and show that ComBat, a popular batch-effect correction tool used in genomics, performs best at modeling and removing the unwanted inter-site variability in FA and MD maps. Using age as a biological phenotype of interest, we show that ComBat both preserves biological variability and removes the unwanted variation introduced by site. Finally, we assess the different harmonization methods in the presence of different levels of confounding between site and age, in addition to test robustness to small sample size studies.

Whole brain white matter connectivity analysis using machine learning: An application to autism

In this paper, we propose an automated white matter connectivity analysis method for machine learning classification and characterization of white matter abnormality via identification of discriminative fiber tracts. The proposed method uses diffusion MRI tractography and a data-driven approach to find fiber clusters corresponding to subdivisions of the white matter anatomy. Features extracted from each fiber cluster describe its diffusion properties and are used for machine learning. The method is demonstrated by application to a pediatric neuroimaging dataset from 149 individuals, including 70 children with autism spectrum disorder (ASD) and 79 typically developing controls (TDC). A classification accuracy of 78.33% is achieved in this cross-validation study. We investigate the discriminative diffusion features based on a two-tensor fiber tracking model. We observe that the mean fractional anisotropy from the second tensor (associated with crossing fibers) is most affected in ASD. We also find that local along-tract (central cores and endpoint regions) differences between ASD and TDC are helpful in differentiating the two groups. These altered diffusion properties in ASD are associated with multiple robustly discriminative fiber clusters, which belong to several major white matter tracts including the corpus callosum, arcuate fasciculus, uncinate fasciculus and aslant tract; and the white matter structures related to the cerebellum, brain stem, and ventral diencephalon. These discriminative fiber clusters, a small part of the whole brain tractography, represent the white matter connections that could be most affected in ASD. Our results indicate the potential of a machine learning pipeline based on white matter fiber clustering.