Resolving heterogeneity in schizophrenia through a novel systems approach to brain structure: individualized structural covariance network analysis

Voxel-based texture similarity networks reveal individual variability and correlate with biological ontologies

The human brain is organized as a complex, hierarchical network. However, the structural covariance patterns among brain regions and the underlying biological substrates of such covariance networks remain to be clarified. The present study proposed a novel individualized structural covariance network termed voxel-based texture similarity networks (vTSNs) based on 76 refined voxel-based textural features derived from structural magnetic resonance images. Validated in three independent longitudinal healthy cohorts (40, 23, and 60 healthy participants, respectively) with two common brain atlases, we found that the vTSN could robustly resolve inter-subject variability with high test-retest reliability. In contrast to the regional-based texture similarity networks (rTSNs) that calculate radiomic features based on region-of-interest information, vTSNs had higher inter- and intra-subject variability ratios and test-retest reliability in connectivity strength and network topological properties. Moreover, the Spearman correlation indicated a stronger association of the gene expression similarity network (GESN) with vTSNs than with rTSNs (vTSN: r = 0.600, rTSN: r = 0.433, z = 39.784, P < 0.001). Hierarchical clustering identified 3 vTSN subnets with differential association patterns with 13 coexpression modules, 16 neurotransmitters, 7 electrophysiology, 4 metabolism, and 2 large-scale structural and 4 functional organization maps. Moreover, these subnets had unique biological hierarchical organization from the subcortex-limbic system to the ventral neocortex and then to the dorsal neocortex. Based on 424 unrelated, qualified healthy subjects from the Human Connectome Project, we found that vTSNs could sensitively represent sex differences, especially for connections in the subcortex-limbic system and between the subcortex-limbic system and the ventral neocortex. Moreover, a multivariate variance component model revealed that vTSNs could explain a significant proportion of inter-subject behavioral variance in cognition (80.0 %) and motor functions (63.4 %). Finally, using 494 healthy adults (aged 19-80 years old) from the Southwest University Adult Lifespan Dataset, the Spearman correlation identified a significant association between aging and vTSN strength, especially within the subcortex-limbic system and between the subcortex-limbic system and the dorsal neocortex. In summary, our proposed vTSN is robust in uncovering individual variability and neurobiological brain processes, which can serve as biologically plausible measures for linking biological processes and human behavior.

Individual Structural Covariance Network Predicts Long-Term Motor Improvement in Parkinson Disease with Subthalamic Nucleus Deep Brain Stimulation

BACKGROUND AND PURPOSE

The efficacy of long-term chronic subthalamic nucleus deep brain stimulation (STN-DBS) in treating Parkinson disease (PD) exhibits substantial variability among individuals. The preoperative identification of suitable deep brain stimulation (DBS) candidates through predictive means becomes crucial. Our study aims to investigate the predictive value of characterizing individualized structural covariance networks for long-term efficacy of DBS, offering patients a precise and cost-effective preoperative screening tool.

MATERIALS AND METHODS

We included 138 patients with PD and 40 healthy controls. We developed individualized structural covariance networks from T1-weighted images utilizing network template perturbation, and computed the networks' topological characteristics. Patients were categorized according to their long-term motor improvement following STN-DBS. Intergroup analyses were conducted on individual network edges and topological indices, alongside correlation analyses with long-term outcomes for the entire patient cohort. Finally, machine learning algorithms were employed for regression and classification to predict post-DBS motor improvement.

RESULTS

Among the patients with PD, 6 edges (left middle frontal and left caudate nucleus, right olfactory and right insula, left superior medial frontal gyrus and right insula, right middle frontal and left paracentral lobule, right middle frontal and cerebellum, left lobule VIIb of the cerebellum and the vermis of the cerebellum) exhibited significant results in intergroup comparisons and correlation analyses. Increased degree centrality and local efficiency of the cerebellum, parahippocampal gyrus, and postcentral gyrus were associated with DBS improvement. A regression model constructed from these 6 edges revealed a significant correlation between predicted and observed changes in the unified PD rating scale (R = 0.671, P < .001) and receiver operating characteristic analysis demonstrated an area under the curve of 0.802, effectively distinguishing between patients with good and moderate improvement post-DBS.

CONCLUSIONS

Our findings reveal the link between individual structural covariance network fingerprints in patients with PD and long-term motor outcome following STN-DBS. Additionally, binary and continuous cerebellum-basal ganglia-frontal structural covariance network edges have emerged as potential predictive biomarkers for DBS motor outcome.

Individualized Texture Similarity Network in Schizophrenia

BACKGROUND

Structural covariance network disruption has been considered an important pathophysiological indicator for schizophrenia. Here, we introduced a novel individualized structural covariance network measure, referred to as a texture similarity network (TSN), and hypothesized that the TSN could reliably reveal unique intersubject heterogeneity and complex dysconnectivity patterns in schizophrenia.

METHODS

The TSN was constructed by measuring the covariance of 180 three-dimensional voxelwise gray-level co-occurrence matrix feature maps between brain areas in each participant. We first tested the validity and reproducibility of the TSN in characterizing the intersubject variability in 2 longitudinal test-retest healthy cohorts. The TSN was further applied to elucidate intersubject variability and dysconnectivity patterns in 10 schizophrenia case-control datasets (609 schizophrenia cases vs. 579 controls) as well as in a first-episode depression dataset (69 patients with depression vs. 69 control participants).

RESULTS

The test-retest analysis demonstrated higher TSN intersubject than intrasubject variability. Moreover, the TSN reliably revealed higher intersubject variability in both chronic and first-episode schizophrenia, but not in depression. The TSN also reproducibly detected coexistent increased and decreased TSN strength in widespread brain areas, increased global small-worldness, and the coexistence of both structural hyposynchronization in the central networks and hypersynchronization in peripheral networks in patients with schizophrenia but not in patients with depression. Finally, aberrant intersubject variability and covariance strength patterns revealed by the TSN showed a missing or weak correlation with other individualized structural covariance network measures, functional connectivity, and regional volume changes.

CONCLUSIONS

These findings support the reliability of a TSN in revealing unique structural heterogeneity and complex dysconnectivity in patients with schizophrenia.

Acute limbic system connectivity predicts chronic cognitive function in mild traumatic brain injury: An individualized differential structural covariance network study

Mild traumatic brain injury (mTBI) is a known risk factor for neurodegenerative diseases, yet the precise pathophysiological mechanisms remain poorly understand, often obscured by group-level analysis in non-invasive neuroimaging studies. Individual-based method is critical to exploring heterogeneity in mTBI. We recruited 80 mTBI patients and 40 matched healthy controls, obtaining high-resolution structural MRI for constructing Individual Differential Structural Covariance Networks (IDSCN). Comparisons were conducted at both the individual and group levels. Connectome-based Predictive Modeling (CPM) was applied to predict cognitive performance based on whole-brain connectivity. During the acute stage of mTBI, patients exhibited significant heterogeneity in the count and direction of altered edges, obscured by group-level analysis. In the chronic stage, the number of altered edges decreased and became more consistent, aligning with clinical observations of acute cognitive impairment and gradual improvement. Subgroup analysis based on loss of consciousness/post-traumatic amnesia revealed distinct patterns of alterations. The temporal lobe, particularly regions related to the limbic system, significantly predicted cognitive function from acute to chronic stage. The use of IDSCN and CPM has provided valuable individual-level insights, reconciling discrepancies from previous studies. Additionally, the limbic system may be an appropriate target for future intervention efforts.

Individualized Structural Perturbations on Normative Brain Connectome Restrict Deep Brain Stimulation Outcomes in Parkinson's Disease

BACKGROUND

Patients with Parkinson's disease (PD) respond to deep brain stimulation (DBS) variably. However, how brain substrates restrict DBS outcomes remains unclear.

OBJECTIVE

In this article, we aim to identify prognostic brain signatures for explaining the response variability.

METHODS

We retrospectively investigated a cohort of patients with PD (n = 141) between 2017 and 2022, and defined DBS outcomes as the improvement ratio of clinical motor scores. We used a deviation index to quantify individual perturbations on a reference structural covariance network acquired with preoperative T1-weighted magnetic resonance imaging. The neurobiological perturbations of patients were represented as z scored indices based on the chronological perturbations measured on a group of normal aging adults.

RESULTS

After applying stringent statistical tests (z > 2.5) and correcting for false discoveries (P < 0.01), we found that accelerated deviations mainly affected the prefrontal cortex, motor strip, limbic system, and cerebellum in PD. Particularly, a negative network within the accelerated deviations, expressed as "more preoperative deviations, less postoperative improvements," could predict DBS outcomes (mean absolute error = 0.09, R2 = 0.15). Moreover, a fusion of personal brain predictors and medical responses significantly improved traditional evaluations of DBS outcomes. Notably, the most important brain predictor, a pathway connecting the cognitive unit (prefrontal cortex) and motor control unit (cerebellum and motor strip), partially mediates DBS outcomes with the age at surgery.

CONCLUSIONS

Our findings suggest that individual structural perturbations on the cognitive motor control circuit are critical for modulating DBS outcomes. Interventions toward the circuit have the potential for additional clinical improvements. © 2024 The Author(s). Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Neural basis underlying the association between thought control ability and happiness: The moderating role of the amygdala

