A technical review of canonical correlation analysis for neuroscience applications

Integrating host and microbiome biology using holo-omics

Holo-omics is the use of omics data to study a host and its inherent microbiomes - a biological system known as a "holobiont". A microbiome that exists in such a space often encounters habitat stability and in return provides metabolic capacities that can benefit their host. Here we present an overview of beneficial host-microbiome systems and propose and discuss several methodological frameworks that can be used to investigate the intricacies of the many as yet undefined host-microbiome interactions that influence holobiont homeostasis. While this is an emerging field, we anticipate that ongoing methodological advancements will enhance the biological resolution that is necessary to improve our understanding of host-microbiome interplay to make meaningful interpretations and biotechnological applications.

Slow and steady: auditory features for discriminating animal vocalizations

We propose that listeners can use temporal regularities - spectro-temporal correlations that change smoothly over time - to discriminate animal vocalizations within and between species. To test this idea, we used Slow Feature Analysis (SFA) to find the most temporally regular components of vocalizations from birds (blue jay, house finch, American yellow warbler, and great blue heron), humans (English speakers), and rhesus macaques. We projected vocalizations into the learned feature space and tested intra-class (same speaker/species) and inter-class (different speakers/species) auditory discrimination by a trained classifier. We found that: 1) Vocalization discrimination was excellent (> 95%) in all cases; 2) Performance depended primarily on the ~10 most temporally regular features; 3) Most vocalizations are dominated by ~10 features with high temporal regularity; and 4) These regular features are highly correlated with the most predictable components of animal sounds.

Data Efficiency, Dimensionality Reduction, and the Generalized Symmetric Information Bottleneck

The symmetric information bottleneck (SIB), an extension of the more familiar information bottleneck, is a dimensionality-reduction technique that simultaneously compresses two random variables to preserve information between their compressed versions. We introduce the generalized symmetric information bottleneck (GSIB), which explores different functional forms of the cost of such simultaneous reduction. We then explore the data set size requirements of such simultaneous compression. We do this by deriving bounds and root-mean-squared estimates of statistical fluctuations of the involved loss functions. We show that in typical situations, the simultaneous GSIB compression requires qualitatively less data to achieve the same errors compared to compressing variables one at a time. We suggest that this is an example of a more general principle that simultaneous compression is more data efficient than independent compression of each of the input variables.

Optimizing spatial specificity and signal quality in fNIRS: an overview of potential challenges and possible options for improving the reliability of real-time applications

The optical brain imaging method functional near-infrared spectroscopy (fNIRS) is a promising tool for real-time applications such as neurofeedback and brain-computer interfaces. Its combination of spatial specificity and mobility makes it particularly attractive for clinical use, both at the bedside and in patients' homes. Despite these advantages, optimizing fNIRS for real-time use requires careful attention to two key aspects: ensuring good spatial specificity and maintaining high signal quality. While fNIRS detects superficial cortical brain regions, consistently and reliably targeting specific regions of interest can be challenging, particularly in studies that require repeated measurements. Variations in cap placement coupled with limited anatomical information may further reduce this accuracy. Furthermore, it is important to maintain good signal quality in real-time contexts to ensure that they reflect the true underlying brain activity. However, fNIRS signals are susceptible to contamination by cerebral and extracerebral systemic noise as well as motion artifacts. Insufficient real-time preprocessing can therefore cause the system to run on noise instead of brain activity. The aim of this review article is to help advance the progress of fNIRS-based real-time applications. It highlights the potential challenges in improving spatial specificity and signal quality, discusses possible options to overcome these challenges, and addresses further considerations relevant to real-time applications. By addressing these topics, the article aims to help improve the planning and execution of future real-time studies, thereby increasing their reliability and repeatability.

dCCA: detecting differential covariation patterns between two types of high-throughput omics data

MOTIVATION

The advent of multimodal omics data has provided an unprecedented opportunity to systematically investigate underlying biological mechanisms from distinct yet complementary angles. However, the joint analysis of multi-omics data remains challenging because it requires modeling interactions between multiple sets of high-throughput variables. Furthermore, these interaction patterns may vary across different clinical groups, reflecting disease-related biological processes.

RESULTS

We propose a novel approach called Differential Canonical Correlation Analysis (dCCA) to capture differential covariation patterns between two multivariate vectors across clinical groups. Unlike classical Canonical Correlation Analysis, which maximizes the correlation between two multivariate vectors, dCCA aims to maximally recover differentially expressed multivariate-to-multivariate covariation patterns between groups. We have developed computational algorithms and a toolkit to sparsely select paired subsets of variables from two sets of multivariate variables while maximizing the differential covariation. Extensive simulation analyses demonstrate the superior performance of dCCA in selecting variables of interest and recovering differential correlations. We applied dCCA to the Pan-Kidney cohort from the Cancer Genome Atlas Program database and identified differentially expressed covariations between noncoding RNAs and gene expressions.

AVAILABILITY AND IMPLEMENTATION

The R package that implements dCCA is available at https://github.com/hwiyoungstat/dCCA.

Neurobiology-based Cognitive Biotypes Using Multi-scale Intrinsic Connectivity Networks in Psychotic Disorders

Objective

Understanding the neurobiology of cognitive dysfunction in psychotic disorders remains elusive, as does developing effective interventions. Limited knowledge about the biological heterogeneity of cognitive dysfunction hinders progress. This study aimed to identify subgroups of patients with psychosis with distinct patterns of functional brain alterations related to cognition (cognitive biotypes).

Methods

B-SNIP consortium data (2,270 participants including participants with psychotic disorders, relatives, and controls) was analyzed. Researchers used reference-informed independent component analysis and the NeuroMark 100k multi-scale intrinsic connectivity networks (ICN) template to obtain subject-specific ICNs and whole-brain functional network connectivity (FNC). FNC features associated with cognitive performance were identified through multivariate joint analysis. K-means clustering identified subgroups of patients based on these features in a discovery set. Subgroups were further evaluated in a replication set and in relatives.

Results

Two biotypes with different functional brain alteration patterns were identified. Biotype 1 exhibited brain-wide alterations, involving hypoconnectivity in cerebellar-subcortical and somatomotor-visual networks and worse cognitive performance. Biotype 2 exhibited hyperconnectivity in somatomotor-subcortical networks and hypoconnectivity in somatomotor-high cognitive processing networks, and better preserved cognitive performance. Demographic, clinical, cognitive, and FNC characteristics of biotypes were consistent in discovery and replication sets, and in relatives. 70.12% of relatives belonged to the same biotype as their affected family members.

Conclusions

These findings suggest two distinctive psychosis-related cognitive biotypes with differing functional brain patterns shared with their relatives. Patient stratification based on these biotypes instead of traditional diagnosis may help to optimize future research and clinical trials addressing cognitive dysfunction in psychotic disorders.

Leveraging multivariate approaches to advance the science of early-life adversity

Since the landmark Adverse Childhood Experiences (ACEs) study, adversity research has expanded to more precisely account for the multifaceted nature of adverse experiences. The complex data structures and interrelated nature of adversity data require robust multivariate statistical methods, and recent methodological and statistical innovations have facilitated advancements in research on childhood adversity. Here, we provide an overview of a subset of multivariate methods that we believe hold particular promise for advancing the field's understanding of early-life adversity, and discuss how these approaches can be practically applied to explore different research questions. This review covers data-driven or unsupervised approaches (including dimensionality reduction and person-centered clustering/subtype identification) as well as supervised/prediction-based approaches (including linear and tree-based models and neural networks). For each, we highlight studies that have effectively applied the method to provide novel insight into early-life adversity. Taken together, we hope this review serves as a resource to adversity researchers looking to expand upon the cumulative approach described in the original ACEs study, thereby advancing the field's understanding of the complexity of adversity and related developmental consequences.

