Multi-organ imaging-derived polygenic indexes for brain and body health.

The UK Biobank (UKB) imaging project is a crucial resource for biomedical research, but is limited to 100,000 participants due to cost and accessibility barriers. Here we used genetic data to predict heritable imaging-derived phenotypes (IDPs) for a larger cohort. We developed and evaluated 4,375 IDP genetic scores (IGS) derived from UKB brain and body images. When applied to UKB participants who were not imaged, IGS revealed links to numerous phenotypes and stratified participants at increased risk for both brain and somatic diseases. For example, IGS identified individuals at higher risk for Alzheimer's disease and multiple sclerosis, offering additional insights beyond traditional polygenic risk scores of these diseases. When applied to independent external cohorts, IGS also stratified those at high disease risk in the All of Us Research Program and the Alzheimer's Disease Neuroimaging Initiative study. Our results demonstrate that, while the UKB imaging cohort is largely healthy and may not be the most enriched for disease risk management, it holds immense potential for stratifying the risk of various brain and body diseases in broader external genetic cohorts.

Merging or ensembling: integrative analysis in multiple neuroimaging studies.

The aim of this paper is to systematically investigate merging and ensembling methods for spatially varying coefficient mixed effects models (SVCMEM) in order to carry out integrative learning of neuroimaging data obtained from multiple biomedical studies. The "merged" approach involves training a single learning model using a comprehensive dataset that encompasses information from all the studies. Conversely, the "ensemble" approach involves creating a weighted average of distinct learning models, each developed from an individual study. We systematically investigate the prediction accuracy of the merged and ensemble learners under the presence of different degrees of interstudy heterogeneity. Additionally, we establish asymptotic guidelines for making strategic decisions about when to employ either of these models in different scenarios, along with deriving optimal weights for the ensemble learner. To validate our theoretical results, we perform extensive simulation studies. The proposed methodology is also applied to 3 large-scale neuroimaging studies.

Source-free unsupervised domain adaptation: A survey.

Unsupervised domain adaptation (UDA) via deep learning has attracted appealing attention for tackling domain-shift problems caused by distribution discrepancy across different domains. Existing UDA approaches highly depend on the accessibility of source domain data, which is usually limited in practical scenarios due to privacy protection, data storage and transmission cost, and computation burden. To tackle this issue, many source-free unsupervised domain adaptation (SFUDA) methods have been proposed recently, which perform knowledge transfer from a pre-trained source model to the unlabeled target domain with source data inaccessible. A comprehensive review of these works on SFUDA is of great significance. In this paper, we provide a timely and systematic literature review of existing SFUDA approaches from a technical perspective. Specifically, we categorize current SFUDA studies into two groups, i.e., white-box SFUDA and black-box SFUDA, and further divide them into finer subcategories based on different learning strategies they use. We also investigate the challenges of methods in each subcategory, discuss the advantages/disadvantages of white-box and black-box SFUDA methods, conclude the commonly used benchmark datasets, and summarize the popular techniques for improved generalizability of models learned without using source data. We finally discuss several promising future directions in this field.

FPLS-DC: functional partial least squares through distance covariance for imaging genetics.

