Cortical thickness analysis in autism with heat kernel smoothing

Graph-based prototype inverse-projection for identifying cortical sulcal pattern abnormalities in congenital heart disease

Examining the altered arrangement and patterning of sulcal folds offers insights into the mechanisms of neurodevelopmental differences in psychiatric and neurological disorders. Previous sulcal pattern analysis used spectral graph matching of sulcal pit-based graph structures to assess deviations from normative sulcal patterns. However, challenges exist, including the absence of a standard criterion for defining a typical reference set, time-consuming cost of graph matching, user-defined feature weight sets, and assumptions about uniform node distribution. We developed a deep learning-based sulcal pattern analysis to address these challenges by adapting prototype-based graph neural networks to sulcal pattern graphs. Additionally, we proposed a prototype inverse-projection for better interpretability. Unlike other prototype-based models, our approach inversely projects prototypes onto individual node representations to calculate the inverse-projection weights, enabling efficient visualization of prototypes and focusing the model on selective regions. We evaluated our method through a classification task between healthy controls (n = 174, age = 15.4 ±1.9 [mean ± standard deviation, years]) and patients with congenital heart disease (n = 345, age = 15.8 ±4.7) from four cohort studies and a public dataset. Our approach demonstrated superior classification performance compared to other state-of-the-art models, supported by extensive ablative studies. Furthermore, we visualized and examined the learned prototypes to enhance understanding. We believe our method has the potential to be a sensitive and understandable tool for sulcal pattern analysis.

Analyzing data in complicated 3D domains: Smoothing, semiparametric regression, and functional principal component analysis

In this work, we introduce a family of methods for the analysis of data observed at locations scattered in three-dimensional (3D) domains, with possibly complicated shapes. The proposed family of methods includes smoothing, regression, and functional principal component analysis for functional signals defined over (possibly nonconvex) 3D domains, appropriately complying with the nontrivial shape of the domain. This constitutes an important advance with respect to the literature, because the available methods to analyze data observed in 3D domains rely on Euclidean distances, which are inappropriate when the shape of the domain influences the phenomenon under study. The common building block of the proposed methods is a nonparametric regression model with differential regularization. We derive the asymptotic properties of the methods and show, through simulation studies, that they are superior to the available alternatives for the analysis of data in 3D domains, even when considering domains with simple shapes. We finally illustrate an application to a neurosciences study, with neuroimaging signals from functional magnetic resonance imaging, measuring neural activity in the gray matter, a nonconvex volume with a highly complicated structure.

Anatomically compliant modes of variations: New tools for brain connectivity

Anatomical complexity and data dimensionality present major issues when analysing brain connectivity data. The functional and anatomical aspects of the connections taking place in the brain are in fact equally relevant and strongly intertwined. However, due to theoretical challenges and computational issues, their relationship is often overlooked in neuroscience and clinical research. In this work, we propose to tackle this problem through Smooth Functional Principal Component Analysis, which enables to perform dimensional reduction and exploration of the variability in functional connectivity maps, complying with the formidably complicated anatomy of the grey matter volume. In particular, we analyse a population that includes controls and subjects affected by schizophrenia, starting from fMRI data acquired at rest and during a task-switching paradigm. For both sessions, we first identify the common modes of variation in the entire population. We hence explore whether the subjects' expressions along these common modes of variation differ between controls and pathological subjects. In each session, we find principal components that are significantly differently expressed in the healthy vs pathological subjects (with p-values < 0.001), highlighting clearly interpretable differences in the connectivity in the two subpopulations. For instance, the second and third principal components for the rest session capture the imbalance between the Default Mode and Executive Networks characterizing schizophrenia patients.

Classification of Alzheimer's disease based on hippocampal multivariate morphometry statistics

BACKGROUND

Alzheimer's disease (AD) is a neurodegenerative disease characterized by progressive cognitive decline, and mild cognitive impairment (MCI) is associated with a high risk of developing AD. Hippocampal morphometry analysis is believed to be the most robust magnetic resonance imaging (MRI) markers for AD and MCI. Multivariate morphometry statistics (MMS), a quantitative method of surface deformations analysis, is confirmed to have strong statistical power for evaluating hippocampus.

AIMS

We aimed to test whether surface deformation features in hippocampus can be employed for early classification of AD, MCI, and healthy controls (HC).

METHODS

We first explored the differences in hippocampus surface deformation among these three groups by using MMS analysis. Additionally, the hippocampal MMS features of selective patches and support vector machine (SVM) were used for the binary classification and triple classification.

RESULTS

By the results, we identified significant hippocampal deformation among the three groups, especially in hippocampal CA1. In addition, the binary classification of AD/HC, MCI/HC, AD/MCI showed good performances, and area under curve (AUC) of triple-classification model achieved 0.85. Finally, positive correlations were found between the hippocampus MMS features and cognitive performances.

CONCLUSIONS

The study revealed significant hippocampal deformation among AD, MCI, and HC. Additionally, we confirmed that hippocampal MMS can be used as a sensitive imaging biomarker for the early diagnosis of AD at the individual level.

A PDE-regularized smoothing method for space-time data over manifolds with application to medical data

We propose an innovative statistical-numerical method to model spatio-temporal data, observed over a generic two-dimensional Riemanian manifold. The proposed approach consists of a regression model completed with a regularizing term based on the heat equation. The model is discretized through a finite element scheme set on the manifold, and solved by resorting to a fixed point-based iterative algorithm. This choice leads to a procedure which is highly efficient when compared with a monolithic approach, and which allows us to deal with massive datasets. After a preliminary assessment on simulation study cases, we investigate the performance of the new estimation tool in practical contexts, by dealing with neuroimaging and hemodynamic data.

Toward evaluation of multiresolution cortical thickness estimation with FreeSurfer, MaCRUISE, and BrainSuite

Advances in Magnetic Resonance Imaging hardware and methodologies allow for promoting the cortical morphometry with submillimeter spatial resolution. In this paper, we generated 3D self-enhanced high-resolution (HR) MRI imaging, by adapting 1 deep learning architecture, and 3 standard pipelines, FreeSurfer, MaCRUISE, and BrainSuite, have been collectively employed to evaluate the cortical thickness. We systematically investigated the differences in cortical thickness estimation for MRI sequences at multiresolution homologously originated from the native image. It has been revealed that there systematically exhibited the preferences in determining both inner and outer cortical surfaces at higher resolution, yielding most deeper cortical surface placements toward GM/WM or GM/CSF boundaries, which directs a consistent reduction tendency of mean cortical thickness estimation; on the contrary, the lower resolution data will most probably provide a more coarse and rough evaluation in cortical surface reconstruction, resulting in a relatively thicker estimation. Although the differences of cortical thickness estimation at the diverse spatial resolution varied with one another, almost all led to roughly one-sixth to one-fifth significant reduction across the entire brain at the HR, independent to the pipelines we applied, which emphasizes on generally coherent improved accuracy in a data-independent manner and endeavors to cost-efficiency with quantitative opportunities.

Quantification of sulcal emergence timing and its variability in early fetal life: Hemispheric asymmetry and sex difference

Human fetal brains show regionally different temporal patterns of sulcal emergence following a regular timeline, which may be associated with spatiotemporal patterns of gene expression among cortical regions. This study aims to quantify the timing of sulcal emergence and its temporal variability across typically developing fetuses by fitting a logistic curve to presence or absence of sulcus. We found that the sulcal emergence started from the central to the temporo-parieto-occipital lobes and frontal lobe, and the temporal variability of emergence in most of the sulci was similar between 1 and 2 weeks. Small variability (< 1 week) was found in the left central and postcentral sulci and larger variability (>2 weeks) was shown in the bilateral occipitotemporal and left superior temporal sulci. The temporal variability showed a positive correlation with the emergence timing that may be associated with differential contributions between genetic and environmental factors. Our statistical analysis revealed that the right superior temporal sulcus emerged earlier than the left. Female fetuses showed a trend of earlier sulcal emergence in the right superior temporal sulcus, lower temporal variability in the right intraparietal sulcus, and higher variability in the right precentral sulcus compared to male fetuses. Our quantitative and statistical approach quantified the temporal patterns of sulcal emergence in detail that can be a reference for assessing the normality of developing fetal gyrification.

