Introduction
Epidemiological studies have consistently shown an association between cannabis use and the emergence of psychosis (Di Forti et al., Reference Di Forti, Quattrone, Freeman, Tripoli, Gayer-Anderson and Quigley2019a; Reference Di Forti, Wu-Choi, Quattrone, Richards, Freeman, Tripoli and Murray2019b; Pries et al., Reference Pries, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2018). Mendelian randomization studies examining the association between genetic predisposition, cannabis use and psychosis have found a bidirectional relationship i.e. genetic vulnerability predisposes to earlier cannabis use and cannabis use predisposes to psychosis among those with genetic vulnerability (Gillespie & Kendler, Reference Gillespie and Kendler2021). Patients with psychosis also have high rates of cannabis use (Boydell et al., Reference Boydell, van Os, Caspi, Kennedy, Giouroukou, Fearon and Murray2006). The reasons for these seemingly counterintuitive findings of greater risk of psychosis with cannabis use and high rates of cannabis use among patients with psychosis are not fully understood. One possible explanation for these discrepant findings is that cannabis use ameliorates anxiety and depressive symptoms while also increasing the risk for psychotic experiences (PE).
While it is difficult to test this hypothesis in patients who have already developed psychosis, a general population, longitudinal sample such as the Netherlands Mental Health Survey and Incidence Study-2 (NEMESIS-2), provides a unique opportunity to examine this question. Approximately 5–8% of the general population may report PE (Linscott & van Os, Reference Linscott and van Os2013). We have previously shown that anxiety and depressive symptoms increase the risk of PE in interaction with environmental risk factors such as cannabis use and childhood adversity, consistent with an ‘affective pathway’ to psychosis (Guloksuz et al., Reference Guloksuz, van Nierop, Lieb, van Winkel, Wittchen and van Os2015; Pries et al., Reference Pries, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2018; Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019; van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020). The important role of an ‘affective pathway’ to psychosis is also supported by experimental and observational studies (Bird, Waite, Rowsell, Fergusson, & Freeman, Reference Bird, Waite, Rowsell, Fergusson and Freeman2017; Freeman et al., Reference Freeman, Dunn, Fowler, Bebbington, Kuipers, Emsley and Garety2013a, Reference Freeman, Startup, Dunn, Cernis, Wingham, Pugh and Kingdon2013b; Garety, Kuipers, Fowler, Freeman, & Bebbington, Reference Garety, Kuipers, Fowler, Freeman and Bebbington2001; Hanssen, Bak, Bijl, Vollebergh, & van Os, Reference Hanssen, Bak, Bijl, Vollebergh and van Os2005) where anxiety, depressive symptoms and worry is associated with paranoia and anomalous perceptual experiences, as well as experience sampling studies showing that increased stress-reactivity and negative affect are associated with subsequent increase in PE (Kramer et al., Reference Kramer, Simons, Wigman, Collip, Jacobs, Derom and Wichers2014; Myin-Germeys & van Os, Reference Myin-Germeys and van Os2007). A network analysis of psychopathology showed that affective symptoms were on the pathway between environmental risk-factors such as cannabis use and psychosis expression (Isvoranu et al., Reference Isvoranu, van Borkulo, Boyette, Wigman, Vinkers and Borsboom2017). More recently, depression was shown to contribute to greater than multiplicative effect to the association between polygenic risk of schizophrenia, psychosis and childhood trauma (van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020). Additionally, it has been shown recently in the NEMESIS-2 data, that a fixed-effects model yields reliable estimates when examining the causal association between cannabis use and psychosis, and that the direction of causality is from cannabis use to psychosis i.e. cannabis use precedes the onset of psychosis (van Os et al., Reference van Os, Pries, ten Have, de Graaf, van Dorsselaer, Bak and Guloksuz2021a, Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Bak and Guloksuz2021b). The question then arises whether the causal association between cannabis use and PE is mediated by anxiety/depressive symptoms or whether the relationship is bidirectional.
The presence of a bidirectional relationship between cannabis use and anxiety/depressive symptoms in the relationship to PE can be examined using mediation analysis. In this regard, we sought to examine the mediational relationship between cannabis use, anxiety/depressive symptoms, and PE. We hypothesized that cannabis use mediates the relationship between anxiety/depressive symptoms and PE. We also examined whether anxiety/depressive symptoms mediated the relationship between cannabis use and PE.
Method
Study population
The Netherlands Mental Health Survey and Incidence Study-2 (NEMESIS-2) is a longitudinal cohort study comprising 4 waves that examines the prevalence, incidence, course, and consequences of psychiatric disorders in the Dutch general population. The first wave (T0) that was conducted from 2007 to 2009, enrolled 6646 participants (response rate 65.1%; average interview duration: 95 min). They were followed up at year 3 (T1), year 6 (T2), and year 9 (T3) with successive response rates of 80.4% (n = 5303; excluding those who deceased; interview duration: 84 min), 87.8% (n = 4618; interview duration: 83 min), and 87.7% (n = 4007; interview duration: 101 min), respectively. The study was approved by the Medical Ethics Review Committee for Institutions on Mental Health Care and written informed consent was obtained from participants at each wave. A multistage random sampling procedure was applied to ensure representativeness of the sample to the general population in terms of age (between the ages of 18 and 65 at baseline), region, and population density. All participants were required to be literate in Dutch. The Composite International Diagnostic Interview (CIDI) version 3.0 [de Graaf, Ormel, ten Have, Burger, & Buist-Bouwman, Reference de Graaf, Ormel, ten Have, Burger and Buist-Bouwman2008; Alonso J, et al., 2004] and additional questionnaires were administered by non-clinician, trained interviewers during home visits. Additional details of NEMESIS-2 are provided elsewhere (de Graaf, Ten Have, & van Dorsselaer, Reference de Graaf, Ten Have and van Dorsselaer2010; de Graaf, ten Have, van Gool, & van Dorsselaer, Reference de Graaf, ten Have, van Gool and van Dorsselaer2012). Data from all four waves were utilized. While there was some attrition between T0 and T3, this was not significantly associated with any of the individual 12-month mental disorders at T0 after controlling for sociodemographic characteristics (Nuyen et al., Reference Nuyen, Tuithof, de Graaf, van Dorsselaer, Kleinjan and Have2020).
Outcome measure: psychotic experiences (PE)
PE were assessed using a 20-item binary-response questionnaire based on CIDI 1.1 as used elsewhere (Pries et al., Reference Pries, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2018; Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019; van Nierop et al., Reference van Nierop, Viechtbauer, Gunther, van Zelst, de Graaf, Ten Have and van Winkel2015), since earlier CIDI versions were found not to be adequate for capturing positive psychotic symptomatology. The primary binary outcome was the presence of any self-report PE conforming to previous studies (Guloksuz et al., Reference Guloksuz, Pries, ten Have, de Graaf, van Dorsselaer, Klingenberg and van Os2020; Pries et al., Reference Pries, Erzin, van Os, ten Have, de Graaf, van Dorsselaer and Guloksuz2021, Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020; Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019).
Assessment of anxiety and depressive symptoms
Symptoms of anxiety disorders including social phobia, specific phobia, panic disorder, generalized anxiety disorder, agoraphobia without panic disorder and symptoms of depressive disorders were assessed using CIDI 3.0 at each visit (de Graaf et al., Reference de Graaf, Ormel, ten Have, Burger and Buist-Bouwman2008). Conforming to previous publications, we used two binary variables for depressive and anxiety symptoms, respectively (e.g. considered present if participants experienced at least one of the CIDI 3.0 core symptoms of depressive (or anxiety) disorder at each time-point) (Pries et al., Reference Pries, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2018; Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019).
Assessment of cannabis use
The frequency of cannabis use over the last year was assessed using the section Illegal Substance Use of the CIDI 3.0 at each visit as used previously (Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019). If subjects reported cannabis use, they were rated on frequency of use on a scale of 1 (never) to 7 (every day). A binary variable (absent = ‘0’ and present = ‘1’) was constructed by using the cut-off value of once per week or more in the previous year, consistent with other studies (Pries et al., Reference Pries, Erzin, van Os, ten Have, de Graaf, van Dorsselaer and Guloksuz2021, Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020; Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019).
Statistical analysis
All analyses were performed using Stata, version 16 (StataCorp., 2019). The level of significance (α) was set at 0.05. To establish mediation, the following assumptions should be met [detailed in Baron and Kenny (Reference Baron and Kenny1986), Gelfand, Mensinger, and Tenhave (Reference Gelfand, Mensinger and Tenhave2009)]: (1) The independent variable should be associated with the outcome; (2) There should be a association between the independent variable and the mediator, and between the mediator and the outcome; (3) The independent variable, mediated through the mediator, should be associated with the outcome. For mediation, the total effect of the independent variable on the outcome should be, at least partly, explained by the mediator (indirect effect). In the mediation model, the direct effect of the independent variable should decrease.
Given that mediation requires temporal precedence from independent variable to mediator to outcome, thus requiring absence of the outcome at the previous wave T–1, PE incident at time point T was modeled as a function of independent variable and mediator assessed over the interval T–1 to T, conforming to previous analyses conducted in the NEMESIS-2 (van Os et al., 2021a, Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Bak and Guloksuz2021b).
Logistic regression models were performed to analyze whether there was an association between the presence of affective dysregulation (any anxiety symptom and any depressive symptom, separately) at the preceding time point (T–1) and PE outcomes. Consequently, the KHB procedure, developed for estimating mediation in logistic regression models, was applied to estimate the mediation effect of cannabis use (Fig. 1a). Logistic regression models were performed to analyze whether there was an association between cannabis use in the previous year and PE outcomes. Subsequently, the KHB procedure was applied to estimate the mediation effect of the presence of affective dysregulation (any anxiety symptom and any depressive symptom, separately) (Fig. 1b). To correct for clustering of multiple observations within subjects, the CLUSTER option was applied in all logistic regression models to estimate cluster-robust standard errors.
All analyses were corrected for sex, age, education (1) primary school, (2) lower secondary education, (3) higher secondary education, (4) higher professional education).
Results
This analysis used data from 6646 participants at four time points (T0, T1, T2, and T3) [n = 6646 at baseline (T0), n = 5303 in the first wave (T1), n = 4618 in the second wave (T2) and n = 4007 in the third wave (T3)]. Table 1 reports the demographics of the NEMESIS-2 at each time point.
Data are given as number (percentage) unless otherwise indicated.
a Baseline covers lifetime, and each follow-up covers the preceding period.
Cannabis use as mediator for the association between anxiety symptoms and psychotic experience
Preceding anxiety symptoms were associated with cannabis use (OR 2.55, 95% CI 1.78–3.66, p < 0.001) and any self-reported PE (OR 1.99, 95% CI 1.65–2.42, p < 0.001); cannabis use was associated with self-reported PE (OR 2.36, 95% CI 1.25–4.46, p = 0.008). The mediation model revealed that cannabis use mediated the association between preceding anxiety symptoms and self-reported PE but explained only 1% of the relationship (Total effect of anxiety OR 1.97, 95% CI 1.63–2.39, p < 0.001; direct effect of anxiety OR 1.96, 95% CI 1.62–2.37, p < 0.001; indirect effect attributable to cannabis use OR 1.01, 95% CI 1.00–1.01, p = 0.045; % of total effect attributable to cannabis use = 1%) (Table 2).
Anxiety symptoms as mediator for the association between cannabis use and psychotic experience
Preceding cannabis use was associated with anxiety (OR 1.55, 95% CI 1.00–2.39, p = 0.05) and, directionally and with greater effect size, any self-reported PE (OR 1.77, 95% CI 0.76–4.12, p = 0.187); anxiety was associated with self-reported PE (OR 3.05, 95% CI 2.44–3.82, p < 0.001). The mediation model revealed that anxiety explained 17% of the total relationship (Total effect of cannabis OR 1.75, 95% CI 0.75–4.06, p = 0.195; direct effect of cannabis OR 1.56, 95% CI 0.68–3.71, p = 0.281; indirect effect attributable to anxiety OR 1.10, 95% CI 0.99–1.12, p = 0.065; % of total effect attributable to anxiety = 17%) (Table 2).
Cannabis use as the mediator for the association between depressive symptoms and psychotic experience
Preceding depression was associated with cannabis use (OR 3.39, 95% CI 2.37–4.86, p < 0.001) and self-reported PE (OR 1.89, 95% CI 1.56–2.28, p < 0.001); cannabis use was associated with self-reported PE (OR 2.02, 95% CI 1.06–3.85, p = 0.033). The mediation model revealed that cannabis use was found to be a significant mediator of the relationship between depression and self-reported PE, but explained only 1.4% of the relationship (total effect of depression OR 1.86, 95% CI 1.54–2.25, p < 0.001; direct effect of depression OR 1.84, 95% CI 1.52–2.23, p < 0.001; indirect effect attributable to cannabis use OR 1.01, 95% CI 1.00–1.02, p = 0.045; % of total effect attributable to cannabis use = 1.4%) (Table 3).
Depressive symptoms as the mediator for the association between cannabis use and psychotic experiences
Preceding cannabis use was associated with depressive symptoms (OR 2.43, 95% CI 1.58–3.74, p < 0.001) and self-reported PE at a trend level (OR 1.84, 95% CI 0.86–3.98, p = 0.118); depressive symptoms was associated with self-reported PE (OR 2.86, 95% CI 2.27–3.60, p < 0.001). The mediation model revealed that depressive symptoms were found to be a significant mediator of the relationship between cannabis use and self-reported PE, and explained 37% of the relationship (Total effect of cannabis use OR 1.77, 95% CI 0.82–3.85, p = 0.146; direct effect of cannabis use OR 1.44, 95% CI 0.66–3.13, p = 0.359; indirect effect attributable to depressive symptoms OR 1.23, 95% CI 1.11–1.37, p < 0.001; % of total effect attributable to depressive symptoms = 37%) (Table 3).
Discussion
In this analysis of the NEMESIS-2 dataset, a large longitudinal cohort among the general population, we examined the mediational relationship between cannabis use, anxiety/depression symptoms (i.e. affective dysregulation) and PE. Our analysis showed that the mediational relationship between cannabis use, anxiety/depressive symptoms (affective dysregulation) and PE was bidirectional i.e. anxiety/depressive symptoms mediated the relationship between cannabis use and PE; and cannabis use mediated the relationship between anxiety/depressive symptoms and PE. However, although the relationship was bidirectional, the contribution of affective dysregulation in mediating the relationship between cannabis use and PE over time was greater than the mediational contribution of cannabis. These results lend further support for the presence of an ‘affective pathway’ to psychosis (Guloksuz et al., Reference Guloksuz, van Nierop, Lieb, van Winkel, Wittchen and van Os2015; Pries et al., Reference Pries, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2018; Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019; van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020).
This is the first study to our knowledge, to show a bidirectional mediational relationship between cannabis use, affective dysregulation and psychosis in a general population sample. The directionality of increased anxiety mediating the association between cannabis use and psychosis is consistent with human laboratory challenge studies using delta-9 tetrahydrocannabinol (THC), the primary psychoactive ingredient, which show that administration of THC results in a significant increase in anxiety and psychosis-like experience (Ganesh et al., Reference Ganesh, Cortes-Briones, Ranganathan, Radhakrishnan, Skosnik and D'Souza2020; Hindley et al., Reference Hindley, Beck, Borgan, Ginestet, McCutcheon, Kleinloog and Howes2020). Interestingly, preceding anxiety (pharmacologically induced using the GABA inverse agonist, Iomazenil) has also been shown to increase the psychosis-like effects of THC (Radhakrishnan et al., Reference Radhakrishnan, Skosnik, Cortes-Briones, Sewell, Carbuto, Schnakenberg and D'Souza2015). In epidemiological studies, cannabis use has been associated with greater depression in the general population (Rabiee et al., Reference Rabiee, Lundin, Agardh, Hensing, Allebeck and Danielsson2020), and among first-episode psychosis (Elowe, Golay, Baumann, Solida-Tozzi, & Conus, Reference Elowe, Golay, Baumann, Solida-Tozzi and Conus2020), consistent with the directionality of cannabis use resulting in greater depression in our analysis. Animal studies show that THC has biphasic effects, i.e. low-dose THC is anxiolytic while high-dose THC is anxiogenic al (Salviato et al., Reference Salviato, Raymundi, Rodrigues Da Silva, Salemme, Batista Sohn, Araujo and Stern2021). While this has not been shown in human studies, an anxiolytic and antidepressant response to cannabis at low doses may explain the directionality of preceding anxiety, depressive symptoms resulting in increased cannabis use. These results are also consistent with the seemingly contradictory clinical observation that cannabis use is associated with psychosis, and people with psychosis are also at high risk of cannabis use.
However, as noted earlier, the contribution of cannabis use as a mediator of the relationship between affective dysregulation and PE was small in our analysis. These results hence do not support a self-medication hypothesis, but rather suggest that cannabis use results in increased affective dysregulation, supporting an ‘affective pathway’ to psychosis (Radhakrishnan et al., Reference Radhakrishnan, Guloksuz, Ten Have, de Graaf, van Dorsselaer, Gunther and van Os2019; van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020).
Our results are consistent with studies examining the interactive effects of psychosis risk using schizophrenia polygenic risk score (PRS-SZ) for psychosis, cannabis use and affective dysregulation (Di Forti et al., Reference Di Forti, Quattrone, Freeman, Tripoli, Gayer-Anderson and Quigley2019a, Reference Di Forti, Wu-Choi, Quattrone, Richards, Freeman, Tripoli and Murray2019b; Guloksuz et al., Reference Guloksuz, Pries, Delespaul, Kenis, Luykx, Lin and van Os2019; van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020). van Os et al., recently found that the relationship between PRS-SZ and psychosis showed significant dependence on the concomitant presence of affective dysregulation in 2 independent cohorts, NEMESIS-2 and EUGEI (van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020). This relationship was particularly evident for delusional ideation but not for hallucinatory experiences (van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020). A similar relationship was also with childhood trauma, such that the relationship between PRS-SZ and childhood trauma showed significant dependence on the concomitant presence of affective dysregulation, and was particularly evident for delusional ideation (van Os et al., Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Delespaul and Guloksuz2020). Cannabis-induced dysphoria was also shown to mediate the relationship between childhood trauma and psychotic experiences (Carlyle et al., Reference Carlyle, Constable, Walter, Wilson, Newland and Hides2021).
The common biological mechanism that may explain this relationship remains to be elucidated. Interestingly, cannabis use, childhood trauma and affective dysregulation result in impairment in hippocampal structure and function (Barch et al., Reference Barch, Shirtcliff, Elsayed, Whalen, Gilbert, Vogel and Luby2020; D'Souza et al., Reference D'Souza, Radhakrishnan, Naganawa, Ganesh, Nabulsi, Najafzadeh and Skosnik2020; Treadway et al., Reference Treadway, Waskom, Dillon, Holmes, Park, Chakravarty and Pizzagalli2015; Yuan et al., Reference Yuan, Rubin-Falcone, Lin, Rizk, Miller, Sublette and Mann2020; Zeredo et al., Reference Zeredo, Quah, Wallis, Alexander, Cockcroft, Santangelo and Clarke2019). Abnormalities in hippocampal microstructure have also been seen in schizophrenia using PET imaging (Radhakrishnan et al., Reference Radhakrishnan, Skosnik, Ranganathan, Naganawa, Toyonaga, Finnema and D'Souza2021). One hypothesis is that hippocampal dysfunction, via the hippocampal-striatal-midbrain network, leads to aberrant signaling of salience (Modinos et al., Reference Modinos, Allen, Zugman, Dima, Azis, Samson and McGuire2020) and hence, may represent a common biological mechanism leading to delusion formation, especially in those with affective dysregulation (Rauschenberg et al., Reference Rauschenberg, Reininghaus, Ten Have, de Graaf, van Dorsselaer, Simons and van Os2021).
Strengths and limitations
These findings should be interpreted in the light of several strengths and limitations. In terms of strengths, the present study has the advantage of a large number of participants representative of the Dutch population (Linscott & van Os, Reference Linscott and van Os2010; van Nierop et al., Reference van Nierop, van Os, Gunther, Myin-Germeys, de Graaf, ten Have and van Winkel2012). The analysis was conducted among 6646 participants at four time points. The study incorporated thorough symptom assessments across diagnostic boundaries providing an opportunity to examine trans-syndromal psychopathology in the general population. One limitation of the study is that the assessment of cannabis use was based on a questionnaire and might be subject to recall-bias and under-reporting. Similarly, the assessment of anxiety and depressive symptoms were based on CIDI 3.0 rather than clinical scales that specifically assess the severity of anxiety and depression, such as. Additionally, the current analysis could not answer the question of whether within-subject effects were different than between-subject effects as noted earlier (van Os et al., Reference van Os, Pries, ten Have, de Graaf, van Dorsselaer, Bak and Guloksuz2021a, Reference van Os, Pries, Ten Have, de Graaf, van Dorsselaer, Bak and Guloksuz2021b). Another limitation is that the average age of participants at baseline in our sample was 44.3 years. This can in part account for the relatively lower contribution of cannabis use on the outcome compared to that of anxiety and depression. The impact of cannabis use on psychosis risk is greater in adolescence and it decreases over time. The information regarding the potency of cannabis was also not available in this sample. Additionally, cannabis use was coded as a dichotomous variable which limits the range of this measure in identifying whether more frequent use has a higher impact on psychosis.
Conclusion
In this large, general population sample, we found that the mediational relationship between cannabis use, anxiety/depressive symptoms (affective dysregulation) and PE is bidirectional. The contribution of affective dysregulation in mediating the relationship between cannabis use and PE was far greater than the mediational contribution of cannabis. This study provides further support for the existence of an ‘affective pathway’ to psychosis. Whether interventions addressing affective dysregulation among cannabis-using individuals can reduce the risk of the emergence of PE remains to be verified.
Financial support
This work was supported by the Ministry of Health, Welfare and Sport (Grant Number 310253), with supplement support from the Netherlands Organization for Health Research and Development (ZonMw) and the Genetic Risk and Outcome of Psychosis investigators. Lotta Pries is supported by the Kootstra Talent Fellowship of Maastricht University. Jim van Os and Sinan Guloksuz are supported by the Ophelia research project, ZonMw grant number: 636340001. Bart Rutten is supported by a VIDI award, grant number 91718336 from the Netherlands Scientific Organisation.
Conflict of interest
Rajiv Radhakrishnan, Sinan Guloksuz, Lotta-Katrin Pries, Saskia van Dorsselaer, Margreet ten Have, Ron de Graaf, Jim van Os, and Maarten Bak do not have any conflicts of interest and have received no funding or compensation from companies. Rajiv Radhakrishnan is supported by Dana Foundation David Mahoney program, NIMH (grants R21 MH115316-02 and R21 MH123870-01A1), CTSA Grant Number UL1 TR001863 from the National Center for Advancing Translational Science, components of the National Institutes of Health (NIH), and NIH roadmap for Medical Research. The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official view of NIH.