MOTIVATION: Imaging genetics integrates imaging and genetic techniques to examine how genetic variations influence the function and structure of organs like the brain or heart, providing insights into their impact on behavior and disease phenotypes. The use of organ-wide imaging endophenotypes has increasingly been used to identify potential genes associated with complex disorders. However, analyzing organ-wide imaging data alongside genetic data presents two significant challenges: high dimensionality and complex relationships. To address these challenges, we propose a novel, nonlinear inference framework designed to partially mitigate these issues. RESULTS: We propose a functional partial least squares through distance covariance (FPLS-DC) framework for efficient genome wide analyses of imaging phenotypes. It consists of two components. The first component utilizes the FPLS-derived base functions to reduce image dimensionality while screening genetic markers. The second component maximizes the distance correlation between genetic markers and projected imaging data, which is a linear combination of the FPLS-basis functions, using simulated annealing algorithm. In addition, we proposed an iterative FPLS-DC method based on FPLS-DC framework, which effectively overcomes the influence of inter-gene correlation on inference analysis. We efficiently approximate the null distribution of test statistics using a gamma approximation. Compared to existing methods, FPLS-DC offers computational and statistical efficiency for handling large-scale imaging genetics. In real-world applications, our method successfully detected genetic variants associated with the hippocampus, demonstrating its value as a statistical toolbox for imaging genetic studies. AVAILABILITY AND IMPLEMENTATION: The FPLS-DC method we propose opens up new research avenues and offers valuable insights for analyzing functional and high-dimensional data. In addition, it serves as a useful tool for scientific analysis in practical applications within the field of imaging genetics research. The R package FPLS-DC is available in Github: https://github.com/BIG-S2/FPLSDC.

A comprehensive patient-specific prediction model for temporomandibular joint osteoarthritis progression.

Temporomandibular joint osteoarthritis (TMJ OA) is a prevalent degenerative disease characterized by chronic pain and impaired jaw function. The complexity of TMJ OA has hindered the development of prognostic tools, posing a significant challenge in timely, patient-specific management. Addressing this gap, our research employs a comprehensive, multidimensional approach to advance TMJ OA prognostication. We conducted a prospective study with 106 subjects, 74 of whom were followed up after 2 to 3 y of conservative treatment. Central to our methodology is the development of an innovative, open-source predictive modeling framework, the Ensemble via Hierarchical Predictions through Nested cross-validation tool (EHPN). This framework synergistically integrates 18 feature selection, statistical, and machine learning methods to yield an accuracy of 0.87, with an area under the ROC curve of 0.72 and an F1 score of 0.82. Our study, beyond technical advancements, emphasizes the global impact of TMJ OA, recognizing its unique demographic occurrence. We highlight key factors influencing TMJ OA progression. Using SHAP analysis, we identified personalized prognostic predictors: lower values of headache, lower back pain, restless sleep, condyle high gray level-GL-run emphasis, articular fossa GL nonuniformity, and long-run low GL emphasis; and higher values of superior joint space, mouth opening, saliva Vascular-endothelium-growth-factor, Matrix-metalloproteinase-7, serum Epithelial-neutrophil-activating-peptide, and age indicate recovery likelihood. Our multidimensional and multimodal EHPN tool enhances clinicians' decision-making, offering a transformative translational infrastructure. The EHPN model stands as a significant contribution to precision medicine, offering a paradigm shift in the management of temporomandibular disorders and potentially influencing broader applications in personalized healthcare.

Leveraging Brain Modularity Prior for Interpretable Representation Learning of fMRI.

Resting-state functional magnetic resonance imaging (rs-fMRI) can reflect spontaneous neural activities in the brain and is widely used for brain disorder analysis. Previous studies focus on extracting fMRI representations using machine/deep learning methods, but these features typically lack biological interpretability. The human brain exhibits a remarkable modular structure in spontaneous brain functional networks, with each module comprised of functionally interconnected brain regions-of-interest (ROIs). However, existing learning-based methods cannot adequately utilize such brain modularity prior. In this paper, we propose a brain modularity-constrained dynamic representation learning framework for interpretable fMRI analysis, consisting of dynamic graph construction, dynamic graph learning via a novel modularity-constrained graph neural network (MGNN), and prediction and biomarker detection. The designed MGNN is constrained by three core neurocognitive modules (i.e., salience network, central executive network, and default mode network), encouraging ROIs within the same module to share similar representations. To further enhance discriminative ability of learned features, we encourage the MGNN to preserve network topology of input graphs via a graph topology reconstruction constraint. Experimental results on 534 subjects with rs-fMRI scans from two datasets validate the effectiveness of the proposed method. The identified discriminative brain ROIs and functional connectivities can be regarded as potential fMRI biomarkers to aid in clinical diagnosis.

