Disorganized Patterns of Sulcal Position in Fetal Brains with Agenesis of Corpus Callosum.

Fetuses with isolated agenesis of the corpus callosum (ACC) are associated with a broad spectrum of neurodevelopmental disability that cannot be specifically predicted in prenatal neuroimaging. We hypothesized that ACC may be associated with aberrant cortical folding. In this study, we determined altered patterning of early primary sulci development in fetuses with isolated ACC using novel quantitative sulcal pattern analysis which measures deviations of regional sulcal features (position, depth, and area) and their intersulcal relationships in 7 fetuses with isolated ACC (27.1 ± 3.8 weeks of gestation, mean ± SD) and 17 typically developing (TD) fetuses (25.7 ± 2.0 weeks) from normal templates. Fetuses with ACC showed significant alterations in absolute sulcal positions and relative intersulcal positional relationship compared to TD fetuses, which were not detected by traditional gyrification index. Our results reveal altered sulcal positional development even in isolated ACC that is present as early as the second trimester and continues throughout the fetal period. It might originate from altered white matter connections and portend functional variances in later life.

Brain extraction in pediatric ADC maps, toward characterizing neuro-development in multi-platform and multi-institution clinical images.

Apparent Diffusion Coefficient (ADC) maps can be used to characterize myelination and to detect abnormalities in the developing brain. However, given the normal variation in regional ADC with myelination, detection of abnormalities is difficult when based on visual assessment. Quantitative and automated analysis of pediatric ADC maps is thus desired but requires accurate brain extraction as the first step. Currently, most existing brain extraction methods are optimized for structural T1-weighted MR images of fully myelinated brains. Due to differences in age and image contrast, these approaches do not translate well to pediatric ADC maps. To address this problem, we present a multi-atlas brain extraction framework that has 1) specificity: designed and optimized specifically for pediatric ADC maps; 2) generality: applicable to multi-platform and multi-institution data, and to subjects at various neuro-developmental stages across the first 6 years of life; 3) accuracy: highly accurate compared to expert annotations; and 4) consistency: consistently accurate regardless of sources of data and ages of subjects. We show how we achieve these goals, via optimizing major components in a multi-atlas brain extraction framework, and via developing and evaluating new criteria for its atlas ranking component. Moreover, we demonstrate that these goals can be achieved with a fixed set of atlases and a fixed set of parameters, which opens doors for our optimized framework to be used in large-scale and multi-institution neuro-developmental and clinical studies. In a pilot study, we use this framework in a dataset containing scanner-generated ADC maps from 308 pediatric patients collected during the course of routine clinical care. Our framework leads to successful quantifications of the changes in whole-brain volumes and mean ADC values across the first 6 years of life.

Quantification and discrimination of abnormal sulcal patterns in polymicrogyria.

Polymicrogyria (PMG) is a malformation of cortical development characterized by an irregular gyral pattern and its diagnosis and severity have been qualitatively judged by visual inspection of imaging features. We aimed to provide a quantitative description of abnormal sulcal patterns for individual PMG brains using our sulcal graph-based analysis and examined the association with language impairment. The sulcal graphs were constructed from magnetic resonance images in 26 typical developing and 18 PMG subjects and the similarity between sulcal graphs was computed by using their geometric and topological features. The similarities between typical and PMG groups were significantly lower than the similarities measured within the typical group. Furthermore, more lobar regions were determined to be abnormal in most patients when compared with the visual diagnosis of PMG involvement, suggesting that PMG may have more global effects on cortical folding than previously expected. Among the PMG, the group with intact language development showed sulcal patterns more closely matched with the typical than the impaired group in the left parietal lobe. Our approach shows the potential to provide a quantitative means for detecting the severity and extent of involvement of cortical malformation and a greater understanding of genotype-phenotype and clinical-imaging features correlations.

Regional infant brain development: an MRI-based morphometric analysis in 3 to 13 month olds.

