The thalamus and its subnuclei-a gateway to obsessive-compulsive disorder

Effects of gene dosage and development on subcortical nuclei volumes in individuals with 22q11.2 copy number variations

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in a total of 385 scans from 22qDel (n = 96, scans = 191, 53.1% female), 22qDup (n = 37, scans = 64, 45.9% female), and TD controls (n = 80, scans = 130, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the linear effects of 22q11.2 gene dosage and non-linear effects of age were characterized with generalized additive mixed models (GAMMs). Positive gene dosage effects (volume increasing with copy number) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Distinct Volume Alterations of Thalamic Nuclei Across the Schizophrenia Spectrum

BACKGROUND AND HYPOTHESIS

Abnormal thalamic nuclei volumes and their link to cognitive impairments have been observed in schizophrenia. However, whether and how this finding extends to the schizophrenia spectrum is unknown. We hypothesized a distinct pattern of aberrant thalamic nuclei volume across the spectrum and examined its potential associations with cognitive symptoms.

STUDY DESIGN

We performed a FreeSurfer-based volumetry of T1-weighted brain MRIs from 137 healthy controls, 66 at-risk mental state (ARMS) subjects, 89 first-episode psychosis (FEP) individuals, and 126 patients with schizophrenia to estimate thalamic nuclei volumes of six nuclei groups (anterior, lateral, ventral, intralaminar, medial, and pulvinar). We used linear regression models, controlling for sex, age, and estimated total intracranial volume, both to compare thalamic nuclei volumes across groups and to investigate their associations with positive, negative, and cognitive symptoms.

STUDY RESULTS

We observed significant volume alterations in medial and lateral thalamic nuclei. Medial nuclei displayed consistently reduced volumes across the spectrum compared to controls, while lower lateral nuclei volumes were only observed in schizophrenia. Whereas positive and negative symptoms were not associated with reduced nuclei volumes across all groups, higher cognitive scores were linked to lower volumes of medial nuclei in ARMS. In FEP, cognition was not linked to nuclei volumes. In schizophrenia, lower cognitive performance was associated with lower medial volumes.

CONCLUSIONS

Results demonstrate distinct thalamic nuclei volume reductions across the schizophrenia spectrum, with lower medial nuclei volumes linked to cognitive deficits in ARMS and schizophrenia. Data suggest a distinctive trajectory of thalamic nuclei abnormalities along the course of schizophrenia.

Exploring volumetric abnormalities in subcortical L-HPA axis structures in pediatric generalized anxiety disorder

BACKGROUND

Pediatric generalized anxiety disorder (GAD) is debilitating and increasingly prevalent, yet its etiology remains unclear. Some believe the disorder to be propagated by chronic dysregulation of the limbic-hypothalamic-pituitary-adrenal (L-HPA) axis, but morphometric studies of implicated subcortical areas have been largely inconclusive. Recognizing that certain subcortical subdivisions are more directly involved in L-HPA axis functioning, this study aims to detect specific abnormalities in these critical areas.

METHODS

Thirty-eight MRI scans of preschool children with (n = 15) and without (n = 23) GAD underwent segmentation and between-group volumetric comparisons of the basolateral amygdala (BLA), ventral hippocampal subiculum (vSC), and mediodorsal medial magnocellular (MDm) area of the thalamus.

RESULTS

Children with GAD displayed significantly larger vSC compared to healthy peers, F(1, 31) = 6.50, pFDR = .048. On average, children with GAD presented with larger BLA and MDm, Fs(1, 31) ≥ 4.86, psFDR ≤ .054. Exploratory analyses revealed right-hemispheric lateralization of all measures, most notably the MDm, F(1, 31) = 8.13, pFDR = .024, the size of which scaled with symptom severity, r = .83, pFDR = .033.

CONCLUSION

The BLA, vSC, and MDm are believed to be involved in the regulation of anxiety and stress, both individually and collectively through the excitation and inhibition of the L-HPA axis. All were found to be enlarged in children with GAD, perhaps reflecting hypertrophy related to hyperexcitability, or early neuronal overgrowth. Longitudinal studies should investigate the relationship between these early morphological differences and the long-term subcortical atrophy previously observed.

Structural and resting-state connection abnormalities of habenula in obsessive-compulsive disorder

BACKGROUND

Previous studies have suggested that the habenula (Hb) may be involved in the mechanism of obsessive-compulsive disorder (OCD). However, the specific role of Hb in OCD remains unclear. This study aimed to explore the structural and functional abnormalities of Hb in OCD and their relationship with the clinical symptoms.

METHODS

Eighty patients with OCD and 85 healthy controls (HCs) were recruited as the primary dataset. The grey matter volume, resting-state functional connectivity (FC), and effective connectivity (EC) of the Hb were calculated and compared between OCD group and HCs. An independent replication dataset was used to verify the stability and robustness of the results.

RESULTS

Patients with OCD exhibited smaller Hb volume and increased FC of right Hb-left hippocampus than HCs. Dynamic causal model revealed an increased EC from left hippocampus to right Hb and a less inhibitory causal influence from the right Hb to left hippocampus in the OCD group compared to HCs. Similar results were found in the replication dataset.

CONCLUSIONS

This study suggested that abnormal structure of Hb and hippocampus-Hb connectivity may contribute to the pathological basis of OCD.

Thalamic nuclei volume partially mediates the effects of aerobic capacity on fatigue in people with multiple sclerosis

BACKGROUND

Fatigue is frequent in people with multiple sclerosis (pwMS) impacting physical and cognitive functions. Lower aerobic capacity and regional thalamic volume may be involved in the pathophysiology of fatigue in pwMS.

OBJECTIVES

To identify associations between thalamic nuclei volumes, aerobic capacity and fatigue and to investigate whether the influence of aerobic capacity on fatigue in pwMS is mediated by thalamic integrity.

METHODS

Eighty-three pwMS underwent a clinical evaluation with assessment of fatigue (Modified Fatigue Impact Scale [MFIS]), including physical (pMFIS) and cognitive (cMFIS) components, and peak of oxygen uptake (VO2peak). PwMS and 63 sex- and age-matched healthy controls (HC) underwent a 3 T brain MRI to quantify volume of the whole thalamus and its nuclei.

RESULTS

Compared to HC, pwMS showed higher global MFIS, pMFIS and cMFIS scores, and lower VO2peak and thalamic volumes (p < 0.001). In pwMS, higher VO2peak was significantly associated with lower MFIS and pMFIS scores (r value = - 0.326 and - 0.356; pFDR ≤ 0.046) and higher laterodorsal thalamic nucleus (Dor) cluster volume (r value = 0.300; pFDR = 0.047). Moreover, lower Dor thalamic cluster volume was significantly associated with higher MFIS, pMFIS and cMFIS scores (r value range = - 0.305; - 0.293; pFDR ≤ 0.049). The volume of Dor thalamic cluster partially mediated the positive effects of VO2peak on both MFIS and cMFIS, with relative indirect effects of 21% and 32% respectively. No mediation was found for pMFIS.

CONCLUSIONS

Higher VO2peak is associated with lower fatigue in pwMS, likely acting on Dor thalamic cluster volume integrity. Such an effect might be different according to the type of fatigue (cognitive or physical).

Aberrant dynamic functional connectivity of thalamocortical circuitry in major depressive disorder

Thalamocortical circuitry has a substantial impact on emotion and cognition. Previous studies have demonstrated alterations in thalamocortical functional connectivity (FC), characterized by region-dependent hypo- or hyper-connectivity, among individuals with major depressive disorder (MDD). However, the dynamical reconfiguration of the thalamocortical system over time and potential abnormalities in dynamic thalamocortical connectivity associated with MDD remain unclear. Hence, we analyzed dynamic FC (dFC) between ten thalamic subregions and seven cortical subnetworks from resting-state functional magnetic resonance images of 48 patients with MDD and 57 healthy controls (HCs) to investigate time-varying changes in thalamocortical FC in patients with MDD. Moreover, dynamic laterality analysis was conducted to examine the changes in functional lateralization of the thalamocortical system over time. Correlations between the dynamic measures of thalamocortical FC and clinical assessment were also calculated. We identified four dynamic states of thalamocortical circuitry wherein patients with MDD exhibited decreased fractional time and reduced transitions within a negative connectivity state that showed strong correlations with primary cortical networks, compared with the HCs. In addition, MDD patients also exhibited increased fluctuations in functional laterality in the thalamocortical system across the scan duration. The thalamo-subnetwork analysis unveiled abnormal dFC variability involving higher-order cortical networks in the MDD cohort. Significant correlations were found between increased dFC variability with dorsal attention and default mode networks and the severity of symptoms. Our study comprehensively investigated the pattern of alteration of the thalamocortical dFC in MDD patients. The heterogeneous alterations of dFC between the thalamus and both primary and higher-order cortical networks may help characterize the deficits of sensory and cognitive processing in MDD.

Role of autonomic, nociceptive, and limbic brainstem nuclei in core autism features

Although multiple theories have speculated about the brainstem reticular formation's involvement in autistic behaviors, the in vivo imaging of brainstem nuclei needed to test these theories has proven technologically challenging. Using methods to improve brainstem imaging in children, this study set out to elucidate the role of the autonomic, nociceptive, and limbic brainstem nuclei in the autism features of 145 children (74 autistic children, 6.0-10.9 years). Participants completed an assessment of core autism features and diffusion- and T1-weighted imaging optimized to improve brainstem images. After data reduction via principal component analysis, correlational analyses examined associations among autism features and the microstructural properties of brainstem clusters. Independent replication was performed in 43 adolescents (24 autistic, 13.0-17.9 years). We found specific nuclei, most robustly the parvicellular reticular formation-alpha (PCRtA) and to a lesser degree the lateral parabrachial nucleus (LPB) and ventral tegmental parabrachial pigmented complex (VTA-PBP), to be associated with autism features. The PCRtA and some of the LPB associations were independently found in the replication sample, but the VTA-PBP associations were not. Consistent with theoretical perspectives, the findings suggest that individual differences in pontine reticular formation nuclei contribute to the prominence of autistic features. Specifically, the PCRtA, a nucleus involved in mastication, digestion, and cardio-respiration in animal models, was associated with social communication in children, while the LPB, a pain-network nucleus, was associated with repetitive behaviors. These findings highlight the contributions of key autonomic brainstem nuclei to the expression of core autism features.

Age of onset of obsessive-compulsive disorder differentially affects white matter microstructure

Previous diffusion MRI studies have reported mixed findings on white matter microstructure alterations in obsessive-compulsive disorder (OCD), likely due to variation in demographic and clinical characteristics, scanning methods, and underpowered samples. The OCD global study was created across five international sites to overcome these challenges by harmonizing data collection to identify consistent brain signatures of OCD that are reproducible and generalizable. Single-shell diffusion measures (e.g., fractional anisotropy), multi-shell Neurite Orientation Dispersion and Density Imaging (NODDI) and fixel-based measures, were extracted from skeletonized white matter tracts in 260 medication-free adults with OCD and 252 healthy controls. We additionally performed structural connectome analysis. We compared cases with controls and cases with early (<18) versus late (18+) OCD onset using mixed-model and Bayesian multilevel analysis. Compared with healthy controls, adult OCD individuals showed higher fiber density in the sagittal stratum (B[SE] = 0.10[0.05], P = 0.04) and credible evidence for higher fiber density in several other tracts. When comparing early (n = 145) and late-onset (n = 114) cases, converging evidence showed lower integrity of the posterior thalamic radiation -particularly radial diffusivity (B[SE] = 0.28[0.12], P = 0.03)-and lower global efficiency of the structural connectome (B[SE] = 15.3[6.6], P = 0.03) in late-onset cases. Post-hoc analyses indicated divergent direction of effects of the two OCD groups compared to healthy controls. Age of OCD onset differentially affects the integrity of thalamo-parietal/occipital tracts and the efficiency of the structural brain network. These results lend further support for the role of the thalamus and its afferent fibers and visual attentional processes in the pathophysiology of OCD.

Exploring the multidimensional nature of repetitive and restricted behaviors and interests (RRBI) in autism: neuroanatomical correlates and clinical implications

BACKGROUND

Repetitive and restricted behaviors and interests (RRBI) are core symptoms of autism with a complex entity and are commonly categorized into 'motor-driven' and 'cognitively driven'. RRBI symptomatology depends on the individual's clinical environment limiting the understanding of RRBI physiology, particularly their associated neuroanatomical structures. The complex RRBI heterogeneity needs to explore the whole RRBI spectrum by integrating the clinical context [autistic individuals, their relatives and typical developing (TD) individuals]. We hypothesized that different RRBI dimensions would emerge by exploring the whole spectrum of RRBI and that these dimensions are associated with neuroanatomical signatures-involving cortical and subcortical areas.

METHOD

A sample of 792 individuals composed of 267 autistic subjects, their 370 first-degree relatives and 155 TD individuals was enrolled in the study. We assessed the whole patterns of RRBI in each individual by using the Repetitive Behavior Scale-Revised and the Yale-Brown Obsessive Compulsive Scale. We estimated brain volumes using MRI scanner for a subsample of the subjects (n = 152, 42 ASD, 89 relatives and 13 TD). We first investigated the dimensionality of RRBI by performing a principal component analysis on all items of these scales and included all the sampling population. We then explored the relationship between RRBI-derived factors with brain volumes using linear regression models.

RESULTS

We identified 3 main factors (with 30.3% of the RRBI cumulative variance): Factor 1 (FA1, 12.7%) reflected mainly the 'motor-driven' RRBI symptoms; Factor 2 and 3 (respectively, 8.8% and 7.9%) gathered mainly Y-BOCS related items and represented the 'cognitively driven' RRBI symptoms. These three factors were significantly associated with the right/left putamen volumes but with opposite effects: FA1 was negatively associated with an increased volume of the right/left putamen conversely to FA2 and FA3 (all uncorrected p < 0.05). FA1 was negatively associated with the left amygdala (uncorrected p < 0.05), and FA2 was positively associated with the left parietal structure (uncorrected p = 0.001).

CONCLUSION

Our results suggested 3 coherent RRBI dimensions involving the putamen commonly and other structures according to the RRBI dimension. The exploration of the putamen's integrative role in RSBI needs to be strengthened in further studies.

Regional homogeneity in patients with obsessive-compulsive disorder and depression: A resting state functional magnetic resonance imaging study

OBJECTIVE

To explore the brain functional impairment of patients with obsessive-compulsive disorder (OCD) with and without depressive symptoms and analyze the correlation between the degree of impairment and the severity of symptoms.

METHOD

Fourteen patients with OCD who met the ICD-10 diagnostic criteria for OCD were included. The group having OCD with depression (OCDd) consisted of 15 patients, and 17 healthy controls (HC) matched for age and education were also included. The Yale-Brown OCD Scale (Y-BOCS) and the 24-item Hamilton Assessment of Depression Scale (HAMD) were administered to the OCD and OCDd groups. Resting-state functional brain magnetic resonance imaging was performed in the three groups of participants.

RESULT

The OCDd group had lower scores on the HAMD, Y-BOCS, and obsessive-compulsive thinking subscales compared with the OCD group (P < 0.05). The scores on the OCDd subscale were negatively correlated with the HAMD scores (R =  - 0.568, P = 0.027). The OCDd group had higher regional homogeneity (ReHo) values in the lingual gyrus than the OCD group. The OCDd group had higher ReHo values in the lingual gyrus than the HC group, and the OCDd group had higher ReHo values than the HC group. These differences were statistically significant (P < 0.05). After correction for multiple comparisons, significant difference was observed between the OCDd and HC groups (P＜0.05). In the OCD group, the ReHo value of the lingual gyrus was negatively correlated with the Y-BOCS total score and the compulsive behavior subscale score (R =  - 0.609, -0.552; P = 0.016, 0.033).

CONCLUSION

Abnormal ReHo values in the lingual gyrus and right medial superior frontal gyrus were found in the patients with OCDd. In the OCDd group, the ReHo values of the lingual gyrus were negatively correlated with the scores on the Y-BOCS total and obsessive-compulsive subscales, suggesting that abnormal local coherence of the lingual gyrus may be related to the severity of OCD.

Effects of Gene Dosage and Development on Subcortical Nuclei Volumes in Individuals with 22q11.2 Copy Number Variations

The 22q11.2 locus contains genes critical for brain development. Reciprocal Copy Number Variations (CNVs) at this locus impact risk for neurodevelopmental and psychiatric disorders. Both 22q11.2 deletions (22qDel) and duplications (22qDup) are associated with autism, but 22qDel uniquely elevates schizophrenia risk. Understanding brain phenotypes associated with these highly penetrant CNVs can provide insights into genetic pathways underlying neuropsychiatric disorders. Human neuroimaging and animal models indicate subcortical brain alterations in 22qDel, yet little is known about developmental differences across specific nuclei between reciprocal 22q11.2 CNV carriers and typically developing (TD) controls. We conducted a longitudinal MRI study in 22qDel (n=96, 53.1% female), 22qDup (n=37, 45.9% female), and TD controls (n=80, 51.2% female), across a wide age range (5.5-49.5 years). Volumes of the thalamus, hippocampus, amygdala, and anatomical subregions were estimated using FreeSurfer, and the effect of 22q11.2 gene dosage was examined using linear mixed models. Age-related changes were characterized with general additive mixed models (GAMMs). Positive gene dosage effects (22qDel < TD < 22qDup) were observed for total intracranial and whole hippocampus volumes, but not whole thalamus or amygdala volumes. Several amygdala subregions exhibited similar positive effects, with bi-directional effects found across thalamic nuclei. Distinct age-related trajectories were observed across the three groups. Notably, both 22qDel and 22qDup carriers exhibited flattened development of hippocampal CA2/3 subfields relative to TD controls. This study provides novel insights into the impact of 22q11.2 CNVs on subcortical brain structures and their developmental trajectories.

Rhesus monkeys exhibiting spontaneous ritualistic behaviors resembling obsessive-compulsive disorder

Obsessive-compulsive disorder (OCD) is a chronic and debilitating psychiatric disorder that affects ∼2%-3% of the population globally. Studying spontaneous OCD-like behaviors in non-human primates may improve our understanding of the disorder. In large rhesus monkey colonies, we found 10 monkeys spontaneously exhibiting persistent sequential motor behaviors (SMBs) in individual-specific sequences that were repetitive, time-consuming and stable over prolonged periods. Genetic analysis revealed severely damaging mutations in genes associated with OCD risk in humans. Brain imaging showed that monkeys with SMBs had larger gray matter (GM) volumes in the left caudate nucleus and lower fractional anisotropy of the corpus callosum. The GM volume of the left caudate nucleus correlated positively with the daily duration of SMBs. Notably, exposure to a stressor (human presence) significantly increased SMBs. In addition, fluoxetine, a serotonergic medication commonly used for OCD, decreased SMBs in these monkeys. These findings provide a novel foundation for developing better understanding and treatment of OCD.

Cortical and subcortical contributions to non-motor inhibitory control: an fMRI study

Inhibition is a core executive cognitive function. However, the neural correlates of non-motor inhibitory control are not well understood. We investigated this question using functional Magnetic Resonance Imaging (fMRI) and a simple Count Go/NoGo task (n = 23), and further explored the causal relationships between activated brain regions. We found that the Count NoGo task activated a distinct pattern in the subcortical basal ganglia, including bilateral ventral anterior/lateral nucleus of thalamus (VA/VL), globus pallidus/putamen (GP/putamen), and subthalamic nucleus (STN). Stepwise regressions and mediation analyses revealed that activations in these region(s) were modulated differently by only 3 cortical regions i.e. the right inferior frontal gyrus/insula (rIFG/insula), along with left IFG/insula, and anterior cingulate cortex/supplementary motor area (ACC/SMA). The activations of bilateral VA/VL were modulated by both rSTN and rIFG/insula (with rGP/putamen as a mediator) independently, and the activation of rGP/putamen was modulated by ACC/SMA, with rIFG/insula as a mediator. Our findings provide the neural correlates of inhibitory control of counting and causal relationships between them, and strongly suggest that both indirect and hyperdirect pathways of the basal ganglia are involved in the Count NoGo condition.

The functional connectome in obsessive-compulsive disorder: resting-state mega-analysis and machine learning classification for the ENIGMA-OCD consortium

Current knowledge about functional connectivity in obsessive-compulsive disorder (OCD) is based on small-scale studies, limiting the generalizability of results. Moreover, the majority of studies have focused only on predefined regions or functional networks rather than connectivity throughout the entire brain. Here, we investigated differences in resting-state functional connectivity between OCD patients and healthy controls (HC) using mega-analysis of data from 1024 OCD patients and 1028 HC from 28 independent samples of the ENIGMA-OCD consortium. We assessed group differences in whole-brain functional connectivity at both the regional and network level, and investigated whether functional connectivity could serve as biomarker to identify patient status at the individual level using machine learning analysis. The mega-analyses revealed widespread abnormalities in functional connectivity in OCD, with global hypo-connectivity (Cohen's d: -0.27 to -0.13) and few hyper-connections, mainly with the thalamus (Cohen's d: 0.19 to 0.22). Most hypo-connections were located within the sensorimotor network and no fronto-striatal abnormalities were found. Overall, classification performances were poor, with area-under-the-receiver-operating-characteristic curve (AUC) scores ranging between 0.567 and 0.673, with better classification for medicated (AUC = 0.702) than unmedicated (AUC = 0.608) patients versus healthy controls. These findings provide partial support for existing pathophysiological models of OCD and highlight the important role of the sensorimotor network in OCD. However, resting-state connectivity does not so far provide an accurate biomarker for identifying patients at the individual level.

Dynamic functional thalamocortical dysconnectivity in schizophrenia correlates to antipsychotics response

Although many studies have showed abnormal thalamocortical networks in patients with schizophrenia (SCZ), the dynamic functional thalamocortical connectivity of individuals with SCZ and the effect of antipsychotics on this connectivity have not been investigated. Drug-naïve first-episode individuals with SCZ and healthy controls were recruited. Patients were treated with risperidone for 12 weeks. Resting-state functional magnetic resonance imaging was acquired at baseline and week 12. We identified six functional thalamic subdivisions. The sliding window strategy was used to determine the dynamic functional connectivity (dFC) of each functional thalamic subdivision. Individuals with SCZ displayed decreased or increased dFC variance in different thalamic subdivisions. The baseline dFC between ventral posterior-lateral (VPL) portions and right dorsolateral superior frontal gyrus (rdSFG) correlated with psychotic symptoms. The dFC variance between VPL and right medial orbital superior frontal gyrus (rmoSFG) or rdSFG decreased after 12-week risperidone treatment. The decreased dFC variance between VPL and rmoSFG correlated with the reduction of PANSS scores. Interestingly, the dFC between VPL and rmoSFG or rdSFG decreased in responders. The dFC variance change of VPL and the averaged whole brain signal correlated with the risperidone efficacy. Our study demonstrates abnormal variability in thalamocortical dFC may be implicated in psychopathological symptoms and risperidone response in individuals with schizophrenia, suggesting that thalamocortical dFC variance may be correlated to the efficacy of antipsychotic treatment.Registration: ClinicalTrials.gov Identifier: NCT00435370. https://www.clinicaltrials.gov/ct2/show/NCT00435370?term=NCT00435370&draw=2&rank=1.

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.

Enhancing neuroimaging genetics through meta-analysis for Tourette syndrome (ENIGMA-TS): A worldwide platform for collaboration

Tourette syndrome (TS) is characterized by multiple motor and vocal tics, and high-comorbidity rates with other neuropsychiatric disorders. Obsessive compulsive disorder (OCD), attention deficit hyperactivity disorder (ADHD), autism spectrum disorders (ASDs), major depressive disorder (MDD), and anxiety disorders (AXDs) are among the most prevalent TS comorbidities. To date, studies on TS brain structure and function have been limited in size with efforts mostly fragmented. This leads to low-statistical power, discordant results due to differences in approaches, and hinders the ability to stratify patients according to clinical parameters and investigate comorbidity patterns. Here, we present the scientific premise, perspectives, and key goals that have motivated the establishment of the Enhancing Neuroimaging Genetics through Meta-Analysis for TS (ENIGMA-TS) working group. The ENIGMA-TS working group is an international collaborative effort bringing together a large network of investigators who aim to understand brain structure and function in TS and dissect the underlying neurobiology that leads to observed comorbidity patterns and clinical heterogeneity. Previously collected TS neuroimaging data will be analyzed jointly and integrated with TS genomic data, as well as equivalently large and already existing studies of highly comorbid OCD, ADHD, ASD, MDD, and AXD. Our work highlights the power of collaborative efforts and transdiagnostic approaches, and points to the existence of different TS subtypes. ENIGMA-TS will offer large-scale, high-powered studies that will lead to important insights toward understanding brain structure and function and genetic effects in TS and related disorders, and the identification of biomarkers that could help inform improved clinical practice.

Obsessive-compulsive symptoms and resting-state functional characteristics in pre-adolescent children from the general population

While functional brain characteristics of obsessive-compulsive disorder have been extensively studied, literature on network topology and subnetwork connectivity related to obsessive-compulsive symptoms (OCS) is sparse. Here we investigated the functional brain characteristics of OCS in children from the general population using a multiscale approach. Since we previously observed OCS-related differences in thalamus morphology, we also focused on the network participation of thalamic subregions. The study included 1701 participants (9-12 years) from the population-based Generation R study. OCS were measured using the Short Obsessive-Compulsive Disorder Screener. We studied the brain network at multiple scales: global network topology, subnetwork connectivity and network participation of thalamic nodes (pre-registration: https://osf.io/azr9c ). Modularity, small-worldness and average participation coefficient were calculated on the global scale. We used a data-driven consensus community approach to extract a partition of five subnetworks involving thalamic subregions and calculate the within- and between-subnetwork functional connectivity and topology. Multiple linear regression models were fitted to model the relationship between OCS and functional brain measures. No significant associations were found when using our preregistered definition of probable OCS. However, post-hoc analyses showed that children endorsing at least one OCS (compared with controls) had higher modularity, lower connectivity between frontoparietal, limbic and visual networks as well as altered participation of the lateral prefrontal thalamus node. Our results suggest that network characteristics of OCS in children from the general population are partly symptom-specific and severity-dependent. Thorough assessment of symptom dimensions can deepen our understanding of OCS-related brain networks.