Cortical thickness systematically varies with curvature and depth in healthy human brains

Cortical thickness varies throughout the cortex in a systematic way. However, it is challenging to investigate the patterns of cortical thickness due to the intricate geometry of the cortex. The cortex has a folded nature both in radial and tangential directions which forms not only gyri and sulci but also tangential folds and intersections. In this article, cortical curvature and depth are used to characterize the spatial distribution of the cortical thickness with much higher resolution than conventional regional atlases. To do this, a computational pipeline was developed that is capable of calculating a variety of quantitative measures such as surface area, cortical thickness, curvature (mean curvature, Gaussian curvature, shape index, intrinsic curvature index, and folding index), and sulcal depth. By analyzing 501 neurotypical adult human subjects from the ABIDE-I dataset, we show that cortex has a very organized structure and cortical thickness is strongly correlated with local shape. Our results indicate that cortical thickness consistently increases along the gyral-sulcal spectrum from concave to convex shape, encompassing the saddle shape along the way. Additionally, tangential folds influence cortical thickness in a similar way as gyral and sulcal folds; outer folds are consistently thicker than inner.

Investigating the Effect of Tau Deposition and Apoe on Hippocampal Morphometry in Alzheimer’s Disease: A Federated Chow Test Model

Alzheimer's disease (AD) affects more than 1 in 9 people age 65 and older and becomes an urgent public health concern as the global population ages. Tau tangle is the specific protein pathological hallmark of AD and plays a crucial role in leading to dementia-related structural deformations observed in magnetic resonance imaging (MRI) scans. The volume loss of hippocampus is mainly related to the development of AD. Besides, apolipoprotein E (APOE) also has significant effects on the risk of developing AD. However, few studies focus on integrating genotypes, MRI, and tau deposition to infer multimodal relationships. In this paper, we proposed a federated chow test model to study the synergistic effects of APOE and tau on hippocampal morphometry. Our experimental results demonstrate our model can detect the difference of tau deposition and hippocampal atrophy among the cohorts with different genotypes and subiculum and cornu ammonis 1 (CA1 subfield) were identified as hippocampal subregions where atrophy is strongly associated with abnormal tau in the homozygote cohort. Our model will provide novel insight into the neural mechanisms about the individual impact of APOE and tau deposition on brain imaging.

Brain network effects by continuous theta burst stimulation in mal de débarquement syndrome: simultaneous EEG and fMRI study

Objective. Heterogeneous clinical responses to treatment with non-invasive brain stimulation are commonly observed, making it necessary to determine personally optimized stimulation parameters. We investigated neuroimaging markers of effective brain targets of treatment with continuous theta burst stimulation (cTBS) in mal de débarquement syndrome (MdDS), a balance disorder of persistent oscillating vertigo previously shown to exhibit abnormal intrinsic functional connectivity.Approach.Twenty-four right-handed, cTBS-naive individuals with MdDS received single administrations of cTBS over one of three stimulation targets in randomized order. The optimal target was determined based on the assessment of acute changes after the administration of cTBS over each target. Repetitive cTBS sessions were delivered on three consecutive days with the optimal target chosen by the participant. Electroencephalography (EEG) was recorded at single-administration test sessions of cTBS. Simultaneous EEG and functional MRI data were acquired at baseline and after completion of 10-12 sessions. Network connectivity changes after single and repetitive stimulations of cTBS were analyzed.Main results.Using electrophysiological source imaging and a data-driven method, we identified network-level connectivity changes in EEG that correlated with symptom responses after completion of multiple sessions of cTBS. We further determined that connectivity changes demonstrated by EEG during test sessions of single administrations of cTBS were signatures that could predict optimal targets.Significance.Our findings demonstrate the effect of cTBS on resting state brain networks and suggest an imaging-based, closed-loop stimulation paradigm that can identify optimal targets during short-term test sessions of stimulation.ClinicalTrials.gov Identifier:NCT02470377.

Federated Morphometry Feature Selection for Hippocampal Morphometry Associated Beta-Amyloid and Tau Pathology

Amyloid-β (Aβ) plaques and tau protein tangles in the brain are now widely recognized as the defining hallmarks of Alzheimer's disease (AD), followed by structural atrophy detectable on brain magnetic resonance imaging (MRI) scans. One of the particular neurodegenerative regions is the hippocampus to which the influence of Aβ/tau on has been one of the research focuses in the AD pathophysiological progress. This work proposes a novel framework, Federated Morphometry Feature Selection (FMFS) model, to examine subtle aspects of hippocampal morphometry that are associated with Aβ/tau burden in the brain, measured using positron emission tomography (PET). FMFS is comprised of hippocampal surface-based feature calculation, patch-based feature selection, federated group LASSO regression, federated screening rule-based stability selection, and region of interest (ROI) identification. FMFS was tested on two Alzheimer's Disease Neuroimaging Initiative (ADNI) cohorts to understand hippocampal alterations that relate to Aβ/tau depositions. Each cohort included pairs of MRI and PET for AD, mild cognitive impairment (MCI), and cognitively unimpaired (CU) subjects. Experimental results demonstrated that FMFS achieves an 89× speedup compared to other published state-of-the-art methods under five independent hypothetical institutions. In addition, the subiculum and cornu ammonis 1 (CA1 subfield) were identified as hippocampal subregions where atrophy is strongly associated with abnormal Aβ/tau. As potential biomarkers for Aβ/tau pathology, the features from the identified ROIs had greater power for predicting cognitive assessment and for survival analysis than five other imaging biomarkers. All the results indicate that FMFS is an efficient and effective tool to reveal associations between Aβ/tau burden and hippocampal morphometry.

Contributions of Sex Chromosomes and Gonadal Hormones to the Male Bias in a Maternal Antibody-Induced Model of Autism Spectrum Disorder

Autism Spectrum Disorder (ASD) is a group of neurodevelopmental conditions that is four times more commonly diagnosed in males than females. While susceptibility genes located in the sex chromosomes have been identified in ASD, it is unclear whether they are sufficient to explain the male bias or whether gonadal hormones also play a key role. We evaluated the sex chromosomal and hormonal influences on the male bias in a murine model of ASD, in which mice are exposed in utero to a maternal antibody reactive to contactin-associated protein-like 2 (Caspr2), which was originally cloned from a mother of a child with ASD (termed C6 mice henceforth). In this model, only male mice are affected. We used the four-core-genotypes (FCG) model in which the Sry gene is deleted from the Y chromosome (Y⁻) and inserted into autosome 3 (TgSry). Thus, by combining the C6 and FCG models, we were able to differentiate the contributions of sex chromosomes and gonadal hormones to the development of fetal brain and adult behavioral phenotypes. We show that the presence of the Y chromosome, or lack of two X chromosomes, irrespective of gonadal sex, increased the susceptibility to C6-induced phenotypes including the abnormal growth of the developing fetal cerebral cortex, as well as a behavioral pattern of decreased open-field exploration in adult mice. Our results indicate that sex chromosomes are the main determinant of the male bias in the maternal C6-induced model of ASD. The less dominant hormonal effect may be due to modulation by sex chromosome genes of factors involved in gonadal hormone pathways in the brain.

The "Narratives" fMRI dataset for evaluating models of naturalistic language comprehension

The "Narratives" collection aggregates a variety of functional MRI datasets collected while human subjects listened to naturalistic spoken stories. The current release includes 345 subjects, 891 functional scans, and 27 diverse stories of varying duration totaling ~4.6 hours of unique stimuli (~43,000 words). This data collection is well-suited for naturalistic neuroimaging analysis, and is intended to serve as a benchmark for models of language and narrative comprehension. We provide standardized MRI data accompanied by rich metadata, preprocessed versions of the data ready for immediate use, and the spoken story stimuli with time-stamped phoneme- and word-level transcripts. All code and data are publicly available with full provenance in keeping with current best practices in transparent and reproducible neuroimaging.

