Multiscale adaptive regression models for neuroimaging data: Multiscale Adaptive Regression Models

Neuroimaging studies aim to analyze imaging data with complex spatial patterns in a large number of locations (called voxels) on a two-dimensional (2D) surface or in a 3D volume. Conventional analyses of imaging data include two sequential steps: spatially smoothing imaging data and then independently fitting a statistical model at each voxel. However, conventional analyses suffer from the same amount of smoothing throughout the whole image, the arbitrary choice of smoothing extent, and low statistical power in detecting spatial patterns. We propose a multiscale adaptive regression model (MARM) to integrate the propagation–separation (PS) approach (Polzehl and Spokoiny, 2000, 2006) with statistical modeling at each voxel for spatial and adaptive analysis of neuroimaging data from multiple subjects. MARM has three features: being spatial, being hierarchical, and being adaptive. We use a multiscale adaptive estimation and testing procedure (MAET) to utilize imaging observations from the neighboring voxels of the current voxel to adaptively calculate parameter estimates and test statistics. Theoretically, we establish consistency and asymptotic normality of the adaptive parameter estimates and the asymptotic distribution of the adaptive test statistics. Our simulation studies and real data analysis confirm that MARM significantly outperforms conventional analyses of imaging data

Family Poverty Affects the Rate of Human Infant Brain Growth

Living in poverty places children at very high risk for problems across a variety of domains, including schooling, behavioral regulation, and health. Aspects of cognitive functioning, such as information processing, may underlie these kinds of problems. How might poverty affect the brain functions underlying these cognitive processes? Here, we address this question by observing and analyzing repeated measures of brain development of young children between five months and four years of age from economically diverse backgrounds (n = 77). In doing so, we have the opportunity to observe changes in brain growth as children begin to experience the effects of poverty. These children underwent MRI scanning, with subjects completing between 1 and 7 scans longitudinally. Two hundred and three MRI scans were divided into different tissue types using a novel image processing algorithm specifically designed to analyze brain data from young infants. Total gray, white, and cerebral (summation of total gray and white matter) volumes were examined along with volumes of the frontal, parietal, temporal, and occipital lobes. Infants from low-income families had lower volumes of gray matter, tissue critical for processing of information and execution of actions. These differences were found for both the frontal and parietal lobes. No differences were detected in white matter, temporal lobe volumes, or occipital lobe volumes. In addition, differences in brain growth were found to vary with socioeconomic status (SES), with children from lower-income households having slower trajectories of growth during infancy and early childhood. Volumetric differences were associated with the emergence of disruptive behavioral problems

Probabilistic Air Segmentation and Sparse Regression Estimated Pseudo CT for PET/MR Attenuation Correction

A probabilistic air segmentation and sparse regression method was developed for PET attenuation correction with a mean whole-brain PET error of 2.42% ± 1.0 by estimating continuous pseudo CT images from T1-weighted MR and atlas CT images

Identification of MCI individuals using structural and functional connectivity networks

Different imaging modalities provide essential complementary information that can be used to enhance our understanding of brain disorders. This study focuses on integrating multiple imaging modalities to identify individuals at risk for mild cognitive impairment (MCI). MCI, often an early stage of Alzheimer’s disease (AD), is difficult to diagnose due to its very mild or insignificant symptoms of cognitive impairment. Recent emergence of brain network analysis has made characterization of neurological disorders at a whole-brain connectivity level possible, thus providing new avenues for brain diseases classification. Employing multiple-kernel Support Vector Machines (SVMs), we attempt to integrate information from diffusion tensor imaging (DTI) and resting-state functional magnetic resonance imaging (rs-fMRI) for improving classification performance. Our results indicate that the multimodality classification approach yields statistically significant improvement in accuracy over using each modality independently. The classification accuracy obtained by the proposed method is 96.3%, which is an increase of at least 7.4% from the single modality-based methods and the direct data fusion method. A cross-validation estimation of the generalization performance gives an area of 0.953 under the receiver operating characteristic (ROC) curve, indicating excellent diagnostic power. The multimodality classification approach hence allows more accurate early detection of brain abnormalities with greater sensitivity