On the stability of canonical correlation analysis and partial least squares with application to brain-behavior associations

Associations between datasets can be discovered through multivariate methods like Canonical Correlation Analysis (CCA) or Partial Least Squares (PLS). A requisite property for interpretability and generalizability of CCA/PLS associations is stability of their feature patterns. However, stability of CCA/PLS in high-dimensional datasets is questionable, as found in empirical characterizations. To study these issues systematically, we developed a generative modeling framework to simulate synthetic datasets. We found that when sample size is relatively small, but comparable to typical studies, CCA/PLS associations are highly unstable and inaccurate; both in their magnitude and importantly in the feature pattern underlying the association. We confirmed these trends across two neuroimaging modalities and in independent datasets with n ≈ 1000 and n = 20,000, and found that only the latter comprised sufficient observations for stable mappings between imaging-derived and behavioral features. We further developed a power calculator to provide sample sizes required for stability and reliability of multivariate analyses. Collectively, we characterize how to limit detrimental effects of overfitting on CCA/PLS stability, and provide recommendations for future studies.

The links between physical activity and prosocial behavior: an fNIRS hyperscanning study

The prevalence of physically inactive lifestyles in modern society raises concerns about the potential association with poor brain health, particularly in the lateral prefrontal cortex, which is crucial for human prosocial behavior. Here, we explored the relationship between physical activity and prosocial behavior, focusing on potential neural markers, including intra-brain functional connectivity and inter-brain synchrony in the lateral prefrontal cortex. Forty participants, each paired with a stranger, completed two experimental conditions in a randomized order: (i) face-to-face and (ii) face stimulus (eye-to-eye contact with a face stimulus of a fictitious person displayed on the screen). Following each condition, participants played economic games with either their partner or an assumed person displayed on the screen. Neural activity in the lateral prefrontal cortex was recorded by functional near-infrared spectroscopy hyperscanning. Sparse multiset canonical correlation analysis showed that a physically inactive lifestyle was covaried with poorer reciprocity, greater trust, shorter decision-making time, and weaker intra-brain connectivity in the dorsal lateral prefrontal cortex and poorer inter-brain synchrony in the ventral lateral prefrontal cortex. These associations were observed exclusively in the face-to-face condition. Our findings suggest that a physically inactive lifestyle may alter human prosocial behavior by impairing adaptable prosocial decision-making in response to social factors through altered intra-brain functional connectivity and inter-brain synchrony.

Reduced covariation between brain morphometry and local spontaneous activity in young children with ASD

While disruptions in brain maturation in the first years of life in ASD are well documented, little is known about how the brain structure and function are related in young children with ASD compared to typically developing peers. We applied a multivariate pattern analysis to examine the covariation patterns between brain morphometry and local brain spontaneous activity in 38 toddlers and preschoolers with ASD and 31 typically developing children using T1-weighted structural MRI and resting-state fMRI data acquired during natural sleep. The results revealed significantly reduced brain structure-function correlations in ASD. The resultant brain structure and function composite indices were associated with age among typically developing children, but not among those with ASD, suggesting mistiming of typical brain maturational trajectories early in life in autism. Additionally, the brain function composite indices were associated with the overall developmental and adaptive behavior skills in the ASD group, highlighting the neurodevelopmental significance of early local brain activity in autism.

Unsupervised Neural Manifold Alignment for Stable Decoding of Movement from Cortical Signals

The stable decoding of movement parameters using neural activity is crucial for the success of brain-machine interfaces (BMIs). However, neural activity can be unstable over time, leading to changes in the parameters used for decoding movement, which can hinder accurate movement decoding. To tackle this issue, one approach is to transfer neural activity to a stable, low-dimensional manifold using dimensionality reduction techniques and align manifolds across sessions by maximizing correlations of the manifolds. However, the practical use of manifold stabilization techniques requires knowledge of the true subject intentions such as target direction or behavioral state. To overcome this limitation, an automatic unsupervised algorithm is proposed that determines movement target intention before manifold alignment in the presence of manifold rotation and scaling across sessions. This unsupervised algorithm is combined with a dimensionality reduction and alignment method to overcome decoder instabilities. The effectiveness of the BMI stabilizer method is represented by decoding the two-dimensional (2D) hand velocity of two rhesus macaque monkeys during a center-out-reaching movement task. The performance of the proposed method is evaluated using correlation coefficient and R-squared measures, demonstrating higher decoding performance compared to a state-of-the-art unsupervised BMI stabilizer. The results offer benefits for the automatic determination of movement intents in long-term BMI decoding. Overall, the proposed method offers a promising automatic solution for achieving stable and accurate movement decoding in BMI applications.

Static and dynamic fMRI-derived functional connectomes represent largely similar information

Functional connectivity (FC) of blood oxygen level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate), and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Fair Canonical Correlation Analysis

This paper investigates fairness and bias in Canonical Correlation Analysis (CCA), a widely used statistical technique for examining the relationship between two sets of variables. We present a framework that alleviates unfairness by minimizing the correlation disparity error associated with protected attributes. Our approach enables CCA to learn global projection matrices from all data points while ensuring that these matrices yield comparable correlation levels to group-specific projection matrices. Experimental evaluation on both synthetic and real-world datasets demonstrates the efficacy of our method in reducing correlation disparity error without compromising CCA accuracy.

Multivariate Associations Among White Matter, Neurocognition, and Social Cognition Across Individuals With Schizophrenia Spectrum Disorders and Healthy Controls

BACKGROUND AND HYPOTHESIS

Neurocognitive and social cognitive abilities are important contributors to functional outcomes in schizophrenia spectrum disorders (SSDs). An unanswered question of considerable interest is whether neurocognitive and social cognitive deficits arise from overlapping or distinct white matter impairment(s).

STUDY DESIGN

We sought to fill this gap, by harnessing a large sample of individuals from the multi-center Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) dataset, unique in its collection of advanced diffusion imaging and an extensive battery of cognitive assessments. We applied canonical correlation analysis to estimates of white matter microstructure, and cognitive performance, across people with and without an SSD.

STUDY RESULTS

Our results established that white matter circuitry is dimensionally and strongly related to both neurocognition and social cognition, and that microstructure of the uncinate fasciculus and the rostral body of the corpus callosum may assume a "privileged role" subserving both. Further, we found that participant-wise estimates of white matter microstructure, weighted by cognitive performance, were largely consistent with participants' categorical diagnosis, and predictive of (cross-sectional) functional outcomes.

CONCLUSIONS

The demonstrated strength of the relationship between white matter circuitry and neurocognition and social cognition underscores the potential for using relationships among these variables to identify biomarkers of functioning, with potential prognostic and therapeutic implications.

Assessing the impacts of climate variables on long-term air quality trends in Peninsular Malaysia