Thought control ability (TCA) plays an important role in individuals' health and happiness. Previous studies demonstrated that TCA was closely conceptually associated with happiness. However, empirical research supporting this relationship was limited. In addition, the neural basis underlying TCA and how this neural basis influences the relationship between TCA and happiness remain unexplored. In the present study, the voxel-based morphometry (VBM) method was adopted to investigate the neuroanatomical basis of TCA in 314 healthy subjects. The behavioral results revealed a significant positive association between TCA and happiness. On the neural level, there was a significant negative correlation between TCA and the gray matter density (GMD) of the bilateral amygdala. Split-half validation analysis revealed similar results, further confirming the stability of the VBM analysis findings. Furthermore, gray matter covariance network and graph theoretical analyses showed positive association between TCA and both the node degree and node strength of the amygdala. Moderation analysis revealed that the GMD of the amygdala moderated the relationship between TCA and happiness. Specifically, the positive association between TCA and self-perceived happiness was stronger in subjects with a lower GMD of the amygdala. The present study indicated the neural basis underlying the association between TCA and happiness and offered a method of improving individual well-being.

Neuroimage Analysis Methods and Artificial Intelligence Techniques for Reliable Biomarkers and Accurate Diagnosis of Schizophrenia: Achievements Made by Chinese Scholars Around the Past Decade

BACKGROUND AND HYPOTHESIS

Schizophrenia (SZ) is characterized by significant cognitive and behavioral disruptions. Neuroimaging techniques, particularly magnetic resonance imaging (MRI), have been widely utilized to investigate biomarkers of SZ, distinguish SZ from healthy conditions or other mental disorders, and explore biotypes within SZ or across SZ and other mental disorders, which aim to promote the accurate diagnosis of SZ. In China, research on SZ using MRI has grown considerably in recent years.

STUDY DESIGN

The article reviews advanced neuroimaging and artificial intelligence (AI) methods using single-modal or multimodal MRI to reveal the mechanism of SZ and promote accurate diagnosis of SZ, with a particular emphasis on the achievements made by Chinese scholars around the past decade.

STUDY RESULTS

Our article focuses on the methods for capturing subtle brain functional and structural properties from the high-dimensional MRI data, the multimodal fusion and feature selection methods for obtaining important and sparse neuroimaging features, the supervised statistical analysis and classification for distinguishing disorders, and the unsupervised clustering and semi-supervised learning methods for identifying neuroimage-based biotypes. Crucially, our article highlights the characteristics of each method and underscores the interconnections among various approaches regarding biomarker extraction and neuroimage-based diagnosis, which is beneficial not only for comprehending SZ but also for exploring other mental disorders.

CONCLUSIONS

We offer a valuable review of advanced neuroimage analysis and AI methods primarily focused on SZ research by Chinese scholars, aiming to promote the diagnosis, treatment, and prevention of SZ, as well as other mental disorders, both within China and internationally.

Using individualized structural covariance networks to analyze the heterogeneity of cerebral small vessel disease with cognitive impairment

BACKGROUND

Cerebral small vessel disease (CSVD) includes vascular disorders characterized by heterogeneous pathomechanisms and different neuropathological clinical manifestations. Cognitive dysfunction in CSVD is associated with reductions in structural covariance networks (SCNs). A majority of research conducted on SCNs focused on group-level analysis. However, it is crucial to investigate the individualized variations in order to gain a better understanding of heterogeneous disorders such as CSVD. Therefore, this study aimed to utilize individualized differential structural covariance network (IDSCN) analysis to detect individualized structural covariance aberration.

METHODS

A total of 35 healthy controls and 33 CSVD patients with cognitive impairment participated in this investigation. Using the regional gray matter volume in their T1 images, the IDSCN was constructed for each participant. Finally, the differential structural covariance edges between the two groups were determined by comparing their IDSCN using paired-sample t-tests. On the basis of these differential edges, the two subtypes of cognitively impaired CSVD patients were identified.

RESULTS

The findings revealed that the differential structural covariance edges in CSVD patients with cognitive impairment showed a highly heterogeneous distribution, with the edges primarily cross-distributed between the occipital lobe (specifically inferior occipital gyrus and cuneus), temporal lobe (specifically superior temporal gyrus), and the cerebellum. To varying degrees, the inferior frontal gyrus and the superior parietal gyrus were also distributed. Subsequently, a correlation analysis was performed between the resulting differential edges and the cognitive scale scores. A significant negative association was observed between the cognitive scores and the differential edges distributed in the inferior frontal gyrus and inferior occipital gyrus, the superior temporal gyrus and inferior occipital gyrus, and within the temporal lobe. Particularly in the cognitive domain of attention, the two subtypes separated by differential edges exhibited differences in cognitive scale scores [Mini-Mental State Examination (MMSE) and Montreal Cognitive Assessment (MoCA)]. The differential edges of the subtype 1, characterized by lower cognitive level, were mainly cross-distributed in the limbic lobe (specifically the cingulate gyrus and hippocampus), the parietal lobe (including the superior parietal gyrus and precuneus), and the cerebellum. In contrast, the differential edges of the subtype 2 with a relatively high level of cognition were distributed between the cuneus and the cerebellum.

CONCLUSIONS

The differential structural covariance was investigated between the healthy controls and the CSVD patients with cognitive impairment, showing that differential structural covariance existed between the two groups. The edge distributions in certain parts of the brain, such as cerebellum and occipital and temporal lobes, verified this. Significant associations were seen between cognitive scale scores and some of those differential edges .The two subtypes that differed in both differential edges and cognitive levels were also identified. The differential edges of subtype 1 with relatively lower cognitive levels were more distributed in the cingulate gyrus, hippocampus, superior parietal gyrus, and precuneus. This could potentially offer significant benefits in terms of accurate diagnosis and targeted treatment of heterogeneous disorders such as CSVD.

Neuroimaging epicenters as potential sites of onset of the neuroanatomical pathology in schizophrenia

Schizophrenia lacks a clear definition at the neuroanatomical level, capturing the sites of origin and progress of this disorder. Using a network-theory approach called epicenter mapping on cross-sectional magnetic resonance imaging from 1124 individuals with schizophrenia, we identified the most likely "source of origin" of the structural pathology. Our results suggest that the Broca's area and adjacent frontoinsular cortex may be the epicenters of neuroanatomical pathophysiology in schizophrenia. These epicenters can predict an individual's response to treatment for psychosis. In addition, cross-diagnostic similarities based on epicenter mapping over of 4000 individuals diagnosed with neurological, neurodevelopmental, or psychiatric disorders appear to be limited. When present, these similarities are restricted to bipolar disorder, major depressive disorder, and obsessive-compulsive disorder. We provide a comprehensive framework linking schizophrenia-specific epicenters to multiple levels of neurobiology, including cognitive processes, neurotransmitter receptors and transporters, and human brain gene expression. Epicenter mapping may be a reliable tool for identifying the potential onset sites of neural pathophysiology in schizophrenia.

Acoustic assessment in mandarin-speaking Parkinson's disease patients and disease progression monitoring and brain impairment within the speech subsystem

Approximately 90% of Parkinson's patients (PD) suffer from dysarthria. However, there is currently a lack of research on acoustic measurements and speech impairment patterns among Mandarin-speaking individuals with PD. This study aims to assess the diagnosis and disease monitoring possibility in Mandarin-speaking PD patients through the recommended speech paradigm for non-tonal languages, and to explore the anatomical and functional substrates. We examined total of 160 native Mandarin-speaking Chinese participants consisting of 80 PD patients, 40 healthy controls (HC), and 40 MRI controls. We screened the optimal acoustic metric combination for PD diagnosis. Finally, we used the objective metrics to predict the patient's motor status using the Naïve Bayes model and analyzed the correlations between cortical thickness, subcortical volumes, functional connectivity, and network properties. Comprehensive acoustic screening based on prosodic, articulation, and phonation abnormalities allows differentiation between HC and PD with an area under the curve of 0.931. Patients with slowed reading exhibited atrophy of the fusiform gyrus (FDR p = 0.010, R = 0.391), reduced functional connectivity between the fusiform gyrus and motor cortex, and increased nodal local efficiency (NLE) and nodal efficiency (NE) in bilateral pallidum. Patients with prolonged pauses demonstrated atrophy in the left hippocampus, along with decreased NLE and NE. The acoustic assessment in Mandarin proves effective in diagnosis and disease monitoring for Mandarin-speaking PD patients, generalizing standardized acoustic guidelines beyond non-tonal languages. The speech impairment in Mandarin-speaking PD patients not only involves motor aspects of speech but also encompasses the cognitive processes underlying language generation.

Aberrant individual structure covariance network in patients with mesial temporal lobe epilepsy

Objective

Mesial temporal lobe epilepsy (mTLE) is a complex neurological disorder that has been recognized as a widespread global network disorder. The group-level structural covariance network (SCN) could reveal the structural connectivity disruption of the mTLE but could not reflect the heterogeneity at the individual level.

Methods