Estimating patient-specific and anatomically correct reference model for craniomaxillofacial deformity via sparse representation: Estimating patient-specific and anatomically correct reference model

A significant number of patients suffer from craniomaxillofacial (CMF) deformity and require CMF surgery in the United States. The success of CMF surgery depends on not only the surgical techniques but also an accurate surgical planning. However, surgical planning for CMF surgery is challenging due to the absence of a patient-specific reference model. Currently, the outcome of the surgery is often subjective and highly dependent on surgeon’s experience. In this paper, the authors present an automatic method to estimate an anatomically correct reference shape of jaws for orthognathic surgery, a common type of CMF surgery

Enriched white matter connectivity networks for accurate identification of MCI patients

Mild cognitive impairment (MCI), often a prodromal phase of Alzheimer’s disease (AD), is frequently considered to be good target for early diagnosis and therapeutic interventions of AD. Recent emergence of reliable network characterization techniques has made it possible to understand neurological disorders at a whole-brain connectivity level. Accordingly, we propose an effective network-based multivariate classification algorithm, using a collection of measures derived from white-matter (WM) connectivity networks, to accurately identify MCI patients from normal controls. An enriched description of WM connections, utilizing six physiological parameters, i.e., fiber count, fractional anisotropy (FA), mean diffusivity (MD), and principal diffusivities (λ1, λ2, λ3), results in six connectivity networks for each subject to account for the connection topology and the biophysical properties of the connections. Upon parcellating the brain into 90 regions-of-interest (ROIs), these properties can be quantified for each pair of regions with common traversing fibers. For building an MCI classifier, clustering coefficient of each ROI in relation to the remaining ROIs is extracted as feature for classification. These features are then ranked according to their Pearson correlation with respect to the clinical labels, and are further sieved to select the most discriminant subset of features using a SVM-based feature selection algorithm. Finally, support vector machines (SVMs) are trained using the selected subset of features. Classification accuracy was evaluated via leave-one-out cross-validation to ensure generalization of performance. The classification accuracy given by our enriched description of WM connections is 88.9%, which is an increase of at least 14.8% from that using simple WM connectivity description with any single physiological parameter. A cross-validation estimation of the generalization performance shows an area of 0.929 under the receiver operating characteristic (ROC) curve, indicating excellent diagnostic power. It was also found, based on the selected features, that portions of the prefrontal cortex, orbitofrontal cortex, parietal lobe and insula regions provided the most discriminant features for classification, in line with results reported in previous studies. Our MCI classification framework, especially the enriched description of WM connections, allows accurate early detection of brain abnormalities, which is of paramount importance for treatment management of potential AD patients

Automated bone segmentation from dental CBCT images using patch-based sparse representation and convex optimization: Segmentation of CBCT image

Purpose: Cone-beam computed tomography (CBCT) is an increasingly utilized imaging modality for the diagnosis and treatment planning of the patients with craniomaxillofacial (CMF) deformities. Accurate segmentation of CBCT image is an essential step to generate three-dimensional (3D) models for the diagnosis and treatment planning of the patients with CMF deformities. However, due to the poor image quality, including very low signal-to-noise ratio and the widespread image artifacts such as noise, beam hardening, and inhomogeneity, it is challenging to segment the CBCT images. In this paper, the authors present a new automatic segmentation method to address these problems

Predicting Future Clinical Changes of MCI Patients Using Longitudinal and Multimodal Biomarkers

