Hostname: page-component-7f64f4797f-42qgm Total loading time: 0 Render date: 2025-11-08T09:03:15.424Z Has data issue: false hasContentIssue false

Mendelian randomisation reveals modifiable pathways and epigenetic markers from childhood maltreatment to neuropsychiatric disorders

Published online by Cambridge University Press:  04 November 2025

Nicole Ying Ting Ng
Affiliation:
Department of Paediatrics and Adolescent Medicine, School of Clinical Medicine, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
John Yen Tang
Affiliation:
Department of Paediatrics and Adolescent Medicine, School of Clinical Medicine, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
Jie V. Zhao
Affiliation:
School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
Christopher Chun Yu Mak
Affiliation:
Department of Paediatrics and Adolescent Medicine, School of Clinical Medicine, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
Brian Hon Yin Chung*
Affiliation:
Department of Paediatrics and Adolescent Medicine, School of Clinical Medicine, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
*
Correspondence: Brian H. Y. Chung. Email: bhychung@hku.hk.
Rights & Permissions [Opens in a new window]

Abstract

Background

Observational studies suggested an association between childhood maltreatment and neuropsychiatric disorders; however, mediators remain disputed.

Aims

We aimed to confirm the relationship between childhood maltreatmen and neuropsychiatric disorders, and to identify addiction-related, biological, behavioural, cognitive, socioeconomic and epigenetic mediators.

Method

We used two-sample Mendelian randomisation and publicly available genome-wide association data to evaluate the effect of genetically predicted childhood maltreatment (N = 143 473) on the risk of six neuropsychiatric disorders (up to N = 500 199). We used two-step Mendelian randomisation to determine the proportion of the effect of childhood maltreatment on disorders that was mediated by mediators. We used multivariable Mendelian randomisation to determine the direct effect of childhood maltreatment on disorders accounting for mediators. We used epigenetic Mendelian randomisation to determine the effect of DNA methylation at childhood maltreatment-associated CpG sites on disorders.

Results

Childhood maltreatment was significantly associated with higher risk of attention-deficit/hyperactivity disorder (ADHD) (odds ratio 10.09, 95% CI: 4.76–21.40), major depressive disorder (MDD) (odds ratio 1.89, 95% CI: 1.32–2.70) and schizophrenia (odds ratio: 5.89, 95% CI: 1.46–23.78). We determined that 4.14–22.17% of the effect of childhood maltreatment was mediated by addiction-related behaviours (smoking initiation, leisure screen time and substance abuse), cognitive traits (executive functioning, intelligence and risk tolerance) and educational attainment. We found that the direct effects of childhood maltreatment on ADHD (odds ratio 2.57) and schizophrenia (odds ratio 5.10) were less than the total effects, while the direct effect on MDD (odds ratio 1.95) remained relatively unchanged. We found altered DNA methylation levels at 3, 4 and 19 CpG sites to be significantly associated with ADHD, MDD and schizophrenia, respectively.

Conclusions

These results emphasise the need for preventative strategies to reduce childhood maltreatment prevalence, including strengthening support for high-risk families and responsive strategies to mitigate consequences for victims, with clinical screening for childhood maltreatment history and holistic approaches addressing addiction-related, cognitive and socioeconomic mediators.

Information

Type
Original Article
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2025. Published by Cambridge University Press on behalf of Royal College of Psychiatrists

Childhood maltreatment involves physical, sexual and emotional abuse, and emotional and physical neglect, to people under 18 years of age within a relationship of responsibility, leading to harm in their health, survival and development. Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1,Reference Baumeister, Akhtar, Ciufolini, Pariante and Mondelli2 Before the COVID-19 pandemic, over 1 billion children experienced childhood maltreatment globally. Reference Park and Walsh3 During the pandemic, this number increased globally despite Target 16.2 of the United Nations’ Sustainable Development Goals (SDGs): ending all forms of abuse against children. Reference Lee and Kim4

Given the high current prevalence of childhood maltreatment, the long-term psychological consequences cannot be ignored. Reference Park and Walsh3 In addition to cross-sectional and longitudinal studies, a recent systematic review found significantly higher risks of depression, anxiety, attention-deficit/hyperactivity disorder (ADHD) and internalising and externalising problems in childhood maltreatment victims. Reference Baldwin, Wang, Karwatowska, Schoeler, Tsaligopoulou and Munafò5Reference Strathearn, Giannotti, Mills, Kisely, Najman and Abajobir10 However, pre-existing environmental and genetic risk factors for psychopathology may confound these findings. As such, previous evidence needs to be triangulated with other methods to confirm causality. Reference Baldwin, Wang, Karwatowska, Schoeler, Tsaligopoulou and Munafò5

Besides confirming causality, identification of causal mechanisms prioritises areas for intervention. Observational studies showed independent associations among stress, chronic inflammation, socioeconomic status, maladaptive behaviours and cognitive and emotional dysregulation with childhood maltreatment and neuropsychiatric outcomes. Reference Carr, Duff and Craddock6,Reference Macpherson, Gray, Ip, McCallum, Hanlon and Welsh7 However, while some studies supported the mediating role of behavioural, social and cognitive traits, Reference Alameda, Rodriguez, Carr, Aas, Trotta and Marino11Reference Keyes, Eaton, Krueger, McLaughlin, Wall and Grant13 others found inconsistent or weak mediating effects Reference Macpherson, Gray, Ip, McCallum, Hanlon and Welsh7,Reference Li, Luyten and Midgley14,Reference Kent, Hopfer, Corley and Stallings15 due to inaccurately measured mediators and inadequate power. Besides, unmeasured confounding and the inability to temporally order investigated traits limit the delineation of causal mediators. Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 These justify further research into behavioural, social and cognitive traits linking childhood maltreatment to neuropsychiatric outcomes.

The role of epigenetics such as DNA methylation (DNAm) has also been of recent interest. Reference Yang, Zhang, Ge, Weder, Douglas-Palumberi and Perepletchikova17 DNAm at genes involved in stress response, serotonin signalling, inflammation and cell growth were previously associated with childhood maltreatment in epigenome-wide association studies (EWAS) of childhood maltreatment. Reference Parade, Huffhines, Daniels, Stroud, Nugent and Tyrka18 However, current research on the mediating effects of DNAm is limited to candidate genes. Thus, the relationship among DNAm, childhood maltreatment and neuropsychiatric disorders remains unclear. Reference Smearman, Almli, Conneely, Brody, Sales and Bradley19,Reference Comtois-Cabana, Barr, Provençal, Ouellet-Morin and Provenzi20 This warrants further exploration of the effect of DNAm at other genes on the risk of neuropsychiatric disorders.

Mendelian randomisation uses genetic variants to proxy the effects of an exposure on an outcome. Reference Bowden, Del Greco, Minelli, Davey Smith, Sheehan and Thompson21 Because genetic variants are randomly allocated at conception, Mendelian randomisation overcomes the limitations of observational studies to support their findings. Recent Mendelian randomisation studies found that childhood maltreatment is causally associated with ADHD, major depressive disorder (MDD) and schizophrenia. Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1,Reference Alkema, Marchi, van der Zaag, van der Sluis, Warrier and Ophoff22,Reference Baltramonaityte, Pingault, Cecil, Choudhary, Järvelin and Penninx23 However, two of these studies utilised outcome data-sets with samples from the UK Biobank, which overlap with the exposure data, Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1,Reference Alkema, Marchi, van der Zaag, van der Sluis, Warrier and Ophoff22 thus introducing sample overlap bias; in addition, one investigated only a composite psychocardiometabolic outcome trait. Reference Baltramonaityte, Pingault, Cecil, Choudhary, Järvelin and Penninx23 Furthermore, none of the existing studies have conducted mediation Mendelian randomisation to delineate causal mediators of this relationship.

Aims

This study uses a mediation Mendelian randomisation study design to first confirm the causal association between childhood maltreatment and neuropsychiatric outcomes, then to evaluate the mediating role of addiction-related, biological, behavioural, cognitive and socioeconomic traits, and the potential role of DNAm.

Method

Study design

This study follows Strengthening the Reporting of Observational Studies in Epidemiology using Mendelian randomisation guidelines (Appendix 1 (found in the Supplementary material)). The Mendelian randomisation assumptions are described in Appendix 2. Figure 1 shows the overall workflow. Univariable Mendelian randomisation (UVMR) investigated the causal effect of childhood maltreatment on neuropsychiatric outcomes. Two-step Mendelian randomisation determined indirect effects acting through mediators. Multivariable Mendelian randomisation (MVMR) determined residual direct effects of childhood maltreatment, adjusting for significant mediators. Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 Epigenetic Mendelian randomisation determined the relationship between DNAm and outcomes. Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 All analyses utilised data from summary statistics of genome-wide association studies (GWAS) or epigenome-wide association studies (EWAS). Different sets of covariates were used for exposures and outcomes. Ethics approval and consent were obtained for each GWAS and EWAS, and can be found in the original publications.

