Multivariate tensor-based morphometry on surfaces: application to mapping ventricular abnormalities in HIV/AIDS

Surface-based morphometry of the corpus callosum in young children of ages 1-5

The corpus callosum (CC) is a large white matter fiber bundle in the brain and is involved in various cognitive, sensory, and motor processes. While implicated in various developmental and psychiatric disorders, much is yet to be uncovered about the normal development of this structure, especially in young children. Additionally, while sexual dimorphism has been reported in prior literature, observations have not necessarily been consistent. In this study, we use morphometric measures including surface tensor-based morphometry (TBM) to investigate local changes in the shape of the CC in children between the ages of 12 and 60 months, in intervals of 12 months. We also analyze sex differences in each of these age groups. We observed larger significant clusters in the earlier ages between 12 v 24 m and between 48 v 60 m and localized differences in the anterior region of the body of the CC. Sex differences were most pronounced in the 12 m group. This study adds to the growing literature of work aiming to understand the developing brain and emphasizes the utility of surface TBM as a useful tool for analyzing regional differences in neuroanatomical morphometry.

Topological Embedding of Human Brain Networks with Applications to Dynamics of Temporal Lobe Epilepsy

We introduce a novel, data-driven topological data analysis (TDA) approach for embedding brain networks into a lower-dimensional space in quantifying the dynamics of temporal lobe epilepsy (TLE) obtained from resting-state functional magnetic resonance imaging (rs-fMRI). This embedding facilitates the orthogonal projection of 0D and 1D topological features, allowing for the visualization and modeling of the dynamics of functional human brain networks in a resting state. We then quantify the topological disparities between networks to determine the coordinates for embedding. This framework enables us to conduct a coherent statistical inference within the embedded space. Our results indicate that brain network topology in TLE patients exhibits increased rigidity in 0D topology but more rapid flections compared to that of normal controls in 1D topology.

Combining Blood-Based Biomarkers and Structural MRI Measurements to Distinguish Persons with and without Significant Amyloid Plaques

Background

Amyloid-β (Aβ) plaques play a pivotal role in Alzheimer's disease. The current positron emission tomography (PET) is expensive and limited in availability. In contrast, blood-based biomarkers (BBBMs) show potential for characterizing Aβ plaques more affordably. We have previously proposed an MRI-based hippocampal morphometry measure to be an indicator of Aβ plaques.

Objective

To develop and validate an integrated model to predict brain amyloid PET positivity combining MRI feature and plasma Aβ42/40 ratio.

Methods

We extracted hippocampal multivariate morphometry statistics from MR images and together with plasma Aβ42/40 trained a random forest classifier to perform a binary classification of participant brain amyloid PET positivity. We evaluated the model performance using two distinct cohorts, one from the Alzheimer's Disease Neuroimaging Initiative (ADNI) and the other from the Banner Alzheimer's Institute (BAI), including prediction accuracy, precision, recall rate, F1 score, and AUC score.

Results

Results from ADNI (mean age 72.6, Aβ+ rate 49.5%) and BAI (mean age 66.2, Aβ+ rate 36.9%) datasets revealed the integrated multimodal (IMM) model's superior performance over unimodal models. The IMM model achieved prediction accuracies of 0.86 in ADNI and 0.92 in BAI, surpassing unimodal models based solely on structural MRI (0.81 and 0.87) or plasma Aβ42/40 (0.73 and 0.81) predictors.

CONCLUSIONS

Our IMM model, combining MRI and BBBM data, offers a highly accurate approach to predict brain amyloid PET positivity. This innovative multiplex biomarker strategy presents an accessible and cost-effective avenue for advancing Alzheimer's disease diagnostics, leveraging diverse pathologic features related to Aβ plaques and structural MRI.

Correlation studies of Hippocampal Morphometry and Plasma NFL Levels in Cognitively Unimpaired Subjects

Alzheimer's disease(AD) is being the burden of society and family. Applying computing-aided strategies to reveal its pathology is one of the research highlights. Plasma neurofilament light (NFL) is an emerging noninvasive and economic biomarker for AD molecular pathology. It is valuable to reveal the correlations between the plasma NFL levels and neurodegeneration, especially hippcampal deformations at the preclinical stage. The negative correlation between plasma NFL levels and hippocampal volumes has been documented. However, the relationship between the plasma NFL levels and the hippocampal morphometry details at the preclinical stage is still elusive. This study seeks to demonstrate the capacity of our proposed surface-based hippocampal morphometry system to discern the plasma NFL positive (NFL+>41.9 pg/L) level and plasma NFL negative (NFL-<41.9pg/L) level and illustrate its superiority to the hippocampal volume measurement by drawing the cohort of 154 CU middle aged and elderly adults. We also apply this morphometry measure and a proposed sparse coding based classification algorithm to classify CU individuals with NFL+ and NFL- levels. Experimental results show that the proposed hippocampal morphometry system offers stronger statistical power to discriminate CU subjects with NFL+ and NFL- levels, comparing with the hippocampal volume measure. Furthermore, this system can discriminate plasma NFL levels in CU individuals (Accuracy=0.86). Both the group level and individual level analysis results indicate that the association between plasma NFL levels and the hippocampal shapes can be mapped at the preclinical stage.

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.

Every individual makes a difference: A trinity derived from linking individual brain morphometry, connectivity and mentalising ability

Mentalising ability, indexed as the ability to understand others' beliefs, feelings, intentions, thoughts and traits, is a pivotal and fundamental component of human social cognition. However, considering the multifaceted nature of mentalising ability, little research has focused on characterising individual differences in different mentalising components. And even less research has been devoted to investigating how the variance in the structural and functional patterns of the amygdala and hippocampus, two vital subcortical regions of the "social brain", are related to inter-individual variability in mentalising ability. Here, as a first step toward filling these gaps, we exploited inter-subject representational similarity analysis (IS-RSA) to assess relationships between amygdala and hippocampal morphometry (surface-based multivariate morphometry statistics, MMS), connectivity (resting-state functional connectivity, rs-FC) and mentalising ability (interactive mentalisation questionnaire [IMQ] scores) across the participants ( N = 24 $$ N=24 $$ ). In IS-RSA, we proposed a novel pipeline, that is, computing patching and pooling operations-based surface distance (CPP-SD), to obtain a decent representation for high-dimensional MMS data. On this basis, we found significant correlations (i.e., second-order isomorphisms) between these three distinct modalities, indicating that a trinity existed in idiosyncratic patterns of brain morphometry, connectivity and mentalising ability. Notably, a region-related mentalising specificity emerged from these associations: self-self and self-other mentalisation are more related to the hippocampus, while other-self mentalisation shows a closer link with the amygdala. Furthermore, by utilising the dyadic regression analysis, we observed significant interactions such that subject pairs with similar morphometry had even greater mentalising similarity if they were also similar in rs-FC. Altogether, we demonstrated the feasibility and illustrated the promise of using IS-RSA to study individual differences, deepening our understanding of how individual brains give rise to their mentalising abilities.

Improved Prediction of Amyloid-β and Tau Burden Using Hippocampal Surface Multivariate Morphometry Statistics and Sparse Coding

BACKGROUND

Amyloid-β (Aβ) plaques and tau protein tangles in the brain are the defining 'A' and 'T' hallmarks of Alzheimer's disease (AD), and together with structural atrophy detectable on brain magnetic resonance imaging (MRI) scans as one of the neurodegenerative ('N') biomarkers comprise the "ATN framework" of AD. Current methods to detect Aβ/tau pathology include cerebrospinal fluid (invasive), positron emission tomography (PET; costly and not widely available), and blood-based biomarkers (promising but mainly still in development).

