Study population and follow-up
The Cape Cod Health Study was approved by the Institutional Review Boards of the Massachusetts Department of Public Health and Boston University Medical Center. Eligible individuals were born during 1969-1983 to married women living in Cape Cod, Massachusetts towns known to have some VL/AC water pipes in their water distribution system. The details of study selection have been previously described (27). In brief, two groups of “index” subjects were selected: (1) those born to women who were exposed to PCE-contaminated drinking water at birth (N=1,910), and (2) those born to women who were unexposed at their birth (N=1,928). These groups were identified by cross-matching maternal addresses on birth certificates with information from local water companies on the location and installation year of the VL/AC pipes. Unexposed subjects were randomly selected and frequency matched to exposed subjects on month and year of birth. In addition, older siblings of exposed and unexposed subjects who were born in Massachusetts during 1969-1983 (N=1,202) were identified. These siblings were initially considered unexposed because they were born before the family moved to an affected residence, but final exposure designations for all subjects were assigned after more extensive exposure assessments were conducted, as described below.
Birth certificates of all individuals were reviewed to obtain information on the study family, including the names of the subject and parents; the subject’s date of birth, birth weight and gestational duration; and the parents’ ages and educational levels when the subject was born. During Phase 1 of the study (2000-2005), individuals were traced and sent invitation letters describing the purpose of the study and requesting that they complete a self-administered questionnaire. Overall 40.5% of successfully located subjects returned the study questionnaire (N=1,689). The Phase 1 questionnaire gathered information on cigarette smoking, alcoholic beverage consumption, and drug use, demographic and medical characteristics; and occupational and non-occupational exposure to solvents. Information was also collected on the family’s residences from the subject’s birth through 1990, including the exact street address and calendar years of residence for all Cape Cod addresses. Lastly, the survey collected information on the subject’s knowledge of the PCE contamination episode and knowledge of their own PCE exposure.
Phase 2 of the study (2017-2020) focused on the 1,512 subjects from Phase 1 who had adequate residential data for assessing their PCE exposure status. A total of 694 subjects completed the entire Phase 2 survey (46.8%). The remainder never responded to several contact attempts (N=718), were not found (N=54), were deceased (N=6), declined to participate (N=26), or completed out only a small portion of the survey (N=14). The Phase 2 questionnaire updated information on cigarette smoking, alcoholic beverage consumption and drug use; and demographic and medical characteristics. Drug information included use of cannabis, cocaine, heroin, hallucinogens, inhalants, and methamphetamine, and misuse of prescription pain relievers, tranquilizers, stimulants and sedatives. The Phase 2 questionnaire also collected new information on the eleven established Diagnostic and Statistical Manual of Mental Disorders Fifth Edition (DSM5) criteria for drug use disorder for any drug. Questions for these criteria, which were adapted from the validated Alcohol Use Disorder and Associated Disabilities Interview Schedule (AUDADIS) (40), ascertained the lifetime prevalence of the following behaviors: ever engaging in hazardous use that increased the chances of getting hurt, ever having social/interpersonal problems due to use, often neglecting major roles due to use, ever developing tolerance such that the usual drug amount had much less of an effect than it once did, often using larger amounts or for a longer period of time than intended, ever making unsuccessful attempts to stop or cut down use, ever spending much time getting or using the drug, ever continuing to use a drug despite having physical and psychological problems related to use, ever giving up activities due to use, ever developing cravings for the drug, and ever having withdrawal symptoms when the drug effects were wearing off. We were unable to obtain physician confirmation of these criteria.
Comparison of Phase 2 participants (N=694) and non-participants (N=818) found that they were quite similar with regard to many characteristics, including PCE exposure status, age, race, maternal age, birth weight, gestational age, and receipt of prenatal care (Table 1). However, compared to Phase 2 non-participants, a larger proportion of Phase 2 participants were female, had mothers who were college graduates and fathers employed in white collar jobs, and a smaller proportion smoked cigarettes regularly, reported heavy episodic drinking and used drugs including crack/cocaine, heroin, psychedelics/hallucinogens, club/designer drugs, inhalants, and Ritalin without a prescription (Table 1).
Table 1 Distribution of Selected Characteristics by Phase 2 Participation Status
|
Characteristic
|
Participant
(N=694)
|
Non-Participant
(N=818)
|
PCE exposure status
|
|
|
Both pre and postnatal exposure
|
363 (52.3%)
|
468 (57.2%)
|
Only postnatal exposure
|
76 (11.0%)
|
58 (7.1%)
|
Unexposed
|
255 (36.7%)
|
292 (35.7%)
|
Age (Phase 2, mean, sd)
|
41 (4)
|
40 (4)
|
White race
|
688 (99.1%)
|
802 (98.0%)
|
Sex
|
|
|
Male
|
242 (34.9%)
|
360 (44.0%)
|
Female
|
452 (65.1%)
|
458 (56.0%)
|
Ever smoked on a regular basis (Phase 1 data)
|
202 (29.4%)
|
343 (42.3%)
|
Missing
|
6
|
8
|
Drank 5+(men)/4+(women) drinks* per day, past 30 days (Phase 1 data)
|
80 (11.8%)
|
147 (18.7%)
|
Missing
|
16
|
32
|
Any drug** use (Phase 1 data)
|
230 (33.1%)
|
351 (42.9%)
|
Mother’s age at subject’s birth (Mean, SD)
|
27 (4)
|
27 (5)
|
Mother’s educational level at subject’s birth
|
|
|
High school graduate or less
|
230 (33.2%)
|
340 (41.7%)
|
Some college
|
210 (30.3%)
|
254 (31.1%)
|
4 year college graduate or higher
|
253 (36.5%)
|
222 (27.2%)
|
Missing
|
1
|
2
|
Father’s occupation at subject’s birth
|
|
|
White collar
|
371 (53.9%)
|
370 (46.0%)
|
Blue collar
|
211 (30.7%)
|
285 (35.4%)
|
Other
|
106 (15.4%)
|
149 (18.5%)
|
Missing
|
6
|
14
|
Subject’s birth weight (Mean, SD)
|
3423 (524)
|
3438 (507)
|
Missing
|
60
|
44
|
Subject’s gestational age (weeks) (Mean, SD)
|
40 (2)
|
40 (2)
|
Missing
|
29
|
60
|
Mother received prenatal care during subject’s pregnancy
|
663 (99.8%)
|
764 (99.6%)
|
Missing
|
30
|
51
|
*A drink was defined as 12-ounce bottle, can, or glass of beer, 4-ounce glass of wine, 12-ounce bottle of wine cooler, hard lemonade, or hard cider, shot of liquor straight or in a mixed drink
*Included crack/cocaine, heroin, psychedelics/hallucinogens, club/designer drugs, inhalants, and Ritalin without a prescription
PCE exposure assessment
PCE exposure assessments were conducted during Phase 1. As a first step, all residential addresses on Cape Cod reported by subjects were geocoded to a latitude and longitude using ArcGIS 8.1. Approximately 95% of reported addresses were successfully geocoded without knowledge of the exposure or outcome status. Visual inspection of the water distribution maps determined initial exposure status and was followed by a more detailed assessment using a leaching and transport model developed by Webler and Brown (41). The model estimates the amount of PCE entering a residence by using the initial amount of PCE in the liner, the age of the pipe, and the leaching rate of PCE from the liner into the water. This rate is modeled as an exponential relationship with a rate constant of 2.25 years based on experimental data (42). Additional model parameters include the water flow rate and direction which were determined using EPANET, water distribution system modeling software developed by the U.S. EPA. This software has been used in other epidemiological studies of drinking water contaminants (43-44). We were able to calculate only annual PCE exposures because only move-in and pipe installation years were available. We estimated PCE exposure during the prenatal period by multiplying the annual mass of PCE that entered the subject’s residence during their birth year by 9/12. We estimated cumulative exposure during early childhood by summing the estimated mass of PCE that entered their residences from the month and year following birth through the month and year of the fifth birthday. Simple percentages were used to account for partial years.
Statistical analysis
We compared the criteria for drug use disorder among subjects with prenatal and early childhood exposure combined (N=363) to unexposed subjects (N=255). Nearly all subjects with prenatal exposure also had childhood exposure and so we were unable to examine the impact of prenatal exposure alone. In addition, the number of subjects exposed only during childhood was too small to provide stable results (N=76). In addition, we examined early life PCE exposure in relation to the lifetime presence of any criteria (>=1), two or more criteria, and finer groupings of the number of criteria (1, 2-3, 4-5, and 6 or more) overall and stratified by sex. The lifetime presence of two or more criteria was examined because the formal diagnosis of a drug use disorder requires the presence of two or more criteria, albeit in a single year. We also examined specific criteria (e.g., neglected major roles due to use, physical and psychological problems related to use) and the impact of PCE exposure levels (>=median and >0-<median) to determine if a dose-response relationship was present. Lastly, we compared the drug use history among individuals with and without any criteria for drug use disorder.
The risk ratio (RR) was used to estimate the strength of the association between PCE exposure and the criteria. Ninety-five percent confidence intervals were used to assess their precision. First, crude analyses were conducted and then generalized estimating equation (GEE) analyses were performed to account for non-independent outcomes arising from several children from the same family (45-46). Twelve percent of Phase 2 subjects were siblings. Lastly, adjusted GEE analyses were conducted to control for confounding variables. Covariates considered for these analyses were factors associated with drug use disorder that, according to our directed acyclic graph (DAG), also influenced their PCE exposure. While many variables were risk factors for having criteria for drug use disorder, only those with plausible causal relationships with PCE exposure were ultimately controlled. These variables, which reflected calendar time and socioeconomic status, were the subject’s age, mother’s educational level and father’s occupation.