Previous research on the changes in resting-state functional connectivity (rsFC) in anorexia nervosa (AN) has been limited by an insufficient sample size, which reduced the reliability of the results and made it difficult to set the whole brain as regions of interest (ROIs).

We analyzed functional magnetic resonance imaging data from 114 female AN patients and 135 healthy controls (HC) and obtained self-reported psychological scales, including eating disorder examination questionnaire 6.0. One hundred sixty-four cortical, subcortical, cerebellar, and network parcellation regions were considered as ROIs. We calculated the ROI-to-ROI rsFCs and performed group comparisons.

Compared to HC, AN patients showed 12 stronger rsFCs mainly in regions containing dorsolateral prefrontal cortex (DLPFC), and 33 weaker rsFCs primarily in regions containing cerebellum, within temporal lobe, between posterior fusiform cortex and lateral part of visual network, and between anterior cingulate cortex (ACC) and thalamus (p < 0.01, false discovery rate [FDR] correction). Comparisons between AN subtypes showed that there were stronger rsFCs between right lingual gyrus and right supracalcarine cortex and between left temporal occipital fusiform cortex and medial part of visual network in the restricting type compared to the binge/purging type (p < 0.01, FDR correction).

Stronger rsFCs in regions containing mainly DLPFC, and weaker rsFCs in regions containing primarily cerebellum, within temporal lobe, between posterior fusiform cortex and lateral part of visual network, and between ACC and thalamus, may represent categorical diagnostic markers discriminating AN patients from HC.

Less than a third of patients with depression achieve successful remission with standard first-step antidepressant monotherapy. The process for determining appropriate second-step care is often based on clinical intuition and involves a protracted course of trial and error, resulting in substantial patient burden and unnecessary delay in the provision of optimal treatment. To address this problem, we adopt an ensemble machine learning approach to improve prediction accuracy of remission in response to second-step treatments.

Data were derived from the Level 2 stage of the STAR*D dataset, which included 1439 patients who were randomized into one of seven different second-step treatment strategies after failing to achieve remission during first-step antidepressant treatment. Ensemble machine learning models, comprising several individual algorithms, were evaluated using nested cross-validation on 155 predictor variables including clinical and demographic measures.

The ensemble machine learning algorithms exhibited differential classification performance in predicting remission status across the seven second-step treatments. For the full set of predictors, AUC values ranged from 0.51 to 0.82 depending on the second-step treatment type. Predicting remission was most successful for cognitive therapy (AUC = 0.82) and least successful for other medication and combined treatment options (AUCs = 0.51–0.66).

Ensemble machine learning has potential to predict second-step treatment. In this study, predictive performance varied by type of treatment, with greater accuracy in predicting remission in response to behavioral treatments than to pharmacotherapy interventions. Future directions include considering more informative predictor modalities to enhance prediction of second-step treatment response.

Dysmyelination could be part of the pathophysiology of schizophrenia spectrum (SCZ) and bipolar disorders (BPD), yet few studies have examined myelination of the cerebral cortex. The ratio of T1- and T2-weighted magnetic resonance images (MRI) correlates with intracortical myelin. We investigated the T1w/T2w-ratio and its age trajectories in patients and healthy controls (CTR) and explored associations with antipsychotic medication use and psychotic symptoms.

Patients with SCZ (n = 64; mean age = 30.4 years, s.d. = 9.8), BPD (n = 91; mean age 31.0 years, s.d. = 10.2), and CTR (n = 155; mean age = 31.9 years, s.d. = 9.1) who participated in the TOP study (NORMENT, University of Oslo, Norway) were clinically assessed and scanned using a General Electric 3 T MRI system. T1w/T2w-ratio images were computed using an optimized pipeline with intensity normalization and field inhomogeneity correction. Vertex-wise regression models were used to compare groups and examine group × age interactions. In regions showing significant differences, we explored associations with antipsychotic medication use and psychotic symptoms.

No main effect of diagnosis was found. However, age slopes of the T1w/T2w-ratio differed significantly between SCZ and CTR, predominantly in frontal and temporal lobe regions: Lower T1w/T2w-ratio values with higher age were found in CTR, but not in SCZ. Follow-up analyses revealed a more positive age slope in patients who were using antipsychotics and patients using higher chlorpromazine-equivalent doses.

While we found no evidence of reduced intracortical myelin in SCZ or BPD relative to CTR, different regional age trajectories in SCZ may suggest a promyelinating effect of antipsychotic medication.

