- Open Access
Nonshared environmental factors in the aetiology of autism and other neurodevelopmental conditions: a monozygotic co-twin control study
Molecular Autism volume 13, Article number: 8 (2022)
A significant proportion of variation in likelihood of neurodevelopmental conditions (NDCs) has been attributed to nonshared environmental (NSE) factors, although it remains unclear which NSE factors pose specific risks for certain NDCs.
A monozygotic co-twin design was applied in a sample of 224 twins (mean age = 17.70 years, SD = 6.28) controlling for confounders such as genes and shared environment. Generalized estimating equation models were fitted, using perinatal and postnatal indications of NSEs as exposure, operationalized both as separate risk factors and as cumulative risk loads. Categorical and dimensional operationalizations of autism spectrum disorder (ASD), attention-deficit/hyperactivity disorder (ADHD), intellectual disability and other NDCs were used as outcomes.
Birth weight discordance was associated with dimensional autism and ADHD for the smaller twin, and medication during infancy was associated with dimensional autism. Among postnatal factors scarlet fever during early childhood was associated with lower IQ. Especially autism was associated with a greater cumulative perinatal or postnatal risk load.
When exploring the associations between each condition and specific NSEs the risk of being statistically underpowered increases. Hence, we limit the reported findings on specific indicators of NSEs to trait levels and present descriptive data for categorical NDCs.
The findings support previous research by indicating an association between exposure to perinatal and postnatal risks and subsequent NDCs within twin pairs and suggest that autism may be especially linked to accumulative early environmental risks. The findings are potentially important for developmental outcomes prognoses and may inform targeted prevention and early interventions.
Neurodevelopmental conditions (NDCs) are childhood onset disabilities that persist into adulthood and are associated with altered architecture, maturation, and functioning of the developing central nervous system . NDCs affect a wide range of mental functions, such as attention, executive control, top-down versus bottom-up processing, speech-language and communication, motivation, motor control, learning, and social awareness [2, 3]. In the Diagnostic and Statistical Manual of Mental Disorders 5th ed. (DSM-5) and the International Classification of Diseases 11th Revision (ICD-11), the umbrella concept of NDCs includes attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), intellectual disabilities (ID), communication disorders, specific learning disorders, and motor disorders, such as tics and developmental coordination disorder (DCD) [4, 5]. Most NDCs are viewed as extreme manifestations of traits that are continuously distributed in the general population . Even though different NDCs are associated with similar neurocognitive alterations, with an overlap of symptoms and frequent coexisting psychiatric and somatic problems, they differ with regard to their genetic vulnerabilities, clinical course, prognosis and response to treatment . Therefore, the need for both grouping NDCs together while retaining diagnostic distinctions, as well as for including both categorical diagnoses and dimensional traits, has been emphasised .
NDCs are common paediatric challenges, with an estimated prevalence of 17.8% in the years 2015 to 2017 in the USA . Although considered highly heritable, a proportion of NDC causality has been attributed to environmental factors , especially nonshared environmental (NSE) effects, i.e., exposures unique to different children within the same family that make them differ from each other. For example, twin–twin differences in autistic traits are substantially influenced by NSEs at a clinical severity level  and estimates of the NSE impact on ADHD and ASD range between 17 and 44% [10,11,12]. Importantly, the identification of valid NSEs in the aetiology of NDCs might facilitate research on prevention of early adverse influences on the developing child.
NSE, or rather indicators of NSE, include prenatal factors such as gestational hypertension and diabetes, gestational bleeding, maternal use of medication and teratogens [13,14,15,16]; perinatal factors, such as breech presentation, foetal distress, preterm birth, caesarean delivery, birth complications, or birth weight [8, 13, 17]; and postnatal factors such as infections and antibiotic exposure [18, 19]. Despite reported associations between NSE and NDCs, the evidence to support a causal relationship remains insufficient  and given the considerable overlap between different types of NDCs  it remains to be determined if the associations are generic to NDCs or if certain environmental factors are linked to specific conditions. For instance, the interaction between the same NSEs with different genetic vulnerabilities may result in an elevated likelihood for specific NDCs, or the association between NSEs and NDCs may be across diagnoses/conditions, where changes in brain development may mediate the association by means of shared neurocognitive alterations.
Thus, research designs are needed that control for the broad range of familial factors (genetics and shared environmental factors) indicated in NDCs . The discordant monozygotic (MZ) co-twin design which automatically controls for both genetic factors (MZ twins, also called identical, are assumed to share essentially 100% of each other’s genes) and shared environment (including the prenatal environment) offers unique opportunities to investigate possible causal pathways to NDCs by focusing on the within-twin pair differences with respect to unique environmental experiences . There is however a shortage of studies investigating MZ twin populations that are sufficiently phenotyped for concordance and discordance with regard to NDCs. In addition, most previous twin-studies were limited by small sample size, a lack of standardized assessments, and by overly focusing on single NDC diagnoses. In spite of these limitations, some support has been found for a within pairs association between lower birth weight and ASD traits . Similarly, within pairs, the twin with lower birthweight had a higher risk of having ADHD  and traits of ADHD . Birthweight discordance has also been associated with cognitive and verbal developmental delay  and lower IQ performance at 3 years of age  in the lower weight twin. Twin studies have also reported an association between perinatal hypoxia and respiratory distress and ASD  and DCD . In addition, the cumulative load of early medical events has been linked to both ASD and autistic trait severity [28, 29].
In sum, while the general importance of NSE for NDCs is undisputed, it remains unclear which NSE factors may pose a specific risk for a certain NDC and whether exposure to several NSEs may further increase the likelihood of a neuroatypical development. Therefore, this study utilized a MZ co-twin design in order to (1) explore if a wide range of NSE are associated with dimensional (traits) operationalizations of NDCs within twin pairs, (2) test whether different NDCs are associated with different NSE factors and (3) investigate the role of cumulative exposure to several NSE factors for categorical (diagnosis) and dimensional operationalizations of NDCs within twin pairs.
Materials and methods
Sample and study design
In this MZ co-twin study, a total of 224 participants (112 MZ same-sex twin pairs) were included (48.2% females; mean age = 17.70 years, SD = 6.28, range: 8–33 years). In total, 49 participants had a diagnosis of ASD (21 females), 43 ADHD (14 females), 13 ID (5 females), 30 other NDCs (e.g., communication disorders, specific learning disorders or motor disorders; 8 females), and 138 typically developed (i.e., with no NDC diagnosis). Thirty-five participants had two or more NDC diagnoses and were included in more than one of the above categories. In total, 52 pairs were concordant for NDCs and 68 pairs were discordant for NDCs (note that concordant pairs may be discordant at the trait level with within-pair differences in dimensional score points). A subsample was included in a previous publication  and the present study adds to the previous research by focusing specifically on MZ twins and by including and comparing risk factors for different types of NDCs. The sample was recruited from the Roots of Autism and ADHD Twin Study Sweden (RATSS), an ongoing study that includes twins from the population-based Child and Adolescent Twin Study in Sweden  and the Young Adult Twins in Sweden Study, where one or both twins positively screened for ASD or ADHD, but also includes typically developing controls . Included twins are comprehensively clinically phenotyped, and zygosity is determined on a panel of 48 single nucleotide polymorphisms . All background and diagnostic information, including retrospective assessment of NSE risk factors were obtained at the same time, during a diagnostic visit at the clinic and complemented by medical records. The study was approved by the National Swedish Ethical Review Board, and written and oral informed consent was obtained from all participants after the nature of the study procedure had been fully explained.
Clinical DSM-5 consensus diagnoses of NDCs were determined by a group of clinicians during the twins’ 2½ day visit at a clinical research unit using medical history records, diagnostic interviews and by first choice standardized diagnostic tools (for more detail and references please see [31, 33]). These include the Schedule for Affective Disorders and Schizophrenia for School-Age Children-Present and Lifetime Version (K-SADS-PL) and the Structured Clinical Interview for DSM-IV-Axis I Disorders (SCID-I), depending on the participant’s age; the Diagnostic Interview for ADHD in adults (DIVA); the Autism Diagnostic Interview—Revised (ADI-R) and the Autism Diagnostic Observation Schedule Second Edition (ADOS-2, modules 3 and 4), neuropsychological tests; and the Adaptive Behavior Assessment System-2 (ABAS-2). When DSM-IV specific instruments were used, a translation to DSM-5 criteria was made by the team of clinicians.
The total raw score points of the parent-report version of the Social Responsiveness Scale Second Edition (SRS-2) were applied to measure autistic traits [34, 35]. The questionnaire was completed during the visit at the clinical research unit. The SRS includes 65 items on interpersonal behavior, communication, inflexible, restricted, and repetitive interests and behaviors. The items are rated on a 4-point Likert-scale ranging from 0 (“not true”) to 3 (“almost always true”), resulting in a total score between 0 and 195, a higher score indicating more autistic traits. Differences of ≥ 10 raw points have in previous research been deemed as meaningful . For this study, SRS was used as a continuous measure.
ADHD traits were measured during the visit at the clinical research unit with the DSM-oriented ADHD subscale of the parent-rated Child Behavior Checklist (CBCL) and the Adult Behavior Checklist (ABCL) depending on the participants age . The items were rated from 0 (“not true”) to 2 (“very often true”), with higher scores indicating more inattention, hyperactivity and impulsivity. The total score could range between 0 and 14 in the CBCL and 0 and 26 in the ABCL. For this study however, CBCL and ABCL score were separately standardized into a Z-score and merged.
The General Ability Index (GAI) of the Wechsler Intelligence Scales for Children or Adults-IV (WISC-IV/WAIS-IV) was used to assess IQ. The GAI is a composite score that is based on three Verbal Comprehension and three Perceptual Reasoning subtests, and the score can range between 40 and 160 where GAI scores between 90 and 109 are considered average.
Nonshared environmental (NSE) factors
Medical and developmental history with focus on the first 5 years of life was collected from the RATSS parent-report questionnaire completed during the visit at the clinical research unit, designed to cover pre-, peri- and postnatal factors and child medical history. The questionnaire has been validated against medical registry data . Based on 13 MZ pairs discordant for NDCs and 13 MZ TD pairs, information on birthweight showed an excellent agreement with the medical records (ICC = 0.93, CI = 0.88–0.96), whereas the agreement was moderate for cumulative load of early medical events (ICC = 0.55, 95% CI = 0.22–0.75) . Information on agreement for specific indicators of NSEs is shown in more detail in Additional file 1: Table S1. Most questions had a yes/no response format with a request to specify for positive answers. When information on any particular risk factor was left blank (less than 5% of all reports) it was considered as non-present. As specific and separate exposure information for each twin is needed for an informative within-pair analysis, which is often impossible to derive in twins that share prenatal milieu, the analyses of prenatal environmental factors were not included in this study. For the analyses of perinatal and postnatal factors, we selected those which have been highlighted in previous research on psychiatric and NDCs [8, 13, 17,18,19].
Birthweight information for each twin was collected with a low birthweight defined as ≤ 2000 g (scored as 1), a definition applied in previous twin studies, e.g., . In addition, considering that twin weight discordance is an indicator of divergent twin growth and that the difference among the twins’ birthweight at 15–25% is associated with an increased risk of morbidity and mortality rates as well as other perinatal complications , the twin birth weight discordance was used as an additional risk factor. The cut-off for the birth weight discordance ([larger twin weight − smaller twin weight]/larger twin weight × 100) indicating the potential risk was set at 18% or greater, using a stringent approach suggested by Breathnach et al. , where the twin with a lower birthweight was coded as 1. Parent-reported convulsions, oxygen therapy (as an indicator for hypoxia and respiratory distress), phototherapy (as an indicator for hyperbilirubinemia), breech presentation, and use of medication (e.g., antibiotics, iron) during infancy, were also coded as 1 if present. A cumulative risk load score was calculated as a sum of the single risk factors, with a possible range from 0 to 7, a higher score indicating more adverse perinatal factors.
Parents’ reports on whether their child had experienced measles, jaundice, mumps, scarlet fever, asthma, frequent ear infections, head injury, convulsions or heart disease during early childhood were collected, and a cumulative load was generated based on the single events ranging between 0 and 9, with a higher score indicating more postnatal risk.
Within-pair distributions of autistic traits, traits of ADHD and IQ in terms of the level of twin-twin discordance are shown in Fig. 1. We examined within-pair association between NSE factors and NDCs using conditional regression analyses within the generalized estimating equations (GEE) framework implemented in the “drgee” package of R which does not make any distributional assumptions [40, 41]. In this co-twin design, the difference in the exposure variable within a pair is correlated to the difference in the outcome variable within the same pair (see Additional file 2: Fig. S1). For the continuous outcome (i.e., traits and number of NDCs), we used the identity link within the drgee function, and for binary outcomes (i.e., diagnosis of NDC) we used the log link within the drgee function and the log risk ratio (RR) are presented in back transformed format. The framework accounts for the use of related individuals in the analysis, including calculation of robust standard errors [40, 42]. In the regressions, all factors shared within clusters (here: twin pairs) are adjusted for. In the case of twins, these include additive and non-additive genetic variants, and environmental exposures that increase familial resemblance (i.e., shared environment). No other factors were adjusted for in the models. Two models were fitted for peri- and postnatal NSEs separately. First, we conducted a set of analyses with the single peri- and postnatal NSEs as exposures and dimensional NDCs as separate outcomes. Since the model reduces the number of informative participants for using categorical outcomes (being discordant for specific NDC and also being discordant for exposure), we only present descriptive data for diagnostic outcome. Second, the cumulative environmental risk loads for perinatal and postnatal factors were used as exposure and number of NDC diagnoses, specific NDCs, i.e., ASD, ADHD, ID, as well as autistic traits, traits of ADHD and IQ, were used as separate outcomes. In addition, all other NDCs were combined into one category. For the cumulative environmental risk loads, all included NSEs were weighted the same (i.e., present or not-present). P-values < 0.05 were considered significant. The Benjamini–Hochberg false discovery rate correction (FDR corr) of 0.05 was additionally applied in order to correct for multiple comparisons . We also present confidence intervals and Z statistics for the reported associations.
Figure 1 depicts the distribution of and differences in the key dimensional constructs in terms of level of twin-twin discordance. In nine pairs twin-to-twin transfusion was reported. One participant had reported hydrocephalus, one cerebral haemorrhage, one craniosynostosis, one encephalitis and one had meningitis during childhood. Among the participants, 27.7% of the twins had one, 12.9% two and 30.9% three or more perinatal risk factors reported. The corresponding numbers for postnatal risk factors were 32.1% for one, 16.1% for two and 7.1% for three or more cumulative factors. Partial correlation, adjusting for sex and age, between the NSEs and NDCs are presented in Additional file 3: Table S2.
Non-shared perinatal factors
Table 1 shows the percentage of twins with single perinatal risk factor exposure and mean value of the cumulative risk load stratified by categorical NDC concordance, while Table 2 shows the models for neurodevelopmental trait outcomes. Analyses showed that in pairs with birth weight discordance ≥ 18%, the lower weighted twin has a higher risk of autistic trait score and higher traits of ADHD. More specifically the lower weighted twin scored almost half a standard deviation higher on the SRS-2 (13.9 raw points) and the CBCL/ABCL ADHD subscale compared to their co-twin. The association between birth weight discordance and traits of ADHD, but not ASD traits, remained when applying FDR corr. As shown in Table 1, all of those reporting ≥ 18% lower birth weight in discordant pairs had either a diagnosis of ASD or ADHD. The association between birth weight discordance and autistic traits and traits of ADHD, as well as lower IQ, for the smaller twin was also present when applying birthweight as a continuous measure, see Additional file 4: Table S3. The twin who received medication during infancy had higher levels of autistic traits compared to their co-twin, a finding that remained after FDR corr.
Tables 2 and 3 show the models for dimensional and categorical neurodevelopmental outcomes with cumulative loads of NSE risk factors as exposure. Cumulative load of perinatal factors was associated with an increase in number of NDC diagnoses. When investigating specific NDCs, the perinatal risk load contributed to ASD and autistic traits but not to any other traits or NDC diagnoses. After applying the FDR corr, the association between the perinatal risk load and increase number of NDC diagnoses and autistic traits remained.
Nonshared postnatal factors
As shown in Table 2, postnatal convulsions, frequent ear infections and head injury were associated with autistic traits, but the association did not survive FDR corr. The twin with reported convulsions scored on average 1 standard deviation higher on the SRS-2 (32 raw points) than their co-twin. In addition, scarlet fever and jaundice were associated with lower IQ, and the association between scarlet fever and IQ remained after FDR corr. As shown in Table 1, all of those reporting convulsions in discordant pairs had either diagnosis of ASD, ADHD or ID.
The postnatal cumulative risk load was associated with an increasing number of NDC diagnoses. The cumulative load also contributed to both ASD and autistic traits, as well as ID and IQ, see Tables 2 and 3. After applying the FDR corr, the association between the postnatal cumulative risk load and autistic traits remained.
We investigated the significance of nonshared peri- and postnatal environmental factors in neurodevelopmental conditions in a monozygotic co-twin design. Cumulative perinatal and postnatal risks were associated with neurodevelopmental conditions, and especially autistic traits, within the pairs. For perinatal NSE, divergent twin weight was associated with traits of ASD and ADHD for the smaller twin, and among discordant pairs, only the twin with ASD or ADHD diagnosis had a birthweight ≥ 18% compared to their co-twin. In addition, the twin who received medication during infancy had more traits of ASD compared to their co-twin. For postnatal NSE, the twin with convulsions had more ASD traits, and at a diagnostic level only those with a diagnosis of ASD, ADHD or ID reported convulsions within discordant pairs. In addition, frequent ear infections and head injury were associated with ASD traits, whereas scarlet fever and jaundice were associated with lower IQ.
Since unmeasured familial confounders, such as genes and shared environment are adjusted for, the discordant MZ co-twin design offers unique opportunities to investigate possible causal NSE pathways to NDCs. In our study, twin birth weight discordance was the perinatal factor contributing most to the likelihood of emerging autistic and ADHD traits, as well as lower IQ when applying birthweight as a continuous measure, for the twin with less weight. In line with this, both ASD traits , ADHD , traits of ADHD  and IQ  have been previously associated with birth-weight difference, with the smaller twin showing increased severity of problems. Moreover, birth weight discordance may double the risk of other adverse perinatal outcomes among twins, including mortality, hypoxic-ischemic encephalopathy, respiratory distress, or sepsis .
Since twins have the same gestational age, a difference in their birth weights occurs due to the factors affecting the individual growth of each twin. Possible mechanisms for discordant growth include twin–twin transfusion syndrome or unequal placental sharing , although the association between weight discordance and adverse outcomes in the study by Breathnach and colleagues  remained even after excluding twins with transfusion syndrome . Other factors that may influence birth weight discordance include maternal factors, such as older age and tobacco use, fetal risk factors, such as viral infections, and placental factors, such as velamentous cord insertion . Hypothetically, a growth discordance might influence brain development in ways specific to ASD and ADHD, as it has been reported that variation in birth weight among MZ twins may be associated with alterations of brain anatomy persisting into late adolescence with lower birth weight being associated with reduced brain volume, and more specifically with altered cortical surface area , a cerebral area of significance to ASD . The use of medication reported in our study, such as antibiotics, during infancy for the twin with higher autistic traits may increase microbiome depletion and bacterial dysbiosis, also reported in children with ASD .
For postnatal NSE factors, convulsions were associated with ASD traits, albeit no longer significant when correcting for multiple comparisons using Benjamini–Hochberg FDR corr. Even though the cause of convulsions was not investigated in this study, and only few participants had reported convulsions, previous research suggested that children with epileptic and non-epileptic convulsions in early childhood have a higher degree of neurodevelopmental and psychopathology scores, as compared to controls, at least nine years later [46, 47]. Epilepsy has also been more commonly reported among individuals with ASD, as compared to peers with ADHD . A history of convulsions has also been linked to structural correlates, such as decreased hippocampal volume and white matter tract changes . In our study, head injury during early childhood and frequent ear infections, was associated with more autistic traits, and scarlet fever and jaundice were associated with lower IQ. The association between scarlet fever and IQ remained significant even after using the FDR corr. Even though streptococcal infection has been associated with a wider range of NDCs  and exposure to infections in early childhood has been associated with lower IQ in adolescence , more research is needed to corroborate our findings. The association with head injury needs to be interpreted with caution because of the moderate degree of agreement between parental reports and medical records and the association lost its significance after the FDR corr.
Generally, our findings are consistent with the previous research on NSE in NDCs , indicating, an association between autism and perinatal and postnatal risks within twin pairs. Our finding thus suggests that autism may be particularly sensitive to cumulative early environmental risk load. Previous research has also reported a potential cumulative effect of the number of environmental adversities on the severity of ASD and autistic traits [28, 29, 51]. Accordingly, it has been suggested that ASD may be triggered by additive or cascading effects of high-level exposure to endogenous and environmental factors, hypothesised to affect the immune system and central nervous system functions during critical developmental stages . In contrast with findings from previous twin and sibling studies , perinatal hypoxia and respiratory distress were not associated with neurodevelopmental outcomes in our study. This negative finding could reflect limited power in our model and that we used a proxy (oxygen therapy) for respiratory distress. We did however observe 11-raw score points more on the SRS-2 for the twin who received oxygen, although not significant.
In order to fully evaluate the validity of our findings, some limitations of the current study need to be discussed. First, data on NSE were based on retrospective parental inquiry, and may be subject to confirmation or recall bias. To maximize data validity, the self-reported medical information from home-stored medical reports and record-keepings was validated against medical registry data. While ratings on birthweight had an excellent agreement with the medical records, other NSEs had a lower agreement, e.g., good agreement for oxygen therapy (hypoxia and respiratory distress) and moderate agreement for breech position, which should be considered when interpreting the results. Second, the NSE factors should be understood as rather indirect indicators of environmental risks for NDCs, e.g., phototherapy, why future studies should seek to assess primary medical events as authentically as possible, such as hyperbilirubinemia. Third, several other potential indicators of NSEs that may be of importance, including Apgar score, were not included. Since twins share prenatal milieu, prenatal environmental factors such as maternal use of medication and teratogens could not be included in the within-pair analysis. Fourth, considering the design of the study, we were unable to draw any firm conclusions about the possible causality. Even though the MZ-twin design provides a unique opportunity to untangle the contribution of NSEs on development, reverse causality could not be excluded where NDCs or other factors related to these conditions may increase the risk for NSEs in early development. Fifth, although an MZ co-twin design implies a high degree of genetic and shared environment control, there may be post-twining de novo mutations that may contribute to a minority of cases of discordance. Sixth, concerns have been raised that findings in twin samples may not be generalizable to non-twin populations. However, previous research suggests that twins do not systematically differ from the general population of non-twins on measures of behavior and development [53, 54], including the prevalence of autism . Our sample is also biased toward twins discordant for NDCs. Although this recruitment strategy makes it inappropriate to use bivariate twin models such as A (additive genetics) C (common/shared environment) and E (unique environment) ACE models, it is appropriate for investigating the association between environmental risk factors and outcomes while controlling for familial confounders. Seventh, when exploring the associations between each condition and specific indicators for NSEs the risk of being statistically underpowered increases. In addition, as the case with birth weight discordance and postnatal convulsions, no NSE was reported among the non-affected twins, why it is not feasible to obtain confidence intervals for the RRs in a small sample size using standard measures. Therefore, we only include descriptive data on the association between specific NSEs and NDCs. Lastly, the validity of applying a specific p-value could be discussed. In this study we conducted numerous statistical calculations which increases the risk for type I errors. At the same time, correction for multiple testing with e.g., Bonferroni correction, has been criticized since the interpretations depend on the number of tests that are performed and the likelihood of type II errors is also increased . In this study, we therefore also present results after correcting for multiple comparisons using Benjamini–Hochberg procedure with FDR corr , which are deemed to be less conservative than the Bonferroni correction using familywise error rate.
On balance, the strengths of the study include a relatively large and rare sample of well phenotyped MZ twins, including a considerable number of MZ twins in Sweden concordant and discordant for NDCs (for a specification of number of concordant and discordant pairs from a nation-wide population-based study in Sweden see , for a more comprehensive review on previous twin research in ASD see ). Furthermore, the applied discordant MZ-twin design is a powerful method for investigating the contribution of NSE to neurodevelopmental conditions.
The current study offers insights into early risk factors associated with neurodevelopmental problems and provides additional perspectives on prognosis and prevention. Weight discordance at birth was an important risk factor for both autistic and ADHD traits. Medication early in life, as well the presence of convulsions and streptococcal infection in the first year of life were also associated with traits of NDCs in the within-pair analyses. Even though our finding needs to be interpreted with caution due to a small number of pairs being discordant for specific NSEs, the findings extend previous research by applying an MZ co-twin control design with a high degree of genetic and shared environment control. Our results support a relationship between early NSE and NDCs, where autism seems to be particularly sensitive to a cumulative impact of negative early environmental events. As previously suggested  and as supported by the present study, autistic traits tend to be continuously distributed even in discordant twins. It may be hypothesized that when the traits approximate a diagnostic threshold, multiple exposures to specific risk factors, such as those described in the present study may act as triggers and aggravate the phenotype to a clinical level. The findings highlight the need for interventions targeting factors contributing to fetal growth reduction including viral infections and maternal health problems. Improving prenatal care might influence a number of risk NSEs. Further identification of valid NSEs in the aetiology of NDCs may facilitate prevention efforts with regard to early adverse influences on the developing child.
Availability of data and materials
The datasets presented in this article are not readily available because of the regulations in the ethical approval and university policies, requiring among others a data sharing agreement. Requests to access the datasets should be directed to SB, firstname.lastname@example.org.
Adaptive behavior assessment system
Adult Behavior Checklist
Autism Diagnostic Inventory-Revised
Autism Diagnostic Observation Schedule—Second Edition
Autism Spectrum Disorder
Child Behavior Checklist
Developmental coordination disorder
Diagnostic Interview for ADHD in Adults
Diagnostic and statistical manual of mental disorders
General Ability Index
Generalized estimating equations
Intraclass correlation coefficient
International classification of diseases
Schedule for affective disorders and schizophrenia for school-age children-present and lifetime version
Roots of autism and ADHD Twin Study Sweden
Structured clinical interview for DSM-IV-Axis I disorders
Social Responsiveness Scale
Wechsler Intelligence Scales for Children or Adults-IV
Thapar A, Cooper M, Rutter M. Neurodevelopmental disorders. Lancet Psychiatry. 2017;4(4):339–46.
Bölte S, Mahdi S, Coghill D, Gau SS, Granlund M, Holtmann M, Karande S, Levy F, Rohde LA, Segerer W, de Vries PJ, Selb M. Standardised assessment of functioning in ADHD: consensus on the ICF core sets for ADHD. Eur Child Adolesc Psychiatry. 2018;27(10):1261–81.
Bölte S, Mahdi S, de Vries PJ, Granlund M, Robison JE, Shulman C, Swedo S, Tonge B, Wong V, Zwaigenbaum L, Segerer W, Selb M. The Gestalt of functioning in autism spectrum disorder: results of the international conference to develop final consensus International Classification of Functioning. Disabil Health Core Sets Autism. 2019;23(2):449–67.
Battle DE. Diagnostic and statistical manual of mental disorders (DSM). Codas. 2013;25(2):191–2.
World Health Organization. International statistical classification of diseases and related health problems, 11th edn. 2019. https://icd.who.int/.
Martin J, Taylor MJ, Lichtenstein P. Assessing the evidence for shared genetic risks across psychiatric disorders and traits. Psychol Med. 2018;48(11):1759–74.
Zablotsky B, Black LI, Maenner MJ, Schieve LA, Danielson ML, Bitsko RH, Blumberg SJ, Kogan MD, Boyle CA. Prevalence and trends of developmental disabilities among children in the United States: 2009–2017. Pediatrics. 2019;144(4):e20190811.
Carlsson T, Molander F, Taylor MJ, Jonsson U, Bölte S. Early environmental risk factors for neurodevelopmental disorders—a systematic review of twin and sibling studies. Dev Psychopathol. 2020;24:1–48.
Castelbaum L, Sylvester CM, Zhang Y, Yu Q, Constantino JN. On the nature of monozygotic twin concordance and discordance for autistic trait severity: a quantitative analysis. Behavi Genet. 2020;50:263–72.
Lichtenstein P, Carlstrom E, Rastam M, Gillberg C, Anckarsater H. The genetics of autism spectrum disorders and related neuropsychiatric disorders in childhood. Am J Psychiatry. 2010;167(11):1357–63.
Sandin S, Lichtenstein P, Kuja-Halkola R, Hultman C, Larsson H, Reichenberg A. The heritability of autism spectrum disorder. JAMA. 2017;318(12):1182–4.
Bai D, Yip BHK, Windham GC, Sourander A, Francis R, Yoffe R, Glasson E, Mahjani B, Suominen A, Leonard H, Gissler M, Buxbaum JD, Wong K, Schendel D, Kodesh A, Breshnahan M, Levine SZ, Parner ET, Hansen SN, Hultman C, Reichenberg A, Sandin S. Association of genetic and environmental factors with autism in a 5-country cohort. JAMA Psychiat. 2019;76(10):1035–43.
Gardener H, Spiegelman D, Buka SL. Prenatal risk factors for autism: comprehensive meta-analysis. Br J Psychiatry. 2009;195(1):7–14.
Wang C, Geng H, Liu W, Zhang G. Prenatal, perinatal, and postnatal factors associated with autism: a meta-analysis. Medicine. 2017;96(18):e6696.
Mill J, Petronis A. Pre- and peri-natal environmental risks for attention-deficit hyperactivity disorder (ADHD): the potential role of epigenetic processes in mediating susceptibility. J Child Psychol Psychiatry. 2008;49(10):1020–30.
Thapar A, Cooper M, Eyre O, Langley K. What have we learnt about the causes of ADHD? J Child Psychol Psychiatry. 2013;54(1):3–16.
Sciberras E, Mulraney M, Silva D, Coghill D. Prenatal risk factors and the etiology of ADHD—review of existing evidence. Curr Psychiatry Rep. 2017;19(1):1.
Bittker SS, Bell KR. Acetaminophen, antibiotics, ear infection, breastfeeding, vitamin D drops, and autism: an epidemiological study. Neuropsychiatr Dis Treat. 2018;14:1399–414.
Lavebratt C, Yang LL, Giacobini M, Forsell Y, Schalling M, Partonen T, Gissler M. Early exposure to antibiotic drugs and risk for psychiatric disorders: a population-based study. Transl Psychiatry. 2019;9(1):317.
Vitaro F, Brendgen M, Arseneault L. The discordant MZ-twin method: One step closer to the holy grail of causality. Int J Behav Dev. 2009;33(4):376–82.
Losh M, Esserman D, Anckarsater H, Sullivan PF, Lichtenstein P. Lower birth weight indicates higher risk of autistic traits in discordant twin pairs. Psychol Med. 2012;42(5):1091–102.
Lehn H, Derks EM, Hudziak JJ, Heutink P, van Beijsterveldt TCEM, Boomsma DI. Attention problems and attention-deficit/hyperactivity disorder in discordant and concordant monozygotic twins: evidence of environmental mediators. J Am Acad Child Adolesc Psychiatry. 2007;46(1):83–91.
Pettersson E, Sjölander A, Almqvist C, Anckarsäter H, D’Onofrio BM, Lichtenstein P, Larsson H. Birth weight as an independent predictor of ADHD symptoms: a within-twin pair analysis. J Child Psychol Psychiatry. 2015;56(4):453–9.
Bellido-González M, Defior-Citoler S, Angel D-López M. Cognitive and verbal development of discordant twins without neurological morbidity. J Reprod Infant Psychol. 2007;25(2):161–8.
Ross GS, Krauss AN, Perlman JM. Physical growth and cognitive abilities in concordant versus discordant birth weight twins at three years old. Early Hum Dev. 2012;88(9):753–6.
Froehlich-Santino W, Londono Tobon A, Cleveland S, Torres A, Phillips J, Cohen B, Torigoe T, Miller J, Fedele A, Collins J, Smith K, Lotspeich L, Croen LA, Ozonoff S, Lajonchere C, Grether JK, O’Hara R, Hallmayer J. Prenatal and perinatal risk factors in a twin study of autism spectrum disorders. J Psychiatr Res. 2014;54:100–8.
Pearsall-Jones JG, Piek JP, Rigoli D, Martin NC, Levy F. An investigation into etiological pathways of DCD and ADHD using a monozygotic twin design. Twin Res Hum Genet. 2009;12(4):381–91.
Willfors C, Carlsson T, Anderlid BM, Nordgren A, Kostrzewa E, Berggren S, Ronald A, Kuja-Halkola R, Tammimies K, Bölte S. Medical history of discordant twins and environmental etiologies of autism. Transl Psychiatry. 2017;7(1):e1014.
Bölte S, Girdler S, Marschik PB. The contribution of environmental exposure to the etiology of autism spectrum disorder. Cell Mol Life Sci. 2019;76(7):1275–97.
Anckarsäter H, Lundström S, Kollberg L, Kerekes N, Palm C, Carlström E, Långström N, Magnusson PK, Halldner L, Bölte S, Gillberg C, Gumpert C, Råstam M, Lichtenstein P. The child and adolescent twin study in Sweden (CATSS). Twin Res Hum Genet. 2011;14(6):495–508.
Bölte S, Willfors C, Berggren S, Norberg J, Poltrago L, Mevel K, Coco C, Fransson P, Borg J, Sitnikov R, Toro R, Tammimies K, Anderlid BM, Nordgren A, Falk A, Meyer U, Kere J, Landén M, Dalman C, Ronald A, Anckarsäter H, Lichtenstein P. The roots of autism and ADHD twin study in Sweden (RATSS). Twin Res Hum Genet. 2014;17(3):164–76.
Hannelius U, Gherman L, Mäkelä VV, Lindstedt A, Zucchelli M, Lagerberg C, Tybring G, Kere J, Lindgren CM. Large-scale zygosity testing using single nucleotide polymorphisms. Twin Res Hum Genet. 2007;10(4):604–25.
Isaksson J, Taylor MJ, Lundin K, Neufeld J, Bölte S. Familial confounding on the ability to read minds: a co-twin control study. Autism. 2019;23(8):1948–56.
Constantino JN, Gruber CP. Social responsiveness scale-2. 2nd ed. Western Psychological Services; 2012.
Zander E, Bölte S. Social Responsiveness Scale (SRS-2), Swedish version. 2nd ed. Hogrefe; 2019.
Achenbach TM, Rescorla LA. Manual for the ASEBA school-age forms & profiles. University of Vermont, Research Centre for Children, Youth, & Families, Burlington; 2001.
Johansson S, Iliadou A, Bergvall N, dé Fairé U, Kramer MS, Pawitan Y, Pedersen NL, Norman M, Lichtenstein P, Cnattingius S. The association between low birth weight and type 2 diabetes: contribution of genetic factors. Epidemiology. 2008;19(5):659–65.
Miller J, Chauhan SP, Abuhamad AZ. Discordant twins: diagnosis, evaluation and management. Am J Obstet Gynecol. 2012;206(1):10–20.
Breathnach FM, McAuliffe FM, Geary M, Daly S, Higgins JR, Dornan J, Morrison JJ, Burke G, Higgins S, Dicker P, Manning F, Mahony R, Malone FD. Perinatal Ireland Research Consortium. Definition of intertwin birth weight discordance. Obstet Gynecol. 2011;118(1):94–103.
Zetterqvist J, Sjölander A. Doubly Robust Estimation with the R Package drgee. Epidemiol Methods. 2015;4(1):69–86.
Goetgeluk S, Vansteelandt S. Conditional generalized estimating equations for the analysis of clustered and longitudinal data. Biometrics. 2008;64(3):772–80.
Neuhaus JM, McCulloch E. Separating between and within-cluster covariate effects by using conditional and partitioning methods. Stat Methodol. 2006;68:859–72.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Roy Stat Soc: Ser B (Methodol). 1995;57(1):289–300.
Raznahan A, Greenstein D, Lee NR, Clasen LS, Giedd JN. Prenatal growth in humans and postnatal brain maturation into late adolescence. PNAS. 2012;109(28):11366–71.
Mensen VT, Wierenga LM, van Dijk S, Rijks Y, Oranje B, Mandl RC, Durston S. Development of cortical thickness and surface area in autism spectrum disorder. Neuroimage Clin. 2016;7(13):215–22.
Chin RFM. The outcomes of childhood convulsive status epilepticus. Epilepsy Behav. 2019;101(Pt B):106286.
Gillberg C, Lundstrom S, Fernell E, Nilsson G, Neville B. Febrile seizures and epilepsy: association with autism and other neurodevelopmental disorders in the child and adolescent twin study in Sweden. Pediatr Neurol. 2017;74:80–6.
Pan PY, Bölte S, Kaur P, Jamil S, Jonsson U. Neurological disorders in autism: a systematic review and meta-analysis. Autism. 2021;25(3):812–30.
Ajdacic-Gross V, Aleksandrowicz A, Rodgers S, Mutsch M, Tesic A, Müller M, Kawohl W, Rössler W, Seifritz E, Castelao E, Strippoli MF, Vandeleur C, von Känel R, Paolicelli R, Landolt MA, Witthauer C, Lieb R, Preisig M. Infectious, atopic and inflammatory diseases, childhood adversities and familial aggregation are independently associated with the risk for mental disorders: results from a large Swiss epidemiological study. World J Psychiatry. 2016;6(4):419–30.
Khandaker GM, Dalman C, Kappelmann N, Stochl J, Dal H, Kosidou K, Jones PB, Karlsson H. Association of childhood infection with IQ and adult nonaffective psychosis in Swedish men: a population-based longitudinal cohort and co-relative study. JAMA Psychiat. 2018;75(4):356–62.
Chien YL, Chou MC, Chou WJ, Wu YY, Tsai WC, Chiu YN, Gau SS. Prenatal and perinatal risk factors and the clinical implications on autism spectrum disorder. Autism. 2019;23(3):783–91.
Crawford S. On the origins of autism: the quantitative threshold exposure hypothesis. Med Hypotheses. 2015;85(6):798–806.
Barnes JC, Boutwell BB. A demonstration of the generalizability of twin-based research on antisocial behavior. Behav Genet. 2013;43(2):120–31.
Kendler KS, Pedersen NL, Farahmand BY, Persson PG. The treated incidence of psychotic and affective illness in twins compared with population expectation: a study in the Swedish Twin and Psychiatric Registries. Psychol Med. 1996;26(6):1135–44.
Sandin S, Lichtenstein P, Kuja-Halkola R, Larsson H, Hultman CM, Reichenberg A. The familial risk of autism. JAMA. 2014;311(17):1770–7.
Perneger TV. What’s wrong with Bonferroni adjustments. BMJ. 1998;316:1236–8.
Lichtenstein P, Carlström E, Råstam M, Gillberg C, Anckarsäter H. The genetics of autism spectrum disorders and related neuropsychiatric disorders in childhood. Am J Psychiatry. 2010;167(11):1357–63.
Willfors C, Tammimies K, Bölte S. Twin research in autism spectrum disorder. In: Casanova MF, El-Baz A, Suri JS, editors. Autism imaging and devices. CRC Press; 2017. p. 15–36.
We express our sincere gratitude to all twins and parents who have participated in this research. We also thank the research team at the Center of Neurodevelopmental Disorders at Karolinska Institutet (KIND), RATSS collaborators. We sincerely acknowledge The Swedish Twin Registry for access to data.
Open access funding provided by Karolinska Institute. Funding for this study was provided by Swedish Research Council, Vinnova, Formas, FORTE, Hjärnfonden, Stockholm Brain Institute, Autism and Asperger Association Stockholm, Queen Silvia Jubilee Fund, Solstickan Foundation, PRIMA Child and Adult Psychiatry, the Pediatric Research Foundation at Astrid Lindgren Children’s Hospital, the Swedish Foundation for Strategic Research, Jerring Foundation, the Swedish Order of Freemasons, Kempe-Carlgrenska Foundation, Sunnderdahls Handikappsfond, The Jeansson Foundation, and EU-AIMS (European Autism Intervention), with support from the Innovative Medicines Initiative Joint Undertaking (Grant Agreement No. 115300), the resources of which are composed of financial contributions from the European Union’s Seventh Framework Programme (grant FP7/2007–2013), from the European Federation of Pharmaceutical Industries and Associations companies’ in-kind contributions, and from Autism Speaks. It was also supported by a new IMI initiative–EU AIMS-2-TRIALS. The Swedish Twin Registry is managed by Karolinska Institutet and receives funding through the Swedish Research Council under the Grant No 2017-00641. The funding organization had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication.
Ethics approval and consent to participate
This study was approved by the National Swedish Ethical Review Board and written and oral informed consent was obtained from all participants after the nature of the study procedure had been fully explained.
Consent for publication
All subjects involved in both methods have provided informed consent for publication.
The authors declare no conflict of interest related to this study. Bölte discloses that he has in the last 3 years acted as an author, consultant or lecturer for Medice and Roche. He receives royalties for textbooks, diagnostic and intervention tools from Hogrefe, Kohlhammer and UTB. Bölte is shareholder in SB Education and Psychological Consulting AB and NeuroSupportSolutions International AB.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Results of the intraclass correlation coefficient based on some of the included measurements, consistency agreement 2-way mixed effects model.
Within-pair association between autistic traits and perinatal or postnatal risk factors in two example pairs.
Partial Correlation Coefficients between the study variables, adjusted for age and sex.
Results from within-pair linear regression analyses with birthweight difference in gram as main predictors of trait Neurodevelopmental conditions.
About this article
Cite this article
Isaksson, J., Ruchkin, V., Aho, N. et al. Nonshared environmental factors in the aetiology of autism and other neurodevelopmental conditions: a monozygotic co-twin control study. Molecular Autism 13, 8 (2022). https://doi.org/10.1186/s13229-022-00487-5
- Early adversities
- Non-shared environment
- Risk factors