Tetrahedral spectral feature-Based bayesian manifold learning for grey matter morphometry: Findings from the Alzheimer's disease neuroimaging initiative

Structural and anatomical analyses of magnetic resonance imaging (MRI) data often require a reconstruction of the three-dimensional anatomy to a statistical shape model. Our prior work demonstrated the usefulness of tetrahedral spectral features for grey matter morphometry. However, most of the current methods provide a large number of descriptive shape features, but lack an unsupervised scheme to automatically extract a concise set of features with clear biological interpretations and that also carries strong statistical power. Here we introduce a new tetrahedral spectral feature-based Bayesian manifold learning framework for effective statistical analysis of grey matter morphology. We start by solving the technical issue of generating tetrahedral meshes which preserve the details of the grey matter geometry. We then derive explicit weak-form tetrahedral discretizations of the Hamiltonian operator (HO) and the Laplace-Beltrami operator (LBO). Next, the Schrödinger's equation is solved for constructing the scale-invariant wave kernel signature (SIWKS) as the shape descriptor. By solving the heat equation and utilizing the SIWKS, we design a morphometric Gaussian process (M-GP) regression framework and an active learning strategy to select landmarks as concrete shape descriptors. We evaluate the proposed system on publicly available data from the Alzheimers Disease Neuroimaging Initiative (ADNI), using subjects structural MRI covering the range from cognitively unimpaired (CU) to full blown Alzheimer's disease (AD). Our analyses suggest that the SIWKS and M-GP compare favorably with seven other baseline algorithms to obtain grey matter morphometry-based diagnoses. Our work may inspire more tetrahedral spectral feature-based Bayesian learning research in medical image analysis.

Automatic classification of autism spectrum disorder in children using cortical thickness and support vector machine

OBJECTIVE

Autism spectrum disorder (ASD) is a neurodevelopmental condition with a heterogeneous phenotype. The role of biomarkers in ASD diagnosis has been highlighted; cortical thickness has proved to be involved in the etiopathogenesis of ASD core symptoms. We apply support vector machine, a supervised machine learning method, in order to identify specific cortical thickness alterations in ASD subjects.

METHODS

A sample of 76 subjects (9.5 ± 3.4 years old) has been selected, 40 diagnosed with ASD and 36 typically developed subjects. All children underwent a magnetic resonance imaging (MRI) examination; T1-MPRAGE sequences were analyzed to extract features for the characterization and parcellation of regions of interests (ROI); average cortical thickness (CT) has been measured for each ROI. For the classification process, the extracted features were used as input for a classifier to identify ASD subjects through a "learning by example" procedure; the features with best performance was then selected by "greedy forward-feature selection." Finally, this model underwent a leave-one-out cross-validation approach.

RESULTS

From the training set of 68 ROIs, five ROIs reached accuracies of over 70%. After this phase, we used a recursive feature selection process in order to identify the eight features with the best accuracy (84.2%). CT resulted higher in ASD compared to controls in all the ROIs identified at the end of the process.

CONCLUSION

We found increased CT in various brain regions in ASD subjects, confirming their role in the pathogenesis of this condition. Considering the brain development curve during ages, these changes in CT may normalize during development. Further validation on a larger sample is required.

Cortical Morphology in Autism: Findings from a Cortical Shape-Adaptive Approach to Local Gyrification Indexing

It has been challenging to elucidate the differences in brain structure that underlie behavioral features of autism. Prior studies have begun to identify patterns of changes in autism across multiple structural indices, including cortical thickness, local gyrification, and sulcal depth. However, common approaches to local gyrification indexing used in prior studies have been limited by low spatial resolution relative to functional brain topography. In this study, we analyze the aforementioned structural indices, utilizing a new method of local gyrification indexing that quantifies this index adaptively in relation to specific sulci/gyri, improving interpretation with respect to functional organization. Our sample included n = 115 autistic and n = 254 neurotypical participants aged 5-54, and we investigated structural patterns by group, age, and autism-related behaviors. Differing structural patterns by group emerged in many regions, with age moderating group differences particularly in frontal and limbic regions. There were also several regions, particularly in sensory areas, in which one or more of the structural indices of interest either positively or negatively covaried with autism-related behaviors. Given the advantages of this approach, future studies may benefit from its application in hypothesis-driven examinations of specific brain regions and/or longitudinal studies to assess brain development in autism.

Attention, awareness, and the right temporoparietal junction

The attention schema theory posits a specific relationship between subjective awareness and attention, in which awareness is the control model that the brain uses to aid in the endogenous control of attention. In previous experiments, we developed a behavioral paradigm in human subjects to manipulate awareness and attention. The paradigm involved a visual cue that could be used to guide attention to a target stimulus. In task 1, subjects were aware of the cue, but not aware that it provided information about the target. The cue measurably drew exogenous attention to itself. In addition, implicitly, the subjects' endogenous attention mechanism used the cue to help shift attention to the target. In task 2, subjects were no longer aware of the cue. The cue still measurably drew exogenous attention to itself, yet without awareness of the cue, the subjects' endogenous control mechanism was no longer able to use the cue to control attention. Thus, the control of attention depended on awareness. Here, we tested the two tasks while scanning brain activity in human volunteers. We predicted that the right temporoparietal junction (TPJ) would be active in relation to the process in which awareness helps control attention. This prediction was confirmed. The right TPJ was active in relation to the effect of the cue on attention in task 1; it was not measurably active in task 2. The difference was significant. In our interpretation, the right TPJ is involved in an interaction in which awareness permits the control of attention.

Brain-wide functional diffuse optical tomography of resting state networks

Objective.Diffuse optical tomography (DOT) has the potential in reconstructing resting state networks (RSNs) in human brains with high spatio-temporal resolutions and multiple contrasts. While several RSNs have been reported and successfully reconstructed using DOT, its full potential in recovering a collective set of distributed brain-wide networks with the number of RSNs close to those reported using functional magnetic resonance imaging (fMRI) has not been demonstrated.Approach.The present study developed a novel brain-wide DOT (BW-DOT) framework that integrates a cap-based whole-head optode placement system with multiple computational approaches, i.e. finite-element modeling, inverse source reconstruction, data-driven pattern recognition, and statistical correlation tomography, to reconstruct RSNs in dual contrasts of oxygenated (HbO) and deoxygenated hemoglobins (HbR).Main results.Our results from the proposed framework revealed a comprehensive set of RSNs and their subnetworks, which collectively cover almost the entire neocortical surface of the human brain, both at the group level and individual participants. The spatial patterns of these DOT RSNs suggest statistically significant similarities to fMRI RSN templates. Our results also reported the networks involving the medial prefrontal cortex and precuneus that had been missed in previous DOT studies. Furthermore, RSNs obtained from HbO and HbR suggest similarity in terms of both the number of RSN types reconstructed and their corresponding spatial patterns, while HbR RSNs show statistically more similarity to fMRI RSN templates and HbO RSNs indicate more bilateral patterns over two hemispheres. In addition, the BW-DOT framework allowed consistent reconstructions of RSNs across individuals and across recording sessions, indicating its high robustness and reproducibility, respectively.Significance.Our present results suggest the feasibility of using the BW-DOT, as a neuroimaging tool, in simultaneously mapping multiple RSNs and its potential values in studying RSNs, particularly in patient populations under diverse conditions and needs, due to its advantages in accessibility over fMRI.

Predicting future cognitive decline with hyperbolic stochastic coding