Identification and validation of supervariants reveal novel loci associated with human white matter microstructure.

As an essential part of the central nervous system, white matter coordinates communications between different brain regions and is related to a wide range of neurodegenerative and neuropsychiatric disorders. Previous genome-wide association studies (GWASs) have uncovered loci associated with white matter microstructure. However, GWASs suffer from limited reproducibility and difficulties in detecting multi-single-nucleotide polymorphism (multi-SNP) and epistatic effects. In this study, we adopt the concept of supervariants, a combination of alleles in multiple loci, to account for potential multi-SNP effects. We perform supervariant identification and validation to identify loci associated with 22 white matter fractional anisotropy phenotypes derived from diffusion tensor imaging. To increase reproducibility, we use United Kingdom (UK) Biobank White British (n = 30,842) data for discovery and internal validation, and UK Biobank White but non-British (n = 1927) data, Europeans from the Adolescent Brain Cognitive Development study (n = 4399) data, and Europeans from the Human Connectome Project (n = 319) data for external validation. We identify 23 novel loci on the discovery set that have not been reported in the previous GWASs on white matter microstructure. Among them, three supervariants on genomic regions 5q35.1, 8p21.2, and 19q13.32 have P-values lower than 0.05 in the meta-analysis of the three independent validation data sets. These supervariants contain genetic variants located in genes that have been related to brain structures, cognitive functions, and neuropsychiatric diseases. Our findings provide a better understanding of the genetic architecture underlying white matter microstructure.

The interaction effects of age, APOE and common environmental risk factors on human brain structure.

Mounting evidence suggests considerable diversity in brain aging trajectories, primarily arising from the complex interplay between age, genetic, and environmental risk factors, leading to distinct patterns of micro- and macro-cerebral aging. The underlying mechanisms of such effects still remain unclear. We conducted a comprehensive association analysis between cerebral structural measures and prevalent risk factors, using data from 36,969 UK Biobank subjects aged 44-81. Participants were assessed for brain volume, white matter diffusivity, Apolipoprotein E (APOE) genotypes, polygenic risk scores, lifestyles, and socioeconomic status. We examined genetic and environmental effects and their interactions with age and sex, and identified 726 signals, with education, alcohol, and smoking affecting most brain regions. Our analysis revealed negative age-APOE-ε4 and positive age-APOE-ε2 interaction effects, respectively, especially in females on the volume of amygdala, positive age-sex-APOE-ε4 interaction on the cerebellar volume, positive age-excessive-alcohol interaction effect on the mean diffusivity of the splenium of the corpus callosum, positive age-healthy-diet interaction effect on the paracentral volume, and negative APOE-ε4-moderate-alcohol interaction effects on the axial diffusivity of the superior fronto-occipital fasciculus. These findings highlight the need of considering age, sex, genetic, and environmental joint effects in elucidating normal or abnormal brain aging.

The Association of Cardiovascular and Neurological Comorbidities in Geriatric Patients Sustaining Ocular Trauma.

Purpose: Falls are associated with ocular trauma in the elderly. However, it is unlikely for a fall to cause ocular injury unless there is a disruption in the protective maneuvers that shield the face. We suspect ocular injury may be an early indicator of systemic or neurologic degeneration. This study investigates the 5-year incidence of cardiovascular and neurodegenerative diseases in older patients who sustained ocular or periorbital injuries. Patients and Methods: This was a retrospective cohort study. The study population included 141 patients over the age of 65 who sustained trauma to the eye, orbit, or eyelid between April 2011 and June 2016. The control population included 141 patients with a similar range of comorbidities who received cataract surgery during the same period. The study measured new diagnoses of various disorders during the 5-year period following presentation. Results: There were a total of 180 females and 102 males in the study. The mean ages of the control and subject group were 76 and 81.8, respectively. Of our twelve tested comorbidity types, patients that suffered a periocular trauma were more likely to develop heart failure (p=0.00244), dementia (p=0.00002), Alzheimer's disease (p=0.00087), and vascular disease (p=0.00037). Conclusion: Geriatric patients who sustained ocular and periocular injuries had a greater incidence of heart failure, dementia, Alzheimer's disease, and atherosclerosis diagnoses in the 5-year period following injury. The findings of this study suggest that periocular trauma may be an early indicator of underlying degenerative or systemic disease. Ophthalmologists should ensure proper primary care follow-up in conjunction with recovery from injury.