This study adopted a recently proposed individual structural covariance network (IDSCN) method to clarify the alternated structural covariance connection mode in mTLE and to associate IDSCN features with the clinical manifestations and regional brain atrophy.

Results

We found significant IDSCN abnormalities in the ipsilesional hippocampus, ipsilesional precentral gyrus, bilateral caudate, and putamen in mTLE patients than in healthy controls. Moreover, the IDSCNs of these areas were positively correlated with the gray matter atrophy rate. Finally, we identified several connectivities with weak associations with disease duration, frequency, and surgery outcome.

Significance

Our research highlights the role of hippo-thalamic-basal-cortical circuits in the pathophysiologic process of disrupted whole-brain morphological covariance networks in mTLE, and builds a bridge between brain-wide covariance network changes and regional brain atrophy.

Altered individual gray matter structural covariance networks in early abstinence patients with alcohol dependence

While alterations in cortical thickness have been widely observed in individuals with alcohol dependence, knowledge about cortical thickness-based structural covariance networks is limited. This study aimed to explore the topological disorganization of structural covariance networks based on cortical thickness at the single-subject level among patients with alcohol dependence. Structural imaging data were obtained from 61 patients with alcohol dependence during early abstinence and 59 healthy controls. The single-subject structural covariance networks were constructed based on cortical thickness data from 68 brain regions and were analyzed using graph theory. The relationships between network architecture and clinical characteristics were further investigated using partial correlation analysis. In the structural covariance networks, both patients with alcohol dependence and healthy controls displayed small-world topology. However, compared to controls, alcohol-dependent individuals exhibited significantly altered global network properties characterized by greater normalized shortest path length, greater shortest path length, and lower global efficiency. Patients exhibited lower degree centrality and nodal efficiency, primarily in the right precuneus. Additionally, scores on the Alcohol Use Disorder Identification Test were negatively correlated with the degree centrality and nodal efficiency of the left middle temporal gyrus. The results of this correlation analysis did not survive after multiple comparisons in the exploratory analysis. Our findings may reveal alterations in the topological organization of gray matter networks in alcoholism patients, which may contribute to understanding the mechanisms of alcohol addiction from a network perspective.

Gray matter atrophy is constrained by normal structural brain network architecture in depression

BACKGROUND

There is growing evidence that gray matter atrophy is constrained by normal brain network (or connectome) architecture in neuropsychiatric disorders. However, whether this finding holds true in individuals with depression remains unknown. In this study, we aimed to investigate the association between gray matter atrophy and normal connectome architecture at individual level in depression.

METHODS

In this study, 297 patients with depression and 256 healthy controls (HCs) from two independent Chinese dataset were included: a discovery dataset (105 never-treated first-episode patients and matched 130 HCs) and a replication dataset (106 patients and matched 126 HCs). For each patient, individualized regional atrophy was assessed using normative model and brain regions whose structural connectome profiles in HCs most resembled the atrophy patterns were identified as putative epicenters using a backfoward stepwise regression analysis.

RESULTS

In general, the structural connectome architecture of the identified disease epicenters significantly explained 44% (±16%) variance of gray matter atrophy. While patients with depression demonstrated tremendous interindividual variations in the number and distribution of disease epicenters, several disease epicenters with higher participation coefficient than randomly selected regions, including the hippocampus, thalamus, and medial frontal gyrus were significantly shared by depression. Other brain regions with strong structural connections to the disease epicenters exhibited greater vulnerability. In addition, the association between connectome and gray matter atrophy uncovered two distinct subgroups with different ages of onset.

CONCLUSIONS

These results suggest that gray matter atrophy is constrained by structural brain connectome and elucidate the possible pathological progression in depression.

Changes in Structural Covariance among Olfactory-related Brain Regions in Anosmia Patients

Anosmia, characterized by the loss of smell, is associated not only with dysfunction in the peripheral olfactory system but also with changes in several brain regions involved in olfactory processing. Specifically, the orbitofrontal cortex is recognized for its pivotal role in integrating olfactory information, engaging in bidirectional communication with the primary olfactory regions, including the olfactory cortex, amygdala, and entorhinal cortex. However, little is known about alterations in structural connections among these brain regions in patients with anosmia. In this study, high-resolution T1-weighted images were obtained from participants. Utilizing the volumes of key brain regions implicated in olfactory function, we employed a structural covariance approach to investigate brain reorganization patterns in patients with anosmia (n=22) compared to healthy individuals (n=30). Our structural covariance analysis demonstrated diminished connectivity between the amygdala and entorhinal cortex, components of the primary olfactory network, in patients with anosmia compared to healthy individuals (z=-2.22, FDR-corrected p=0.039). Conversely, connectivity between the orbitofrontal cortex-a major region in the extended olfactory network-and amygdala was found to be enhanced in the anosmia group compared to healthy individuals (z=2.32, FDR-corrected p=0.039). However, the structural connections between the orbitofrontal cortex and entorhinal cortex did not differ significantly between the groups (z=0.04, FDR-corrected p=0.968). These findings suggest a potential structural reorganization, particularly of higher-order cortical regions, possibly as a compensatory effort to interpret the limited olfactory information available in individuals with olfactory loss.

Analysis of brain structural covariance network in Cushing disease

Background

Cushing disease (CD) is a rare clinical neuroendocrine disease. CD is characterized by abnormal hypercortisolism induced by a pituitary adenoma with the secretion of adrenocorticotropic hormone. Individuals with CD usually exhibit atrophy of gray matter volume. However, little is known about the alterations in topographical organization of individuals with CD. This study aimed to investigate the structural covariance networks of individuals with CD based on the gray matter volume using graph theory analysis.

Methods

High-resolution T1-weighted images of 61 individuals with CD and 53 healthy controls were obtained. Gray matter volume was estimated and the structural covariance network was analyzed using graph theory. Network properties such as hubs of all participants were calculated based on degree centrality.

Results

No significant differences were observed between individuals with CD and healthy controls in terms of age, gender, and education level. The small-world features were conserved in individuals with CD but were higher than those in healthy controls. The individuals with CD showed higher global efficiency and modularity, suggesting higher integration and segregation as compared to healthy controls. The hub nodes of the individuals with CD were Short insular gyri (G_insular_short_L), Anterior part of the cingulate gyrus and sulcus (G_and_S_cingul-Ant_R), and Superior frontal gyrus (G_front_sup_R).

Conclusions

Significant differences in the structural covariance network of patients with CD were found based on graph theory. These findings might help understanding the pathogenesis of individuals with CD and provide insight into the pathogenesis of this CD.

Altered cerebellar and caudate gray-matter volumes and structural covariance networks preceding dual cognitive and mobility impairments in older people

INTRODUCTION

The neuroanatomical changes driving both cognitive and mobility impairments, an emerging preclinical dementia syndrome, are not fully understood. We examined gray-matter volumes (GMVs) and structural covariance networks (SCNs) abnormalities in community-based older people preceding the conversion to physio-cognitive decline syndrome (PCDS).

METHODS

Voxel-wise brain GMV and established SCNs were compared between PCDS and non-PCDS converters.

RESULTS

The study included 343 individuals (60.2 ± 6.9 years, 49.6% men) with intact cognitive and mobility functions. Over an average 5.6-year follow-up, 116 transitioned to PCDS. Identified regions with abnormal GMVs in PCDS converters were over cerebellum and caudate, which served as seeds for SCNs establishment. Significant differences in cerebellum-based (to right frontal pole and left middle frontal gyrus) and caudate-based SCNs (to right caudate putamen, right planum temporale, left precentral gyrus, right postcentral gyrus, and left parietal operculum) between converters and nonconverters were observed.

DISCUSSION

This study reveals early neuroanatomic changes, emphasizing the cerebellum's role, in dual cognitive and mobility impairments.

HIGHLIGHTS

Neuroanatomic precursors of dual cognitive and mobility impairments are identified. Cerebellar GMV reductions and increased right caudate GMV precede the onset of PCDS. Altered cerebellum- and caudate-based SCNs drive PCDS transformation. This research establishes a foundation for understanding PCDS as a specific dementia syndrome.

Detecting differential transcript usage in complex diseases with SPIT

Differential transcript usage (DTU) plays a crucial role in determining how gene expression differs among cells, tissues, and developmental stages, contributing to the complexity and diversity of biological systems. In abnormal cells, it can also lead to deficiencies in protein function and underpin disease pathogenesis. Analyzing DTU via RNA sequencing (RNA-seq) data is vital, but the genetic heterogeneity in populations with complex diseases presents an intricate challenge due to diverse causal events and undetermined subtypes. Although the majority of common diseases in humans are categorized as complex, state-of-the-art DTU analysis methods often overlook this heterogeneity in their models. We therefore developed SPIT, a statistical tool that identifies predominant subgroups in transcript usage within a population along with their distinctive sets of DTU events. This study provides comprehensive assessments of SPIT's methodology and applies it to analyze brain samples from individuals with schizophrenia, revealing previously unreported DTU events in six candidate genes.

Associations Between Structural Covariance Network and Antipsychotic Treatment Response in Schizophrenia

BACKGROUND AND HYPOTHESIS