Climate change is thought to influence the composition of atmospheric air, but little is known about the direct relationship between these variables, especially in a hot tropical climate like that of Malaysia. This work summarizes and analyzes the climate state and air quality of Peninsular Malaysia based on selected ground-based observations of the temperature, precipitation, relative humidity, wind speed, wind direction and concentrations of PM10, O3, CO, NO2, and SO2 over the last 20 years (2000-2019). The relationship between the climate state and air quality is analyzed using the Pearson correlation and canonical correlation analysis (CCA) methods is employed to predict the degree of change in the future air quality under different warming scenarios. It is found that the Peninsular Malaysia mainly experienced strong precipitation in the central and mountainous regions, while air pollutants are primarily concentrated in densely populated areas. Throughout the period of study (interannual, monthly, and diurnal time series analyses), Peninsular Malaysia became warmer and drier, with a significant increase in temperature (+4.2 %), decrease in the relative humidity (-4.5 %), and greater fluctuation in precipitation amount. The pollution conditions have worsened; there has been an increase in the PM10 (+16.4 %), O3 (+39.5 %), and NO2 (+2.1 %) concentration over the last 20 years. However, the amount of SO2 (-53.6 %) and CO (-20.6 %) decreased significantly. The analysis of the monthly variation shows a strong bimodality of the PM10 and O3 concentrations that corresponds to the monsoon transition. Intensive diurnal fluctuations and correlations are observed for all the variables in this study. According to the CCA, the air quality factors are strongly correlated with meteorological factors; in particular, the CO, O3, and PM10 concentrations interact strongly with the air temperature. These findings show that the future air quality in Peninsular Malaysia has high possibility to deteriorate under warming condition.

Multi-view information fusion using multi-view variational autoencoder to predict proximal femoral fracture load

Background

Hip fracture occurs when an applied force exceeds the force that the proximal femur can support (the fracture load or "strength") and can have devastating consequences with poor functional outcomes. Proximal femoral strengths for specific loading conditions can be computed by subject-specific finite element analysis (FEA) using quantitative computerized tomography (QCT) images. However, the radiation and availability of QCT limit its clinical usability. Alternative low-dose and widely available measurements, such as dual energy X-ray absorptiometry (DXA) and genetic factors, would be preferable for bone strength assessment. The aim of this paper is to design a deep learning-based model to predict proximal femoral strength using multi-view information fusion.

Results

We developed new models using multi-view variational autoencoder (MVAE) for feature representation learning and a product of expert (PoE) model for multi-view information fusion. We applied the proposed models to an in-house Louisiana Osteoporosis Study (LOS) cohort with 931 male subjects, including 345 African Americans and 586 Caucasians. We performed genome-wide association studies (GWAS) to select 256 genetic variants with the lowest p-values for each proximal femoral strength and integrated whole genome sequence (WGS) features and DXA-derived imaging features to predict proximal femoral strength. The best prediction model for fall fracture load was acquired by integrating WGS features and DXA-derived imaging features. The designed models achieved the mean absolute percentage error of 18.04%, 6.84% and 7.95% for predicting proximal femoral fracture loads using linear models of fall loading, nonlinear models of fall loading, and nonlinear models of stance loading, respectively.

Conclusion

The proposed models are capable of predicting proximal femoral strength using WGS features and DXA-derived imaging features. Though this tool is not a substitute for predicting FEA using QCT images, it would make improved assessment of hip fracture risk more widely available while avoiding the increased radiation exposure from QCT.

Characterising the covariance pattern between lifestyle factors and structural brain measures: a multivariable replication study of two independent ageing cohorts

Modifiable lifestyle factors have been shown to promote healthy brain ageing. However, studies have typically focused on a single factor at a time. Given that lifestyle factors do not occur in isolation, multivariable analyses provide a more realistic model of the lifestyle-brain relationship. Here, canonical correlation analyses (CCA) examined the relationship between nine lifestyle factors and seven MRI-derived indices of brain structure. The resulting covariance pattern was further explored with Bayesian regressions. CCA analyses were first conducted on a Danish cohort of older adults (n = 251) and then replicated in a British cohort (n = 668). In both cohorts, the latent factors of lifestyle and brain structure were positively correlated (UK: r = .37, p < 0.001; Denmark: r = .27, p < 0.001). In the cross-validation study, the correlation between lifestyle-brain latent factors was r = .10, p = 0.008. However, the pattern of associations differed between datasets. These findings suggest that baseline characterisation and tailoring towards the study sample may be beneficial for achieving targeted lifestyle interventions.

Role of caveolin-1 in metabolic programming of fetal brain

Mice lacking caveolin-1 (Cav1), a key protein of plasma membrane, exhibit brain aging at an early adult stage. Here, integrative analyses of metabolomics, transcriptomics, epigenetics, and single-cell data were performed to test the hypothesis that metabolic deregulation of fetal brain due to the ablation of Cav1 is linked to brain aging in these mice. The results of this study show that lack of Cav1 caused deregulation in the lipid and amino acid metabolism in the fetal brain, and genes associated with these deregulated metabolites were significantly altered in the brain upon aging. Moreover, ablation of Cav1 deregulated several metabolic genes in specific cell types of the fetal brain and impacted DNA methylation of those genes in coordination with mouse epigenetic clock. The findings of this study suggest that the aging program of brain is confounded by metabolic abnormalities in the fetal stage due to the absence of Cav1.

Evaluation of the Hierarchical Correspondence between the Human Brain and Artificial Neural Networks: A Review

Artificial neural networks (ANNs) that are heavily inspired by the human brain now achieve human-level performance across multiple task domains. ANNs have thus drawn attention in neuroscience, raising the possibility of providing a framework for understanding the information encoded in the human brain. However, the correspondence between ANNs and the brain cannot be measured directly. They differ in outputs and substrates, neurons vastly outnumber their ANN analogs (i.e., nodes), and the key algorithm responsible for most of modern ANN training (i.e., backpropagation) is likely absent from the brain. Neuroscientists have thus taken a variety of approaches to examine the similarity between the brain and ANNs at multiple levels of their information hierarchy. This review provides an overview of the currently available approaches and their limitations for evaluating brain-ANN correspondence.

Spatial association between distributed β-amyloid and tau varies with cognition

Several PET studies have explored the relationship between β-amyloid load and tau uptake at the early stages of Alzheimer's disease (AD) progression. Most of these studies have focused on the linear relationship between β-amyloid and tau at the local level and their synergistic effect on different AD biomarkers. We hypothesize that patterns of spatial association between β-amyloid and tau might be uncovered using alternative association metrics that account for linear as well as more complex, possible nonlinear dependencies. In the present study, we propose a new Canonical Distance Correlation Analysis (CDCA) to generate distinctive spatial patterns of the cross-correlation structure between tau, as measured by [18F]flortaucipir PET, and β-amyloid, as measured by [18F]florbetapir PET, from the Alzheimer's Disease Neuroimaging Initiative (ADNI) study. We found that the CDCA-based β-amyloid scores were not only maximally distance-correlated to tau in cognitively normal (CN) controls and mild cognitive impairment (MCI), but also differentiated between low and high levels of β-amyloid uptake. The most distinctive spatial association pattern was characterized by a spread of β-amyloid covering large areas of the cortex and localized tau in the entorhinal cortex. More importantly, this spatial dependency varies according to cognition, which cannot be explained by the uptake differences in β-amyloid or tau between CN and MCI subjects. Hence, the CDCA-based scores might be more accurate than the amyloid or tau SUVR for the enrollment in clinical trials of those individuals on the path of cognitive deterioration.

Deep learning integrates histopathology and proteogenomics at a pan-cancer level

We introduce a pioneering approach that integrates pathology imaging with transcriptomics and proteomics to identify predictive histology features associated with critical clinical outcomes in cancer. We utilize 2,755 H&E-stained histopathological slides from 657 patients across 6 cancer types from CPTAC. Our models effectively recapitulate distinctions readily made by human pathologists: tumor vs. normal (AUROC = 0.995) and tissue-of-origin (AUROC = 0.979). We further investigate predictive power on tasks not normally performed from H&E alone, including TP53 prediction and pathologic stage. Importantly, we describe predictive morphologies not previously utilized in a clinical setting. The incorporation of transcriptomics and proteomics identifies pathway-level signatures and cellular processes driving predictive histology features. Model generalizability and interpretability is confirmed using TCGA. We propose a classification system for these tasks, and suggest potential clinical applications for this integrated human and machine learning approach. A publicly available web-based platform implements these models.

