Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-05T14:55:47.874Z Has data issue: false hasContentIssue false

Development and initial validation of a multivariable predictive Early Adversity Scale for Schizophrenia (EAS-Sz) using register data to quantify environmental risk for adult schizophrenia diagnosis after childhood exposure to adversity

Published online by Cambridge University Press:  12 July 2022

Patsy Di Prinzio*
Affiliation:
Neuropsychiatric Epidemiology Research Unit, School of Population and Global Health, The University of Western Australia, Perth, Australia
Jonas Björk
Affiliation:
Division of Occupational and Environmental Medicine, Lund University, Lund, Sweden
Giulietta Valuri
Affiliation:
School of Population and Global Health, The University of Western Australia, Perth, Australia
Taryn Ambrosi
Affiliation:
Neuropsychiatric Epidemiology Research Unit, School of Population and Global Health, The University of Western Australia, Perth, Australia
Maxine Croft
Affiliation:
Neuropsychiatric Epidemiology Research Unit, School of Population and Global Health, The University of Western Australia, Perth, Australia
Vera A. Morgan
Affiliation:
Neuropsychiatric Epidemiology Research Unit, School of Population and Global Health, The University of Western Australia, Perth, Australia
*
Author for correspondence: Patsy Di Prinzio, E-mail: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Background

Additional to a child's genetic inheritance, environmental exposures are associated with schizophrenia. Many are broadly described as childhood adversity; modelling the combined impact of these is complex. We aimed to develop and validate a scale on childhood adversity, independent of genetic and other environmental liabilities, for use in schizophrenia risk analysis models, using data from cross-linked electronic health and social services registers.

Method

A cohort of N = 428 970 Western Australian children born 1980–2001 was partitioned into three samples: scale development sample (N = 171 588), and two scale validation samples (each N = 128 691). Measures of adversity were defined before a child's 10th birthday from five domains: discontinuity in parenting, family functioning, family structure, area-level socioeconomic/demographic environment and family-level sociodemographic status. Using Cox proportional hazards modelling of follow-up time from 10th birthday to schizophrenia diagnosis or censorship, weighted combinations of measures were firstly developed into scales for each domain, then combined into a final global scale. Discrimination and calibration performance were validated using Harrell's C and graphical assessment respectively.

Results

A weighted combination of 42 measures of childhood adversity was derived from the development sample. Independent application to identical measures in validation samples produced Harrell's Concordance statistics of 0.656 and 0.624. Average predicted time to diagnosis curves corresponded with 95% CI limits of observed Kaplan–Meier curves in five prognostic categories.

Conclusions