Schizophrenia is associated with widespread cortical thinning and abnormality in the structural covariance network, which may reflect connectome alterations due to treatment effect or disease progression. Notably, patients with treatment-resistant schizophrenia (TRS) have stronger and more widespread cortical thinning, but it remains unclear whether structural covariance is associated with treatment response in schizophrenia.

STUDY DESIGN

We organized a multicenter magnetic resonance imaging study to assess structural covariance in a large population of TRS and non-TRS, who had been resistant and responsive to non-clozapine antipsychotics, respectively. Whole-brain structural covariance for cortical thickness was assessed in 102 patients with TRS, 77 patients with non-TRS, and 79 healthy controls (HC). Network-based statistics were used to examine the difference in structural covariance networks among the 3 groups. Moreover, the relationship between altered individual differentiated structural covariance and clinico-demographics was also explored.

STUDY RESULTS

Patients with non-TRS exhibited greater structural covariance compared with HC, mainly in the fronto-temporal and fronto-occipital regions, while there were no significant differences in structural covariance between TRS and non-TRS or HC. Higher individual differentiated structural covariance was associated with lower general scores of the Positive and Negative Syndrome Scale in the non-TRS group, but not in the TRS group.

CONCLUSIONS

These findings suggest that reconfiguration of brain networks via coordinated cortical thinning is related to treatment response in schizophrenia. Further longitudinal studies are warranted to confirm if greater structural covariance could serve as a marker for treatment response in this disease.

Resolving heterogeneity in Alzheimer's disease based on individualized structural covariance network

The heterogeneity of Alzheimer's disease (AD) poses a challenge to precision medicine. We aimed to identify distinct subtypes of AD based on the individualized structural covariance network (IDSCN) analysis and to research the underlying neurobiology mechanisms. In this study, 187 patients with AD (age = 73.57 ± 6.00, 50% female) and 143 matched normal controls (age = 74.30 ± 7.80, 44% female) were recruited from the Alzheimer's Disease Neuroimaging Initiative (ADNI) project database, and T1 images were acquired. We utilized the IDSCN analysis to generate individual-level altered structural covariance network and performed k-means clustering to subtype AD based on structural covariance network. Cognition, disease progression, morphological features, and gene expression profiles were further compared between subtypes, to characterize the heterogeneity in AD. Two distinct AD subtypes were identified in a reproducible manner, and we named the two subtypes as slow progression type (subtype 1, n = 104, age = 76.15 ± 6.44, 42% female) and rapid progression type (subtype 2, n = 83, age = 71.98 ± 8.72, 47% female), separately. Subtype 1 had better baseline visuospatial function than subtype 2 (p < 0.05), whereas subtype 2 had better baseline memory function than subtype 1 (p < 0.05). Subtype 2 showed worse progression in memory (p = 0.003), language (p = 0.003), visuospatial function (p = 0.020), and mental state (p = 0.038) than subtype 1. Subtype 1 often shared increased structural covariance network, mainly in the frontal lobe and temporal lobe regions, whereas subtype 2 often shared increased structural covariance network, mainly in occipital lobe regions and temporal lobe regions. Functional annotation further revealed that all differential structural covariance network between the two AD subtypes were mainly implicated in memory, learning, emotion, and cognition. Additionally, differences in gray matter volume (GMV) between AD subtypes were identified, and genes associated with GMV differences were found to be enriched in the terms potassium ion transport, synapse organization, and histone modification and the pathways viral infection, neurodegeneration-multiple diseases, and long-term depression. The two distinct AD subtypes were identified and characterized with neuroanatomy, cognitive trajectories, and gene expression profiles. These comprehensive results have implications for neurobiology mechanisms and precision medicine.

Toward individualized connectomes of brain morphology

The morphological brain connectome (MBC) delineates the coordinated patterns of local morphological features (such as cortical thickness) across brain regions. While classically constructed using population-based approaches, there is a growing trend toward individualized modeling. Currently, the methods for individualized MBCs are varied, posing challenges for method selection and cross-study comparisons. Here, we summarize how individualized MBCs are modeled through low-order methods (correlation-, divergence-, distance-, and deviation-based methods) describing relations in brain morphology, as well as high-order methods capturing similarities in these low-order relations. We discuss the merits and limitations of different methods, examining them in the context of robustness, reproducibility, and reliability. We highlight the importance of elucidating the cellular and molecular mechanisms underlying the individualized connectome, and establishing normative benchmarks to assess individual variation in development, aging, and neuropsychiatric disorders.

Identifying psychosis subtypes use individualized covariance structural differential networks and multi-site clustering

BACKGROUND

Similarities among schizophrenia (SZ), schizoaffective disorder (SAD) and bipolar disorder (BP) including clinical phenotypes, brain alterations and risk genes, make it challenging to perform reliable separation among them. However, previous subtype identification that transcend traditional diagnostic boundaries were based on group-level neuroimaging features, ignoring individual-level inferences.

METHODS

455 psychoses (178 SZs, 134 SADs and 143 BPs), their first-degree relatives (N = 453) and healthy controls (HCs, N = 220) were collected from Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP I) consortium. Individualized covariance structural differential networks (ICSDNs) were constructed for each patient and multi-site clustering was used to identify psychosis subtypes. Group differences between subtypes in clinical phenotypes and voxel-wise fractional amplitude of low frequency fluctuation (fALFF) were calculated, as well as between the corresponding relatives.

RESULTS

Two psychosis subtypes were identified with increased whole brain structural covariance, with decreased connectivity between amygdala-hippocampus and temporal-occipital cortex in subtype I (S-I) compared to subtype II (S-II), which was replicated under different clustering methods, number of edges and across datasets (B-SNIP II) and different brain atlases. S-I had higher emotional-related symptoms than S-II and showed significant fALFF decrease in temporal and occipital cortex, while S-II was more similar to HC. This pattern was consistently validated on relatives of S-I and S-II in both fALFF and clinical symptoms.

CONCLUSIONS

These findings reconcile categorical and dimensional perspectives of psychosis neurobiological heterogeneity, indicating that relatives of S-I might have greater predisposition in developing psychosis, while relatives of S-II are more likely to be healthy. This study contributes to the development of neuroimaging informed diagnostic classifications within psychosis spectrum.

Effects of PRRT2 mutation on brain gray matter networks in paroxysmal kinesigenic dyskinesia

Although proline-rich transmembrane protein 2 is the primary causative gene of paroxysmal kinesigenic dyskinesia, its effects on the brain structure of paroxysmal kinesigenic dyskinesia patients are not yet clear. Here, we explored the influence of proline-rich transmembrane protein 2 mutations on similarity-based gray matter morphological networks in individuals with paroxysmal kinesigenic dyskinesia. A total of 51 paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations, 55 paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation, and 80 healthy controls participated in the study. We analyzed the structural connectome characteristics across groups by graph theory approaches. Relative to paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation and healthy controls, paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations exhibited a notable increase in characteristic path length and a reduction in both global and local efficiency. Relative to healthy controls, both patient groups showed reduced nodal metrics in right postcentral gyrus, right angular, and bilateral thalamus; Relative to healthy controls and paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation, paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations showed almost all reduced nodal centralities and structural connections in cortico-basal ganglia-thalamo-cortical circuit including bilateral supplementary motor area, bilateral pallidum, and right caudate nucleus. Finally, we used support vector machine by gray matter network matrices to classify paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 mutations and paroxysmal kinesigenic dyskinesia patients possessing proline-rich transmembrane protein 2 non-mutation, achieving an accuracy of 73%. These results show that proline-rich transmembrane protein 2 related gray matter network deficits may contribute to paroxysmal kinesigenic dyskinesia, offering new insights into its pathophysiological mechanisms.

Integrative omics analysis reveals epigenomic and transcriptomic signatures underlying brain structural deficits in major depressive disorder

Several lines of evidence support the involvement of transcriptomic and epigenetic mechanisms in the brain structural deficits of major depressive disorder (MDD) separately. However, research in these two areas has remained isolated. In this study, we proposed an integrative strategy that combined neuroimaging, brain-wide gene expression, and peripheral DNA methylation data to investigate the genetic basis of gray matter abnormalities in MDD. The MRI T1-weighted images and Illumina 850 K DNA methylation microarrays were obtained from 269 patients and 416 healthy controls, and brain-wide transcriptomic data were collected from Allen Human Brain Atlas. The between-group differences in gray matter volume (GMV) and differentially methylated CpG positions (DMPs) were examined. The genes with their expression patterns spatially related to GMV changes and genes with DMPs were overlapped and selected. Using principal component regression, the associations between DMPs in overlapped genes and GMV across individual patients were investigated, and the region-specific correlations between methylation status and gene expression were examined. We found significant associations between the decreased GMV and DMPs methylation status in the anterior cingulate cortex, inferior frontal cortex, and fusiform face cortex regions. These DMPs genes were primarily enriched in the neurodevelopmental and synaptic transmission process. There was a significant negative correlation between DNA methylation and gene expression in genes associated with GMV changes of the frontal cortex in MDD. Our findings suggest that GMV abnormalities in MDD may have a transcriptomic and epigenetic basis. This imaging-transcriptomic-epigenetic integrative analysis provides spatial and biological links between cortical morphological deficits and peripheral epigenetic signatures in MDD.