Homogeneous-Multiset-CCA-Based Brain Covariation and Contravariance Connectivity Network Modeling

Brain connectivity networks based on functional magnetic resonance imaging (fMRI) have expanded our understanding of brain functions in both healthy and diseased states. However, most current studies construct connectivity networks using averaged regional time courses with the strong assumption that the activities of voxels contained in each brain region are similar, ignoring their possible variations. Additionally, pairwise correlation analysis is often adopted with more attention to positive relationships, while joint interactions at the network level as well as anti-correlations are less investigated. In this paper, to provide a new strategy for regional activity representation and brain connectivity modeling, a novel homogeneous multiset canonical correlation analysis (HMCCA) model is proposed, which enforces sign constraints on the weights of voxels to guarantee homogeneity within each brain region. It is capable of obtaining regional representative signals and constructing covariation and contravariance networks simultaneously, at both group and subject levels. Validations on two sessions of fMRI data verified its reproducibility and reliability when dealing with brain connectivity networks. Further experiments on subjects with and without Parkinson's disease (PD) revealed significant alterations in brain connectivity patterns, which were further associated with clinical scores and demonstrated superior prediction ability, indicating its potential in clinical practice.

Evaluation of Classification Systems and Their Correlation With Clinical and Quality-of-life Parameters in Patients With Surgically Treated Spondylodiskitis

INTRODUCTION

Surgical classifications for spondylodiscitis (SD) typically include radiologic features and the status of neurologic impairment. Clinical factors such as preoperative pain, function/disability, overall quality of life (QoL), and risk of recurrence and mortality, which are necessary for a comprehensive assessment of SD and measurement of treatment success, are not considered. There is a lack of external validation of SD classifications. The aim of this study was to validate classifications of SD and to correlate these classifications with the above clinical factors.

METHODS

One hundred fourteen patients from a prospective SD register (2008-2020) with available imaging, preoperative neurologic status, backpain, function/disability data (Oswestry Disability Index and Core Outcome Measures Index), QoL data (Short Form 36, European Quality-of-life Questionnaire), and a 1-year follow-up were retrospectively classified according to Akbar, Homagk, and Pola classifications. Interrater reliability, correlation among classifications, and correlation between classifications and QoL were calculated.

RESULTS

Interrater reliability was κ = 0.83 for Akbar, κ = 0.94 for Homagk, and κ = 0.99 for Pola. The correlation of Akbar with Pola and Homagk was moderate (ρ s = 0.47; ρ s = 0.46) and high between Pola and Homagk (ρ s = 0.7). No notable correlation was observed between any of the classifications and preoperative Oswestry Disability Index, Core Outcome Measures Index, QoL, mortality, and recurrence within 1 year. Only a weak correlation was observed between Homagk and preoperative leg pain and back pain.

CONCLUSION

Available SD classifications have a very good interrater reliability and moderate-to-high correlation with each other but lack correlation with preoperative pain, function/disability, and overall QoL. Because these factors are important for a comprehensive assessment of SD in severity, decision making, and prognosis, they should be included in future SD classifications. This could allow for more comprehensive treatment algorithms.

LEVEL OF EVIDENCE

Level II. Diagnostic study = prospective cohort study; development of diagnostic criteria.

DATA AVAILABILITY

The data sets used and analyzed during this study are available from the corresponding author on reasonable request.

Multivariate analytical approaches for investigating brain-behavior relationships

Background

Many studies of brain-behavior relationships rely on univariate approaches where each variable of interest is tested independently, which does not allow for the simultaneous investigation of multiple correlated variables. Alternatively, multivariate approaches allow for examining relationships between psychopathology and neural substrates simultaneously. There are multiple multivariate methods to choose from that each have assumptions which can affect the results; however, many studies employ one method without a clear justification for its selection. Additionally, there are few studies illustrating how differences between methods manifest in examining brain-behavior relationships. The purpose of this study was to exemplify how the choice of multivariate approach can change brain-behavior interpretations.

Method

We used data from 9,027 9- to 10-year-old children from the Adolescent Brain Cognitive DevelopmentSM Study (ABCD Study®) to examine brain-behavior relationships with three commonly used multivariate approaches: canonical correlation analysis (CCA), partial least squares correlation (PLSC), and partial least squares regression (PLSR). We examined the associations between psychopathology dimensions including general psychopathology, attention-deficit/hyperactivity symptoms, conduct problems, and internalizing symptoms with regional brain volumes.

Results

The results of CCA, PLSC, and PLSR showed both consistencies and differences in the relationship between psychopathology symptoms and brain structure. The leading significant component yielded by each method demonstrated similar patterns of associations between regional brain volumes and psychopathology symptoms. However, the additional significant components yielded by each method demonstrated differential brain-behavior patterns that were not consistent across methods.

Conclusion

Here we show that CCA, PLSC, and PLSR yield slightly different interpretations regarding the relationship between child psychopathology and brain volume. In demonstrating the divergence between these approaches, we exemplify the importance of carefully considering the method's underlying assumptions when choosing a multivariate approach to delineate brain-behavior relationships.

Examining Factors Influencing Cognitive Load of Computer Programmers

In this study, the factors influencing the cognitive load of computer programmers during the perception of different code tasks were investigated. The eye movement features of computer programmers were used to provide a significant relationship between the perceptual processes of the sample codes and cognitive load. Thanks to the relationship, the influence of various personal characteristics of programmers on cognitive load was examined. Various personal parameters such as programming experience, age, native language, and programming frequency were used in the study. The study was performed on the Eye Movements in Programming (EMIP) dataset containing 216 programmers with different characteristics. Eye movement information recorded during two different code comprehension tasks was decomposed into sub-information, such as pupil movement speed and diameter change. Rapid changes in eye movement signals were adaptively detected using the z-score peak detection algorithm. Regarding the cognitive load calculations, canonical correlation analysis was used to build a statistically significant and efficient mathematical model connecting the extracted eye movement features and the different parameters of the programmers, and the results were statistically significant. As a result of the analysis, the factors affecting the cognitive load of computer programmers for the related database were converted into percentages, and it was seen that linguistic distance is an essential factor in the cognitive load of programmers and the effect of gender on cognitive load is quite limited.

Penalized decomposition using residuals (PeDecURe) for feature extraction in the presence of nuisance variables

Neuroimaging data are an increasingly important part of etiological studies of neurological and psychiatric disorders. However, mitigating the influence of nuisance variables, including confounders, remains a challenge in image analysis. In studies of Alzheimer's disease, for example, an imbalance in disease rates by age and sex may make it difficult to distinguish between structural patterns in the brain (as measured by neuroimaging scans) attributable to disease progression and those characteristic of typical human aging or sex differences. Concerningly, when not properly accounted for, nuisance variables pose threats to the generalizability and interpretability of findings from these studies. Motivated by this critical issue, in this work, we examine the impact of nuisance variables on feature extraction methods and propose Penalized Decomposition Using Residuals (PeDecURe), a new method for obtaining nuisance variable-adjusted features. PeDecURe estimates primary directions of variation which maximize covariance between partially residualized imaging features and a variable of interest (e.g., Alzheimer's diagnosis) while simultaneously mitigating the influence of nuisance variation through a penalty on the covariance between partially residualized imaging features and those variables. Using features derived using PeDecURe's first direction of variation, we train a highly accurate and generalizable predictive model, as evidenced by its robustness in testing samples with different underlying nuisance variable distributions. We compare PeDecURe to commonly used decomposition methods (principal component analysis (PCA) and partial least squares) as well as a confounder-adjusted variation of PCA. We find that features derived from PeDecURe offer greater accuracy and generalizability and lower correlations with nuisance variables compared with the other methods. While PeDecURe is primarily motivated by challenges that arise in the analysis of neuroimaging data, it is broadly applicable to data sets with highly correlated features, where novel methods to handle nuisance variables are warranted.

