Abstract
Conceptualizations that distinguish systems-level stress exposures are lacking; the stimulation (lack of safety and high attentional demands), discrepancy (social exclusion and lack of belonging), and deprivation (SDD; lack of environmental enrichment) theory of psychosis and stressors occurring at the systems level has not been directly tested. Exploratory factor analysis was conducted on 3,207 youths, and associations with psychotic-like experiences (PLEs) were explored. Although model fit was suboptimal, five factors were defined, and four were consistent with the SDD theory and related to PLEs. Objective and subjective or self-report exposures for deprivation showed significantly stronger PLE associations compared with discrepancy and objective stimulation factors. Objective and subjective or self-report measures converged overall, although self-report stimulation exhibited a significantly stronger association with PLEs compared with objective stimulation. Considering distinct systems-level exposures could help clarify putative mechanisms and psychosis vulnerability. The preliminary approach potentially informs health policy efforts aimed at psychopathology prevention and intervention.
Keywords
Stress exposure has been widely implicated to play a causal role in the etiology of psychosis (Mayo et al., 2017; McEwen, 2004; Mittal & Walker, 2019; Pruessner et al., 2017; Shah & Malla, 2015). A wide array of work has examined individual-level stressors and their relation to psychosis risk. More recently, however, the field has looked toward examining systems-level factors, such as neighborhood features. A recent review synthesizing the literature on systemic environmental risk factors for psychotic disorders hypothesized that many previously explored exposures fall within three domains (Vargas et al., 2020). These hypothesized environmental risk domains include stimulation (systemic factors conferring lack of safety and high attentional demands), discrepancy (systemic factors conferring lack of belonging and social exclusion), and deprivation (systemic factors conferring lack of needed environmental enrichment). To date, relevant research on systemic environmental risk has focused primarily on adult populations, although late childhood and early adolescent periods could represent a highly informative time span for early detection and prevention efforts (Anniko et al., 2019; Felner et al., 1995; Zinzow et al., 2009).
Furthermore, although a body of literature supports the important roles for each of these three domains, they have yet to be directly tested together. The present investigation directly tested this theory in three stages: (a) first by examining whether risk factor exposures would separate into the hypothesized domains, (b) then relating the resulting domains to psychotic-like experiences (PLEs) to assess relevance to psychosis risk, and (c) finally, by comparing the magnitude of effects for observed associations by domain to explore differences in degrees of vulnerability. Examining systems-level stress by synthesizing these factors into environmental domains of influence is a critical priority. Despite the considerable research attention dedicated to psychotic spectrum disorders, the field also lacks a clear understanding of Environment × Liability interactions, particularly in the pediatric developmental periods. Given the challenging prognosis of the condition, a better understanding of causal factors and systems-level prevention efforts are paramount.
Psychotic disorders are chronic in nature, difficult to treat, and highly debilitating, constituting one of the top 15 leading causes of disability worldwide (GBD 2016 Disease and Injury Incidence and Prevalence Collaborators, 2017). After psychotic illness onset, confounds related to factors such as medication use and functional decline make it difficult to distill factors driving illness onset. Thus, assessment of associated symptoms or experiences on the psychosis spectrum, such as PLEs, provides a promising alternative for identifying factors relating to psychosis etiology (Van Os et al., 2009). PLEs, including experiences such as unusual beliefs, suspiciousness, and perceptual abnormalities, have been associated with pathogenic factors implicated in formal psychosis (Kelleher & Cannon, 2011; Morgan et al., 2009; Olin & Mednick, 1996; Orr et al., 2014; Papanastasiou et al., 2020; Yung et al., 2006). Furthermore, PLEs are experienced by 13% to 15% of children (Laurens et al., 2007; Poulton et al., 2000), and childhood experiences of PLEs have been shown to increase later risk for psychotic disorder onset (Kline et al., 2014; Poulton et al., 2000; Welham et al., 2009). Investigating PLEs in childhood provides an opportunity to understand environmental risk factors early in development before illness-related confounds. PLEs may also be of key importance during a critical developmental period in which prevention could be particularly effective.
Stress exposure is a central contributing factor in the development of psychotic disorders (Green et al., 2010; McLaughlin et al., 2012). Robust and varied evidence suggests a wide host of stressors can cumulatively contribute to psychosis risk (Mayo et al., 2017; McEwen, 2004; Mittal & Walker, 2019; Pruessner et al., 2017; Shah & Malla, 2015). The neural diathesis-stress model, for example, posits that the accumulation of stressors can comprise “multiple hits,” which, acting on a vulnerable system, can increase the likelihood of developing a psychotic disorder (Pruessner et al., 2017; Walker et al., 2008). Although this theory highlights the importance of conceptualizing stressors collectively, isolating qualitatively distinct stressors could also be uniquely informative to psychosis etiology. Work on trauma exposure, for instance, has illuminated differing neurodevelopmental consequences depending on the type of trauma (Gibson et al., 2016; McGrath et al., 2017; McLaughlin, Sheridan, & Lambert, 2014). Beyond individual-level factors such as trauma, distinct domains of systemic environmental risk factors may also contribute to vulnerability for developing a psychotic disorder. Systemic environmental risk factors have received less attention compared with individual-level factors (e.g., trauma, bullying, family environment). Conceptualizing systemic environmental factors into domains could ultimately aid in understanding the complex and multifaceted nature of psychotic illness presentation. Isolating the impacts of different dimensions of experience could be informative to psychotic disorder etiology. Indeed, the nonpsychiatric literature has successfully spearheaded conceptualizing neurodevelopment in this manner with individual-level stressors (McLaughlin, Sheridan, & Lambert, 2014). Despite this fact, the psychotic disorder literature has largely conceptualized stressors collectively, and less attention has been given to understanding specific kinds of exposures or systemic environmental factors.
Individuals operate within a larger environmental, social, and cultural context (i.e., structural or systems-level factors); thus, stress can also occur at the systems level (Bronfenbrenner, 1994; Glass & McAtee, 2006). Local structural characteristics (e.g., neighborhood socioeconomic status or cultural integration) could systematically affect well-being and risk for psychopathology. Increasing understanding of structural barriers to mental health has strong potential to inform health policy initiatives and prevention and intervention efforts at the societal level. To this end, our group has developed a literature-backed theoretical model of different types of structural exposures along with distinct intermediary mechanisms of impact and proposed relevant neural systems (Vargas et al., 2020)—the stimulation, discrepancy, and deprivation (SDD) model of psychosis. Although each hypothesized SDD domain benefits from support from the psychotic disorder literature, the distinctiveness of the domains has yet to be tested as they relate to the psychosis continuum. As a result, the degree to which each domain contributes to psychosis risk is unclear. Understanding each domain as it relates to degrees of risk for psychopathology could be immensely useful in identifying and prioritizing treatment targets. Furthermore, existing evidence in support of each domain has often homed in on adolescent or young adult populations. Determining whether structural, environmental risk exposures are also meaningful in earlier development (during childhood and early adolescence) in a way that is relevant to subthreshold psychotic symptoms would lend granularity to the understanding of environmental risk exposure across developmental periods.
The Adolescent Brain Cognitive Development (ABCD) data set is the largest investigation of brain development and child health in the United States (Garavan et al., 2018; Volkow et al., 2018). The study constitutes a nationally representative collaboration across 21 sites aiming to understand child and adolescent development. Previous investigations on the ABCD data set have indeed explored similar questions, lending unique insights in the process. Some of these have targeted systems-level features—Karcher and colleagues found urbanicity, deprivation, and lead exposure risk related to PLEs (Karcher et al., 2021) as well as systems-level environmental risk and neural features (Karcher et al., 2021; Marshall et al., 2020). At the individual level, adverse childhood experiences were also found to relate to PLEs (Karcher et al., 2020). Although these investigations have provided an invaluable perspective, several outstanding questions remain. First, although the existing evidence suggests that systems-level factors are influential, the theoretical understanding is more limited. It is unclear to what extent types of systemic stressors relate to psychosis vulnerability—examining differences in the magnitude of such associations would expand and inform the SDD theory, aiding understanding of degrees of impact. Furthermore, although it is clear that there are converging mechanisms through which systemic stressors could be influential, efforts to distinguish qualitatively distinct stressors have been limited. Distinguishing types of systemic stressors is a first and necessary preliminary step toward identifying and understanding how structural or systemic factors can contribute to stress, psychosis etiology, and symptoms.
In the current study, we used a nationally representative sample of youths 9 to 11 years old to further understand exposure to environmental stressors in relation to PLEs. The first aim was to directly test the SDD theory by exploring whether relevant items would load into factors consistent with the three hypothesized domains. The second aim was to determine whether the environmental stress domains would relate to PLEs, consistent with the SDD theory. The third aim was to then compare relative strengths of existing associations between environmental exposures and PLEs (to see whether certain exposures would show greater associations with psychosis risk than others). Finally, given investigations that have found divergences in self-reported compared with objective environmental and neighborhood measures (Gallagher et al., 2016; Hidalgo et al., 2015), a final exploratory aim sought to determine whether associations between objective, census-derived neighborhood metrics of environmental exposures and self-report measures indexing the same exposure would exhibit relations similar in magnitude.
Method
Participants
The ABCD data set includes a large representative sample of children 9 to 11 years old across 21 centers in the United States (for demographic information, see Table S1 in the Supplemental Material available online; Barch et al., 2018; Garavan et al., 2018). All centers obtained the parents’ informed consent as well as the children’s assent. Research procedures followed ethical guidelines laid out by respective institutional review boards (https://dx-doi-org-s.web.bisu.edu.cn/10.15154/1519171). The current sample included baseline data for participants who had available data for items in the final factor solution and PLE data. In the case that there were two or three siblings that completed the study, one child per family was randomly chosen for inclusion, which resulted in 1,020 participants being excluded. A total of 6,415 children were used to do the random sample split for exploratory factor analysis (EFA) and confirmatory factor analysis (CFA) because they had available data on all self-report items that were initially considered, and a subset of 6,072 children also had available data on “objective” neighborhood measures, PLEs, and self-report items used in the factor analysis. Group demeaning per site was conducted to account for possible effects of nesting within sites according to recommendations (Bear et al., 2011; Huang, 2016; Huang & Cornell, 2016). PLE analyses were conducted with site demeaned values.
PLEs
The Prodromal Questionnaire–Brief Child (PPS) version was used to assess psychotic-like experiences (Cicero et al., 2019; Karcher et al., 2018; Loewy et al., 2011). The 21-item self-report questionnaire was previously validated in the ABCD study sample (Karcher et al., 2018). The questionnaire asked participants about specific PLEs that were endorsed with a binary response (yes/no). Participants also indicated whether there was distress related to endorsed symptoms on a five-item Likert scale. Consistent with prior research, PLE scores accounting for distress were calculated such that the total number of endorsed symptoms was weighed by the level of distress (0 = zero endorsement, 1 = endorsement without distress, 2–6 = endorsement with incremental distress levels; Karcher et al., 2018; Loewy et al., 2011). PLE scores accounting for distress were used for all analyses. For the current sample, the average item endorsement accounting for distress was 6.49 (SD = 10.77). Of our total sample, 57.4% of participants had a rating of more than zero on the PPS.
Self-report questionnaires
Self-report measures relevant to theoretical interest in the three domains of deprivation, discrepancy, and stimulation (for full item prompts, see Table S2 in the Supplemental Material) were chosen across numerous administered scales (Vargas et al., 2020). These included the ABCD Parent Multi-Group Ethnic Identity–Revised Survey (MEIM; Phinney & Ong, 2007), which separates into “ethnic identity search” and “affirmation, belonging, and commitment” (Phinney & Ong, 2007). The ABCD Parent Vancouver Index of Acculturation (VIA)—Short Survey (Ryder et al., 2000), which subdivides into “heritage” and “American” subscores (Ryder et al., 2000). Other scales administered included the ABCD Parent Neighborhood Safety/Crime (NSC) survey modified from PhenX (Echeverria et al., 2004; Mujahid et al., 2007), ABCD Parent Acculturation (ACC) survey modified from PhenX, ABCD Youth ACC (Alegria et al., 2004; Marin et al., 1987), and the ABCD Parent Demographics survey. Primary guardians or parents of the youths completed the ABCD Parent MEIM, ABCD Parent NSC, ABCD Parent VIA, and ACC. Items related to English proficiency were omitted because of the lack of theoretical relevance to the current study. Measures were developed by the ABCD team to index environmental and cultural factors that could be relevant to development (Alegria et al., 2004; Zucker et al., 2018). Thus, these measures index exposures that occur at the systems level that have been shown to increase vulnerability to experiencing chronic stress (Vargas et al., 2020).
Objective neighborhood features
Residential history was collected through addresses where participants had lived across their lifetime. Addresses were used to determine census tracts corresponding to each location. Each tract represents census-delineated neighborhoods. Census and Federal Bureau of Investigation data were used to calculate neighborhood population density, total crimes occurring in a certain neighborhood, and the area deprivation index (ADI). The ADI metric has been successfully adapted to measure neighborhood deprivation; it is calculated using the American Community Survey 2015 5-year summary (Kind et al., 2014). Because these metrics are compiled using government data, they will be referred to as objective neighborhood features, drawing a contrast from neighborhood features of interest that are also assessed through self-report (e.g., the NSC).
Theorized systemic environmental exposure domains
As mentioned above, in the present investigation, we were not interested in exposure to individual-level stressors, which have traditionally received more exposure in the clinical literature (e.g., bullying, family conflict, and other individual-level stress exposures). Rather, the study sought to home in on environmental risk exposure occurring at the systems level, building on a broader literature of various systemic environmental exposures and their relations to psychotic disorder incidence. Thus, the variables that were chosen for inclusion reflect only exposure to systems-level factors. The ABCD data set provided us with a valuable opportunity to pull as many relevant variables that were theoretically consistent with the three previously hypothesized domains of systemic exposures: discrepancy, deprivation, and stimulation (Vargas et al., 2020). For discrepancy, the MEIM, VIA, and ACC scales were used, consistent with evidence that a lack of sense of belonging within one’s culture and lack of participation and engagement with the majority culture and with one’s culture are cultural or systems-level factors that can confer a lack of social capital and social exclusion (Emerson et al., 2018; Veling et al., 2008; Yang et al., 2018). For the deprivation domain, the ABCD parent demographic survey was used to index lack of access to environmental enrichment (through questions indexing access to resources such as access to doctors if needed, food, and utilities). For objective measures, the ADI was used as a measure of neighborhood deprivation. Finally, for the stimulation domain, high-crime regions and urban areas with high population density have been theorized to comprise high attentional demands by engaging threat neural correlates and conferring higher arousal of stress systems (Freeman et al., 2015; Gong et al., 2016; Newbury et al., 2017). Thus, the NSC survey was chosen, which assesses neighborhood safety. For objective measures, total neighborhood crimes and population density were chosen as part of the stimulation domain.
EFA
To determine whether environmental risk factors would fall within hypothesized domains (Vargas et al., 2020), we conducted an EFA 1 on the self-report scales using the minimum residuals method (Comrey, 1962) with the psych package (Version 2.0.12; Revelle, 2017) for the R software environment (Version 3.6.3; R Core Team, 2020). An EFA was not conducted for the objective measures because they spanned two theorized domains (stimulation and deprivation) with only three items or objective measures, likely not comprising enough items for a stable multifactor EFA (Raubenheimer, 2004). Given the theoretical expectation that some factors would correlate, an oblimin rotation was chosen. The number of factors was decided according to inspection of the scree plot as well as theoretical consistency and interpretability. A cutoff value of 0.4 was chosen for factor loadings; items falling beneath this threshold were excluded (Peterson, 2000). The total sample of 6,415 was randomly split in half to create two samples, one for the EFA (n = 3,207) and the other for the CFA (n = 3,208; for factor correlations and correlation matrix, see Tables S3 and S4 in the Supplemental Material).
CFA
The solution found in the first sample using EFA was tested in the second sample using CFA with R packages psych (Revelle, 2017), lavaan (Version 0.6-7; Rosseel, 2012), and semTools (Version 0.5-0; Jorgensen et al., 2018). In the case that the model did not achieve adequate fit according to conventional thresholds, modification indices were used to make theoretically consistent modifications to improve the final model fit (Brown & Moore, 2012). Conventional thresholds include model χ2 p value > .05, comparative fit index (CFI) ≥ .90, Tucker-Lewis index (TLI) ≥ .90, root mean square error of approximation (RMSEA) < 0.08, and standardized root mean square residual (SRMR) < 0.08 (Hooper et al., 2008). Because of the χ2 test’s sensitivity to large sample sizes, this index was deemphasized when assessing model fit given the size of our sample (Hu & Bentler, 1999).
Associations between self-report factors, objective neighborhood features, and PLEs
To determine whether PLEs and self-report factors (from CFA solution) related to each other, we conducted nonparametric Spearman’s correlations, adjusting for age and sex. Spearman’s correlations adjusting for age and sex were also conducted to test the association between PLEs and objective neighborhood features. A central aim was to compare the strength of associations between PLEs and self-report items and between objective neighborhood features and PLEs. Differences between correlations between PLEs, self-report, and objective neighborhood features indexing the same construct were also tested (i.e., between self-report stimulation and neighborhood safety and between stimulation and total number of crimes in neighborhood, neighborhood deprivation index, and between self-reported deprivation and lack of resources). To test whether associations were significantly different, we converted correlation coefficients into a z score using Fisher’s r to z transformation (Meng et al., 1992). Then, the asymptotic covariance of the estimates was computed and used in an asymptotic z test to determine whether one correlation was significantly greater than the other (Lee & Preacher, 2013; Steiger, 1980). All analyses were Bonferroni corrected by dividing α = .05 by the number of tests conducted (Bonferroni, 1935; Shaffer, 1995).
Results
EFA for self-report items
The Kaiser-Meyer-Olkin (KMO) and Bartlett’s test of sphericity indicated the data were appropriate to analyze. The KMO measure of sampling adequacy was 0.89 (Kaiser, 1970, 1974; Kaiser & Rice, 1974). Bartlett’s test of sphericity was significant, χ2(435) = 59,580.6, p < .05 (Bartlett, 1951). The five-factor solution was chosen (a) after inspection of the scree plot, (b) given previous theoretical support, and (c) considering that six-and seven-factor solutions were difficult to interpret and had an insufficient number of primary loadings. One item was eliminated from the original variable set because of failure to meet the minimum criteria of having a primary factor loading of 0.4 or above (“In the past 12 months, were evicted from your home for not paying the rent or mortgage”). After item removal, the five-factor EFA solution using minimum residuals with an oblimin rotation explained 55% of the variance. All items in the analysis had primary loadings greater than 0.4 (Table 1). Reliability across items was adequate, total ω = .94 (Revelle & Zinbarg, 2009).
Exploratory Factor Analysis Loadings, Uniqueness, and Communality Metrics
Note: h2 = communalities metrics; u2 = uniqueness metrics; com = the complexity of component loadings per item.
Factor loadings were partially consistent with theoretical predictions. Factor 5 was theoretically consistent with the deprivation domain. Factor 3, in turn, was theoretically consistent with the stimulation domain. The discrepancy domain comprised the remaining three factors. Factors 1 and 4 index participation in heritage and American culture, respectively. Factor 2, in turn, indexes an individual’s sense of belonging with ethnic group. Factors were named according to corresponding items and theoretical hypotheses. Factor 1 was named Heritage Culture Participation, Factor 2 was deemed Sense of Belonging With Ethnic Group, Factor 3 was named Neighborhood Safety, Factor 4 was deemed American Culture Participation, and Factor 5 was deemed Deprivation (see Fig. 1, Table 1).

Exploratory factor analysis of self-report items. Representation of the distribution and interquartile range for endorsement of items for each respective factor appears at left. Neighborhood Safety is part of the stimulation domain. Sense of Belonging With Ethnic Group and American Culture Participation constitute the discrepancy domain. Deprivation constitutes the deprivation domain. Heritage Culture Participation did not relate to psychotic-like experiences (PLEs) and was not included within the three stimulation, discrepancy, and deprivation (SDD) domains. Factor correlations with medium effect sizes were included. For the rest of the factor correlations, see Table S4 in the Supplemental Material available online.
CFA and post hoc modification indices
The initial CFA of the five-factor model found using EFA did not reach adequate fit thresholds, robust χ2 = 7,927.997, p < .0001, robust CFI = .820, robust TLI = .802, robust RMSEA = 0.087, robust SRMR = 0.056. Modification indices were run to identify contributors to the inadequate fit. These included likely cross-loadings between the American Culture Participation and Heritage Culture Participation factors, including in items that indexed participation in heritage versus American or mainstream cultural traditions, comfortability interacting with individuals of same heritage culture versus typical American people, and interest in having friends from heritage versus typical American friends. It is also likely that there was some redundancy within the Sense of Belonging With Ethnic Group factor given that residuals correlated across several items, including items that assessed strong sense of belonging and strong attachment to ethnic group, doing things to understand ethnic background and talking to others to learn more about ethnic group, and participating in heritage versus American cultural traditions. Within the American Culture Participation factor, items related to maintaining or developing American mainstream cultural practices, believing in mainstream American values, and comfortability interacting with typical American people and enjoying typical American entertainment showed correlations. Finally, within the Deprivation factor, items related to limited doctor availability versus limited dentist availability showed correlated residuals. After modifications allowing residual correlations for these items, the model achieved adequate fit thresholds (except for the χ2 test, which was expected given the sample size): robust χ2 = 4,520.459, p < .0001, robust CFI = .909, robust TLI = .897, robust RMSEA = 0.06, robust SRMR = 0.049 (for individual item loadings, see Table S5 in the Supplemental Material).
Associations between self-report factors and PLEs
Greater endorsement of the Deprivation factor or deprivation domain related to more PLEs (r = .10; Table 2, Fig. 2a). Likewise, endorsement of safer neighborhoods or the stimulation domain related to fewer PLEs (r = −.09), consistent with predictions. Contrary to predictions, the Heritage Cultural Participation factor was not found to significantly relate to PLEs. However, the other hypothesized discrepancy domain factors (Sense of Belonging With Ethnic Group and American Culture Participation) related to PLEs such that less sense of belonging with ethnic group (r = −.05) and greater American culture participation (r = −.03) predicted fewer PLEs. The strength of the correlations was subsequently compared to gauge the relative strength of observed associations; the Deprivation factor association with PLEs was significantly greater than that of the discrepancy domain American Culture Participation and Sense of Belonging With Ethnic Group factors, which survived Bonferroni correction (see Table 2, Fig. 2a). Likewise, the association between the stimulation domain/Neighborhood Safety factor and PLEs was significantly stronger than that of the discrepancy domain/American Culture Participation and Sense of Belonging With Ethnic Group factors. However, the association between self-report deprivation and PLEs was not significantly stronger than the association between self-report stimulation domain/Neighborhood Safety and PLEs.
Spearman’s Partial Correlations for Self-Report and Objective Measures Controlling for Sex and Age and Comparisons of the Strength of Association
Note: Factor 3 (Neighborhood Safety), neighborhood population density, and crime are part of the stimulation domain. Factor 2 (Sense of Belonging With Ethnic Group), Factor 1 (Heritage Culture Participation), and Factor 4 (American Culture Participation) constitute the discrepancy domain. Factor 5 (Deprivation) and neighborhood deprivation constitute the deprivation domain. CI = confidence interval.
This column indicates whether the association survived correction for multiple comparisons. The threshold was determined by dividing the α level (.05) by the number of comparisons (five for self-report measures, three for objective measures). The self-report value is .01. The objective measure value is .01667. For self-report comparisons (six), the Bonferroni-corrected threshold was .00833. For self-report/objective measure comparisons, the Bonferroni-corrected threshold was .025. For objective measure comparisons, the Bonferroni-corrected threshold was .01667.

Associations between self-report factors discrepancy domain/Sense of Belonging With Ethnic Group, discrepancy domain/American Culture Participation, Deprivation, stimulation domain/Neighborhood Safety, and Heritage Culture Participation. (a) Self-report factors and PLEs. (b) Objective neighborhood features and PLEs.
Associations between objective neighborhood features and PLEs
Regarding objective neighborhood measures, increased stimulation domain/neighborhood population density (r = .07) related to increased experience of PLEs. This was not the case for stimulation domain/neighborhood total crimes, which did not observe a significant association. Deprivation domain/neighborhood deprivation (r = .14) related to increased experience of PLEs. The strength of the correlations was compared; the association between deprivation domain/neighborhood deprivation and PLEs was significantly stronger than the associations observed for both stimulation domain/total crimes and stimulation domain/population density (Table 2, Fig. 2b).
Exploratory comparison of objective and self-report measures
Finally, as an exploratory aim, self-report measures and their corresponding objective measures were compared in terms of strength of the association observed. The association between deprivation domain/neighborhood deprivation and PLEs was significantly stronger than that between self-reported Deprivation factor and PLEs. However, the association between self-reported stimulation domain/Neighborhood Safety and PLEs was significantly stronger than that of stimulation domain/neighborhood total crimes, for which we did not observe a significant association with PLEs.
Discussion
In the current study, we used a large representative sample to examine whether structural environmental exposures could be distinguished as theorized in the SDD model of psychosis. Furthermore, exposures (self-report and objective) were explored in relation to PLEs. Although we expected three factors to emerge, the EFA of self-report data identified five factors, corresponding to Heritage Culture Participation, Sense of Belonging With Ethnic Group, American Culture Participation, Neighborhood Safety, and Deprivation. Although more factors were fit than anticipated and fit indices were suboptimal before modification indices were applied, these factors were partially consistent with the SDD model. Critically, the present investigation is among the first to compare the relative strength of the association between these distinct domains of environmental exposure to stress and psychosis risk. Environmental exposure to stress across the three domains related to PLEs (|rs| range = .03–.14): The Deprivation factor’s association with PLEs was significantly greater than that of the discrepancy domain/American Culture Participation and discrepancy domain/Sense of Belonging With Ethnic Group factors, suggesting deprivation exposures could relate more strongly to psychosis vulnerability. Consistent with this interpretation, in terms of objective neighborhood measures, the association between neighborhood deprivation and PLEs was also significantly stronger than that between stimulation domain/neighborhood population density and stimulation domain/neighborhood total crime and PLEs. Self-report stimulation exposures (Neighborhood Safety) also showed a stronger association compared with discrepancy domain/American Culture Participation and discrepancy domain/Sense of Belonging With Ethnic Group. Taken together, results aid in refining and building on the SDD theory of psychosis, which could in time be informative to relevant public policy conceptualizations of prevention and intervention.
EFA was used for self-report scales measuring systemic factors applicable to the SDD theory. Using self-report items completed by the parents alleviated concerns common to self-report scales regarding state effects (e.g., mood and fatigue) influencing ratings of exposure (the inclusion of objective measures was also helpful in this regard). As expected, items taken from the NSC survey loaded onto a Neighborhood Safety factor (Echeverria et al., 2004). These items reflected feelings of safety in one’s living environment, consistent with the stimulation SDD domain. Items inquiring about income or resource availability, indexing degrees of environmental enrichment, consistent with the deprivation domain, loaded onto a Deprivation factor. Finally, items related to culture separated into factors relating to Heritage Culture Participation, American Culture Participation, and Sense of Belonging With Ethnic Group. Of these factors, Sense of Belonging With Ethnic Group and American Culture Participation are conceptually consistent with the discrepancy domain because they index current feelings of belonging and participation. The fact that these items did separate into three factors could indicate that the discrepancy domain requires further granularity in conceptualization—perhaps separate domains or distinct subdomains more closely represent these exposures. Future investigations will benefit from these insights in further refining and modifying the SDD theory.
It is necessary to highlight that CFA did not show adequate model fit for the five-factor solution supported by the EFA in Sample 1. As a result, any interpretations relating results to the SDD theory should be considered to be preliminary. Future investigations will be essential in further refining the model and identifying ideal self-report measures for each of the domains. Modification indices suggest that cross-loadings between the American Culture Participation and Heritage Culture Participation factors contributed to the degree of model fit (possibly at least partially because of similar wording across items). Furthermore, items for the American Culture Participation and Heritage Culture Participation factors did not explicitly assess for relative degrees of participation or comfortability; the items did not directly assess to what degree someone participates in mainstream American culture relative to heritage culture, which could have resulted in correlated residuals between items.
Modification indices also revealed the presence of correlated residuals among items within the Sense of Belonging With Ethnic Group factor and for items within the Deprivation factor, possibly because of redundancy across items within each factor. Although the ABCD data set provided an invaluable opportunity to test the SDD theory in a preliminary fashion, it also included items that were not originally designed to test the SDD theory. Thus, some “noise” is likely attributable to the fact that the items do not fully capture the theoretical domains. Perhaps the domains represent formative, rather than reflective, latent variables. Future refining of current measures to improve granularity and theoretical consistency to the SDD theory will be helpful, as will further refining and revising of the SDD theory itself. Nonetheless, given the limited existing literature on qualitatively distinct systemic stressors and their relation to psychosis vulnerability, the current investigation offers a preliminary starting point to understanding these exposures through self-report items.
Furthermore, the SDD theory concerns itself with stressors occurring at the systems level but also recognizes the multitude of individual-level factors that could also contribute to indicators—thus, cross-loadings would be expected. Although cross-loadings could likely be present and the items were not uniquely designed for measuring the intended latent constructs, the observed factors could nonetheless be useful in conceptualizing systemic stressors and predicting psychosis vulnerability. For example, Big 5 personality models, which have often shown cross-loadings and suboptimal fit with a five-factor CFA structure, are nonetheless highly useful models and reliably predict outcomes of interest (Gurven et al., 2013; Marsh et al., 2010; McCrae & Costa, 1997). Thus, despite suboptimal CFA fit, in the current study, we proceeded to correlate the self-report factors to PLEs, as we did with objective measures of the domains (neighborhood crime, population density, and deprivation), to increase the understanding of whether the factors are indeed useful in predicting psychosis vulnerability. Associations with PLEs should be treated as preliminary and interpreted with caution.
Collectively, findings showed that both objective (deprivation domain/neighborhood deprivation and stimulation domain/neighborhood population density) and subjective endorsement of stimulation and deprivation exposures contribute to the association with PLEs. Decreased reports of stimulation/Neighborhood Safety (along with higher stimulation/neighborhood population density) predicted greater PLE endorsement. Results are consistent with previous investigations on adolescents and adults that found that stimulation exposures relate to increased risk for developing a psychotic disorder (Bhavsar et al., 2014; Freeman et al., 2015; Gong et al., 2016; Kirkbride et al., 2012; Newbury et al., 2018; Wilson-Genderson & Pruchno, 2013). Increased reporting of deprivation and higher objective neighborhood deprivation predicted greater endorsement of PLEs. Similar to stimulation exposures, results are consistent with a strong body of literature suggesting that exposure to both individual and neighborhood deprivation can confer risk for psychosis (Bhavsar et al., 2014, 2018; Kirkbride et al., 2012; Lasalvia et al., 2014; O’Donoghue et al., 2015; Omer et al., 2014) and adversely affect an individual’s physical health and functional outcomes (Akman et al., 2004; Beckett et al., 2006; Gee et al., 2013; Kobayashi et al., 2017; Lang et al., 2009; Mackes et al., 2020; McCann et al., 2018; McLaughlin, Sheridan, Winter, et al., 2014; Mensah & Hobcraft, 2008; Richards et al., 2015; Uysal et al., 2005; Wiesel & Hubel, 1965). Furthermore, results suggest that the relation between these environmental exposures and psychosis spectrum symptoms extends to nonclinical psychosis and is evident as early as late childhood to late adolescence. Note that self-report factors for stimulation/Neighborhood Safety and Deprivation were correlated (r = −.3)—which was not the case for self-report factors relating to the discrepancy domain. This is consistent with our theory that the systemic exposures share common underlying mechanisms as well as phenomenologically distinct effects and manifestations.
Of the self-report factors relevant to the discrepancy domain, decreased Sense of Belonging With Ethnic Group and increased American Culture Participation related to less endorsement of PLEs. Observed associations are congruent with evidence that high ethnic density (Schofield et al., 2017; Termorshuizen et al., 2014; Veling et al., 2008) and social cohesion (Crush et al., 2018) can serve as protective factors. Results extend the existing literature by providing evidence that discrepancy systemic exposures could relate to PLEs as early as late childhood to early adolescence (Allardyce et al., 2005; Crush et al., 2018; Lasalvia et al., 2014; Schofield et al., 2017; Silver et al., 2002; Termorshuizen et al., 2014; Van Os et al., 2000; Veling et al., 2008). The Heritage Culture Participation factor, on the other hand, was not associated with PLEs. The lack of significant association between Heritage Culture Participation and PLEs could be due to insufficient factor specificity regarding feelings of belonging. That is, one could participate in one’s heritage culture and yet not feel a sense of belonging with one’s surroundings more broadly or with the majority culture. Psychosis environmental risk factors theorized by the SDD theory discrepancy domain include ethnic minority status (Lasalvia et al., 2014; Termorshuizen et al., 2014), low ethnic density (Schofield et al., 2017; Veling et al., 2008), and social fragmentation (Allardyce et al., 2005; Silver et al., 2002; Van Os et al., 2000). Perhaps the Heritage Culture Participation factor does not fully capture these experiences of social exclusion or lack of belonging. American Culture Participation, on the other hand, did relate to PLEs. The association makes sense because American culture in the United States would be the “majority” culture.
Given the vast evidence of low ethnic density and minority status conferring psychosis risk (Lasalvia et al., 2014; Schofield et al., 2017; Termorshuizen et al., 2014; Veling et al., 2008), perhaps American Culture Participation indexes comfortability within the majority culture, which could directly affect overall social capital and sense of social cohesion, effecting psychosis risk through this mechanism (Butler & Muir, 2017; Crush et al., 2018; Schellenberg et al., 2018; Verhaeghe & Tampubolon, 2012). The ABCD data set provided an opportunity to test some systemic exposures that fall within the discrepancy domain. Future investigations assessing other environmental exposures theorized to fall under the discrepancy domain and further establishing specificity could aid in offering more nuance to current results.
The literature thus far has been limited in comparing the magnitude of associations between distinct environmental exposures and psychosis vulnerability, and the proposed SDD model originally was more or less agnostic regarding the relative magnitude of the associations between differing systemic exposures. The three domains showed differential associations with PLEs, which demonstrates that their relative contribution could be informative to consider and highlights the importance of considering individual domains (as opposed to an aggregate of systemic environmental stress exposure). In the current sample, observed associations between deprivation (objective and self-report) exposures and PLEs were significantly stronger than associations between PLEs and discrepancy exposures; effect sizes for deprivation were twice as large (deprivation: rs range = .09–.14; discrepancy: rs range = .03–.05). The observed difference is consistent with animal and human literature that suggests that lack of neurodevelopmentally appropriate enrichment can have widespread consequences, affecting a host of key systems necessary for general functioning as well as overall health and well-being (Akman et al., 2004; Beckett et al., 2006; Gee et al., 2013; Kobayashi et al., 2017; Lang et al., 2009; Mackes et al., 2020; McCann et al., 2018; McLaughlin, Sheridan, Winter, et al., 2014; Mensah & Hobcraft, 2008; Richards et al., 2015; Uysal et al., 2005; Wiesel & Hubel, 1965). Note that the deprivation exposure objective measures also exhibited significantly stronger associations with PLEs (r = .14) compared with objective stimulation/population density (r = .07) and stimulation/neighborhood crimes (r = .02). Future investigations will be necessary to clarify possible mechanisms through which systemic exposures could differ in magnitude.
PLEs are complex and include a multitude of putative contributing factors. Likewise, systemic environmental exposures are complex and multifaceted and include a rich variety of complex protective and exacerbating factors (including potent individual-level stressors) to moderate and mediate relationships. This can lead to measures of systemic effects and neighborhood effects tending toward what could be considered small effect sizes (Barrington et al., 2014; Crump et al., 2011; Cubbin & Winkleby, 2005; Cummins et al., 2005; Forsberg et al., 2018; Gale et al., 2011; Jaffe et al., 2005; Kirkbride et al., 2012; Lang et al., 2009). However, it is critical to interpret effect sizes in the context of the relationships they are depicting (Funder & Ozer, 2019). Effect sizes could be small in the face of single events and prove more ultimately consequential as effects accumulate over the medium and long terms (Funder & Ozer, 2019). Thus, effects observed during childhood could accumulate over many years and have a nontrivial impact across the lifetime. Thus, observed effects, although small, can be meaningful when considered in context, especially when assessing aggregate effects for communities, regions, or even countries as a whole; indeed, a host of rather small effects has served as springboards for effective health policy initiatives (Arnett et al., 2019; Funder & Ozer, 2019; Schwingshackl et al., 2015).
An exploratory aim sought to compare the magnitude of PLE associations between self-report and objective environmental exposures. The PLE association with objective neighborhood deprivation was significantly stronger in magnitude than the PLE association with self-report deprivation. However, this was not the case with stimulation exposures. The PLE association with self-report stimulation/Neighborhood Safety was significantly stronger in magnitude compared with objective stimulation/neighborhood total crimes (which did not show a significant association with PLEs). Although preliminary, results suggest that perhaps the association between PLEs and stimulation/crime exposure is at least partially contingent on conscious awareness: In this case, self-reported lack of safety could represent heightened awareness of the objective circumstances, or it could represent inaccurate reporting due to heightened vigilance. Alternatively, it could be the case that collection sites had different average crime rates, which could have blunted an existing effect when the variance related to site was partialled out in analyses. Future investigations are needed to parse out these possibilities.
Current results constitute a promising, although preliminary, start to broader questions with the potential to inform models of psychosis vulnerability and highlight the value of targeting specific environmental factors in preventive health policy efforts for psychotic disorders. Results ought to be treated as preliminary and interpreted with caution. Initial findings aid in refining and modifying the SDD theory moving forward and must be contextualized within certain limitations. There is the problem of intracategory variability that characterizes many, if not most, stressful life event measures (Dohrenwend, 2006). The self-report measures used in this study showed robust reliability (Alegria et al., 2004; Echeverria et al., 2004; Marin et al., 1987; Phinney & Ong, 2007; Ryder et al., 2000). Yet future studies may benefit from further building on practices recommended for limiting intracategory variability, such as implementing more closed probes as well as more stringent inclusion and exclusion criteria for considering an experience endorsed (Dohrenwend, 2006). Thankfully, our concerns related to self-report and intracategory variability were tempered by the inclusion of objective measures of the constructs of interest. Future investigations will benefit from continuing to benefit from self-report data and complement it with corroborating evidence or objective measures. Likewise, future investigations may benefit from increasing granularity and more fully incorporating other facets of systems-level environmental exposures and integrating them to build on and improve conceptual frameworks of systems-level environmental stressors.
The present work focused on systems-level environmental exposures, aiming to distinguish distinct components of systemic exposures that could confer stress. The theoretical interest was identifying qualitatively distinct systemic exposures with theorized intermediary mechanisms (based on prior research). We chose not to include race in our models given the overwhelming evidence of systemic disadvantage that disproportionally affects Black, indigenous, and people of color (BIPOC) and ethnic minorities (Adler & Stewart, 2009; Braveman, 2014; Braveman et al., 2011; Brondolo et al., 2009; Jackson et al., 2010; Joynt et al., 2011). If the effect of race is partialled out, one may miss meaningful information, given that ethnic minorities and BIPOC are disproportionately more likely to experience systemic environmental risk factors. Thus, the aim of the present study was to detail differing components of systemic exposures such that these may be better understood and conceptualized as a whole when considering systemic disadvantage. Although race was beyond the scope of the study questions, future investigations could further investigate interactions of types of systemic exposures with race and other factors and build on this knowledge. Beyond this point, it is necessary to consider that the current investigation used one (albeit well-powered) sample and thus did not replicate the established factor structure in an independent sample. Therefore, it will be critical for future work to confirm the observed factor structure in independent samples, thus establishing replicability and increasing generalizability of the work. The current results ought to be interpreted as preliminary until replicability is established.
Future work would also benefit from exploring neural and biological mechanisms that could underlie the associations observed between environmental exposures and psychosis vulnerability. Collecting information on the precise timing of exposure would further add richness to neurodevelopmental conceptualizations of risk and resilience. Examining interactions between individual-level factors (e.g., trauma, family environment, exposure to bullying), systemic factors, and relations to psychosis vulnerability in children and adolescents would also be a worthwhile line of inquiry. Finally, longitudinal investigations would aid our ability to predict the directionality of the associations and account for confounds such as social drift.
Supplemental Material
sj-pdf-1-cpx-10.1177_21677026211016415 – Supplemental material for Differentiating Kinds of Systemic Stressors With Relation to Psychotic-Like Experiences in Late Childhood and Early Adolescence: The Stimulation, Discrepancy, and Deprivation Model of Psychosis
Supplemental material, sj-pdf-1-cpx-10.1177_21677026211016415 for Differentiating Kinds of Systemic Stressors With Relation to Psychotic-Like Experiences in Late Childhood and Early Adolescence: The Stimulation, Discrepancy, and Deprivation Model of Psychosis by Teresa Vargas, Katherine S. F. Damme, K. Juston Osborne and Vijay A. Mittal in Clinical Psychological Science
Footnotes
Acknowledgements
Data used in the preparation of this article were obtained from the Adolescent Brain Cognitive Development (ABCD) Study (https://abcdstudy.org), held in the National Institute of Mental Health (NIMH) Data Archive. This is a multisite, longitudinal study designed to recruit more than 10,000 children ages 9 to 10 and follow them over 10 years into early adulthood. A full list of supporters is available at https://abcdstudy.org/federal-partners.html. A listing of participating sites and a complete listing of the study investigators can be found at https://abcdstudy.org/Consortium_Members.pdf. ABCD consortium investigators designed and implemented the study and/or provided data but did not necessarily participate in analysis or writing of this article. This article reflects the views of the authors and may not reflect the opinions or views of the National Institutes of Health or ABCD consortium investigators. The ABCD data repository grows and changes over time. The ABCD data used in this report can be found at
.
Transparency
Action Editor: Andrew Littlefield
Editor: Kenneth J. Sher
Author Contributions
T. Vargas and V. A. Mittal conceptualized the study. T. Vargas did data analyses and wrote the manuscript draft with supervision from V. A. Mittal. K. S. F. Damme and K. J. Osborne provided analysis suggestions as well as conceptual feedback. All of the authors approved the final manuscript for submission.
Notes
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