Bayesian bi-level variable selection for genome-wide survival study.

Mild cognitive impairment (MCI) is a clinical syndrome characterized by the onset and evolution of cognitive impairments, often considered a transitional stage to Alzheimer's disease (AD). The genetic traits of MCI patients who experience a rapid progression to AD can enhance early diagnosis capabilities and facilitate drug discovery for AD. While a genome-wide association study (GWAS) is a standard tool for identifying single nucleotide polymorphisms (SNPs) related to a disease, it fails to detect SNPs with small effect sizes due to stringent control for multiple testing. Additionally, the method does not consider the group structures of SNPs, such as genes or linkage disequilibrium blocks, which can provide valuable insights into the genetic architecture. To address the limitations, we propose a Bayesian bi-level variable selection method that detects SNPs associated with time of conversion from MCI to AD. Our approach integrates group inclusion indicators into an accelerated failure time model to identify important SNP groups. Additionally, we employ data augmentation techniques to impute censored time values using a predictive posterior. We adapt Dirichlet-Laplace shrinkage priors to incorporate the group structure for SNP-level variable selection. In the simulation study, our method outperformed other competing methods regarding variable selection. The analysis of Alzheimer's Disease Neuroimaging Initiative (ADNI) data revealed several genes directly or indirectly related to AD, whereas a classical GWAS did not identify any significant SNPs.

The pivotal role of the X-chromosome in the genetic architecture of the human brain.

Genes on the X-chromosome are extensively expressed in the human brain, resulting in substantial influences on brain development, intellectual disability, and other brain-related disorders. To comprehensively investigate the X-chromosome's impact on the cerebral cortex, white matter tract microstructures, and intrinsic and extrinsic brain functions, we examined 2,822 complex brain imaging traits obtained from n=34,000 subjects in the UK Biobank. We unveiled potential autosome-X-chromosome interaction, while proposing an atlas of dosage compensation (DC) for each set of traits. We observed a pronounced X-chromosome impact on the corticospinal tract and the functional amplitude and connectivity of visual networks. In association studies, we identified 50 genome-wide significant trait-locus pairs enriched in Xq28, 22 of which replicated in independent datasets (n=4,900). Notably, 13 newly identified pairs were in the X-chromosome's non-pseudo-autosomal regions (NPR). The volume of the right ventral diencephalon shared genetic architecture with schizophrenia and educational attainment in a locus indexed by rs2361468 (located ~3kb upstream of PJA1, a conserved and ubiquitously expressed gene implicated in multiple psychiatric disorders). No significant associations were identified in the pseudo-autosomal regions (PAR) or the Y-chromosome. Finally, we explored sex-specific associations on the X-chromosome and compared differing genetic effects between sexes. We found much more associations can be identified in males (33 versus 9) given a similar sample size. In conclusion, our research provides invaluable insights into the X-chromosome's role in the human brain, contributing to the observed sex differences in brain structure and function.

Charting Aging Trajectories of Knee Cartilage Thickness for Early Osteoarthritis Risk Prediction: An MRI Study from the Osteoarthritis Initiative Cohort.