OBJECTIVE

To develop a non-invasive and widely available structural MRI-based framework to quantitatively predict the amyloid and tau measurements.

METHODS

With MRI-based hippocampal multivariate morphometry statistics (MMS) features, we apply our Patch Analysis-based Surface Correntropy-induced Sparse coding and max-pooling (PASCS-MP) method combined with the ridge regression model to individual amyloid/tau measure prediction.

RESULTS

We evaluate our framework on amyloid PET/MRI and tau PET/MRI datasets from the Alzheimer's Disease Neuroimaging Initiative. Each subject has one pair consisting of a PET image and MRI scan, collected at about the same time. Experimental results suggest that amyloid/tau measurements predicted with our PASCP-MP representations are closer to the real values than the measures derived from other approaches, such as hippocampal surface area, volume, and shape morphometry features based on spherical harmonics.

CONCLUSION

The MMS-based PASCP-MP is an efficient tool that can bridge hippocampal atrophy with amyloid and tau pathology and thus help assess disease burden, progression, and treatment effects.

A Surface-Based Federated Chow Test Model for Integrating APOE Status, Tau Deposition Measure, and Hippocampal Surface Morphometry

BACKGROUND

Alzheimer's disease (AD) is the most common type of age-related dementia, affecting 6.2 million people aged 65 or older according to CDC data. It is commonly agreed that discovering an effective AD diagnosis biomarker could have enormous public health benefits, potentially preventing or delaying up to 40% of dementia cases. Tau neurofibrillary tangles are the primary driver of downstream neurodegeneration and subsequent cognitive impairment in AD, resulting in structural deformations such as hippocampal atrophy that can be observed in magnetic resonance imaging (MRI) scans.

OBJECTIVE

To build a surface-based model to 1) detect differences between APOE subgroups in patterns of tau deposition and hippocampal atrophy, and 2) use the extracted surface-based features to predict cognitive decline.

METHODS

Using data obtained from different institutions, we develop a surface-based federated Chow test model to study the synergistic effects of APOE, a previously reported significant risk factor of AD, and tau on hippocampal surface morphometry.

RESULTS

We illustrate that the APOE-specific morphometry features correlate with AD progression and better predict future AD conversion than other MRI biomarkers. For example, a strong association between atrophy and abnormal tau was identified in hippocampal subregion cornu ammonis 1 (CA1 subfield) and subiculum in e4 homozygote cohort.

CONCLUSION

Our model allows for identifying MRI biomarkers for AD and cognitive decline prediction and may uncover a corner of the neural mechanism of the influence of APOE and tau deposition on hippocampal morphology.

LESA: Longitudinal Elastic Shape Analysis of Brain Subcortical Structures

Over the past 30 years, magnetic resonance imaging has become a ubiquitous tool for accurately visualizing the change and development of the brain's subcortical structures (e.g., hippocampus). Although subcortical structures act as information hubs of the nervous system, their quantification is still in its infancy due to many challenges in shape extraction, representation, and modeling. Here, we develop a simple and efficient framework of longitudinal elastic shape analysis (LESA) for subcortical structures. Integrating ideas from elastic shape analysis of static surfaces and statistical modeling of sparse longitudinal data, LESA provides a set of tools for systematically quantifying changes of longitudinal subcortical surface shapes from raw structure MRI data. The key novelties of LESA include: (i) it can efficiently represent complex subcortical structures using a small number of basis functions and (ii) it can accurately delineate the spatiotemporal shape changes of the human subcortical structures. We applied LESA to analyze three longitudinal neuroimaging data sets and showcase its wide applications in estimating continuous shape trajectories, building life-span growth patterns, and comparing shape differences among different groups. In particular, with the Alzheimer's Disease Neuroimaging Initiative (ADNI) data, we found that the Alzheimer's Disease (AD) can significantly speed the shape change of ventricle and hippocampus from 60 to 75 years old compared with normal aging.

The impact of Alzheimer's disease susceptibility loci on lateral ventricular surface morphology in older adults

The enlargement of ventricular volume is a general trend in the elderly, especially in patients with Alzheimer's disease (AD). Multiple susceptibility loci have been reported to have an increased risk for AD and the morphology of brain structures are affected by the variations in the risk loci. Therefore, we hypothesized that genes contributed significantly to the ventricular surface, and the changes of ventricular surface were associated with the impairment of cognitive functions. After the quality controls (QC) and genotyping, a lateral ventricular segmentation method was employed to obtain the surface features of lateral ventricle. We evaluated the influence of 18 selected AD susceptibility loci on both volume and surface morphology across 410 subjects from Alzheimer's Disease Neuroimaging Initiative (ADNI). Correlations were conducted between radial distance (RD) and Montreal Cognitive Assessment (MoCA) subscales. Only the C allele at the rs744373 loci in BIN1 gene significantly accelerated the atrophy of lateral ventricle, including the anterior horn, body, and temporal horn of left lateral ventricle. No significant effect on lateral ventricle was found at other loci. Our results revealed that most regions of the bilateral ventricular surface were significantly negatively correlated with cognitive scores, particularly in delayed recall. Besides, small areas of surface were negatively correlated with language, orientation, and visuospatial scores. Together, our results indicated that the genetic variation affected the localized areas of lateral ventricular surface, and supported that lateral ventricle was an important brain structure associated with cognition in the elderly.

Altered Glymphatic System in Middle-Aged cART-Treated Patients With HIV: A Diffusion Tensor Imaging Study

Objective:

The brain relies on the glymphatic system to clear metabolic wastes and maintain brain homeostasis to fulfill its functions better. Yet, the complexity of the glymphatic flow and clearance and its changes in HIV infection and its role in neurocognitive dysfunction remain poorly understood. This study aims to explore the impact of HIV and combination antiretroviral therapy (cART) on the glymphatic system and establish a potential biomarker of HIV-associated neurocognitive disorders (HAND).

Methods

Here, we examined the glymphatic profiles of middle-aged virosuppressed patients with HIV (n = 27) receiving cART over 1–6 years and healthy controls (n = 28) along the perivascular space (PVS) using diffusion tensor image analysis along the perivascular space (ALPS) with guided and unguided approaches. We later combined data from these analyses to investigate MRI glymphatic correlates of cognitive impairment and other clinical tests of HIV (CD4+ T-cell counts and CD4+/CD8+ ratio).

Results

We found that glymphatic function as measured by the ALPS index increased significantly in the right and left PVSs of patients with HIV having cART. On antiretroviral therapy, a changing pattern in glymphatic clearance function in patients with HIV having cART correlated with attention and working memory. Duration on cART was also associated with cognitive performances of abstract and executive function and learning and memory.

Conclusion

These findings provide MRI evidence of the presence of HIV-induced changes in the glymphatic flow and clearance, which might underlie cognitive impairment among patients with HIV having cART. An increase in the glymphatic activity might reflect a compensatory mechanism to regulate microenvironment homeostasis compromised by HIV. This compensation might be necessary to maintain the proper functioning of the brain while coping with HIV pathology. These findings also shed light on the clinical importance of evaluating glymphatic function based on the ALPS index and suggest that improving the glymphatic system may serve as an alternative therapeutic strategy for HAND.