Comorbid depression substantially affects the management of glycemia and diabetes-related complications among patients with type 2 diabetes mellitus. In this study, we sought to determine the association between weight change over 4 years and depression risk among patients with type 2 diabetes mellitus.

This population-based retrospective cohort study from the National Health Insurance Services of Korea included 1 111 345 patients with type 2 diabetes who were divided into groups according to body weight change over 4 years. Body weight changes were compared with the preceding 4-year period (2005–2008). Depression was defined according to the International Classification of Diseases 10th revision code for depression (F32 and F33) on one or more inpatient or outpatient claims.

During a median follow-up of 7.4 years, 244 081 cases of depression were identified. We observed a U-shaped association between body weight change and depression risk with a higher risk among both groups of weight loss (hazard ratio (HR) 1.17, 95% CI 1.15–1.19 for ⩾ −10%; HR 1.07, 95% CI 1.06–1.08 for −10 to −5%) and weight gain (HR 1.06, 95% CI 1.04–1.08 for ⩾10%; HR 1.02, 95% CI 1.01–1.04 for 5–10%) compared with the stable weight group (−5 to 5%).

A U-shaped association between body weight change and depression risk was observed in this large nationwide cohort study. Our study suggests that patients with type 2 diabetes and weight change, either gain or loss, could be considered a high-risk group for depression.

Major depressive disorder (MDD) is highly prevalent and burdensome for individuals and society. While there are psychological interventions able to prevent and treat MDD, uptake remains low. To overcome structural and attitudinal barriers, an indirect approach of using online insomnia interventions seems promising because insomnia is less stigmatized, predicts MDD onset, is often comorbid and can outlast MDD treatment. This individual-participant-data meta-analysis evaluated the potential of the online insomnia intervention GET.ON Recovery as an indirect treatment to reduce depressive symptom severity (DSS) and potential MDD onset across a range of participant characteristics.

Efficacy on depressive symptom outcomes was evaluated using multilevel regression models controlling for baseline severity. To identify potential effect moderators, clinical, sociodemographic, and work-related variables were investigated using univariable moderation and random-forest methodology before developing a multivariable decision tree.

IPD were obtained from four of seven eligible studies (N = 561); concentrating on workers with high work-stress. DSS was significantly lower in the intervention group both at post-assessment (d = −0.71 [95% CI−0.92 to −0.51]) and at follow-up (d = −0.84 [95% CI −1.11 to −0.57]). In the subsample (n = 121) without potential MDD at baseline, there were no significant group differences in onset of potential MDD. Moderation analyses revealed that effects on DSS differed significantly across baseline severity groups with effect sizes between d = −0.48 and −0.87 (post) and d = − 0.66 to −0.99 (follow-up), while no other sociodemographic, clinical, or work-related characteristics were significant moderators.

An online insomnia intervention is a promising approach to effectively reduce DSS in a preventive and treatment setting.

Growing evidence points to the pivotal role of vitamin D in the pathophysiology and treatment of major depressive disorder (MDD). However, there is a paucity of longitudinal research investigating the effects of vitamin D supplementation on the brain of MDD patients.

We conducted a double-blind randomized controlled trial in 46 MDD patients, who were randomly allocated into either VD (antidepressant medication + vitamin D supplementation) or NVD (antidepressant medication + placebos) groups. Data from diffusion tensor imaging, resting-state functional MRI, serum vitamin D concentration, and clinical symptoms were obtained at baseline and after an average of 7 months of intervention.

Both VD and NVD groups showed significant improvement in depression and anxiety symptoms but with no significant differences between the two groups. However, a greater increase in serum vitamin D concentration was found to be associated with greater improvement in depression and anxiety symptoms in VD group. More importantly, neuroimaging data demonstrated disrupted white matter integrity of right inferior fronto-occipital fasciculus along with decreased functional connectivity between right frontoparietal and medial visual networks after intervention in NVD group, but no changes in VD group.

These findings suggest that vitamin D supplementation as adjunctive therapy to antidepressants may not only contribute to improvement in clinical symptoms but also help preserve brain structural and functional connectivity in MDD patients.

The relationship between childhood trauma (CT) and psychotic symptoms in patients with schizophrenia (SCZ), and subthreshold psychotic experiences in non-clinical populations is well-established. However, little is known about the relationship between subtypes of trauma and specific symptoms in patients, their siblings, and controls. It is also not clear which variables mediate the relationship between trauma and psychotic symptoms.