Elucidation of infant brain development is a critically important goal given the enduring impact of these early processes on various domains including later cognition and language. Although infants' whole-brain growth rates have long been available, regional growth rates have not been reported systematically. Accordingly, relatively less is known about the dynamics and organization of typically developing infant brains. Here we report global and regional volumetric growth of cerebrum, cerebellum, and brainstem with gender dimorphism, in 33 cross-sectional scans, over 3 to 13 months, using T1-weighted 3-dimensional spoiled gradient echo images and detailed semi-automated brain segmentation. Except for the midbrain and lateral ventricles, all absolute volumes of brain regions showed significant growth, with 6 different patterns of volumetric change. When normalized to the whole brain, the regional increase was characterized by 5 differential patterns. The putamen, cerebellar hemispheres, and total cerebellum were the only regions that showed positive growth in the normalized brain. Our results show region-specific patterns of volumetric change and contribute to the systematic understanding of infant brain development. This study greatly expands our knowledge of normal development and in future may provide a basis for identifying early deviation above and beyond normative variation that might signal higher risk for neurological disorders.

Morphological Features of Language Regions in Individuals with Tuberous Sclerosis Complex.

A significant number of individuals with tuberous sclerosis complex (TSC) exhibit language difficulties. Here, we examined the language-related brain morphometry in 59 participants (7 participants with TSC and comorbid autism spectrum disorder (ASD) (TSC + ASD), 13 with TSC but no ASD (TSC-ASD), 10 with ASD-only (ASD), and 29 typically developing (TD) controls). A hemispheric asymmetry was noted in surface area and gray matter volume of several cortical language areas in TD, ASD, and TSC-ASD groups, but not in TSC + ASD group. TSC + ASD group demonstrated increased cortical thickness and curvature values in multiple language regions for both hemispheres, compared to other groups. After controlling for tuber load in the TSC groups, within-group differences stayed the same but the differences between TSC-ASD and TSC + ASD were no longer statistically significant. These preliminary findings suggest that comorbid ASD in TSC as well as tuber load in TSC is associated with changes in the morphometry of language regions. Future studies with larger sample sizes will be needed to confirm these findings.

A quantitative method for correlating observations of decreased apparent diffusion coefficient with elevated cerebral blood perfusion in newborns presenting cerebral ischemic insults.

In patients presenting with cerebral ischemic injury, the outcome of injured brain tissue quantified as decreased apparent diffusion coefficient (ADC) may depend on associated alterations in cerebral blood perfusion (CBP). This study proposes a non-biased method to quantify associations between ADC and CBP in newborns with global or focal cerebral ischemia. The study population consisted of nine neonates (age: 0 to 3 days) presenting with clinical and imaging evidence of ischemia (seven with global hypoxic ischemia, and two with focal arterial ischemic stroke) with decreased ADC. Six newborns without diffusion abnormalities on magnetic resonance (MR) imaging served as a comparative cohort (age: 0 days to 4 weeks). All patients underwent MR imaging including diffusion weighted imaging (DWI) to determine ADC and axial arterial spin labeling (ASL) to determine CBP. An algorithm was developed that uses the B0 volume from the DWI raw data as a reference, co-registers the ADC and ASL-CBP data to the B0, generates mask filters, and finally performs a statistical analysis to automatically select regions of interest (ROIs) with ADC or ASL-CBP values that deviate significantly from the rest of the brain. If ROIs are identified in this analysis, the algorithm then evaluates correlation based on ROI location and volume. A significant correlation was found between decreased ADC and elevated ASL-CBP with regions of elevated ASL-CBP typically larger than the corresponding ADC abnormality. The association between decreased diffusivity and increased ASL-CBP suggests that, for this cohort, cerebral ischemia is associated with hyperperfusion.

Quantitative comparison and analysis of sulcal patterns using sulcal graph matching: a twin study.