Metabolic interactions between organs in overweight and obesity using total-body positron emission tomography

BACKGROUND AND OBJECTIVES

Overweight and obesity is a complex condition resulting from unbalanced energy homeostasis among various organs. However, systemic abnormalities in overweight and obese people are seldom explored in vivo by metabolic imaging techniques. The aim of this study was to determine metabolic abnormities throughout the body in overweight and obese adults using total-body positron emission tomography (PET) glucose uptake imaging.

METHODS

Thirty normal weight subjects [body mass index (BMI) < 25 kg/m2, 55.47 ± 13.94 years, 16 men and 14 women], and 26 overweight and obese subjects [BMI ≥ 25 kg/m2, 52.38 ± 9.52 years, 21 men and 5 women] received whole-body 18F-fluorodeoxyglucose PET imaging using the uEXPLORER. Whole-body standardized uptake value normalized by lean body mass (SUL) images were calculated. Metabolic networks were constructed based on the whole-body SUL images using covariance network approach. Both group-level and individual-level network differences between normal weight and overweight/obese subjects were evaluated. Correlation analysis was conducted between network properties and BMI for the overweight/obese subjects.

RESULTS

Compared with normal weight subjects, overweight/obese subjects exhibited altered network connectivity strength in four network nodes, namely the pancreas (p = 0.033), spleen (p = 0.021), visceral adipose tissue (VAT) (p = 1.12 × 10-5) and bone (p = 0.021). Network deviations of overweight/obese subjects from the normal weight were positively correlated with BMI (r = 0.718, p = 3.64 × 10-5). In addition, overweight/obese subjects experienced altered connections between organs, and some of the altered connections, including pancreas-right lung and VAT-bilateral lung connections were significantly correlated with BMI.

CONCLUSION

Overweight/obese individuals exhibit metabolic alterations in organ level, and altered metabolic interactions at the systemic level. The proposed approach using total-body PET imaging can reveal individual metabolic variability and metabolic deviations between organs, which would open up a new path for understanding metabolic alterations in overweight and obesity.

Individual-level brain morphological similarity networks: Current methodologies and applications

AIMS

The human brain is an extremely complex system in which neurons, clusters of neurons, or regions are connected to form a complex network. With the development of neuroimaging techniques, magnetic resonance imaging (MRI)-based brain networks play a key role in our understanding of the intricate architecture of human brain. Among them, the structural MRI-based brain morphological network approach has attracted increasing attention due to the advantages in data acquisition, image quality, and in revealing the structural organizing principles intrinsic to the brain. This review is to summarize the methodology and related applications of individual-level morphological networks.

BACKGROUND

There have been a growing number of studies related to brain morphological similarity networks. Conventional morphological networks are intersubject covariance networks constructed using a certain morphological indicator of a group of subjects; individual-level morphological networks, on the other hand, measure the morphological similarity between brain regions for individual brains and can reflect the morphological information of single subjects. In recent years, individual morphological networks have demonstrated significant worth in exploring the topological changes of the human brain under both normal and disease conditions. Such studies provided novel perspectives for understanding human brain development and exploring the pathological mechanisms of neuropsychiatric disorders.

CONCLUSION

This paper mainly focuses on the studies of brain morphological networks at the individual level, introduces several ways for network construction, reviews representative work in this field, and finally points out current problems and future directions.

Disrupted brain gray matter connectome in social anxiety disorder: a novel individualized structural covariance network analysis

Phenotyping approaches grounded in structural network science can offer insights into the neurobiological substrates of psychiatric diseases, but this remains to be clarified at the individual level in social anxiety disorder (SAD). Using a recently developed approach combining probability density estimation and Kullback-Leibler divergence, we constructed single-subject structural covariance networks (SCNs) based on multivariate morphometry (cortical thickness, surface area, curvature, and volume) and quantified their global/nodal network properties using graph-theoretical analysis. We compared network metrics between SAD patients and healthy controls (HC) and analyzed the relationship to clinical characteristics. We also used support vector machine analysis to explore the ability of graph-theoretical metrics to discriminate SAD patients from HC. Globally, SAD patients showed higher global efficiency, shorter characteristic path length, and stronger small-worldness. Locally, SAD patients showed abnormal nodal centrality mainly involving left superior frontal gyrus, right superior parietal lobe, left amygdala, right paracentral gyrus, right lingual, and right pericalcarine cortex. Altered topological metrics were associated with the symptom severity and duration. Graph-based metrics allowed single-subject classification of SAD versus HC with total accuracy of 78.7%. This finding, that the topological organization of SCNs in SAD patients is altered toward more randomized configurations, adds to our understanding of network-level neuropathology in SAD.

Resolving heterogeneity in depression using individualized structural covariance network analysis

BACKGROUND

Elucidating individual aberrance is a critical first step toward precision medicine for heterogeneous disorders such as depression. The neuropathology of depression is related to abnormal inter-regional structural covariance indicating a brain maturational disruption. However, most studies focus on group-level structural covariance aberrance and ignore the interindividual heterogeneity. For that reason, we aimed to identify individualized structural covariance aberrance with the help of individualized differential structural covariance network (IDSCN) analysis.

METHODS

T1-weighted anatomical images of 195 first-episode untreated patients with depression and matched healthy controls (n = 78) were acquired. We obtained IDSCN for each patient and identified subtypes of depression based on shared differential edges.

RESULTS

As a result, patients with depression demonstrated tremendous heterogeneity in the distribution of differential structural covariance edges. Despite this heterogeneity, altered edges within subcortical-cerebellum network were often shared by most of the patients. Two robust neuroanatomical subtypes were identified. Specifically, patients in subtype 1 often shared decreased motor network-related edges. Patients in subtype 2 often shared decreased subcortical-cerebellum network-related edges. Functional annotation further revealed that differential edges in subtype 2 were mainly implicated in reward/motivation-related functional terms.

CONCLUSIONS

In conclusion, we investigated individualized differential structural covariance and identified that decreased edges within subcortical-cerebellum network are often shared by patients with depression. The identified two subtypes provide new insights into taxonomy and facilitate potential clues to precision diagnosis and treatment of depression.

Altered topology of individual brain structural covariance networks in major depressive disorder

BACKGROUND

The neurobiological pathogenesis of major depression disorder (MDD) remains largely controversial. Previous literatures with limited sample size utilizing group-level structural covariance networks (SCN) commonly generated mixed findings regarding the topology of brain networks.

METHODS

We analyzed T1 images from a high-powered multisite sample including 1173 patients with MDD and 1019 healthy controls (HCs). We used regional gray matter volume to construct individual SCN by utilizing a novel approach based on the interregional effect size difference. We further investigated MDD-related structural connectivity alterations using topological metrics.

RESULTS

Compared to HCs, the MDD patients showed a shift toward randomization characterized by increased integration. Further subgroup analysis of patients in different stages revealed this randomization pattern was also observed in patients with recurrent MDD, while the first-episode drug naïve patients exhibited decreased segregation. Altered nodal properties in several brain regions which have a key role in both emotion regulation and executive control were also found in MDD patients compared with HCs. The abnormalities in inferior temporal gyrus were not influenced by any specific site. Moreover, antidepressants increased nodal efficiency in the anterior ventromedial prefrontal cortex.

CONCLUSIONS

The MDD patients at different stages exhibit distinct patterns of randomization in their brain networks, with increased integration during illness progression. These findings provide valuable insights into the disruption in structural brain networks that occurs in patients with MDD and might be useful to guide future therapeutic interventions.

Detecting differential transcript usage in complex diseases with SPIT

Differential transcript usage (DTU) plays a crucial role in determining how gene expression differs among cells, tissues, and different developmental stages, thereby contributing to the complexity and diversity of biological systems. In abnormal cells, it can also lead to deficiencies in protein function, potentially leading to pathogenesis of diseases. Detecting such events for single-gene genetic traits is relatively uncomplicated; however, the heterogeneity of populations with complex diseases presents an intricate challenge due to the presence of diverse causal events and undetermined subtypes. SPIT is the first statistical tool that quantifies the heterogeneity in transcript usage within a population and identifies predominant subgroups along with their distinctive sets of DTU events. We provide comprehensive assessments of SPIT's methodology in both single-gene and complex traits and report the results of applying SPIT to analyze brain samples from individuals with schizophrenia. Our analysis reveals previously unreported DTU events in six candidate genes.

A reliable subtyping of de novo Parkinson disease: biomarkers, medication effects and longitudinal progression

Divergent clinical symptoms and pathological progression suggest multiple subtypes of Parkinson disease (PD). Here, we proposed a reliable PD subtyping approach that quantifies the disturbance of an individual patient to the reference structural covariance networks derived from healthy controls. We revealed two subtypes of de novo PD patients by using longitudinal data from the PPMI dataset. Compared to the conventional clinical TD/PIGD phenotypes, our subtyping was highly stable in 5 years' visits. The two subtypes of PD showed significant differences in motor symptoms, medication effects, CSF biomarkers, and longitudinal progression. Moreover, patients of subtype 2 showed widespread lower cortical-to-dorsal raphe nucleus (DRN) connections and higher medication effects on motor symptoms which was regulated by 5-HT neurons in DRN. Our results suggest distinct neuropathological pathways underlying the two subtypes, such that, in contrast to the typical PD subtype, patients of subtype 2 may be affected by serotonergic modulation on dopaminergic neurons in striatum. Our study opens new avenue to precision medicine and personalized treatments in PD and may be applicable to other neurodegenerative diseases.