Fig. 1 Overall study design with the aims, data-sets, statistical methods and expected outcomes. CM, childhood maltreatment; EWAS, epigenome-wide association study; GWAS, genome-wide association study; PGC, Psychiatric Genomics Consortium; ADHD, attention-deficit/hyperactivity disorder; MR, Mendelian randomisation; MR-PRESSO, Mendelian randomisation pleiotropy residual sum and outlier; MR-Egger, Mendelian randomisation-Egger; IL6, interleukin-6; IL6R, Interleukin-6 receptor; CRP, C-reactive protein; BMI, body mass index; IVW, inverse-variance weighted; MVMR, multivariable Mendelian randomisation; DNAm, DNA methylation; mQTL, methylation quantitative trait loci; GoDMC, Genetics of DNA Methylation Consortium.

Data sources

Exposure

Genetically predicted childhood maltreatment information was obtained from a GWAS of 143 473 participants of European ancestry in the UK Biobank. Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1 Childhood maltreatment was retrospectively self-reported through a five-item Childhood Trauma Screener, which rates each trauma subtype out of four, with zero indicating never experienced and four indicating often experienced. Scores for each subtype were summed to obtain a composite score for childhood maltreatment. The quantitative phenotypic trait in this GWAS is a log-transformed score of the Childhood Trauma Screener, with effect sizes representing a log-score increase per effect allele. Additional questionnaire and GWAS details are reported in the original study. Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1

Outcomes

Six neuropsychiatric outcomes were included: ADHD (N = 55 374), anxiety disorder (N = 21 761), autism spectrum disorder (N = 46 350), bipolar disorder (N = 413 466), MDD (N = 500 199) and schizophrenia (N = 275 249). Summary statistics from European-ancestry participants were obtained from the Psychiatric Genomics Consortium, and excluded UK Biobank data. Reference Sullivan, Agrawal, Bulik, Andreassen, Borglum and Breen24 Controls were participants without neuropsychiatric disorders or associated symptoms. Study details are shown in Supplementary Table 1 available at https://doi.org/10.1192/bjp.2025.10433.

Mediators

Each mediator was proxied by genetic instruments obtained from their respective GWAS summary statistics (see details in Supplementary Table 1). A total of 27 mediators were selected based on previously hypothesised pathways from observational studies. Reference Baldwin, Wang, Karwatowska, Schoeler, Tsaligopoulou and Munafò5,Reference Strathearn, Giannotti, Mills, Kisely, Najman and Abajobir10,Reference Scott, Malacova, Mathews, Haslam, Pacella and Higgins25 These mediators have also been reported to be associated with childhood maltreatment and/or neuropsychiatric outcomes in observational studies, indicating potential mediating pathways. To evaluate the role of systemic inflammation and chronic stress, four biological mediators were included (serum C-reactive protein (CRP), plasma cortisol, plasma interleukin-6 and interleukin-6 receptor), along with body mass index (BMI) as a general health indicator. Reference Baumeister, Akhtar, Ciufolini, Pariante and Mondelli2 For behavioural traits, seven lifestyle and behavioural traits relating to sedentary lifestyle and sleep problems were evaluated (self-reported: moderate/vigorous physical activity, TV watching hours, sleep duration, leisure screen time, propensity for speeding; accelerometer measured: average acceleration, sleep duration). Childhood maltreatment has known negative effects on different domains of cognitive functioning in children, which further manifests as risk-taking behaviours later in life, including engaging in addiction-related behaviours. As such, three cognitive traits (executive functioning, intelligence and internalising behaviour in childhood/adolescence), two risk-taking measures (risk tolerance and propensity for speeding) and nine risky health behaviours (addiction-related traits: smoking initiation, age of smoking initiation, pack years smoked, smoking cessation, cigarettes per day, cannabis use disorder, alcoholic drinks per week, alcohol dependency and ICD-coded substance abuse) were also assessed as mediators in this study, because these traits have been observationally associated with both childhood maltreatment and mental disorders. Reference Lunding, Ueland, Aas, Høegh, Werner and Rødevand26,Reference Aas, Etain, Bellivier, Henry, Lagerberg and Ringen27 Finally, to consider the potential protective effect of socioeconomic factors, household income and educational attainment were also included as mediators in the analyses. Additional details and references supporting the inclusion of these mediators are shown in Appendix 2.

CpG sites and mQTLs

On 13 June 2023, five childhood maltreatment-related phenotypes were searched in the EWAS Atlas (adverse childhood experience, child abuse, childhood adversity, childhood sexual victimisation and childhood stress) to identify 2929 childhood maltreatment-associated CpG sites at significance P < 1 × 10−4 from previous EWAS. Reference Li, Zou, Li, Gao, Sang and Zhang28 These were measured in blood, buccal cells, leukocytes, saliva and sperm samples. Methylation quantitative trait loci (mQTLs) were obtained from the Genetics of DNA Methylation Consortium (GoDMC) to proxy DNAm levels at these CpG sites. Reference Min, Hemani, Hannon, Dekkers, Castillo-Fernandez and Luijk29 DNAm data were generated using the Illumina Infinium HumanMethylation450 BeadChip from whole-blood samples of 27 750 European participants. Summary statistics of 36 cohorts were combined in an inverse-variance, fixed-effect meta-analysis to determine mQTLs for each CpG site. Reference Min, Hemani, Hannon, Dekkers, Castillo-Fernandez and Luijk29 Details can be found in Supplementary Table 1.

Instrument selection

Single-nucleotide polymorphisms (SNPs) were used as genetic instruments for childhood maltreatment at genome-wide significance (P < 5 × 10 −8). Linkage disequilibrium clumping (R 2 < 0.01, 1000 kb window), using the European 1000 Genomes reference panel, Reference Auton, Abecasis, Altshuler, Durbin, Abecasis and Bentley30 obtained independent SNPs. F-statistics were calculated to evaluate instrument strength, with F > 10 indicating strong instruments.

Independent cis-mQTLs (mQTL < 1 Mb from CpG site) at P < 1 × 10−8 were extracted from the GoDMC database to minimise potential horizontal pleiotropy. Reference Min, Hemani, Hannon, Dekkers, Castillo-Fernandez and Luijk29,Reference Richardson, Zheng, Davey Smith, Timpson, Gaunt and Relton31 This resulted in 1091 CpG sites with viable cis-mQTLs for epigenetic Mendelian randomisation.

Statistical analysis

Mediation Mendelian randomisation leverages two-step Mendelian randomisation and MVMR to separate the total causal effect of childhood maltreatment on each outcome into indirect and direct effects, respectively (Supplementary Fig. 1). Statistical analyses were conducted in R (version 4.3.1 for Windows) using these R packages: TwoSampleMR (version 0.6.8) Reference Hemani, Zheng, Elsworth, Wade, Haberland and Baird32 and MendelianRandomization (version 0.10.0) Reference Broadbent, Foley, Grant, Mason, Staley and Burgess33 for UVMR, ieugwasr (version 1.0.3) Reference Elsworth and Rasteiro34 for linkage disequilibrium clumping, MVMR (version 0.4) Reference Sanderson, Spiller and Bowden35 for MVMR, RMediation (version 1.2.2) Reference Tofighi and MacKinnon36 for calculation of indirect effect 95% confidence intervals and coloc (version 5.2.3) Reference Giambartolomei, Vukcevic, Schadt, Franke, Hingorani and Wallace37 for colocalisation. Total and direct effects are presented as odds ratio per unit increase in the Childhood Trauma Screener score. Indirect effects acting through mediators are presented as a percentage of the total effect of childhood maltreatment on each outcome.