Knee osteoarthritis (OA), a prevalent joint disease in the U.S., poses challenges in terms of predicting of its early progression. Although high-resolution knee magnetic resonance imaging (MRI) facilitates more precise OA diagnosis, the heterogeneous and multifactorial aspects of OA pathology remain significant obstacles for prognosis. MRI-based scoring systems, while standardizing OA assessment, are both time-consuming and labor-intensive. Current AI technologies facilitate knee OA risk scoring and progression prediction, but these often focus on the symptomatic phase of OA, bypassing initial-stage OA prediction. Moreover, their reliance on complex algorithms can hinder clinical interpretation. To this end, we make this effort to construct a computationally efficient, easily-interpretable, and state-of-the-art approach aiding in the radiographic OA (rOA) auto-classification and prediction of the incidence and progression, by contrasting an individual's cartilage thickness with a similar demographic in the rOA-free cohort. To better visualize, we have developed the toolset for both prediction and local visualization. A movie demonstrating different subtypes of dynamic changes in local centile scores during rOA progression is available at https://tli3.github.io/KneeOA/. Specifically, we constructed age-BMI-dependent reference charts for knee OA cartilage thickness, based on MRI scans from 957 radiographic OA (rOA)-free individuals from the Osteoarthritis Initiative cohort. Then we extracted local and global centiles by contrasting an individual's cartilage thickness to the rOA-free cohort with a similar age and BMI. Using traditional boosting approaches with our centile-based features, we obtain rOA classification of KLG ≤ 1 versus KLG = 2 (AUC = 0.95, F1 = 0.89), KLG ≤ 1 versus KLG ≥ 2 (AUC = 0.90, F1 = 0.82) and prediction of KLG2 progression (AUC = 0.98, F1 = 0.94), rOA incidence (KLG increasing from < 2 to ≥ 2; AUC = 0.81, F1 = 0.69) and rOA initial transition (KLG from 0 to 1; AUC = 0.64, F1 = 0.65) within a future 48-month period. Such performance in classifying KLG ≥ 2 matches that of deep learning methods in recent literature. Furthermore, its clinical interpretation suggests that cartilage changes, such as thickening in lateral femoral and anterior femoral regions and thinning in lateral tibial regions, may serve as indicators for prediction of rOA incidence and early progression. Meanwhile, cartilage thickening in the posterior medial and posterior lateral femoral regions, coupled with a reduction in the central medial femoral region, may signify initial phases of rOA transition.

Interactions between Bifidobacterium and Bacteroides and human milk oligosaccharides and their associations with infant cognition.

While ample research on independent associations between infant cognition and gut microbiota composition and human milk (HM) oligosaccharides (HMOs) has been reported, studies on how the interactions between gut microbiota and HMOs may yield associations with cognitive development in infancy are lacking. We aimed to determine how HMOs and species of Bacteroides and Bifidobacterium genera interact with each other and their associations with cognitive development in typically developing infants. A total of 105 mother-infant dyads were included in this study. The enrolled infants [2.9-12 months old (8.09 ± 2.48)] were at least predominantly breastfed at 4 months old. A total of 170 HM samples from the mothers and fecal samples of the children were collected longitudinally. Using the Mullen Scales of Early Learning to assess cognition and the scores as the outcomes, linear mixed effects models including both the levels of eight HMOs and relative abundance of Bacteroides and Bifidobacterium species as main associations and their interactions were employed with adjusting covariates; infant sex, delivery mode, maternal education, site, and batch effects of HMOs. Additionally, regression models stratifying infants based on the A-tetrasaccharide (A-tetra) status of the HM they received were also employed to determine if the associations depend on the A-tetra status. With Bacteroides species, we observed significant associations with motor functions, while Bif. catenulatum showed a negative association with visual reception in the detectable A-tetra group both as main effect (value of p = 0.012) and in interaction with LNFP-I (value of p = 0.007). Additionally, 3-FL showed a positive association with gross motor (p = 0.027) and visual reception (p = 0.041). Furthermore, significant associations were observed with the interaction terms mainly in the undetectable A-tetra group. Specifically, we observed negative associations for Bifidobacterium species and LNT [breve (p = 0.011) and longum (p = 0.022)], and positive associations for expressive language with 3'-SL and Bif. bifidum (p = 0.01), 6'-SL and B. fragilis (p = 0.019), and LNFP-I and Bif. kashiwanohense (p = 0.048), respectively. Our findings suggest that gut microbiota and HMOs are both independently and interactively associated with early cognitive development. In particular, the diverse interactions between HMOs and Bacteroides and Bifidobacterium species reveal different candidate pathways through which HMOs, Bifidobacterium and Bacteroides species potentially interact to impact cognitive development in infancy.