ICA-based Individualized Differential Structure Similarity Networks for Predicting Symptom Scores in Adolescents with Major Depressive Disorder

Major depressive disorder (MDD) is a complex mood disorder characterized by persistent and overwhelming depression. Previous studies have identified large scale structural brain alterations in MDD, yet most are group analyses with atlas-parcellated anatomical regions. Here we proposed a method to measure individual difference by independent component analysis (ICA)-based individual difference structural similarity network (IDSSN). This approach provided a data-adaptive, atlas-free solution that can be applied to new individual subjects. Specifically, we constructed individualized whole-brain structural covariance networks based on network perturbation approach using spatially constrained ICA. First, a set of benchmark independent components (ICs) were generated using gray matter volume (GMV) from all healthy controls. Then individual heterogeneity was obtained by calculating differences and other similarity metrics between ICs derived from "each one patient + all controls" and the benchmark ICs, resulting in 32 imaging features and structural similarity networks for each patient, which can be used for predicting multiple clinical symptoms. We estimated IDSSN for 189 adolescent MDD patients aged 10-20 years and compared them to 112 healthy adolescents. We tested their predictability of the Hamilton Anxiety Scale , the 17-item Hamilton Depression Scale and six clinical syndromes using connectome-based predictive modeling. The prediction results showed that ICA-based IDSSN features reveal more disease-specific information than those using other brain templates. We also found that depression-associated networks mainly involved the default-mode network and visual network. In conclusion, our study proposed an adaptive method that improves the ability to detect GMV deviations and specificity between one individual patient and healthy groups, providing a new perspectives and insights for evaluating individual-level clinical heterogeneity based on brain structures.

Structural covariances of prefrontal subregions selectively associate with dopamine-related gene coexpression and schizophrenia

Evidence highlights that dopamine (DA) system dysregulation and prefrontal cortex (PFC) dysfunction may underlie the pathophysiology of schizophrenia. However, the associations among DA genes, PFC morphometry, and schizophrenia have not yet been fully clarified. Based on the brain gene expression dataset from Allen Human Brain Atlas and structural magnetic resonance imaging data (NDIS = 1727, NREP = 408), we first identified 10 out of 22 PFC subregions whose gray matter volume (GMV) covariance profiles were reliably associated with their DA genes coexpression profiles, then four out of the identified 10 PFC subregions demonstrated abnormally increased GMV covariance with the hippocampus, insula, and medial frontal areas in schizophrenia patients (NCASE = 100; NCONTROL = 102). Moreover, based on a schizophrenia postmortem expression dataset, we found that the DA genes coexpression of schizophrenia was significantly reduced between the middle frontal gyrus and hippocampus, in which 21 DA genes showed significantly unsynchronized expression changes, and the 21 genes' brain expression were enriched in brain activity invoked by working memory, reward, speech production, and episodic memory. Our findings indicate the DA genes selectively regulate the structural covariance of PFC subregions by their coexpression profiles, which may underlie the disrupted GMV covariance and impaired cognitive functions in schizophrenia.

A multicentre study on grey matter morphometric biomarkers for classifying early schizophrenia and parkinson's disease psychosis

Psychotic symptoms occur in a majority of schizophrenia patients and in ~50% of all Parkinson's disease (PD) patients. Altered grey matter (GM) structure within several brain areas and networks may contribute to their pathogenesis. Little is known, however, about transdiagnostic similarities when psychotic symptoms occur in different disorders, such as in schizophrenia and PD. The present study investigated a large, multicenter sample containing 722 participants: 146 patients with first episode psychosis, FEP; 106 individuals in at-risk mental state for developing psychosis, ARMS; 145 healthy controls matching FEP and ARMS, Con-Psy; 92 PD patients with psychotic symptoms, PDP; 145 PD patients without psychotic symptoms, PDN; 88 healthy controls matching PDN and PDP, Con-PD. We applied source-based morphometry in association with receiver operating curves (ROC) analyses to identify common GM structural covariance networks (SCN) and investigated their accuracy in identifying the different patient groups. We assessed group-specific homogeneity and variability across the different networks and potential associations with clinical symptoms. SCN-extracted GM values differed significantly between FEP and Con-Psy, PDP and Con-PD, PDN and Con-PD, as well as PDN and PDP, indicating significant overall grey matter reductions in PD and early schizophrenia. ROC analyses showed that SCN-based classification algorithms allow good classification (AUC ~0.80) of FEP and Con-Psy, and fair performance (AUC ~0.72) when differentiating PDP from Con-PD. Importantly, the best performance was found in partly the same networks, including the thalamus. Alterations within selected SCNs may be related to the presence of psychotic symptoms in both early schizophrenia and PD psychosis, indicating some commonality of underlying mechanisms. Furthermore, results provide evidence that GM volume within specific SCNs may serve as a biomarker for identifying FEP and PDP.

Modular architecture and resilience of structural covariance networks in first-episode antipsychotic-naive psychoses

Structural covariance network (SCN) studies on first-episode antipsychotic-naïve psychosis (FEAP) have examined less granular parcellations on one morphometric feature reporting lower network resilience among other findings. We examined SCNs of volume, cortical thickness, and surface area using the Human Connectome Project atlas-based parcellation (n = 358 regions) from 79 FEAP and 68 controls to comprehensively characterize the networks using a descriptive and perturbational network neuroscience approach. Using graph theoretical methods, we examined network integration, segregation, centrality, community structure, and hub distribution across the small-worldness threshold range and correlated them with psychopathology severity. We used simulated nodal "attacks" (removal of nodes and all their edges) to investigate network resilience, calculated DeltaCon similarity scores, and contrasted the removed nodes to characterize the impact of simulated attacks. Compared to controls, FEAP SCN showed higher betweenness centrality (BC) and lower degree in all three morphometric features and disintegrated with fewer attacks with no change in global efficiency. SCNs showed higher similarity score at the first point of disintegration with ≈ 54% top-ranked BC nodes attacked. FEAP communities consisted of fewer prefrontal, auditory and visual regions. Lower BC, and higher clustering and degree, were associated with greater positive and negative symptom severity. Negative symptoms required twice the changes in these metrics. Globally sparse but locally dense network with more nodes of higher centrality in FEAP could result in higher communication cost compared to controls. FEAP network disintegration with fewer attacks suggests lower resilience without impacting efficiency. Greater network disarray underlying negative symptom severity possibly explains the therapeutic challenge.

Mapping the neuroanatomical heterogeneity of OCD using a framework integrating normative model and non-negative matrix factorization

Obsessive-compulsive disorder (OCD) is a spectrum disorder with high interindividual heterogeneity. We propose a comprehensible framework integrating normative model and non-negative matrix factorization (NMF) to quantitatively estimate the neuroanatomical heterogeneity of OCD from a dimensional perspective. T1-weighted magnetic resonance images of 98 first-episode untreated patients with OCD and matched healthy controls (HCs, n = 130) were acquired. We derived individualized differences in gray matter morphometry using normative model and parsed them into latent disease factors using NMF. Four robust disease factors were identified. Each patient expressed multiple factors and exhibited a unique factor composition. Factor compositions of patients were significantly correlated with severity of symptom, age of onset, illness duration, and exhibited sex differences, capturing sources of clinical heterogeneity. In addition, the group-level morphological differences obtained with two-sample t test could be quantitatively derived from the identified disease factors, reconciling the group-level and subject-level findings in neuroimaging studies. Finally, we uncovered two distinct subtypes with opposite morphological differences compared with HCs from factor compositions. Our findings suggest that morphological differences of individuals with OCD are the unique combination of distinct neuroanatomical patterns. The proposed framework quantitatively estimating neuroanatomical heterogeneity paves the way for precision medicine in OCD.

Spatial navigation is associated with subcortical alterations and progression risk in subjective cognitive decline

BACKGROUND

Subjective cognitive decline (SCD) may serve as a symptomatic indicator for preclinical Alzheimer's disease; however, SCD is a heterogeneous entity regarding clinical progression. We aimed to investigate whether spatial navigation could reveal subcortical structural alterations and the risk of progression to objective cognitive impairment in SCD individuals.

METHODS

One hundred and eighty participants were enrolled: those with SCD (n = 80), normal controls (NCs, n = 77), and mild cognitive impairment (MCI, n = 23). SCD participants were further divided into the SCD-good (G-SCD, n = 40) group and the SCD-bad (B-SCD, n = 40) group according to their spatial navigation performance. Volumes of subcortical structures were calculated and compared among the four groups, including basal forebrain, thalamus, caudate, putamen, pallidum, hippocampus, amygdala, and accumbens. Topological properties of the subcortical structural covariance network were also calculated. With an interval of 1.5 years ± 12 months of follow-up, the progression rate to MCI was compared between the G-SCD and B-SCD groups.

RESULTS