Seven hundred and forty-two patients with SCZ, 718 of their unaffected siblings and 1039 controls from three EU-GEI sites were assessed for CT, symptom severity, and cognitive schemas about self/others. CT was assessed with the Childhood Trauma Questionnaire, and cognitive schemas were assessed by The Brief Core Schema Scale.

Patients with psychosis were affected by CT more than their siblings and controls in all domains. Childhood emotional abuse and neglect were more common in siblings than controls. CT was related to negative cognitive schemas toward self/others in patients, siblings, and controls. We found that negative schemas about self-mediated the relationship between emotional abuse and thought withdrawal and thought broadcasting. Approximately 33.9% of the variance in these symptoms was explained by the mediator. It also mediated the relationship between sexual abuse and persecutory delusions in SCZ.

Our findings suggest that childhood abuse and neglect are more common in patients with schizophrenia than their siblings and healthy controls, and have different impacts on clinical domains which we searched. The relationship between CT and positive symptoms seems to be mediated by negative cognitive schemas about self in schizophrenia.

Cardiac vagal tone is an indicator of parasympathetic nervous system functioning, and there is increasing interest in its relation to antisocial behavior. It is unclear however whether antisocial individuals are characterized by increased or decreased vagal tone, and whether increased vagal tone is the source of the low heart rate frequently reported in antisocial populations.

Participants consisted of four groups of community-dwelling adolescent boys aged 15.7 years: (1) controls, (2) childhood-only antisocial, (3) adolescent-only antisocial, and (4) persistently antisocial. Heart rate and vagal tone were assessed in three different conditions: rest, cognitive stressor, and social stressor.

All three antisocial groups had both lower resting heart rates and increased vagal tone compared to the low antisocial controls across all three conditions. Low heart rate partially mediated the relationship between vagal tone and antisocial behavior.

Results indicate that increased vagal tone and reduced heart rate are relatively broad risk factors for different developmental forms of antisocial behavior. Findings are the first to implicate vagal tone as an explanatory factor in understanding heart rate – antisocial behavior relationships. Future experimental work using non-invasive vagus nerve stimulation or heart rate variability biofeedback is needed to more systematically evaluate this conclusion.

First-episode psychotic disorders comprise a heterogeneous phenotype with a complex etiology involving numerous common small-effect genetic variations and a wide range of environmental exposures. We examined whether a family of schizophrenia spectrum disorder (FH-Sz) interacts with an environmental risk score (ERS-Sz) regarding the outcome of patients with non-affective first episode psychosis (NAFEP).

We included 288 patients with NAFEP who were evaluated after discharge from an intensive 2-year program. We evaluated three outcome measures: symptomatic remission, psychosocial functioning, and personal recovery. We analyzed the main and joint associations of a FH-Sz and the ERS-Sz on the outcomes by using the relative excess risk due to interaction (RERI) approach.

A FH-Sz showed a significant association with poor symptomatic remission and psychosocial functioning outcomes, although there was no significant interaction between a FH-Sz and the ERS-Sz on these outcomes. The ERS-Sz did not show a significant association with poor symptomatic remission and psychosocial functioning outcomes, even though the magnitude of the interaction between ERS-Sz and FH-Sz with the later outcome was moderate (RERI = 6.89, 95% confidence interval −16.03 to 29.81). There was no association between a FH-Sz and the ERS-Sz and personal recovery.

Our results provide further empirical support regarding the contribution of FH-Sz to poor symptomatic remission and poor psychosocial functioning outcomes in patients with NAFEP.

Class and social disadvantage have long been identified as significant factors in the etiology and epidemiology of psychosis. Few studies have explicitly examined the impact of intersecting social disadvantage on long-term employment and financial independence.

We applied latent class analysis (LCA) to 20-year longitudinal data from participants with affective and non-affective psychosis (n = 256) within the Chicago Longitudinal Research. LCA groups were modeled using multiple indicators of pre-morbid disadvantage (parental social class, educational attainment, race, gender, and work and social functioning prior to psychosis onset). The comparative longitudinal work and financial functioning of LCA groups were then examined.

We identified three distinct latent classes: one comprised entirely of White participants, with the highest parental class and highest levels of educational attainment; a second predominantly working-class group, with equal numbers of Black and White participants; and a third with the lowest parental social class, lowest levels of education and a mix of Black and White participants. The latter, our highest social disadvantage group experienced significantly poorer employment and financial outcomes at all time-points, controlling for diagnosis, symptoms, and hospitalizations prior to baseline. Contrary to our hypotheses, on most measures, the two less disadvantaged groups did not significantly differ from each other.