Integrating Transcriptomics, Genomics, and Imaging in Alzheimer's Disease: A Federated 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. In clinical practice, structural magnetic resonance imaging (sMRI) is the most accessible and widely used diagnostic imaging modality. Additionally, genome-wide association studies (GWAS) and transcriptomics-the study of gene expression-also play an important role in understanding AD etiology and progression. Sophisticated imaging genetics systems have been developed to discover genetic factors that consistently affect brain function and structure. However, most studies to date focused on the relationships between brain sMRI and GWAS or brain sMRI and transcriptomics. To our knowledge, few methods have been developed to discover and infer multimodal relationships among sMRI, GWAS, and transcriptomics. To address this, we propose a novel federated model, Genotype-Expression-Imaging Data Integration (GEIDI), to identify genetic and transcriptomic influences on brain sMRI measures. The relationships between brain imaging measures and gene expression are allowed to depend on a person's genotype at the single-nucleotide polymorphism (SNP) level, making the inferences adaptive and personalized. We performed extensive experiments on publicly available Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset. Experimental results demonstrated our proposed method outperformed state-of-the-art expression quantitative trait loci (eQTL) methods for detecting genetic and transcriptomic factors related to AD and has stable performance when data are integrated from multiple sites. Our GEIDI approach may offer novel insights into the relationship among image biomarkers, genotypes, and gene expression and help discover novel genetic targets for potential AD drug treatments.

Quantitative characterization of the human retinotopic map based on quasiconformal mapping

The retinotopic map depicts the cortical neurons' response to visual stimuli on the retina and has contributed significantly to our understanding of human visual system. Although recent advances in high field functional magnetic resonance imaging (fMRI) have made it possible to generate the in vivo retinotopic map with great detail, quantifying the map remains challenging. Existing quantification methods do not preserve surface topology and often introduce large geometric distortions to the map. In this study, we developed a new framework based on computational conformal geometry and quasiconformal Teichmüller theory to quantify the retinotopic map. Specifically, we introduced a general pipeline, consisting of cortical surface conformal parameterization, surface-spline-based cortical activation signal smoothing, and vertex-wise Beltrami coefficient-based map description. After correcting most of the violations of the topological conditions, the result was a "Beltrami coefficient map" (BCM) that rigorously and completely characterizes the retinotopic map by quantifying the local quasiconformal mapping distortion at each visual field location. The BCM provided topological and fully reconstructable retinotopic maps. We successfully applied the new framework to analyze the V1 retinotopic maps from the Human Connectome Project (n=181), the largest state of the art retinotopy dataset currently available. With unprecedented precision, we found that the V1 retinotopic map was quasiconformal and the local mapping distortions were similar across observers. The new framework can be applied to other visual areas and retinotopic maps of individuals with and without eye diseases, and improve our understanding of visual cortical organization in normal and clinical populations.

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.

Predicting Tau Accumulation in Cerebral Cortex with Multivariate MRI Morphometry Measurements, Sparse Coding, and Correntropy

Biomarker-assisted diagnosis and intervention in Alzheimer's disease (AD) may be the key to prevention breakthroughs. One of the hallmarks of AD is the accumulation of tau plaques in the human brain. However, current methods to detect tau pathology are either invasive (lumbar puncture) or quite costly and not widely available (Tau PET). In our previous work, structural MRI-based hippocampal multivariate morphometry statistics (MMS) showed superior performance as an effective neurodegenerative biomarker for preclinical AD and Patch Analysis-based Surface Correntropy-induced Sparse coding and max-pooling (PASCS-MP) has excellent ability to generate low-dimensional representations with strong statistical power for brain amyloid prediction. In this work, we apply this framework together with ridge regression models to predict Tau deposition in Braak12 and Braak34 brain regions separately. We evaluate our framework on 925 subjects from the Alzheimer's Disease Neuroimaging Initiative (ADNI). Each subject has one pair consisting of a PET image and MRI scan which were collected at about the same times. Experimental results suggest that the representations from our MMS and PASCS-MP have stronger predictive power and their predicted Braak12 and Braak34 are closer to the real values compared to the measures derived from other approaches such as hippocampal surface area and volume, and shape morphometry features based on spherical harmonics (SPHARM).

Cortical Surface Shape Analysis Based on Alexandrov Polyhedra

Shape analysis has been playing an important role in early diagnosis and prognosis of neurodegenerative diseases such as Alzheimer's diseases (AD). However, obtaining effective shape representations remains challenging. This paper proposes to use the Alexandrov polyhedra as surface-based shape signatures for cortical morphometry analysis. Given a closed genus-0 surface, its Alexandrov polyhedron is a convex representation that encodes its intrinsic geometry information. We propose to compute the polyhedra via a novel spherical optimal transport (OT) computation. In our experiments, we observe that the Alexandrov polyhedra of cortical surfaces between pathology-confirmed AD and cognitively unimpaired individuals are significantly different. Moreover, we propose a visualization method by comparing local geometry differences across cortical surfaces. We show that the proposed method is effective in pinpointing regional cortical structural changes impacted by AD.

Predicting Brain Amyloid Using Multivariate Morphometry Statistics, Sparse Coding, and Correntropy: Validation in 1,101 Individuals From the ADNI and OASIS Databases

Biomarker assisted preclinical/early detection and intervention in Alzheimer’s disease (AD) may be the key to therapeutic breakthroughs. One of the presymptomatic hallmarks of AD is the accumulation of beta-amyloid (Aβ) plaques in the human brain. However, current methods to detect Aβ pathology are either invasive (lumbar puncture) or quite costly and not widely available (amyloid PET). Our prior studies show that magnetic resonance imaging (MRI)-based hippocampal multivariate morphometry statistics (MMS) are an effective neurodegenerative biomarker for preclinical AD. Here we attempt to use MRI-MMS to make inferences regarding brain Aβ burden at the individual subject level. As MMS data has a larger dimension than the sample size, we propose a sparse coding algorithm, Patch Analysis-based Surface Correntropy-induced Sparse-coding and Max-Pooling (PASCS-MP), to generate a low-dimensional representation of hippocampal morphometry for each individual subject. Then we apply these individual representations and a binary random forest classifier to predict brain Aβ positivity for each person. We test our method in two independent cohorts, 841 subjects from the Alzheimer’s Disease Neuroimaging Initiative (ADNI) and 260 subjects from the Open Access Series of Imaging Studies (OASIS). Experimental results suggest that our proposed PASCS-MP method and MMS can discriminate Aβ positivity in people with mild cognitive impairment (MCI) [Accuracy (ACC) = 0.89 (ADNI)] and in cognitively unimpaired (CU) individuals [ACC = 0.79 (ADNI) and ACC = 0.81 (OASIS)]. These results compare favorably relative to measures derived from traditional algorithms, including hippocampal volume and surface area, shape measures based on spherical harmonics (SPHARM) and our prior Patch Analysis-based Surface Sparse-coding and Max-Pooling (PASS-MP) methods.

Neuroimaging PheWAS (Phenome-Wide Association Study): A Free Cloud-Computing Platform for Big-Data, Brain-Wide Imaging Association Studies

Large-scale, case-control genome-wide association studies (GWASs) have revealed genetic variations associated with diverse neurological and psychiatric disorders. Recent advances in neuroimaging and genomic databases of large healthy and diseased cohorts have empowered studies to characterize effects of the discovered genetic factors on brain structure and function, implicating neural pathways and genetic mechanisms in the underlying biology. However, the unprecedented scale and complexity of the imaging and genomic data requires new advanced biomedical data science tools to manage, process and analyze the data. In this work, we introduce Neuroimaging PheWAS (phenome-wide association study): a web-based system for searching over a wide variety of brain-wide imaging phenotypes to discover true system-level gene-brain relationships using a unified genotype-to-phenotype strategy. This design features a user-friendly graphical user interface (GUI) for anonymous data uploading, study definition and management, and interactive result visualizations as well as a cloud-based computational infrastructure and multiple state-of-art methods for statistical association analysis and multiple comparison correction. We demonstrated the potential of Neuroimaging PheWAS with a case study analyzing the influences of the apolipoprotein E (APOE) gene on various brain morphological properties across the brain in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. Benchmark tests were performed to evaluate the system's performance using data from UK Biobank. The Neuroimaging PheWAS system is freely available. It simplifies the execution of PheWAS on neuroimaging data and provides an opportunity for imaging genetics studies to elucidate routes at play for specific genetic variants on diseases in the context of detailed imaging phenotypic data.

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.