Hyperbolic geometry has been successfully applied in modeling brain cortical and subcortical surfaces with general topological structures. However, such approaches, similar to other surface-based brain morphology analysis methods, usually generate high dimensional features. It limits their statistical power in cognitive decline prediction research, especially in datasets with limited subject numbers. To address the above limitation, we propose a novel framework termed as hyperbolic stochastic coding (HSC). We first compute diffeomorphic maps between general topological surfaces by mapping them to a canonical hyperbolic parameter space with consistent boundary conditions and extracts critical shape features. Secondly, in the hyperbolic parameter space, we introduce a farthest point sampling with breadth-first search method to obtain ring-shaped patches. Thirdly, stochastic coordinate coding and max-pooling algorithms are adopted for feature dimension reduction. We further validate the proposed system by comparing its classification accuracy with some other methods on two brain imaging datasets for Alzheimer's disease (AD) progression studies. Our preliminary experimental results show that our algorithm achieves superior results on various classification tasks. Our work may enrich surface-based brain imaging research tools and potentially result in a diagnostic and prognostic indicator to be useful in individualized treatment strategies.

Anatomical bases of fast parietal grasp control in humans: A diffusion-MRI tractography study

The dorso-posterior parietal cortex (DPPC) is a major node of the grasp/manipulation control network. It is assumed to act as an optimal forward estimator that continuously integrates efferent outflows and afferent inflows to modulate the ongoing motor command. In agreement with this view, a recent per-operative study, in humans, identified functional sites within DPPC that: (i) instantly disrupt hand movements when electrically stimulated; (ii) receive short-latency somatosensory afferences from intrinsic hand muscles. Based on these results, it was speculated that DPPC is part of a rapid grasp control loop that receives direct inputs from the hand-territory of the primary somatosensory cortex (S1) and sends direct projections to the hand-territory of the primary motor cortex (M1). However, evidence supporting this hypothesis is weak and partial. To date, projections from DPPC to M1 grasp zone have been identified in monkeys and have been postulated to exist in humans based on clinical and transcranial magnetic studies. This work uses diffusion-MRI tractography in two samples of right- (n = 50) and left-handed (n = 25) subjects randomly selected from the Human Connectome Project. It aims to determine whether direct connections exist between DPPC and the hand control sectors of the primary sensorimotor regions. The parietal region of interest, related to hand control (hereafter designated DPPC_hand), was defined permissively as the 95% confidence area of the parietal sites that were found to disrupt hand movements in the previously evoked per-operative study. In both hemispheres, irrespective of handedness, we found dense ipsilateral connections between a restricted part of DPPC_hand and focal sectors within the pre and postcentral gyrus. These sectors, corresponding to the hand territories of M1 and S1, targeted the same parietal zone (spatial overlap > 92%). As a sensitivity control, we searched for potential connections between the angular gyrus (AG) and the pre and postcentral regions. No robust pathways were found. Streamline densities identified using AG as the starting seed represented less than 5 % of the streamline densities identified from DPPC_hand. Together, these results support the existence of a direct sensory-parietal-motor loop suited for fast manual control and more generally, for any task requiring rapid integration of distal sensorimotor signals.

Intergenerational Transmission of Cortical Sulcal Patterns from Mothers to their Children

Intergenerational effects are described as the genetic, epigenetic, as well as pre- and postnatal environmental influence parents have on their offspring's behavior, cognition, and brain. During fetal brain development, the primary cortical sulci emerge with a distinctive folding pattern that are under strong genetic influence and show little change of this pattern throughout postnatal brain development. We examined intergenerational transmission of cortical sulcal patterns by comparing primary sulcal patterns between children (N = 16, age 5.5 ± 0.81 years, 8 males) and their biological mothers (N = 15, age 39.72 ± 4.68 years) as well as between children and unrelated adult females. Our graph-based sulcal pattern comparison method detected stronger sulcal pattern similarity for child-mother pairs than child-unrelated pairs, where higher similarity between child-mother pairs was observed mostly for the right lobar regions. Our results also show that child-mother versus child-unrelated pairs differ for daughters and sons with a trend toward significance, particularly for the left hemisphere lobar regions. This is the first study to reveal significant intergenerational transmission of cortical sulcal patterns, and our results have important implications for the study of the heritability of complex behaviors, brain-based disorders, the identification of biomarkers, and targets for interventions.

Zero-crossing patterns reveal subtle epileptiform discharges in the scalp EEG

Clinical diagnosis of epilepsy depends heavily on the detection of interictal epileptiform discharges (IEDs) from scalp electroencephalographic (EEG) signals, which by purely visual means is far from straightforward. Here, we introduce a simple signal analysis procedure based on scalp EEG zero-crossing patterns which can extract the spatiotemporal structure of scalp voltage fluctuations. We analyzed simultaneous scalp and intracranial EEG recordings from patients with pharmacoresistant temporal lobe epilepsy. Our data show that a large proportion of intracranial IEDs manifest only as subtle, low-amplitude waveforms below scalp EEG background and could, therefore, not be detected visually. We found that scalp zero-crossing patterns allow detection of these intracranial IEDs on a single-trial level with millisecond temporal precision and including some mesial temporal discharges that do not propagate to the neocortex. Applied to an independent dataset, our method discriminated accurately between patients with epilepsy and normal subjects, confirming its practical applicability.

A gyrification analysis approach based on Laplace Beltrami eigenfunction level sets

An accurate measure of the complexity of patterns of cortical folding or gyrification is necessary for understanding normal brain development and neurodevelopmental disorders. Conventional gyrification indices (GIs) are calculated based on surface curvature (curvature-based GI) or an outer hull surface of the cortex (outer surface-based GI). The latter is dependent on the definition of the outer hull surface and a corresponding function between surfaces. In the present study, we propose the Laplace Beltrami-based gyrification index (LB-GI). This is a new curvature-based local GI computed using the first three Laplace Beltrami eigenfunction level sets. As with outer surface-based GI methods, this method is based on the hypothesis that gyrification stems from a flat surface during development. However, instead of quantifying gyrification with reference to corresponding points on an outer hull surface, LB-GI quantifies the gyrification at each point on the cortical surface with reference to their surrounding gyral points, overcoming several shortcomings of existing methods. The LB-GI was applied to investigate the cortical maturation profile of the human brain from preschool to early adulthood using the PING database. The results revealed more detail in patterns of cortical folding than conventional curvature-based methods, especially on frontal and posterior tips of the brain, such as the frontal pole, lateral occipital, lateral cuneus, and lingual. Negative associations of cortical folding with age were observed at cortical regions, including bilateral lingual, lateral occipital, precentral gyrus, postcentral gyrus, and superior frontal gyrus. The results also indicated positive significant associations between age and the LB-GI of bilateral insula, the medial orbitofrontal, frontal pole and rostral anterior cingulate regions. It is anticipated that the LB-GI will be advantageous in providing further insights in the understanding of brain development and degeneration in large clinical neuroimaging studies.

Apolipoprotein E4 effects on topological brain network organization in mild cognitive impairment

The Apolipoprotein E isoform E4 (ApoE4) is consistently associated with an elevated risk of developing late-onset Alzheimer's Disease (AD); however, less is known about the potential genetic modulation of the brain networks organization during prodromal stages like Mild Cognitive Impairment (MCI). To investigate this issue during this critical stage, we used a dataset with a cross-sectional sample of 253 MCI patients divided into ApoE4-positive (‛Carriers') and ApoE4-negative ('non-Carriers'). We estimated the cortical thickness (CT) from high-resolution T1-weighted structural magnetic images to calculate the correlation among anatomical regions across subjects and build the CT covariance networks (CT-Nets). The topological properties of CT-Nets were described through the graph theory approach. Specifically, our results showed a significant decrease in characteristic path length, clustering-index, local efficiency, global connectivity, modularity, and increased global efficiency for Carriers compared to non-Carriers. Overall, we found that ApoE4 in MCI shaped the topological organization of CT-Nets. Our results suggest that in the MCI stage, the ApoE4 disrupting the CT correlation between regions may be due to adaptive mechanisms to sustain the information transmission across distant brain regions to maintain the cognitive and behavioral abilities before the occurrence of the most severe symptoms.

Cortical Morphological Changes in Congenital Amusia: Surface-Based Analyses