Our analyses add to a growing literature on the impact of multiple forms of social disadvantage on long-term functional trajectories, underscoring the importance of proactive attention to sociostructural disadvantage early in treatment, and the development and evaluation of interventions designed to mitigate ongoing social stratification.

Post-traumatic stress symptoms (PTSS) were the most frequently reported mental health concern for youth during COVID-19, yet variations in youth's PTSS responses warrant empirical consideration. Features of the caregiving environment influence youth's responses to environmental stressors, and youth's parasympathetic nervous system regulation may qualify the magnitude and/or direction of these effects. This prospective investigation evaluated diathesis stress and differential susceptibility models of caregiving and parasympathetic influences on youth's PTSS responses to COVID-19.

Participants were 225 caregiver-youth dyads (youth 49.8% female at birth; 88.4% non-white) followed from childhood through adolescence and COVID-19. Youth's resting respiratory sinus arrhythmia (RSA; Mage = 6.11, s.d. = 0.21), caregiving features (i.e. attachment security [youth Mage = 12.24, s.d. = 0.35] and caregiver internalizing psychopathology [caregiver Mage = 39.29, s.d. = 6.78]) were assessed pre-pandemic. Youth's PTSS was assessed one year prior to the US COVID-19 pandemic (Mage = 14.24, s.d. = 0.50) and during the spring of 2020 at the height of the pandemic (Mage = 15.23, s.d. = 0.57).

Youth's PTSS increased during COVID-19. Youth with relatively high resting RSA evidenced the lowest PTSS when their caregiving environment featured high attachment security or low caregiver internalizing problems, but the highest PTSS when their caregiving environment featured low attachment security or high caregiver internalizing problems. In contrast, PTSS levels of youth with relatively low or average resting RSA did not differ significantly depending on attachment security or caregiver internalizing.

Results are consistent with a differential susceptibility hypothesis, wherein relatively high resting RSA conferred heightened sensitivity to caregiving environments in a for-better-and-for-worse manner during COVID-19.

Attention-deficit hyperactivity disorder (ADHD) is often comorbid with other medical conditions in adult patients. However, ADHD is extremely underdiagnosed in adults and little is known about the medical comorbidities in undiagnosed adult individuals with high ADHD liability. In this study we investigated associations between ADHD genetic liability and electronic health record (EHR)-based ICD-10 diagnoses across all diagnostic categories, in individuals without ADHD diagnosis history.

We used data from the Estonian Biobank cohort (N = 111 261) and generated polygenic risk scores (PRS) for ADHD (PRSADHD) based on the ADHD genome-wide association study. We performed a phenome-wide association study (PheWAS) to test for associations between standardized PRSADHD and 1515 EHR-based ICD-10 diagnoses in the full and sex-stratified sample. We compared the observed significant ICD-10 associations to associations with (1) ADHD diagnosis and (2) questionnaire-based high ADHD risk analyses.

After Bonferroni correction (p = 3.3 × 10−5) we identified 80 medical conditions associated with PRSADHD. The strongest evidence was seen with chronic obstructive pulmonary disease (OR 1.15, CI 1.11–1.18), obesity (OR 1.13, CI 1.11–1.15), and type 2 diabetes (OR 1.11, CI 1.09–1.14). Sex-stratified analysis generally showed similar associations in males and females. Out of all identified associations, 40% and 78% were also observed using ADHD diagnosis or questionnaire-based ADHD, respectively, as the predictor.

Overall our findings indicate that ADHD genetic liability is associated with an increased risk of a substantial number of medical conditions in undiagnosed individuals. These results highlight the need for timely detection and improved management of ADHD symptoms in adults.

Cerebral microvascular dysfunction may contribute to depression via disruption of brain structures involved in mood regulation, but evidence is limited. We investigated the association of retinal microvascular function, a proxy for microvascular function in the brain, with incidence and trajectories of clinically relevant depressive symptoms.