A novel pipeline leveraging surface-based features of small subcortical structures to classify individuals with autism spectrum disorder

Autism spectrum disorder (ASD) is accompanied with widespread impairment in social-emotional functioning. Classification of ASD using sensitive morphological features derived from structural magnetic resonance imaging (MRI) of the brain may help us to better understand ASD-related mechanisms and improve related automatic diagnosis. Previous studies using T1 MRI scans in large heterogeneous ABIDE dataset with typical development (TD) controls reported poor classification accuracies (around 60%). This may because they only considered surface-based morphometry (SBM) as scalar estimates (such as cortical thickness and surface area) and ignored the neighboring intrinsic geometry information among features. In recent years, the shape-related SBM achieves great success in discovering the disease burden and progression of other brain diseases. However, when focusing on local geometry information, its high dimensionality requires careful treatment in its application to machine learning. To address the above challenges, we propose a novel pipeline for ASD classification, which mainly includes the generation of surface-based features, patch-based surface sparse coding and dictionary learning, Max-pooling and ensemble classifiers based on adaptive optimizers. The proposed pipeline may leverage the sensitivity of brain surface morphometry statistics and the efficiency of sparse coding and Max-pooling. By introducing only the surface features of bilateral hippocampus that derived from 364 male subjects with ASD and 381 age-matched TD males, this pipeline outperformed five recent MRI-based ASD classification studies with >80% accuracy in discriminating individuals with ASD from TD controls. Our results suggest shape-related SBM features may further boost the classification performance of MRI between ASD and TD.

Neuroimaging Advances in Diagnosis and Differentiation of HIV, Comorbidities, and Aging in the cART Era

In the "cART era" of more widely available and accessible treatment, aging and HIV-related comorbidities, including symptoms of brain dysfunction, remain common among HIV-infected individuals on suppressive treatment. A better understanding of the neurobiological consequences of HIV infection is essential for developing thorough treatment guidelines and for optimizing long-term neuropsychological outcomes and overall brain health. In this chapter, we first summarize magnetic resonance imaging (MRI) methods used in over two decades of neuroHIV research. These methods evaluate brain volumetric differences and circuitry disruptions in adults living with HIV, and help map clinical correlations with brain function and tissue microstructure. We then introduce and discuss aging and associated neurological complications in people living with HIV, and processes by which infection may contribute to the risk for late-onset dementias. We describe how new technologies and large-scale international collaborations are helping to disentangle the effect of genetic and environmental risk factors on brain aging and neurodegenerative diseases. We provide insights into how these advances, which are now at the forefront of Alzheimer's disease research, may advance the field of neuroHIV. We conclude with a summary of how we see the field of neuroHIV research advancing in the decades to come and highlight potential clinical implications.

Computing Univariate Neurodegenerative Biomarkers with Volumetric Optimal Transportation: A Pilot Study

Changes in cognitive performance due to neurodegenerative diseases such as Alzheimer's disease (AD) are closely correlated to the brain structure alteration. A univariate and personalized neurodegenerative biomarker with strong statistical power based on magnetic resonance imaging (MRI) will benefit clinical diagnosis and prognosis of neurodegenerative diseases. However, few biomarkers of this type have been developed, especially those that are robust to image noise and applicable to clinical analyses. In this paper, we introduce a variational framework to compute optimal transportation (OT) on brain structural MRI volumes and develop a univariate neuroimaging index based on OT to quantify neurodegenerative alterations. Specifically, we compute the OT from each image to a template and measure the Wasserstein distance between them. The obtained Wasserstein distance, Wasserstein Index (WI) for short to specify the distance to a template, is concise, informative and robust to random noise. Comparing to the popular linear programming-based OT computation method, our framework makes use of Newton's method, which makes it possible to compute WI in large-scale datasets. Experimental results, on 314 subjects (140 Aβ + AD and 174 Aβ- normal controls) from the Alzheimer's Disease Neuroimaging Initiative (ADNI) baseline dataset, provide preliminary evidence that the proposed WI is correlated with a clinical cognitive measure (the Mini-Mental State Examination (MMSE) score), and it is able to identify group difference and achieve a good classification accuracy, outperforming two other popular univariate indices including hippocampal volume and entorhinal cortex thickness. The current pilot work suggests the application of WI as a potential univariate neurodegenerative biomarker.

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.

Brain structural correlates of trajectories to cognitive impairment in men with and without HIV disease

There are distinct trajectories to cognitive impairment among participants in the Multicenter AIDS Cohort Study (MACS). Here we analyzed the relationship between regional brain volumes and the individual trajectories to impairment in a subsample (n = 302) of the cohort. 302 (167 HIV-infected; mean age = 55.7 yrs.; mean education: 16.2 yrs.) of the men enrolled in the MACS MRI study contributed data to this analysis. We used voxel-based morphometry (VBM) to segment the brain images to analyze gray and white matter volume at the voxel-level. A Mixed Membership Trajectory Model had previously identified three distinct profiles, and each study participant had a membership weight for each of these three trajectories. We estimated VBM model parameters for 100 imputations, manually performed the post-hoc contrasts, and pooled the results. We examined the associations between brain volume at the voxel level and the MMTM membership weights for two profiles: one considered "unhealthy" and the other considered "Premature aging." The unhealthy profile was linked to the volume of the posterior cingulate gyrus/precuneus, the inferior frontal cortex, and the insula, whereas the premature aging profile was independently associated with the integrity of a portion of the precuneus. Trajectories to cognitive impairment are the result, in part, of atrophy in cortical regions linked to normal and pathological aging. These data suggest the possibility of predicting cognitive morbidity based on patterns of CNS atrophy.

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.

Subcortical shape and neuropsychological function among U.S. service members with mild traumatic brain injury

In a recent manuscript, our group demonstrated shape differences in the thalamus, nucleus accumbens, and amygdala in a cohort of U.S. Service Members with mild traumatic brain injury (mTBI). Given the significant role these structures play in cognitive function, this study directly examined the relationship between shape metrics and neuropsychological performance. The imaging and neuropsychological data from 135 post-deployed United States Service Members from two groups (mTBI and orthopedic injured) were examined. Two shape features modeling local deformations in thickness (RD) and surface area (JD) were defined vertex-wise on parametric mesh-representations of 7 bilateral subcortical gray matter structures. Linear regression was used to model associations between subcortical morphometry and neuropsychological performance as a function of either TBI status or, among TBI patients, subjective reporting of initial concussion severity (CS). Results demonstrated several significant group-by-cognition relationships with shape metrics across multiple cognitive domains including processing speed, memory, and executive function. Higher processing speed was robustly associated with more dilation of caudate surface area among patients with mTBI who reported more than one CS variables (loss of consciousness (LOC), alteration of consciousness (AOC), and/or post-traumatic amnesia (PTA)). These significant patterns indicate the importance of subcortical structures in cognitive performance and support a growing functional neuroanatomical literature in TBI and other neurologic disorders. However, prospective research will be required before exact directional evolution and progression of shape can be understood and utilized in predicting or tracking cognitive outcomes in this patient population.