Background: Congenital amusia (CA) is a rare disorder characterized by deficits in pitch perception, and many structural and functional magnetic resonance imaging studies have been conducted to better understand its neural bases. However, a structural magnetic resonance imaging analysis using a surface-based morphology method to identify regions with cortical features abnormalities at the vertex-based level has not yet been performed. Methods: Fifteen participants with CA and 13 healthy controls underwent structural magnetic resonance imaging. A surface-based morphology method was used to identify anatomical abnormalities. Then, the surface parameters' mean value of the identified clusters with statistically significant between-group differences were extracted and compared. Finally, Pearson's correlation analysis was used to assess the correlation between the Montreal Battery of Evaluation of Amusia (MBEA) scores and surface parameters. Results: The CA group had significantly lower MBEA scores than the healthy controls (p = 0.000). The CA group exhibited a significant higher fractal dimension in the right caudal middle frontal gyrus and a lower sulcal depth in the right pars triangularis gyrus (p < 0.05; false discovery rate-corrected at the cluster level) compared to healthy controls. There were negative correlations between the mean fractal dimension values in the right caudal middle frontal gyrus and MBEA score, including the mean MBEA score (r = -0.5398, p = 0.0030), scale score (r = -0.5712, p = 0.0015), contour score (r = -0.4662, p = 0.0124), interval score (r = -0.4564, p = 0.0146), rhythmic score (r = -0.5133, p = 0.0052), meter score (r = -0.3937, p = 0.0382), and memory score (r = -0.3879, p = 0.0414). There was a significant positive correlation between the mean sulcal depth in the right pars triangularis gyrus and the MBEA score, including the mean score (r = 0.5130, p = 0.0052), scale score (r = 0.5328, p = 0.0035), interval score (r = 0.4059, p = 0.0321), rhythmic score (r = 0.5733, p = 0.0014), meter score (r = 0.5061, p = 0.0060), and memory score (r = 0.4001, p = 0.0349). Conclusion: Individuals with CA exhibit cortical morphological changes in the right hemisphere. These findings may indicate that the neural basis of speech perception and memory impairments in individuals with CA is associated with abnormalities in the right pars triangularis gyrus and middle frontal gyrus, and that these cortical abnormalities may be a neural marker of CA.

Developing univariate neurodegeneration biomarkers with low-rank and sparse subspace decomposition

Cognitive decline due to Alzheimer's disease (AD) is closely associated with brain structure alterations captured by structural magnetic resonance imaging (sMRI). It supports the validity to develop sMRI-based univariate neurodegeneration biomarkers (UNB). However, existing UNB work either fails to model large group variances or does not capture AD dementia (ADD) induced changes. We propose a novel low-rank and sparse subspace decomposition method capable of stably quantifying the morphological changes induced by ADD. Specifically, we propose a numerically efficient rank minimization mechanism to extract group common structure and impose regularization constraints to encode the original 3D morphometry connectivity. Further, we generate regions-of-interest (ROI) with group difference study between common subspaces of Aβ+AD and Aβ-cognitively unimpaired (CU) groups. A univariate morphometry index (UMI) is constructed from these ROIs by summarizing individual morphological characteristics weighted by normalized difference between Aβ+AD and Aβ-CU groups. We use hippocampal surface radial distance feature to compute the UMIs and validate our work in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. With hippocampal UMIs, the estimated minimum sample sizes needed to detect a 25% reduction in the mean annual change with 80% power and two-tailed P=0.05are 116, 279 and 387 for the longitudinal Aβ+AD, Aβ+mild cognitive impairment (MCI) and Aβ+CU groups, respectively. Additionally, for MCI patients, UMIs well correlate with hazard ratio of conversion to AD (4.3, 95% CI = 2.3-8.2) within 18 months. Our experimental results outperform traditional hippocampal volume measures and suggest the application of UMI as a potential UNB.

Structural abnormalities of cingulate cortex in patients with first-episode drug-naïve schizophrenia comorbid with depressive symptoms

Depressive symptoms are common in patients with first-episode psychosis. However, the neural mechanisms underlying the comorbid depression in schizophrenia are still unknown. The main purpose of this study was to characterize the structural abnormalities of first-episodes drug-naïve (FEDN) schizophrenia comorbid with depression by utilizing both volume-based and surface-based morphometric measurements. Forty-two patients with FEDN schizophrenia and 29 healthy controls were recruited. The 24-item Hamilton Depression Rating Scale (HAMD-24) was administrated to divide all patients into depressive patients (DP) and non-depressive patients (NDP). Compared with NDP, DP had a significantly larger volume and surface area in the left isthmus cingulate cortex and also had a greater volume in the left posterior cingulate cortex. Correlation analysis showed that HAMD total score was positively correlated with the surface area of the left isthmus cingulate and gray matter volume of the left isthmus cingulate cortex. In addition, gray matter volume of the left isthmus cingulate was also correlated with the PANSS general psychopathology or total score. The findings suggest that prominent structural abnormalities of gray matter are mainly concentrated on the cingulate cortex in FEDN schizophrenia patients comorbid with depression, which may contribute to depressive symptoms and psychopathological symptoms.

Patch-Based Surface Morphometry Feature Selection with Federated Group Lasso Regression

Collectively, vast quantities of brain imaging data exist across hospitals and research institutions, providing valuable resources to study brain disorders such as Alzheimer's disease (AD). However, in practice, putting all these distributed datasets into a centralized platform is infeasible due to patient privacy concerns, data restrictions and legal regulations. In this study, we propose a novel federated feature selection framework that can analyze the data at each individual institution without data-sharing or accessing private patient information. In this framework, we first propose a federated group lasso optimization method based on block coordinate descent. We employ stability selection to determine statistically significant features, by solving the group lasso problem with a sequence of regularization parameters. To accelerate the stability selection, we further propose a federated screening rule, which can identify and exclude the irrelevant features before solving the group lasso. Here, we use this framework for patch based feature selection on hippocampal morphometry. Shape is characterized through two different kinds of local measures, the radial distance and the surface area determined via tensor-based morphometry (TBM). The method is tested on 1,127 T1-weighted brain magnetic resonance images (MRI) of AD, mild cognitive impairment (MCI) and elderly control subjects, randomly assigned to five independent hypothetical institutions for testing purpose. We examine the association of MRI-based anatomical measures with general cognitive assessment and amyloid burden to identify the morphometry changes related to AD deterioration and plaque accumulation. Finally, we visualize the significance of the association on the hippocampal surfaces. Our experimental results successfully demonstrate the efficiency and effectiveness of our method.

A Novel Method for High-Dimensional Anatomical Mapping of Extra-Axial Cerebrospinal Fluid: Application to the Infant Brain

Cerebrospinal fluid (CSF) plays an essential role in early postnatal brain development. Extra-axial CSF (EA-CSF) volume, which is characterized by CSF in the subarachnoid space surrounding the brain, is a promising marker in the early detection of young children at risk for neurodevelopmental disorders. Previous studies have focused on global EA-CSF volume across the entire dorsal extent of the brain, and not regionally-specific EA-CSF measurements, because no tools were previously available for extracting local EA-CSF measures suitable for localized cortical surface analysis. In this paper, we propose a novel framework for the localized, cortical surface-based analysis of EA-CSF. The proposed processing framework combines probabilistic brain tissue segmentation, cortical surface reconstruction, and streamline-based local EA-CSF quantification. The quantitative analysis of local EA-CSF was applied to a dataset of typically developing infants with longitudinal MRI scans from 6 to 24 months of age. There was a high degree of consistency in the spatial patterns of local EA-CSF across age using the proposed methods. Statistical analysis of local EA-CSF revealed several novel findings: several regions of the cerebral cortex showed reductions in EA-CSF from 6 to 24 months of age, and specific regions showed higher local EA-CSF in males compared to females. These age-, sex-, and anatomically-specific patterns of local EA-CSF would not have been observed if only a global EA-CSF measure were utilized. The proposed methods are integrated into a freely available, open-source, cross-platform, user-friendly software tool, allowing neuroimaging labs to quantify local extra-axial CSF in their neuroimaging studies to investigate its role in typical and atypical brain development.

Altered White Matter Organization in the TUBB3 E410K Syndrome