Accurate prediction of clinical changes of mild cognitive impairment (MCI) patients, including both qualitative change (i.e., conversion to Alzheimer's disease (AD)) and quantitative change (i.e., cognitive scores) at future time points, is important for early diagnosis of AD and for monitoring the disease progression. In this paper, we propose to predict future clinical changes of MCI patients by using both baseline and longitudinal multimodality data. To do this, we first develop a longitudinal feature selection method to jointly select brain regions across multiple time points for each modality. Specifically, for each time point, we train a sparse linear regression model by using the imaging data and the corresponding clinical scores, with an extra ‘group regularization’ to group the weights corresponding to the same brain region across multiple time points together and to allow for selection of brain regions based on the strength of multiple time points jointly. Then, to further reflect the longitudinal changes on the selected brain regions, we extract a set of longitudinal features from the original baseline and longitudinal data. Finally, we combine all features on the selected brain regions, from different modalities, for prediction by using our previously proposed multi-kernel SVM. We validate our method on 88 ADNI MCI subjects, with both MRI and FDG-PET data and the corresponding clinical scores (i.e., MMSE and ADAS-Cog) at 5 different time points. We first predict the clinical scores (MMSE and ADAS-Cog) at 24-month by using the multimodality data at previous time points, and then predict the conversion of MCI to AD by using the multimodality data at time points which are at least 6-month ahead of the conversion. The results on both sets of experiments show that our proposed method can achieve better performance in predicting future clinical changes of MCI patients than the conventional methods

Hierarchical Anatomical Brain Networks for MCI Prediction: Revisiting Volumetric Measures

Owning to its clinical accessibility, T1-weighted MRI (Magnetic Resonance Imaging) has been extensively studied in the past decades for prediction of Alzheimer's disease (AD) and mild cognitive impairment (MCI). The volumes of gray matter (GM), white matter (WM) and cerebrospinal fluid (CSF) are the most commonly used measurements, resulting in many successful applications. It has been widely observed that disease-induced structural changes may not occur at isolated spots, but in several inter-related regions. Therefore, for better characterization of brain pathology, we propose in this paper a means to extract inter-regional correlation based features from local volumetric measurements. Specifically, our approach involves constructing an anatomical brain network for each subject, with each node representing a Region of Interest (ROI) and each edge representing Pearson correlation of tissue volumetric measurements between ROI pairs. As second order volumetric measurements, network features are more descriptive but also more sensitive to noise. To overcome this limitation, a hierarchy of ROIs is used to suppress noise at different scales. Pairwise interactions are considered not only for ROIs with the same scale in the same layer of the hierarchy, but also for ROIs across different scales in different layers. To address the high dimensionality problem resulting from the large number of network features, a supervised dimensionality reduction method is further employed to embed a selected subset of features into a low dimensional feature space, while at the same time preserving discriminative information. We demonstrate with experimental results the efficacy of this embedding strategy in comparison with some other commonly used approaches. In addition, although the proposed method can be easily generalized to incorporate other metrics of regional similarities, the benefits of using Pearson correlation in our application are reinforced by the experimental results. Without requiring new sources of information, our proposed approach improves the accuracy of MCI prediction from (of conventional volumetric features) to (of hierarchical network features), evaluated using data sets randomly drawn from the ADNI (Alzheimer's Disease Neuroimaging Initiative) dataset

Infant Brain Atlases from Neonates to 1- and 2-Year-Olds

Background: Studies for infants are usually hindered by the insufficient image contrast, especially for neonates. Prior knowledge, in the form of atlas, can provide additional guidance for the data processing such as spatial normalization, label propagation, and tissue segmentation. Although it is highly desired, there is currently no such infant atlas which caters for all these applications. The reason may be largely due to the dramatic early brain development, image processing difficulties, and the need of a large sample size. Methodology: To this end, after several years of subject recruitment and data acquisition, we have collected a unique longitudinal dataset, involving 95 normal infants (56 males and 39 females) with MRI scanned at 3 ages, i.e., neonate, 1-yearold, and 2-year-old. State-of-the-art MR image segmentation and registration techniques were employed, to construct which include the templates (grayscale average images), tissue probability maps (TPMs), and brain parcellation maps (i.e., meaningful anatomical regions of interest) for each age group. In addition, the longitudinal correspondences between agespecific atlases were also obtained. Experiments of typical infant applications validated that the proposed atlas outperformed other atlases and is hence very useful for infant-related studies. Conclusions: We expect that the proposed infant 0–1–2 brain atlases would be significantly conducive to structural and functional studies of the infant brains. These atlases are publicly available in our website