Volumes of the basal forebrain, the right hippocampus, and their respective subfields differed significantly among the four groups (p < 0.05, false discovery rate corrected). The B-SCD group showed lower volumes in the basal forebrain than the G-SCD group, especially in the Ch4p and Ch4a-i subfields. Furthermore, the structural covariance network of the basal forebrain and right hippocampal subfields showed that the B-SCD group had a larger Lambda than the G-SCD group, which suggested weakened network integration in the B-SCD group. At follow-up, the B-SCD group had a significantly higher conversion rate to MCI than the G-SCD group.

CONCLUSION

Compared to SCD participants with good spatial navigation performance, SCD participants with bad performance showed lower volumes in the basal forebrain, a reorganized structural covariance network of subcortical nuclei, and an increased risk of progression to MCI. Our findings indicated that spatial navigation may have great potential to identify SCD subjects at higher risk of clinical progression, which may contribute to making more precise clinical decisions for SCD individuals who seek medical help.

Large-scale analysis of structural brain asymmetries in schizophrenia via the ENIGMA consortium

Left-right asymmetry is an important organizing feature of the healthy brain that may be altered in schizophrenia, but most studies have used relatively small samples and heterogeneous approaches, resulting in equivocal findings. We carried out the largest case-control study of structural brain asymmetries in schizophrenia, with MRI data from 5,080 affected individuals and 6,015 controls across 46 datasets, using a single image analysis protocol. Asymmetry indexes were calculated for global and regional cortical thickness, surface area, and subcortical volume measures. Differences of asymmetry were calculated between affected individuals and controls per dataset, and effect sizes were meta-analyzed across datasets. Small average case-control differences were observed for thickness asymmetries of the rostral anterior cingulate and the middle temporal gyrus, both driven by thinner left-hemispheric cortices in schizophrenia. Analyses of these asymmetries with respect to the use of antipsychotic medication and other clinical variables did not show any significant associations. Assessment of age- and sex-specific effects revealed a stronger average leftward asymmetry of pallidum volume between older cases and controls. Case-control differences in a multivariate context were assessed in a subset of the data (N = 2,029), which revealed that 7% of the variance across all structural asymmetries was explained by case-control status. Subtle case-control differences of brain macrostructural asymmetry may reflect differences at the molecular, cytoarchitectonic, or circuit levels that have functional relevance for the disorder. Reduced left middle temporal cortical thickness is consistent with altered left-hemisphere language network organization in schizophrenia.

Identification of shared and distinct patterns of brain network abnormality across mental disorders through individualized structural covariance network analysis

BACKGROUND

Mental disorders, including depression, obsessive compulsive disorder (OCD), and schizophrenia, share a common neuropathy of disturbed large-scale coordinated brain maturation. However, high-interindividual heterogeneity hinders the identification of shared and distinct patterns of brain network abnormalities across mental disorders. This study aimed to identify shared and distinct patterns of altered structural covariance across mental disorders.

METHODS

Subject-level structural covariance aberrance in patients with mental disorders was investigated using individualized differential structural covariance network. This method inferred structural covariance aberrance at the individual level by measuring the degree of structural covariance in patients deviating from matched healthy controls (HCs). T1-weighted anatomical images of 513 participants (105, 98, 190 participants with depression, OCD and schizophrenia, respectively, and 130 age- and sex-matched HCs) were acquired and analyzed.

RESULTS

Patients with mental disorders exhibited notable heterogeneity in terms of altered edges, which were otherwise obscured by group-level analysis. The three disorders shared high difference variability in edges attached to the frontal network and the subcortical-cerebellum network, and they also exhibited disease-specific variability distributions. Despite notable variability, patients with the same disorder shared disease-specific groups of altered edges. Specifically, depression was characterized by altered edges attached to the subcortical-cerebellum network; OCD, by altered edges linking the subcortical-cerebellum and motor networks; and schizophrenia, by altered edges related to the frontal network.

CONCLUSIONS

These results have potential implications for understanding heterogeneity and facilitating personalized diagnosis and interventions for mental disorders.

Infections During Pregnancy and Risks for Adult Psychosis: Findings from the New England Family Study

For the past 40 years, our team has conducted a unique program of research investigating the prenatal risks for schizophrenia and related adult psychiatric disorders. The New England Family Study is a long-term prospective cohort study of over 16,000 individuals followed from the prenatal period for over 50 years. This chapter summarizes several major phases and findings from this work, highlighting recent results on maternal prenatal bacterial infections and brain imaging. Implications regarding the causes and potential prevention of major psychotic disorders are discussed.

Schizophrenia and psychedelic state: Dysconnection versus hyper-connection. A perspective on two different models of psychosis stemming from dysfunctional integration processes

Psychotic symptoms are a cross-sectional dimension affecting multiple diagnostic categories, despite schizophrenia represents the prototype of psychoses. Initially, dopamine was considered the most involved molecule in the neurobiology of schizophrenia. Over the next years, several biological factors were added to the discussion helping to constitute the concept of schizophrenia as a disease marked by a deficit of functional integration, contributing to the formulation of the Dysconnection Hypothesis in 1995. Nowadays the notion of dysconnection persists in the conceptualization of schizophrenia enriched by neuroimaging findings which corroborate the hypothesis. At the same time, in recent years, psychedelics received a lot of attention by the scientific community and astonishing findings emerged about the rearrangement of brain networks under the effect of these compounds. Specifically, a global decrease in functional connectivity was found, highlighting the disintegration of preserved and functional circuits and an increase of overall connectivity in the brain. The aim of this paper is to compare the biological bases of dysconnection in schizophrenia with the alterations of neuronal cyto-architecture induced by psychedelics and the consequent state of cerebral hyper-connection. These two models of psychosis, despite diametrically opposed, imply a substantial deficit of integration of neural signaling reached through two opposite paths.

Leveraging Machine Learning for Gaining Neurobiological and Nosological Insights in Psychiatric Research

Much attention is currently devoted to developing diagnostic classifiers for mental disorders. Complementing these efforts, we highlight the potential of machine learning to gain biological insights into the psychopathology and nosology of mental disorders. Studies to this end have mainly used brain imaging data, which can be obtained noninvasively from large cohorts and have repeatedly been argued to reveal potentially intermediate phenotypes. This may become particularly relevant in light of recent efforts to identify magnetic resonance imaging-derived biomarkers that yield insight into pathophysiological processes as well as to refine the taxonomy of mental illness. In particular, the accuracy of machine learning models may be used as dependent variables to identify features relevant to pathophysiology. Moreover, such approaches may help disentangle the dimensional (within diagnosis) and often overlapping (across diagnoses) symptomatology of psychiatric illness. We also point out a multiview perspective that combines data from different sources, bridging molecular and system-level information. Finally, we summarize recent efforts toward a data-driven definition of subtypes or disease entities through unsupervised and semisupervised approaches. The latter, blending unsupervised and supervised concepts, may represent a particularly promising avenue toward dissecting heterogeneous categories. Finally, we raise several technical and conceptual aspects related to the reviewed approaches. In particular, we discuss common pitfalls pertaining to flawed input data or analytic procedures that would likely lead to unreliable outputs.

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.

Structural covariance network of the hippocampus-amygdala complex in medication-naïve patients with first-episode major depressive disorder

Background

The hippocampus and amygdala are densely interconnected structures that work together in multiple affective and cognitive processes that are important to the etiology of major depressive disorder (MDD). Each of these structures consists of several heterogeneous subfields. We aim to explore the topologic properties of the volume-based intrinsic network within the hippocampus-amygdala complex in medication-naïve patients with first-episode MDD.

Methods

High-resolution T1-weighted magnetic resonance imaging scans were acquired from 123 first-episode, medication-naïve, and noncomorbid MDD patients and 81 age-, sex-, and education level-matched healthy control participants (HCs). The structural covariance network (SCN) was constructed for each group using the volumes of the hippocampal subfields and amygdala subregions; the weights of the edges were defined by the partial correlation coefficients between each pair of subfields/subregions, controlled for age, sex, education level, and intracranial volume. The global and nodal graph metrics were calculated and compared between groups.

Results

Compared with HCs, the SCN within the hippocampus-amygdala complex in patients with MDD showed a shortened mean characteristic path length, reduced modularity, and reduced small-worldness index. At the nodal level, the left hippocampal tail showed increased measures of centrality, segregation, and integration, while nodes in the left amygdala showed decreased measures of centrality, segregation, and integration in patients with MDD compared with HCs.

Conclusion

Our results provide the first evidence of atypical topologic characteristics within the hippocampus-amygdala complex in patients with MDD using structure network analysis. It provides more delineate mechanism of those two structures that underlying neuropathologic process in MDD.

Stage-dependent differential influence of metabolic and structural networks on memory across Alzheimer's disease continuum

Background

Large-scale neuronal network breakdown underlies memory impairment in Alzheimer's disease (AD). However, the differential trajectories of the relationships between network organisation and memory across pathology and cognitive stages in AD remain elusive. We determined whether and how the influences of individual-level structural and metabolic covariance network integrity on memory varied with amyloid pathology across clinical stages without assuming a constant relationship.

Methods