Lysobacter changpingensis sp. nov., a novel species of the genus Lysobacter isolated from a rhizosphere soil of strawberry in China.

In the present work, we characterized in detail strain CM-3-T8T, which was isolated from the rhizosphere soil of strawberries in Beijing, China, in order to elucidate its taxonomic position. Cells of strain CM-3-T8T were Gram-negative, non-spore-forming, aerobic, short rod. Growth occurred at 25-37 °C, pH 5.0-10.0, and in the presence of 0-8% (w/v) NaCl. Phylogenetic analysis based on 16S rRNA gene sequences revealed that strain CM-3-T8T formed a stable clade with Lysobacter soli DCY21T and Lysobacter panacisoli CJ29T, with the 16S rRNA gene sequence similarities of 98.91% and 98.50%. The average nucleotide identity and digital DNA-DNA hybridization values between strain SG-8 T and the two reference type strains listed above were 76.3%, 79.6%, and 34.3%, 27%, respectively. The DNA G + C content was 68.4% (mol/mol). The major cellular fatty acids were comprised of C15:0 iso (36.15%), C17:0 iso (8.40%), and C11:0 iso 3OH (8.28%). The major quinone system was ubiquinone Q-8. The major polar lipids were phosphatidylethanolamine (PE), phosphatidylethanolamine (PME), diphosphatidylglycerol (DPG), and aminophospholipid (APL). On the basis of phenotypic, genotypic, and phylogenetic evidence, strain CM-3-T8T (= ACCC 61714 T = JCM 34576 T) represents a new species within the genus Lysobacter, for which the name Lysobacter changpingensis sp. nov. is proposed.

Heart-brain connections: Phenotypic and genetic insights from magnetic resonance images.

Cardiovascular health interacts with cognitive and mental health in complex ways, yet little is known about the phenotypic and genetic links of heart-brain systems. We quantified heart-brain connections using multiorgan magnetic resonance imaging (MRI) data from more than 40,000 subjects. Heart MRI traits displayed numerous association patterns with brain gray matter morphometry, white matter microstructure, and functional networks. We identified 80 associated genomic loci (P < 6.09 × 10-10) for heart MRI traits, which shared genetic influences with cardiovascular and brain diseases. Genetic correlations were observed between heart MRI traits and brain-related traits and disorders. Mendelian randomization suggests that heart conditions may causally contribute to brain disorders. Our results advance a multiorgan perspective on human health by revealing heart-brain connections and shared genetic influences.

Mapping Genetic Topography of Cortical Thickness and Surface Area in Neonatal Brains.

Adult twin neuroimaging studies have revealed that cortical thickness (CT) and surface area (SA) are differentially influenced by genetic information, leading to their spatially distinct genetic patterning and topography. However, the postnatal origins of the genetic topography of CT and SA remain unclear, given the dramatic cortical development from neonates to adults. To fill this critical gap, this study unprecedentedly explored how genetic information differentially regulates the spatial topography of CT and SA in the neonatal brain by leveraging brain magnetic resonance (MR) images from 202 twin neonates with minimal influence by the complicated postnatal environmental factors. We capitalized on infant-dedicated computational tools and a data-driven spectral clustering method to parcellate the cerebral cortex into a set of distinct regions purely according to the genetic correlation of cortical vertices in terms of CT and SA, respectively, and accordingly created the first genetically informed cortical parcellation maps of neonatal brains. Both genetic parcellation maps exhibit bilaterally symmetric and hierarchical patterns, but distinct spatial layouts. For CT, regions with closer genetic relationships demonstrate an anterior-posterior (A-P) division, while for SA, regions with greater genetic proximity are typically within the same lobe. Certain genetically informed regions exhibit strong similarities between neonates and adults, with the most striking similarities in the medial surface in terms of SA, despite their overall substantial differences in genetic parcellation maps. These results greatly advance our understanding of the development of genetic influences on the spatial patterning of cortical morphology.SIGNIFICANCE STATEMENT Genetic influences on cortical thickness (CT) and surface area (SA) are complex and could evolve throughout the lifespan. However, studies revealing distinct genetic topography of CT and SA have been limited to adults. Using brain structural magnetic resonance (MR) images of twins, we unprecedentedly discovered the distinct genetically-informed parcellation maps of CT and SA in neonatal brains, respectively. Each genetic parcellation map comprises a distinct spatial layout of cortical regions, where vertices within the same region share high genetic correlation. These genetic parcellation maps of CT and SA of neonates largely differ from those of adults, despite their highly remarkable similarities in the medial cortex of SA. These discoveries provide important insights into the genetic organization of the early cerebral cortex development.