Applying surface-based hippocampal morphometry to study APOE-E4 allele dose effects in cognitively unimpaired subjects

Apolipoprotein E (APOE) e4 is the major genetic risk factor for late-onset Alzheimer's disease (AD). The dose-dependent impact of this allele on hippocampal volumes has been documented, but its influence on general hippocampal morphology in cognitively unimpaired individuals is still elusive. Capitalizing on the study of a large number of cognitively unimpaired late middle aged and older adults with two, one and no APOE-e4 alleles, the current study aims to characterize the ability of our automated surface-based hippocampal morphometry algorithm to distinguish between these three levels of genetic risk for AD and demonstrate its superiority to a commonly used hippocampal volume measurement. We examined the APOE-e4 dose effect on cross-sectional hippocampal morphology analysis in a magnetic resonance imaging (MRI) database of 117 cognitively unimpaired subjects aged between 50 and 85 years (mean = 57.4, SD = 6.3), including 36 heterozygotes (e3/e4), 37 homozygotes (e4/e4) and 44 non-carriers (e3/e3). The proposed automated framework includes hippocampal surface segmentation and reconstruction, higher-order hippocampal surface correspondence computation, and hippocampal surface deformation analysis with multivariate statistics. In our experiments, the surface-based method identified APOE-e4 dose effects on the left hippocampal morphology. Compared to the widely-used hippocampal volume measure, our hippocampal morphometry statistics showed greater statistical power by distinguishing cognitively unimpaired subjects with two, one, and no APOE-e4 alleles. Our findings mirrored previous studies showing that APOE-e4 has a dose effect on the acceleration of brain structure deformities. The results indicated that the proposed surface-based hippocampal morphometry measure is a potential preclinical AD imaging biomarker for cognitively unimpaired individuals.

•

Applied surface-based hippocampal morphometry on cognitively unimpaired subjects.

•

Our study identified APOE-e4 dose effects on cognitively unimpaired subjects.

•

Surface-based hippocampal morphometry outperformed the hippocampal volume measure.

•

Surface-based hippocampal morphometry may be a potential preclinical AD biomarker.

Combined morphologic and metabolic pipeline for Positron emission tomography/computed tomography based radiotherapy response evaluation in locally advanced pancreatic adenocarcinoma

•

A novel morphologic and metabolic combined pipeline for PA response evaluation.

•

The derived metric outperformed traditional imaging metrics in risk stratification.

•

May serve as a new image biomarker to characterize heterogeneous tumor response.

Background and purpose

Adaptive radiation planning for pancreatic adenocarcinoma (PA) relies on accurate treatment response assessment, while traditional response evaluation criteria inefficiently characterize tumors with complex morphological features or intrinsically low metabolism. To better assess treatment response of PA, we quantify and compare regional morphological and metabolic features of the 3D pre- and post-radiation therapy (RT) tumor models.

Materials and methods

Thirty-one PA patients with pre and post-RT Positron emission tomography/computed tomography (PET/CT) scans were evaluated. 3D meshes of pre- and post-RT tumors were generated and registered to establish vertex-wise correspondence. To assess tumor response, Mahalanobis distances (M_dist|_Fusion) between pre- and post-RT tumor surfaces with anatomic and metabolic fused vectors were calculated for each patient. M_dist|_Fusion was evaluated by overall survival (OS) prediction and survival risk classification. As a comparison, the same analyses were conducted on traditional imaging/physiological predictors, and distances measurements based on metabolic and morphological features only.

Results

Among all the imaging/physiological parameters, M_dist|_Fusion was shown to be the best predictor of OS (HR = 0.52, p = 0.008), while other parameters failed to reach significance. Moreover, M_dist|_Fusion outperformed traditional morphologic and metabolic measurements in patient risk stratification, either alone (HR = 11.51, p < 0.001) or combined with age (HR = 9.04, p < 0.001).

Conclusions

We introduced a PET/CT-based novel morphologic and metabolic pipeline for response evaluation in locally advanced PA. The fused M_dist|_Fusion outperformed traditional morphologic, metabolic, and physiological measurements in OS prediction and risk stratification. The novel fusion model may serve as a new imaging-marker to more accurately characterize the heterogeneous tumor RT response.

Discriminating lung adenocarcinoma from lung squamous cell carcinoma using respiration-induced tumor shape changes

Based on 4D-CT, we aimed to characterize the pattern of morphological changes in lung tumors during respiration, and investigated its potential in non-invasively differentiating lung adenocarcinoma (AC) and squamous cell carcinoma (SCC). We applied a 3D surface analysis on 22 tumors (13 AC, 9 SCC) to investigate the tumor regional morphological fluctuations in response to respiration phases. Tumor surface vertices among ten respiratory phases were matched using surface-based registration, and the shape descriptors (ρ and detJ) were calculated and tracked across respiration stages in a regionally aligned scenario. Pair-wise group comparisons were performed between lung AC and SCC subtypes, in terms of ratios of maximal shape changes as well as correlation coefficients between tumor shape and respiratory stage indicators from the lung. AC type tumors had averaged larger surface measurements at exhale than at inhale, and these surface measurements were negatively correlated with lung volumes across respiratory stages. In contrast, SCC type tumors had averaged smaller surface measurements at exhale than at inhale, and the correlations with lung volumes were positive. The group differences in maximal shape changes as well as correlations were both statistically significant (p  <  0.05). We developed a non-invasive lung tumor sub-type detection pipeline based on respiration-induced tumor surface deformation. Significant differences in deformation patterns were detected between lung AC and SCC. The derived surface measurements may potentially serve as a new non-invasive imaging biomarker of lung cancer subtypes.

HIPPOCAMPUS MORPHOMETRY STUDY ON PATHOLOGY-CONFIRMED ALZHEIMER'S DISEASE PATIENTS WITH SURFACE MULTIVARIATE MORPHOMETRY STATISTICS

Alzheimer's disease (AD) is one of the most prevalent neurodegenerative diseases in elderly and the incidence of this disease is increasing with older ages. One of the hallmarks of AD is the accumulation of beta-amyloid plaques (Aβ) in human brains. Most of prior brain imaging researchers used the clinical symptom based diagnosis without the confirmation of imaging or fluid Aβ information. In this work, we study hippocampus morphometry on a cohort consisting of Aβ positive AD (N = 151) and matched Aβ negative cognitively unimpaired subjects (N = 271) with Aβ positivity determined via florbetapir PET. The brain images are obtained from publicly available Alzheimer's Disease Neuroimaging Initiative (ADNI). We compute our surface multivariate morphometry statistics from segmented hippocampus structure in structural MR images. With these features, we find statistically significant difference by using Hotelling's T² tests. Meanwhile, we apply a patch-based analysis of sparse coding system for binary group classification and achieve an accuracy rate of 90.48%. Our results demonstrate that our surface multivariate morphometry statistics (MMS) perform better than traditional hippocampal volume measures in classification and it may be applied as a potential biomarker for distinguishing dementia due to AD from age matched normal aging individuals.

Contribution to speech development of the right anterior putamen revealed with multivariate tensor-based morphometry

In our previous study1, we suggested that the difference between tensor-based metrics in the anterior part of the right putamen between 21 and 18 months age groups associated with speech development during this ages. Here we used a correlational analysis between verbal scores and determinant of the Jacobian matrix to confirm our hypothesis. Significant correlations in anterior part of the right putamen between verbal scores and surface metric were revealed in the 18 and 21 age groups.

Feature selective temporal prediction of Alzheimer's disease progression using hippocampus surface morphometry

INTRODUCTION