The global pattern of cortical sulci provides important information on brain development and functional compartmentalization. Sulcal patterns are routinely used to determine fetal brain health and detect cerebral malformations. We present a quantitative method for automatically comparing and analyzing the sulcal pattern between individuals using a graph matching approach. White matter surfaces were reconstructed from volumetric T1 MRI data and sulcal pits, the deepest points in local sulci, were identified on this surface. The sulcal pattern was then represented as a graph structure with sulcal pits as nodes. The similarity between graphs was computed with a spectral-based matching algorithm by using the geometric features of nodes (3D position, depth and area) and their relationship. In particular, we exploited the feature of graph topology (the number of edges and the paths between nodes) to highlight the interrelated arrangement and patterning of sulcal folds. We applied this methodology to 48 monozygotic twins and showed that the similarity of the sulcal graphs in twin pairs was significantly higher than in unrelated pairs for all hemispheres and lobar regions, consistent with a genetic influence on sulcal patterning. This novel approach has the potential to provide a quantitative and reliable means to compare sulcal patterns.

Multi-channel attention-fusion neural network for brain age estimation: Accuracy, generality, and interpretation with 16,705 healthy MRIs across lifespan.

Brain age estimated by machine learning from T1-weighted magnetic resonance images (T1w MRIs) can reveal how brain disorders alter brain aging and can help in the early detection of such disorders. A fundamental step is to build an accurate age estimator from healthy brain MRIs. We focus on this step, and propose a framework to improve the accuracy, generality, and interpretation of age estimation in healthy brain MRIs. For accuracy, we used one of the largest sample sizes (N = 16,705). For each subject, our proposed algorithm first explicitly splits the T1w image, which has been commonly treated as a single-channel 3D image in other studies, into two 3D image channels representing contrast and morphometry information. We further proposed a "fusion-with-attention" deep learning convolutional neural network (FiA-Net) to learn how to best fuse the contrast and morphometry image channels. FiA-Net recognizes varying contributions across image channels at different brain anatomy and different feature layers. In contrast, multi-channel fusion does not exist for brain age estimation, and is mostly attention-free in other medical image analysis tasks (e.g., image synthesis, or segmentation), where treating channels equally may not be optimal. For generality, we used lifespan data 0-97 years of age for real-world utility; and we thoroughly tested FiA-Net for multi-site and multi-scanner generality by two phases of cross-validations in discovery and replication data, compared to most other studies with only one phase of cross-validation. For interpretation, we directly measured each artificial neuron's correlation with the chronological age, compared to other studies looking at the saliency of features where salient features may or may not predict age. Overall, FiA-Net achieved a mean absolute error (MAE) of 3.00 years and Pearson correlation r=0.9840 with known chronological ages in healthy brain MRIs 0-97 years of age, comparing favorably with state-of-the-art algorithms and studies for accuracy and generality across sites and datasets. We also provided interpretations on how different artificial neurons and real neuroanatomy contribute to the age estimation.

Neuropsychiatric disease classification using functional connectomics - results of the connectomics in neuroimaging transfer learning challenge.

Large, open-source datasets, such as the Human Connectome Project and the Autism Brain Imaging Data Exchange, have spurred the development of new and increasingly powerful machine learning approaches for brain connectomics. However, one key question remains: are we capturing biologically relevant and generalizable information about the brain, or are we simply overfitting to the data? To answer this, we organized a scientific challenge, the Connectomics in NeuroImaging Transfer Learning Challenge (CNI-TLC), held in conjunction with MICCAI 2019. CNI-TLC included two classification tasks: (1) diagnosis of Attention-Deficit/Hyperactivity Disorder (ADHD) within a pre-adolescent cohort; and (2) transference of the ADHD model to a related cohort of Autism Spectrum Disorder (ASD) patients with an ADHD comorbidity. In total, 240 resting-state fMRI (rsfMRI) time series averaged according to three standard parcellation atlases, along with clinical diagnosis, were released for training and validation (120 neurotypical controls and 120 ADHD). We also provided Challenge participants with demographic information of age, sex, IQ, and handedness. The second set of 100 subjects (50 neurotypical controls, 25 ADHD, and 25 ASD with ADHD comorbidity) was used for testing. Classification methodologies were submitted in a standardized format as containerized Docker images through ChRIS, an open-source image analysis platform. Utilizing an inclusive approach, we ranked the methods based on 16 metrics: accuracy, area under the curve, F1-score, false discovery rate, false negative rate, false omission rate, false positive rate, geometric mean, informedness, markedness, Matthew's correlation coefficient, negative predictive value, optimized precision, precision, sensitivity, and specificity. The final rank was calculated using the rank product for each participant across all measures. Furthermore, we assessed the calibration curves of each methodology. Five participants submitted their method for evaluation, with one outperforming all other methods in both ADHD and ASD classification. However, further improvements are still needed to reach the clinical translation of functional connectomics. We have kept the CNI-TLC open as a publicly available resource for developing and validating new classification methodologies in the field of connectomics.