Seven unrelated individuals (four pediatric, three adults) with the TUBB3 E410K syndrome, harboring identical de novo heterozygous TUBB3 c.1228 G>A mutations, underwent neuropsychological testing and neuroimaging. Despite the absence of cortical malformations, they have intellectual and social disabilities. To search for potential etiologies for these deficits, we compared their brain's structural and white matter organization to 22 controls using structural and diffusion magnetic resonance imaging. Diffusion images were processed to calculate fractional anisotropy (FA) and perform tract reconstructions. Cortical parcellation-based network analysis and gyral topology-based FA analyses were performed. Major interhemispheric, projection and intrahemispheric tracts were manually segmented. Subjects had decreased corpus callosum volume and decreased network efficiency. While only pediatric subjects had diffuse decreases in FA predominantly affecting mid- and long-range tracts, only adult subjects had white matter volume loss associated with decreased cortical surface area. All subjects showed aberrant corticospinal tract trajectory and bilateral absence of the dorsal language network long segment. Furthermore, pediatric subjects had more tracts with decreased FA compared with controls than did adult subjects. These findings define a TUBB3 E410K neuroimaging endophenotype and lead to the hypothesis that the age-related changes are due to microscopic intrahemispheric misguided axons that are pruned during maturation.

Early-Emerging Sulcal Patterns Are Atypical in Fetuses with Congenital Heart Disease

Fetuses with congenital heart disease (CHD) have third trimester alterations in cortical development on brain magnetic resonance imaging (MRI). However, the intersulcal relationships contributing to global sulcal pattern remain unknown. This study applied a novel method for examining the geometric and topological relationships between sulci to fetal brain MRIs from 21-30 gestational weeks in CHD fetuses (n = 19) and typically developing (TD) fetuses (n = 17). Sulcal pattern similarity index (SI) to template fetal brain MRIs was determined for the position, area, and depth for corresponding sulcal basins and intersulcal relationships for each subject. CHD fetuses demonstrated altered global sulcal patterns in the left hemisphere compared with TD fetuses (TD [SI, mean ± SD]: 0.822 ± 0.023, CHD: 0.795 ± 0.030, P = 0.002). These differences were present in the earliest emerging sulci and were driven by differences in the position of corresponding sulcal basins (TD: 0.897 ± 0.024, CHD: 0.878 ± 0.019, P = 0.006) and intersulcal relationships (TD: 0.876 ± 0.031, CHD: 0.857 ± 0.018, P = 0.033). No differences in cortical gyrification index, mean curvature, or surface area were present. These data suggest our methods may be more sensitive than traditional measures for evaluating cortical developmental alterations early in gestation.

Atypical age-related changes in cortical thickness in autism spectrum disorder

Recent longitudinal neuroimaging and neurophysiological studies have shown that tracking relative age-related changes in neural signals, rather than a static snapshot of a neural measure, could offer higher sensitivity for discriminating typically developing (TD) individuals from those with autism spectrum disorder (ASD). It is not clear, however, which aspects of age-related changes (trajectories) would be optimal for identifying atypical brain development in ASD. Using a large cross-sectional data set (Autism Brain Imaging Data Exchange [ABIDE] repository; releases I and II), we aimed to explore age-related changes in cortical thickness (CT) in TD and ASD populations (age range 6–30 years old). Cortical thickness was estimated from T1-weighted MRI images at three scales of spatial coarseness (three parcellations with different numbers of regions of interest). For each parcellation, three polynomial models of age-related changes in CT were tested. Specifically, to characterize alterations in CT trajectories, we compared the linear slope, curvature, and aberrancy of CT trajectories across experimental groups, which was estimated using linear, quadratic, and cubic polynomial models, respectively. Also, we explored associations between age-related changes with ASD symptomatology quantified as the Autism Diagnostic Observation Schedule (ADOS) scores. While no overall group differences in cortical thickness were observed across the entire age range, ASD and TD populations were different in terms of age-related changes, which were located primarily in frontal and tempo-parietal areas. These atypical age-related changes were also associated with ADOS scores in the ASD group and used to predict ASD from TD development. These results indicate that the curvature is the most reliable feature for localizing brain areas developmentally atypical in ASD with a more pronounced effect with symptomatology and is the most sensitive in predicting ASD development.

Applying surface-based morphometry to study ventricular abnormalities of cognitively unimpaired subjects prior to clinically significant memory decline

•

A completely automated surface-based ventricular morphometry system.

•

Generate a whole connected 3D ventricular shape model.

•

Test-retest the system in two independent CU subject cohorts.

•

Subregional ventricular abnormalities prior to clinically memory decline.

Ventricular volume (VV) is a widely used structural magnetic resonance imaging (MRI) biomarker in Alzheimer’s disease (AD) research. Abnormal enlargements of VV can be detected before clinically significant memory decline. However, VV does not pinpoint the details of subregional ventricular expansions. Here we introduce a ventricular morphometry analysis system (VMAS) that generates a whole connected 3D ventricular shape model and encodes a great deal of ventricular surface deformation information that is inaccessible by VV. VMAS contains an automated segmentation approach and surface-based multivariate morphometry statistics. We applied VMAS to two independent datasets of cognitively unimpaired (CU) groups. To our knowledge, it is the first work to detect ventricular abnormalities that distinguish normal aging subjects from those who imminently progress to clinically significant memory decline. Significant bilateral ventricular morphometric differences were first shown in 38 members of the Arizona APOE cohort, which included 18 CU participants subsequently progressing to the clinically significant memory decline within 2 years after baseline visits (progressors), and 20 matched CU participants with at least 4 years of post-baseline cognitive stability (non-progressors). VMAS also detected significant differences in bilateral ventricular morphometry in 44 Alzheimer’s Disease Neuroimaging Initiative (ADNI) subjects (18 CU progressors vs. 26 CU non-progressors) with the same inclusion criterion. Experimental results demonstrated that the ventricular anterior horn regions were affected bilaterally in CU progressors, and more so on the left. VMAS may track disease progression at subregional levels and measure the effects of pharmacological intervention at a preclinical stage.

Cortical Microstructural Alterations in Mild Cognitive Impairment and Alzheimer's Disease Dementia

In Alzheimer's disease (AD), neurodegenerative processes are ongoing for years prior to the time that cortical atrophy can be reliably detected using conventional neuroimaging techniques. Recent advances in diffusion-weighted imaging have provided new techniques to study neural microstructure, which may provide additional information regarding neurodegeneration. In this study, we used neurite orientation dispersion and density imaging (NODDI), a multi-compartment diffusion model, in order to investigate cortical microstructure along the clinical continuum of mild cognitive impairment (MCI) and AD dementia. Using gray matter-based spatial statistics (GBSS), we demonstrated that neurite density index (NDI) was significantly lower throughout temporal and parietal cortical regions in MCI, while both NDI and orientation dispersion index (ODI) were lower throughout parietal, temporal, and frontal regions in AD dementia. In follow-up ROI analyses comparing microstructure and cortical thickness (derived from T1-weighted MRI) within the same brain regions, differences in NODDI metrics remained, even after controlling for cortical thickness. Moreover, for participants with MCI, gray matter NDI-but not cortical thickness-was lower in temporal, parietal, and posterior cingulate regions. Taken together, our results highlight the utility of NODDI metrics in detecting cortical microstructural degeneration that occurs prior to measurable macrostructural changes and overt clinical dementia.

An Atypical Sulcal Pattern in Children with Disorders of the Corpus Callosum and Its Relation to Behavioral Outcomes

Hypogenesis (hCC) and dysgenesis (dCC) of the corpus callosum (CC) are characterized by its smaller size or absence. The outcomes of these patients vary considerably and are unrelated to the size of the CC abnormality. The aim of the current study was to characterize the sulcal pattern in children with hCC and dCC and to explore its relation to clinical outcome. We used quantitative sulcal pattern analysis that measures deviation (similarity index, SI) of the composite or individual sulcal features (position, depth, area, and graph topology) compared to the control group. We calculated SI for each hemisphere and lobe in 11 children with CC disorder (hCC = 4, dCC = 7) and 15 controls. hCC and dCC had smaller hemispheric SI compared to controls. dCC subjects had smaller regional SI in the frontal and occipital lobes, which were driven by a smaller SI in a position or a graph topology. The significantly decreased SI gradient was found across groups only in the sulcal graph topology of the temporal lobes (controls > hCC > dCC) and was related to clinical outcome. Our results suggest that careful examination of sulcal pattern in hCC and dCC patients could be a useful biomarker of outcome.