Longitudinal data are from The Maastricht Study of 5952 participants (59.9 ± 8.5 years/49.7% women) without clinically relevant depressive symptoms at baseline (2010–2017). Central retinal arteriolar equivalent and central retinal venular equivalent (CRAE and CRVE) and a composite score of flicker light-induced retinal arteriolar and venular dilation were assessed at baseline. We assessed incidence and trajectories of clinically relevant depressive symptoms (9-item Patient Health Questionnaire score ⩾10). Trajectories included continuously low prevalence (low, n = 5225 [87.8%]); early increasing, then chronic high prevalence (early-chronic, n = 157 [2.6%]); low, then increasing prevalence (late-increasing, n = 247 [4.2%]); and remitting prevalence (remitting, n = 323 [5.4%]).

After a median follow-up of 7.0 years (range 1.0–11.0), 806 (13.5%) individuals had incident clinically relevant depressive symptoms. After full adjustment, a larger CRAE and CRVE were each associated with a lower risk of clinically relevant depressive symptoms (hazard ratios [HRs] per standard deviation [s.d.]: 0.89 [95% confidence interval (CI) 0.83–0.96] and 0.93 [0.86–0.99], respectively), while a lower flicker light-induced retinal dilation was associated with a higher risk of clinically relevant depressive symptoms (HR per s.d.: 1.10 [1.01–1.20]). Compared to the low trajectory, a larger CRAE was associated with lower odds of belonging to the early-chronic trajectory (OR: 0.83 [0.69–0.99]) and a lower flicker light-induced retinal dilation was associated with higher odds of belonging to the remitting trajectory (OR: 1.23 [1.07–1.43]).

These findings support the hypothesis that cerebral microvascular dysfunction contributes to the development of depressive symptoms.

The nature of the pathway from conduct disorder (CD) in adolescence to antisocial behavior in adulthood has been debated and the role of certain mediators remains unclear. One perspective is that CD forms part of a general psychopathology dimension, playing a central role in the developmental trajectory. Impairment in reflective functioning (RF), i.e., the capacity to understand one's own and others' mental states, may relate to CD, psychopathology, and aggression. Here, we characterized the structure of psychopathology in adult male-offenders and its role, along with RF, in mediating the relationship between CD in their adolescence and current aggression.

A secondary analysis of pre-treatment data from 313 probation-supervised offenders was conducted, and measures of CD symptoms, general and specific psychopathology factors, RF, and aggression were evaluated through clinical interviews and questionnaires.

Confirmatory factor analyses indicated that a bifactor model best fitted the sample's psychopathology structure, including a general psychopathology factor (p factor) and five specific factors: internalizing, disinhibition, detachment, antagonism, and psychoticism. The structure of RF was fitted to the data using a one-factor model. According to our mediation model, CD significantly predicted the p factor, which was positively linked to RF impairments, resulting in increased aggression.

These findings highlight the critical role of a transdiagnostic approach provided by RF and general psychopathology in explaining the link between CD and aggression. Furthermore, they underscore the potential utility of treatments focusing on RF, such as mentalization-based treatment, in mitigating aggression in offenders with diverse psychopathologies.

Emotional problems, especially anxiety, have become increasingly common in recent generations. Few population-based studies have examined trajectories of emotional problems from early childhood to late adolescence or investigated differences in psychiatric and functional outcomes.

Using the Avon Longitudinal Study of Parents and Children (ALSPAC, n = 8286, 50.4% male), we modeled latent class growth trajectories of emotional problems, using the parent-reported Strength and Difficulties Questionnaire emotional scale (SDQ-E) on seven occasions (4–17 years). Psychiatric outcomes in young adulthood (21–25 years) were major depressive disorder (MDD), generalized anxiety disorder (GAD), and self-harm. Functional outcomes were exam attainment, educational/occupational status, and social relationship quality.

We identified four classes of emotional problems: low (67.0%), decreasing (18.4%), increasing (8.9%), and persistent (5.7%) problems. Compared to those in the low class, individuals with decreasing emotional problems were not at elevated risk of any poor adult outcome. Individuals in the increasing and persistent classes had a greater risk of adult MDD (RR: 1.59 95% CI 1.13–2.26 and RR: 2.25 95% CI 1.49–3.41) and self-harm (RR: 2.37 95% CI 1.91–2.94 and RR: 1.87 95% CI 1.41–2.48), and of impairment in functional domains. Childhood sleep difficulties, irritability, conduct and neurodevelopmental problems, and family adversity were associated with a persistent course of emotional problems.

Childhood emotional problems were common, but those whose symptoms improved over time were not at increased risk for adverse adult outcomes. In contrast, individuals with persistent or adolescent-increasing emotional problems had a higher risk of mental ill-health and social impairment in young adulthood which was especially pronounced for those with persistent emotional problems.