Multi-View Learning to Unravel the Different Levels Underlying Hepatitis B Vaccine Response

The immune system acts as an intricate apparatus that is dedicated to mounting a defense and ensures host survival from microbial threats. To engage this faceted immune response and provide protection against infectious diseases, vaccinations are a critical tool to be developed. However, vaccine responses are governed by levels that, when interrogated, separately only explain a fraction of the immune reaction. To address this knowledge gap, we conducted a feasibility study to determine if multi-view modeling could aid in gaining actionable insights on response markers shared across populations, capture the immune system's diversity, and disentangle confounders. We thus sought to assess this multi-view modeling capacity on the responsiveness to the Hepatitis B virus (HBV) vaccination. Seroconversion to vaccine-induced antibodies against the HBV surface antigen (anti-HBs) in early converters (n = 21; <2 months) and late converters (n = 9; <6 months) and was defined based on the anti-HBs titers (>10IU/L). The multi-view data encompassed bulk RNA-seq, CD4+ T-cell parameters (including T-cell receptor data), flow cytometry data, and clinical metadata (including age and gender). The modeling included testing single-view and multi-view joint dimensionality reductions. Multi-view joint dimensionality reduction outperformed single-view methods in terms of the area under the curve and balanced accuracy, confirming the increase in predictive power to be gained. The interpretation of these findings showed that age, gender, inflammation-related gene sets, and pre-existing vaccine-specific T-cells could be associated with vaccination responsiveness. This multi-view dimensionality reduction approach complements clinical seroconversion and all single modalities. Importantly, this modeling could identify what features could predict HBV vaccine response. This methodology could be extended to other vaccination trials to identify the key features regulating responsiveness.

Canonical correlation analysis in high dimensions with structured regularization

Canonical correlation analysis (CCA) is a technique for measuring the association between two multivariate data matrices. A regularized modification of canonical correlation analysis (RCCA) which imposes an ℓ2 penalty on the CCA coefficients is widely used in applications with high-dimensional data. One limitation of such regularization is that it ignores any data structure, treating all the features equally, which can be ill-suited for some applications. In this article we introduce several approaches to regularizing CCA that take the underlying data structure into account. In particular, the proposed group regularized canonical correlation analysis (GRCCA) is useful when the variables are correlated in groups. We illustrate some computational strategies to avoid excessive computations with regularized CCA in high dimensions. We demonstrate the application of these methods in our motivating application from neuroscience, as well as in a small simulation example.

Static and dynamic functional connectomes represent largely similar information

Functional connectivity (FC) of blood-oxygen-level-dependent (BOLD) fMRI time series can be estimated using methods that differ in sensitivity to the temporal order of time points (static vs. dynamic) and the number of regions considered in estimating a single edge (bivariate vs. multivariate). Previous research suggests that dynamic FC explains variability in FC fluctuations and behavior beyond static FC. Our aim was to systematically compare methods on both dimensions. We compared five FC methods: Pearson's/full correlation (static, bivariate), lagged correlation (dynamic, bivariate), partial correlation (static, multivariate) and multivariate AR model with and without self-connections (dynamic, multivariate). We compared these methods by (i) assessing similarities between FC matrices, (ii) by comparing node centrality measures, and (iii) by comparing the patterns of brain-behavior associations. Although FC estimates did not differ as a function of sensitivity to temporal order, we observed differences between the multivariate and bivariate FC methods. The dynamic FC estimates were highly correlated with the static FC estimates, especially when comparing group-level FC matrices. Similarly, there were high correlations between the patterns of brain-behavior associations obtained using the dynamic and static FC methods. We conclude that the dynamic FC estimates represent information largely similar to that of the static FC.

Dysconnection and cognition in schizophrenia: A spectral dynamic causal modeling study

Schizophrenia (SZ) is a severe mental disorder characterized by failure of functional integration (aka dysconnection) across the brain. Recent functional connectivity (FC) studies have adopted functional parcellations to define subnetworks of large-scale networks, and to characterize the (dys)connection between them, in normal and clinical populations. While FC examines statistical dependencies between observations, model-based effective connectivity (EC) can disclose the causal influences that underwrite the observed dependencies. In this study, we investigated resting state EC within seven large-scale networks, in 66 SZ and 74 healthy subjects from a public dataset. The results showed that a remarkable 33% of the effective connections (among subnetworks) of the cognitive control network had been pathologically modulated in SZ. Further dysconnection was identified within the visual, default mode and sensorimotor networks of SZ subjects, with 24%, 20%, and 11% aberrant couplings. Overall, the proportion of discriminative connections was remarkably larger in EC (24%) than FC (1%) analysis. Subsequently, to study the neural correlates of impaired cognition in SZ, we conducted a canonical correlation analysis between the EC parameters and the cognitive scores of the patients. As such, the self-inhibitions of supplementary motor area and paracentral lobule (in the sensorimotor network) and the excitatory connection from parahippocampal gyrus to inferior temporal gyrus (in the cognitive control network) were significantly correlated with the social cognition, reasoning/problem solving and working memory capabilities of the patients. Future research can investigate the potential of whole-brain EC as a biomarker for diagnosis of brain disorders and for neuroimaging-based cognitive assessment.

Identifying vulnerable brain networks associated with Alzheimer's disease risk

The selective vulnerability of brain networks in individuals at risk for Alzheimer's disease (AD) may help differentiate pathological from normal aging at asymptomatic stages, allowing the implementation of more effective interventions. We used a sample of 72 people across the age span, enriched for the APOE4 genotype to reveal vulnerable networks associated with a composite AD risk factor including age, genotype, and sex. Sparse canonical correlation analysis (CCA) revealed a high weight associated with genotype, and subgraphs involving the cuneus, temporal, cingulate cortices, and cerebellum. Adding cognitive metrics to the risk factor revealed the highest cumulative degree of connectivity for the pericalcarine cortex, insula, banks of the superior sulcus, and the cerebellum. To enable scaling up our approach, we extended tensor network principal component analysis, introducing CCA components. We developed sparse regression predictive models with errors of 17% for genotype, 24% for family risk factor for AD, and 5 years for age. Age prediction in groups including cognitively impaired subjects revealed regions not found using only normal subjects, i.e. middle and transverse temporal, paracentral and superior banks of temporal sulcus, as well as the amygdala and parahippocampal gyrus. These modeling approaches represent stepping stones towards single subject prediction.

Targeted rhythmic visual stimulation at individual participants' intrinsic alpha frequency causes selective increase of occipitoparietal BOLD-fMRI and EEG functional connectivity

