Introduction
Loneliness is a distressing emotional state that arises from the discrepancy between one's perceived and desired levels of social connection (Perlman & Peplau, Reference Perlman, Peplau, Gilmour and Duck1981). As such, perceived loneliness can exist for individuals who are not socially isolated (Cacioppo et al., Reference Cacioppo, Hawkley, Ernst, Burleson, Berntson, Nouriani and Spiegel2006). Loneliness is a major risk factor for psychological disturbance and poor health outcomes in adulthood. In population-based studies, loneliness in adolescents has been prospectively associated with social anxiety (Maes et al., Reference Maes, Nelemans, Danneel, Fernández-Castilla, Van den Noortgate, Goossens and Vanhalst2019; Teo, Lerrigo, & Rogers, Reference Teo, Lerrigo and Rogers2013) and depression (Qualter, Brown, Munn, & Rotenberg, Reference Qualter, Brown, Munn and Rotenberg2010).
Loneliness is relatively common among both children and adults. For example, between 11% and 27% of children (aged 10–15 years) reported feeling lonely in the UK (Snape & Manclossi, Reference Snape and Manclossi2018; Yang, Petersen, & Qualter, Reference Yang, Petersen and Qualter2020), whereas more than 40% of Americans of all ages reported feeling lonely (Cacioppo, Grippo, London, Goossens, & Cacioppo, Reference Cacioppo, Grippo, London, Goossens and Cacioppo2015). Perceived loneliness is also associated with a substantial increase in the risk of premature mortality, independent of income, education, sex, and ethnicity (Cacioppo, Cacioppo, Capitanio, & Cole, Reference Cacioppo, Cacioppo, Capitanio and Cole2015; Holt-Lunstad & Smith, Reference Holt-Lunstad and Smith2016; Miller, Reference Miller2011). A longitudinal analysis of four nationally representative US samples showed the influence of social isolation on several biomarkers of cardiovascular heart disease including hypertension, body mass index, and waist circumference and inflammation (High-sensitivity C-reactive protein) across the lifespan (Yang et al., Reference Yang, Boen, Gerken, Li, Schorpp and Harris2016). The magnitude of this effect is comparable to that of smoking and obesity or physical inactivity.
Much less is known about the long-term effects of childhood loneliness. A recent meta-analysis of 102 cross-sectional studies of childhood loneliness found associations with social anxiety (M age below 21 years) (Maes et al., Reference Maes, Nelemans, Danneel, Fernández-Castilla, Van den Noortgate, Goossens and Vanhalst2019). The same meta-analysis looked at 10 longitudinal studies and reported small but positive associations between loneliness and social anxiety symptoms both within and across time, and across childhood and adolescence. To our knowledge, there is one longitudinal study on the association between childhood loneliness and depression in pre-school and school-aged children (Qualter et al., Reference Qualter, Brown, Munn and Rotenberg2010). No studies to date have addressed the key question of whether childhood loneliness is associated with a long-term disruption in mental health that extends into adulthood.
Parents and children often disagree about the presence and severity of child symptoms and psychopathology (Achenbach, McConaughy, & Howell, Reference Achenbach, McConaughy and Howell1987; De Los Reyes et al., Reference De Los Reyes, Youngstrom, Pabon, Youngstrom, Feeny and Findling2011). Thus, parental and child ratings may capture overlapping but largely distinct information about a child's experience of loneliness and associated risk for psychiatric disorder outcomes. Loneliness is, by definition, a subjective, internal state of mind. Yet, it is reasonable to expect that raters who know the child well (e.g. parents or friends) may have access to information about the individual's loneliness and therefore provide accurate judgments of internal traits.
To address these gaps, we conducted a large population-based study of the association between childhood loneliness and psychiatric disorders in adulthood. We focused on loneliness measures assessed eight times from children and a parent from ages 9 to 16 years. Participants were then followed up four times in adulthood between ages 19 and 30 years, and assessed for anxiety, depression, and substance use symptoms and disorders. The present study had three main aims. First, we aimed to examine whether repeated measures of loneliness in childhood are associated with a broad range of psychiatric problems that include anxiety, depression, and substance use disorders in adults. We also examined the associations between the groups of childhood loneliness trajectories and these psychiatric outcomes. Second, we aimed to examine whether the associations between childhood loneliness and adult psychiatric disorders persist after accounting for other co-occurring adverse childhood experiences. Finally, we examined whether parent-reported loneliness predicts adult outcomes similar to child-reported loneliness.
Method
Participants
This report follows the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guideline for cohort studies (Von Elm et al., Reference Von Elm, Altman, Egger, Pocock, Gøtzsche and Vandenbroucke2007). Data were drawn from the Great Smoky Mountains Study (GSMS), a longitudinal, community-representative cohort in which participants were followed up from age 9 years into adulthood. The GSMS, which enrolled children in 11 predominantly rural North Carolina counties, was originally designed to estimate the prevalence of mental illness and service (Costello, Mustillo, Erkanli, Keeler, & Angold, Reference Costello, Mustillo, Erkanli, Keeler and Angold2003). Initially, three cohorts of children, aged 9, 11, and 13 years, were recruited from a pool of approximately **12 000 children using a two-stage sampling design, resulting in 1420 participants [630 girls (weighted percentage, 49.0%)] (Costello et al., Reference Costello, Mustillo, Erkanli, Keeler and Angold2003). Sampling weights were applied to adjust for differential probability of selection. An ascertainment figure appears in Fig. 1 in the online Supplementary material, and the original study articles (Copeland, Angold, Shanahan, & Costello, Reference Copeland, Angold, Shanahan and Costello2014; Costello et al., Reference Costello, Angold, Burns, Stangl, Tweed, Erkanli and Worthman1996, Reference Costello, Mustillo, Erkanli, Keeler and Angold2003) provide additional detail on sampling and derivation of sample weights.
Annual assessments were completed until participants were 16 years old and then they were assessed again at ages 19, 21, 25, and 30 years, for a total of 11 230 total assessments from January 1993 to December 2015. Across all assessments, 83.0% of possible interviews (11 230 of 13 530) were completed. Race/ethnicity was determined based on parent report. Before all interviews, parents and children signed institutional review board-approved informed consent and assent forms. All procedures and protocols for the present study were approved by the Duke University Institutional Review Board.
Measures
Childhood loneliness
Questions about perceived loneliness were collected as part of self or parental-report interviews using the Child and Adolescent Psychiatric Assessment (CAPA) from ages 9 to 16 (Angold & Costello, Reference Angold and Costello2000). The CAPA defines loneliness as ‘a feeling of being alone and/or friendless, regardless of the justification for the feeling; total daily duration of at least 1 h’. The CAPA focuses on the 3 months immediately preceding the interview as its primary period. To assess this item, participants and their parents were asked first ‘Do you think s/he feels lonely?’ and then ‘Sometimes children feel that they have no one who would help them. Does s/he ever feel like that?’. These two questions are followed by a series of secondary prompts, if necessary, to clarify whether the child met the operational definition of lonely. Secondary prompts comprised of items such as: ‘How often is that?’, ‘When was the last time?’, ‘Does s/he feel cared for by friends?’, ‘Does s/he feel lonely even though s/he has some friends?’. More details on how loneliness is defined and assessed through primary and follow-up prompts can be found in online Supplementary Fig. 2.
These items index both a perceived quantitative lack of contacts in one's social network and a perceived qualitative deficit in existing relations. As such, this construct is best considered perceived social isolation (feeling lonely) as opposed to objective social isolation (being lonely). Different dimensions of social functioning (e.g. subjective loneliness, network quality, and network size) have been found to have different associations with health (Victor, Scambler, Bond, & Bowling, Reference Victor, Scambler, Bond and Bowling2000). Subjective loneliness has been linked to psychological consequences e.g. Weiss (Reference Weiss1973) and Peplau and Perlman (Reference Peplau, Perlman, Mark and Glenn1979). Other studies of loneliness have used a similar construct (Qualter et al., Reference Qualter, Brown, Munn and Rotenberg2010, Reference Qualter, Brown, Rotenberg, Vanhalst, Harris, Goossens and Munn2013).
Loneliness is scored on a three-point rating scale 0 = ‘not at all’, 1 = ‘the interviewee definitely feels lonely in a way that interferes with at least two activities and is uncontrollable’, and 2 = ‘the interviewee feels lonely most of the time’. For the current study, loneliness was scored dichotomously (‘0’ v. ‘1’ and ‘2’) and we aggregated the loneliness measure between ages 9 and 16 into a single measure of childhood loneliness.
Adult psychiatric disorders and functioning
All outcomes except where noted were assessed using the Young Adult Psychiatric Assessment (Angold et al., Reference Angold, Cox, Prendergast, Rutter, Simonoff, Costello and Ascher1999), an upward extension of the CAPA interview administered to the participants. The assessment of adult psychiatric disorders resembled that of childhood disorders, but with only self- (and not parent-) reports. A 3-month ‘primary period’ was selected because longer recall periods are associated with forgetting and recall bias (Copeland et al., Reference Copeland, Angold, Shanahan and Costello2014; Costello et al., Reference Costello, Angold, Burns, Stangl, Tweed, Erkanli and Worthman1996, Reference Costello, Mustillo, Erkanli, Keeler and Angold2003). Disorders included any DSM generalized anxiety disorder, social phobia, post-traumatic stress disorder, panic disorder, agoraphobia, depressive disorder, and substance use disorder (including, nicotine, alcohol, cannabis, and other illicit drugs). For the current study, we examined anxiety, depression, as well as substance use disorders. Each psychiatric disorder measure between ages 19 and 30 is aggregated into a single measure of adult psychiatric status. The participant was positive for diagnosis if criteria were met at any adult observation. In addition, we examined symptoms scores of the total anxiety, depressive, and substance use symptoms.
Covariates
Parent and child characteristics examined as potential confounders are depicted in the Appendix in the online Supplementary material. These included the following: sex of the child, rural v. urban area, family hardships (including low socio-economic status, single parent, change in parent structure, maltreatment, and depression of the mother), and childhood psychiatric comorbidities (including anxiety, depression, and disruptive behavior). For psychiatric symptoms, the CAPA focuses on the 3 months immediately preceding the interview to minimize recall bias. Scoring programs written in SAS statistical software (SAS Institute Inc.) combine information about the date of onset, duration, and intensity of each symptom to create DSM diagnoses. Test–retest reliability and validity of the CAPA diagnoses are similar to other psychiatric interviews (Angold & Costello, Reference Angold and Costello1995, Reference Angold and Costello2000). Psychiatric disorders assessed included anxiety disorders, depression disorder, and substance use disorders. The categories of family hardships or childhood adversities were assessed at each observation. A full description of these variables is available in a previous publication (Copeland, Wolke, Angold, & Costello, Reference Copeland, Wolke, Angold and Costello2013). In our study, we aggregated childhood covariates between ages 9 and 16 into a single measure of each covariate.
Statistical analysis
First, we computed descriptive statistics for self- and parent-reported loneliness at different time points. Then, we tested prospective associations of childhood loneliness (between ages 9 and 16 years) reported by parent and child with adult anxiety, depression, and substance use diagnostic status (between ages 19 and 30 years) with separate weighted logistic regressions. Then, we tested prospective associations of childhood loneliness with adult anxiety, depressive, and substance use symptom scores with separate linear regressions. All analyses applied sampling weights; therefore, results provide estimates of the original representative population from which the sample was drawn.
As a follow-up to primary analysis, a latent class growth analysis was used to test for the associations of childhood loneliness trajectories with adult psychiatric disorders. Trajectories of loneliness were modeled in children from whom data were available for two or more time points (N = 1334). Models with two to five trajectories were assessed and compared. We used the full information maximum likelihood ratio to account for missing data. To assess model fit, we evaluated the bootstrap likelihood ratio test, the Bayesian information criterion, and the Akaike information criterion (Nylund et al., Reference Nylund, Asparouhov and Muthén2007). Subsequently, linear regression analyses were used to assess childhood loneliness trajectory in relation to adult outcomes.
Across all assessments, 83% of possible interviews were completed (online Supplementary Table S1). All 1420 participants were interviewed at least once in childhood (ages 9–16); 1260 participants (88.7%) had 3+ childhood observations. Of the total sample, 1334 (94.0%) were followed up at least once in adulthood at ages 19, 21, 25, or 30. All analyses were performed using SAS 9.4 software. All missing values of the potential confounding factors were imputed using multiple imputations. With the Markov Chain Monte Carlo multiple imputation techniques, 10 complete data sets were created. Multivariate analyses were performed on each imputed data set, and effect estimates were pooled (Schunk, Reference Schunk2008).
Results
Prevalence of loneliness and the associations with sociodemographic factors
Descriptive information is provided in Table 1. The prevalence of childhood/adolescent loneliness was 13.4%, meaning that >1 in 10 children reported feeling lonely at some point by age 16. Overall, the prevalence of loneliness differed by sex of the child (more common in girls than boys), but not by race/ethnicity and rural v. urban area. In our study, the within-individual correlation of loneliness measure over time was r = 0.42. Concurrent analyses showed that childhood loneliness was associated with a change in parent structure, as well as maternal depression, after adjustment for sex of the child. We observed no associations between childhood loneliness and low socio-economic status, single-parent family status, and child maltreatment.
Numbers denote children included in one or more analyses. Numbers are unweighted, and percentages are weighted.
a p value from binary logistic regression of childhood loneliness and childhood adversities outcome. The models [ORs] are adjusted for child sex.
Next, we tested the associations of child- and parent-reported loneliness with anxiety, depression, and substance use disorder within childhood. Child- and parent-reported loneliness was concurrently associated with childhood anxiety, depression, and substance use disorder (Table 2). The association of loneliness with psychiatric disorder outcomes was consistently stronger if loneliness was rated by the child than by the parent.
Numbers denote children included in one or more analyses. Values are frequencies for categorical (numbers are unweighted, and percentages are weighted). Means and standard deviations (M ± s.d.) for continuous measures.
a p value from binary logistic regression of childhood loneliness and childhood psychiatric disorder outcomes. The models [ORs] are adjusted for child sex and adversities.
Furthermore, we tested whether childhood loneliness was associated with the number of arguments with peers and teased or bullied during childhood (online Supplementary Table S1). Concurrent associations showed that childhood loneliness was associated with the number of arguments with peers and teased or bullied in childhood. In contrast, we observed no associations of childhood loneliness and frequency of contact with peers, confidante with family and peers, as well as shyness with peers. These associations and non-associations suggest that our measure best approximates subjective loneliness.
Childhood loneliness and adult psychiatric disorders
We then examined the associations of child- and parent-report loneliness and psychiatric disorders in adulthood adjusted for demographic variables and childhood adversities. As shown in Table 3, self-reported loneliness was associated with an anxiety disorder ( Odds ratio (OR) = 3.53, 95% Confidence interval (CI) 1.55–8.04, p = 0.002) but not with depressive and substance use disorders. Next, to test whether childhood loneliness was associated with adult psychiatric symptoms. Self-reported loneliness was associated with anxiety (B = 1.20, 95% CI 0.43–1.97, p = 0.002), and depressive symptoms scores (B = 0.76, 95% CI 0.27–1.25, p = 0.002), but not with substance use symptoms (B = 0.20, 95% CI, −0.14 to 0.54, p = 0.246). Effect estimates were modestly attenuated when we accounted for childhood psychiatric status (model 2).
Regression analysis of childhood loneliness and adult psychiatric outcomes. Values are frequencies for categorical (numbers are unweighted, and percentages are weighted) and means and standard deviations (M ± s.d.) for continuous measures. Model 1 is adjusted for child sex and childhood adversities. Childhood adversities include low family socioeconomic status, change in parent structure, mother depression. Model 2 is additionally adjusted for childhood psychiatric disorders. Child psychiatric disorders include anxiety, depression, and disruptive disorders.
OR indicate binary logistic regression coefficients for psychiatric disorders.
B statistics indicate linear regression coefficients for psychiatric symptoms.
a Model 2 is additionally adjusted for childhood psychiatric symptoms.
In fully adjusted models (model 2), we observed no associations between parent-reported loneliness and any measure of adult diagnostic status. Parent-reported loneliness was, however, associated with adult anxiety (B = 1.12, 95% CI 0.49–1.75, p = 0.001) and depressive (B = 1.12, 95% CI 0.49–1.75, p = 0.001) psychiatric symptoms.
Childhood loneliness trajectories and adult psychiatric disorders
We tested associations of childhood loneliness trajectories with adult psychiatric disorders and symptoms. Figure 1 illustrates the mean scores for the total sample of the three trajectory groups. The first class (N = 1260, 94.4%), termed ‘low’, consisted of children who reported no or very few feelings of loneliness over time. The second class (N = 65, 5.1%), labeled ‘moderate’ included children with a moderate increasing level of loneliness. Finally, there was a small group of children (N = 9, 0.5%), with slightly higher levels of lonely feelings than the moderate group. The three-class model was found to be the optimal model, with the lowest Bayesian information criterion scores and p values = 0.05 for the bootstrap likelihood ratio test (online Supplementary Table S2). Because of the sample size for the high group, it will be combined with the moderate group for all analyses.
As shown in Table 4, children in the moderate/high groups had higher levels of anxiety disorders than the children who were in the low group. There were no differences in depression and substance use disorder between groups. Children in the moderate/high groups had higher levels of anxiety and depressive symptoms than those in the low loneliness group. No differences in substance use symptoms between groups were found.
Regression analysis of childhood loneliness trajectories and adult psychiatric outcomes. Model 1 is adjusted for child sex and childhood adversities. Childhood adversities include low family socioeconomic status, change in parent structure, mother depression. Model 2 is additionally adjusted for childhood psychiatric disorders. Child psychiatric disorders include anxiety, depression, and disruptive disorders.
Odds ratio (OR) indicate binary logistic regression coefficients for psychiatric disorders.
B statistics indicate linear regression coefficients for psychiatric symptoms.
a Model 2 is additionally adjusted for childhood psychiatric symptoms.
Discussion
This prospective population-based study examined the associations of childhood loneliness and adult psychiatric disorders while carefully controlling other common childhood adversities and childhood psychiatric functioning. We highlight three key findings. First, childhood loneliness was prospectively associated with adult self-reported anxiety and depression but not substance use outcomes. Moreover, children exposed to persistently moderate and high levels of loneliness trajectories showed more symptoms of anxiety and depression. Second, these associations remained when we account for childhood-assessed adversities and psychiatric comorbidities (e.g. anxiety, depression, and conduct disorder). Third, the associations were stronger for self- than for parent-reported childhood loneliness, although there was evidence of cross-informant associations for psychiatric symptoms. Overall, our findings suggest that children's experience of loneliness is associated with increased risk for psychiatric disorders, and has the potential to have lifelong effects on one's social and emotional functioning.
Prior studies examining developmental trajectories of loneliness from childhood to early adulthood have indicated that between 3% and 22% of people experience persistent loneliness (Qualter et al., Reference Qualter, Brown, Munn and Rotenberg2010, Reference Qualter, Brown, Rotenberg, Vanhalst, Harris, Goossens and Munn2013; Vanhalst, Luyckx, & Goossens, Reference Vanhalst, Luyckx and Goossens2014; Yang et al., Reference Yang, Petersen and Qualter2020). With our data spanning 20+ years and multiple informants, we were able to extend previous findings across developmental periods indicating that loneliness experienced in childhood had particularly robust associations with adult self-reported anxiety and depression. Importantly, our results suggest that the associations of childhood loneliness with adult psychiatric outcomes were independent of childhood sociodemographic factors, adversities, and psychiatric functioning, indicating a unique contribution of loneliness to mental health outcomes later in life. No such associations were observed for adult substance disorders or symptoms.
The trajectory analysis largely confirmed the primary analysis with evidence that moderate to high loneliness is likely to affect anxiety and depression outcomes. Within these groups, loneliness tended to peak between ages 14 and 16 during early adolescence when peer groups and influences are taking on increasing importance. At the same time, it is notable that the rank order of the trajectories remained the same from childhood through adolescence suggesting the potential for early identification of children at risk. These results suggest that childhood loneliness may be a potential risk factor (Offord & Kraemer, Reference Offord and Kraemer2000) to emerging mental health problems and that age-appropriate interventions for loneliness may alleviate later suffering.
The finding that early childhood loneliness precedes mental health outcomes suggests that loneliness may be attributed to myriad processes including neural (Cacioppo, Capitanio, & Cacioppo, Reference Cacioppo, Capitanio and Cacioppo2014), neuroendocrine (Cacioppo et al., Reference Cacioppo, Cacioppo, Capitanio and Cole2015), and genetic mechanisms (Goossens et al., Reference Goossens, Van Roekel, Verhagen, Cacioppo, Cacioppo, Maes and Boomsma2015), as well as physiological stress state. The evolutionary theory of loneliness posits that loneliness increases the motivation to attend to and approach social stimuli for potential relief from the aversive state (e.g. like hunger, thirst, and pain promotes behavior change to increase the likelihood of the survival of one's genes) (Boomsma, Cacioppo, Muthén, Asparouhov, & Clark, Reference Boomsma, Cacioppo, Muthén, Asparouhov and Clark2007). The evolutionary theory of loneliness suggests that such experiences may contribute to (a) increased vigilance for social threats along with increased anxiety, hostility, and social withdrawal to avoid predation, (b) increased sleep fragmentation, (c) elevated vascular activity, increase extended periods of hypothalamo–pituitary–adrenal activations, and decreased gene expression and immunity to deal with potential assaults that may arise, (d) decreased impulse control (e.g. prepotent responding), and (e) increased depressive symptomatology (Cacioppo et al., Reference Cacioppo, Capitanio and Cacioppo2014). Tests of these mechanisms are beyond the scope of this study but may explain the long-term associations observed here.
Our study extended prior findings by using a multi-informant approach to address the potential for a bias when relying only on self-report data. Because loneliness is a subjective experience, and an internal state of mind, it is often examined using self-reports (Heinrich & Gullone, Reference Heinrich and Gullone2006). However, as close others (e.g. parents, teacher, or friends) can observe behavioral changes resulting from loneliness, they could provide additional information on children's loneliness (Cacioppo, Cacioppo, & Boomsma, Reference Cacioppo, Cacioppo and Boomsma2014). In our study, both self- and parent-loneliness ratings were associated with adult self-reported anxiety and depression outcomes, although higher associations were observed by self-reported loneliness. Self-ratings may be the best indicators of internal traits such as loneliness (Vazire, Reference Vazire2010), but parental reports also pick up on problems that suggest long-term distress. Discrepancies between informants might arise from substantial changes in individuals' social experiences and expectations across adolescent development (Qualter et al., Reference Qualter, Vanhalst, Harris, Van Roekel, Lodder, Bangee and Verhagen2015). These differences between self- and parent-report loneliness could also be explained attentional processes in adults, either directly or indirectly (e.g. mothers or fathers) (Cacioppo et al., Reference Cacioppo, Cacioppo, Capitanio and Cole2015).
The current study has several limitations. First, this study was representative of a rural area in the Southeast USA but not the US population. Second, loneliness was measured using a series of primary and secondary prompts but only coded into a single item. Insofar as the measure is limited, it may underestimate the effect of loneliness on adult outcomes. Next, this study did not include the potential influence of genetic variation on the young adult's mental health. Genetic and environmental determinants of loneliness (e.g. predispositional vulnerability and exposures to specific life experiences), could contribute to differences across individuals in which pathways operate (Goossens et al., Reference Goossens, Van Roekel, Verhagen, Cacioppo, Cacioppo, Maes and Boomsma2015). For example, McGuire & Clifford (Reference McGuire and Clifford2000) showed significant heritability and non-shared environmental influences for children's loneliness using adoption and twin studies (McGuire & Clifford, Reference McGuire and Clifford2000). The strengths of the present study are large samples, repeated assessments of loneliness across childhood and adolescence, multiple informants, and assessment of a broad spectrum of measured childhood covariates.
Our findings may have implications for future research and clinical practice. First, increased opportunities for social contact and social support, and improved social skills may reduce the risk of future psychiatric disorders in lonely children. However, such study needs to be targeted at the subjective experience of loneliness rather than merely at increasing objective social contacts. Overall, early interventions targeting children's maladaptive social cognitions may be an efficient way to alleviate such subjective feelings of loneliness. Such interventions would have to be implemented in a developmentally appropriate way given the social, cognitive, and emotional changes from childhood to adolescence.
In summary, our findings show that loneliness is relatively common and is observed by parents and children early in life. The question is whether this is a transient dysphoric state that affects current health only or it has the potential to compromise emotional health long term. The current study suggests that loneliness is not transient and the effects are long term. It is also important to remember that even if childhood loneliness is a risk factor for adult distress, most children with childhood loneliness do not experience adult distress. Future research should identify who is at risk for such long-term effects of loneliness and how this risk is propagated across significant developmental transitions.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S0033291721001422.
Acknowledgements
William E. Copeland was supported by a grant from the National Institute of Mental Health (R01MH117559, R01HD093651). Yllza Xerxa was supported by a grant from the Royal Netherlands Academy of Arts and Sciences (KNAW) Ter Meulen (KNAWWF/DA/2061/TMB367). The authors would like to thank the contribution of participating parents and their children for being part of the study for the past 25 years.
Conflict of interest
The authors declare no competing interest.