Our Early Adversity Scale for Schizophrenia (EAS-Sz), the first using routinely collected register data, predicts schizophrenia diagnosis above chance, and has potential to help untangle contributions of genetic and environmental liability to schizophrenia risk.

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 (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2022. Published by Cambridge University Press

Introduction

A child's genetic inheritance contributes to their risk of developing schizophrenia, a highly polygenic disorder, with heritability estimated at 60–80% (Owen, Sawa, & Mortensen, Reference Owen, Sawa and Mortensen2016). Lifetime morbid risk ranges from approximately 1% within the general population to a concordance of 48% for monozygotic twins where one twin is affected (Gottesman, Reference Gottesman1994; Gottesman & Erlenmeyer-Kimling, Reference Gottesman and Erlenmeyer-Kimling2001; Gottesman, McGuffin, & Farmer, Reference Gottesman, McGuffin and Farmer1987). Environmental exposures are also implicated in schizophrenia risk, either exclusively or in interaction with genetic risk (Morgan & Gayer-Anderson, Reference Morgan and Gayer-Anderson2016; van Os, Kenis, & Rutten, Reference van Os, Kenis and Rutten2010; Varese et al., Reference Varese, Smeets, Drukker, Lieverse, Lataster, Viechtbauer and Bentall2012). Environmental exposures previously studied include physical trauma such as obstetric complications (Cannon, Jones, & Murray, Reference Cannon, Jones and Murray2002), cannabis use (Marconi, Di Forti, Lewis, Murray, & Vassos, Reference Marconi, Di Forti, Lewis, Murray and Vassos2016) and exposure to a wide range of experiences termed ‘childhood adversity’ (Morgan & Gayer-Anderson, Reference Morgan and Gayer-Anderson2016; van Os et al., Reference van Os, Kenis and Rutten2010; Varese et al., Reference Varese, Smeets, Drukker, Lieverse, Lataster, Viechtbauer and Bentall2012). The latter include discontinuity in parenting either by death or separation (Morgan et al., Reference Morgan, Kirkbride, Leff, Craig, Hutchinson, McKenzie and Fearon2007; Paksarian, Eaton, Mortensen, Merikangas, & Pedersen, Reference Paksarian, Eaton, Mortensen, Merikangas and Pedersen2015a); childhood abuse (Read, van Os, Morrison, & Ross, Reference Read, van Os, Morrison and Ross2005; Varese et al., Reference Varese, Smeets, Drukker, Lieverse, Lataster, Viechtbauer and Bentall2012); social disadvantage (Wicks, Hjern, & Dalman, Reference Wicks, Hjern and Dalman2010); being born or raised in an highly urbanised environment (Vassos, Pedersen, Murray, Collier, & Lewis, Reference Vassos, Pedersen, Murray, Collier and Lewis2012); migration; and minority ethnicity (Cantor-Graae & Pedersen, Reference Cantor-Graae and Pedersen2013; Morgan, Charalambides, Hutchinson, & Murray, Reference Morgan, Charalambides, Hutchinson and Murray2010), including Indigenous status (Tapsell, Hallett, & Mellsop, Reference Tapsell, Hallett and Mellsop2018). Childhood exposure to maladaptive family functioning has been associated with a broad range of psychopathology in adulthood (Kessler et al., Reference Kessler, McLaughlin, Green, Gruber, Sampson, Zaslavsky and Williams2010). While many factors appear to operate at the individual level, neighbourhood (or area) level characteristics such as disadvantage, residential mobility and crime rates have also been found to impact on risk (Veling, Susser, Selten, & Hoek, Reference Veling, Susser, Selten and Hoek2015).

Modelling adversities and their interaction with genetic liability in risk factor analysis is complex. Severity of impact will vary depending on the adversity and may be related to age at exposure (Agid et al., Reference Agid, Shapira, Zislin, Ritsner, Hanin, Murad and Lerer1999; Paksarian et al., Reference Paksarian, Eaton, Mortensen, Merikangas and Pedersen2015a; Paksarian, Eaton, Mortensen, & Pedersen, Reference Paksarian, Eaton, Mortensen and Pedersen2015b). Risk may increase in a dose–response relationship as the number of adversities increases (Morgan et al., Reference Morgan, Kirkbride, Hutchinson, Craig, Morgan, Dazzan and Fearon2008a; Wicks, Hjern, Gunnell, Lewis, & Dalman, Reference Wicks, Hjern, Gunnell, Lewis and Dalman2005; Zammit, Lewis, Dalman, & Allebeck, Reference Zammit, Lewis, Dalman and Allebeck2010) and, because adversities cluster (Dong et al., Reference Dong, Anda, Felitti, Dube, Williamson, Thompson and Giles2004; Kessler et al., Reference Kessler, McLaughlin, Green, Gruber, Sampson, Zaslavsky and Williams2010), a comprehensive measure of early life adversity exposure should account for the interactive effects of co-occurring adversities on the total risk for outcome (Evans, Li, & Whipple, Reference Evans, Li and Whipple2013). To date, we have only simple tools for measuring relative contributions in risk attributable to given sets of childhood adversities, and for advancing our understanding of how they may interact with genetic liability.

A common approach to modelling adversities in the study of schizophrenia has been to include them as separate terms in multivariable explanatory models (Morgan et al., Reference Morgan, Di Prinzio, Valuri, Croft, McNeil and Jablensky2019; Veling et al., Reference Veling, Susser, Selten and Hoek2015) with some studies allowing for the possibility of interactions between risk factors (Zammit et al., Reference Zammit, Lewis, Dalman and Allebeck2010). Some summary exposure scores have also been proposed. One approach is to sum the number of risk exposures (Veling et al., Reference Veling, Susser, Selten and Hoek2015), which implicitly assumes that the risk contribution of each exposure is equal. Other approaches apply weightings reflecting the severity of impact (determined from estimates in literature) to exposures before summing (Padmanabhan, Shah, Tandon, & Keshavan, Reference Padmanabhan, Shah, Tandon and Keshavan2017; Vassos et al., Reference Vassos, Sham, Kempton, Trotta, Stilo, Gayer-Anderson and Morgan2020). More complex summaries that take interactions into account have been developed to predict several more common adult mental disorders, although not schizophrenia (Green et al., Reference Green, McLaughlin, Berglund, Gruber, Sampson, Zaslavsky and Kessler2010; Kessler et al., Reference Kessler, McLaughlin, Green, Gruber, Sampson, Zaslavsky and Williams2010). One prognostic scale recently developed is the Exposome Score for Schizophrenia (Pries et al., Reference Pries, Lage-Castellanos, Delespaul, Kenis, Luykx, Lin and Guloksuz2019). Here, candidate models, incorporating severity weightings of eight exposures (cannabis use, winter birth, hearing impairment, childhood neglect or psychological, physical or sexual abuse, and childhood peer victimisation), were alternately derived via techniques including logistic regression. Testing and assessment conducted in an independent data set (Pries et al., Reference Pries, Erzin, van Os, Ten Have, de Graaf, van Dorsselaer and Guloksuz2020) indicated that models derived via logistic regression predicted schizophrenia at least as well as, or more accurately, than other techniques trialled, including those using machine learning.

Absent in the literature is a comprehensive adversity scale for risk of schizophrenia that considers correlations between variables, clustering of observations, possible interactions between effects, differences in severity of impact and differential impact depending on age at exposure. Critically, there is a need for a scale that can capture adverse exposures in childhood when, according to the neurodevelopmental/developmental perspectives on schizophrenia, it is assumed that the impact of stress on a rapidly developing brain may be greatest (Murray, Bhavsar, Tripoli, & Howes, Reference Murray, Bhavsar, Tripoli and Howes2017). Our aim was to fill this gap by developing and validating a childhood adversity scale covering multiple domains of adversity with the goal of delivering such as a robust tool which could be further applied to risk analysis models. Moreover, because childhood exposures may occur up to two decades or more before schizophrenia onset, a tool was needed that could be applied to existing prospective data registers. The use of prospectively collected data has great advantages over retrospective, possibly biased recall data (Vargas & Mittal, Reference Vargas and Mittal2018). Hence, our focus was on developing a scale for use with cross-linked data from multiple electronic health and social services registers. Finally, to maximise translational capacity, the methodology needed to be transparent, easy to replicate and adaptable to variations in specific register data.

Methods

Cohort identification

We used a whole-population cohort of children (N = 428 970) [henceforth referred to as the e-Cohort (Morgan et al., Reference Morgan, Valuri, Croft, Griffith, Shah, Young and Jablensky2011)] which included all children born alive in Western Australia between 1 January 1980 and 31 December 2001, born to two sets of mothers: (i) those with no recorded contact with mental health services (N = 407 639 children, N = 217 929 mothers); and (ii) those with a diagnosis of psychotic illness (N = 21 331 children, N = 10 605 mothers) before 31 December 2015. No data were provided to the research team on mothers on the registers with a non-psychotic psychiatric illness diagnosed before 31 December 2015, nor on their children [estimated to be approximately 61 000 mothers and 115 000 children (Justice, 2022)], as they were excluded from the original ethics agreement covering the project. Children of the e-Cohort were identified on the Midwives' Notification System (Gee & Dawes, Reference Gee and Dawes1994), which includes mandatory, prospectively collected data on all infants born in Western Australia at 20 weeks gestation or greater or weighing at least 400 g, including home births, and not restricted to live births. Mothers with a psychotic illness were sourced through linkage of records on the Midwives' Notification System to the Hospital Morbidity Data Collection and Mental Health Information System which covers records of all public and private inpatient hospital admissions, as well as public outpatient and ambulatory care contacts with mental health services across the State, dating back to 1966. Psychotic illness was broadly defined using ICD-9 codes 295–298 to include schizophrenia (ICD-9 295; N = 1220 mothers), affective psychoses (ICD-9 296; N = 8365 mothers), paranoid states and other non-organic psychoses (ICD-9 297 and 298; N = 1020 mothers). ICD-8 and ICD-10 equivalents, used in a minority of records, were mapped to ICD-9. For those fathers named on the Midwives' Notification System, relevant covariate data were also extracted from linked State registers. The full cohort of children, their mothers and fathers named on the Midwives' Notification System were linked to electronic mental and physical health data, mortality, child protection and corrective services records. Linkage was carried out by the Data Linkage Branch of the Western Australian Department of Health. Full details on the linked registers used to establish and characterise the e-Cohort have been published previously (Di Prinzio et al., Reference Di Prinzio, Morgan, Bjork, Croft, Lin, Jablensky and McNeil2018; Morgan et al., Reference Morgan, Valuri, Croft, Griffith, Shah, Young and Jablensky2011).

This study was approved by the Western Australian Department of Health Human Research Ethics Committee (2011/75) and The University of Western Australia Human Research Ethics Committee (RA/4/1/1322).

Cohort partitioning

We randomly generated three non-overlapping partitions of children from the e-Cohort: (i) a 40% (N = 171 588) random sample to develop the computer algorithm, (ii) a 30% (N = 128 691) random sample to validate it and explore calibration, and (iii) a final 30% (N = 128 691) partition to repeat validation and calibration on a fresh sample.

Outcome

Children's follow up and schizophrenia diagnosis: Children were followed up from their 10th birthday to 30 June 2015, their date of death or date of onset of schizophrenia, whichever came first. Children were defined as having schizophrenia if they had any recorded diagnosis of schizophrenia (ICD-9 295-all or ICD-10 F20, F21, F23.1, F23.2 or F25) in their mental health inpatient admissions or ambulatory/outpatient contacts; the date of the first record of schizophrenia was used as a proxy for date of onset. Information on child mortality was obtained from the Western Australian Death Registry.

Measures of adversity exposure

Guided by literature, we considered the range of childhood adversities associated with schizophrenia. We calculated measures of such adversities where the scope and accuracy of register data provided sufficient information to do so. Complete descriptions are in the Covariate Dictionary in the online Supplementary materials. These measures were grouped into five broad domains of adversity, A1 to A5:

A1 – Discontinuity in parenting: Before the age of 10 years covered any separation of a child from a parent including hospitalisation, parental death, placement in foster care and parental incarceration. Separations were further categorised by the child's developmental epoch: (i) under 1 year of age, (ii) ages 1–4 and (iii) ages 5–9.

A2 – Family functioning: Parental corrective services contacts included detention, diversionary programmes and non-custodial orders. Child protection contacts were available from 1989 onwards and covered any notification of allegations of child abuse (sexual, physical and emotional) and neglect made to child protection services before a child's 10th birthday concerning the index child or any of their maternal siblings; notifications were further classified as substantiated or unsubstantiated.

A3 – Family structure: Mother's age and paternal age at the time of child's birth were calculated from birth records, as was maternal marital status and child's birth order. Calculation of the size of the family a child belonged to at the time of their 10th birthday was possible using Midwives Notification data up to 2001 and hospital morbidity records for maternal obstetric events after 2001.

A4 – Area-level socioeconomic/demographic environment: Area-level socioeconomic disadvantage and geographical remoteness/level of urbanicity of the mother's neighbourhood were determined from census data using the mother's address at the time of a child's birth. Additional area-level measures from census information reflected the proportion of persons within an area who were Indigenous (of Aboriginal or Torres Strait Islander descent), Australian born, never married, living in one parent families, unemployed, living in a different residence to 5 years previously, living in a different residence to 1 year previously, living in semi-detached dwellings, living in flats, living in rented dwellings, who had no post school qualifications or who spoke a language other than English at home and did not speak English well. An area-level measure of ethnic heterogeneity was calculated (Morgan et al., Reference Morgan, Morgan, Clare, Valuri, Woodman, Ferrante and Jablensky2008b). Area-level crime rates based on the arrest rate per 1000 residents in 2002 were determined. The distribution of socioeconomic disadvantage within a postcode provided a measure of the geographical variability of disadvantage within a localised area.

A5 – Family-level sociodemographic status: A child's father was classified as known if he was named on the birth record. The child's Indigenous status was scored positive if the child and/or either parent was identified as Indigenous (of Aboriginal or Torres Strait Islander descent) in any of the data sources available. Parental Australian state/territory of origin and migrant status were recorded; the latter was categorised according to the affluence of the country of origin using World Bank rankings. Father's occupation was classified according to skill level, as defined by the Australian Bureau of Statistics.

Online Supplementary Table S1 details missing data. Child protection records were incomplete for children born prior to July 1989, affecting 70 474 children across the whole e-Cohort, and resulting in underestimation of derived exposures. Father's occupation status was missing for 3.7% of known fathers (6186/165 754). They were grouped with unemployed fathers for analysis. The level of missing data for remaining exposures was low (<1%) and considered unlikely to introduce bias.

Statistical analysis

We used (i) the largest of our three non-overlapping partitions (40%) to develop the algorithm, (ii) one of the 30% partitions to validate and explore calibration, and (iii) the final 30% partition to repeat validation and calibration, and thus assess generalisation error (Hastie & Fieldman, Reference Hastie and Fieldman2009).

Development: We used a hierarchical approach to model associations between adversity exposures and schizophrenia diagnosis, using the 40% development sample of the e-Cohort. For each domain of adversity, we developed multivariable Cox models considering exposures specific to that domain only. This process was repeated independently for each domain in turn. The resultant models generated five adversity scores for each child, one for each domain (Equation 1). Next, these domain-level scores were used as input to a global multivariable Cox model which estimated their optimal combination to best explain the risk of schizophrenia diagnosis. Hence, a single score – the Early Adversity Scale for Schizophrenia (EAS-Sz) score – could be determined for each child (Equation 2).

To determine domain-specific models, we first considered the parsimonious parameterisation of potential covariates into discrete measures. This was guided by the assessment of estimated bivariate hazard ratios, and graphical assessment of the hazard rate shapes as illustrated by Nelson Aalen plots. Measures that displayed bivariate association with schizophrenia diagnosis of p ⩽ 0.2 were considered as candidates in further model development. Within a domain, all candidate exposures were evaluated in the presence of each other using augmented backwards elimination (Dunkler, Plischke, Leffondre, & Heinze, Reference Dunkler, Plischke, Leffondre and Heinze2014) and a significance level of 5%. Clearly influential subjects based on visual assessment of likelihood displacement [v-martingale residual graphs (Hosmer, Lemeshow, & May, Reference Hosmer, Lemeshow and May2011)] were excluded. From the optimal domain model, each child's domain score was calculated as the sum of the estimated log hazard ratios for the exposures experienced (Equation 1).

(1)$$\eqalign{& score_{child\;j,\, domain\;i} = \mathop \sum \limits_{k_i = 1}^{n_i} {\rm log}( Hazard\;Ratio_{k_i}) } \times [ 1\;if\;child\;j\;exposed\;to\;adversity\;k_i; \;\;0\;if\;not\;exposed] $$

where

i = adversity domain = 1–5,

j = child = 1 to n development sample,

ni = number of terms in the optimal multivariable model for domain i.

The global model considered the five domain scores as continuous variables, each associated with risk of schizophrenia diagnosis. The model also considered the ability of any of the 10 two-way multiplicative interactions to improve the explanation of variation. Augmented backwards elimination and a significance level of 5% were used to determine which minimal set of interactions and raw domain scores were of optimal predictive benefit. Screening and exclusion of influential observations mirrored the domain-level approach. The EAS-Sz score is defined as the global score – a weighted combination of domain scores as per the preferred global model and is detailed in Equation 2.

(2)$$\eqalign{\mathop {EAS-Sz}\limits_{child\;j} & = \mathop \sum \limits_{domain\;i = 1}^5 \matrix{ {{\rm log}( HR_i) \times [ {score_{child\;j, \;domain\;i}} ] } } \cr & \quad + \mathop \sum \limits_{all\;pairs\;ab\;a = 1..4, \;b = 2..5, \;a < b} \;\log ( {HR_{domain\;ab\;interaction}} ) \cr & \quad \times [ {score_{child\;j, \;domain\;a}} ] \times [ {score_{child\;j, \;domain\;b}} ] } $$

where

i = adversity domain = 1–5,

j = child = 1 to n development sample,

score child j, domain i is as defined in Equation 1,

HR is estimated hazard ratio.

Throughout domain-level and global modelling processes, proportionality assumptions were assessed by χ2 tests using Schoenfeld residuals. The robustness of parameter estimates was investigated by comparing them with corresponding estimates determined from equivalent logistic regression models with log(time) offsets.

Validations: We used two separate 30% random samples of our cohort (each N = 128 691), which did not overlap with our development dataset or each other, to conduct validation assessments on domain and global models. To assess discrimination ability of our Cox model, we calculated Harrell's Concordance (Harrell, Califf, Pryor, Lee, & Rosati, Reference Harrell, Califf, Pryor, Lee and Rosati1982). For calibration assessment, we prepared graphical comparisons of predicted v. observed diagnosis prevalence over follow-up, across five risk stratification categories, using the methods of Royston (Royston, Reference Royston2015). The categories were defined based on a pragmatic partitioning of the global adversity scores where their definition balanced the need for sufficient categories to provide useful discrimination, while maintaining stability of parameter estimates derived from such categories. Observed natural cut-points of the score's distribution were exploited where possible.

Results

Our development dataset comprised N = 171 588 births. Children who died before age 10 years (N = 1258) were treated as not at risk of developing schizophrenia and were excluded. Of the remaining 170 330 children, 656 had at least one schizophrenia diagnosis before the end of follow-up. A further 967 children died before the end of follow-up and before receiving a diagnosis. Adversity exposure profiles for children excluded or censored due to death are detailed in online Supplementary Table S2. Loss to follow-up due to outmigration is assumed to be low (Morgan et al., Reference Morgan, Valuri, Croft, Griffith, Shah, Young and Jablensky2011).

After assessment and relevant consolidation of 70 candidate measures, including some re-parametrisation of original versions to comprise fewer category classifications, 42 were found to have potential to be predictive of schizophrenia diagnosis: 15 in the discontinuity in parenting domain, two in family functioning, five in family structure, 16 in area-level socioeconomic/demographic environment and four in family-level sociodemographic status. The covariate dictionary in online Supplementary materials provides more details.

The proportion of children exposed to adversities ranged from 0.4% (having a mother incarcerated before the child's 10th birthday) to 97.3% (being born in an area with all but the lowest of crime rates – total annual crime arrests ⩾50 per 1000 persons) (Table 1). Corresponding proportions for validation datasets can be viewed in online Supplementary Table S3. Associations between adversity exposures and schizophrenia diagnosis are summarised in Table 1 and range from modest (either or both parents deceased before child's 10th birthday, HR = 1.39, 95% CI 0.77–2.53) to large (any placement in foster care before a child's 10th birthday, HR = 9.56, 95% CI 5.80–15.7).

Table 1. Exposure to adversity in development sample (number and percentage) and bivariate association with schizophrenia diagnosis (hazard ratios and 95% confidence intervals)

Italics highlight that numbers so formatted are percentages.

a In area of mother's residence at child's birth, at 2001 census.

Model development for five adversity domains

For each of the five adversity domains, the Cox model chosen to explain variation in schizophrenia diagnosis is summarised in Table 2. Adversity exposures for excluded influential observations can be found in online Supplementary Table S4, and changes in model parameter estimates due to their exclusion can be viewed in online Supplementary Table S5. Tests of proportionality assumptions did not show any substantial violations. Parameter estimates were considered robust to model specification after comparison with those corresponding to logistic regression formulation (online Supplementary Table S6).

Table 2. Risk of developing schizophrenia following exposure to adversity in each of five domains (hazard ratios and 95% confidence intervals)

a In area of mother's residence at child's birth, at 2001 census.

The exposure profile of each child in the development sample was combined with the model parameters of Table 2 using Equation 1 to calculate each child's five continuous domain scores. Summaries of the distributions of these scores are presented in online Supplementary Fig. S1. The great majority of children were allocated scores in the lower ends of the distributions, with only small proportions allocated scores in the upper ranges.

Global model development: combining domain models

Table 3 summarises the chosen global model which includes the five domain scores, together with interaction terms. This was chosen due to its combination of relative high concordance, good calibration and simplicity of application. Parameter estimates were considered robust (online Supplementary Table S6) and proportionality assumptions were considered valid.

Table 3. Risk of developing schizophrenia following exposure to adversity across five domains (hazard ratios and 95% confidence intervals)

a Exposure to adversity in each domain is measured on a continuous scale as a weighted sum of exposures defined in Table 2 and calculated using Equation 1. The hazard ratio of a domain is the estimated factor increase in the hazard for a 1 unit increase in score in that domain.

The global score – the EAS-Sz score – was calculated for each child in the development sample by applying Equation 2 to the model parameters of Table 3. A distributional summary of these scores is contained in online Supplementary Fig. S1.

Model validation: calibration

Five risk stratification categories were defined in the development sample, representing increasing levels of exposure to adversity: (i) negligible: 0 ⩽ EAS-Sz < 0.1 [N = 20 304, with Sz diagnosis = 33 (0.16%)]; (ii) mild: 0.1 ⩽ EAS-Sz < 0.5 [N = 75 815, with Sz diagnosis = 179 (0.24%)]; (iii) moderate: 0.5 ⩽ EAS-Sz < 1.0 [N = 52 639, with Sz diagnosis = 229 (0.44%)]; (iv) elevated: 1.0 ⩽ EAS-Sz < 2 [N = 19 889, with Sz diagnosis = 152 (0.76%)]; and (v) substantial: EAS-Sz ⩾ 2.0 [N = 2941, with Sz diagnosis = 63 (2.14%)]. The corresponding distributions for validation sets 1 and 2 are presented in Table 4 and display dose response. The proportion of children in validation sets 1 and 2 who received a schizophrenia diagnosis before the end of their follow-up increased from 0.19% and 0.26% respectively for those with EAS-Sz defined ‘negligible’ risk to 2.19% and 1.85% respectively for those with ‘substantial’ risk. In Fig. 1, the model assumed average predicted time-to-diagnosis curve for all children in each of the risk categories is presented, together with 95% CI limits of the observed Kaplan–Meier time-to-diagnosis curves. The illustration is repeated for the development sample and each validation sample. The observed Kaplan–Meier curves displayed consistency with model predictions.

Fig. 1. Calibration illustration for five prognostic risk strata of exposure to childhood adversity. Observed 95% limits for Kaplan–Meier curves and predicted time to diagnosis curves.

Table 4. Validation of Cox models describing risk of schizophrenia following exposure to adversity (distribution of schizophrenia diagnoses across risk categories in validation samples for global model and Harrell's Concordance statistics for global and domain models)

Italics highlight that numbers so formatted are percentages.

Model validation: discrimination

Predictive performance for selected models is summarised in Table 4. The global model produced Harrell's Concordance values of 0.656 and 0.624 for validation sets 1 and 2 respectively.

For models describing individual domains, the highest values of Harrell's Concordance were for A3 Family structure (0.609 in validation set 1 and 0.591 in validation set 2).

Discussion

Principal findings

To the best of our knowledge, we have derived and validated the first comprehensive early adversity scale for schizophrenia [EAS-Sz] describing the risk of schizophrenia diagnosis after exposure to broadly defined childhood adversity, using routinely (and prospectively) collected register data. The scale alone accurately predicts schizophrenia above chance, and so has good potential to be combined with other risk factors for schizophrenia as a tool to help untangle the relative contributions of genetic and environmental liability.

Our method estimated the relative impact of various components of adversity exposure, simultaneously adjusting for other exposures. The scale components, all derived from electronic register data, cover a wide range of exposures across multiple broad domains up to a child's 10th birthday and are far greater in number than any sets of exposures previously investigated. Using a large population cohort, we developed the scale using a partitioned 40% sample and evaluated its internal validity using two separate 30% samples. The main scale was a complex global scale which incorporated five domains: (i) discontinuity in parenting, (ii) family functioning, (iii) family structure, (iv) area-level socioeconomic/demographic environment and (v) family-level sociodemographic status, and its values of Harrell's Concordance were well above chance, ranging from 0.624 to 0.656. Although superficially modest, we highlight that these concordance values are achieved for a partial prediction model. We note that the score does not include any proxy measures for genetic liability to schizophrenia (familial diagnosis history) or other risk factors that have been reported in the literature such as cognitive function, obstetric complications and cannabis use, and only includes exposures up to age 10 years. Predicted time-to-diagnosis curves for five risk strata calculated from the development sample were consistent with 95% CI bounds of the corresponding observed Kaplan–Meier curves in both validation samples. In addition to the global scale, we produced five subscales, each based on a single domain of adversity. These also indicated prediction above chance (Harrell's Concordance values range 0.551–0.609). These results indicate that, while some predictive capacity can be achieved from examining adversity exposure in a single domain, more can be gained by considering a broad measure of adversity across many domains. We consider this a particular strength of our scale.

Strengths and limitations

Although simple in application, our scale is unique in the level of sophistication of data used in its derivation, drawing on high-quality State electronic registers of excellent coverage for a large population. The information was prospectively collected, avoiding recall bias. In addition, by accessing data retrospectively, we covered a large time span efficiently. Importantly, the scale allows for different severity weightings for different exposures, some differing by age at occurrence. Finally, by limiting exposure time from birth to a child's 10th birthday, well before the outcome of schizophrenia diagnosis (in all but the rarest of cases), we have avoided potential confounding with premorbid behaviours associated with attenuated or at-risk mental states in early psychosis.

Of great statistical advantage, this scale allows researchers to incorporate many measures of a complex adversity profile into their analyses as a single metric – leaving degrees of freedom free to estimate other contributions to risk of schizophrenia such as genetic liability. See, for example, Vargas and Mittal (Reference Vargas and Mittal2018). As its derivation was achieved through the application of straightforward statistical techniques, other researchers can emulate the process.

The validation of our scale using large datasets which were randomly selected from our cohort gives confidence to the scale's ability to accurately represent risk of schizophrenia in Western Australia, and augers well for its generalisability to other similar populations. However, we believe further validation of EAS-Sz by testing its application with registers from other jurisdictions would confirm whether the relationship between early childhood adversity exposure and schizophrenia outcome that we have quantified holds true when applied to electronic register records elsewhere in the world. This remains an important future step.

Our results have not yet been assessed for specificity to schizophrenia compared to other mental health outcomes, many of which have also been associated with exposure to early adversity.

Our measures of adversity exposure used in constructing EAS-Sz have been derived from records in State registers. We acknowledge that some measures of exposure may not always correspond to an expected increased risk of an undesirable outcome at an individual level and can reasonably be labelled ‘blunt’ measures. However, the substantial difference in prevalence of schizophrenia diagnosis between children classified with ‘substantial’ risk compared to lower risk levels gives support to the ability of EAS-Sz to identify childhood adversity at levels of severity that have greatest impact on the risk of schizophrenia diagnosis.

Despite the breadth of our information, some important exposures could not be considered because they are not recorded in registers, for example, experience of bullying, or family rearing style. Also we did not have access to information for children of mothers on the psychiatric register with psychiatric disorders other than psychotic illness. This was consistent with the ethics approval for the original e-Cohort design; however, we acknowledge that the omission of these children could exclude some patterns in adversity exposure and subsequent schizophrenia risk. In addition, follow-up time for those children born late in the study period was short and, due to their age, they had not passed through the risk period for developing schizophrenia.

Despite minor indications that the model predicted rates of diagnosis at the lower end of the observed 95% Kaplan–Meier estimates for younger ages and towards the higher end of those observed for older ages, we observe that our calibration is generally very good.

Our Department of Child Protection data covering abuse notifications and any subsequent substantiations only include records from 1989 onwards, which reduces the data coverage for around 40% of the children in our cohort. We estimate this translates to an undercount of approximately 20% of children exposed to a substantiated allegation of child abuse in their maternal family (i.e. we estimate we have missed identifying 0.4% of the cohort population as being exposed to a substantiated allegation).

Results in relation to previous studies

There are few reports on the validation and performance of adversity scales. The data-driven method adopted by Pries et al. (Reference Pries, Lage-Castellanos, Delespaul, Kenis, Luykx, Lin and Guloksuz2019) to develop their Exposome Score is one of the closest in approach to ours, although major differences remain: they developed their predictive model using far fewer exposure measures, most of which were determined by recall, and in the case of patients, after diagnosis. In general, their predictive exposures were more chronologically concurrent with outcome than our own, which we deliberately limited to exposures before 10 years of age. Due to the observational structure of their training sample, just under 40% had a diagnosis, time of follow-up was not relevant, and due to the age of the participants (mean = 34 years, s.d. = 10), the likelihood of censored outcome in the controls could be assumed non-existent. The predictive model they developed using logistic regression was found to have an area under the receiver operating curve (AU-ROC) of 0.73 when tested in a separate observational sample (Pries et al., Reference Pries, Lage-Castellanos, Delespaul, Kenis, Luykx, Lin and Guloksuz2019). Pries et al. further validated their Exposome Score, reporting an AU-ROC of 0.84 after application to a third cross-sectional sample (Pries et al., Reference Pries, Erzin, van Os, Ten Have, de Graaf, van Dorsselaer and Guloksuz2020).

In comparison, we measured baseline exposure, prospectively recorded from objective registers, at 10 years old, and incorporated variable follow-up time (mean age 24 years, range 10–35.5 years) with high levels of censoring using Cox regression. For those who developed schizophrenia, the mean age at diagnosis was 22 years. We report Harrell's Concordance (the AU-ROC equivalent for Cox regression) of 0.656 and 0.624 for our two internal validation samples.

Conclusions and future directions

The causal pathway to schizophrenia is complex. Environmental factors have long been implicated in schizophrenia risk, but much remains to be understood about the role of these factors in schizophrenia liability (van Os et al., Reference van Os, Kenis and Rutten2010). Capturing these environmental risks adequately and in a replicable form has been a fraught exercise (Guloksuz, van Os, & Rutten, Reference Guloksuz, van Os and Rutten2018). However, perseverance in developing more comprehensive measures of multidimensional environmental exposures is critical: as Vassos et al. (Reference Vassos, Sham, Kempton, Trotta, Stilo, Gayer-Anderson and Morgan2020) highlight, quantifying cumulative environmental risk is vital, given its potential with respect to intervention and prevention. Our scale offers an objective method of efficiently combining multiple exposures to adversity that, in the literature, have been associated with risk of schizophrenia. The method affords different weightings to different exposures and considers potential interactions, to create a comprehensive global metric of cumulative risk that also includes an option for parsing adversity risk into its underlying sub-domains. To further increase the portability of the scale across platforms and research groups, future developments may include remodelling it as a nominal scale, where each increment in risk is characterised by a smaller set of necessary and sufficient exposures. In developing the nominal scale, approximations with differing levels of detail would be assessed in terms of the relative loss in predictive capacity incurred, against the benefit of improved simplicity of risk measurement calculation.

Our scale provides a nuanced tool for furthering our understanding of the role of cumulative adversity in schizophrenia risk, including its interplay with genetic liability and other exposures such as cognitive impairment and obstetric complications. In current work we are studying these interplays with the global EAS-Sz score and separate domain-level scores, including the specificity of EAS-Sz and its components to schizophrenia, in contrast with other mental health outcomes. Of great importance, its application to whole-population registers allows for retrospective analysis of data from large cohorts to identify modifiable factors along the developmental pathway to schizophrenia that may also influence genetic vulnerability and become potential targets for intervention and prevention.

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/S0033291722001945.

Acknowledgements

Special thanks are extended to Professor Tom McNeil for the original motivation for this work, and for his continued support and advice. We thank the Data Linkage Branch of the Western Australia (WA) Department of Health for data linkage and extraction of the WA data, and for client support. We also thank the custodians of the WA Hospital Morbidity Data Collection, Mental Health Information System, the WA Midwives Notification System, WA Registry of Births, Deaths and Marriages, WA Department of Communities and the WA Department of Justice for the provision of data.

Financial support

This work was supported by grants from the Australian National Health and Medical Research Council (#303235; #458702; APP1002259; APP1080606), March of Dimes (#12-FY04-48; #12-FY07-224), Stanley Foundation and Swedish Research Council (VR; dnr 2019-00198).

Conflict of interest

None.

References

Agid, O., Shapira, B., Zislin, J., Ritsner, M., Hanin, B., Murad, H., … Lerer, B. (1999). Environment and vulnerability to major psychiatric illness: A case control study of early parental loss in major depression, bipolar disorder and schizophrenia. Molecular Psychiatry, 4, 163172.CrossRefGoogle ScholarPubMed
Cannon, M., Jones, P. B., & Murray, R. M. (2002). Obstetric complications and schizophrenia: Historical and meta-analytic review. American Journal of Psychiatry, 159, 10801092.CrossRefGoogle ScholarPubMed
Cantor-Graae, E., & Pedersen, C. B. (2013). Full spectrum of psychiatric disorders related to foreign migration: A Danish population-based cohort study. JAMA Psychiatry, 70, 427435.CrossRefGoogle ScholarPubMed
Di Prinzio, P., Morgan, V. A., Bjork, J., Croft, M., Lin, A., Jablensky, A., & McNeil, T. F. (2018). Intellectual disability and psychotic disorders in children: Association with maternal severe mental illness and exposure to obstetric complications in a whole-population cohort. American Journal of Psychiatry, 175, 12321242.CrossRefGoogle Scholar
Dong, M., Anda, R. F., Felitti, V. J., Dube, S. R., Williamson, D. F., Thompson, T. J., … Giles, W. H. (2004). The interrelatedness of multiple forms of childhood abuse, neglect, and household dysfunction. Child Abuse & Neglect, 28, 771784.CrossRefGoogle ScholarPubMed
Dunkler, D., Plischke, M., Leffondre, K., & Heinze, G. (2014). Augmented backward elimination: A pragmatic and purposeful way to develop statistical models. PLoS ONE, 9, e113677.CrossRefGoogle ScholarPubMed
Evans, G. W., Li, D., & Whipple, S. S. (2013). Cumulative risk and child development. Psychological Bulletin, 139, 13421396.CrossRefGoogle ScholarPubMed
Gee, V., & Dawes, V. (1994). Validation study of the Western Australian Midwives’ Notification System 1992. Occasional Paper 59. Perth, Western Australia: Health Department of Western Australia, p. 18.Google Scholar
Gottesman, I. I. (1994). Complications to the complex inheritance of schizophrenia. Clinical Genetics, 46, 116123.CrossRefGoogle Scholar
Gottesman, I. I., & Erlenmeyer-Kimling, L. (2001). Family and twin strategies as a head start in defining prodromes and endophenotypes for hypothetical early-interventions in schizophrenia. Schizophrenia Research, 51, 93102.CrossRefGoogle ScholarPubMed
Gottesman, I. I., McGuffin, P., & Farmer, A. E. (1987). Clinical genetics as clues to the ‘real’ genetics of schizophrenia (a decade of modest gains while playing for time). Schizophrenia Bulletin, 13, 2347.CrossRefGoogle Scholar
Green, J. G., McLaughlin, K. A., Berglund, P. A., Gruber, M. J., Sampson, N. A., Zaslavsky, A. M., & Kessler, R. C. (2010). Childhood adversities and adult psychiatric disorders in the national comorbidity survey replication I: Associations with first onset of DSM-IV disorders. Archives of General Psychiatry, 67, 113123.CrossRefGoogle ScholarPubMed
Guloksuz, S., van Os, J., & Rutten, B. P. F. (2018). The Exposome paradigm and the complexities of environmental research in psychiatry. JAMA Psychiatry, 75, 985986.CrossRefGoogle ScholarPubMed
Harrell, F. E., Califf, R. M., Pryor, D. B., Lee, K. L., & Rosati, R. A. (1982). Evaluating the yield of medical tests. JAMA-Journal of the American Medical Association, 247, 25432546.CrossRefGoogle ScholarPubMed
Hastie, T. T. R, & Fieldman, J. H. (2009). The elements of statistical learning: Data mining, inference, and prediction. New York: Springer.CrossRefGoogle Scholar
Hosmer, D., Lemeshow, S., & May, S. (2011). Applied survival analysis: Regression modeling of time-to-event data (2nd ed., pp. 169206). Hoboken, NJ: John Wiley & Sons, Inc.Google Scholar
Kessler, R. C., McLaughlin, K. A., Green, J. G., Gruber, M. J., Sampson, N. A., Zaslavsky, A. M., … Williams, D. R. (2010). Childhood adversities and adult psychopathology in the WHO World Mental Health Surveys. The British Journal of Psychiatry, 197, 378385.CrossRefGoogle ScholarPubMed
Marconi, A., Di Forti, M., Lewis, C. M., Murray, R. M., & Vassos, E. (2016). Meta-analysis of the association between the level of cannabis use and risk of psychosis. Schizophrenia Bulletin, 42, 12621269.CrossRefGoogle ScholarPubMed
Morgan, C., Charalambides, M., Hutchinson, G., & Murray, R. M. (2010). Migration, ethnicity, and psychosis: Toward a sociodevelopmental model. Schizophrenia Bulletin, 36, 655664.CrossRefGoogle Scholar
Morgan, C., & Gayer-Anderson, C. (2016). Childhood adversities and psychosis: Evidence, challenges, implications. World Psychiatry, 15, 93102.CrossRefGoogle ScholarPubMed
Morgan, C., Kirkbride, J., Hutchinson, G., Craig, T., Morgan, K., Dazzan, P., … Fearon, P. (2008a). Cumulative social disadvantage, ethnicity and first-episode psychosis: A case-control study. Psychological Medicine, 38, 17011715.CrossRefGoogle ScholarPubMed
Morgan, C., Kirkbride, J., Leff, J., Craig, T., Hutchinson, G., McKenzie, K., … Fearon, P. (2007). Parental separation, loss and psychosis in different ethnic groups: A case-control study. Psychological Medicine, 37, 495503.CrossRefGoogle ScholarPubMed
Morgan, V. A., Di Prinzio, P., Valuri, G., Croft, M., McNeil, T., & Jablensky, A. (2019). Are familial liability for schizophrenia and obstetric complications independently associated with risk of psychotic illness, after adjusting for other environmental stressors in childhood? Australian and New Zealand Journal of Psychiatry, 53, 11051115.CrossRefGoogle ScholarPubMed
Morgan, F., Morgan, V. A., Clare, J., Valuri, G., Woodman, R., Ferrante, A., … Jablensky, A. (2008b). Schizophrenia and offending: Area of residence and the impact of social disorganisation and urbanicity. Trends and Issues in Crime and Criminal Justice, 365, 16.Google Scholar
Morgan, V. A., Valuri, G. M., Croft, M. L., Griffith, J. A., Shah, S., Young, D. J., & Jablensky, A. V. (2011). Cohort profile: Pathways of risk from conception to disease: The Western Australian schizophrenia high-risk e-Cohort. International Journal of Epidemiology, 40, 14771485.CrossRefGoogle ScholarPubMed
Murray, R. M., Bhavsar, V., Tripoli, G., & Howes, O. (2017). 30 Years on: How the neurodevelopmental hypothesis of schizophrenia morphed into the developmental risk factor model of psychosis. Schizophrenia Bulletin, 43, 11901196.CrossRefGoogle ScholarPubMed
Owen, M. J., Sawa, A., & Mortensen, P. B. (2016). Schizophrenia. Lancet, 388, 8697.CrossRefGoogle ScholarPubMed
Padmanabhan, J. L., Shah, J. L., Tandon, N., & Keshavan, M. S. (2017). The ‘polyenviromic risk score’: Aggregating environmental risk factors predicts conversion to psychosis in familial high-risk subjects. Schizophrenia Research, 181, 1722.CrossRefGoogle ScholarPubMed
Paksarian, D., Eaton, W. W., Mortensen, P. B., Merikangas, K. R., & Pedersen, C. B. (2015a). A population-based study of the risk of schizophrenia and bipolar disorder associated with parent–child separation during development. Psychological Medicine, 45, 28252837.CrossRefGoogle ScholarPubMed
Paksarian, D., Eaton, W. W., Mortensen, P. B., & Pedersen, C. B. (2015b). Childhood residential mobility, schizophrenia, and bipolar disorder: A population-based study in Denmark. Schizophrenia Bulletin, 41, 346354.CrossRefGoogle ScholarPubMed
Pries, L. K., Erzin, G., van Os, J., Ten Have, M., de Graaf, R., van Dorsselaer, S., … Guloksuz, S. (2020). Predictive performance of Exposome score for schizophrenia in the general population. Schizophrenia Bulletin, 47(2), 277283.CrossRefGoogle Scholar
Pries, L. K., Lage-Castellanos, A., Delespaul, P., Kenis, G., Luykx, J. J., Lin, B. D., … Guloksuz, S. (2019). Estimating Exposome score for schizophrenia using predictive modeling approach in two independent samples: The results from the EUGEI study. Schizophrenia Bulletin, 45, 960965.CrossRefGoogle ScholarPubMed
Read, J., van Os, J., Morrison, A. P., & Ross, C. A. (2005). Childhood trauma, psychosis and schizophrenia: A literature review with theoretical and clinical implications. Acta Psychiatrica Scandinavica, 112, 330350.CrossRefGoogle ScholarPubMed
Royston, P. (2015). Tools for checking calibration of a Cox model in external validation: Prediction of population-averaged survival curves based on risk groups. Stata Journal, 15, 275291.CrossRefGoogle Scholar
Tapsell, R., Hallett, C., & Mellsop, G. (2018). The rate of mental health service use in New Zealand as analysed by ethnicity. Australasian Psychiatry, 26, 290293.CrossRefGoogle ScholarPubMed
van Os, J., Kenis, G., & Rutten, B. P. (2010). The environment and schizophrenia. Nature, 468, 203212.CrossRefGoogle ScholarPubMed
Varese, F., Smeets, F., Drukker, M., Lieverse, R., Lataster, T., Viechtbauer, W., … Bentall, R. P. (2012). Childhood adversities increase the risk of psychosis: A meta-analysis of patient-control, prospective- and cross-sectional cohort studies. Schizophrenia Bulletin, 38, 661671.CrossRefGoogle ScholarPubMed
Vargas, T., & Mittal, V. (2018). Issues affecting reliable and valid assessment of early life stressors in psychosis. Schizophrenia Research, 192, 465466.CrossRefGoogle ScholarPubMed
Vassos, E., Pedersen, C. B., Murray, R. M., Collier, D. A., & Lewis, C. M. (2012). Meta-analysis of the association of urbanicity with schizophrenia. Schizophrenia Bulletin, 38, 11181123.CrossRefGoogle ScholarPubMed
Vassos, E., Sham, P., Kempton, M., Trotta, A., Stilo, S. A., Gayer-Anderson, C., … Morgan, C. (2020). The Maudsley environmental risk score for psychosis. Psychological Medicine, 50, 22132220.CrossRefGoogle ScholarPubMed
Veling, W., Susser, E., Selten, J. P., & Hoek, H. W. (2015). Social disorganization of neighborhoods and incidence of psychotic disorders: A 7-year first-contact incidence study. Psychological Medicine, 45, 17891798.CrossRefGoogle ScholarPubMed
WA Department of Justice (2022). Statistics – births, deaths and marriages registered: Number of births, deaths and marriages registered in Western Australia 1981–2021. Perth: WA Department of Justice.Google Scholar
Wicks, S., Hjern, A., & Dalman, C. (2010). Social risk or genetic liability for psychosis? A study of children born in Sweden and reared by adoptive parents. American Journal of Psychiatry, 167, 12401246.CrossRefGoogle ScholarPubMed
Wicks, S., Hjern, A., Gunnell, D., Lewis, G., & Dalman, C. (2005). Social adversity in childhood and the risk of developing psychosis: A national cohort study. American Journal of Psychiatry, 162, 16521657.CrossRefGoogle ScholarPubMed
Zammit, S., Lewis, G., Dalman, C., & Allebeck, P. (2010). Examining interactions between risk factors for psychosis. British Journal of Psychiatry, 197, 207211.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Exposure to adversity in development sample (number and percentage) and bivariate association with schizophrenia diagnosis (hazard ratios and 95% confidence intervals)

Figure 1

Table 2. Risk of developing schizophrenia following exposure to adversity in each of five domains (hazard ratios and 95% confidence intervals)

Figure 2

Table 3. Risk of developing schizophrenia following exposure to adversity across five domains (hazard ratios and 95% confidence intervals)

Figure 3

Fig. 1. Calibration illustration for five prognostic risk strata of exposure to childhood adversity. Observed 95% limits for Kaplan–Meier curves and predicted time to diagnosis curves.

Figure 4

Table 4. Validation of Cox models describing risk of schizophrenia following exposure to adversity (distribution of schizophrenia diagnoses across risk categories in validation samples for global model and Harrell's Concordance statistics for global and domain models)

Supplementary material: File

Di Prinzio et al. supplementary material

Di Prinzio et al. supplementary material 1

Download Di Prinzio et al. supplementary material(File)
File 171.8 KB
Supplementary material: PDF

Di Prinzio et al. supplementary material

Di Prinzio et al. supplementary material 2

Download Di Prinzio et al. supplementary material(PDF)
PDF 736.1 KB