Seven hundred and eight participants from the Alzheimer's Disease Neuroimaging Initiative were studied. Individual-level structural and metabolic covariance scores in higher-level cognitive and hippocampal networks were derived from magnetic resonance imaging and [18F] fluorodeoxyglucose positron emission tomography using seed-based partial least square analyses. The non-linear associations between network scores and memory across cognitive stages in each pathology group were examined using sparse varying coefficient modelling.

Results

We showed that the associations of memory with structural and metabolic networks in the hippocampal and default mode regions exhibited pathology-dependent differential trajectories across cognitive stages using sparse varying coefficient modelling. In amyloid pathology group, there was an early influence of hippocampal structural network deterioration on memory impairment in the preclinical stage, and a biphasic influence of the angular gyrus-seeded default mode metabolic network on memory in both preclinical and dementia stages. In non-amyloid pathology groups, in contrast, the trajectory of the hippocampus-memory association was opposite and weaker overall, while no metabolism covariance networks were related to memory. Key findings were replicated in a larger cohort of 1280 participants.

Conclusions

Our findings highlight potential windows of early intervention targeting network breakdown at the preclinical AD stage.

Funding

Data collection and sharing for this project was funded by the Alzheimer's Disease Neuroimaging Initiative (ADNI) (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). We also acknowledge the funding support from the Duke NUS/Khoo Bridge Funding Award (KBrFA/2019-0020) and NMRC Open Fund Large Collaborative Grant (OFLCG09May0035).

SPAMRI: A MATLAB Toolbox for Surface-Based Processing and Analysis of Magnetic Resonance Imaging

Structural magnetic resonance imaging (MRI) has elicited increasing attention in morphological surface studies due to its stability and sensitivity to neurodegenerative processes, particularly in exploring brain aging and psychiatric disease. However, a user-friendly toolbox for the surface-based analysis of structural MRI is still lacking. On the basis of certain software functions in FreeSurfer, CAT and ANTs, a MATLAB toolbox called "surface-based processing and analysis of MRI" (SPAMRI) has been developed, which can be performed in Windows, Linux and Mac-OS. SPAMRI contains several features as follows: (1) open-source MATLAB-based package with a graphical user interface (GUI); (2) a set of images that can be generated for quality checking, such as Talairach transform, skull strip, and surface reconstruction; (3) user-friendly GUI with capabilities on statistical analysis, multiple comparison corrections, reporting of results, and surface measurement extraction; and (4) provision of a conversion tool between surface files (e.g., mesh files) and volume files (e.g., NIFTI files). SPAMRI is applied to a publicly released structural MRI dataset of 44 healthy young adults and 39 old adults. Findings showed that old people have decreased cortical thickness, especially in prefrontal cortex, relative to those of young adults, thereby suggesting a cognitive decline in the former. SPAMRI is anticipated to substantially simplify surface-based image processing and MRI dataset analyses and subsequently open new opportunities to investigate structural morphologies.

Two distinct subtypes of obsessive compulsive disorder revealed by a framework integrating multimodal neuroimaging information

Patients with obsessive compulsive disorder (OCD) exhibit tremendous heterogeneity in structural and functional neuroimaging aberrance. However, most previous studies just focus on group‐level aberrance of a single modality ignoring heterogeneity and multimodal features. On that account, we aimed to uncover OCD subtypes integrating structural and functional neuroimaging features with the help of a multiview learning method and examined multimodal aberrance for each subtype. Ninety‐nine first‐episode untreated patients with OCD and 104 matched healthy controls (HCs) undergoing structural and functional MRI were included in this study. Voxel‐based morphometric and amplitude of low‐frequency fluctuation (ALFF) were adopted to assess gray matter volumes (GMVs) and the spontaneous neuronal fluctuations respectively. Structural/functional distance network was obtained by calculating Euclidean distance between pairs of regional GMVs/ALFF values across patients. Similarity network fusion, one of multiview learning methods capturing shared and complementary information from multimodal data sources, was used to fuse multimodal distance networks into one fused network. Then spectral clustering was adopted to categorize patients into subtypes. As a result, two robust subtypes were identified. These two subtypes presented opposite GMV aberrance and distinct ALFF aberrance compared with HCs while shared indistinguishable clinical and demographic features. In addition, these two subtypes exhibited opposite structure–function difference correlation reflecting distinct adaptive modifications between multimodal aberrance. Altogether, these results uncover two objective subtypes with distinct multimodal aberrance and provide a new insight into taxonomy of OCD.

Identifying the individual metabolic abnormities from a systemic perspective using whole-body PET imaging

Introduction

Distinct physiological states arise from complex interactions among the various organs present in the human body. PET is a non-invasive modality with numerous successful applications in oncology, neurology, and cardiology. However, while PET imaging has been applied extensively in detecting focal lesions or diseases, its potential in detecting systemic abnormalities is seldom explored, mostly because total-body imaging was not possible until recently.

Methods

In this context, the present study proposes a framework capable of constructing an individual metabolic abnormality network using a subject’s whole-body ¹⁸F-FDG SUV image and a normal control database. The developed framework was evaluated in the patients with lung cancer, the one discharged after suffering from Covid-19 disease, and the one that had gastrointestinal bleeding with the underlying cause unknown.

Results

The framework could successfully capture the deviation of these patients from healthy subjects at the level of both system and organ. The strength of the altered network edges revealed the abnormal metabolic connection between organs. The overall deviation of the network nodes was observed to be highly correlated to the organ SUV measures. Therefore, the molecular connectivity of glucose metabolism was characterized at a single subject level.

Conclusion

The proposed framework represents a significant step toward the use of PET imaging for identifying metabolic dysfunction from a systemic perspective. A better understanding of the underlying biological mechanisms and the physiological interpretation of the interregional connections identified in the present study warrant further research.

Two distinct subtypes of obsessive compulsive disorder revealed by heterogeneity through discriminative analysis

Neurobiological heterogeneity in obsessive compulsive disorder (OCD) is understudied leading to conflicting neuroimaging findings. Therefore, we investigated objective neuroanatomical subtypes of OCD by adopting a newly proposed method based on gray matter volumes (GMVs). GMVs were derived from T1‐weighted anatomical images of patients with OCD (n = 100) and matched healthy controls (HCs; n = 106). We first inquired whether patients with OCD presented higher interindividual variability HCs in terms of GMVs. Then, we identified distinct subtypes of OCD by adopting heterogeneity through discriminative analysis (HYDRA), where regional GMVs were treated as features. Patients with OCD presented higher interindividual variability than HCs, suggesting a high structural heterogeneity of OCD. HYDRA identified two distinct robust subtypes of OCD presenting opposite neuroanatomical aberrances compared with HCs, while sharing indistinguishable clinical and demographic features. Specifically, Subtype 1 exhibited widespread increased GMVs in cortical and subcortical regions, including the orbitofrontal gyrus, right anterior insula, bilateral hippocampus, and bilateral parahippocampus and cerebellum. Subtype 2 demonstrated overall decreased GMVs in regions such as the orbitofrontal gyrus, right anterior insula, and precuneus. When mixed together, none of patients presented significant differences compared with HCs. In addition, the total intracranial volume of Subtype 2 was significantly correlated with the total score of the Yale–Brown Obsessive Compulsive Scale while that of Subtype 1 was not. These results identified two distinct neuroanatomical subtypes, providing a possible explanation for conflicting neuroimaging findings, and proposed a potential objective taxonomy contributing to precise clinical diagnosis and treatment in OCD.

Cell type-specific manifestations of cortical thickness heterogeneity in schizophrenia

Brain morphology differs markedly between individuals with schizophrenia, but the cellular and genetic basis of this heterogeneity is poorly understood. Here, we sought to determine whether cortical thickness (CTh) heterogeneity in schizophrenia relates to interregional variation in distinct neural cell types, as inferred from established gene expression data and person-specific genomic variation. This study comprised 1849 participants in total, including a discovery (140 cases and 1267 controls) and a validation cohort (335 cases and 185 controls). To characterize CTh heterogeneity, normative ranges were established for 34 cortical regions and the extent of deviation from these ranges was measured for each individual with schizophrenia. CTh deviations were explained by interregional gene expression levels of five out of seven neural cell types examined: (1) astrocytes; (2) endothelial cells; (3) oligodendrocyte progenitor cells (OPCs); (4) excitatory neurons; and (5) inhibitory neurons. Regional alignment between CTh alterations with cell type transcriptional maps distinguished broad patient subtypes, which were validated against genomic data drawn from the same individuals. In a predominantly neuronal/endothelial subtype (22% of patients), CTh deviations covaried with polygenic risk for schizophrenia (sczPRS) calculated specifically from genes marking neuronal and endothelial cells (r = -0.40, p = 0.010). Whereas, in a predominantly glia/OPC subtype (43% of patients), CTh deviations covaried with sczPRS calculated from glia and OPC-linked genes (r = -0.30, p = 0.028). This multi-scale analysis of genomic, transcriptomic, and brain phenotypic data may indicate that CTh heterogeneity in schizophrenia relates to inter-individual variation in cell-type specific functions. Decomposing heterogeneity in relation to cortical cell types enables prioritization of schizophrenia subsets for future disease modeling efforts.