Shared or common environment (C), along with genetic effects and individual-specific environment, has long been a key parameter in the quantitative genetic studies of individual differences in both model organisms and humans (Falconer, Reference Falconer1989; Lynch & Walsh, Reference Lynch and Walsh1998; Mather & Jinks, Reference Mather and Jinks1982). Years ago, traditional quantitative methods proposed a decomposition of genetic effects into additive, dominance, and epistatic components (Falconer, Reference Falconer1989; Mather & Jinks, Reference Mather and Jinks1982). With modern molecular methods, it is now possible to further sub-divide these genetic effects into the proportion that reflects common single-nucleotide variants, rare variants, and various genomic abnormalities (Lee et al., Reference Lee, DeCandia, Ripke, Yang, Sullivan, Goddard and Wray2012; Purcell et al., Reference Purcell, Wray, Stone, Visscher, O'Donovan, Sullivan and Scolnick2009, Reference Purcell, Moran, Fromer, Ruderfer, Solovieff, Roussos and Sklar2014). By contrast, potential sub-components of the shared environment have received less attention. The only subcomponents that are occasionally used in modeling are the special twin environment (T) and the proportion of variance in shared environment resulting from parental behaviors (W; Truett et al., Reference Truett, Eaves, Walters, Heath, Hewitt, Meyer and Kendler1994).
There is an obvious reason for these developments. Examining only reared together twin pairs — by far the most popular design in human behavior genetics — permits only the estimation of a single parameter for the shared environment. Including full siblings permits the addition of the special twin environment. Studying parents in twin-family designs can allow the estimation of a shared environment resulting from parental behaviors. So, estimation of sub-components of the shared environment depends critically on the availability of informative samples.
In a thoughtful review paper, Burt effectively countered the claims of the influential article by Plomin and Daniels (Reference Plomin and Daniels1987) that shared environmental effects are of limited import for psychological, psychiatric, and behavioral traits. She conducted a meta-analysis of twin and adoption studies of internalizing and externalizing psychopathology through adolescence, concluding that c2 estimates for most traits were in the range of 10–20%. She also noted that shared environmental influences are typically stable, unlikely to be artefactual, and may be particularly important in high-risk environments. In a subsequent review, Burt (Reference Burt2014) noted the limitations of the traditional twin design for the study of the shared environment and makes the important point that ‘shared environmental influences are not restricted to familial influences’ (p. 307). Community and school effects, she observed, both directly and as proxies for other influences like peer deviance, likely contribute to c2 estimates. Finally, in a recent empirical study, Burt et al. (Reference Burt, Klump, Kashy, Gorman-Smith and Neiderhiser2015) showed, using community samples, a substantial impact of neighborhoods on the key externalizing dimension of rule breaking.
In this report, congruent with the research perspective articulated by Burt, we attempt to sub-divide the shared environment into four subcomponents termed (1) close family, (2) family, (3) household, and (4) community. This has been made possible by using Swedish population registries that contain the informative relationships necessary to estimate these parameters.
We examine three externalizing syndromes that we have previously examined in Sweden using both adoption and traditional twin or twin-sibling designs: drug abuse (DA; Kendler, Maes et al., Reference Kendler, Maes, Sundquist, Ohlsson and Sundquist2013; Kendler et al., Reference Kendler, Sundquist, Ohlsson, Palmer, Maes, Winkleby and Sundquist2012), alcohol use disorder (AUD; Kendler, Ji et al., Reference Kendler, Ji, Edwards, Ohlsson, Sundquist and Sundquist2015; Kendler, Lonn et al., Reference Kendler, Lonn, Maes, Lichtenstein, Sundquist and Sundquist2016) and criminal behavior (CB; Kendler, Larsson Lönn et al., Reference Kendler, Larsson Lönn, Morris, Sundquist, Långström and Sundquist2014; Kendler, Maes et al., Reference Kendler, Maes, Lönn, Morris, Lichtenstein, Sundquist and Sundquist2015). We study these disorders because, unlike more traditional psychiatric disorders where shared environmental effects are difficult to detect, we have, for all of them, found both evidence for substantial shared environmental effects in standard twin or twin-siblings studies (Kendler, Lonn et al., Reference Kendler, Lonn, Maes, Lichtenstein, Sundquist and Sundquist2016; Kendler, Maes et al., Reference Kendler, Maes, Lönn, Morris, Lichtenstein, Sundquist and Sundquist2015; Kendler, Maes et al., Reference Kendler, Maes, Sundquist, Ohlsson and Sundquist2013) and for an impact on risk in adoptees of adoptive family characteristics (Kendler, Ji et al., Reference Kendler, Ji, Edwards, Ohlsson, Sundquist and Sundquist2015; Kendler, Larsson Lönn et al., Reference Kendler, Larsson Lönn, Morris, Sundquist, Långström and Sundquist2014; Kendler et al., Reference Kendler, Sundquist, Ohlsson, Palmer, Maes, Winkleby and Sundquist2012). Evidence for shared environmental effects for these syndromes has also emerged from other twin samples for DA (Kendler et al., Reference Kendler, Jacobson, Prescott and Neale2003) and AUD (Kendler et al., Reference Kendler, Prescott, Neale and Pedersen1997; McGue et al., Reference McGue, Pickens and Svikis1992), and from meta-analyses of twin and adoption samples for externalizing/antisocial behavior (Rhee & Waldman, Reference Rhee and Waldman2002) and AUD (Verhulst et al., Reference Verhulst, Neale and Kendler2015). Furthermore, with a common pathway model fitted to monozygotic (MZ) and dizygotic (DZ) twin pairs, we have found evidence of both shared and unique genetic and common environmental effects (Kendler, Lonn et al., Reference Kendler, Lonn, Maes, Lichtenstein, Sundquist and Sundquist2016) suggesting that we might gain further insight into the operation of shared environmental effects by undertaking a further multivariate analysis of these three syndromes. Therefore, in this study, we also apply a common factor model separately in males and females with the goal of understanding the sources of the shared environmental influences that are common to, versus unique to, DA, AUD, and CB, and the degree to which the pattern of results are similar or different in men and women.
Methods
As outlined previously for DA, AUD, and CB, we used linked data from multiple Swedish nationwide registries and healthcare data (Kendler, Ji et al., Reference Kendler, Ji, Edwards, Ohlsson, Sundquist and Sundquist2015; Kendler, Larsson Lönn et al., Reference Kendler, Larsson Lönn, Morris, Sundquist, Långström and Sundquist2014; Kendler et al., Reference Kendler, Sundquist, Ohlsson, Palmer, Maes, Winkleby and Sundquist2012). Linking was achieved via the unique individual 10-digit personal identification number assigned at birth or immigration to all Swedish residents. To preserve confidentiality, this identification number was replaced by a serial number. We secured ethical approval for this study from the Regional Ethical Review Board of Lund University. For a definition of DA, AUD, and CB see Appendix A.
Among all individuals born in Sweden between 1960 and 1990 (N = 3,257,987) and alive at age 15, we created seven datasets based on their genetic and geographical/household relationship. All analyses were stratified by gender, suggesting that we created 7×2 = 14 datasets. See Table 1 for a description of the datasets. AUD, CB, and DA were measured from age 15 until end of 2010 and were considered as binary variables in the models. For a definition of household and geographical status, see Appendix B.
TABLE 1 Description of Datasets and Assumptions for the Structural Equation Model
Note: *A = additive genetics; C = shared environment. **SAMS = Small area market statistics — there are approximately 9,200 SAMS throughout Sweden, their average population being around 1,000. Note that we could not model 249,00,000 or 344,000,000 pairs in OpenMx. We therefore took a random sample of 500,000 pairs to be included in the models. The correlations were almost identical in the random sample as in the original sample.
Statistical Analysis
Our structural model divides the sources of individual differences in liability to DA, AUD, and CB into the following: additive genetic (A); four different components of the shared environment: close family (C_close), family (C_family), household (C_household), community (C_community); and unique environment (E). We extended the standard models used in twin research to include six variance sources for each phenotype. The four components of shared environment (C_close, C_family, C_household, C_community) are supposed to reflect different levels of the shared environment, ranging from close family experiences (C_close) to more distant community effects (C_community). In the models, we assumed that individuals within the same twin (MZ and DZ) and full-sibling (0–2 years apart — termed full-siblings 1) pairs share all forms of the environment — from close family experiences to community effects (C_close + C_family + C_household + C_community). This assumption is based on several empirical results from our prior analyses showing that, for these externalizing phenotypes, age differences in siblings are substantially and inversely correlated with resemblance (Kendler, Maes et al., Reference Kendler, Maes, Sundquist, Ohlsson and Sundquist2013; Kendler, Morris et al., Reference Kendler, Morris, Lönn, Sundquist and Sundquist2014; Kendler et al., Reference Kendler, Sundquist, Ohlsson, Palmer, Maes, Winkleby and Sundquist2012). These results suggest that twins and close-aged siblings share more of their rearing environment than do siblings further apart in age. For full- (termed full-siblings 2) and half-siblings more distant in age (3–10 years apart), but sharing the same household for at least part of their childhood, we assume that they share the same environment as twins and full-siblings close in age, except for the close family experiences (C_family + C_household + C_community). For step-siblings, we assume that they share the same household and community experiences (C_household + C_community). This assumption is based on results from our previous studies. Using the heritability estimate (a2) for DA, CB, and AUD from our previous studies on those comparisons between pairs where the C_family component could easily be extracted (i.e., the comparison between full-siblings 2 and step-siblings as well as between half-siblings and step-siblings), we see that 10 out of 12 comparisons indicate an extra environmental similarity between pairs that are genetically related (Kendler, Ohlsson et al., Reference Kendler, Ohlsson, Edwards, Lichtenstein, Sundquist and Sundquist2016). Furthermore, we assume that non-related individuals within the same SAMS only share the community effects (C_community). For all individuals residing in the same SAMS from at least 6 years until age 15 and with a maximum of 5 years’ age difference, we created all possible pairs of non-related individuals. As it was not possible to model 249,000,000 male or 344,000,000 female pairs in OpenMx, we took a random sample of 500,000 male–male pairs and 500,000 female–female pairs to be included in the models. Unique environment (E) includes random developmental effects, environmental experiences not shared by individuals within the different types of pairs, and random error.
In the next step, we wanted to investigate to what extent genetic and environmental factors are common for the three syndromes, with a particular focus on the different components of the shared environment. In the first model, an independent pathway model, we assumed that each of the six variance components (A, C_close, C_family, C_household, C_community, and E) consists of two parts: one that is common to all three phenotypes (denoted Ac, C_closec, C_familyc, C_householdc, C_communityc, and Ec) and one that is specific to each one of them (As, C_closes, C_familys, C_households, C_communitys, and Es). In the second model, a common pathway model, we assumed that the common variance components (Ac, C_closec, C_familyc, C_householdc, C_communityc, and Ec) are mediated by a latent phenotype represented by a common factor (called ‘externalizing behavior’ in Figures 1 and 2). The common pathway model is nested within the independent pathway model. This latent phenotype mediates the common genetic and environmental effects, as the paths from Ac, C_closec, C_familyc, C_householdc, C_communityc, and Ec run via the latent phenotype.
Although our sample size is considerable, the models include many parameters and have limited statistical power. We therefore followed the recommendations based on simulations that showed that in such situations, parameter estimates from the full model are typically more accurate than those from sub-models even if the latter provides a better model fit (Sullivan & Eaves, Reference Sullivan and Eaves2002). Unfortunately, we were not able to fit 95% confidence intervals for the parameters in the common pathway model. We therefore only present the standard errors for the estimated path coefficients in males only (see limitations section for further details). Models were fit in the OpenMx software (Neale et al., Reference Neale, Hunter, Pritikin, Zahery, Brick, Kirkpatrick and Boker2015).
Results
Descriptive Findings Univariate Analyses
The prevalence of registration for DA, AUD, and CB in our entire sample of Swedish males and females are seen in Table 1. CB was most common followed by AUD and DA. The prevalence of all three syndromes were considerably higher in the step-sibling and half-sibling groups, and lower in the twin groups, likely because the twins were screened for co-operativeness in having to return questionnaires about zygosity. For females, the full-siblings had lower prevalence. The overall prevalence in the entire population was 5.0% (DA), 4.6% (AUD), and 17.7% (CB) for males, and 2.1%, 2.0%, and 6.5%, respectively, for females.
The tetrachoric correlations within pairs for DA, AUD, and CB are seen in Figure 1a (males) and Figure 1b (females). Focusing on males and the different C components, the following findings are noteworthy. First, the correlation among unrelated individuals within the same SAMS is above zero for all phenotypes, suggesting that parts of the shared environment could be attributed to community experiences. Second, the correlation among unrelated individuals within the same household (step-siblings) is higher than the correlation among unrelated individuals within the same SAMS, suggesting that parts of the shared environment could be attributed to household experiences. Third, for AUD and CB, the correlation between full-sibling pairs close in age is higher than the correlation between full-sibling pairs further apart in age, giving evidence for a C_close component. Also, the DZ and full-sibling (close in age) correlation is more than half the MZ correlation, suggesting a total C effect (C_close + C_family + C_household + C_community) different from zero.
FIGURE 1a Tetrachoric correlations within the different types of pairs (males).
FIGURE 1b Tetrachoric correlations within the different types of pairs (females).
Focusing on females and the different C components, it is noteworthy that the correlation among unrelated individuals within the same SAMS is above zero for all phenotypes, but substantially lower than among males. Also, the within-household correlation (step-siblings) is low, and for CB and AUD actually lower than among individuals in the same SAMS, suggesting no appreciable household effect. On the other hand, the correlation between full-sibling pairs close in age is higher than the correlation between full-sibling pairs further apart in age, giving evidence for a C_close component. Also, except for AUD, the DZ and full-sibling (close in age) correlation is more than half the MZ correlation, suggesting a total C effect (C_close + C_family + C_household + C_community) different from zero.
Model Fitting: Univariate Models
In Table 2, the six variance components (a2+c_close2+c family 2+c household2+c community2+e2) for the three syndromes are shown for males. Focusing on shared environmental effects for DA, the C_community and C_household were most important, capturing more than 85% of the total effect of the shared environment. For CB and AUD, the patterns were different; the community and household captured less than 50% of the total effect of the shared environment, while the close family effect was more important. Note, however, that the 95% confidence intervals for the variance components are broad (except for the community effects). For females (Table 2), the pattern differed. For all syndromes, the C_close component was most important. For AUD and CB, all the other C components were small. For DA, the C component increased as the nearness to the individual increased.
TABLE 2 Variance Components from the Univariate Models
Model Fitting: Multivariate Models
We first fitted an independent pathway model, which resulted in a log-likelihood and AIC of 3110534 and -8589022, respectively. The common pathway model, which included seven fewer parameters, had the following log likelihood and AIC: Likelihood: 3110392; estimated parameters 44; AIC: -8589076. The fit of the common pathway model was superior, and we therefore focused on the results from this model.
The parameter estimates for the four different C components as well as their loadings on the common latent factor for the common pathway model are seen in Figure 2a (males) and 2b (females). The estimated total contributions of additive genetic, shared (four parts) and unique environmental effects on the three forms of externalizing behaviors, and the proportion of each that is common versus specific in its effect are seen in Table 3a (males) and 3b (females).
FIGURE 2a Parameter estimates and standard errors for the four C components from the best fit common pathway model (males). Note: C1 = C_close; C2 = C_family; C3 = C_household; and C4 = C_community. Paths are standardized partial regression coefficients. *Suggest that the path is statistically significant at the 5% level.
FIGURE 2b Parameter estimates and standard errors for the four C components from the best fit common pathway model (females). Note: C1 = C_close; C2 = C_family; C3 = C_household; and C4 = C_community. Paths are standardized partial regression coefficients.
TABLE 3a Variance Components from the Best-Fitting Model (Males)
TABLE 3b Variance Components from the Best-Fitting Model (Females)
Males
First, all behaviors loaded quite strongly on the common factor, with the loading for DA (+0.85) being slightly higher than that seen for AUD (+0.80) and CB (+0.74). Similar to but somewhat more pronounced than with the univariate model, estimates for total shared environmental contributions (C_close + C_family + C_household + C_community) were higher for DA than for AUD and CB (c_close2+c family2+c household2+c community2 = 0.21 [DA]; 0.17 [CB]; and 0.16 [AUD]). Also similar to the univariate models, among the four components of C, community experiences were most important for DA (c_community2 = 0.10), while for CB and AUD, the close family environment were of more importance (c_close2 = 0.08 [CB] and c_close2 = 0.05 [AUD]). More similar estimates, and also more modest, were seen across syndromes for the two c components in between community experiences and close family experiences (c_family, c_household). This was slightly different from what was seen in the univariate models. Third, syndrome-specific community effects were mainly seen for DA and partly for CB, while syndrome-specific close family effects were mainly seen for AUD and CB. The two C components in between community experiences and close family experiences (c_family, c_household) were estimated to almost entirely derive from the common latent factor. Note that for all syndromes, the effect of additive genetics was 2–4 times the size of the total shared environmental effect.
Females
As for males, we first fitted an independent pathway model that resulted in a log-likelihood and AIC of 1554073 and -9569501, respectively. The common pathway model had a lower AIC (-9569773) and we therefore focused on the results from this model.
All behaviors loaded strongly on the common factor, but the loading for DA (+0.92) was substantially higher than that seen for CB (+0.63) and AUD (+0.74). Second, as seen in the univariate models, estimates for total shared environmental contributions (C_close + C_family + C_household + C_community) were substantially lower than for males. They were higher for DA (c_close2+c family2+c household2+c community2= 0.18 [DA] than for either CB [0.10] or AUD [0.14]). Among the four components of c, the community experiences was just slightly above zero while the c components referred to as the household effect were almost zero. The total close family experiences were similar and most important across syndromes and also divided into common and specific components. For c_family, the total effect was strongest for DA (0.09) and it was only for DA that there was a syndrome specific effect of c_family. For all syndromes, the effects of additive genetic factors were 2–4 times the size of the total effect of the shared environment.
Discussion
We sought in these analyses, using information available to us from a range of Swedish national registries, to sub-divide the latent concept of ‘shared environment’ or C into four subcomponents that we termed (1) close family, (2) family, (3) household, and (4) community. We applied our model to three externalizing syndromes that we had previously studied and for which we had detected substantial evidence for shared environment. We first fitted univariate models and then a common pathway model to males and females separately, which permitted us to determine the degree to which these C sub-components were common across the three syndromes versus syndrome specific and were similar in magnitude in males and females.
Examining first the influences common to all three externalizing syndromes, in both males and females, all the shared environmental components had only modest effects. But among them, in males, family environment had the largest effect, followed by household and community with close family effects being the weakest. The pattern was different in females. Family effects remained the strongest but was then followed by close family with household and community effects being considerably weaker.
With regard to syndrome-specific, shared environmental effects, in males, family and household effects were consistently either absent or very small in magnitude. By contrast, close family effects were relatively strong for CB and AUD, and weaker for DA. However, community effects were relatively strong for DA and weaker for CB and AUD. The pattern was different in women. By far the strongest impact was seen for close family effects specific to CB, followed by family effects also for CB. Family effects also impacted modestly on AUD and DA with all other effects being quite small in magnitude.
These patterns were similar in the univariate models. However, as seen from the univariate models, where confidence intervals could be estimated, the uncertainty of the estimates for the different C components was rather high.
These results are broadly consistent with several prior results we have obtained for Swedish samples using different analytic methods. First, congruent with our current findings, in a prior traditional common pathway twin analysis of DA, CB, and AUD (Kendler, Lonn et al., Reference Kendler, Lonn, Maes, Lichtenstein, Sundquist and Sundquist2016), we found shared environmental influences on the common liability and specific liabilities to CB and DA to be considerably stronger in males than females. Second, using multilevel modeling controlling for familial factors, we found that community level influences had a considerably stronger influence on DA in males than females (Kendler, Ohlsson et al., Reference Kendler, Ohlsson, Sundquist and Sundquist2015). We found exactly this pattern for our DA-specific community effects in these analyses. Third, using conditional logistic models, resemblance in male and female sibling pairs for both DA (Kendler, Ohlsson et al., Reference Kendler, Ohlsson, Sundquist and Sundquist2013) and CB (Kendler, Morris et al., Reference Kendler, Morris, Lönn, Sundquist and Sundquist2014) were greater for those closer versus more distant in age. These results reflect what we would here call our ‘close family’ effects — biological sibling pairs of close age raised together — that we see to be especially prominently for CB in these analyses. Fourth, using Cox proportional hazard regression models, we found in male and female sib-pairs that the syndromes of both DA and CB were ‘transmitted’ considerably more strongly from older to younger than younger to older siblings (Kendler, Morris et al., Reference Kendler, Morris, Lönn, Sundquist and Sundquist2014; Kendler, Ohlsson et al., Reference Kendler, Ohlsson, Sundquist and Sundquist2013). In our current conceptual framework, these effects would likely be reflected in both ‘close family’ and ‘family’ effects. Fifth, we found in male–male full- and half-sibling pairs resemblance in risk for DA was predicted by years of living together in the same household (Kendler, Maes et al., Reference Kendler, Maes, Sundquist, Ohlsson and Sundquist2013). These results would also reflect what we call close family and family effects. These same analyses found slightly weaker but significant effects on pair resemblance for years of residing in the same community but not the same household — what we would here call community effects (Kendler, Maes et al., Reference Kendler, Maes, Sundquist, Ohlsson and Sundquist2013). Finally, congruent with our findings, we examined the same two effects in female–female half- and full-sibling pairs and found that they were present but weaker than those observed in the male–male pairs (Kendler, Maes et al., Reference Kendler, Maes, Sundquist, Ohlsson and Sundquist2013). Our results are also consistent with findings from Burt et al. (Reference Burt, Klump, Kashy, Gorman-Smith and Neiderhiser2015), who sampled communities across Michigan and found that 11% of the population variance in self-reported non-aggressive rule-breaking — a broad index of an externalizing predisposition — was predicated from the participant's neighborhood.
We found that community effects on externalizing syndromes, especially DA, were larger in males than females. This is consistent with evidence that compared to females, males are more motivated to use psychoactive substances to conform to subgroup values and are more peer influenced in their intake of drugs (Borsari & Carey, Reference Borsari and Carey2001; Graziano et al., Reference Graziano, Bina, Giannotta and Ciairano2012). In a relevant Swedish high-school student survey, Svensson (Reference Svensson2003) found that males compared to females had higher exposure rates to deviant peers. This arose because parents of girls compared to parents of boys monitored their behavior and friends more closely and that the probability of drug use was more strongly predicted by peer deviance in males than in females (Svensson, Reference Svensson2003).
Our results are at the latent level and so do not illuminate the more specific social mechanisms through which our components of shared environment actually impact on resemblance in relatives. A few speculations, however, are in order. First, we suspect that the close family effects are likely driven largely through peer effects. Close-aged siblings are much more likely to socialize together and share peer groups than are siblings distant in age with one another. It would also be likely that close-aged siblings would use drugs or alcohol together or one would provide the substance (i.e., ‘turn on’) the other. For CB, close-age siblings are more likely to collaborate in criminal activity, which might be responsible for the substantial close family-specific effects for CB in our analyses. Second, we suspect that family effects for externalizing syndromes likely reflect parental characteristics that have broadly consistent effects across all of their offspring including SES, parenting styles, and levels of parental monitoring, all of which have been shown to impact on risk for externalizing syndromes in children (Bleuler, Reference Bleuler and Diethelm1955; Farrington, Reference Farrington2005; Hawkins et al., Reference Hawkins, Herrenkohl, Farrington, Brewer, Catalano, Harachi, Loeber and Farrington1998; Kendler et al., Reference Kendler, Myers and Prescott2000; Otowa et al., Reference Otowa, Gardner, Kendler and Hettema2013). Third, our distinction between family and household effects was perhaps the most speculative part of our model. This was meant to capture the difference between biological and step-siblings living in the same household with the expectation that biological siblings would experience more similar environments. Our analyses provided some evidence that our hypothesis was correct. Fourth, community effects could result from differences in the general peer environment in the community and particularly in the schools and the degree to which deviant behaviors in adolescence were encouraged. For example, a study in nearby Finland documented substantially different rates of smoking and drinking in different schools (Rose et al., Reference Rose, Viken, Dick, Bates, Pulkkinen and Kaprio2003). For DA, an additional potentially important feature of community effects is likely availability and price of illicit substances, which might be driving the DA-specific community effects.
In summary, in line with recent reviews by Burt (Reference Burt2009, Reference Burt2014), our results confirm both the overall importance of the shared environment for broadly defined psychopathology and the potential progress that can be made in understanding the nature of these influences by decomposing c2 into its constituent parts.
Limitations
These results should be interpreted in the context of several potential methodological limitations. First, optimization issues in OpenMx made it infeasible to obtain confidence intervals for the parameter estimates in the common pathway model. Still, for males, we presented standard errors of the path coefficients. For females, we encountered problems in finding the global likelihood maximum. Focusing on finding the best fit for the models, we ran each model from a variety of starting values, and reran them from the previously fitted values, to ensure that a global likelihood maximum had been reached. This made the standard errors for some parameters unreliable; hence, we did not report them. Note, however, that we are using complete population data; consequently, the precision of the estimates cannot be improved. Second, the definition and the assumptions in the models for the four C components in the model could be questioned. The geographical boundaries of the SAMS, while generally set to reflect natural community boundaries, were set for administrative reasons. It is possible that other criteria, such as municipalities, school districts, or geographic units centered on individuals’ residences, might have done a superior job of capturing community-level influences. Despite these shortcomings, the SAMS units have previously been used successfully to capture a collective influence and, as seen here particularly for males, it also captures a substantial collective influence. While the household is a clearly defined unit, the separation between close family (C1) and family (C2) could be questioned. However, for all six within-syndrome correlations, the full-siblings closer in age had a higher correlation than full-siblings further apart in age. Also, we did not include a T parameter in the model, which would capture the difference between close-in-age full-siblings and DZ pairs. Nevertheless, comparing the correlations between DZ pairs and full-sibling pairs do not support a T effect as for three out of six within-syndromes correlations the full-sibling correlation were higher than the DZ correlations. Third, we defined the shared environments as any environmental, that is, non-genetic, influence that increases similarity regardless of the proportion of genes shared. Hence, we then assume that pairs do not create or provoke more similar environments for themselves because they are more similar genetically. However, for household, we allowed for a difference between biological and step-siblings living in the same household with the expectation that biological siblings would experience more similar environments. We did not allow this for the community effect (C4). We partly tested this by only analyzing non-related individuals and estimated a model with only the household (C3) and the community (C4) effect. The community effect (measured as the share of total variance) was identical over the two models. Fourth, we used several different registers to construct the three syndromes, and the time frame during which data for CB, AUD, and DA are available is different for different registries, which might explain why full-siblings and half-siblings further apart in age are less similar. However, the difference in the within-syndrome correlations when comparing full-siblings close in age and full-siblings further apart in age is similar when we investigate those born 1980–1990 and are exposed during the period when we have the most information from all registers included. Fifth, while using registry data has the important advantage of not requiring accurate respondent recall and reporting, the risk for time-dependent bias due to registry effects remains. In addition, only a subset of individuals who are drug abusers (or alcohol abusers) would ever require medical treatment or commit an alcohol or drug-related crime, which means that our registries alone do not capture the true prevalence in the population. However, such data likely contains both false negative and false positive diagnoses, the frequency of which we cannot estimate.
Acknowledgments
This project was supported by the grants R01AA023534 and F32AA02226 from the National Institute of Alcohol Abuse and Alcoholism, the Swedish Research Council (K2012-70X-15428-08-3), the Swedish Research Council for Health, Working Life, and Welfare (In Swedish: Forte; Reg.nr: 2013-1836), the Swedish Research Council (2012-2378; 2014–10134) and FORTE (2014-0804) as well as ALF funding from Region Skåne awarded. None of the authors have any conflicts of interest to report.
Ethical Standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant national and institutional committees on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008. We secured ethical approval for this study from the Regional Ethical Review Board of Lund University (No. 2008/409).
Appendix A
DA was identified in the Swedish medical (the Swedish Hospital Discharge Register, containing all hospitalizations for all Swedish inhabitants from 1973–2010 and the Outpatient Care Register, containing information from outpatient clinics from 2001 to 2010) and mortality registries by ICD codes (ICD-8: drug dependence [304]; ICD-9: drug-induced psychoses [292] and drug dependence [304]; ICD-10: mental and behavioral disorders due to psychoactive substance use [F10–F19], except those due to alcohol [F10] or tobacco [F17]); in the Suspicion Register, which records arrests by codes 3070 (driving under the influence of drugs), 5010 (drug possession), 5011 (drug use), and 5012 (possession and use), which reflect crimes related to drug abuse; and in the Crime Register, which reflects convictions by references to laws covering narcotics (law 1968:64, paragraph 1, point 6) and drug-related driving offenses (law 1951:649, paragraph 4, subsection 2 and paragraph 4A, subsection 2). Drug abuse was identified in individuals (excluding cancer patients) in the Prescribed Drug Register who had filled prescriptions of hypnotics and sedatives (Anatomical Therapeutic Chemical [ATC] Classification System N05C and N05BA) or opioids (ATC: N02A) in dosages consistent with drug abuse, that is, on average more than four defined daily doses a day for 12 months.
CB was defined from the Swedish Crime Register based on the following convictions, law and chapter in parentheses: (aggravated) assault (3:5, 3:6); illegal threat (4:5); threats and violence against an officer (17:1, 17:2); intimidation (4:7); (Gross) violation of a person's/woman's integrity (4:4a); kidnapping (4:1); illegal confinement or restraint (4:2); (aggravated) robbery (8:5, 8:6); illegal coercion (4:4); (aggravated) Arson (13:1, 13:2); murder, manslaughter or filicide (3:1, 3:2, 3:3); sexual crimes (excluding prostitution and the buying of sexual services but including child pornography) (6:1–6:10, 6:12, 16:10A); theft of a vehicle (8:1–2, 8:4, 8:7–8); theft (including burglary) (8:1–2, 8:4); vandalism (12:1–4); vandalism causing danger to the public, sabotage, hijacking (13:3–10 (5a–b)); unlawful entering of a person's home, trespassing (4:6); fraud (9:1–10); embezzlement (10:1–8 (5a–e)); dishonesty/crime towards a creditor (includes forged book-keeping in companies) (11:1–5); and forgery (14:1–10); and thereby excluding convictions for minor crimes like traffic infractions.
AUD was defined by ICD codes for main and secondary diagnoses from Swedish medical registries for the following diagnoses: ICD8 and 9: alcohol-related psychiatric disorders (291), alcohol dependence (303), alcohol abuse (305A), alcohol-related polyneuropathy (357F), alcohol-related cardiomyopathy (425F), alcohol-related gastritis (535D), alcoholic fatty liver, alcohol hepatitis, alcoholic cirrhosis, unspecified liver damage caused by alcohol (571A-D), toxic effects of alcohol (980), alcoholism (V79B); ICD10: alcohol6 related psychiatric and behavioral disorders (F10, excluding acute alcohol intoxication: F10.0), rehabilitation of a person with alcohol abuse (Z50.2), guidance and medical advice to a person with alcohol abuse (Z71.4), alcohol-related pseudo-Cushing syndrome (E24.4), alcohol-related degeneration of the nervous system and brain (G31.2), alcohol-related polyneuropathy (G62.1), alcohol-related myopathy (G72.1), alcohol-related cardiomyopathy (I42.6), alcohol-related gastritis (K29.2), liver diseases caused by alcohol (K70.0-K70.9), acute pancreatitis caused by alcohol (K85.2), chronic pancreatitis caused by alcohol (K86.0), treatment of pregnant alcoholic woman (O35.4), toxic effects of alcohol (T51.0-T51.9), and based on Anatomical Therapeutic Chemical (ATC) codes in the Prescribed Drug Register (containing all prescriptions in Sweden picked up by patients from 2005 to 2010): disulfiram (N07BB01), acamprosate (N07BB03), or naltrexone (N07BB04). Additionally, we identified individuals with at least two convictions of drunk driving (law 1951:649) or drunk in charge of maritime vessel (law 1994:1009) in the Crime registers insuring that we did not count arrests in the suspicion register that descried the same event contained in the conviction register.
Appendix B
Household was ascertained as follows: From 1960 to 1985 (every 5th year), we used the household identification number from the Population and Housing Census. The household identification number includes all individuals living in the same dwelling. For the years we did not have data, we approximated the household identification number with the available information from the year closest in time. From 1986 onward (every year), we used the family identification number from the Total Population Register. Family identification numbers are assigned to individuals who are related and who are registered as residing at the same property (each individual can be part of only one family). In addition, adults who are registered at the same property and have children together but are not married are registered as being in the same family.
Geographical status was defined in terms of Small Areas for Market Statistics (SAMS), which are small geographical units defined by Statistics Sweden, the Swedish government statistics bureau. There are approximately 9,200 SAMS throughout Sweden, with an average population of around 1,000. From 1960 to 1970, we had no information on SAMS areas, and therefore we used parishes as a proxy for SAMS for these years. The parishes serve as districts for the Swedish census and elections, and have approximately the same number of inhabitants as SAMS areas. From 1960 to 1985, we only had information for every 5th year, and for that reason we approximated the geographical status with the information from the year closest in time.