BRAIN AGE ESTIMATION USING LSTM ON CHILDREN'S BRAIN MRI.

Brain age prediction based on children's brain MRI is an important biomarker for brain health and brain development analysis. In this paper, we consider the 3D brain MRI volume as a sequence of 2D images and propose a new framework using the recurrent neural network for brain age estimation. The proposed method is named as 2D-ResNet18+Long short-term memory (LSTM), which consists of four parts: 2D ResNet18 for feature extraction on 2D images, a pooling layer for feature reduction over the sequences, an LSTM layer, and a final regression layer. We apply the proposed method on a public multisite NIH-PD dataset and evaluate generalization on a second multisite dataset, which shows that the proposed 2D-ResNet18+LSTM method provides better results than traditional 3D based neural network for brain age estimation.

Cortical surface shape analysis based on spherical wavelets.

In vivo quantification of neuroanatomical shape variations is possible due to recent advances in medical imaging and has proven useful in the study of neuropathology and neurodevelopment. In this paper, we apply a spherical wavelet transformation to extract shape features of cortical surfaces reconstructed from magnetic resonance images (MRIs) of a set of subjects. The spherical wavelet transformation can characterize the underlying functions in a local fashion in both space and frequency, in contrast to spherical harmonics that have a global basis set. We perform principal component analysis (PCA) on these wavelet shape features to study patterns of shape variation within normal population from coarse to fine resolution. In addition, we study the development of cortical folding in newborns using the Gompertz model in the wavelet domain, which allows us to characterize the order of development of large-scale and finer folding patterns independently. Given a limited amount of training data, we use a regularization framework to estimate the parameters of the Gompertz model to improve the prediction performance on new data. We develop an efficient method to estimate this regularized Gompertz model based on the Broyden-Fletcher-Goldfarb-Shannon (BFGS) approximation. Promising results are presented using both PCA and the folding development model in the wavelet domain. The cortical folding development model provides quantitative anatomic information regarding macroscopic cortical folding development and may be of potential use as a biomarker for early diagnosis of neurologic deficits in newborns.

Reusable Client-Side JavaScript Modules for Immersive Web-Based Real-Time Collaborative Neuroimage Visualization.

In this paper we present a web-based software solution to the problem of implementing real-time collaborative neuroimage visualization. In both clinical and research settings, simple and powerful access to imaging technologies across multiple devices is becoming increasingly useful. Prior technical solutions have used a server-side rendering and push-to-client model wherein only the server has the full image dataset. We propose a rich client solution in which each client has all the data and uses the Google Drive Realtime API for state synchronization. We have developed a small set of reusable client-side object-oriented JavaScript modules that make use of the XTK toolkit, a popular open-source JavaScript library also developed by our team, for the in-browser rendering and visualization of brain image volumes. Efficient realtime communication among the remote instances is achieved by using just a small JSON object, comprising a representation of the XTK image renderers' state, as the Google Drive Realtime collaborative data model. The developed open-source JavaScript modules have already been instantiated in a web-app called MedView, a distributed collaborative neuroimage visualization application that is delivered to the users over the web without requiring the installation of any extra software or browser plugin. This responsive application allows multiple physically distant physicians or researchers to cooperate in real time to reach a diagnosis or scientific conclusion. It also serves as a proof of concept for the capabilities of the presented technological solution.

ChRIS--A web-based neuroimaging and informatics system for collecting, organizing, processing, visualizing and sharing of medical data.