The brain can be represented as a network, with nodes as brain regions and edges as region-to-region connections. Nodes with the most connections (hubs) are central to efficient brain function. Current findings on structural differences in Major Depressive Disorder (MDD) identified using network approaches remain inconsistent, potentially due to small sample sizes. It is still uncertain at what level of the connectome hierarchy differences may exist, and whether they are concentrated in hubs, disrupting fundamental brain connectivity.

We utilized two large cohorts, UK Biobank (UKB, N = 5104) and Generation Scotland (GS, N = 725), to investigate MDD case–control differences in brain network properties. Network analysis was done across four hierarchical levels: (1) global, (2) tier (nodes grouped into four tiers based on degree) and rich club (between-hub connections), (3) nodal, and (4) connection.

In UKB, reductions in network efficiency were observed in MDD cases globally (d = −0.076, pFDR = 0.033), across all tiers (d = −0.069 to −0.079, pFDR = 0.020), and in hubs (d = −0.080 to −0.113, pFDR = 0.013–0.035). No differences in rich club organization and region-to-region connections were identified. The effect sizes and direction for these associations were generally consistent in GS, albeit not significant in our lower-N replication sample.

Our results suggest that the brain's fundamental rich club structure is similar in MDD cases and controls, but subtle topological differences exist across the brain. Consistent with recent large-scale neuroimaging findings, our findings offer a connectomic perspective on a similar scale and support the idea that minimal differences exist between MDD cases and controls.

Although dopaminergic disturbances are well-known in schizophrenia, the understanding of dopamine-related brain dynamics remains limited. This study investigates the dynamic coactivation patterns (CAPs) associated with the substantia nigra (SN), a key dopaminergic nucleus, in first-episode treatment-naïve patients with schizophrenia (FES).

Resting-state fMRI data were collected from 84 FES and 94 healthy controls (HCs). Frame-wise clustering was implemented to generate CAPs related to SN activation or deactivation. Connectome features of each CAP were derived using an edge-centric method. The occurrence for each CAP and the balance ratio for antagonistic CAPs were calculated and compared between two groups, and correlations between temporal dynamic metrics and symptom burdens were explored.

Functional reconfigurations in CAPs exhibited significant differences between the activation and deactivation states of SN. During SN activation, FES more frequently recruited a CAP characterized by activated default network, language network, control network, and the caudate, compared to HCs (F = 8.54, FDR-p = 0.030). Moreover, FES displayed a tilted balance towards a CAP featuring SN-coactivation with the control network, caudate, and thalamus, as opposed to its antagonistic CAP (F = 7.48, FDR-p = 0.030). During SN deactivation, FES exhibited increased recruitment of a CAP with activated visual and dorsal attention networks but decreased recruitment of its opposing CAP (F = 6.58, FDR-p = 0.034).

Our results suggest that neuroregulatory dysfunction in dopaminergic pathways involving SN potentially mediates aberrant time-varying functional reorganizations in schizophrenia. This finding enriches the dopamine hypothesis of schizophrenia from the perspective of brain dynamics.

Epigenetic changes are plausible molecular sources of clinical heterogeneity in schizophrenia. A subgroup of schizophrenia patients with elevated inflammatory or immune-dysregulation has been reported by previous studies. However, little is known about epigenetic changes in genes related to immune activation in never-treated first-episode patients with schizophrenia (FES) and its consistency with that in treated long-term ill (LTS) patients.

In this study, epigenome-wide profiling with a DNA methylation array was applied using blood samples of both FES and LTS patients, as well as their corresponding healthy controls. Non-negative matrix factorization (NMF) and k -means clustering were performed to parse heterogeneity of schizophrenia, and the consistency of subtyping results from two cohorts. was tested.

This study identified a subtype of patients in FES participants (47.5%) that exhibited widespread methylation level alterations of genes enriched in immune cell activity and a significantly higher proportion of neutrophils. This clustering of FES patients was validated in LTS patients, with high correspondence in epigenetic and clinical features across two cohorts

In summary, this study demonstrated a subtype of schizophrenia patients across both FES and LTS cohorts, defined by widespread alterations in methylation profile of genes related to immune function and distinguishing clinical features. This finding illustrates the promise of novel treatment strategies targeting immune dysregulation for a subpopulation of schizophrenia patients.