Abstract
Objective Examine whether the relationship between childhood conduct problems and substance use in adulthood is explained through potential snares (police arrest, gang membership, school non-completion) in a Brazilian population-based birth cohort.
Method Data were analyzed from 4,599 young people from the 1993 Pelotas Birth Cohort in Brazil. The exposure was conduct problems at age 11 years. Outcomes included hazardous alcohol consumption and illicit drug use (age 22 years). Mediators included police arrest (by age 18 years), gang membership (ages 18 and 22 years), and school non-completion (by age 22 years). We performed counterfactual mediation using the parametric g-computation formula to estimate the indirect effect via all three mediators simultaneously.
Results After adjusting for confounders (including comorbidity), conduct problems were weakly associated with police arrest [OR (95% CI)=1.45 (0.97, 2.16)] and school non-completion [OR (95% CI)=1.46 (1.22, 1.74)], but not with gang membership. Police arrest and gang membership were associated with illicit drug use [OR (95% CI)=3.84 (2.46, 5.99); OR (95% CI)=7.78 (4.30, 14.10), respectively] and with hazardous alcohol use [OR (95% CI)=1.60 (1.08, 2.38); OR (95% CI)=1.88 (1.07, 3.30)]. However, there was no association between school non-completion and either outcome after confounder adjustment. There was little evidence for an indirect effect of conduct problems on hazardous alcohol use and illicit drug use via all three mediators after confounder adjustment.
Conclusion The indirect effects of childhood conduct problems on substance use in early adulthood via police arrest, gang membership, and school non-completion were very small, after accounting for confounders and comorbidity.
Introduction
Conduct problems refer to behaviors related to conduct disorder (norm-breaking behaviors and violations of the rights of others) and behaviors related to oppositional defiant disorder (noncompliant, angry and defiant behaviours). In Brazil, conduct problems are almost twice as prevalent as in high-income countries (HIC), with an average prevalence of 21% compared to 13% in the UK and 11% in the US (measured using the Strengths and Difficulties Questionnaire)1, with a higher prevalence even after accounting for measurement differences2. There is robust evidence that childhood conduct problems are associated with various adverse outcomes in early adulthood, including problem behaviors such as alcohol and drug misuse3. Although the majority of evidence comes from HIC, a previous study comparing a British and Brazilian birth cohort found that associations for childhood conduct problems with hazardous alcohol use and illicit drug use in early adulthood were stronger in Brazil compared to the UK2.
Problematic substance use places a significant burden on society as it is linked with alcohol dependency, road traffic accidents, injury, violence, and crime4. To reduce problematic substance use in Brazil, it is important to understand factors that are on the causal pathway from childhood conduct problems. An influential developmental taxonomy published by Moffitt in 1993 introduced the concept of snares that can trap adolescents into long-term problematic behaviors beyond an age when most otherwise desist from these behaviors5. Previous literature has suggested that interrupted education, gang membership and contact with the criminal justice system can act as snares trapping some adolescents into persistent problem behavior, including problematic substance use, when desistance is normative3,6–8.
Very little research has examined the importance of potential snares in low- and middle-income countries (LMIC). This is particularly important to address in countries, such as Brazil, with a higher prevalence of both conduct problems1 and potential snares such as gang membership and school non-completion9,10, compared to HIC. In LMIC settings, the processes for the association between conduct problems and later problematic substance use might differ substantially from HIC due to fewer supportive systems and different social contexts for children with conduct problems1,11.
Considering school non-completion as a possible “snare” linking conduct problems and later substance use, in Brazil, repetitive grade retention is a very common practice12 and is a strong predictor of dropping out of school13. Longitudinal data show Brazilian children repeating a school grade by age 11 had increased odds of elevated conduct problems into mid-adolescence14, and those with an externalising disorder in adolescence had higher odds of grade repetition three years later15. The association between childhood conduct problems and school non-completion is also well-documented in HIC16–20. However, it is not clear if there is a causal relationship between conduct problems and school non-completion, with some studies suggesting it is confounded by inattention19. There is evidence that school non-completion is associated with problematic behaviors in early adulthood, such as substance use21,22, and completing school by age 22 years was associated with lower odds of violent and non-violent crime, even after accounting for school grade repetitions in Brazil23.
Considering another possible “snare” following child conduct problems, youth gangs are more prevalent in Brazil than in most HIC and may offer a form of social capital and protection to adolescents who have failed to comply with societal norms9,24. Excessive alcohol, drug use and interpersonal violence are common among gang members and becoming a member increases the chances of early initiation and frequent consumption of alcohol and illicit drugs25. There is evidence that childhood conduct problems are a risk factor for gang membership in both HIC and LMIC9,24,26. Finally, contact with the criminal justice system can be a consequence of childhood conduct problems, given the strong links between childhood conduct problems and later crime in LMIC2,27. According to labelling theory, stigmatising effects of obtaining a criminal record are internalised by young people, resulting in self-stigma, low self-efficacy and self-esteem, increasing susceptibility to long-term problematic alcohol and drug use28.
One of the challenges for understanding relationships between child conduct problems, potential snares (such as gang membership) and later substance use, is the possibility of bidirectional relationships and shared or distinct sets of confounders for each relationship. Addressing these challenges requires data from a study where young people are followed longitudinally across childhood, adolescence, and early adulthood, with potential confounders assessed across multiple domains. In the current study, we examine whether potential snares in late adolescence (school non-completion, police arrest, and gang membership) mediate the association between childhood conduct problems and hazardous alcohol use and illicit drug use in early adulthood using a population-based, prospective birth cohort in Brazil: the 1993 Pelotas Birth Cohort.
In a previous study using the 1993 Pelotas Birth Cohort, conduct problems at age 11 years were associated with hazardous alcohol use and illicit drug use at age 22 years, after accounting for sociodemographic, individual, peer, family, and community confounders2. In this study, we hypothesise that associations between conduct problems at age 11 years and hazardous alcohol use and illicit drug use at age 22 years will be partly explained through potential snares between age 18 and 22 years (including school non-completion, police arrest, and gang membership). We hypothesise that potential baseline and intermediate confounders will weaken these indirect effects, particularly when accounting for childhood hyperactivity on the associations for conduct problems2,19,29,30.
Methods
Sample
The 1993 Pelotas Birth Cohort Study is an ongoing population-based study designed to investigate the effects of a wide range of influences on health and development. Pelotas is a city located in the extreme south of Brazil, with an estimated population of 345,000 inhabitants, 93% of whom live in the urban area. All births occurring in the five maternity clinics in the town were monitored in 1993 (99% of births in Pelotas occurred in hospital). For the 5,265 children born alive, only 16 mothers could not be interviewed or refused to participate in the study. The 5,249 newborns, whose mothers lived in the urban area, were included in the cohort (81 were either twins or triplets). The detailed methodology of this study can be found elsewhere31,32. During the perinatal study, mothers were interviewed to collect demographic, health and socioeconomic information about the family. Follow-up visits were conducted in 2004–2005 (age 11; N = 4,452 mothers; retention rate of 87.5%, after accounting for N = 141 deaths), 2008 (age 15; 4,349 mothers; retention rate of 85.7%, after accounting for N = 147 deaths), 2011-2012 (age 18; 4,106 mothers; retention rate of 81.4%, after accounting for N = 164 deaths), and 2015-2016 (age 22; N = 3,810 young people; retention rate of 76.3%, after accounting for N = 193 deaths)31,33. Study data at age 22 years was collected and managed using REDCap electronic data capture tools34. The perinatal study and each follow-up were approved by the Research Ethics Committee of the Federal University of Pelotas School of Medicine. After being informed of the details of the study, participants signed a term of informed consent.
Measures
Exposure: Conduct problems at age 11 years
Conduct problems were assessed at age 11 during an interview with the primary caregiver (usually mothers) using the Strengths and Difficulties Questionnaire (SDQ)35. The SDQ is a screening questionnaire for children’s mental health problems that occurred in the previous six months, and has been validated in Brazil36,37 using independently diagnosed psychiatric disorders. Five items form the conduct problems subscale, including: ‘often has temper tantrums or hot tempers’, ‘often fights with other children or bullies them’, ‘often lies or cheats’, ‘steals from home, school or elsewhere’ and ‘generally obedient, usually does what adults request’. Responses for each item are ‘Not true’, ‘Somewhat true’ and ‘Certainly true’. The Item ‘generally obedient’ was reverse coded. All items were summed together (range 0 to 10) and the scale was dichotomised to create a binary exposure of low conduct problems (less than four) versus high conduct problems (four or more) according to the cut point on the newer 4-band categorisation (https://www.sdqinfo.org/py/sdqinfo/c0.py).
Outcomes: Hazardous alcohol use and illicit drug use at age 22 years
The 10-item self-report Alcohol Use Disorder Identification Test (AUDIT38) was used to assess hazardous alcohol consumption at age 22 years. The AUDIT is a brief screening tool aiming to detect risky drinking and risk for alcohol dependence with high validity and reliability39 and has been validated using diagnostic interview, physical examinations, and laboratory testing40. The AUDIT was dichotomised at a cut-point of eight and treated as a binary variable. The cut-point of eight or more represents hazardous levels of drinking and has been validated in Brazil41 using psychiatric diagnoses of alcohol use disorders.
Illicit drug use was assessed using the same self-report questionnaire at age 22 years which included questions about lifetime use of cannabis, cocaine, crack, amphetamine-type stimulants, nitrous oxide or other inhalants, hallucinogens, opioids, and other injected illegal drugs. A binary variable was created representing current drug use with “no” representing ‘never used’ (61%), ‘I just tried it’ (17%), and ‘I used to use it but I don’t anymore’ (8%) and “yes” representing ‘I use it occasionally’ (7%), ‘I only use it on weekends’ (1%), and ‘I use it every day or almost every day’ (5%).
Mediators: Police arrest (age 18), gang membership (age 18-22), and school non-completion (22 years)
Police arrest was assessed by self-report at the age of 18, asking if the young person had ever been arrested or detained (“yes” if ever arrested or detained and “no” otherwise).
Gang membership was assessed by self-report at ages 18 and 22 years with questions asking whether the participant had been a gang member in the previous year. Both time points were combined into a binary variable representing gang membership at either time point due to the low prevalence of gang membership (approximately 1% at both time points).
School non-completion (assessed at age 22) was a dichotomous variable indicating if the young person has finished high school (passed all school grades) by age 22 (“yes” if not completed school and “no” if completed). In Brazil, young people typically finish school by age 17, assuming they do not fail any grades; however, a student may continue studying into adulthood until they complete all grades23.
Baseline and intermediate confounders
Baseline confounders and intermediate confounders were drawn from five domains (sociodemographic, individual, peer, family, and community) that have previously been identified as key factors influencing a young person’s behavior (see Supplement 1, available online). Baseline confounders included sex, a score of sociodemographic (e.g., low maternal education) and health (e.g., mother smoking in pregnancy) risk factors for conduct problems (assessed by maternal report or observation during the perinatal period), maternal depression, the quality of the father-child relationship, the quality of the mother-child relationship, parental alcohol consumption, parental smoking status, parental separation, and neighbourhood safety (all assessed at child age 11 years). Intermediate confounders included peer deviance (assessed at child age 11 years), frequency of cigarette smoking and alcohol consumption, and peer drug use (all assessed at child age 15 years).
In secondary analyses, hyperactivity problems (at age 11 years) were considered as an additional baseline confounder and emotional problems (at age 15 years) were considered as an additional intermediate confounder.
Statistical analyses
Prior to conducting any statistical analysis, we constructed a Directed Acyclic Graph (DAG) to present hypothesised causal pathways based on previously published articles (see Figure 1). All statistical analyses were carried out using Stata version 1742. First, we performed descriptive statistics to report characteristics of the population for the whole sample. Second, we conducted unadjusted and adjusted logistic regressions in stages to examine exposure-mediator and mediator-outcome associations. Third, we performed a counterfactual approach to mediation using the parametric g-computation formula to estimate indirect effects using binary mediators and common binary outcomes, and incorporate intermediate confounders43. The package -gformula-in Stata44 was used to estimate the total causal effect, natural direct effect and natural indirect effect via all three mediators simultaneously. All mediators were included simultaneously as we did not consider them to be independent. However, given that we were not able to confidently specify a causal direction between, we were only able to estimate one indirect effect via all three mediators together. -gformula-uses Monte Carlo simulations to simulate the mediator, outcome, and intermediate confounders under each hypothetical “counter to the fact” scenario. Standard errors were estimated using 50 bootstrap samples and normal-based 95% confidence intervals were calculated. More detail on the statistical analyses can be found in Supplement 2, available online.
Sensitivity analyses
First, we re-ran the mediation models making the assumption that police arrest (by age 18 years) and gang membership (at age 18 or 22 years) caused school non-completion (by age 22 years) i.e. we performed a single mediator model for school non-completion with police arrest and gang membership as additional intermediate confounders. Second, we re-ran all analyses using complete case rather than imputed data (further details provided below).
Missing data
Of the 5,249 newborns included in the cohort, 45% (N = 2,342) had complete data on all analysis variables (see Figure S1, available online, for a flow chart of retention). Table S1, available online, provides a comparison of the complete case sample to the sample with missing information in at least one of the analysis variables (n = 2,907) on the exposure and baseline confounders. As shown in Table S1, high conduct problems at age 11 was associated with missingness, as was female sex, a higher sociodemographic and health risk score, higher maternal depression symptoms, no parental separation, no parental alcohol use, poorer father-child and mother-child relationship, and higher hyperactivity problems.
To increase sample size, we performed multivariate imputation by chained equations to impute missing data on all variables, across all assessments up to our starting sample (those that provided data on alcohol use and drug use from at least one time point between age 11 and age 22 years; N = 4,599). This starting sample was chosen to ensure we had strong auxiliary variables available to impute substance use at age 22 years. For a detailed description of multiple imputation, see Supplement 3, available online. Imputation and subsequent analysis code for the mediation models can be found on GitHub: https://github.com/gemmahammerton/gformula_1993_Pelotas. All analyses in the main text were performed on imputed data.
Results
Descriptive statistics
Of 4,599 children included in the imputed analyses, 51% were female, 77% had White mothers, and 23% had Black or mixed race mothers. Table 1 shows descriptive statistics for all analysis variables for the whole sample, using both imputed (N = 4,599) and complete case data (N = 2,342). Based on imputed data, 31% of children at age 11 years had high conduct problems, 4% of young people had been arrested or detained by the police by age 18 years, 2% had been members of a gang in the previous year at age 18 or 22 years, and 43% had not completed high school (passed all school grades) by age 22 years. 22% of young people reported hazardous alcohol consumption at age 22 years, and 14% reported current use of illicit drugs.
Associations between conduct problems (exposure) and police arrest, gang membership and school non-completion (mediators)
Figure 2 shows associations between the exposure (high conduct problems) and the mediators (police arrest, gang membership, and school non-completion) before and after adjusting for confounders. After adjusting for baseline confounders and hyperactivity problems, those with high childhood conduct problems had slightly higher odds of being arrested by the police by age 18 [OR (95% CI) = 1.45 (0.97, 2.16)] and not completing school by age 22 [OR (95% CI) = 1.46 (1.22, 1.74)] compared to those with low childhood conduct problems. There was no evidence of an association between conduct problems and gang membership after adjusting for baseline confounders and hyperactivity problems [OR (95% CI) = 1.17 (0.70, 1.95)].
Associations between police arrest, gang membership and school non-completion (mediators) and hazardous alcohol use and illicit drug use (outcomes)
Univariable associations for baseline confounders with mediators and outcomes are shown in Table S3, available online, and univariable associations for intermediate confounders are shown in Table S4.
Table 2 shows associations between the mediators and the outcomes before and after adjusting for confounders (including hyperactivity and emotional problems). After adjusting for all confounders, those who had been arrested [OR (95% CI) = 1.60 (1.08, 2.38)] and those in a gang [OR (95% CI) = 1.88 (1.07, 3.30)] had slightly higher odds of hazardous alcohol use at age 22. There was no evidence of an association between school non-completion and hazardous alcohol use after adjusting for confounders [OR (95% CI) = 1.10 (0.90, 1.33)].
After adjusting for confounders, those who had been arrested [OR (95% CI) = 3.84 (2.46, 5.99)] and those in a gang [OR (95% CI) = 7.78 (4.30, 14.10)] had higher odds of illicit drug use at age 22. There was no evidence of an association between school non-completion and illicit drug use after adjusting for confounders [OR (95% CI) = 1.17 (0.92, 1.49)].
Indirect effects of conduct problems (exposure) on hazardous alcohol use and illicit drug use (outcomes) via police arrest, gang membership and school non-completion (mediators)
Table 3 shows the total causal effect, natural direct effect and natural indirect effect of conduct problems on hazardous alcohol use and illicit drug use via all three mediators simultaneously. After adjusting for baseline and intermediate confounders, there was weak evidence for a small total effect of conduct problems on hazardous alcohol use [OR (95% CI) = 1.18 (1.00, 1.39)] but little evidence for a total effect of conduct problems on illicit drug use [OR (95% CI) = 1.05 (0.86, 1.30)]. There was weak evidence for a small natural indirect effect of conduct problems on hazardous alcohol use [OR (95% CI) = 1.03 (1.00, 1.06)] and illicit drug use [OR (95% CI) = 1.08 (1.01, 1.17)] via all three mediators.
After additionally adjusting for hyperactivity (at age 11) as a baseline confounder and emotional problems (at age 15) as an intermediate confounder, all effects weakened. There was no longer evidence for a total effect of conduct problems on hazardous alcohol use [OR (95% CI) = 1.02 (0.84, 1.24)], or a natural indirect effect of conduct problems on hazardous alcohol use [OR (95% CI) = 1.01 (0.99, 1.04)] or illicit drug use [OR (95% CI) = 1.03 (0.96, 1.10)] via all three mediators.
Sensitivity analyses
Table S5, available online, shows the total causal effect, natural direct effect and natural indirect effect of conduct problems on hazardous alcohol use and illicit drug use via school non-completion, treating police arrest and gang membership as intermediate confounders. There was no evidence for an indirect effect via school non-completion, for either hazardous alcohol use or illicit drug use after adjusting for confounders.
Tables S6 to S8, available online, show results for all analyses using complete case data (N = 2,342). In general, odds ratios were slightly stronger but confidence intervals were wider using the complete case sample; however, conclusions were unchanged.
Discussion
Summary of findings
Using a large, population-based birth cohort in Brazil, we found that high childhood conduct problems were weakly associated with being arrested by the police by age 18 years and strongly associated with not completing school by age 22 years after adjusting for potential confounders, including childhood hyperactivity problems. We also found an association between childhood conduct problems and gang membership in late adolescence, but this weakened substantially after accounting for childhood hyperactivity problems. Being arrested by the police and gang membership were strongly associated with illicit drug use and weakly associated with hazardous alcohol use at age 22 years, after adjusting for confounders (including childhood conduct, hyperactivity, and emotional problems). However, there was no association between school non-completion and either hazardous alcohol use or illicit drug use after adjusting for confounders. Results from the mediation analysis showed weak evidence for a total effect of conduct problems on hazardous alcohol use but little evidence for a total effect of conduct problems on illicit drug use after adjusting for confounders. We found weak evidence for a natural indirect effect of conduct problems on hazardous alcohol use and illicit drug use via all three mediators; however, after additionally adjusting for hyperactivity and emotional problems, there was no longer evidence of an indirect effect on either outcome.
Strengths and limitations
Due to the longitudinal design and wealth of data in the 1993 Pelotas Birth Cohort, we were able to examine mediators of the association between conduct problems at age 11 and substance use over 10 years later, whilst adjusting for potential baseline and intermediate confounders across multiple domains, including sociodemographic, individual, peer, family, and the community. However, our results need to be interpreted in light of several limitations. First, in our population-based sample, few people joined a gang or were arrested by the police (2% and 4%, respectively); therefore, the weak evidence for associations could be explained by a lack of power and other factors more common among children with conduct problems might be more important mediators for the relationship between conduct problems and substance use.
Second, even though we adjusted for a wide range of confounders, there might be some residual confounding present in our analysis due to either unmeasured confounders (such as genetic risk) or confounders measured with error. Third, as with most cohort studies, there was selective attrition over time, with missingness associated with conduct and hyperactivity problems, alongside sociodemographic and family-based confounders.
However, we used multiple imputation to address missing data incorporating a wide range of strong auxiliary variables to make the missing-at-random assumption plausible and conclusions using complete case and imputed data were similar.
Fourth, although the exposure (conduct problems) was assessed using parent-report, the mediators and outcomes were assessed using self-report meaning that results could be subject to reporting bias and shared-rater bias. Additionally, although the exposure and outcomes were assessed using established scales that have been validated in Brazil, each mediator was assessed using a single item. Finally, the measure of school non-completion captures both dropping out of school and grade repetition. It may be that school dropout is a more important risk factor for substance use than grade repetition, but we were not able to distinguish between these related constructs.
Comparison with previous literature
We found that childhood conduct problems were strongly associated with not completing school by age 22 years after adjusting for confounders, including childhood hyperactivity problems. This finding partially supports a previous study using a school-based community cohort of Brazilian children and adolescents (aged 6 to 14 years) which found that externalising disorders (including attention deficit and hyperactivity, conduct and oppositional-defiant disorders) were associated with grade repetition and age-grade distortion three years later15, although not with school dropout. This finding also supports evidence from HIC showing an association between childhood conduct problems and school non-completion16–20. Contrary to the conclusion in a recent systematic review19, we found evidence of an association between conduct problems and school non-completion, even after adjusting for hyperactivity problems. However, although high childhood conduct problems were associated with both police arrest and gang membership after adjusting for confounders, these associations weakened after adjusting for childhood hyperactivity problems, particularly for gang membership. This supports previous reviews of the literature which highlight the importance of risk-taking, impulsivity, and low self-control as risk factors for joining a gang9,24,45. However, any one risk factor in isolation is unlikely to be sufficient to identify young people at high or low risk of joining a gang; alternatively, the total accumulation of risk factors or exposure to risk factors across multiple domains is important24,46,47.
In the current study, there was no association between school non-completion and either hazardous alcohol use or illicit drug use after adjusting for confounders. This is in contrast to a study of American adolescents which found an association between school failure and alcohol use one year later21. The lack of association in the current study could be explained by the wide range of confounders adjusted for (given that there was evidence for an association in the unadjusted analyses). Alternatively, it could be due to the age of the young people, given that those who have completed school may be attending university at age 22, where levels of substance use are generally high. Despite the lack of association between school non-completion and substance use, it is still likely that school non-completion is an important snare in Brazil, that may trap those with childhood conduct problems into involvement with violent and non-violent crime in early adulthood23.
In contrast to the findings for school non-completion, we found that being arrested by the police and gang membership were weakly associated with hazardous alcohol use, and strongly associated with illicit drug use at age 22 years, even after adjusting for confounders (including childhood conduct, hyperactivity, and emotional problems). This supports results from previous studies that found associations between gang membership9,25,48,49 and police arrest28 with problematic substance use.
Conclusions and future directions
The findings from the current study highlight the importance of school professionals being aware of the risk for school non-completion for those with childhood conduct problems, and collaboration between education and health sectors in Brazil. Early intervention to support those with childhood conduct problems to continue with their education is a priority. However, it is also important to consider that many behavioral problems have to do with not adapting to the school system, which is not a problem in itself for the child or adolescent, but our inability to think about education in other ways. Joining gangs can be a way of protesting, finding peers who may also not agree with the ideas of others, of the institutions.
The findings also highlight the importance of early intervention and treatment for hyperactivity (in addition to conduct problems) to prevent long-term negative outcomes such as gang membership, criminal activity, and substance use in Brazil. Given the high rates of comorbidity of conduct and hyperactivity problems, it is difficult to examine the independent effects and preventing or treating both is a priority17,29. Finally, given the strong associations observed between gang membership and police arrest with substance use (particularly illicit drug use), programmes and policies designed to reduce substance use in Brazil should focus on young people involved in gangs, and in the criminal justice system.
In the current study, we found little evidence for an indirect effect of childhood conduct problems on substance use in early adulthood via school non-completion, gang membership, and police arrest after accounting potential confounders including comorbidity. The weak evidence for mediating effects might be because the causes of problematic substance use in children with conduct problems operate early in life6. Alternatively, there may be other mediators, such exposure to trauma or the influence of deviant peers, which might have a stronger influence on later substance use50. Future research is needed to examine mediators of the associations between childhood hyperactivity problems and later substance use in Brazil, and to examine the potential snares that may trap children with conduct problems into other negative long-term outcomes such as criminal behavior and unemployment.
Data Availability
Access to data from the 1993 Pelotas Birth Cohort is through a system of managed open access (https://epidemio-ufpel.org.br/cohort-1993/).
Funding
Dr. Hammerton is supported by the National Institute for Health Research (NIHR) Health Protection Research Unit in Behavioural Science and Evaluation (HPRU BSE) at the University of Bristol, in partnership with the UK Health Security Agency (UKHSA). This research was also funded by a Sir Henry Wellcome Postdoctoral Fellowship (209138/Z/17/Z) awarded to Dr. Hammerton. Drs. Hammerton and Heron are members of the MRC Integrative Epidemiology Unit at the University of Bristol (MC_UU_00011/7). This research was funded in whole, or in part, by the Wellcome Trust (209138/Z/17/Z; 210735_A_18_Z). For the purpose of Open Access, the authors have applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission. Professor Murray is supported by a Wellcome Trust Investigator Award (210735_A_18_Z). Professors Menezes, Gonçalves, Wehrmeister, and Murray are supported by the Brazilian National Research Council (CNPq).
Drs. Heron and Hammerton, PhD, served as the statistical experts for this research.
Supplement 1. Details on the baseline and intermediate confounders
Baseline and intermediate confounders
Baseline confounders included sex, a score of sociodemographic and health risk factors for conduct problems, maternal depression, the quality of the father-child and mother-child relationship, parental alcohol consumption, parental smoking status, parental separation, and neighbourhood safety. These were chosen based on evidence from previous literature of associations with conduct problems 1–4, the potential snares 2,4–8, and substance use 9–13.
Sociodemographic and health risk factors were assessed by maternal report during the perinatal period. Sociodemographic factors included maternal age (< 20 years/≥20 years), low maternal education (yes/no; referring to 0–8 vs≥ 9 years of schooling), marital status (single mother/ with partner), three or more siblings (yes/no) and family income (lowest quintile/ second-fifth quintiles). The cumulative number of sociodemographic risk factors was summed, up to five, for each child as has been done previously 2, with five indicating the highest risk for conduct problems and substance use. Health risk factors included unplanned pregnancy (yes/no), mother ever smoked in pregnancy (yes/no), mother used alcohol in pregnancy (yes/no), maternal urinary infection during pregnancy (yes/no), intrauterine growth restriction (yes/no; referring to < 10th percentile/≥ 10th percentile for gestational age and sex, according to the reference curve developed by Kramer and colleagues (Kramer et al, 2001) and premature birth < 37 weeks (yes/no). Again, the cumulative number of health risk factors was summed, up to six, for each child as has been done previously 2.
All other baseline confounders were assessed when the child was 11 years old. Maternal depression was assessed with mothers using the Self Report Questionnaire, which has been validated in Brazil 14. The quality of the mother-child and father-child relationship was assessed with child-report questions including the following response options: bad, regular, good, very good, excellent. Parental alcohol consumption was assessed using two maternal-report questions, capturing whether she and/or her husband never, not currently or currently drink alcohol. Likewise, parental smoking was assessed using two maternal-report questions, one addressing the mother’s lifetime smoking status and one addressing the father’s lifetime smoking status (ever versus never). Parental separation was assessed using a maternal-report question about the child’s biological father living in the same household, and a child-report question about parental divorce. A binary variable was created, were either child-report or maternal-report of the father not living in the same household were classed as parental separation. Neighbourhood safety was assessed using a maternal-report question asking whether she was afraid of living in the neighbourhood (yes/no).
Intermediate confounders included peer deviance, frequency of cigarette smoking and alcohol consumption, and peer drug use. Intermediate confounders were chosen based on evidence that they confound the associations between the potential snares and substance use 5–8,10,12,15,16, and are also causally downstream from childhood conduct problems (our exposure) 5,12,17–19.
Peer deviance was assessed at age 11 by asking the young person’s mother if she considered her child’s friends as having a bad influence or being bad company (yes/no). Frequency of cigarette smoking and frequency of alcohol consumption were assessed at age 15 by child-report. The questionnaire item for smoking had five response categories: ‘Never smoked cigarettes, not smoked in the last month, 1 to 5 days, 6 to 9 days, ten or more days and every day in the last month’ and the questionnaire item for alcohol had five response categories: ‘Never drank alcohol, not drank alcohol in the last month, one to five days, six to nine days, ten or more days and every day in last month’. For both items, the last three response categories were combined to the category “current smoking” or “current drinking”, respectively, due to small cell counts. Peer drug use was assessed by child-report at age 15. A binary variable was created indicating if friends or anyone in the same class as the young person used cannabis, glue, solvents, or cocaine.
Comorbidity: hyperactivity problems (age 11) and emotional problems (age 15)
In secondary analyses, hyperactivity was considered as an additional baseline confounder based on evidence for associations with conduct problems 1,3,4,20, the potential snares 6,7,19,20, and substance use 19,21–23. Hyperactivity problems were assessed by maternal-report at child age 11 years using the five-item hyperactivity problems subscale from the Strengths and Difficulties Questionnaire (SDQ) 24. The five items in the questionnaire include ‘restless, overactive, cannot stay still for long’, ‘constantly fidgeting or squirming’, ‘easily distracted, concentration wanders’, ‘thinks things out before acting’ and ‘sees tasks through to the end, good attention span’. Responses for each item are ‘Not true’, ‘Somewhat true’ and ‘Certainly true’. The items ‘thinks things out’ and ‘sees tasks through’ were reverse coded. All items were summed together (range 0 to 10).
In secondary analyses, emotional problems were considered as an additional intermediate confounder given that they may be a precursor for the potential snares and substance use 23,25–28, and are also causally downstream from childhood conduct problems (our exposure) 19,29. Emotional problems were assessed by maternal-report at age 15 years using the five-item emotional problems subscale from the SDQ. The five items in the questionnaire include ‘often complains of headaches, stomach-aches or sickness’, ‘many worries, often seems worried’, ‘often unhappy, down-hearted or tearful’, ‘nervous or clingy in new situations, easily loses confidence’ and ‘many fears, easily scared’. Responses for each item are ‘Not true’, ‘Somewhat true’ and ‘Certainly true’. All items were summed together (range 0 to 10).
Supplement 2. Further details on the statistical analyses including i) logistic regression models and ii) the counterfactual approach to mediation
Logistic regression models
We conducted unadjusted and adjusted logistic regressions in stages to examine exposure-mediator and mediator-outcome associations. For exposure-mediator regressions, we first adjusted for baseline confounders, followed by baseline confounders including hyperactivity problems. For mediator-outcome regressions, we first adjusted for baseline and intermediate confounders, followed by these confounders including hyperactivity and emotional problems, followed by a regression model simultaneously including all three mediators, and all confounders. All regression models were governed by the same Directed Acyclic Graph, hence for mediator-outcome regressions, the exposure (conduct problems) becomes another confounder. Caution is required when directly comparing across these logistic regression models, each with different sets of confounders, given that the non-collapsibility of the odds ratio with a common binary outcome (22% for hazardous alcohol use and 14% for illicit drug use) means that the difference between an unadjusted and adjusted odds ratio can be due to both confounder adjustment and non-collapsibility 30.
Counterfactual mediation
The counterfactual approach is based on conceptualising ‘potential outcomes’ for each individual [Y(x)] that would have been observed if particular conditions were met (i.e. had the exposure X been set to the value x through some intervention) – regardless of the conditions that were in fact met for each individual 31. In mediation analyses, four assumptions are made with respect to confounding. These include no unmeasured confounders for any of the paths (exposure-outcome, exposure-mediator and mediator-outcome) and no measured or unmeasured confounders for the association between mediator and outcome which lie on the causal pathway from the exposure. In the current study, we assume that the same set of baseline confounders (sex, a score of sociodemographic and health risk factors, maternal depression, the quality of the father-child and mother-child relationship, parental alcohol consumption, and smoking status, parental separation, neighbourhood safety, and hyperactivity problems) confound all paths, and these were assessed before or during the same assessment as the exposure.
The counterfactual approach to mediation using the parametric g-computation formula was performed to relax the assumption of no measured intermediate confounding and to include binary mediators and outcomes. This counterfactual approach simulates the mediator, outcome and intermediate confounders under each hypothetical intervention or “counter to the fact” scenario. Again, mediation models were performed in stages. First we performed unadjusted models. Next, we adjusted for baseline and intermediate confounders, and finally, we adjusted for all confounders including hyperactivity and emotional problems. Unlike the logistic regression models, -gformula-presents mediation effects as marginal odds ratios meaning that non-collapsibility of the odds ratio is no longer a consideration. We did not hypothesise an exposure-mediator interaction and therefore did not include it in our models. The proportion of the total causal effect that was mediated was calculated by dividing the log odds ratio for the natural indirect effect by the log odds ratio for the total causal effect, and multiplying by 100. We used a Monte Carlo sample size of 10,000 to minimise fluctuations in effect estimates, and a seed of 79.
The total causal effect (TCE) is the value the outcome would take if everybody had been exposed to high conduct problems versus everyone having low conduct problems. The natural direct effect (NDE) is the direct (unmediated) effect of the exposure on the outcome when the mediator takes the value it would take in the absence of the exposure. It is thus modelled as the direct effect of exposure X = 1 (high conduct problems) versus exposure X = 0 (low conduct problems) on outcome Y (e.g., hazardous alcohol use) if mediator M (e.g., gang membership) were set to whatever it would be for X = 0. The NIE captures the effect of the exposure on the outcome that operates by changing the mediator. It is thus modelled as the effect on outcome Y (e.g., hazardous alcohol use) if the exposure were fixed at X = 1 and mediator M (e.g., gang membership) were changed from the level it would take if X = 0 to the level it would take if X = 1.
Supplement 3. Details on multiple imputation
To address missing data, we used multivariate imputation by chained equations32 to impute all incomplete analysis variables up to our starting sample (those that provided data on alcohol use and drug use from at least one time point between age 11 and age 22 years; N = 4,599). Despite the univariable association observed between conduct problems and missingness (suggesting that the exposure could be missing-not-at-random), we decided to impute the small amount of missing data on conduct problems within our starting sample (<5% missing) given the strong auxiliary data and wide range of confounders available for our imputation model, making the missing-at-random assumption plausible. First, we investigated the patterns of missing data, the percentage of missingness for each analysis variable and potential auxiliary variables and determined the distribution of missing variables (see Supplementary Table 2). Next, we performed multivariable regression analyses with each incomplete analysis variable as the outcome and potential auxiliary variables (such as earlier assessments of the mediators and outcomes) as the exposures to determine which auxiliary variables should be used to impute which analysis variables. In the imputation model, all analysis variables (exposure, mediators, outcomes, and confounders) were included in the imputation equation for every imputed variable. Additionally, we used tailored imputation equations to include only the auxiliary variables which were associated with the variable to be imputed in the regression models described above.
We used the Stata command -mi impute chained-with 10 cycles of regression switching and generated 40 imputed datasets. We imputed binary variables using binary logistic regression models, ordinal variables using ordinal logistic regression models, and all continuous variables using predicitve mean matching (with 10 nearest-neighbours) given evidence for a skewed distribution (Supplementary Table 2). Given that gang membership was a composite of two time points, gang membership at age 18 and at age 22 were imputed seperately and then combined to passively impute the composite measure of gang membership at either age, which was then included in the imputation equations for all analysis variables. After data were imputed, we compared the summary statistics (mean or proportion) for each analysis variable across complete case and imputed data to identify if there were any large differences. Convergence plots were also used to ensure that 10 cycles of regression switching was sufficient, and Monte Carlo errors were examined in the regression models for the outcomes and mediators to ensure that 40 imputed datasets were sufficient32.
After imputation, all analyses were combined across the 40 datasets using Rubin’s rules. For descriptive statistics and regression models, this was done using the Stata command -mi estimate-. Mediation models were estimated across each imputed dataset using a loop, with the results from each imputed dataset saved into a logfile. The log odds ratios and bootstrapped standard errors for the mediation effects (total causal effect, natural indirect effect, and natural direct effect) across each imputed dataset were extracted from the logfile. The mean of the log odds ratio for each mediation effect across the 40 imputed datasets was calculated, and the standard error was calculated using Rubin’s rules. Imputation and subsequent analysis code for the mediation models can be found on GitHub: https://github.com/gemmahammerton/gformula_1993_Pelotas
Acknowledgements
This article is based on data from the study "Pelotas Birth Cohort, 1993" conducted by Postgraduate Program in Epidemiology at Universidade Federal de Pelotas with the collaboration of the Brazilian Public Health Association (ABRASCO). From 2004 to 2013, the Wellcome Trust supported the 1993 birth cohort study. The European Union, National Support Program for Centers of Excellence (PRONEX), the Brazilian National Research Council (CNPq), and the Brazilian Ministry of Health supported previous phases of the study. The 22-year follow-up was supported by the Science and Technology Department / Brazilian Ministry of Health, with resources transferred through the Brazilian National Council for Scientific and Technological Development (CNPq), grant 400943/2013-1.
Dr. Hammerton and Professor Hickman would like to acknowledge support from the National Institute for Health Research (NIHR) Health Protection Research Unit in Behavioural Science and Evaluation (HPRU BSE) at the University of Bristol. The views expressed are those of the authors and not necessarily those of the NIHR, the Department of Health and Social Care or UK Health Security Agency (UKHSA).
All authors have reported no biomedical financial interests or potential conflicts of interest.