Prediction of Alzheimer's disease (AD) progression based on baseline measures allows us to understand disease progression and has implications in decisions concerning treatment strategy. To this end, we combine a predictive multi-task machine learning method (cFSGL) with a novel MR-based multivariate morphometric surface map of the hippocampus (mTBM) to predict future cognitive scores of patients.

METHODS

Previous work has shown that a multi-task learning framework that performs prediction of all future time points simultaneously (cFSGL) can be used to encode both sparsity as well as temporal smoothness. The authors showed that this method is able to predict cognitive outcomes of ADNI subjects using FreeSurfer-based baseline MRI features, MMSE score demographic information and ApoE status. Whilst volumetric information may hold generalized information on brain status, we hypothesized that hippocampus specific information may be more useful in predictive modeling of AD. To this end, we applied a multivariate tensor-based parametric surface analysis method (mTBM) to extract features from the hippocampal surfaces.

RESULTS

We combined mTBM features with traditional surface features such as middle axis distance, the Jacobian determinant as well as 2 of the Jacobian principal eigenvalues to yield 7 normalized hippocampal surface maps of 300 points each. By combining these 7 × 300 = 2100 features together with the previous ~350 features, we illustrate how this type of sparsifying method can be applied to an entire surface map of the hippocampus that yields a feature space that is 2 orders of magnitude larger than what was previously attempted.

CONCLUSIONS

By combining the power of the cFSGL multi-task machine learning framework with the addition of AD sensitive mTBM feature maps of the hippocampus surface, we are able to improve the predictive performance of ADAS cognitive scores 6, 12, 24, 36 and 48 months from baseline.

Enhancing Diffusion MRI Measures By Integrating Grey and White Matter Morphometry With Hyperbolic Wasserstein Distance

In order to improve the preclinical diagnose of Alzheimer's disease (AD), there is a great deal of interest in analyzing the AD related brain structural changes with magnetic resonance image (MRI) analyses. As the major features, variation of the structural connectivity and the cortical surface morphometry provide different views of structural changes to determine whether AD is present on presymptomatic patients. However, the large scale tensor-valued information and relatively low imaging resolution in diffusion MRI (dMRI) have created huge challenges for analysis. In this paper, we propose a novel framework that improves dMRI analysis power by fusing cortical surface morphometry features from structural MRI (sMRI). We first compute the hyperbolic harmonic maps between cortical surfaces with the landmark constraints thus to precisely evaluate surface tensor-based morphometry. Meanwhile, the graph-based analysis of structural connectivity derived from dMRI is conducted. Next, we fuse these two features via the optimal mass transportation (OMT) and eventually the Wasserstein distance (WD) based single image index is computed as a potential clinical multimodality imaging score. We apply our framework to brain images of 20 AD patients and 20 matched healthy controls, randomly chosen from the Alzheimer's Disease Neuroimaging Initiative (AD-NI2) dataset. Our preliminary experimental results of group classification outperformed those of some other single dMRI-based features, such as regional hippocampal volume, mean scores of fractional anisotropy (FA) and mean axial (MD). The novel image fusion pipeline and simple imaging score of structural changes may benefit the preclinical AD and AD prevention research.

Ventricular shape and relative position abnormalities in preterm neonates

Recent neuroimaging findings have highlighted the impact of premature birth on subcortical development and morphological changes in the deep grey nuclei and ventricular system. To help characterize subcortical microstructural changes in preterm neonates, we recently implemented a multivariate tensor-based method (mTBM). This method allows to precisely measure local surface deformation of brain structures in infants. Here, we investigated ventricular abnormalities and their spatial relationships with surrounding subcortical structures in preterm neonates. We performed regional group comparisons on the surface morphometry and relative position of the lateral ventricles between 19 full-term and 17 preterm born neonates at term-equivalent age. Furthermore, a relative pose analysis was used to detect individual differences in translation, rotation, and scale of a given brain structure with respect to an average. Our mTBM results revealed broad areas of alterations on the frontal horn and body of the left ventricle, and narrower areas of differences on the temporal horn of the right ventricle. A significant shift in the rotation of the left ventricle was also found in preterm neonates. Furthermore, we located significant correlations between morphology and pose parameters of the lateral ventricles and that of the putamen and thalamus. These results show that regional abnormalities on the surface and pose of the ventricles are also associated with alterations on the putamen and thalamus. The complementarity of the information provided by the surface and pose analysis may help to identify abnormal white and grey matter growth, hinting toward a pattern of neural and cellular dysmaturation.

Hyperbolic Harmonic Mapping for Surface Registration

Automatic computation of surface correspondence via harmonic map is an active research field in computer vision, computer graphics and computational geometry. It may help document and understand physical and biological phenomena and also has broad applications in biometrics, medical imaging and motion capture industries. Although numerous studies have been devoted to harmonic map research, limited progress has been made to compute a diffeomorphic harmonic map on general topology surfaces with landmark constraints. This work conquers this problem by changing the Riemannian metric on the target surface to a hyperbolic metric so that the harmonic mapping is guaranteed to be a diffeomorphism under landmark constraints. The computational algorithms are based on Ricci flow and nonlinear heat diffusion methods. The approach is general and robust. We employ our algorithm to study the constrained surface registration problem which applies to both computer vision and medical imaging applications. Experimental results demonstrate that, by changing the Riemannian metric, the registrations are always diffeomorphic and achieve relatively high performance when evaluated with some popular surface registration evaluation standards.

Generalized Scalar-on-Image Regression Models via Total Variation

The use of imaging markers to predict clinical outcomes can have a great impact in public health. The aim of this paper is to develop a class of generalized scalar-on-image regression models via total variation (GSIRM-TV), in the sense of generalized linear models, for scalar response and imaging predictor with the presence of scalar covariates. A key novelty of GSIRM-TV is that it is assumed that the slope function (or image) of GSIRM-TV belongs to the space of bounded total variation in order to explicitly account for the piecewise smooth nature of most imaging data. We develop an efficient penalized total variation optimization to estimate the unknown slope function and other parameters. We also establish nonasymptotic error bounds on the excess risk. These bounds are explicitly specified in terms of sample size, image size, and image smoothness. Our simulations demonstrate a superior performance of GSIRM-TV against many existing approaches. We apply GSIRM-TV to the analysis of hippocampus data obtained from the Alzheimers Disease Neuroimaging Initiative (ADNI) dataset.

Longitudinal Analysis of Pre-Term Neonatal Cerebral Ventricles From 3D Ultrasound Images Using Spatial-Temporal Deformable Registration

Preterm neonates with a very low birth weight of less than 1,500 grams are at increased risk for developing intraventricular hemorrhage (IVH), which is a major cause of brain injury in preterm neonates. Quantitative measurements of ventricular dilatation or shrinkage play an important role in monitoring patients and evaluating treatment options. 3D ultrasound (US) has been developed to monitor ventricle volume as a biomarker for ventricular changes. However, ventricle volume as a global indicator does not allow for precise analysis of local ventricular changes, which could be linked to specific neurological problems often seen in the patient population later in life. In this work, a 3D+t spatial-temporal deformable registration approachis proposed, which is applied to the analysis of the detailed local changes of preterm IVH neonatal ventricles from 3D US images. In particular, a novel sequential convex/dual optimization algorithm is introduced to extract the optimal 3D+t spatial-temporal deformable field, which simultaneously optimizes the sequence of 3D deformation fieldswhile enjoying both efficiencyand simplicity in numerics. The developed registration technique was evaluated by comparing two manually extracted ventricle surfaces from the baseline and the registered follow-up images using the metrics of Dice similarity coefficient (DSC), mean absolute surface distance (MAD), and maximum absolute surface distance (MAXD). The performed experiments using 14 patients with 5 time-point images per patient show that the proposed 3D+t registration approach accurately recovered the longitudinal deformation of ventricle surfaces from 3D US images. The proposed approach may be potentially used to analyse the change pattern of cerebral ventricles of IVH patients, their response to different treatment options, and to elucidate the deficiencies that a patient could have later in life. To the best of our knowledge, this paper reports the first study on the longitudinalanalysis of neonatal ventricular system from 3D US images.