A priori statistical power was calculated for Mendelian randomisation analyses with an online tool (https://sb452.shinyapps.io/power/; Supplementary Table 2). Calculation details can be found in Appendix 3. For ADHD as an example outcome, we would have sufficient power (>80%) if the expected odds ratio is ≤0.92 or ≥1.12. The power estimates for all other UVMR and epigenetic Mendelian randomisation associations can be found in Supplementary Table 2.

UVMR: total effects

The random-effects inverse-variance weighted (IVW) method was used as the main way to determine total effect (Supplementary Fig. 1A) to investigate the effect of childhood maltreatment on neuropsychiatric outcomes. Various tests were conducted to determine potential violations of the Mendelian randomisation assumptions. Reference Burgess, Davey Smith, Davies, Dudbridge, Gill and Glymour38 Steiger testing assessed whether genetic instruments explained more of the variance in childhood maltreatment than the outcome, and evaluated violations of assumption 1. This supported the use of selected instruments to investigate the direction of the causal association between childhood maltreatment and neuropsychiatric outcomes. F-statistics are a test for instrument strength, and were used to evaluate violations of assumption 1. Cochran’s Q test for instrument heterogeneity and the Mendelian randomisation-Egger intercept test evaluated violations of assumptions 2 and 3. Further pleiotropy-robust estimates were derived using the weighted median and Mendelian randomisation pleiotropy residual sum and outlier methods. Additional details of each sensitivity analysis are shown in Appendix 2.

Given that outcomes were selected based on prior knowledge, and therefore considered hypothesis-based analyses, UVMR effects under the IVW method that passed the Benjamini–Hochberg-corrected P-value threshold (P < 0.05) were considered significant. Of these, those with directionally concordant sensitivity analyses at P < 0.05 were considered robust and were included in subsequent mediation Mendelian randomisation analyses.

To address the key limitation in previous studies, Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1 we further calculated pairwise genetic correlation between childhood maltreatment and each outcome, and their corresponding linkage disequilibrium score regression intercepts. The intercept estimates the inflation in genetic correlation between two traits due to sample overlap. Reference Bulik-Sullivan, Finucane, Anttila, Gusev, Day and Loh39,Reference Lee, McGue, Iacono and Chow40 Bivariate genetic analyses were carried out using linkage disequilibrium score (LDSC) regression software (https://github.com/bulik/ldsc). Intercepts for each pairwise regression were considered. Those with confidence intervals that do not overlap with zero suggest that there may be minimal or residual confounding of the genetic relationship between the two traits because of sample overlap.

Two-step Mendelian randomisation: indirect effects

Indirect effects through each mediator were estimated using two-step Mendelian randomisation (Supplementary Fig. 1B). Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 In step one, UVMR was conducted to determine the causal effect of childhood maltreatment on each mediator. In step two, MVMR was conducted to determine the causal effect of mediators on each outcome after adjusting for childhood maltreatment. Standard errors and corresponding 95% confidence intervals were calculated using the Delta method. Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 The proportion mediated by each mediator is the indirect effect divided by the total effect. IVW estimates for each step were multiplied to obtain indirect effects when nominally significant associations in the same direction were obtained in both steps, with significant and directionally concordant sensitivity analyses. Significant mediators were defined as traits with indirect effects with confidence intervals excluding the null Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 (Supplementary Fig. 2).

MVMR: direct effects

MVMR simultaneously estimates the direct effects of multiple exposures on an outcome (Supplementary Fig. 1C). Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 MVMR serves two purposes: first, to generate direct effects of mediators on the outcome for step two of two-step Mendelian randomisation; second, to obtain the direct effect of childhood maltreatment accounting for mediators.

As the main method, MVMR-IVW was used to obtain the direct effects of each exposure. Modified Q-statistic (Qa min), conditional Sanderson–Windmeijer conditional F-statistic (F SW) and the MVMR-Egger intercept test were conducted to evaluate MVMR assumptions. MVMR-weighted median was conducted to obtain direct effects under relaxed assumptions on horizontal pleiotropy. Qhet_mvmr was conducted to consider the effects of weak instruments. Reference Sanderson, Spiller and Bowden35 Assumptions and details of MVMR are provided in Appendix 2.

Robust findings are indicated through directional concordance between MVMR sensitivity analyses and MVMR-IVW estimates, and significant MVMR sensitivity analyses with confidence intervals excluding the null.

Epigenetic Mendelian randomisation and colocalisation

UVMR analysis was conducted to determine the causal effect of DNAm at each childhood maltreatment-associated CpG site on significant outcomes (Supplementary Fig. 1D). Statistical significance was determined by Bonferroni-corrected P-values (ADHD: P < 6.68 × 10−5 (748 sites); MDD: P < 6.68 × 10−5 (748 sites)); schizophrenia: P < 6.59 × 10−5 (759 sites)). The stringent Bonferroni-corrected P-value thresholds were used to reduce false positives because this was a hypothesis-free analysis.

Colocalisation was conducted as a sensitivity analysis for associations passing the Bonferroni-corrected P-value threshold. Reference Zuber, Grinberg, Gill, Manipur, Slob and Patel41 colocalisation determines the probability that CpG sites and outcomes share causal variants (PP.H4). This was done within a ±200-kilobase (kb) window of either the nearest gene mapped using the Infinium HumanMethylation450 Manifest file 42 or the CpG site position. Default prior probabilities (Prior.H3 = 1 × 10−4, Prior.H4 = 1 × 10−5) were used to calculate posterior probabilities with the approximate Bayes factor. Reference Giambartolomei, Vukcevic, Schadt, Franke, Hingorani and Wallace37 Strong evidence of colocalisation was defined as PP.H4 > 0.8. Reference Rasooly, Peloso, Pereira, Dashti, Giambartolomei and Wheeler43 CpG sites with statistically significant epigenetic Mendelian randomisation estimates and passing colocalisation (PP.H4 > 0.8) were included in downstream pathway analysis. CpG sites were mapped to the nearest gene, combined into a gene set and searched for enriched pathways in https://string-db.org (version 12.0). Reference Szklarczyk, Kirsch, Koutrouli, Nastou, Mehryary and Hachilif44

Results

Genetic liability to childhood maltreatment on the risk of neuropsychiatric disorders

There were six genetic variants used as instruments for childhood maltreatment, with an average F-statistic of 32.63 (Supplementary Table 3), suggesting strong instruments. The genetic liability to childhood maltreatment on each outcome risk is shown in Supplementary Fig. 3 and Supplementary Table 4. Using the IVW method, genetic liability to childhood maltreatment was associated with higher risk of ADHD (odds ratio 10.09, 95% CI: 4.76–21.40), MDD (odds ratio 1.89, 95% CI: 1.32–2.70) and schizophrenia (odds ratio 5.89, 95% CI: 1.46–23.78) at Benjamini–Hochberg-corrected P-values. There was no causal evidence for the effect of childhood maltreatment on the risks of anxiety disorder, autism spectrum disorder and bipolar disorder. All associations passed Steiger testing. Generally, sensitivity estimates were consistent with the IVW estimate, and sensitivity tests did not suggest violated Mendelian randomisation assumptions. Detailed UVMR results are given in Appendix 4. To assess sample overlap, LDSC regression analysis was conducted for all childhood maltreatment to outcome pairs (Supplementary Table 5). Across all bivariate LDSC regression analyses, only childhood maltreatment and depression had an intercept of 0.065 (95% CI: 0.052–0.078) that did not overlap with zero, suggesting residual confounding possibly due to sample overlap. This suggests that all UVMR analyses, except childhood maltreatment and depression, were robust to residual confounding due to sample overlap.

Two-step Mendelian randomisation: indirect effect of mediators

In step one, 11 mediators were found to have causal associations under the IVW method (Supplementary Fig. 4 and Supplementary Table 6). Among biological traits, childhood maltreatment was associated with higher BMI and serum interleukin-6 levels. Among addiction-related traits, childhood maltreatment was associated with increased risks of smoking initiation, smoking cessation and substance abuse. Among cognitive and socioeconomic traits, childhood maltreatment was associated with increased internalising behaviour in childhood and adolescence, risk tolerance and leisure screen time, and negatively associated with educational attainment, executive functioning and intelligence. All associations passed Steiger testing. Seven mediators passed sensitivity analyses: smoking initiation, smoking cessation, educational attainment, executive functioning, intelligence, leisure screen time and risk tolerance (Supplementary Table 6).

In step two, the MVMR-IVW method determined the effect of mediators, independent of childhood maltreatment, on each outcome (Supplementary Table 7).

Smoking initiation, substance abuse and increased leisure screen time were significantly associated with higher ADHD risk, while higher educational attainment, executive functioning and intelligence were significantly associated with lower ADHD risk (Table 1). These effects are directionally consistent with step 1. Smoking initiation (18.98%) and substance abuse (20.38%) mediated the largest proportion of the effect of childhood maltreatment on ADHD.

Table 1 Indirect effects of significant mediators for childhood maltreatment on neuropsychiatric outcome relationship

ADHD, attention-deficit/hyperactivity disorder; MDD, major depressive disorder; NA, not available.

The effect size in step one is the causal effect of childhood maltreatment on the investigated mediator (reported in beta). The effect size in step two is the causal effect of the mediator on the neuropsychiatric outcome (reported in beta). The indirect effect, the product of effect sizes from steps one and two, was calculated if nominally significant associations in the same direction were obtained in both steps, with significant and directionally concordant sensitivity analyses (reported in beta). The total effect is the causal effect of childhood maltreatment on each neuropsychiatric outcome (reported in beta). The proportion mediated by the investigated mediator is the indirect effect divided by the total effect.

Increased leisure screen time was significantly associated with higher risk of MDD, while higher educational attainment and executive functioning were significantly associated with lower MDD risk (Table 1). The directions of these associations were consistent with step 1. These three mediators mediated a similar proportion of the effect of childhood maltreatment on MDD (educational attainment 6.51%; executive functioning 9.27%; leisure screen time 9.90%).

Higher risk tolerance and educational attainment were significantly associated with higher schizophrenia risk, while increased executive functioning and intelligence were significantly associated with lower schizophrenia risk (Table 1). The directions of these associations for all mediators, except educational attainment, were consistent with those in step 1. Risk tolerance (22.17%) and executive functioning (19.26%) mediated the largest proportions of the relationship between childhood maltreatment and schizophrenia.

MVMR sensitivity analyses assessed the robustness of these findings in the context of Mendelian randomisation assumptions. MVMR-weighted median sensitivity analyses for all models were consistent with MVMR-IVW estimates, suggesting that the findings were robust to unbalanced horizontal pleiotropy. Models with smoking initiation and childhood maltreatment as exposures had conditional F-statistics of 11.55 and 10.27, respectively, suggesting that weak instruments would not bias the results (Supplementary Table 7). For other models, at least one exposure had low conditional F-statistics (F sw < 10). Qhet_mvmr estimates were similar to respective MVMR-IVW estimates, suggesting that weak instruments did not bias the MVMR results (Supplementary Table 7).

MVMR: direct effects of childhood maltreatment

To evaluate the direct effects of childhood maltreatment, each mediator was included in the MVMR model independently, then statistically significant mediators were added to a full MVMR model.

Figure 2 shows the total effect of childhood maltreatment from UVMR, and the direct effects of childhood maltreatment from MVMR on each outcome (Supplementary Table 8). For ADHD, adjustment for the combined effects of all statistically significant mediators resulted in a smaller direct effect (odds ratio 2.57, 95% CI: 1.76–3.77) than when individually including mediators, and compared with the total effect (Fig. 2a). For schizophrenia, only the inclusion of all significant mediators in a MVMR model resulted in a slightly lower direct effect of childhood maltreatment on schizophrenia (odds ratio 5.10, 95% CI: 2.71–9.60), compared with the total effect (Fig. 2c). By contrast, regardless of whether individual or all mediators were included in the MVMR model, the direct effect of childhood maltreatment on MDD (odds ratio 1.95, 95% CI: 1.68–2.26) remained largely unchanged from the total effect (Fig. 2b).

Fig. 2 Direct effects of childhood maltreatment (CM) adjusting for individual or all significant mediators. ADHD, attention-deficit/hyperactivity disorder; MDD, major depressive disorder; SNPs, single-nucleotide polymorphisms; OR, odds ratio.

F-statistics for the MVMR model that included all significant mediators ranged from 1.27 to 4.26 with ADHD as the outcome, 1.91–6.93 with MDD as the outcome and 1.77–7.02 with schizophrenia as the outcome (Supplementary Table 8), suggesting potential weak instrument bias. However, Qhet_mvmr estimates were directionally consistent across all analyses, suggesting that weak instrument bias did not affect MVMR results (Supplementary Table 8).

Epigenetic Mendelian randomisation and colocalisation

There were 26 statistically significant associations between CpG sites and ADHD (3), MDD (4) and schizophrenia (19) at Bonferroni-corrected P-value thresholds (Supplementary Table 9). Most CpG sites were associated with only one outcome, except for two sites. Cg23337620 (not near any gene) was associated with ADHD (odds ratio 0.55, 95% CI: 0.42–0.71) and schizophrenia (odds ratio 0.69, 95% CI: 0.58–0.82). Cg15790214 (HCG11) was associated with MDD (odds ratio 0.91, 95% CI: 0.88–0.93) and schizophrenia (odds ratio 0.72, 95% CI: 0.68–0.77). Significant CpG sites were included in downstream colocalisation and pathway analysis.

Colocalisation analysis was performed to determine whether CpG sites and outcomes shared the same causal variant. There was strong evidence of colocalisation with 10 CpG sites (PP.H4 > 0.8) (Supplementary Table 10). For these, the mQTLs used to proxy DNAm of two CpG sites in Mendelian randomisation were found to have 100% posterior probability of sharing causal variants with schizophrenia: cg04585390 (SLC5A10) and cg11329058 (RPTOR).

For these ten CpG sites with strong evidence of colocalisation, nine were situated near one gene whereas one was not near any (Supplementary Table 10). These genes are ANGPLT2, CLU, H2BC6, HNRNPK, MAPT, NDFIP2, RPTOR, SLC5A10 and TSSC1. Although no enrichment of protein–protein interactions was observed, three genes – CLU, MAPT and HNRNPK – were implicated in the formation of neurofibrillary tangles.

Discussion

Main findings

To our knowledge, this is the first study to delineate mediators of the relationship between childhood maltreatment and neuropsychiatric disorders. We replicated causal associations between childhood maltreatment and ADHD, MDD and schizophrenia from previous Mendelian randomisation studies. Reference Warrier, Kwong, Luo, Dalvie, Croft and Sallis1,Reference Alkema, Marchi, van der Zaag, van der Sluis, Warrier and Ophoff22,Reference Baltramonaityte, Pingault, Cecil, Choudhary, Järvelin and Penninx23 We established addiction-related, cognitive and socioeconomic traits as mediators (Fig. 3). The direct effect of childhood maltreatment persisted, particularly for MDD. Epigenetic Mendelian randomisation revealed ten childhood maltreatment-associated CpG sites associated with outcomes, highlighting novel pathways for further investigation.

Fig. 3 Breakdown of the mediators of the relationship between childhood maltreatment (CM) and neuropsychiatric disorders, potential mechanisms and suggested interventions. Mediators in red are addiction-related traits, mediators in green are cognitive-related traits and the mediator in blue is a socioeconomic trait. Solid arrows represent causal links and the dashed arrow represents associations, not causality. DNAm, DNA methylation; ADHD, attention-deficit/hyperactivity disorder; MDD, major depressive disorder; OR, odds ratio.

Interpretation of our findings and comparison with the literature

Our results support the associations between childhood maltreatment and mediators reported in observational studies. Childhood maltreatment was reported to reduce cognitive functioning, intelligence Reference Matte-Landry, Grisé Bolduc, Tanguay-Garneau, Collin-Vézina and Ouellet-Morin45 and academic achievement, Reference Fry, Fang, Elliott, Casey, Zheng and Li46 and to increase tobacco and illicit drug use. Reference Petruccelli, Davis and Berman8 The magnitude of these differed slightly compared with our Mendelian randomisation estimates. Genetic instruments used in Mendelian randomisation studies affect the usual levels of exposure over long periods of time, corresponding to a lifetime effect of the exposure. Reference Burgess, Davey Smith, Davies, Dudbridge, Gill and Glymour38 For example, educational attainment in observational studies has larger effect sizes compared with our Mendelian randomisation estimates. This may reflect childhood maltreatment impact during short yet sensitive periods, when brain development is highly susceptible to external stimuli. However, childhood maltreatment occurs during sensitive periods in childhood development, specifically when brain development is highly susceptible to external stimuli. Reference Schaefer, Cheng and Dunn47 This could explain the larger effects on educational outcomes in observational studies as compared with our Mendelian randomisation estimates. Nevertheless, the directionality of our results was consistent with observational associations, validating our findings.

The results in step two of mediation Mendelian randomisation align with literature that has reported a significant causal effect of executive functioning, educational attainment and smoking initiation on ADHD, MDD and schizophrenia. Reference Ahn, Norman, Justice and Shaw48Reference Treur, Demontis, Smith, Sallis, Richardson and Wiers51

In applying novel Mendelian randomisation methods to delineate mediators, we validated hypothesised mechanisms linking childhood maltreatment and outcomes. Previous studies delineated childhood maltreatment as a predictor of reduced executive functioning, Reference Letkiewicz, Funkhouser and Shankman52,Reference Masson, Bussières, East-Richard, R-Mercier and Cellard53 and Mendelian randomisation studies showed that this increased the risks of ADHD, MDD and schizophrenia. Reference Ahn, Norman, Justice and Shaw48,Reference Burton, Sallis, Hatoum, Munafò and Reed49 We add to this by establishing temporal ordering between childhood maltreatment, executive functioning and neuropsychiatric outcomes. Early life stressors impair executive functioning through atypical development of brain regions including the prefrontal cortex, which influences emotional regulation and reward processing. Reference Kasai, Koike, Matsuoka, Okada, Sakakibara and Sakurada54 Reduced educational attainment exacerbates these cognitive deficits. Reference Fry, Fang, Elliott, Casey, Zheng and Li46,Reference Masson, Bussières, East-Richard, R-Mercier and Cellard53,Reference Puetz and McCrory55 These mechanisms contribute to reduced stress tolerance and the inability to disengage from rumination, resulting in increased risk of adult neuropsychiatric disorders. Reference Rock, Roiser, Riedel and Blackwell56

In observational studies, childhood maltreatment was reported to lower educational attainment. Reference Fry, Fang, Elliott, Casey, Zheng and Li46,Reference Romano, Babchishin and Marquis57 Adding to this, we found that educational attainment mediated the effect of childhood maltreatment on ADHD and MDD. At the individual level, exposure to childhood maltreatment impacts brain and cognitive development, causing mental health problems. Reference Fry, Fang, Elliott, Casey, Zheng and Li46 Challenges associated with mental health further limit educational attainment and, consequently, cognitive impairments persist. At a social level, childhood maltreatment-exposed children are neither provided with social and emotional support nor exposed to environments fostering a sense of agency or self-efficacy. Reference Romano, Babchishin and Marquis57,Reference Slade and Wissow58 As such, they are unable to internalise the value of education, leading to poorer educational outcomes. These suggest that educational, child welfare and care-giving sectors contribute to sustaining limited educational attainment and increased risk of neuropsychiatric outcomes. Reference Romano, Babchishin and Marquis57

Although smoking and alcohol dependency were commonly reported as consequences of childhood maltreatment, Reference Petruccelli, Davis and Berman8 our analyses supported only smoking initiation and substance abuse as mediators of the relationship between childhood maltreatment and ADHD. Non-significant indirect effects of other smoking traits may be due to limited power from smaller sample sizes. Alternatively, the self-medication hypothesis states that victims who experienced trauma use illicit substances as a coping mechanism to reduce awareness and maintain a dissociated state. Reference Schimmenti, Billieux, Santoro, Casale and Starcevic59 Thus, without healthy coping mechanisms, these behaviours can elicit unexpected effects – for example, reducing the risk of neuropsychiatric disorders. Reference Brady, Phalen, Roche, Cowan and Bennett60

We formally delineated novel mediators such as leisure screen time and risk tolerance. The role of risk tolerance can be explained by impaired cognitive functioning, whereas leisure screen time is more complex. Leisure screen time may be a form of avoidant coping hindering emotional regulation development. Reference Westerlund, Alvarsson, Carli and Hadlaczky61 However, its effect on outcomes remains uncertain due to the wide content variety and influences on coping strategies. Reference Westerlund, Alvarsson, Carli and Hadlaczky61,Reference Tang, Werner-Seidler, Torok, Mackinnon and Christensen62 Therefore, further research into the implication of leisure screen time as a mediator is required.

Another novelty of our study is the application of epigenetic Mendelian randomisation to establish outcome-specific DNAm signatures. Significant CpG sites were near genes responsible for biological pathways such as apoptosis (HNRNPK, CLU, RPTOR) and transcription regulation (H2BC6, HNRNPK).

Pathway analysis of these mapped genes identified a functionally enriched in neurofibrillary tangle formation pathways, involving CLU, MAPT and HNRNPK. Neurofibrillary tangles are a hallmark of neurodegenerative disorders, including Alzheimer’s and Parkinson’s diseases, acting as a precursor to microtubule destabilisation and dysregulated dynamics of the neuronal skeleton. Reference Metaxas and Kempf63 Our results showed that hypomethylation at the MAPT gene, corresponding to increased expression of the tau protein, was associated with increased schizophrenia risk. Reference Coupland, Mellick, Silburn, Mather, Armstrong and Sachdev64 HNRNPK, an RNA protein, influences the development of early tau pathology, Reference Kavanagh, Halder and Drummond65 and CLU (clusterin) binds and stabilises tau to prevent degradation, thereby accelerating tau aggregation and its neurotoxic properties. Reference Yuste-Checa, Trinkaus, Riera-Tur, Imamoglu, Schaller and Wang66 Interestingly, increased tau levels were previously associated with depression, anxiety and ADHD. Reference Hall, Petersen and Johnson67,Reference Demirci68 Contrastingly, both total and phosphorylated tau levels were lower in patients with schizophrenia. Reference Demirel, Cetin, Turan, Yıldız, Sağlam and Duran69 Nevertheless, tau pathology has a potential mechanistic role in the pathogenesis of neuropsychiatric disorders. We further identified one additional signal at the RPTOR (regulatory associated protein of MTOR complex 1) gene, which is part of the mTORC1 pathway. This stress-related pathway is reported to be downregulated in the prefrontal cortex of humans and rats with schizophrenia. Reference Chen, Guan, Wang and Lin70,Reference Sanacora, Yan and Popoli71 This aligns with our results, which show that childhood maltreatment-induced hypomethylation at RPTOR, triggering the MTORC1 pathway, was associated with increased schizophrenia risk. Considering the two key findings from our epigenetic Mendelian randomisation analysis, we hypothesise that early childhood adversities such as childhood maltreatment are linked to neuropsychiatric outcomes, possibly through changes in DNAm at key genes involved in neurobiological and/or stress pathways.

However, these epigenetic Mendelian randomisation results should be interpreted with caution. First, the full GWAS summary statistics for DNAm at each CpG site were not available, precluding formal two-step epigenetic mediation Mendelian randomisation analyses and delineation of CpG sites as true causal mediators. Nevertheless, the specificity of genetic instruments ensured specificity of the DNAm effect and allowed postulation of biological mechanisms. Significant CpG sites were near genes responsible for biological pathways such as apoptosis (HNRNPK, CLU, RPTOR) and transcription regulation (H2BC6, HNRNPK). Second, DNAm levels were proxied by at most three mQTLs, which explains only a small proportion of the variation in DNAm. This may overlook any true causal associations resulting from larger methylation changes that are not captured by the selected mQTLs. Third, we considered DNAm levels only in whole blood, potentially missing DNAm signatures in other relevant tissues, including the brain. Nevertheless, we considered all CpG sites associated with any adverse childhood exposures in a hypothesis-free approach to identify potential novel epigenetic mediators. We recognise that epigenetic Mendelian randomisation alone is insufficient to establish biological mechanisms. Therefore, these findings warrant further research to confirm the biological role of each gene in the pathogenesis of ADHD, MDD and schizophrenia in victims. This may identify new markers to screen for high-risk victims, or novel therapeutic targets to reduce the risk of neuropsychiatric disorders.

Clinical implications

Despite SDG Target 16.2 aiming to end abuse against children by 2030, Reference Lee and Kim4,Reference Clark, Coll-Seck, Banerjee, Peterson, Dalglish and Ameratunga72 a large proportion of the adult population is expected to have a history of childhood maltreatment following the COVID-19 pandemic. The direct effect on neuropsychiatric outcomes is indisputable, and therefore strengthening care-giving and economic support for high-risk families would prevent childhood maltreatment by creating a safe environment. Reference Malvaso, Pilkington, Montgomerie, Delfabbro and Lynch73 The consequences of childhood maltreatment can be mitigated through early detection and intervention. Screening for childhood maltreatment history in clinical settings may be an effective strategy to identify high-risk individuals for adverse neuropsychiatric outcomes. 74 Through this, timely provision of trauma-informed therapy and increased access to community trauma services can prevent further childhood maltreatment and mitigate negative consequences.

Addiction-related, cognitive and socioeconomic mediators may be interdependent. For example, without healthy coping strategies, successful prevention of addiction in the childhood maltreatment victim may lead them to engage in other unhealthy behaviours. Therefore, policies and interventions need to target the multidimensional consequences of childhood maltreatment. Simultaneous addiction prevention and provision of psychoeducation and skills to cope with trauma was proposed as an effective intervention. Reference Bailey, Newton, Perry, Grummitt, Baams and Barrett75 ‘Rise Above’ is an example of an evidence-based trauma-informed prevention programme. This is a comprehensive intervention curriculum for youths covering trauma, education on e-cigarette and drug use, psychosocio-emotional development and prosocial behaviour. Reference Shin76 In providing the necessary coping tools for childhood maltreatment victims, the detrimental impacts on mental health can be mitigated. Moreover, psychotherapies, including trauma-focused cognitive behavioural therapy, facilitate the development of cognitive coping strategies to improve self-regulation, while cognitive remediation focuses on strategies that improves cognitive functioning. Reference Masson, Bussières, East-Richard, R-Mercier and Cellard53,Reference Rhoades, Mitnick, Heyman, Slep, Del Vecchio, Wampler and McWey77 Furthermore, close coordination between educational, child welfare and care-giving sectors is required to enable access to education, and to foster a positive attitude towards education in childhood maltreatment victims. Reference Romano, Babchishin and Marquis57 Therefore, a holistic approach is required to reach the therapeutic potential of trauma-informed interventions.

Strengths and limitations

This study used mediation and epigenetic Mendelian randomisation to provide novel insight into the pathways between childhood maltreatment and neuropsychiatric outcomes. It overcomes the limitations of unmeasured confounding and reverse causation in traditional observational study designs to delineate key causal mediators, effectively prioritising areas for intervention.

Our findings should be taken in the context of several limitations. First, while the Mendelian randomisation framework generates causal evidence, one must satisfy the Mendelian randomisation assumptions. Nevertheless, sensitivity analyses and instrument strength were satisfactory for UVMR analyses, and MVMR analyses accounting for weak instruments showed unlikely biased results. Second, childhood maltreatment and some mediators were based on self-reported measurements. Incorrect group assignment or reporting bias in the original study may have biased our results. Third, childhood maltreatment was reported as a composite measure of childhood maltreatment subtypes. We could not evaluate the effects of childhood maltreatment subtypes on outcomes. Fourth, the direct and indirect effects of binary mediators may be inaccurate given the non-collapsibility of odds ratios. However, a recent study showed that log(odds ratio) could be used as an outcome unit to calculate indirect effects. Reference Carter, Sanderson, Hammerton, Richmond, Davey Smith and Heron16 Fifth, our analyses were limited by the type and availability of data. GWAS data-sets used in Mendelian randomisation analyses were limited to European ancestry, and caution is required when generalising these findings to other populations. As discussed, formal mediation epigenetic Mendelian randomisation analyses were not possible with the current data. Furthermore, mQTLs were also limited to regions within 1 Mb from the CpG site, resulting in very few instruments to proxy DNAm levels. Although this may have limited instrument strength, the specificity of genetic instruments enabled postulation of biological mechanisms. Sixth, the effect sizes presented in this study represent a lifetime average effect of exposures on outcomes. Thus, mediators with large effects acting within shorter time frames may have been overlooked in this study. Finally, the heritability of childhood maltreatment is moderate, and genetic instruments represent gene–environmental interactions such as transgenerational transmission and passive or evocative heritability. Reference Espinosa Dice, Lawn, Ratanatharathorn, Roberts, Denckla and Kim78 Because genetics are more stable and predictable than environmental factors, our findings are based on childhood maltreatment variance, as determined by the myriad of childhood maltreatment risk factors with which the variants are associated. Reference ter Kuile, Hübel, Cheesman, Coleman, Peel and Levey79 While the genetic and environmental causes of childhood maltreatment are not distinguished, our findings determine the consequences of childhood maltreatment regardless of the underlying cause. Such analyses reflect how childhood maltreatment is perpetrated, and present crucial evidence to guide interventions.

In conclusion, our results support the role of addiction-related, cognitive and socioeconomic traits in mediating the effect of childhood maltreatment on neuropsychiatric disorders. Responsive interventions in these areas may be more relevant for ADHD and schizophrenia, but not for MDD. These include addiction prevention and treatment, improving cognitive functioning and supporting socioeconomic attainment. Ultimately, preventative strategies are vital to reducing childhood maltreatment prevalence and its long-term impacts. Reference Ng80

Supplementary material

The supplementary material is available online at https://doi.org/10.1192/bjp.2025.10433

Data availability

All genetic instruments can be obtained from summary statistics from the respective genome-wide association study (GWAS). Outcome estimates were obtained using publicly accessible data made accessible by the GWAS consortia.

Acknowledgements

The authors thank the researchers associated with each GWAS for providing publicly available summary statistics used for the analyses in this study. The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975 as revised in 2013. All procedures involving human participants can be found in the original studies. This manuscript includes material adapted from a Master of Philosophy thesis submitted to the University of Hong Kong in 2024. Necessary permissions for reuse have been obtained.

Author contributions

N.Y.T.N. and B.H.Y.C. conceived the study. N.Y.T.N., J.V.Z., C.C.Y.M. and B.H.Y.C. contributed to the study design. N.Y.T.N. and J.Y.T. conducted the main analysis and drafted the manuscript. J.Y.T. conducted sensitivity analyses and generated the figures. All authors were involved in the interpretation of results, helped refine the manuscript and approved its final version. B.H.Y.C. is the guarantor and attests that all authors meet authorship criteria.

Funding

This work was supported by Seed Funding for Basic Research (no. 2202100680) from the University of Hong Kong.

Declaration of interest

None.

Transparency declaration

N.Y.T.N. and B.H.Y.C. affirm that the manuscript is an honest, accurate and transparent account of the reported study; that no important aspects of the study has been omitted; and that any discrepancies from the study as planned have been explained.

Analytic code availability

Code used to carry out statistical analyses can be found in the respective packages and are available upon reasonable request.

References

Warrier, V, Kwong, ASF, Luo, M, Dalvie, S, Croft, J, Sallis, HM, et al. Gene-environment correlations and causal effects of childhood maltreatment on physical and mental health: a genetically informed approach. Lancet Psychiatry 2021; 8: 373–86.10.1016/S2215-0366(20)30569-1CrossRefGoogle ScholarPubMed
Baumeister, D, Akhtar, R, Ciufolini, S, Pariante, CM, Mondelli, V. Childhood trauma and adulthood inflammation: a meta-analysis of peripheral C-reactive protein, interleukin-6 and tumour necrosis factor-α. Mol Psychiatry 2016; 21: 642–9.10.1038/mp.2015.67CrossRefGoogle ScholarPubMed
Park, WJ, Walsh, KA. COVID-19 and the unseen pandemic of child abuse. BMJ Paediatr Open 2022; 6: e001553.10.1136/bmjpo-2022-001553CrossRefGoogle ScholarPubMed
Lee, H, Kim, E. Global prevalence of physical and psychological child abuse during COVID-19: a systematic review and meta-analysis. Child Abuse Negl 2023; 135: 105984.10.1016/j.chiabu.2022.105984CrossRefGoogle ScholarPubMed
Baldwin, JR, Wang, B, Karwatowska, L, Schoeler, T, Tsaligopoulou, A, Munafò, MR, et al. Childhood maltreatment and mental health problems: a systematic review and meta-analysis of quasi-experimental studies. Am J Psychiatry 2023; 180: 117–26.10.1176/appi.ajp.20220174CrossRefGoogle ScholarPubMed
Carr, A, Duff, H, Craddock, F. A systematic review of reviews of the outcome of noninstitutional child maltreatment. Trauma Viol Abuse 2020; 21: 828–43.10.1177/1524838018801334CrossRefGoogle ScholarPubMed
Macpherson, JM, Gray, SR, Ip, P, McCallum, M, Hanlon, P, Welsh, P, et al. Child maltreatment and incident mental disorders in middle and older ages: a retrospective UK Biobank cohort study. Lancet Reg Health Eur 2021; 11: 100224.10.1016/j.lanepe.2021.100224CrossRefGoogle Scholar
Petruccelli, K, Davis, J, Berman, T. Adverse childhood experiences and associated health outcomes: a systematic review and meta-analysis. Child Abuse Negl 2019; 97: 104127.10.1016/j.chiabu.2019.104127CrossRefGoogle ScholarPubMed
Power, C, Li, L, Pinto Pereira, SM. An overview of child maltreatment (neglect and abuse) associations with developmental trajectories and long-term outcomes in the 1958 British birth cohort. Longit Life Course Stud 2020; 11: 431–58.10.1332/175795920X15891281805890CrossRefGoogle ScholarPubMed
Strathearn, L, Giannotti, M, Mills, R, Kisely, S, Najman, J, Abajobir, A. Long-term cognitive, psychological, and health outcomes associated with child abuse and neglect. Pediatrics 2020; 146: e20200438.10.1542/peds.2020-0438CrossRefGoogle ScholarPubMed
Alameda, L, Rodriguez, V, Carr, E, Aas, M, Trotta, G, Marino, P, et al. A systematic review on mediators between adversity and psychosis: potential targets for treatment. Psychol Med 2020; 50: 1966–76.10.1017/S0033291720002421CrossRefGoogle ScholarPubMed
Klumparendt, A, Nelson, J, Barenbrügge, J, Ehring, T. Associations between childhood maltreatment and adult depression: a mediation analysis. BMC Psychiatry 2019; 19: 36.10.1186/s12888-019-2016-8CrossRefGoogle ScholarPubMed
Keyes, KM, Eaton, NR, Krueger, RF, McLaughlin, KA, Wall, MM, Grant, BF, et al. Childhood maltreatment and the structure of common psychiatric disorders. Br J Psychiatry 2012; 200: 107–15.10.1192/bjp.bp.111.093062CrossRefGoogle ScholarPubMed
Li, ET, Luyten, P, Midgley, N. Psychological mediators of the association between childhood emotional abuse and depression: a systematic review. Front Psychiatry 2020; 11: 559213.10.3389/fpsyt.2020.559213CrossRefGoogle ScholarPubMed
Kent, KR, Hopfer, CJ, Corley, RP, Stallings, MC. Executive functions as a mediator of childhood maltreatment on adult psychopathology: a longitudinal mediation analysis comparing maltreatment factor models. Child Abuse Negl 2022; 123: 105369.10.1016/j.chiabu.2021.105369CrossRefGoogle ScholarPubMed
Carter, AR, Sanderson, E, Hammerton, G, Richmond, RC, Davey Smith, G, Heron, J, et al. Mendelian randomisation for mediation analysis: current methods and challenges for implementation. Eur J Epidemiol 2021; 36: 465–78.10.1007/s10654-021-00757-1CrossRefGoogle ScholarPubMed
Yang, B, Zhang, H, Ge, W, Weder, N, Douglas-Palumberi, H, Perepletchikova, F, et al. Child abuse and epigenetic mechanisms of disease risk. Am J Prev Med 2013; 44: 101–7.10.1016/j.amepre.2012.10.012CrossRefGoogle ScholarPubMed
Parade, SH, Huffhines, L, Daniels, TE, Stroud, LR, Nugent, NR, Tyrka, AR. A systematic review of childhood maltreatment and DNA methylation: candidate gene and epigenome-wide approaches. Transl Psychiatry 2021; 11: 134.10.1038/s41398-021-01207-yCrossRefGoogle ScholarPubMed
Smearman, EL, Almli, LM, Conneely, KN, Brody, GH, Sales, JM, Bradley, B, et al. Oxytocin receptor genetic and epigenetic variations: association with child abuse and adult psychiatric symptoms. Child Dev 2016; 87: 122–34.10.1111/cdev.12493CrossRefGoogle ScholarPubMed
Comtois-Cabana, M, Barr, E, Provençal, N, Ouellet-Morin, I, Provenzi, L. Association between child maltreatment and depressive symptoms in emerging adulthood: the mediating and moderating roles of DNA methylation. PLOS One 2023; 18: e0280203.10.1371/journal.pone.0280203CrossRefGoogle ScholarPubMed
Bowden, J, Del Greco, M F, Minelli, C, Davey Smith, G, Sheehan, NA, Thompson, JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. Int J Epidemiol 2016; 45: 1961–74.Google ScholarPubMed
Alkema, A, Marchi, M, van der Zaag, JAJ, van der Sluis, D, Warrier, V, Ophoff, RA, et al. Childhood abuse v. neglect and risk for major psychiatric disorders. Psychol Med 2023; 54: 1598–609.10.1017/S0033291723003471CrossRefGoogle ScholarPubMed
Baltramonaityte, V, Pingault, J-B, Cecil, CAM, Choudhary, P, Järvelin, M-R, Penninx, BWJH, et al. A multivariate genome-wide association study of psycho-cardiometabolic multimorbidity. PLOS Genet 2023; 19: e1010508.10.1371/journal.pgen.1010508CrossRefGoogle ScholarPubMed
Sullivan, PF, Agrawal, A, Bulik, CM, Andreassen, OA, Borglum, AD, Breen, G, et al. Psychiatric genomics: an update and an agenda. Am J Psychiatry 2018; 175: 1527.10.1176/appi.ajp.2017.17030283CrossRefGoogle Scholar
Scott, JG, Malacova, E, Mathews, B, Haslam, DM, Pacella, R, Higgins, DJ, et al. The association between child maltreatment and mental disorders in the Australian Child Maltreatment Study. Med J Aust 2023; 218(suppl 6): S2633.10.5694/mja2.51870CrossRefGoogle ScholarPubMed
Lunding, SH, Ueland, T, Aas, M, Høegh, MC, Werner, MCF, Rødevand, L, et al. Tobacco smoking related to childhood trauma mediated by cognitive control and impulsiveness in severe mental disorders. Schizophr Res 2023; 261: 236–44.10.1016/j.schres.2023.09.041CrossRefGoogle ScholarPubMed
Aas, M, Etain, B, Bellivier, F, Henry, C, Lagerberg, T, Ringen, A, et al. Additive effects of childhood abuse and cannabis abuse on clinical expressions of bipolar disorders. Psychol Med 2014; 44: 1653–62.10.1017/S0033291713002316CrossRefGoogle ScholarPubMed
Li, M, Zou, D, Li, Z, Gao, R, Sang, J, Zhang, Y, et al. EWAS Atlas: a curated knowledgebase of epigenome-wide association studies. Nucleic Acids Res 2019; 47: D9838.10.1093/nar/gky1027CrossRefGoogle Scholar
Min, JL, Hemani, G, Hannon, E, Dekkers, KF, Castillo-Fernandez, J, Luijk, R, et al. Genomic and phenotypic insights from an atlas of genetic effects on DNA methylation. Nat Genet 2021; 53: 1311–21.10.1038/s41588-021-00923-xCrossRefGoogle ScholarPubMed
Auton, A, Abecasis, GR, Altshuler, DM, Durbin, RM, Abecasis, GR, Bentley, DR, et al. A global reference for human genetic variation. Nature 2015; 526: 6874.Google ScholarPubMed
Richardson, TG, Zheng, J, Davey Smith, G, Timpson, NJ, Gaunt, TR, Relton, CL, et al. Mendelian randomization analysis identifies cpG sites as putative mediators for genetic influences on cardiovascular disease risk. Am J Hum Genet 2017; 101: 590602.10.1016/j.ajhg.2017.09.003CrossRefGoogle ScholarPubMed
Hemani, G, Zheng, J, Elsworth, B, Wade, KH, Haberland, V, Baird, D, et al. The MR-Base platform supports systematic causal inference across the human phenome. eLife 2018; 7: e34408.10.7554/eLife.34408CrossRefGoogle ScholarPubMed
Broadbent, JR, Foley, CN, Grant, AJ, Mason, AM, Staley, JR, Burgess, S. MendelianRandomization v0.5.0: updates to an R package for performing Mendelian randomization analyses using summarized data. Wellcome Open Res 2020; 5: 252.10.12688/wellcomeopenres.16374.2CrossRefGoogle Scholar
Elsworth, BP, Rasteiro, R. ieugwasr: Interface to the ‘OpenGWAS’ Database API. Github, 2025 (https://github.com/MRCIEU/ieugwasr).Google Scholar
Sanderson, E, Spiller, W, Bowden, J. Testing and correcting for weak and pleiotropic instruments in two-sample multivariable Mendelian randomization. Stat Med 2021; 40: 5434–52.10.1002/sim.9133CrossRefGoogle ScholarPubMed
Tofighi, D, MacKinnon, DP. RMediation: an R package for mediation analysis confidence intervals. Behav Res Methods 2011; 43: 692700.10.3758/s13428-011-0076-xCrossRefGoogle Scholar
Giambartolomei, C, Vukcevic, D, Schadt, EE, Franke, L, Hingorani, AD, Wallace, C, et al. Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLOS Genet 2014; 10: e1004383.10.1371/journal.pgen.1004383CrossRefGoogle ScholarPubMed
Burgess, S, Davey Smith, G, Davies, NM, Dudbridge, F, Gill, D, Glymour, MM, et al. Guidelines for performing Mendelian randomization investigations: update for summer 2023. Wellcome Open Res 2019; 4: 186.10.12688/wellcomeopenres.15555.1CrossRefGoogle ScholarPubMed
Bulik-Sullivan, B, Finucane, HK, Anttila, V, Gusev, A, Day, FR, Loh, PR, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet 2015; 47: 1236–41.10.1038/ng.3406CrossRefGoogle ScholarPubMed
Lee, JJ, McGue, M, Iacono, WG, Chow, CC. The accuracy of LD score regression as an estimator of confounding and genetic correlations in genome-wide association studies. Genet Epidemiol 2018; 42: 783–95.10.1002/gepi.22161CrossRefGoogle Scholar
Zuber, V, Grinberg, NF, Gill, D, Manipur, I, Slob, EAW, Patel, A, et al. Combining evidence from Mendelian randomization and colocalization: review and comparison of approaches. Am J Hum Genet 2022; 109: 767–82.10.1016/j.ajhg.2022.04.001CrossRefGoogle ScholarPubMed
Illumina. HumanMethylation 450 v1.2 Manifest File (CSV Format). Illumina, 2013 (https://sapac.support.illumina.com/downloads/infinium_humanmethylation450_product_files.html [cited 13 Jun 2023]).Google Scholar
Rasooly, D, Peloso, GM, Pereira, AC, Dashti, H, Giambartolomei, C, Wheeler, E, et al. Genome-wide association analysis and Mendelian randomization proteomics identify drug targets for heart failure. Nat Commun 2023; 14: 3826.10.1038/s41467-023-39253-3CrossRefGoogle ScholarPubMed
Szklarczyk, D, Kirsch, R, Koutrouli, M, Nastou, K, Mehryary, F, Hachilif, R, et al. The STRING database in 2023: protein-protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res 2023; 51: D63846.10.1093/nar/gkac1000CrossRefGoogle ScholarPubMed
Matte-Landry, A, Grisé Bolduc, M-È, Tanguay-Garneau, L, Collin-Vézina, D, Ouellet-Morin, I. Cognitive outcomes of children with complex trauma: a systematic review and meta-analyses of longitudinal studies. Trauma Violence Abuse 2022; 24: 2743–57.10.1177/15248380221111484CrossRefGoogle Scholar
Fry, D, Fang, X, Elliott, S, Casey, T, Zheng, X, Li, J, et al. The relationships between violence in childhood and educational outcomes: a global systematic review and meta-analysis. Child Abuse Negl 2018; 75: 628.10.1016/j.chiabu.2017.06.021CrossRefGoogle ScholarPubMed
Schaefer, JD, Cheng, TW, Dunn, EC. Sensitive periods in development and risk for psychiatric disorders and related endpoints: a systematic review of child maltreatment findings. Lancet Psychiatry 2022; 9: 978–91.10.1016/S2215-0366(22)00362-5CrossRefGoogle ScholarPubMed
Ahn, K, Norman, LJ, Justice, CM, Shaw, P. ADHD and its neurocognitive substrates: a two sample Mendelian randomization study. Transl Psychiatry 2022; 12: 378.10.1038/s41398-022-02139-xCrossRefGoogle ScholarPubMed
Burton, SMI, Sallis, HM, Hatoum, AS, Munafò, MR, Reed, ZE. Is there a causal relationship between executive function and liability to mental health and substance use? A Mendelian randomization approach. R Soc Open Sci 2022; 9: 220631.10.1098/rsos.220631CrossRefGoogle Scholar
Demange, PA, Boomsma, DI, van Bergen, E, Nivard, MG. Evaluating the causal relationship between educational attainment and mental health. medRxiv [Preprint] 2024. Available from: https://www.medrxiv.org/content/10.1101/2023.01.26.23285029v2 [cited 10 Jun 2024].Google ScholarPubMed
Treur, JL, Demontis, D, Smith, GD, Sallis, H, Richardson, TG, Wiers, RW, et al. Investigating causality between liability to ADHD and substance use, and liability to substance use and ADHD risk, using Mendelian randomization. Addict Biol 2021; 26: e12849.10.1111/adb.12849CrossRefGoogle ScholarPubMed
Letkiewicz, AM, Funkhouser, CJ, Shankman, SA. Childhood maltreatment predicts poorer executive functioning in adulthood beyond symptoms of internalizing psychopathology. Child Abuse Negl 2021; 118: 105140.10.1016/j.chiabu.2021.105140CrossRefGoogle ScholarPubMed
Masson, M, Bussières, E-L, East-Richard, C, R-Mercier, A, Cellard, C. Neuropsychological profile of children, adolescents and adults experiencing maltreatment: a meta-analysis. Clin Neuropsychol 2015; 29: 573–94.10.1080/13854046.2015.1061057CrossRefGoogle ScholarPubMed
Kasai, K, Koike, S, Matsuoka, J, Okada, N, Sakakibara, E, Sakurada, H, et al. Shared functional impairment in the prefrontal cortex affects symptom severity across psychiatric disorders. Psychol Med 2022; 52: 2661–70.Google Scholar
Puetz, VB, McCrory, E. Exploring the relationship between childhood maltreatment and addiction: a review of the neurocognitive evidence. Curr Addict Rep 2015; 2: 318–25.10.1007/s40429-015-0073-8CrossRefGoogle ScholarPubMed
Rock, PL, Roiser, JP, Riedel, WJ, Blackwell, AD. Cognitive impairment in depression: a systematic review and meta-analysis. Psychol Med 2014; 44: 2029–40.10.1017/S0033291713002535CrossRefGoogle ScholarPubMed
Romano, E, Babchishin, L, Marquis, R. Childhood maltreatment and educational outcomes. Trauma Violence Abuse 2015; 16: 418–37.10.1177/1524838014537908CrossRefGoogle ScholarPubMed
Slade, EP, Wissow, LS. The influence of childhood maltreatment on adolescents’ academic performance. Econ Educ Rev 2007; 26: 604–14.10.1016/j.econedurev.2006.10.003CrossRefGoogle ScholarPubMed
Schimmenti, A, Billieux, J, Santoro, G, Casale, S, Starcevic, V. A trauma model of substance use: elaboration and preliminary validation. Addict Behav 2022; 134: 107431.10.1016/j.addbeh.2022.107431CrossRefGoogle ScholarPubMed
Brady, DJ, Phalen, PL, Roche, DJO, Cowan, T, Bennett, ME. A reduction in cigarette smoking improves health-related quality of life and does not worsen psychiatric symptoms in individuals with serious mental illness. Addict Behav 2024; 151: 107949.10.1016/j.addbeh.2023.107949CrossRefGoogle Scholar
Hökby S, Westerlund, J, Alvarsson, J, Carli, V, Hadlaczky, G. Longitudinal effects of screen time on depressive symptoms among Swedish adolescents: the moderating and mediating role of coping engagement behavior. Int J Environ Res Public Health 2023; 3771, 20.Google Scholar
Tang, S, Werner-Seidler, A, Torok, M, Mackinnon, AJ, Christensen, H. The relationship between screen time and mental health in young people: a systematic review of longitudinal studies. Clin Psychol Rev 2021; 86: 102021.10.1016/j.cpr.2021.102021CrossRefGoogle Scholar
Metaxas, A, Kempf, SJ. Neurofibrillary tangles in Alzheimer′s disease: elucidation of the molecular mechanism by immunohistochemistry and tau protein phospho-proteomics. Neural Regen Res 2016; 11: 1579–81.Google ScholarPubMed
Coupland, KG, Mellick, GD, Silburn, PA, Mather, K, Armstrong, NJ, Sachdev, PS, et al. DNA methylation of the MAPT gene in Parkinson’s disease cohorts and modulation by vitamin E in vitro. Mov Disord 2014; 29: 1606–14.10.1002/mds.25784CrossRefGoogle ScholarPubMed
Kavanagh, T, Halder, A, Drummond, E. Tau interactome and RNA binding proteins in neurodegenerative diseases. Mol Neurodegen 2022; 17: 66.10.1186/s13024-022-00572-6CrossRefGoogle ScholarPubMed
Yuste-Checa, P, Trinkaus, VA, Riera-Tur, I, Imamoglu, R, Schaller, TF, Wang, H, et al. The extracellular chaperone Clusterin enhances Tau aggregate seeding in a cellular model. Nat Commun 2021; 12: 4863.10.1038/s41467-021-25060-1CrossRefGoogle Scholar
Hall, JR, Petersen, M, Johnson, L. Plasma total tau and neurobehavioral symptoms of cognitive decline in cognitively normal older adults. Front Psychol 2021; 774049, 12.10.3389/fpsyg.2021.774049CrossRefGoogle ScholarPubMed
Saraçoğlu H, Kılıç E, Demirci, E. The study of tau and phospho tau protein levels in attention deficit and hyperactivity disorder. Turk J Med Sci 2021; 51: 2107–11.Google Scholar
Demirel, ÖF, Cetin, I, Turan, Ş, Yıldız, N, Sağlam, T, Duran, A. Total tau and phosphorylated tau protein serum levels in patients with schizophrenia compared with controls. Psychiatr Q 2017; 88: 921–8.10.1007/s11126-017-9507-xCrossRefGoogle ScholarPubMed
Chen, Y, Guan, W, Wang, ML, Lin, XY. PI3K-AKT/mTOR signaling in psychiatric disorders: a valuable target to stimulate or suppress? Int J Neuropsychopharmacol 2024; 27: pyae010.10.1093/ijnp/pyae010CrossRefGoogle ScholarPubMed
Sanacora, G, Yan, Z, Popoli, M. The stressed synapse 2.0: pathophysiological mechanisms in stress-related neuropsychiatric disorders. Nat Rev Neurosci 2022; 23: 86103.10.1038/s41583-021-00540-xCrossRefGoogle ScholarPubMed
Clark, H, Coll-Seck, AM, Banerjee, A, Peterson, S, Dalglish, SL, Ameratunga, S, et al. A future for the worlds children? A WHO–UNICEF–Lancet Commission. Lancet 2020; 395: 605–58.10.1016/S0140-6736(19)32540-1CrossRefGoogle ScholarPubMed
Malvaso, C, Pilkington, R, Montgomerie, A, Delfabbro, P, Lynch, J. A public health approach to preventing child maltreatment: an intelligent information infrastructure to help us know what works. Child Abuse Negl 2020; 106: 104466–9.10.1016/j.chiabu.2020.104466CrossRefGoogle ScholarPubMed
Substance Abuse and Mental Health Services Administration (US). Trauma-informed care: a sociocultural perspective. In Trauma-Informed Care in Behavioral Health Services. Substance Abuse and Mental Health Services Administration (US), 2014.Google Scholar
Bailey, S, Newton, N, Perry, Y, Grummitt, L, Baams, L, Barrett, E. Trauma-informed prevention programmes for depression, anxiety, and substance use among young people: protocol for a mixed-methods systematic review. Syst Rev 2023; 12: 203.10.1186/s13643-023-02365-4CrossRefGoogle ScholarPubMed
Shin, SH. Preventing e-cigarette use among high-risk adolescents: a trauma-informed prevention approach. Addict Behav 2021; 115: 106795.10.1016/j.addbeh.2020.106795CrossRefGoogle ScholarPubMed
Rhoades, KA, Mitnick, DM, Heyman, RE, Slep, AMS, Del Vecchio, T. Systemic approaches to child maltreatment. In The Handbook of Systemic Family: Therapy (eds Wampler, KS, McWey, LM): 163–89. Wiley Blackwell, 2020.Google Scholar
Espinosa Dice, AL, Lawn, RB, Ratanatharathorn, A, Roberts, AL, Denckla, CA, Kim, AH, et al. Childhood maltreatment and health in the UK Biobank: triangulation of outcome-wide and polygenic risk score analyses. BMC Med 2024; 22: 135.10.1186/s12916-024-03360-9CrossRefGoogle ScholarPubMed
ter Kuile, AR, Hübel, C, Cheesman, R, Coleman, JRI, Peel, AJ, Levey, DF, et al. Genetic decomposition of the heritable component of reported childhood maltreatment. Biol Psychiatry Glob Open Sci 2023; 3: 716–24.10.1016/j.bpsgos.2023.03.003CrossRefGoogle ScholarPubMed
Ng, NYT. Modifiable Pathways and Epigenetic Markers for the Effects of Childhood Maltreatment on Neuropsychiatric Disorders: A Mendelian Randomization Study. University of Hong Kong, 2024 (https://hub.hku.hk/handle/10722/352702 [cited 23 Dec 2024]).Google Scholar
Figure 0

Fig. 1 Overall study design with the aims, data-sets, statistical methods and expected outcomes. CM, childhood maltreatment; EWAS, epigenome-wide association study; GWAS, genome-wide association study; PGC, Psychiatric Genomics Consortium; ADHD, attention-deficit/hyperactivity disorder; MR, Mendelian randomisation; MR-PRESSO, Mendelian randomisation pleiotropy residual sum and outlier; MR-Egger, Mendelian randomisation-Egger; IL6, interleukin-6; IL6R, Interleukin-6 receptor; CRP, C-reactive protein; BMI, body mass index; IVW, inverse-variance weighted; MVMR, multivariable Mendelian randomisation; DNAm, DNA methylation; mQTL, methylation quantitative trait loci; GoDMC, Genetics of DNA Methylation Consortium.

Figure 1

Table 1 Indirect effects of significant mediators for childhood maltreatment on neuropsychiatric outcome relationship

Figure 2

Fig. 2 Direct effects of childhood maltreatment (CM) adjusting for individual or all significant mediators. ADHD, attention-deficit/hyperactivity disorder; MDD, major depressive disorder; SNPs, single-nucleotide polymorphisms; OR, odds ratio.

Figure 3

Fig. 3 Breakdown of the mediators of the relationship between childhood maltreatment (CM) and neuropsychiatric disorders, potential mechanisms and suggested interventions. Mediators in red are addiction-related traits, mediators in green are cognitive-related traits and the mediator in blue is a socioeconomic trait. Solid arrows represent causal links and the dashed arrow represents associations, not causality. DNAm, DNA methylation; ADHD, attention-deficit/hyperactivity disorder; MDD, major depressive disorder; OR, odds ratio.

Supplementary material: File

Ng et al. supplementary material 1

Ng et al. supplementary material
Download Ng et al. supplementary material 1(File)
File 694.2 KB
Supplementary material: File

Ng et al. supplementary material 2

Ng et al. supplementary material
Download Ng et al. supplementary material 2(File)
File 135.5 KB
Submit a response

eLetters

No eLetters have been published for this article.