The utility of web browsers for general purpose computing, long anticipated, is only now coming into fruition. In this paper we present a web-based medical image data and information management software platform called ChRIS ([Boston] Children's Research Integration System). ChRIS' deep functionality allows for easy retrieval of medical image data from resources typically found in hospitals, organizes and presents information in a modern feed-like interface, provides access to a growing library of plugins that process these data - typically on a connected High Performance Compute Cluster, allows for easy data sharing between users and instances of ChRIS and provides powerful 3D visualization and real time collaboration.

Longitudinal changes in diffusion properties in white matter pathways of children with tuberous sclerosis complex.

BACKGROUND: Abnormal white matter development in patients with tuberous sclerosis complex, a multisystem hamartomatous disorder caused by aberrant neural proliferation and axonal maturation, may be associated with poorer neurocognitive outcomes. The purpose of this study is to identify predictors of longitudinal changes in diffusion properties of white matter tracts in patients with tuberous sclerosis complex. METHODS: Diffusion magnetic resonance imaging was carried out in 17 subjects with tuberous sclerosis complex (mean age, 7.2 ± 4.4 years) with at least two magnetic resonance imaging scans (mean number of days between scans, 419.4 ± 105.4). There were 10 males; 5 of 17 had autism spectrum disorder and 10 of 17 had epilepsy. Regions of interest were placed to delineate the internal capsule/corona radiata, cingulum, and corpus callosum. The outcomes were mean change in apparent diffusion coefficient and fractional anisotropy. Data were analyzed using Pearson's correlation and multiple linear regression analyses. RESULTS: Gender was a significant predictor of mean change in apparent diffusion coefficient in the left internal capsule, right and left cingulum bundles, and corpus callosum and a significant predictor of mean change in fractional anisotropy in the corpus callosum. Epilepsy was a significant predictor of mean change in apparent diffusion coefficient in the left internal capsule. Autism spectrum disorder was not predictive of diffusion changes in any of the studied pathways. CONCLUSION: Clinical variables, including gender and epilepsy, have an effect on the development of white matter pathways. These variables should be taken into consideration when counseling tuberous sclerosis complex patients and in future imaging studies in this population.

Assessment of the frequency-domain multi-distance method to evaluate the brain optical properties: Monte Carlo simulations from neonate to adult.

The near infrared spectroscopy (NIRS) frequency-domain multi-distance (FD-MD) method allows for the estimation of optical properties in biological tissue using the phase and intensity of radiofrequency modulated light at different source-detector separations. In this study, we evaluated the accuracy of this method to retrieve the absorption coefficient of the brain at different ages. Synthetic measurements were generated with Monte Carlo simulations in magnetic resonance imaging (MRI)-based heterogeneous head models for four ages: newborn, 6 and 12 month old infants, and adult. For each age, we determined the optimal set of source-detector separations and estimated the corresponding errors. Errors arise from different origins: methodological (FD-MD) and anatomical (curvature, head size and contamination by extra-cerebral tissues). We found that the brain optical absorption could be retrieved with an error between 8-24% in neonates and infants, while the error increased to 19-44% in adults over all source-detector distances. The dominant contribution to the error was found to be the head curvature in neonates and infants, and the extra-cerebral tissues in adults.

Intrinsic curvature: a marker of millimeter-scale tangential cortico-cortical connectivity?

In this paper, we draw a link between cortical intrinsic curvature and the distributions of tangential connection lengths. We suggest that differential rates of surface expansion not only lead to intrinsic curvature of the cortical sheet, but also to differential inter-neuronal spacing. We propose that there follows a consequential change in the profile of neuronal connections: specifically an enhancement of the tendency towards proportionately more short connections. Thus, the degree of cortical intrinsic curvature may have implications for short-range connectivity.

Cortical Surface Shape Analysis Based on Spherical Wavelet Transformation.

Shape analysis of neuroanatomical structures has proven useful in the study of neuropathology and neurodevelopment. Advances in medical imaging have made it possible to study this shape variation in vivo. In this paper, we propose the use of a spherical wavelet transformation to extract cortical surface shape features, as wavelets can characterize the underlying functions in a local fashion in both space and frequency. Our results demonstrate the utility of the wavelet approach in both detecting the spatial scale and pattern of shape variation in synthetic data, and for quantifying and visualizing shape variations of cortical surface models in subject populations.