An Optimal Transportation based Univariate Neuroimaging Index

The alterations of brain structures and functions have been considered closely correlated to the change of cognitive performance due to neurodegenerative diseases such as Alzheimer's disease. In this paper, we introduce a variational framework to compute the optimal transformation (OT) in 3D space and propose a univariate neuroimaging index based on OT to measure such alterations. We compute the OT from each image to a template and measure the Wasserstein distance between them. By comparing the distances from all the images to the common template, we obtain a concise and informative index for each image. Our framework makes use of the Newton's method, which reduces the computational cost and enables itself to be applicable to large-scale datasets. The proposed work is a generic approach and thus may be applicable to various volumetric brain images, including structural magnetic resonance (sMR) and fluorodeoxyglucose positron emission tomography (FDG-PET) images. In the classification between Alzheimer's disease patients and healthy controls, our method achieves an accuracy of 82.30% on the Alzheimers Disease Neuroimaging Initiative (ADNI) baseline sMRI dataset and outperforms several other indices. On FDG-PET dataset, we boost the accuracy to 88.37% by leveraging pairwise Wasserstein distances. In a longitudinal study, we obtain a 5% significance with p-value = 1.13×10⁵ in a t-test on FDG-PET. The results demonstrate a great potential of the proposed index for neuroimage analysis and the precision medicine research.

Multivariate surface-based analysis of corpus callosum in patients with sickle cell disease

Sickle cell disease (SCD) is a genetic hematological disease in which the hemoglobin molecule in red blood cells is abnormal. It is closely associated with many symptoms, including pain, anemia, chest syndrome and neurocognitive impairment. One of the most debilitating symptoms is elevated risk for cerebro-vascular accidents. The corpus callosum (CC), as the largest and most prominent white matter (WM) structure in the brain, can reflect the chronic cerebrovascular damage resulting from silent strokes or infarctions in asymptomatic SCD patients. While a lot of studies have reported WM alterations in this cohort, little is known about the shape deformation of the CC. Here we perform the first surface morphometry analysis of the CC in SCD patients using four different shape metrics on T1-weighted magnetic resonance images. We detect regional surface morphological differences in the CC between 11 patients and 10 healthy control subjects. Differences are located in the genu, posterior midbody and splenium, potentially casting light on the anatomical substrates underlying neuropsychological test differences between the SCD and control groups.

Putamen Development in Children 12 to 21 Months Old

We studied the developmental trajectory of the putamen in 13-21 months old children using multivariate surface tensor-based morphometry. Our results indicate surface changes between 12 and 15 months' age groups in the middle superior part the left putamen. The growth of the left putamen at earlier ages slows down after 15 months. The most important surface changes were detected in the right putamen between 18 and 21 months and were located in the anterior part of the structure. Our results demonstrate the heterochronic growth of the right and left putamen related to different functional subregions within putamen. Our results are compatible with previous studies devoted to total putamen volume changes during normal development.

Predicting Alzheimer's Disease Using Combined Imaging-Whole Genome SNP Data

The growing public threat of Alzheimer's disease (AD) has raised the urgency to discover and validate prognostic biomarkers in order to predicting time to onset of AD. It is anticipated that both whole genome single nucleotide polymorphism (SNP) data and high dimensional whole brain imaging data offer predictive values to identify subjects at risk for progressing to AD. The aim of this paper is to test whether both whole genome SNP data and whole brain imaging data offer predictive values to identify subjects at risk for progressing to AD. In 343 subjects with mild cognitive impairment (MCI) enrolled in the Alzheimer's Disease Neuroimaging Initiative (ADNI-1), we extracted high dimensional MR imaging (volumetric data on 93 brain regions plus a surface fluid registration based hippocampal subregion and surface data), and whole genome data (504,095 SNPs from GWAS), as well as routine neurocognitive and clinical data at baseline. MCI patients were then followed over 48 months, with 150 participants progressing to AD. Combining information from whole brain MR imaging and whole genome data was substantially superior to the standard model for predicting time to onset of AD in a 48-month national study of subjects at risk. Our findings demonstrate the promise of combined imaging-whole genome prognostic markers in people with mild memory impairment.

Influence of APOE Genotype on Hippocampal Atrophy over Time - An N=1925 Surface-Based ADNI Study

The apolipoprotein E (APOE) e4 genotype is a powerful risk factor for late-onset Alzheimer’s disease (AD). In the Alzheimer’s Disease Neuroimaging Initiative (ADNI) cohort, we previously reported significant baseline structural differences in APOE e4 carriers relative to non-carriers, involving the left hippocampus more than the right—a difference more pronounced in e4 homozygotes than heterozygotes. We now examine the longitudinal effects of APOE genotype on hippocampal morphometry at 6-, 12- and 24-months, in the ADNI cohort. We employed a new automated surface registration system based on conformal geometry and tensor-based morphometry. Among different hippocampal surfaces, we computed high-order correspondences, using a novel inverse-consistent surface-based fluid registration method and multivariate statistics consisting of multivariate tensor-based morphometry (mTBM) and radial distance. At each time point, using Hotelling’s T² test, we found significant morphological deformation in APOE e4 carriers relative to non-carriers in the full cohort as well as in the non-demented (pooled MCI and control) subjects at each follow-up interval. In the complete ADNI cohort, we found greater atrophy of the left hippocampus than the right, and this asymmetry was more pronounced in e4 homozygotes than heterozygotes. These findings, combined with our earlier investigations, demonstrate an e4 dose effect on accelerated hippocampal atrophy, and support the enrichment of prevention trial cohorts with e4 carriers.

MORPHOMETRIC ANALYSIS OF HIPPOCAMPUS AND LATERAL VENTRICLE REVEALS REGIONAL DIFFERENCE BETWEEN COGNITIVELY STABLE AND DECLINING PERSONS

Alzheimers disease (AD) is a progressive neurodegenerative disease most prevalent in the elderly. Distinguishing disease-related memory decline from normal age-related memory decline has been clinically difficult due to the subtlety of cognitive change during the preclinical stage of AD. In contrast, sensitive biomarkers derived from in vivo neuroimaging data could improve the early identification of AD. In this study, we employed a morphometric analysis in the hippocampus and lateral ventricle. A novel group-wise template-based segmentation algorithm was developed for ventricular segmentation. Further, surface multivariate tensor-based morphometry and radial distance on each surface point were computed. Using Hotellings T² test, we found significant morphometric differences in both hippocampus and lateral ventricle between stable and clinically declining subjects. The left hemisphere was more severely affected than the right during this early disease stage. Hippocampal and ventricular morphometry has significant potential as an imaging biomarker for onset prediction and early diagnosis of AD.

APPLYING SPARSE CODING TO SURFACE MULTIVARIATE TENSOR-BASED MORPHOMETRY TO PREDICT FUTURE COGNITIVE DECLINE