Neural oscillations in distinct frequency bands are ubiquitous in the brain and play a role in many cognitive processes. The "communication by coherence" hypothesis, poses that the synchronization through phase coupling of frequency-specific neural oscillations regulate information flow across distribute brain regions. Specifically, the posterior alpha frequency band (7-12 Hz) is thought to gate bottom-up visual information flow by inhibition during visual processing. Evidence shows that increased alpha phase coherency positively correlates with functional connectivity in resting state connectivity networks, supporting alpha mediates neural communication through coherency. However, these findings have mainly been derived from spontaneous changes in the ongoing alpha rhythm. In this study, we experimentally modulate the alpha rhythm by targeting individuals' intrinsic alpha frequency with sustained rhythmic light to investigate alpha-mediated synchronous cortical activity in both EEG and fMRI. We hypothesize increased alpha coherency and fMRI connectivity should arise from modulation of the intrinsic alpha frequency (IAF) as opposed to control frequencies in the alpha range. Sustained rhythmic and arrhythmic stimulation at the IAF and at neighboring frequencies within the alpha band range (7-12 Hz) was implemented and assessed in a separate EEG and fMRI study. We observed increased cortical alpha phase coherency in the visual cortex during rhythmic stimulation at the IAF as in comparison to rhythmic stimulation of control frequencies. In the fMRI, we found increased functional connectivity for stimulation at the IAF in visual and parietal areas as compared to other rhythmic control frequencies by correlating time courses from a set of regions of interest for the different stimulation conditions and applying network-based statistics. This suggests that rhythmic stimulation at the IAF frequency induces a higher degree of synchronicity of neural activity across the occipital and parietal cortex, which supports the role of the alpha oscillation in gating information flow during visual processing.

Prenatal exposure to air pollution is associated with structural changes in the neonatal brain

BACKGROUND

Prenatal exposure to air pollution is associated with adverse neurologic consequences in childhood. However, the relationship between in utero exposure to air pollution and neonatal brain development is unclear.

METHODS

We modelled maternal exposure to nitrogen dioxide (NO2) and particulate matter (PM2.5 and PM10) at postcode level between date of conception to date of birth and studied the effect of prenatal air pollution exposure on neonatal brain morphology in 469 (207 male) healthy neonates, with gestational age of ≥36 weeks. Infants underwent MR neuroimaging at 3 Tesla at 41.29 (36.71-45.14) weeks post-menstrual age (PMA) as part of the developing human connectome project (dHCP). Single pollutant linear regression and canonical correlation analysis (CCA) were performed to assess the relationship between air pollution and brain morphology, adjusting for confounders and correcting for false discovery rate.

RESULTS

Higher exposure to PM10 and lower exposure to NO2 was strongly canonically correlated to a larger relative ventricular volume, and moderately associated with larger relative size of the cerebellum. Modest associations were detected with higher exposure to PM10 and lower exposure to NO2 and smaller relative cortical grey matter and amygdala and hippocampus, and larger relaive brainstem and extracerebral CSF volume. No associations were found with white matter or deep grey nuclei volume.

CONCLUSIONS

Our findings show that prenatal exposure to air pollution is associated with altered brain morphometry in the neonatal period, albeit with opposing results for NO2 and PM10. This finding provides further evidence that reducing levels of maternal exposure to particulate matter during pregnancy should be a public health priority and highlights the importance of understanding the impacts of air pollution on this critical development window.

EEGformer: A transformer-based brain activity classification method using EEG signal

Background

The effective analysis methods for steady-state visual evoked potential (SSVEP) signals are critical in supporting an early diagnosis of glaucoma. Most efforts focused on adopting existing techniques to the SSVEPs-based brain-computer interface (BCI) task rather than proposing new ones specifically suited to the domain.

Method

Given that electroencephalogram (EEG) signals possess temporal, regional, and synchronous characteristics of brain activity, we proposed a transformer-based EEG analysis model known as EEGformer to capture the EEG characteristics in a unified manner. We adopted a one-dimensional convolution neural network (1DCNN) to automatically extract EEG-channel-wise features. The output was fed into the EEGformer, which is sequentially constructed using three components: regional, synchronous, and temporal transformers. In addition to using a large benchmark database (BETA) toward SSVEP-BCI application to validate model performance, we compared the EEGformer to current state-of-the-art deep learning models using two EEG datasets, which are obtained from our previous study: SJTU emotion EEG dataset (SEED) and a depressive EEG database (DepEEG).

Results

The experimental results show that the EEGformer achieves the best classification performance across the three EEG datasets, indicating that the rationality of our model architecture and learning EEG characteristics in a unified manner can improve model classification performance.

Conclusion

EEGformer generalizes well to different EEG datasets, demonstrating our approach can be potentially suitable for providing accurate brain activity classification and being used in different application scenarios, such as SSVEP-based early glaucoma diagnosis, emotion recognition and depression discrimination.

Multivariate brain-based dimensions of child psychiatric problems: degrees of generalizability

Multivariate machine learning techniques are a promising set of tools for identifying complex brain-behavior associations. However, failure to replicate results from these methods across samples has hampered their clinical relevance. This study aimed to delineate dimensions of brain functional connectivity that are associated with child psychiatric symptoms in two large and independent cohorts: the Adolescent Brain Cognitive Development (ABCD) Study and the Generation R Study (total n =8,605). Using sparse canonical correlations analysis, we identified three brain-behavior dimensions in ABCD: attention problems, aggression and rule-breaking behaviors, and withdrawn behaviors. Importantly, out-of-sample generalizability of these dimensions was consistently observed in ABCD, suggesting robust multivariate brain-behavior associations. Despite this, out-of-study generalizability in Generation R was limited. These results highlight that the degree of generalizability can vary depending on the external validation methods employed as well as the datasets used, emphasizing that biomarkers will remain elusive until models generalize better in true external settings.

GATE: Graph CCA for Temporal Self-Supervised Learning for Label-Efficient fMRI Analysis

In this work, we focus on the challenging task, neuro-disease classification, using functional magnetic resonance imaging (fMRI). In population graph-based disease analysis, graph convolutional neural networks (GCNs) have achieved remarkable success. However, these achievements are inseparable from abundant labeled data and sensitive to spurious signals. To improve fMRI representation learning and classification under a label-efficient setting, we propose a novel and theory-driven self-supervised learning (SSL) framework on GCNs, namely Graph CCA for Temporal sElf-supervised learning on fMRI analysis (GATE). Concretely, it is demanding to design a suitable and effective SSL strategy to extract formation and robust features for fMRI. To this end, we investigate several new graph augmentation strategies from fMRI dynamic functional connectives (FC) for SSL training. Further, we leverage canonical-correlation analysis (CCA) on different temporal embeddings and present the theoretical implications. Consequently, this yields a novel two-step GCN learning procedure comprised of (i) SSL on an unlabeled fMRI population graph and (ii) fine-tuning on a small labeled fMRI dataset for a classification task. Our method is tested on two independent fMRI datasets, demonstrating superior performance on autism and dementia diagnosis. Our code is available at https://github.com/LarryUESTC/GATE.

Examining the intersection of cognitive and physical function measures: Results from the brain networks and mobility (B-NET) study

Background and objectives

Although evidence exists that measures of mobility and cognition are correlated, it is not known to what extent they overlap, especially across various domains. This study aimed to investigate the intersection of 18 different objective cognitive and physical function measures from a sample of unimpaired adults aged 70 years and older.

Research design and methods

Canonical correlation analysis was utilized to explore the joint cross-sectional relationship between 13 cognitive and 6 physical function measures in the baseline visit of the Brain Networks and Mobility Function (B-NET) Study (n = 192).

Results