Four-dimensional mapping of dynamic longitudinal brain subcortical development and early learning functions in infants.

Brain subcortical structures are paramount in many cognitive functions and their aberrations during infancy are predisposed to various neurodevelopmental and neuropsychiatric disorders, making it highly essential to characterize the early subcortical normative growth patterns. This study investigates the volumetric development and surface area expansion of six subcortical structures and their associations with Mullen scales of early learning by leveraging 513 high-resolution longitudinal MRI scans within the first two postnatal years. Results show that (1) each subcortical structure (except for the amygdala with an approximately linear increase) undergoes rapid nonlinear volumetric growth after birth, which slows down at a structure-specific age with bilaterally similar developmental patterns; (2) Subcortical local area expansion reveals structure-specific and spatiotemporally heterogeneous patterns; (3) Positive associations between thalamus and both receptive and expressive languages and between caudate and putamen and fine motor are revealed. This study advances our understanding of the dynamic early subcortical developmental patterns.

Statistical Learning Methods for Neuroimaging Data Analysis with Applications.

The aim of this review is to provide a comprehensive survey of statistical challenges in neuroimaging data analysis, from neuroimaging techniques to large-scale neuroimaging studies and statistical learning methods. We briefly review eight popular neuroimaging techniques and their potential applications in neuroscience research and clinical translation. We delineate four themes of neuroimaging data and review major image processing analysis methods for processing neuroimaging data at the individual level. We briefly review four large-scale neuroimaging-related studies and a consortium on imaging genomics and discuss four themes of neuroimaging data analysis at the population level. We review nine major population-based statistical analysis methods and their associated statistical challenges and present recent progress in statistical methodology to address these challenges.

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.

Addressing multi-site functional MRI heterogeneity through dual-expert collaborative learning for brain disease identification.

Several studies employ multi-site rs-fMRI data for major depressive disorder (MDD) identification, with a specific site as the to-be-analyzed target domain and other site(s) as the source domain. But they usually suffer from significant inter-site heterogeneity caused by the use of different scanners and/or scanning protocols and fail to build generalizable models that can well adapt to multiple target domains. In this article, we propose a dual-expert fMRI harmonization (DFH) framework for automated MDD diagnosis. Our DFH is designed to simultaneously exploit data from a single labeled source domain/site and two unlabeled target domains for mitigating data distribution differences across domains. Specifically, the DFH consists of a domain-generic student model and two domain-specific teacher/expert models that are jointly trained to perform knowledge distillation through a deep collaborative learning module. A student model with strong generalizability is finally derived, which can be well adapted to unseen target domains and analysis of other brain diseases. To the best of our knowledge, this is among the first attempts to investigate multi-target fMRI harmonization for MDD diagnosis. Comprehensive experiments on 836 subjects with rs-fMRI data from 3 different sites show the superiority of our method. The discriminative brain functional connectivities identified by our method could be regarded as potential biomarkers for fMRI-related MDD diagnosis.