Abstract
Importance Recent research has demonstrated that domains of social determinants of health (SDOH) (e.g., air pollution and social context) are associated with psychosis. However, SDOHs have often been studied in isolation.
Objective To identify distinct exposure profiles, estimate their associations with persistent distressing psychotic-like experiences (PLE), and evaluate whether involvement with physical activities partially explains this association.
Design, Setting, and Participants This population-based study used data from the Adolescent Brain and Cognitive Development (ABCD) Study. Participants were recruited from 22 US sites between September 2016 and January 2022. Data from baseline and three follow-ups were included.
Exposures Area-level geocoded variables spanning various domains of SDOH, including socioeconomic status (SES), education, crime, built environment, social context, and crime, were clustered using a self-organizing map method to identify exposure profiles.
Main Outcomes and Measures Persistent distressing PLE was derived from the Prodromal Questionnaire-Brief Child Version across four years. Generalized linear mixed modeling tested the association between exposure profiles and persistent distressing PLE as well as physical activities (i.e., team and individual sports), adjusting for individual-level covariates including age, sex, race/ethnicity, highest level of parent education, family-relatedness, and study sites.
Results Among 8,145 participants (baseline mean [SD] age, 9.92 [0.63] years; 3,868 (47.5%) females; 5,566 (68.3%) White, 956 (11.7%) Black, 159 (2.0%) Asian, and 1,480 (18.4%) Hispanic participants), five exposure profiles were identified. Compared to the reference Profile 1 (suburban affluent areas, 2521 children, 30.9%), Profile 3 (rural areas with low walkability and high ozone; 1459 children, 17.9%; adjusted OR: 1.34, 95% CI: 1.09—1.64) and Profile 4 (urban areas with high SES deprivation, high crime, and high pollution; 715 children, 8.8%; adjusted OR: 1.40, 95% CI: 1.08—1.81), were associated with persistent distressing PLE. Team sports mediated 6.14% of the association for Profile 3.
Conclusion and Relevance This study found that neighborhoods characterized by rural areas with low walkability and urban areas with high socioeconomic deprivation, air pollutants, and crime were associated with persistent distressing PLE. Further research is needed to explore the pathways through which different environmental factors may impact the development of psychosis.
Introduction
Psychotic-like experiences (PLEs), also referred to as subclinical psychotic symptoms or psychotic experiences, are odd or unreal perceptions, thoughts, or beliefs. 1 It is one of the earliest manifestations of psychotic disorders and can be common among children. 2,3 Although PLEs have been considered mild expressions of psychosis liability, PLEs that are persistent and distressing may be more indicative of psychopathology later on in life with greater functional impairment, greater cognitive impairment, and increased mental health service utilization. 4,5 Given the high prevalence of PLE in early adolescence and its clinical significance to later risk for psychosis and general psychopathology, exploring the risk factors and mechanisms underlying persistent and distressful PLE is crucial to improve psychiatric risk assessments and more effective prevention strategies in children and adolescents.
It has been long known that urban upbringing is a major risk factor for developing psychotic disorders, including schizophrenia. 6–8 Under the social determinants of health model (SDOH) for psychosis, certain environmental and social conditions (e.g., air pollutants, social fragmentation, and socioeconomic deprivation) may impact downstream psychosocial and biological processes and, in turn, lead to psychosis. 9,10 Early exposure to air pollution has been shown to be associated with increased odds of psychotic-like experiences in adolescence, potentially through inflammation and excessive oxidative stress. 11–13 The Social Disorganization Theory posits that the wider community plays an important role in shaping social interactions, whether through sports or extracurricular activities, which may be crucial for mental health, especially for children. 14,15 For example, youth living in neighborhoods that were more walkable, had less crime, and greater social cohesion participated in more physical activities, particularly team sports. 16–18 And greater involvement with team sports, as opposed to individual sports, was shown to be associated with reduced risk for future psychopathology, including persistent distressing PLE. 19–21 Perhaps structural factors may facilitate youth coming together to play team sports, and certain aspects of these social ties may be relevant to offset the risk of psychosis. 21
Prior literature has largely focused on analyzing each neighborhood characteristic as it relates to psychosis separately. 22–24 Existing measurement tools, such as the Area Deprivation Index (ADI) and Child Opportunity Index (COI), have been limited in their ability to capture the multifaceted impacts of various neighborhood characteristics simultaneously while also understanding individual contributions of each neighborhood characteristic. These indices are constructed from an arbitrarily chosen set of SDOH variables, and little is currently known about the combinations of these environmental and social factors as well as the mechanisms through which they lead to psychosis. Without such information, it may be difficult to devise effective, targeted policies and interventions.
This study aimed to fill this knowledge gap by using innovative methodologies that capture and analyze SDOH’s multidimensional nature and uncover its associations with distressing PLE and physical activities. We clustered area-level exposures based on various indices of air pollutants and psychosocial factors into exposure profiles. Then we investigated whether these exposure profiles would be associated with persistent distressing. We also tested team sports as a potential mediator of these associations. We hypothesized that certain urban exposure profiles would be associated with persistent distressing PLE and that team sports would partially explain this association.
Methods
The data were collected from a population-based sample of 9-to-10-year-olds in the Adolescent Brain and Cognitive Development Study (ABCD) 5.0 release, which included visits collected between September 1, 2016, and January 15, 2022. 25 Data from baseline and three follow-up time points were included.
The ABCD Study is a nationwide longitudinal study of brain, behavioral, and cognitive development in adolescents and involves 22 research sites throughout the United States, with more than 11,868 children recruited at baseline. Within these sites, public, private, and public charter schools within a 50-mile radius of the data-collecting site were randomly selected for recruitment. 26 General inclusion and exclusion criteria for the ABCD study have been described elsewhere. 26–28 Informed written consent was obtained from children and their parents, and ethical approval was obtained from each site’s institutional review boards. We followed the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) reporting guidelines.
Participants
This study included 8,145 participants based on the availability of environmental exposures, outcome variables, and sociodemographic covariates. eFigure 1 shows a flowchart of included and excluded participants. eTable 1 compares those included and excluded based on missing data.
Sociodemographic and clinical characteristics
Sociodemographic characteristics, including age, sex, race, ethnicity, parental education, household income, and family history of psychosis, were collected through parent reports and interviews during baseline assessment. Race and ethnicity were aggregated into 5 categories: non-Hispanic White, non-Hispanic Black, non-Hispanic Asian, non-Hispanic other races, and Hispanic. High parental education was defined as having at least one parent or caregiver who obtained a bachelor’s degree or higher. The income-to-needs ratio was calculated by the median value of the income band divided by the federal poverty line for the respective household size. 29 A value greater or less than one would denote above or below the poverty threshold, respectively. Family history of psychosis in first-degree and second-degree relatives was assessed using the parent-rated Family History Assessment Module Screener. 30
Area-level exposures
Area-level data were derived from participants’ primary home addresses at baseline, which were geocoded to retrieve information at the census tract level. These geospatial location data were then linked to external environmental constructs, which were part of the ABCD 5.0 release. 25 Exposure variables were used to construct five domains, including the Area Deprivation Index (ADI), 31,32 Child Opportunity Index 2.0 (COI), 33 Crime, 34 Environmental Quality, 35,36 and Social Vulnerability Index (SVI), 37 based on their respective sources and prior literature. 25 A comprehensive list of variables included in our study, exposure definition, and years measured are described in Table 1 and eTable 2. Distance to major roadways, population density, and percentage of households without a car, were used to classify the exposure profiles as urban, suburban, and rural as done in prior literature. 38,39
A. Study Sample Characteristics. B. Area-level Deprivation Index. C. Child Opportunity Index. D. Crime. E. Environmental Quality. F. Social Vulnerability Index. G. Urbanicity Features. H. Physical Activities. I. Outcome Variable. a
Persistent distressing psychotic-like experiences (PLE)
The Prodromal Questionnaire-Brief Child Version (PQBC) was completed by participants to assess psychotic-like experiences (PLE) over the past month. 40 Initially, participants were categorized as either having experienced distressing PLE (rating of at least one PLE>= 3 on a five-point distress scale) or not having experienced distress, following criteria based on prior literature on distressing PLE. 4 Participants with distressing PLEs were grouped into persistent and non-persistent categories, with “persistent” indicating distress at two or more follow-ups out of 4 waves (baseline, 1-year, 2-year, and 3-year). This binary categorization was based on prior literature demonstrating such persistency (as opposed to transiency) was specifically associated with a polygenic risk score for schizophrenia, 21 greater cognitive deficits, worse social functioning, and more health service utilization4 as well as a greater risk for future diagnosis of schizophrenia. 41 Because psychotic symptoms may develop throughout adolescence, we conducted a sensitivity analysis and redefined persistent distressing PLE to incorporate individuals who reported distressing PLE only at the last follow-up.
Physical activities
Physical activities were derived from the Sports and Activities Involvement Questionnaire (SAIQ). Physical activities were further subcategorized into team and individual sports (eTable 3). The study focused on the level of involvement in 23 physical activities, with caregivers reporting (1) the time spent per session in minutes, (2) the number of days per week of participation, and (3) the number of months per year of engagement. The mean participation hours per week for each sport endorsed in the past 12 months were calculated using the formula: average hours per week per sport in the past year = (time spent × days per week × months per year) / 12 / 60. 42
Statistical analysis
The self-organizing map (SOM) is an unsupervised machine learning technique for multi-dimensional data reduction and visualization, and identifying profiles of correlated characteristics (e.g., air pollutants and socioeconomic deprivation). SOM aims to map complex, multi-dimensional data into a simpler, easy-to-visualize format where multiple exposures are grouped into an exposure profile. SOM was used to cluster 29 area-level exposure profiles in this study. See Table 1 and eTable 2 for a complete list of exposure variables. The optimal number of profiles was determined using within-cluster sum of squares, between-cluster sum of square statistics, and visual inspection of the exposure profile star plot. 43 These methods seek to maximize homogeneity within profiles and heterogeneity between profiles. The reference profile had the highest socioeconomic status and lowest hazardous environmental exposure (e.g., air pollution). Once participants were assigned to an exposure profile, unadjusted and adjusted logistic mixed models were used to estimate the association between exposure profiles, modeled as a categorical exposure, and persistent distressing PLE as well as physical activities (i.e., team and individual sports). The model adjusted for age, sex, family history of psychosis, race and ethnicity, parental education, income-to-needs ratio, as individual-level fixed effects, and family groups and recruiting sites as random effects.
We then assessed whether team and individual sports would mediate the association between exposure profile and persistent distressing PLE. Mediation analysis was only conducted for exposure profiles that were statistically significantly associated with persistent distressing PLE by modeling each significant exposure profile as a dichotomous exposure variable in the mediation analysis. Detailed definitions and methodology for the mediation analyses were provided in the eMethods.
Subsequently, we used another method, weighted quantile sum (WQS) regression, to understand the overall effect of the exposure mixture and how much of the total mixture effect on the outcome is explained by each area-level characteristic within the mixture. WQS is a method that estimates the effect on the outcome of increasing all exposures simultaneously by one quantile, as well as the weighted contribution of each exposure variable to the association between the overall exposure mixture and the outcome while accounting for the complex correlation structure among exposures. 44–46 All SOM analyses were performed using code within the ECM package (https://github.com/johnlpearce/ECM). The other analyses used the following R packages, lme4, 47 lmerTest, 48 mediation, 49 gWQS. 50 All analyses were conducted using R version 4.2.1. Statistical significance was determined using an alpha-level of 0.05.
Results
Descriptive statistics
The study included 8,145 participants aged 9 to 10 years at baseline followed until 13 to 14 years, with 3,868 (47.5%) females, 5,566 (68.3%) non-Hispanic White, 956 (11.7%) non-Hispanic Black, 159 (2.0%) non-Hispanic Asian, and 1,480 (18.4%) Hispanic participants. Among the participants, 5,401 (66.3%) had parents with at least a bachelor’s degree and 668 (8.2%) had a family history of psychosis. Across the four years, 1605 (19.7%) participants had persistent distressing PLE (Table 1). All individual neighborhood-level characteristics correlate with each other in the expected direction (eFigure 2).
Five distinct exposure profiles
Five exposure profiles were identified upon visual inspection. Each participant was assigned to one of the five profiles, and their demographics stratified by the exposure profiles are shown in Table 1. Figure 1A presents a detailed description of each profile.
Associations between exposure profiles and persistent distressing PLE
Profile 3, characterized by rural areas with low walkability and high ozone, and Profile 4, characterized by urban areas with high ADI, crime, and pollution, were significantly associated with persistent distressing PLE compared to the reference Profile 1 (Profile 3: adjusted odds ratio (OR): 1.34, 95% CI: 1.09—1.64, p=.006; Profile 4: adjusted OR: 1.40, 95% CI: 1.08—1.81, p=.01) (Table 2, Figure 1B). Sensitivity analyses of (1) an alternative version of persistent distressing PLE and (2) only including participants who lived in their addresses for more than one year showed consistent results (eTable 4). Variance inflation factors (VIF) for exposure profiles and individual-level covariates ruled out multicollinearity in our models (eTable 5).
Mediation of the relationship between Profile 3 and persistent distressing PLE by team sports
Team sports explained 6.14% of the relationship between Profile 3 and persistent distressing PLE (indirect effects adjusted β: 0.002, 95% CI: <.001—.005, p =.013) (Figure 2). However, team sports did not significantly mediate the association between Profile 4 and persistent distressing PLE (Figure 2). Total physical activities or individual sports did not significantly mediate the association between Profile 3 or Profile 4 and persistent distressing PLE (eTables 6 and 7).
Individual exposure contribution assessed by WQS regression weights
The positively constrained model demonstrated that the overall effect of increasing all exposures simultaneously by one decile corresponded to 12% greater odds for endorsement of persistent distressing PLE (adjusted OR: 1.12, 95% CI: 1.02—1.23, p=.03). The tree plot presented in Figure 3 demonstrates the contribution of each exposure to a higher odds of persistent distressing PLE with greater total crime having the highest mean weight at 0.130, followed by lower walkability at 0.106, greater ethnoracial minority concentration at 0.090, higher ozone at 0.085, and greater industrial pollutants at 0.077 (eTable 8). The negatively constrained model was not significant (adjusted estimate: < .01, 95% CI: -0.10—0.09, p=.95) (eTable 9).
Discussion
In this study, we identified five unique exposure profiles and investigated their associations with persistent distressing PLE and physical activities. We found two exposure profiles to be associated with greater odds of developing persistent distressing PLE. Rural areas that were less walkable with high ozone (Profile 3) and urban areas with high socioeconomic disadvantage, crime, and pollution (Profile 4) were both uniquely associated with greater odds for persistent distressing PLE even after adjusting for individual-level covariates. It is likely that a combination of factors, including financial stress, food insecurity, limited opportunities for physical activities, and the presence of air pollutants, might heighten the likelihood of PLE. 13,51,52 Chronic stress has been shown to disrupt the hypothalamic-pituitary-adrenal axis and elevate cortisol levels, which has been shown to be associated with neuroanatomical changes linked to psychotic illnesses. 53–55 Subsequently, our mediation analyses suggested that there may be distinct biopsychosocial pathways through which the environment may play a role in the development of psychosis. For example, the degree of involvement in team sports (but not individual sports) modestly mediated the association between Profile 3 and PLE by more than 6%. However, none of the subcomponents of physical activities mediate such relationships for Profile 4. Although the association between Profile 3 and PLE was mainly driven by other factors rather than team sports, it is possible that involvement in team sports (e.g., pick-up basketball) and perhaps other informal social interactions may be less likely to occur in more rural and less walkable places. 56 This lack of social engagement with peers may play a role in the future risk of psychotic experiences. 57,58
Our findings demonstrating that the association between Profile 4 and PLE was not explained by physical activity suggest that there may be an alternative mechanism through which urban areas characterized by high crime and pollution may be associated with PLE. It is possible that threat from exposure to crime, 59 air pollution, 60 and/or financial stress61 may play a more important role for this subgroup of youth with these specific combinations of exposures. In fact, recent literature has suggested that deprivation and threat may be two environmental factors that underlie partially distinct biological pathways in the development of psychosis. 62,63 Our findings may partially reflect this dimensional model of adversity in psychosis, and it is possible that different environmental factors may impact psychosis through distinct mechanisms.
In addition to the findings using SOM, the WQS regression model also points to similar factors that may drive the development of PLEs. High crime rates and low walkability carried the greatest weight in driving the effects of persistent distressing PLE, which align with high crime rates in Profile 4 and low walkability in Profile 3. Prior research has shown that living in high-crime neighborhoods was associated with subclinical psychotic symptoms, including suspiciousness and paranoia among help-seeking adolescents, 64,65 as well as a higher incidence of first-onset schizophrenia. 59 Although walkability has not been widely studied in psychosis, places with less walkability tend to have less access to community services and recreational centers. 16,66,67 Access to these resources may have downstream effects on exercise and involvement with team sports, which has been shown to be inversely associated with psychopathology and PLE. 19,20,68,69
One of the key strengths of this study is our approach to identifying multidimensional exposure profiles as measured by area-level characteristics, as opposed to traditional single-indexing methods. In addition, we used another exposure mixture method (i.e., WQS) to enhance the rigor of our findings and demonstrated consistent factors that may be driving the association with PLE. Future prospective studies should further investigate the biopsychosocial mechanisms through which neighborhood characteristics, including crime and walkability (etc. types of crime, park access, street connectivity, mixed land use) 70 may drive the development of psychotic disorders.
Limitation
This study has several limitations. First, we excluded several participants due to missing data. Participants excluded were from lower SES households, which may impact the findings’ generalizability. Second, PLEs were self-reported, and physical activities were parent-reported, which may be biased. Third, this study did not test whether there was a sensitive period in which environmental factors may have a more pronounced effect on long-term psychosis risk. Prospective studies should collect this data and analyze whether environmental factors at various developmental periods may differentially impact psychosis risk.
Conclusion
This study identified five exposure profiles, of which two were associated with persistent distressing psychotic-like experiences among children and adolescents across four years. These two exposure profiles were characterized by (1) rural areas with low walkability with high ozone (Profile 3), and (2) urban areas with high socioeconomic deprivation and high pollutants (Profile 4). Moreover, we found that less involvement with team sports partially mediated the positive association between Profile 3 (but not Profile 4) and greater odds of developing persistent distressing psychotic-like experiences. These findings suggest that there may be distinct mechanisms through which environmental factors may impact psychosis. Future studies should further explore the mechanisms through which living in certain neighborhoods may confer risk for the development of psychosis.
Data Availability
All data from the Adolescent Brain Cognitive Development (ABCD) Study (https://nda.nih.gov/abcd/request-access) are made available to researchers from universities and other institutions with research inquiries following institutional review board and National Institute of Mental Health Data Archive approval.
Key Points
Question
What are the underlying complex patterns of the environmental and social determinants of health, and what are their associations with persistent distressing psychotic-like experiences among children in the US?
Findings
Two distinct exposure profiles characterized by rural areas with low walkability and urban areas with high socioeconomic deprivation, air pollutants, and crime were associated with persistent distressing psychotic-like experiences.
Meaning
Several adverse environmental and social factors were jointly associated with psychotic-like experiences. Further research is needed to investigate the mechanisms through which these factors impact the development of psychosis.
Disclosures
There are no conflicts to disclose.
Conflicts of Interest
There are no conflicts of interest for any authors concerning the data or the study.
Data Sharing Statement
The ABCD study anonymized data are released annually and are publicly available via the NIMH Data Archive (NHA). All data from the Adolescent Brain Cognitive Development (ABCD) Study (https://nda.nih.gov/abcd/request-access) are made available to researchers from universities and other institutions with research inquiries following institutional review board and National Institute of Mental Health Data Archive approval.
Acknowledgments
This work was supported by the National Institute of Mental Health (NIMH K23-MH129684; Dr. Ku), the National Institute on Aging (NIA R01AG079170; Huels/Wingo), the National Institute of Mental Health (NIMH R01 MH129855; Risk), the Emory Constructive Collision Grant (Ku/Risk/Huels), and Emory HERCULES Pilot Award (Ku/Risk/Huels). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health, the National Institute of Mental Health, or Emory University.
Footnotes
The title has been revised. Figure 2 has been updated.