Mean age of participants was 76.4 years. Two synthetic functions were identified. Function 1 explained 26.3% of the shared variability between the cognition and physical function variables, whereas Function 2 explained 19.5%. Function 1 termed "cognitive and physical speed" related the expanded Short Physical Performance Battery (eSPPB), 400-m walk speed, and Dual Task gait speed measures of physical function to semantic fluency animals scores, Digit Symbol Coding (DSC), and Trail Making Test B. Function 2 termed "complex motor tasks and cognitive tasks" related the Force Plate Postural Sway Foam Task and Dual Task to the following cognitive variables: MoCA Adjusted Score, Verbal Fluency L words, Craft story immediate and delayed recall, and Trail Making Test B.

Discussion and implications

We identified groups of cognitive and physical functional abilities that were linked in cross-sectional analyses, which may suggest shared underlying neural network pathway(s) related to speed (Function 1) or complexity (Function 2).

Translational significance

Whether such neural processes decline before measurable functional losses or may be important targets for future interventions that aim to prevent disability also remains to be determined.

Assessment of heavy metals contamination and human health risk assessment of the commonly consumed medicinal herbs in China

This study investigates heavy metal contamination of commonly consumed medicinal herbs and human health risks to the Chinese population arising from the consumption of herbs that contain potentially harmful elements. Food safety standards for Chinese residents are becoming stricter, and much work in this field needs to be performed. This study examines Co, Ba, Fe, Cr, Mn, Ni, Zn, As, Cd, Pb, Cu, Be, Sb, and Bi concentrations in four regularly consumed Chinese herb species: Radix Paeoniae Alba (RPA), Radix Angelicae Dahuricae (RAD), Rhizoma Atractylodis Macrocephalae (RAM), and Radix Puerariae (RP). A pollution status examination and evaluation of heavy metals in RPA, RAD, RAM, and RP were performed. The human health risk assessment associated with the intake of potentially harmful elements in herbs was calculated in terms of the estimated daily intake (EDI), the target hazard quotient (THQ), the estimated hazard index (HI), and the lifetime cancer risk (CR). The mean single-factor pollution index (PI) showed that in the RPA, RAD, RAM, and RP samples, approximately 10.0%, 10.0%, 30.0%, and 10.0%, respectively, were polluted by Cd. The present study indicated that the pattern of consumption of the studied herbs in China does not seem to suggest an excessive health hazard associated with any of the toxic elements studied.

Brain network changes in adult victims of violence

INTRODUCTION

Stressful experiences such as violence can affect mental health severely. The effects are associated with changes in structural and functional brain networks. The current study aimed to investigate brain network changes in four large-scale brain networks, the default mode network, the salience network, the fronto-parietal network, and the dorsal attention network in self-identified victims of violence and controls who did not identify themselves as victims.

MATERIALS AND METHODS

The control group (n = 32) was matched to the victim group (n = 32) by age, gender, and primary psychiatric disorder. Sparse inverse covariance maps were derived from functional resting-state measurements and from T1 weighted structural data for both groups.

RESULTS

Our data underlined that mostly the salience network was affected in the sample of self-identified victims. In self-identified victims with a current psychiatric diagnosis, the dorsal attention network was mostly affected underlining the potential role of psychopathological alterations on attention-related processes.

CONCLUSION

The results showed that individuals who identify themselves as victim demonstrated significant differences in all considered networks, both within- and between-network.

Environmental neuroscience linking exposome to brain structure and function underlying cognition and behavior

Individual differences in human brain structure, function, and behavior can be attributed to genetic variations, environmental exposures, and their interactions. Although genome-wide association studies have identified many genetic variants associated with brain imaging phenotypes, environmental exposures associated with these phenotypes remain largely unknown. Here, we propose that environmental neuroscience should pay more attention on exploring the associations between lifetime environmental exposures (exposome) and brain imaging phenotypes and identifying both cumulative environmental effects and their vulnerable age windows during the life course. Exposome-neuroimaging association studies face several challenges including the accurate measurement of the totality of environmental exposures varied in space and time, the highly correlated structure of the exposome, and the lack of standardized approaches for exposome-wide association studies. By agnostically scanning the effects of environmental exposures on brain imaging phenotypes and their interactions with genomic variations, exposome-neuroimaging association analyses will improve our understanding of causal factors associated with individual differences in brain structure and function as well as their relations with cognitive abilities and neuropsychiatric disorders.

Multimodal and multidomain lesion network mapping enhances prediction of sensorimotor behavior in stroke patients

Beyond the characteristics of a brain lesion, such as its etiology, size or location, lesion network mapping (LNM) has shown that similar symptoms after a lesion reflects similar dis-connectivity patterns, thereby linking symptoms to brain networks. Here, we extend LNM by using a multimodal strategy, combining functional and structural networks from 1000 healthy participants in the Human Connectome Project. We apply multimodal LNM to a cohort of 54 stroke patients with the aim of predicting sensorimotor behavior, as assessed through a combination of motor and sensory tests. Results are two-fold. First, multimodal LNM reveals that the functional modality contributes more than the structural one in the prediction of sensorimotor behavior. Second, when looking at each modality individually, the performance of the structural networks strongly depended on whether sensorimotor performance was corrected for lesion size, thereby eliminating the effect that larger lesions generally produce more severe sensorimotor impairment. In contrast, functional networks provided similar performance regardless of whether or not the effect of lesion size was removed. Overall, these results support the extension of LNM to its multimodal form, highlighting the synergistic and additive nature of different types of network modalities, and their corresponding influence on behavioral performance after brain injury.

Accounting for Bias in the Estimation of r 2 between Two Sets of Noisy Neural Responses

The Pearson correlation coefficient squared, r 2, is an important tool used in the analysis of neural data to quantify the similarity between neural tuning curves. Yet this metric is biased by trial-to-trial variability; as trial-to-trial variability increases, measured correlation decreases. Major lines of research are confounded by this bias, including those involving the study of invariance of neural tuning across conditions and the analysis of the similarity of tuning across neurons. To address this, we extend an estimator, [Formula: see text], that was recently developed for estimating model-to-neuron correlation, in which a noisy signal is compared with a noise-free prediction, to the case of neuron-to-neuron correlation, in which two noisy signals are compared with each other. We compare the performance of our novel estimator to a prior method developed by Spearman, commonly used in other fields but widely overlooked in neuroscience, and find that our method has less bias. We then apply our estimator to demonstrate how it avoids drastic confounds introduced by trial-to-trial variability using data collected in two prior studies (macaque, both sexes) that examined two different forms of invariance in the neural encoding of visual inputs-translation invariance and fill-outline invariance. Our results quantify for the first time the gradual falloff with spatial offset of translation-invariant shape selectivity within visual cortical neuronal receptive fields and offer a principled method to compare invariance in noisy biological systems to that in noise-free models.SIGNIFICANCE STATEMENT Quantifying the similarity between two sets of averaged neural responses is fundamental to the analysis of neural data. A ubiquitous metric of similarity, the correlation coefficient, is attenuated by trial-to-trial variability that arises from many irrelevant factors. Spearman recognized this problem and proposed corrected methods that have been extended over a century. We show this method has large asymptotic biases that can be overcome using a novel estimator. Despite the frequent use of the correlation coefficient in neuroscience, consensus on how to address this fundamental statistical issue has not been reached. We provide an accurate estimator of the correlation coefficient and apply it to gain insight into visual invariance.

Age-related brain atrophy is not a homogenous process: Different functional brain networks associate differentially with aging and blood factors