Whole-brain electrophysiological functional connectivity dynamics in resting-state EEG

OBJECTIVE

Functional connectivity (FC) dynamics have been studied in functional magnetic resonance imaging (fMRI) data, while it is largely unknown in electrophysiological data, e.g. EEG.

APPROACH

The present study proposed a novel analytic framework to study spatiotemporal dynamics of FC (dFC) in resting-state human EEG data, including independent component analysis, cortical source imaging, sliding-window correlation analysis, and k-means clustering.

MAIN RESULTS

Our results confirm that major fMRI intrinsic connectivity networks (ICNs) can be successfully reconstructed from EEG using our analytic framework. Prominent spatial and temporal variability were revealed in these ICNs. The mean dFC spatial patterns of individual ICNs resemble their corresponding static FC (sFC) patterns but show fewer cross-talks among distinct ICNs. Our investigation unveils evidences of time-domain variations in individual ICNs comparable to their mean FC level in terms of magnitude. The major contributors to these variations are from the frequency below 0.0156 Hz, in the similar range of FC dynamics from fMRI data. Among different ICNs, larger temporal variabilities are observed in the frontal attention and auditory/visual ICNs, while sensorimotor, salience, and default model networks showed less. Our analytic framework for the first time revealed quasi-stable states within individual EEG ICNs, with various strengths or spatial patterns that were reliably detected at both group and individual levels. These states all together reveal a more complete picture of EEG ICNs: (1) quasi-stable state spatial patterns as a whole for each EEG ICN are more consistent with the corresponding fMRI ICN in terms of the bilateral distribution and multi-nodes structure; (2) EEG ICNs reveal more transient patterns about within-ICN between-node communications than fMRI ICNs.

SIGNIFICANCE

The present findings highlight the fact that rich temporal and spatial dynamics exist in ICN that can be detected from EEG data. Future studies might extend investigations towards spectral dynamics of EEG ICNs.

ENRICHING STATISTICAL INFERENCES ON BRAIN CONNECTIVITY FOR ALZHEIMER'S DISEASE ANALYSIS VIA LATENT SPACE GRAPH EMBEDDING

We develop a graph node embedding Deep Neural Network that leverages statistical outcome measure and graph structure given in the data. The objective is to identify regions of interests (ROIs) in the brain that are affected by topological changes of brain connectivity due to specific neurodegenerative diseases by enriching statistical group analysis. We tackle this problem by learning a latent space where statistical inference can be made more effectively. Our experiments on a large-scale Alzheimer's Disease dataset show promising result identifying ROIs that show statistically significant group differences separating even early and late Mild Cognitive Impairment (MCI) groups whose effect sizes are very subtle.

Integrating Convolutional Neural Networks and Multi-Task Dictionary Learning for Cognitive Decline Prediction with Longitudinal Images

BACKGROUND

Disease progression prediction based on neuroimaging biomarkers is vital in Alzheimer's disease (AD) research. Convolutional neural networks (CNN) have been proved to be powerful for various computer vision research by refining reliable and high-level feature maps from image patches.

OBJECTIVE

A key challenge in applying CNN to neuroimaging research is the limited labeled samples with high dimensional features. Another challenge is how to improve the prediction accuracy by joint analysis of multiple data sources (i.e., multiple time points or multiple biomarkers). To address these two challenges, we propose a novel multi-task learning framework based on CNN.

METHODS

First, we pre-trained CNN on the ImageNet dataset and transferred the knowledge from the pre-trained model to neuroimaging representation. We used this deep model as feature extractor to generate high-level feature maps of different tasks. Then a novel unsupervised learning method, termed Multi-task Stochastic Coordinate Coding (MSCC), was proposed for learning sparse features of multi-task feature maps by using shared and individual dictionaries. Finally, Lasso regression was performed on these multi-task sparse features to predict AD progression measured by the Mini-Mental State Examination (MMSE) and the Alzheimer's Disease Assessment Scale cognitive subscale (ADAS-Cog).

RESULTS

We applied this novel CNN-MSCC system on the Alzheimer's Disease Neuroimaging Initiative dataset to predict future MMSE/ADAS-Cog scales. We found our method achieved superior performances compared with seven other methods.

CONCLUSION

Our work may add new insights into data augmentation and multi-task deep model research and facilitate the adoption of deep models in neuroimaging research.

Neural correlates of binocular depth inversion illusion in antipsychotic-naïve first-episode schizophrenia patients

OBJECTIVES

Binocular depth inversion illusion (BDII), a visual, 'top-down'-driven information process, is impaired in schizophrenia and particularly in its early stages. BDII is a sensitive measure of impaired visual information processing and represents a valid diagnostic tool for schizophrenia and other psychotic disorders. However, neurobiological underpinnings of aberrant BDII in first-episode schizophrenia are largely unknown at present.

METHODS

In this study, 22 right-handed, first-episode, antipsychotic-naïve schizophrenia patients underwent BDII assessment and MRI scanning at 1.5 T. The surface-based analysis via new version of Freesurfer (6.0) enabled calculation of cortical thickness and surface area. BDII total and faces scores were related to the two distinct cortical measurements.

RESULTS

We found a significant correlation between BDII performance and cortical thickness in the inferior frontal gyrus and middle temporal gyrus (p < 0.003, Bonferroni corr.), as well as superior parietal gyrus, postcentral gyrus, supramarginal gyrus, and precentral gyrus (p < 0.05, CWP corr.), respectively. BDII performance was significantly correlated with surface area in the superior parietal gyrus and right postcentral gyrus (p < 0.003, Bonferroni corr.).

CONCLUSION

BDII performance may be linked to cortical thickness and surface area variations in regions involved in "adaptive" or "top-down" modulation and stimulus processing, i.e., frontal and parietal lobes. Our results suggest that cortical features of distinct evolutionary and genetic origin differently contribute to BDII performance in first-episode, antipsychotic-naïve schizophrenia patients.

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.

Altered Cortical Morphometry in Pre-manifest Huntington's Disease: Cross-sectional Data from the IMAGE-HD Study

Huntington's disease (HD) is an inherited progressive neurodegenerative disease mainly associated with subcortical striatal atrophy. There is also strong evidence showing cerebral atrophy and cortical thinning; however, limited research has investigated altered patterns of cortical folding in this disease. Here, we investigated cortical morphometry via both gyrification index (GI, a measure of cortical folding) and cortical thinning. The localized GI was examined using a novel GI, namely LB-GI. As part of a cross-sectional study, pre-manifest (pre-HD) individuals (n = 29) and matched controls (n = 29) underwent T1-MRI using data from the IMAGE-HD study. Compared to controls, pre-HD individuals demonstrated significantly lower GI in the left superior parietal and the right superior temporal regions and greater cortical thinning in the bilateral pre-central and the superior frontal gyri and left caudal middle frontal gyrus, as well as the superior parietal region. For the first time, we report evidence of abnormal localized cortical folding in pre-HD. We also provide evidence that cortical folding impacts different regions of the cortical surface more so than cortical thickness. As a result, we propose a potential new biological marker that may increase our understanding of the neuropathology of HD. Greater understanding of brain changes could inform new therapeutic approaches and target points for clinical trials.

Large-scale brain modes reorganize between infant sleep states and carry prognostic information for preterms