Alzheimer's disease (AD) is a progressive brain disease. Accurate diagnosis of AD and its prodromal stage, mild cognitive impairment, is crucial for clinical trial design. There is also growing interests in identifying brain imaging biomarkers that help evaluate AD risk presymptomatically. Here, we applied a recently developed multivariate tensor-based morphometry (mTBM) method to extract features from hippocampal surfaces, derived from anatomical brain MRI. For such surface-based features, the feature dimension is usually much larger than the number of subjects. We used dictionary learning and sparse coding to effectively reduce the feature dimensions. With the new features, an Adaboost classifier was employed for binary group classification. In tests on publicly available data from the Alzheimers Disease Neuroimaging Initiative, the new framework outperformed several standard imaging measures in classifying different stages of AD. The new approach combines the efficiency of sparse coding with the sensitivity of surface mTBM, and boosts classification performance.

Impact of Early and Late Visual Deprivation on the Structure of the Corpus Callosum: A Study Combining Thickness Profile with Surface Tensor-Based Morphometry

Blindness represents a unique model to study how visual experience may shape the development of brain organization. Exploring how the structure of the corpus callosum (CC) reorganizes ensuing visual deprivation is of particular interest due to its important functional implication in vision (e.g., via the splenium of the CC). Moreover, comparing early versus late visually deprived individuals has the potential to unravel the existence of a sensitive period for reshaping the CC structure. Here, we develop a novel framework to capture a complete set of shape differences in the CC between congenitally blind (CB), late blind (LB) and sighted control (SC) groups. The CCs were manually segmented from T1-weighted brain MRI and modeled by 3D tetrahedral meshes. We statistically compared the combination of local area and thickness at each point between subject groups. Differences in area are found using surface tensor-based morphometry; thickness is estimated by tracing the streamlines in the volumetric harmonic field. Group differences were assessed on this combined measure using Hotelling's T(2) test. Interestingly, we observed that the total callosal volume did not differ between the groups. However, our fine-grained analysis reveals significant differences mostly localized around the splenium areas between both blind groups and the sighted group (general effects of blindness) and, importantly, specific dissimilarities between the LB and CB groups, illustrating the existence of a sensitive period for reorganization. The new multivariate statistics also gave better effect sizes for detecting morphometric differences, relative to other statistics. They may boost statistical power for CC morphometric analyses.

BFLCRM: A BAYESIAN FUNCTIONAL LINEAR COX REGRESSION MODEL FOR PREDICTING TIME TO CONVERSION TO ALZHEIMER'S DISEASE

The aim of this paper is to develop a Bayesian functional linear Cox regression model (BFLCRM) with both functional and scalar covariates. This new development is motivated by establishing the likelihood of conversion to Alzheimer's disease (AD) in 346 patients with mild cognitive impairment (MCI) enrolled in the Alzheimer's Disease Neuroimaging Initiative 1 (ADNI-1) and the early markers of conversion. These 346 MCI patients were followed over 48 months, with 161 MCI participants progressing to AD at 48 months. The functional linear Cox regression model was used to establish that functional covariates including hippocampus surface morphology and scalar covariates including brain MRI volumes, cognitive performance (ADAS-Cog), and APOE status can accurately predict time to onset of AD. Posterior computation proceeds via an efficient Markov chain Monte Carlo algorithm. A simulation study is performed to evaluate the finite sample performance of BFLCRM.

A framework for quantification and visualization of segmentation accuracy and variability in 3D lateral ventricle ultrasound images of preterm neonates

PURPOSE

Intraventricular hemorrhage (IVH) is a major cause of brain injury in preterm neonates. Three dimensional ultrasound (US) imaging systems have been developed to visualize 3D anatomical structure of preterm neonatal intracranial ventricular system with IVH and ventricular dilation. To allow quantitative analysis, the ventricle system is required to be segmented accurately and efficiently from 3D US images. Although semiautomatic segmentation algorithms have been developed, local segmentation accuracy and variability associated with these algorithms should be evaluated statistically before they can be applied in clinical settings. This work proposes a statistical framework to quantify the local accuracy and variability and performs statistical tests to identify locations where the semiautomatically segmented surfaces are significantly different from manually segmented surfaces.

METHODS

Three dimensional lateral ventricle US images of preterm neonates were each segmented six times manually and using a semiautomated segmentation algorithm. The local difference between manually and algorithmically segmented surfaces as well as the segmentation variability for each method was computed and superimposed on the ventricular surface of each subject. To summarize the segmentation performance for a whole group of subjects, the subject-specific local difference and standard deviation maps were registered onto a 3D template ventricular surface using a nonrigid registration algorithm. Pointwise, intersubject average accuracy and pooled variability for the whole group of subjects can be computed and visualized on the template surface, providing a summary of performance of the segmentation algorithm for the whole group of ventricles with highly variable geometry. In addition to pointwise statistical analysis performed on the template surface, statistical conclusion regarding the accuracy of the segmentation algorithm was made for subregions and the whole ventricle with the spatial correlation of pointwise accuracy taken into account.

RESULTS

Ten 3D US images were involved in this study. Pointwise local difference, ΔS, its absolute value |ΔS| as well as the standard deviations of the manual and algorithm segmentations were computed and superimposed on the each ventricle surface. Regions with lower segmentation accuracy and higher segmentation variability can be identified from these maps, and the localized information was applied to improve the accuracy of the algorithm. Intersubject average ΔS and |ΔS| as well as pooled standard deviations was computed on the template surface. Intersubject average ΔS and |ΔS| indicated that the algorithm underestimated regions in the neighborhood of the tips of anterior, inferior, and posterior horns. Intersubject pooled standard deviations indicated that manual segmentation had a higher segmentation variability than algorithm segmentation over the whole ventricle. Statistical analysis on the template surface showed that there was significant difference between algorithm and manual methods for segmenting the right lateral ventricle but not for the left lateral ventricle.

CONCLUSIONS

A framework was proposed for evaluating, visualizing, and summarizing the local accuracy and variability of a segmentation algorithm. This framework can be used for improving the accuracy of segmentation algorithms, as well as providing useful feedback to improve the manual segmentation performance. More importantly, this framework can be applied for longitudinal monitoring of local ventricular changes of neonates with IVH.

Novel Neuroimaging Methods to Understand How HIV Affects the Brain

In much of the developed world, the HIV epidemic has largely been controlled by antiretroviral treatment. Even so, there is growing concern that HIV-infected individuals may be at risk for accelerated brain aging and a range of cognitive impairments. What promotes or resists these changes is largely unknown. There is also interest in discovering factors that promote resilience to HIV and combat its adverse effects in children. Here, we review recent developments in brain imaging that reveal how the virus affects the brain. We relate these brain changes to changes in blood markers, cognitive function, and other patient outcomes or symptoms, such as apathy or neuropathic pain. We focus on new and emerging techniques, including new variants of brain MRI. Diffusion tensor imaging, for example, can map the brain's structural connections, while fMRI can uncover functional connections. Finally, we suggest how large-scale global research alliances, such as ENIGMA, may resolve controversies over effects where evidence is now lacking. These efforts pool scans from tens of thousands of individuals and offer a source of power not previously imaginable for brain imaging studies.

Unified heat kernel regression for diffusion, kernel smoothing and wavelets on manifolds and its application to mandible growth modeling in CT images

We present a novel kernel regression framework for smoothing scalar surface data using the Laplace-Beltrami eigenfunctions. Starting with the heat kernel constructed from the eigenfunctions, we formulate a new bivariate kernel regression framework as a weighted eigenfunction expansion with the heat kernel as the weights. The new kernel method is mathematically equivalent to isotropic heat diffusion, kernel smoothing and recently popular diffusion wavelets. The numerical implementation is validated on a unit sphere using spherical harmonics. As an illustration, the method is applied to characterize the localized growth pattern of mandible surfaces obtained in CT images between ages 0 and 20 by regressing the length of displacement vectors with respect to a surface template.