Aging is characterized by a progressive loss of brain volume at an estimated rate of 5% per decade after age 40. While these morphometric changes, especially those affecting gray matter and atrophy of the temporal lobe, are predictors of cognitive performance, the strong association with aging obscures the potential parallel, but more specific role, of individual subject physiology. Here, we studied a cohort of 554 human subjects who were monitored using structural MRI scans and blood immune protein concentrations. Using machine learning, we derived a cytokine clock (CyClo), which predicted age with good accuracy (Mean Absolute Error = 6 y) based on the expression of a subset of immune proteins. These proteins included, among others, Placenta Growth Factor (PLGF) and Vascular Endothelial Growth Factor (VEGF), both involved in angiogenesis, the chemoattractant vascular cell adhesion molecule 1 (VCAM-1), the canonical inflammatory proteins interleukin-6 (IL-6) and tumor necrosis factor alpha (TNFα), the chemoattractant IP-10 (CXCL10), and eotaxin-1 (CCL11), previously involved in brain disorders. Age, sex, and the CyClo were independently associated with different functionally defined cortical networks in the brain. While age was mostly correlated with changes in the somatomotor system, sex was associated with variability in the frontoparietal, ventral attention, and visual networks. Significant canonical correlation was observed for the CyClo and the default mode, limbic, and dorsal attention networks, indicating that immune circulating proteins preferentially affect brain processes such as focused attention, emotion, memory, response to social stress, internal evaluation, and access to consciousness. Thus, we identified immune biomarkers of brain aging which could be potential therapeutic targets for the prevention of age-related cognitive decline.

Algebraic approach for subspace decomposition and clustering of neural activity

We developed an approach to decompose neuronal signals into disjoint components, corresponding to task- or event-based epochs. This protocol describes how to project behavioral templates onto a low-dimensional subspace of neuronal responses to derive neuronal templates, then how to decompose and cluster neuronal responses using these derived templates. We outline these steps on complementary datasets of calcium imaging and spiking activity. Our approach relies on fundamental, linear algebraic principles and is adaptive to the temporal structure of the neural data.

For complete details on the use and execution of this protocol, please refer to Adam et al. (2022).¹

•

Apply singular value decomposition to derive a low-dimensional neural response subspace

•

Project behavioral templates onto the subspace to derive neuronal templates

•

Every neuronal response can be decomposed as a weighted sum of neuronal templates

•

The derived weights can be used for clustering the neurons

Publisher’s note: Undertaking any experimental protocol requires adherence to local institutional guidelines for laboratory safety and ethics.

Deriving psychiatric symptom-based biomarkers from multivariate relationships between psychophysiological and biochemical measures

Identification of biomarkers for psychiatric disorders remains very challenging due to substantial symptom heterogeneity and diagnostic comorbidity, limiting the ability to map symptoms to underlying neurobiology. Dimensional symptom clusters, such as anhedonia, hyperarousal, etc., are complex and arise due to interactions of a multitude of complex biological relationships. The primary aim of the current investigation was to use multi-set canonical correlation analysis (mCCA) to derive biomarkers (biochemical, physiological) linked to dimensional symptoms across the anxiety and depressive spectrum. Active-duty service members (N = 2,592) completed standardized depression, anxiety and posttraumatic stress questionnaires and several psychophysiological and biochemical assays. Using this approach, we identified two phenotype associations between distinct physiological and biological phenotypes. One was characterized by symptoms of dysphoric arousal (anhedonia, anxiety, hypervigilance) which was associated with low blood pressure and startle reactivity. This finding is in line with previous studies suggesting blunted physiological reactivity is associated with subpopulations endorsing anxiety with comorbid depressive features. A second phenotype of anxious fatigue (high anxiety and reexperiencing/avoidance symptoms coupled with fatigue) was associated with elevated blood levels of norepinephrine and the inflammatory marker C-reactive protein in conjunction with high blood pressure. This second phenotype may describe populations in which inflammation and high sympathetic outflow might contribute to anxious fatigue. Overall, these findings support the growing consensus that distinct neuropsychiatric symptom patterns are associated with differential physiological and blood-based biological profiles and highlight the potential of mCCA to reveal important psychiatric symptom biomarkers from several psychophysiological and biochemical measures.

A cross-cohort replicable and heritable latent dimension linking behaviour to multi-featured brain structure

Identifying associations between interindividual variability in brain structure and behaviour requires large cohorts, multivariate methods, out-of-sample validation and, ideally, out-of-cohort replication. Moreover, the influence of nature vs nurture on brain-behaviour associations should be analysed. We analysed associations between brain structure (grey matter volume, cortical thickness, and surface area) and behaviour (spanning cognition, emotion, and alertness) using regularized canonical correlation analysis and a machine learning framework that tests the generalisability and stability of such associations. The replicability of brain-behaviour associations was assessed in two large, independent cohorts. The load of genetic factors on these associations was analysed with heritability and genetic correlation. We found one heritable and replicable latent dimension linking cognitive-control/executive-functions and positive affect to brain structural variability in areas typically associated with higher cognitive functions, and with areas typically associated with sensorimotor functions. These results revealed a major axis of interindividual behavioural variability linking to a whole-brain structural pattern.

Canonical Correlation Analysis and Partial Least Squares for Identifying Brain-Behavior Associations: A Tutorial and a Comparative Study

Canonical correlation analysis (CCA) and partial least squares (PLS) are powerful multivariate methods for capturing associations across 2 modalities of data (e.g., brain and behavior). However, when the sample size is similar to or smaller than the number of variables in the data, standard CCA and PLS models may overfit, i.e., find spurious associations that generalize poorly to new data. Dimensionality reduction and regularized extensions of CCA and PLS have been proposed to address this problem, yet most studies using these approaches have some limitations. This work gives a theoretical and practical introduction into the most common CCA/PLS models and their regularized variants. We examine the limitations of standard CCA and PLS when the sample size is similar to or smaller than the number of variables. We discuss how dimensionality reduction and regularization techniques address this problem and explain their main advantages and disadvantages. We highlight crucial aspects of the CCA/PLS analysis framework, including optimizing the hyperparameters of the model and testing the identified associations for statistical significance. We apply the described CCA/PLS models to simulated data and real data from the Human Connectome Project and Alzheimer's Disease Neuroimaging Initiative (both of n > 500). We use both low- and high-dimensionality versions of these data (i.e., ratios between sample size and variables in the range of ∼1-10 and ∼0.1-0.01, respectively) to demonstrate the impact of data dimensionality on the models. Finally, we summarize the key lessons of the tutorial.

The mathematical relationship between COVID-19 cases and socio-economic indicators of OECD countries

In most aspects, the COVID-19 pandemic has been the most impactful disease in the last few decades. Although two years have passed since the start of the COVID-19 pandemic, it is clear that the pandemic is far from being brought under full control. This pandemic, which was kept under control in part by the discovery of vaccines, has had a negative impact on people's daily lives as well as global economies. The transmission and death rate of this disease, which was detected in Wuhan, China in December 2019, has differed between countries. The reason for these different rates according to countries has been investigated while developing a mathematical model in this comprehensive study. The effect of 16 different socioeconomic indicators of countries on the parameters of the COVID-19 pandemic was analyzed in this study, which was performed on 35 countries, the majority of which are OECD countries. The main objective of this research is to minimize potential negative effects by improving or taking these indicators into account in future similar outbreaks. In this study, canonical correlation analysis was used to obtain the appropriate mathematical model and the performance of the obtained model has been tested statistically. According to the mathematical model obtained, the age factor comes to the fore as the most important indicator in the COVID-19 pandemic and the next important indicator is the health expenditures of the countries. The influence of 16 different socioeconomic indicators on the COVID-19 pandemic parameters are presented as a percentage with details.