Sleep architecture carries vital information about brain health across the lifespan. In particular, the ability to express distinct vigilance states is a key physiological marker of neurological wellbeing in the newborn infant although systems-level mechanisms remain elusive. Here, we demonstrate that the transition from quiet to active sleep in newborn infants is marked by a substantial reorganization of large-scale cortical activity and functional brain networks. This reorganization is attenuated in preterm infants and predicts visual performance at two years. We find a striking match between these empirical effects and a computational model of large-scale brain states which uncovers fundamental biophysical mechanisms not evident from inspection of the data. Active sleep is defined by reduced energy in a uniform mode of neural activity and increased energy in two more complex anteroposterior modes. Preterm-born infants show a deficit in this sleep-related reorganization of modal energy that carries novel prognostic information.

Local Cortical Gyrification is Increased in Children With Autism Spectrum Disorders, but Decreases Rapidly in Adolescents

Extensive MRI evidence indicates early brain overgrowth in autism spectrum disorders (ASDs). Local gyrification may reflect the distribution and timing of aberrant cortical expansion in ASDs. We examined MRI data from (Study 1) 64 individuals with ASD and 64 typically developing (TD) controls (7-19 years), and from (Study 2) an independent sample from the Autism Brain Imaging Data Exchange (n = 31/group). Local Gyrification Index (lGI), cortical thickness (CT), and surface area (SA) were measured. In Study 1, differences in lGI (ASD > TD) were found in left parietal and temporal and right frontal and temporal regions. lGI decreased bilaterally with age, but more steeply in ASD in left precentral, right lateral occipital, and middle frontal clusters. CT differed between groups in right perisylvian cortex (TD > ASD), but no differences were found for SA. Partial correlations between lGI and CT were generally negative, but associations were weaker in ASD in several clusters. Study 2 results were consistent, though less extensive. Altered gyrification may reflect unique information about the trajectory of cortical development in ASDs. While early overgrowth tends to be undetectable in later childhood in ASDs, findings may indicate that a trace of this developmental abnormality could remain in a disorder-specific pattern of gyrification.

Altered White Matter Connectivity Associated with Intergyral Brain Disorganization in Hemiplegic Cerebral Palsy

Despite extensive literature showing damages in the sensorimotor projection fibers of children with hemiplegic cerebral palsy (HCP), little is known about how these damages affect the global brain network. In this study, we assess the relationship between the structural integrity of sensorimotor projection fibers and the integrity of intergyral association white matter connections in children with HCP. Diffusion tensor imaging was performed in 10 children with HCP and 16 typically developing children. We estimated the regional and global white-matter connectivity using a region-of-interest (ROI)-based approach and a whole-brain gyrus-based parcellation method. Using the ROI-based approach, we tracked the spinothalamic (STh), thalamocortical (ThC), corticospinal (CST), and sensorimotor U- (SMU) fibers. Using the whole-brain parcellation method, we tracked the short-, middle-, and long-range association fibers. We observed for the more affected hemisphere of children with HCP: (i) an increase in axial diffusivity (AD), mean diffusivity (MD), and radial diffusivity (RD) for the STh and ThC fibers; (ii) a decrease in fractional anisotropy (FA) and an increase in MD and RD for the CST and SMU fibers; in (iii) a decrease in FA and an increase in AD, MD, and RD for the middle- and long-range association fibers; and (iv) an association between the integrity of sensorimotor projection and intergyral association fibers. Our findings indicate that altered structural integrity of the sensorimotor projection fibers disorganizes the intergyral association white matter connections among local and distant regions in children with HCP.

Intracortical smoothing of small-voxel fMRI data can provide increased detection power without spatial resolution losses compared to conventional large-voxel fMRI data

Continued improvement in MRI acquisition technology has made functional MRI (fMRI) with small isotropic voxel sizes down to 1 mm and below more commonly available. Although many conventional fMRI studies seek to investigate regional patterns of cortical activation for which conventional voxel sizes of 3 mm and larger provide sufficient spatial resolution, smaller voxels can help avoid contamination from adjacent white matter (WM) and cerebrospinal fluid (CSF), and thereby increase the specificity of fMRI to signal changes within the gray matter. Unfortunately, temporal signal-to-noise ratio (tSNR), a metric of fMRI sensitivity, is reduced in high-resolution acquisitions, which offsets the benefits of small voxels. Here we introduce a framework that combines small, isotropic fMRI voxels acquired at 7 T field strength with a novel anatomically-informed, surface mesh-navigated spatial smoothing that can provide both higher detection power and higher resolution than conventional voxel sizes. Our smoothing approach uses a family of intracortical surface meshes and allows for kernels of various shapes and sizes, including curved 3D kernels that adapt to and track the cortical folding pattern. Our goal is to restrict smoothing to the cortical gray matter ribbon and avoid noise contamination from CSF and signal dilution from WM via partial volume effects. We found that the intracortical kernel that maximizes tSNR does not maximize percent signal change (ΔS/S), and therefore the kernel configuration that optimizes detection power cannot be determined from tSNR considerations alone. However, several kernel configurations provided a favorable balance between boosting tSNR and ΔS/S, and allowed a 1.1-mm isotropic fMRI acquisition to have higher performance after smoothing (in terms of both detection power and spatial resolution) compared to an unsmoothed 3.0-mm isotropic fMRI acquisition. Overall, the results of this study support the strategy of acquiring voxels smaller than the cortical thickness, even for studies not requiring high spatial resolution, and smoothing them down within the cortical ribbon with a kernel of an appropriate shape to achieve the best performance-thus decoupling the choice of fMRI voxel size from the spatial resolution requirements of the particular study. The improvement of this new intracortical smoothing approach over conventional surface-based smoothing is expected to be modest for conventional resolutions, however the improvement is expected to increase with higher resolutions. This framework can also be applied to anatomically-informed intracortical smoothing of higher-resolution data (e.g. along columns and layers) in studies with prior information about the spatial structure of activation.

Social Responsiveness Scale (SRS) in Relation to Longitudinal Cortical Thickness Changes in Autism Spectrum Disorder

The relationship between brain development and clinical heterogeneity in autism (ASD) is unknown. This study examines the Social Responsiveness Scale (SRS) in relation to the longitudinal development of cortical thickness. Participants (N = 91 ASD, N = 56 TDC; 3-39 years at first scan) were scanned up to three times over a 7-year period. Mixed-effects models examined cortical thickness in relation to SRS score. ASD participants with higher SRS scores showed regionally increased age-related cortical thinning. Regional thickness differences and reduced age-related cortical thinning were found in predominantly right lateralized regions in ASD with decreasing SRS scores over time. Our findings emphasize the importance of examining clinical phenotypes in brain-based studies of ASD.

Discrete Heat Kernel Smoothing in Irregular Image Domains

We present the discrete version of heat kernel smoothing on graph data structure. The method is used to smooth data in an irregularly shaped domains in 3D images. New statistical properties of heat kernel smoothing are derived. As an application, we show how to filter out noisy data in the lung blood vessel trees obtained from computed tomography. The method can be further used in representing the complex vessel trees parametrically as a linear combination of basis functions and extracting the skeleton representation of the trees.

Scale-specific analysis of fMRI data on the irregular cortical surface

To fully characterize the activity patterns on the cerebral cortex as measured with fMRI, the spatial scale of the patterns must be ascertained. Here we address this problem by constructing steerable bandpass filters on the discrete, irregular cortical mesh, using an improved Gaussian smoothing in combination with differential operators of directional derivatives. We demonstrate the utility of the algorithm in two ways. First, using modelling we show that our algorithm yields superior results in numerical precision and spatial uniformity of filter kernels compared to the most widely adopted approach for cortical smoothing. As the effective scales of information differ from the nominal filter sizes applied to extract them, we evaluated the effective scales empirically for different filters to make subsequent comparisons well calibrated. Second, we applied the algorithm to an fMRI dataset to assess the scale and pattern form of cortical encoding of information about visual objects in the ventral visual pathway. We found that filtering by our method improved the detection of discriminant information about experimental conditions over previous methods, that the level of categorization (subordinate versus superordinate) of objects was differentially related to the spatial scale of fMRI patterns, and that the spatial scale at which information was encoded increased along the ventral visual pathway. In sum, our results indicate that the proposed algorithm is particularly suited to assess and detect scale-specific information encoding in cortex, and promises further insight into the topography of cortical encoding in the human brain.