Heat-related fatalities in North Carolina 1999-2017

OBJECTIVES

Research shows the highest rates of occupational heat-related fatalities among farm laborers and among Black and Hispanic workers in North Carolina (NC). The Hispanic population and workforce in NC have grown substantially in the past 20 years. We describe the epidemiology of heat-related fatal injuries in the general population and among workers in NC.

METHODS

We reviewed North Carolina death records and records of the North Carolina Office of the Chief Medical Examiner to identify heat-related deaths (primary International Classification of Diseases, Tenth Revision diagnosis code: X30 or T67.0-T67.9) that occurred between January 1, 1999, and December 31, 2017. Decedent age, sex, race, and ethnicity were extracted from both the death certificate and the medical examiner's report as well as determinations of whether the death occurred at work.

RESULTS

In NC between 1999 and 2017, there were 225 deaths from heat-related injuries, and 25 occurred at work. The rates of occupational heat-related deaths were highest among males, workers of Hispanic ethnicity, workers of Black, multiple, or unknown race, and in workers aged 55-64. The highest rate of occupational heat-related deaths occurred in the agricultural industry.

CONCLUSIONS

Since the last report (2001), the number of heat-related fatalities has increased, but fewer were identified as workplace fatalities. Rates of occupational heat-related deaths are highest among Hispanic workers. NC residents identifying as Black are disproportionately burdened by heat-related fatalities in general, with a wider apparent disparity in occupational deaths.

Forty years of struggle in North Carolina: Workplace segregation and fatal occupational injury rates

OBJECTIVE

To assess workplace segregation in fatal occupational injury from 1992 to 2017 in North Carolina.

METHODS

We calculated occupational fatal injury rates within categories of occupation, industry, race, age, and sex; and estimated expected numbers of fatalities among Black and Hispanic male workers had they experienced the rates of White male workers. We also estimated the contribution of workforce segregation to disparities by estimating the expected number of fatalities among Black and Hispanic male workers had they experienced the industry and occupation patterns of White male workers. We assessed person-years of life-lost, using North Carolina life expectancy estimates.

RESULTS

Hispanic workers contributed 32% of their worker-years and experienced 58% of their fatalities in construction. Black workers were most overrepresented in the food manufacturing industry. Hispanic males experienced 2.11 (95% CI: 1.86-2.40) times the mortality rate of White males. The Black-White and Hispanic-White disparities were widest among workers aged 45 and older, and segregation into more dangerous industries and occupations played a substantial role in driving disparities. Hispanic workers who suffered occupational fatalities lost a median 47 life-years, compared to 37 among Black workers and 36 among White workers.

CONCLUSIONS

If Hispanic and Black workers experienced the workplace safety of their White counterparts, fatal injury rates would be substantially reduced. Workforce segregation reflects structural racism, which also contributes to mortality disparities. Root causes must be addressed to eliminate disparities.

Single proxy control

Negative control variables are sometimes used in nonexperimental studies to detect the presence of confounding by hidden factors. A negative control outcome (NCO) is an outcome that is influenced by unobserved confounders of the exposure effects on the outcome in view, but is not causally impacted by the exposure. Tchetgen Tchetgen (2013) introduced the Control Outcome Calibration Approach (COCA) as a formal NCO counterfactual method to detect and correct for residual confounding bias. For identification, COCA treats the NCO as an error-prone proxy of the treatment-free counterfactual outcome of interest, and involves regressing the NCO on the treatment-free counterfactual, together with a rank-preserving structural model, which assumes a constant individual-level causal effect. In this work, we establish nonparametric COCA identification for the average causal effect for the treated, without requiring rank-preservation, therefore accommodating unrestricted effect heterogeneity across units. This nonparametric identification result has important practical implications, as it provides single-proxy confounding control, in contrast to recently proposed proximal causal inference, which relies for identification on a pair of confounding proxies. For COCA estimation we propose 3 separate strategies: (i) an extended propensity score approach, (ii) an outcome bridge function approach, and (iii) a doubly-robust approach. Finally, we illustrate the proposed methods in an application evaluating the causal impact of a Zika virus outbreak on birth rate in Brazil.

Lifetime excess absolute risk for lung cancer due to exposure to radon: results of the pooled uranium miners cohort study PUMA

The Pooled Uranium Miners Analysis (PUMA) study is the largest uranium miners cohort with 119,709 miners, 4.3 million person-years at risk and 7754 lung cancer deaths. Excess relative rate (ERR) estimates for lung cancer mortality per unit of cumulative exposure to radon progeny in working level months (WLM) based on the PUMA study have been reported. The ERR/WLM was modified by attained age, time since exposure or age at exposure, and exposure rate. This pattern was found for the full PUMA cohort and the 1960 + sub-cohort, i.e., miners hired in 1960 or later with chronic low radon exposures and exposure rates. The aim of the present paper is to calculate the lifetime excess absolute risk (LEAR) of lung cancer mortality per WLM using the PUMA risk models, as well as risk models derived in previously published smaller uranium miner studies, some of which are included in PUMA. The same methods were applied for all risk models, i.e., relative risk projection up to <95 years of age, an exposure scenario of 2 WLM per year from age 18-64 years, and baseline mortality rates representing a mixed Euro-American-Asian population. Depending upon the choice of model, the estimated LEAR per WLM are 5.38 × 10-4 or 5.57 × 10-4 in the full PUMA cohort and 7.50 × 10-4 or 7.66 × 10-4 in the PUMA 1960 + sub-cohort, respectively. The LEAR per WLM estimates derived from risk models reported for previously published uranium miners studies range from 2.5 × 10-4 to 9.2 × 10-4. PUMA strengthens knowledge on the radon-related lung cancer LEAR, a useful way to translate models for policy purposes.

25-Year fatal workplace suicide trends in North Carolina: 1992-2017

BACKGROUND

Suicide is a serious public health problem in the United States, but limited evidence is available investigating fatal suicides at work. There is a substantial need to characterize workplace suicides to inform suicide prevention interventions and target high-risk settings. This study aims to examine workplace suicide rates in North Carolina (NC) by worker characteristics, means of suicide used, and industry between 1992 and 2017.

METHODS

Fatal workplace suicides were identified from records of the NC Office of the Chief Medical Examiner system and the NC death certificate. Sex, age, race, ethnicity, class of worker, manner of death, and industry were abstracted. Crude and age-standardized homicide rates were calculated as the number of suicides that occurred at work divided by an estimate of worker-years (w-y). Rate ratios and 95% confidence intervals (CIs) were calculated, and trends over calendar time for fatal workplace suicides were examined overall and by industry.

RESULTS

81 suicides over 109,464,430 w-y were observed. Increased rates were observed in workers who were male, self-employed, and 65+ years old. Firearms were the most common means of death (63%) followed by hanging (16%). Gas service station workers experienced the highest fatal occupational suicide rate, 11.5 times (95% CI: 3.62-36.33) the overall fatal workplace suicide rate, followed by Justice, Public Order, and Safety workers at 3.23 times the overall rate (95% CI: 1.31-7.97).

CONCLUSION

Our findings identify industries and worker demographics that were vulnerable to workplace suicides. Targeted and tailored mitigation strategies for vulnerable industries and workers are recommended.

Forty-year trends in fatal occupational injuries in North Carolina

BACKGROUND

We describe progress in the control of deaths on-the-job due to fatal occupational injury in North Carolina over the period 1978-2017.

METHODS

Forty years of information on fatal occupational injuries in North Carolina has been assembled from medical examiners' reports and death certificates, supplemented by newspaper and police reports. Cases were defined as unintentional fatal occupational injuries among adults. Annual estimates of the population at risk were derived from US Census data, and rates were quantified using Poisson regression methods.

RESULTS

There were 4434 eligible deaths. The unintentional fatal occupational injury rate at the beginning of the study period was more than threefold the rate at the end of the study. The fatal occupational injury rate among men declined from 9.6 per 100,000 worker-years in the period 1978-1982 to 3.1 per 100,000 worker-years in the period 2013-2017. The fatal occupational injury rate among women declined from 0.3 per 100,000 worker-years in the period 1978-1981 to 0.1 per 100,000 worker-years in the period 2013-2017. Declines in rates were observed for young adults as well as older workers and were observed across all major industry categories. Average annual declines in rates were greatest in those industries and occupations that had the highest fatal injury rates at the start of the study period.

CONCLUSIONS

The substantial decline in fatal injury rates underscores the importance of injury prevention and demonstrates the ability to make meaningful reductions in unintentional fatal injury.

Universal Difference-in-Differences for Causal Inference in Epidemiology

Difference-in-differences is undoubtedly one of the most widely used methods for evaluating the causal effect of an intervention in observational (i.e., nonrandomized) settings. The approach is typically used when pre- and postexposure outcome measurements are available, and one can reasonably assume that the association of the unobserved confounder with the outcome has the same absolute magnitude in the two exposure arms and is constant over time; a so-called parallel trends assumption. The parallel trends assumption may not be credible in many practical settings, for example, if the outcome is binary, a count, or polytomous, as well as when an uncontrolled confounder exhibits nonadditive effects on the distribution of the outcome, even if such effects are constant over time. We introduce an alternative approach that replaces the parallel trends assumption with an odds ratio equi-confounding assumption under which an association between treatment and the potential outcome under no treatment is identified with a well-specified generalized linear model relating the pre-exposure outcome and the exposure. Because the proposed method identifies any causal effect that is conceivably identified in the absence of confounding bias, including nonlinear effects such as quantile treatment effects, the approach is aptly called universal difference-in-differences. We describe and illustrate both fully parametric and more robust semiparametric universal difference-in-differences estimators in a real-world application concerning the causal effects of a Zika virus outbreak on birth rate in Brazil. A supplementary digital video is available at: http://links.lww.com/EDE/C90.

Higher-order evidence

Higher-order evidence is evidence about evidence. Epidemiologic examples of higher-order evidence include the settings where the study data constitute first-order evidence and estimates of misclassification comprise the second-order evidence (e.g., sensitivity, specificity) of a binary exposure or outcome collected in the main study. While sampling variability in higher-order evidence is typically acknowledged, higher-order evidence is often assumed to be free of measurement error (e.g., gold standard measures). Here we provide two examples, each with multiple scenarios where second-order evidence is imperfectly measured, and this measurement error can either amplify or attenuate standard corrections to first-order evidence. We propose a way to account for such imperfections that requires third-order evidence. Further illustrations and exploration of how higher-order evidence impacts results of epidemiologic studies is warranted.

Fatal occupational injuries in North Carolina, 1992-2017

OBJECTIVES

After declining for several decades, fatal occupational injury rates have stagnated in the USA since 2009. To revive advancements in workplace safety, interventions targeting at-risk worker groups must be implemented. Our study aims to identify these at-risk populations by evaluating disparities in unintentional occupational fatalities occurring in North Carolina (NC) from 1992 to 2017.

METHODS

Our retrospective cohort study drew on both the NC Office of the Chief Medical Examiner system and the NC death certificate data system to identify unintentional fatal occupational injuries occurring from 1992 to 2017. Unintentional fatal occupational injury rates were reported across industries, occupations and demographic groups, and rate ratios were calculated to assess disparities.

RESULTS

Among those aged 18 and older, 2645 unintentional fatal occupational injuries were identified. Fatal occupational injury rates declined by 0.82 injuries/100 000 person-years over this period, falling consistently from 2004 to 2009 and increasing from 2009 to 2017. Fatal injury rates were highest among Hispanic workers, who experienced 2.75 times the fatal injury rate of non-Hispanic White workers (95% CI 2.42 to 3.11) and self-employed workers, who experienced 1.44 times the fatal injury rate of private workers (95% CI 1.29 to 1.60). We also observed that fatal injury rates increased with age group and were higher among male relative to female workers even after adjustment for differential distributions across occupations.

CONCLUSIONS

The decline in unintentional fatal occupational injury rates over this period is encouraging, but the increase in injury rate after 2009 and the large disparities between occupations, industries and demographic groups highlight the need for additional targeted safety interventions.

Estimating the Effect of a Treatment When There Is Nonadherence in a Trial

Randomized trials offer a powerful strategy for estimating the effect of a treatment on an outcome. However, interpretation of trial results can be complicated when study subjects do not take the treatment to which they were assigned; this is referred to as nonadherence. Prior authors have described instrumental variable approaches to analyze trial data with nonadherence; under their approaches, the initial assignment to treatment is used as an instrument. However, their approaches require the assumption that initial assignment to treatment has no direct effect on the outcome except via the actual treatment received (i.e., the exclusion restriction), which may be implausible. We propose an approach to identification of a causal effect of treatment in a trial with 1-sided nonadherence without assuming exclusion restriction. The proposed approach leverages the study subjects initially assigned to control status as an unexposed reference population; we then employ a bespoke instrumental variable analysis, where the key assumption is "partial exchangeability" of the association between a covariate and an outcome in the treatment and control arms. We provide a formal description of the conditions for identification of causal effects, illustrate the method using simulations, and provide an empirical application.

Dysfunction of the Salivary and Lacrimal Glands After Radioiodine Therapy for Thyroid Cancer: Results of the START Study After 6-Months of Follow-Up

Background: Understanding of changes in salivary and lacrimal gland functions after radioactive iodine therapy (131I-therapy) remains limited, and, to date, no studies have evaluated dose-response relationships between absorbed dose from 131I-therapy and dysfunctions of these glands. This study investigates salivary/lacrimal dysfunctions in differentiated thyroid cancer (DTC) patients six months after 131I-therapy, identifies 131I-therapy-related risk factors for salivary/lacrimal dysfunctions, and assesses the relationships between 131I-therapy radiation dose and these dysfunctions. Methods: A cohort study was conducted involving 136 DTC patients treated by 131I-therapy of whom 44 and 92 patients received 1.1 and 3.7 GBq, respectively. Absorbed dose to the salivary glands was estimated using a dosimetric reconstruction method based on thermoluminescent dosimeter measurements. Salivary and lacrimal functions were assessed at baseline (T0, i.e., immediately before 131I-therapy) and six months later (T6) using validated questionnaires and salivary samplings, with and without stimulation of the salivary glands. Statistical analyses included descriptive analyses and random-effects multivariate logistic and linear regressions. Results: There was no difference between T0 and T6 in the level of parotid gland pain, nor was there difference in the number of patients with hyposalivation, but there were significantly more patients with dry mouth sensation and dry eyes after therapy compared with baseline. Age, menopause, depression and anxiety symptoms, history of systemic disease, and not taking painkillers in the past three months were found to be significantly associated with salivary or lacrimal disorders. Significant associations were found between 131I-exposure and salivary disorders adjusted on the previous variables: for example, per 1-Gy increase in mean dose to the salivary glands, odds ratio = 1.43 [CI 1.02 to 2.04] for dry mouth sensation, ß = -0.08 [CI -0.12 to -0.02] mL/min for stimulated saliva flow, and ß = 1.07 [CI 0.42 to 1.71] mmol/L for salivary potassium concentration. Conclusions: This study brings new knowledge on the relationship between the absorbed dose to the salivary glands from 131I-therapy and salivary/lacrimal dysfunctions in DTC patients six months after 131I-therapy. Despite the findings of some dysfunctions, the results do not show any obvious clinical disorders after the 131I-therapy. Nevertheless, this study raises awareness of the risk factors for salivary disorders, and calls for longer follow-up. Clinical Trials Registration: Number NCT04876287 on the public website (ClinicalTrials.gov).

Cancer mortality after low dose exposure to ionising radiation in workers in France, the United Kingdom, and the United States (INWORKS): cohort study

OBJECTIVE

To evaluate the effect of protracted low dose, low dose rate exposure to ionising radiation on the risk of cancer.

DESIGN

Multinational cohort study.

SETTING

Cohorts of workers in the nuclear industry in France, the UK, and the US included in a major update to the International Nuclear Workers Study (INWORKS).

PARTICIPANTS

309 932 workers with individual monitoring data for external exposure to ionising radiation and a total follow-up of 10.7 million person years.

MAIN OUTCOME MEASURES

Estimates of excess relative rate per gray (Gy) of radiation dose for mortality from cancer.

RESULTS

The study included 103 553 deaths, of which 28 089 were due to solid cancers. The estimated rate of mortality due to solid cancer increased with cumulative dose by 52% (90% confidence interval 27% to 77%) per Gy, lagged by 10 years. Restricting the analysis to the low cumulative dose range (0-100 mGy) approximately doubled the estimate of association (and increased the width of its confidence interval), as did restricting the analysis to workers hired in the more recent years of operations when estimates of occupational external penetrating radiation dose were recorded more accurately. Exclusion of deaths from lung cancer and pleural cancer had a modest effect on the estimated magnitude of association, providing indirect evidence that the association was not substantially confounded by smoking or occupational exposure to asbestos.

CONCLUSIONS

This major update to INWORKS provides a direct estimate of the association between protracted low dose exposure to ionising radiation and solid cancer mortality based on some of the world's most informative cohorts of radiation workers. The summary estimate of excess relative rate solid cancer mortality per Gy is larger than estimates currently informing radiation protection, and some evidence suggests a steeper slope for the dose-response association in the low dose range than over the full dose range. These results can help to strengthen radiation protection, especially for low dose exposures that are of primary interest in contemporary medical, occupational, and environmental settings.

Ionizing radiation and solid cancer mortality among US nuclear facility workers

BACKGROUND

The risk of solid cancers from low-level protracted ionizing radiation is not well characterized. Nuclear workers provide valuable information on the effects of ionizing radiation in contemporary exposure scenarios relevant to workers and the public.

METHODS

We evaluated the association between penetrating ionizing radiation exposure and solid cancer mortality among a pooled cohort of nuclear workers in the USA, with extended follow-up to examine cancers with long latencies. This analysis includes 101 363 workers from five nuclear facilities, with 12 069 solid cancer deaths between 1944 and 2016. The association between cumulative equivalent dose measured in sieverts (Sv) and solid cancer subtypes were modelled as the excess relative rate per Sv (ERR Sv-1) using Cox regression.

RESULTS

For the association between ionizing radiation exposure and all solid cancer mortality we observed an elevated rate (ERR Sv-1=0.19; 95% CI: -0.10, 0.52), which was higher among a contemporary sub-cohort of workers first hired in 1960 or later (ERR Sv-1= 2.23; 95% CI: 1.13, 3.49). Similarly, we observed an elevated rate for lung cancer mortality (ERR Sv-1= 0.65; 95% CI: 0.09, 1.30) that was higher among contemporary hires (ERR Sv-1= 2.90; 95% CI: 1.00, 5.26).

CONCLUSIONS

Although concerns remain about confounding, measurement error and precision, this analysis strengthens the evidence base indicating there are radiogenic risks for several solid cancer types.

Radon and lung cancer in the pooled uranium miners analysis (PUMA): highly exposed early miners and all miners

OBJECTIVES

Radon is a ubiquitous occupational and environmental lung carcinogen. We aim to quantify the association between radon progeny and lung cancer mortality in the largest and most up-to-date pooled study of uranium miners.

METHODS

The pooled uranium miners analysis combines 7 cohorts of male uranium miners with 7754 lung cancer deaths and 4.3 million person-years of follow-up. Vital status and lung cancer deaths were ascertained between 1946 and 2014. The association between cumulative radon exposure in working level months (WLM) and lung cancer was modelled as the excess relative rate (ERR) per 100 WLM using Poisson regression; variation in the association by temporal and exposure factors was examined. We also examined analyses restricted to miners first hired before 1960 and with <100 WLM cumulative exposure.

RESULTS

In a model that allows for variation by attained age, time since exposure and annual exposure rate, the ERR/100 WLM was 4.68 (95% CI 2.88 to 6.96) among miners who were less than 55 years of age and were exposed in the prior 5 to <15 years at annual exposure rates of <0.5 WL. This association decreased with older attained age, longer time since exposure and higher annual exposure rate. In analyses restricted to men first hired before 1960, we observed similar patterns of association but a slightly lower estimate of the ERR/100 WLM.

CONCLUSIONS

This new large, pooled study confirms and supports a linear exposure-response relationship between cumulative radon exposure and lung cancer mortality which is jointly modified by temporal and exposure factors.

Prenatal organophosphate ester exposure and executive function in Norwegian preschoolers

Organophosphate esters (OPEs) are ubiquitous chemicals, used as flame retardants and plasticizers. OPE usage has increased over time as a substitute for other controlled compounds. This study investigates the impact of prenatal OPE exposure on executive function (EF) in preschoolers.

Methods

We selected 340 preschoolers from the Norwegian Mother, Father, and Child Cohort Study. Diphenyl-phosphate (DPhP), di-n-butyl-phosphate (DnBP), bis(2-butoxyethyl) phosphate (BBOEP), and bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) were measured in maternal urine. EF was measured using the Behavior Rating Inventory of Executive Functioning-Preschool (BRIEF-P) and the Stanford-Binet fifth edition (SB-5). EF scores were scaled so a higher score indicated worse performance. We estimated exposure-outcome associations and evaluated modification by child sex using linear regression.

Results

Higher DnBP was associated with lower EF scores across multiple rater-based domains. Higher DPhP and BDCIPP were associated with lower SB-5 verbal working memory (β = 0.49, 95% CI = 0.12, 0.87; β = 0.53, 95% CI = 0.08, 1.02), and higher BBOEP was associated with lower teacher-rated inhibition (β = 0.34, 95% CI = 0.01, 0.63). DPhP was associated with lower parent-reported BRIEF-P measures in boys but not girls [inhibition: boys: 0.37 (95% CI = 0.03, 0.93); girls: -0.48 (95% CI = -1.27, 0.19); emotional control: boys: 0.44 (95% CI = -0.13, 1.26); girls: -0.83 (95% CI = -1.73, -0.00); working memory: boys: 0.49 (95% CI = 0.03, 1.08); girls: -0.40 (95% CI = -1.11, 0.36)]. Fewer sex interactions were observed for DnBP, BBOEP, and BDCIPP, with irregular patterns observed across EF domains.

Conclusions

We found some evidence prenatal OPE exposure may impact EF in preschoolers and variation in associations by sex.

Disparities in Fatal Occupational Injury Rates in North Carolina, 1978-2017: Comparing Non-Managerial Employees to Managers

BACKGROUND

We examined fatal occupational injuries among private-sector workers in North Carolina during the 40-year period 1978-2017, comparing the occurrence of fatal injuries among non-managerial employees to that experienced by managers.

METHODS

We estimated a standardized fatal occupational injury ratio by inverse probability of exposure weighting, taking the non-managerial workers as the target population. When this ratio measure takes a value greater than unity it signals settings in which non-managerial employees are not provided as safe a work environment as that provided for managers.

RESULTS

Across all industries, non-managerial workers in North Carolina experienced fatal occupational injury rates 8.2 (95%CI: 7.0, 10.0) times the rate experienced by managers. Disparities in fatal injury rates between managers and the employees they supervise were greatest in forestry, rubber, and metal manufacturing, wholesale trade, fishing and extractive industries, and construction.

CONCLUSIONS

The results may help focus discussion about workplace safety between labor and management upon equity, with a goal of providing a work environment for non-managerial employees as safe as the one provided for managers.

Inverse Probability Weighting to Estimate Exposure Effects on the Burden of Recurrent Outcomes in the Presence of Competing Events

Recurrent events-outcomes that an individual can experience repeatedly over the course of follow-up-are common in epidemiologic and health services research. Studies involving recurrent events often focus on time to first occurrence or on event rates, which assume constant hazards over time. In this paper, we contextualize recurrent event parameters of interest using counterfactual theory in a causal inference framework and describe an approach for estimating a target parameter referred to as the mean cumulative count. This approach leverages inverse probability weights to control measured confounding with an existing (and underutilized) nonparametric estimator of recurrent event burden first proposed by Dong et al. in 2015. We use simulations to demonstrate the unbiased estimation of the mean cumulative count using the weighted Dong-Yasui estimator in a variety of scenarios. The weighted Dong-Yasui estimator for the mean cumulative count allows researchers to use observational data to flexibly estimate and contrast the expected number of cumulative events experienced per individual by a given time point under different exposure regimens. We provide code to ease application of this method.

Sulfur dioxide reduction at coal-fired power plants in North Carolina and associations with preterm birth among surrounding residents

Coal-fired power plants (CFPP) are major contributors of air pollution, including the majority of anthropogenic sulfur dioxide (SO2) emissions, which have been associated with preterm birth (PTB). To address a 2002 North Carolina (NC) policy, 14 of the largest NC CFPPs either installed desulfurization equipment (scrubbers) or retired coal units, resulting in substantial reductions of SO2 air emissions. We investigated whether SO2 air emission reduction strategies at CFPPs in NC were associated with changes in prevalence of PTB in nearby communities.

Methods

We used US EPA Air Markets Program Data to track SO2 emissions and determine the implementation dates of intervention at CFPPs and geocoded 2003-2015 NC singleton live births. We conducted a difference-in-difference analysis to estimate change in PTB associated with change in SO2 reduction strategies for populations living 0-<4 and 4-<10 miles from CFPPs pre- and postintervention, with a comparison of those living 10-<15 miles from CFPPs.

Results

With the spatial-temporal exposure restrictions applied, 42,231 and 41,218 births were within 15 miles of CFPP-scrubbers and CFPP-retired groups, respectively. For residents within 4-<10 miles from a CFPP, we estimated that the absolute prevalence of PTB decreased by -1.5% [95% confidence interval (CI): -2.6, -0.4] associated with scrubber installation and -0.5% (95% CI: -1.6, 0.6) associated with the retirement of coal units at CFPPs. Our findings were imprecise and generally null-to-positive among those living within 0-<4 miles regardless of the intervention type.

Conclusions

Results suggest a reduction of PTB among residents 4-<10 miles of the CFPPs that installed scrubbers.

The carcinogenicity of opium consumption: a systematic review and meta-analysis

The carcinogenicity of opium consumption was recently evaluated by a Working Group convened by the International Agency for Research on Cancer (IARC). We supplement the recent IARC evaluation by conducting an extended systematic review as well as a quantitative meta-analytic assessment of the role of opium consumption and risk for selected cancers, evaluating in detail various aspects of study quality on meta-analytic findings. We searched the published literature to identify all relevant studies on opium consumption and risk of selected cancers in humans through 31 October, 2022. Meta-relative risks (mRRs) and associated 95% confidence intervals (CIs) were estimated using random-effects models for studies of cancer of the urinary bladder, larynx, lung, oesophagus, pancreas, and stomach. Heterogeneity among studies was assessed using the I2 statistic. We assessed study quality and conducted sensitivity analyses to evaluate the impact of potential reverse causation, protopathic bias, selection bias, information bias, and confounding. In total, 2 prospective cohort studies and 33 case-control studies were included. The overall pooled mRR estimated for 'ever or regular' versus 'never' use of opium ranged from 1.50 (95% CI 1.13-1.99, I2 = 0%, 6 studies) for oesophageal cancer to 7.97 (95% CI 4.79-13.3, I2 = 62%, 7 studies) for laryngeal cancer. Analyses of cumulative opium exposure suggested greater risk of cancer associated with higher opium consumption. Findings were robust in sensitivity analyses excluding studies prone to potential methodological sources of biases and confounding. Findings support an adverse association between opium consumption and cancers of the urinary bladder, larynx, lung, oesophagus, pancreas and stomach.

Ionising radiation and cardiovascular disease: systematic review and meta-analysis

OBJECTIVE

To systematically review and perform a meta-analysis of radiation associated risks of cardiovascular disease in all groups exposed to radiation with individual radiation dose estimates.

DESIGN

Systematic review and meta-analysis.

MAIN OUTCOME MEASURES

Excess relative risk per unit dose (Gy), estimated by restricted maximum likelihood methods.

DATA SOURCES

PubMed and Medline, Embase, Scopus, Web of Science Core collection databases.

ELIGIBILITY CRITERIA FOR SELECTING STUDIES

Databases were searched on 6 October 2022, with no limits on date of publication or language. Animal studies and studies without an abstract were excluded.

RESULTS

The meta-analysis yielded 93 relevant studies. Relative risk per Gy increased for all cardiovascular disease (excess relative risk per Gy of 0.11 (95% confidence interval 0.08 to 0.14)) and for the four major subtypes of cardiovascular disease (ischaemic heart disease, other heart disease, cerebrovascular disease, all other cardiovascular disease). However, interstudy heterogeneity was noted (P<0.05 for all endpoints except for other heart disease), possibly resulting from interstudy variation in unmeasured confounders or effect modifiers, which is markedly reduced if attention is restricted to higher quality studies or those at moderate doses (<0.5 Gy) or low dose rates (<5 mGy/h). For ischaemic heart disease and all cardiovascular disease, risks were larger per unit dose for lower dose (inverse dose effect) and for fractionated exposures (inverse dose fractionation effect). Population based excess absolute risks are estimated for a number of national populations (Canada, England and Wales, France, Germany, Japan, USA) and range from 2.33% per Gy (95% confidence interval 1.69% to 2.98%) for England and Wales to 3.66% per Gy (2.65% to 4.68%) for Germany, largely reflecting the underlying rates of cardiovascular disease mortality in these populations. Estimated risk of mortality from cardiovascular disease are generally dominated by cerebrovascular disease (around 0.94-1.26% per Gy), with the next largest contribution from ischaemic heart disease (around 0.30-1.20% per Gy).

CONCLUSIONS

Results provide evidence supporting a causal association between radiation exposure and cardiovascular disease at high dose, and to a lesser extent at low dose, with some indications of differences in risk between acute and chronic exposures, which require further investigation. The observed heterogeneity complicates a causal interpretation of these findings, although this heterogeneity is much reduced if only higher quality studies or those at moderate doses or low dose rates are considered. Studies are needed to assess in more detail modifications of radiation effect by lifestyle and medical risk factors.

SYSTEMATIC REVIEW REGISTRATION

PROSPERO CRD42020202036.

Gestational organophosphate ester exposure and preschool attention-deficit/hyperactivity disorder in the Norwegian Mother, Father, and Child cohort study

BACKGROUND

Attention-deficit/hyperactivity-disorder (ADHD) is a leading neurodevelopmental disorder in children worldwide; however, few modifiable risk factors have been identified. Organophosphate esters (OPEs) are ubiquitous chemical compounds that are increasingly prevalent as a replacement for other regulated chemicals. Current research has linked OPEs to neurodevelopmental deficits. The purpose of this study was to assess gestational OPE exposure on clinically-assessed ADHD in children at age 3 years.

METHODS

In this nested case-control study within the Norwegian Mother, Father, and Child Cohort study, we evaluated the impact of OPE exposure at 17 weeks' gestation on preschool-age ADHD. Between 2007 and 2011, 260 ADHD cases were identified using the Preschool Age Psychiatric Assessment and compared to a birth-year-stratified control group of 549 children. We categorized bis(2-butoxyethyl) phosphate (BBOEP) and bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as values < limit of detection (LOD) (BBOEP N = 386, BDCIPP N = 632), ≥LOD but < limit of quantification (LOQ) (BBOEP N = 413; BDCIPP N = 75), or above LOQ (BBOEP N = 70; BDCIPP N = 102). Diphenyl phosphate (DPhP) and di-n-butyl phosphate (DnBP) were categorized as quartiles and also modeled with a log10 linear term. We estimated multivariable adjusted odds ratios (ORs) using logistic regression and examined modification by sex using an augmented product term approach.

RESULTS

Mothers in the 3rd DnBP quartile had 1.71 times the odds of having a child with ADHD compared to the 1st quartile (95%CI: 1.13, 2.58); a similar trend was observed for log10 DnBP and ADHD. Mothers with BDCIPP ≥ LOD but < LOQ had 1.39 times the odds of having a child with ADHD compared to those with BDCIPP < LOD (95%CI: 0.83, 2.31). Girls had lower odds of ADHD with increasing BBOEP exposure (log10 OR: 0.55 (95%CI: 0.37, 0.93), however boys had a weakly increased odds (log10 OR: 1.25 (95%CI: 0.74, 2.11) p-interaction = 0.01].

CONCLUSIONS

We found modest increased odds of preschool ADHD with higher DnBP and BDCIPP exposure.

Trends in fatal occupational injury rates among older workers before and after the Great Recession of 2008

BACKGROUND

Older workers experience higher rates of fatal occupational injury than younger workers worldwide. In North Carolina, the population of older workers more than doubled between 2000 and 2017. In 2008, the Great Recession changed occupational patterns among all age groups. We examined annual rates and distribution of fatal occupational injuries experienced by older workers, comparing the pre-recession period (2000-2007) to the post-recession period (2009-2017).

METHODS

Detailed information on all fatal occupational injuries during the period between 1 January 2000 and 31 December 2017 were abstracted from the records of the North Carolina Office of the Chief Medical Examiner and the office of vital records. The decennial Census and American Community Survey were used to estimate the population at risk and derive annual rates of fatal occupational injury.

RESULTS

During the study period, 537 occupational fatalities occurred among workers 55+ years of age. The rate of fatal occupational injury among older workers declined 2.8% per year, with a 7.7% yearly decline in the pre-recession period compared with a 1.4% increase per year in the post-recession period. Workers 65+ years of age experienced rate increases in both periods. The highest rates of unintentional fatal occupational injury (injuries that were not purposefully inflicted) were observed in forestry, fishing hunting and trapping, and wood building manufacturing. Intentional fatal occupational injury rates (homicide, suicide) were highest in transportation, gas/service stations and grocery/food stores.

CONCLUSIONS

Older workers have persistently high rates of fatal occupational injury in North Carolina before and after the Great Recession.

Generalized Difference-in-Differences

Difference-in-differences (DID) analyses are used in a variety of research areas as a strategy for estimating the causal effect of a policy, program, intervention, or environmental hazard (hereafter, treatment). The approach offers a strategy for estimating the causal effect of a treatment using observational (i.e., nonrandomized) data in which outcomes on each study unit have been measured both before and after treatment. To identify a causal effect, a DID analysis relies on an assumption that confounding of the treatment effect in the pretreatment period is equivalent to confounding of the treatment effect in the post treatment period. We propose an alternative approach that can yield identification of causal effects under different identifying conditions than those usually required for DID. The proposed approach, which we refer to as generalized DID, has the potential to be used in routine policy evaluation across many disciplines, as it essentially combines two popular quasiexperimental designs, leveraging their strengths while relaxing their usual assumptions. We provide a formal description of the conditions for identification of causal effects, illustrate the method using simulations, and provide an empirical example based on Card and Krueger's landmark study of the impact of an increase in minimum wage in New Jersey on employment.

Associations of Air Pollution and Serum Biomarker Abnormalities in Individuals with Hemodialysis-Dependent Kidney Failure

BACKGROUND

Ambient particles with a median aerodynamic diameter of <2.5 µm (PM2.5) is a ubiquitous air pollutant with established adverse health consequences. While postulated to promote a systemic inflammatory response, limited studies have demonstrated changes in serum biomarkers related to PM2.5 exposure. We aim to examine associations between short-term PM2.5 exposure and commonly measured biomarkers known to be affected by inflammation among patients receiving maintenance in-center hemodialysis.

METHODS

We conducted a retrospective open cohort study from January 1, 2008, to December 31, 2014. Adult hemodialysis patients were identified from the United States Renal Data System and linked at the patient level to laboratory data from a large dialysis organization. Daily ambient PM2.5 was estimated on a 1-km grid and assigned to cohort patients based on the ZIP codes of dialysis clinics. Serum albumin, serum ferritin, transferrin saturation (TSAT), and serum hemoglobin were ascertained from the dialysis provider organization database. Mixed-effect models were used to assess the changes in biomarker levels associated with PM2.5 exposure.

RESULTS

The final cohort included 173,697 hemodialysis patients. Overall, the daily ZIP-level ambient PM2.5 averages were 8.4-8.5 µg/m3. A 10-µg/m3 increase in same-day ambient PM2.5 exposure was associated with higher relative risks of lower albumin (relative risk [RR], 1.01; 95% confidence interval [95% CI], 1.01 to 1.02) and lower hemoglobin (RR, 1.02; 95% CI, 1.01 to 1.03). Associations of same-day ambient PM2.5 exposure and higher ferritin and lower TSAT did not reach statistical significance.

CONCLUSIONS

Short-term PM2.5 exposure was associated with lower serum hemoglobin and albumin among patients receiving in-center hemodialysis. These findings lend support to the role of inflammation in PM2.5 exposure-outcome associations.

Association Between Long-term Ambient PM2.5 Exposure and Cardiovascular Outcomes Among US Hemodialysis Patients

RATIONALE & OBJECTIVE

Ambient PM_2.5 (particulate matter with a diameter of 2.5 microns) is a ubiquitous air pollutant with established adverse cardiovascular (CV) effects. However, quantitative estimates of the association between PM_2.5 exposure and CV outcomes in the setting of kidney disease are limited. This study assessed the association of long-term PM_2.5 exposure with CV events and cardiovascular disease (CVD)-specific mortality among patients receiving maintenance in-center hemodialysis (HD).

STUDY DESIGN

Retrospective cohort study.

SETTINGS & PARTICIPANTS

314,079 adult kidney failure patients initiating HD between 2011 and 2016 identified from the US Renal Data System.

EXPOSURE

Estimated daily ZIP code-level PM_2.5 concentrations were used to calculate each participant's annual average PM_2.5 exposure based on the dialysis clinics visited during the 365 days before the outcome.

OUTCOME

CV event and CVD-specific mortality were ascertained based on ICD-9/ICD-10 diagnostic codes and recorded cause of death from Centers for Medicare & Medicaid Services form 2746.

ANALYTICAL APPROACH

Discrete time hazards models were used to estimate hazards ratios per 1 μg/m³ greater annual average PM_2.5, adjusting for temperature, humidity, day of the week, season, age at baseline, race, employment status, and geographic region. Effect measure modification was assessed for age, sex, race, and baseline comorbidities.

RESULTS

Each 1 μg/m³ greater annual average PM_2.5 was associated with a greater rate of CV events (HR, 1.02 [95% CI, 1.01-1.02]) and CVD-specific mortality (HR, 1.02 [95% CI, 1.02-1.03]). The association was more pronounced for people who initiated dialysis at an older age, had chronic obstructive pulmonary disease (COPD) at baseline, or were Asian. Evidence of effect modification was also observed across strata of race, and other baseline comorbidities.

LIMITATIONS

Potential exposure misclassification and unmeasured confounding.

CONCLUSIONS

Long-term ambient PM_2.5 exposure was associated with CVD outcomes among patients receiving maintenance in-center HD. Stronger associations between long-term PM_2.5 exposure and adverse effects were observed among patients who were of advanced age, had COPD, or were Asian.

PLAIN-LANGUAGE SUMMARY

Long-term exposure to air pollution, also called PM_2.5, has been linked to adverse cardiovascular outcomes. However, little is known about the association of PM_2.5 and outcomes among patients receiving dialysis, who are individuals with high cardiovascular disease burdens. We conducted an epidemiological study to assess the association between the annual PM_2.5 exposure and cardiovascular events and death among patients receiving regular outpatient hemodialysis in the United States between 2011 and 2016. We found a higher risk of heart attacks, strokes, and related events in patients exposed to higher levels of air pollution. Stronger associations between air pollution and adverse health events were observed among patients who were older at the start of dialysis, had chronic obstructive pulmonary disease, or were Asian. These findings bolster the evidence base linking air pollution and adverse health outcomes and may inform policy makers and clinicians.

Challenges with misclassification of American Indian/Alaska Native race and Hispanic ethnicity on death records in North Carolina occupational fatalities surveillance

As frequently segregated and exploitative environments, workplaces are important sites in driving health and mortality disparities by race and ethnicity. Because many worksites are federally regulated, US workplaces also offer opportunities for effectively intervening to mitigate these disparities. Development of policies for worker safety and equity should be informed by evidence, including results from research studies that use death records and other sources of administrative data. North Carolina has a long history of Black/white disparities in work-related mortality and evidence of such disparities is emerging in Hispanic and American Indian/Alaska Native (AI/AN) worker populations. The size of Hispanic and AI/AN worker populations have increased in North Carolina over the last decade, and North Carolina has the largest AI/AN population in the eastern US. Previous research indicates that misidentification of Hispanic and AI/AN identities on death records can lead to underestimation of race/ethnicity-specific mortality rates. In this commentary, we describe problems and complexities involved in determining AI/AN and Hispanic identities from North Carolina death records. We provide specific examples of misidentification that are likely introducing bias to occupational mortality disparity documentation, and offer recommendations for improved data collection, analysis, and interpretation. Our primary recommendation is to build and maintain relationships with local community leadership, so that improvements in the ascertainment of race and ethnicity are grounded in the lived experience of workers from communities of color.

A Bespoke Instrumental Variable Approach to Correction for Exposure Measurement Error

A covariate-adjusted estimate of an exposure-outcome association may be biased if the exposure variable suffers measurement error. We propose an approach to correct for exposure measurement error in a covariate-adjusted estimate of the association between a continuous exposure variable and outcome of interest. Our proposed approach requires data for a reference population in which the exposure was a priori set to some known level (e.g., 0, and is therefore unexposed); however, our approach does not require an exposure validation study or replicate measures of exposure, which are typically needed when addressing bias due to exposure measurement error. A key condition for this method, which we refer to as "partial population exchangeability," requires that the association between a measured covariate and outcome in the reference population equals the association between that covariate and outcome in the target population in the absence of exposure. We illustrate the approach using simulations and an example.

Nonmalignant respiratory disease mortality in male Colorado Plateau uranium miners, 1960-2016

BACKGROUND

To evaluate trends of nonmalignant respiratory disease (NMRD) mortality among US underground uranium miners on the Colorado Plateau, and to estimate the exposure-response association between cumulative radon progeny exposure and NMRD subtype mortality.

METHODS

Standardized mortality ratios (SMRs) and excess relative rates per 100 working level months (excess relative rate [ERR]/100 WLM) were estimated in a cohort of 4021 male underground uranium miners who were followed from 1960 through 2016.

RESULTS

We observed elevated SMRs for all NMRD subtypes. Silicosis had the largest SMR (n = 52, SMR = 41.4; 95% confidence interval [CI]: 30.9, 54.3), followed by other pneumoconiosis (n = 49, SMR = 39.6; 95% CI: 29.6, 52.3) and idiopathic pulmonary fibrosis (IPF) (n = 64, SMR = 4.77; 95% CI 3.67, 6.09). SMRs for silicosis increased with duration of employment; SMRs for IPF increased with duration of employment and calendar period. There was a positive association between cumulative radon exposure and silicosis with evidence of modification by smoking (ERR/100 WLM≥10 pack-years  = 0.78; 95% CI: 0.05, 24.6 and ERR/100 WLM<10 pack-years  = 0.01; 95% CI: -0.03, 0.52), as well as a small positive association between radon and IPF (ERR/100 WLM = 0.06, 95% CI: 0.00, 0.24); these associations were driven by workers with prior employment in hard rock mining.

CONCLUSIONS

Uranium mining workers had excess NMRD mortality compared with the general population; this excess persisted throughout follow-up. Exposure-response analyses indicated a positive association between radon exposure and IPF and silicosis, but these analyses have limitations due to outcome misclassification and missing information on occupational co-exposures such as silica dust.

Drinking Water-Associated PFAS and Fluoroethers and Lipid Outcomes in the GenX Exposure Study

BACKGROUND

Residents of Wilmington, North, Carolina, were exposed to drinking water contaminated by fluoroethers and legacy per- and polyfluoroalkyl substances (PFAS), such as perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS), with fluoroether exposure occurring from 1980 to 2017. PFOA and PFOS have previously been associated with metabolic dysfunction; however, few prior studies have examined associations between other PFAS and lipid levels.

OBJECTIVES

We measured the association between serum fluoroether and legacy PFAS levels and various cholesterol outcomes.

METHODS

Participants in the GenX Exposure Study contributed nonfasting blood samples in November 2017 and May 2018 that were analyzed for 20 PFAS (10 legacy, 10 fluoroethers) and serum lipids [total cholesterol, low-density lipoprotein (LDL), high-density lipoprotein (HDL), triglycerides] and calculated non-HDL cholesterol. We estimated covariate-adjusted associations between quartiles of exposure to each of the PFAS measures (as well as the summed concentrations of legacy PFAS, fluoroethers, and all 10 targeted PFAS) and lipid outcomes by fitting inverse probability of treatment weighted linear regressions.

RESULTS

In this cross-sectional study of 326 participants (age range 6-86 y), eight PFAS were detected in > 50 % of the population. For PFOS and perfluorononanoic acid (PFNA), non-HDL cholesterol was approximately 5 mg / dL higher per exposure quartile increase: [PFOS: 4.89; 95% confidence interval (CI): 0.10, 9.68 and PFNA: 5.25 (95% CI: 0.39, 10.1)], whereas total cholesterol was approximately 6 mg / dL higher per quartile [PFOS: 5.71 (95% CI: 0.38, 11.0), PFNA: 5.92 (95% CI: 0.19, 11.7)]. In age-stratified analyses, associations were strongest among the oldest participants. Two fluoroethers were associated with higher HDL, whereas other fluoroether compounds were not associated with serum lipid levels.

DISCUSSION

PFNA and PFOS were associated with higher levels of total and non-HDL cholesterol, with associations larger in magnitude among older adults. In the presence of these legacy PFAS, fluoroethers appeared to be associated with HDL but not non-HDL lipid measures. https://doi.org/10.1289/EHP11033.

Exposure to industrial hog operations and gastrointestinal illness in North Carolina, USA

With 9 million hogs, North Carolina (NC) is the second leading hog producer in the United States. Most hogs are housed at concentrated animal feeding operations (CAFOs), where millions of tons of hog waste can pollute air and water with fecal pathogens that can cause diarrhea, vomiting, and/or nausea (known as acute gastrointestinal illness (AGI)). We used NC's ZIP code-level emergency department (ED) data to calculate rates of AGI ED visits (2016-2019) and swine permit data to estimate hog exposure. Case exposure was estimated as the inverse distances from each hog CAFO to census block centroids, weighting with Gaussian decay and by manure amount per CAFO, then aggregated to ZIP code using population weights. We compared ZIP codes in the upper quartile of hog exposure ("high hog exposed") to those without hog exposure. Using inverse probability of treatment weighting, we created a control with similar demographics to the high hog exposed population and calculated rate ratios using quasi-Poisson models. We examined effect measure modification of rurality and race using adjusted models. In high hog exposed areas compared to areas without hog exposure, we observed a 11% increase (95% CI: 1.06, 1.17) in AGI rate and 21% increase specifically in rural areas (95% CI: 0.98, 1.43). When restricted to rural areas, we found an increased AGI rate among American Indian (RR = 4.29, 95% CI: 3.69, 4.88) and Black (RR = 1.45, 95% CI: 0.98, 1.91) residents. The association was stronger during the week after heavy rain (RR = 1.41, 95% CI: 1.19, 1.62) and in areas with both poultry and swine CAFOs (RR = 1.52, 95% CI: 1.48, 1.57). Residing near CAFOs may increase rates of AGI ED visits. Hog CAFOs are disproportionally built near rural Black and American Indian communities in NC and are associated with increased AGI most strongly in these populations.

North Carolina's Changing Energy Generation Profile and Reductions in Key Air Pollutants, 2000-2019

BACKGROUND Coal combustion releases a number of airborne toxins. The North Carolina Clean Smokestacks Act (CSA) of 2002 required North Carolina coal-fired power plants (CFPP) to reduce nitrogen oxides (NO_X) emissions by 2009 and sulfur dioxide (SO₂) emissions to 2 benchmarks by 2009 and 2013.METHODS We utilized publicly available databases from the Energy Information Administration and the Environmental Protection Agency to characterize North Carolina's electricity generation profile from 2000 until 2019 and evaluate corresponding NO_x and SO₂ emissions by sector over the same time period.RESULTS Between 2000 and 2008 in North Carolina, approximately 60% of electric power was generated by CFPPs. Since then, North Carolina's electric power generation has transformed from predominant dependence on coal to approximately equal dependence on natural gas and nuclear power (each at ~ 30%), with coal close behind (~ 25%). Renewables have increased, although marginally relative to the rapid increase in natural gas. Despite the stark drop in reliance on CFPPs for energy in North Carolina and subsequent drop in emissions, CFPPs still contribute ~ 60% of SO₂ air pollution as of 2017.LIMITATIONS This analysis relies upon electricity generation and emissions data self-reported by utilities and publicly available from federal agenciesCONCLUSION North Carolina's electric utilities met the 2009 and 2013 regulatory benchmarks set by the CSA, which resulted in substantial reductions in SO₂ emissions from the fuel combustion electric generation sector. Still, CFPPs remain the primary utility-related and overall anthropogenic contributor of SO₂ air pollution in North Carolina.

Patterns of Children's Blood Lead Screening and Blood Lead Levels in North Carolina, 2011-2018-Who Is Tested, Who Is Missed?

BACKGROUND

No safe level of lead in blood has been identified. Blood lead testing is required for children on Medicaid, but it is at the discretion of providers and parents for others. Elevated blood lead levels (EBLLs) cannot be identified in children who are not tested.

OBJECTIVES

The aims of this research were to identify determinants of lead testing and EBLLs among North Carolina children and estimate the number of additional children with EBLLs among those not tested.

METHODS

We linked geocoded North Carolina birth certificates from 2011-2016 to 2010 U.S. Census data and North Carolina blood lead test results from 2011-2018. We estimated the probability of being screened for lead and created inverse probability (IP) of testing weights. We evaluated the risk of an EBLL of ≥3μg/dL at <30  months of age, conditional on characteristics at birth, using generalized linear models and then applied IP weights to account for missing blood lead results among unscreened children. We estimated the number of additional children with EBLLs of all North Carolina children using the IP-weighted population and bootstrapping to produce 95% credible intervals (CrI).

RESULTS

Mothers of the 63.5% of children (402,002 of 633,159) linked to a blood lead test result were disproportionately young, Hispanic, Black, American Indian, or on Medicaid. In full models, maternal age ≤20y [risk ratio (RR)=1.10; 95% confidence interval (CI): 1.13, 1.20] or smoking (RR=1.14; 95% CI: 1.12, 1.17); proximity to a major roadway (RR=1.10; 95% CI: 1.05, 1.15); proximity to a lead-releasing Toxics Release Inventory site (RR=1.08; 95% CI: 1.03, 1.14) or a National Emissions Inventory site (RR=1.11; 95% CI: 1.07, 1.14); and living in neighborhoods with more housing built before 1950 (RR=1.10; 95% CI: 1.05, 1.14) or before 1940 (RR=1.18; 95% CI: 1.11, 1.25) or more vacant housing (RR=1.14; 95% CI: 1.11, 1.17) were associated with an increased risk of EBLL, whereas overlap with a public water service system was associated with a decreased risk of EBLL (RR=0.85; 95% CI: 0.83, 0.87). Children of Black mothers were no more likely than children of White mothers to have EBLLs (RR=0.98; 95% CI: 0.96, 1.01). Complete blood lead screening in 2011-2018 may have identified an additional 17,543 (95% CrI: 17,462, 17,650) children with EBLLs ≥3μg/dL.

DISCUSSION

Our results indicate that current North Carolina lead screening strategies fail to identify over 30% (17,543 of 57,398) of children with subclinical lead poisoning and that accounting for characteristics at birth alters the conclusions about racial disparities in children's EBLLs. https://doi.org/10.1289/EHP10335.

Lung cancer mortality in the European cohort of titanium dioxide workers: a reanalysis of the exposure-response relationship

OBJECTIVES

Animal bioassays have demonstrated convincing evidence of the potential carcinogenicity to humans of titanium dioxide (TiO2), but limitations in cohort studies have been identified, among which is the healthy worker survivor effect (HWSE). We aimed to address this bias in a pooled study of four cohorts of TiO2 workers.

METHODS

We reanalysed data on respirable TiO2 dust exposure and lung cancer mortality among 7341 male workers employed in TiO2 production in Finland, France, UK and Italy using the parametric g-formula, considering three hypothetical interventions: setting annual exposures at 2.4 (U.S. occupational exposure limit), 0.3 (German limit) and 0 mg/m3 for 25 and 35 years.

RESULTS

The HWSE was evidenced. Taking this into account, we observed a positive association between lagged cumulative exposure to TiO2 and lung cancer mortality. The estimated number of lung cancer deaths at each age group decreased across increasingly stringent intervention levels. At age 70 years, the estimated number of lung cancer deaths expected in the cohort after 35-year exposure was 293 for exposure set at 2.4 mg/m3, 235 for exposure set at 0.3 mg/m3, and 211 for exposure set at 0 mg/m3.

CONCLUSION

This analysis shows that HWSE can hide an exposure-response relationship. It also shows that TiO2 epidemiological data could demonstrate an exposure-effects relationship if analysed appropriately. More epidemiological studies and similar reanalyses of existing cohort studies are warranted to corroborate the human carcinogenicity of TiO2. This human evidence, when combined with the animal evidence, strengthens the overall evidence of carcinogenicity of TiO2.

Lung Cancer and Radon: Pooled Analysis of Uranium Miners Hired in 1960 or Later

BACKGROUND

Despite reductions in exposure for workers and the general public, radon remains a leading cause of lung cancer. Prior studies of underground miners depended heavily upon information on deaths among miners employed in the early years of mine operations when exposures were high and tended to be poorly estimated.

OBJECTIVES

To strengthen the basis for radiation protection, we report on the follow-up of workers employed in the later periods of mine operations for whom we have more accurate exposure information and for whom exposures tended to be accrued at intensities that are more comparable to contemporary settings.

METHODS

We conducted a pooled analysis of cohort studies of lung cancer mortality among 57,873 male uranium miners in Canada, Czech Republic, France, Germany, and the United States, who were first employed in 1960 or later (thereby excluding miners employed during the periods of highest exposure and focusing on miners who tend to have higher quality assessments of radon progeny exposures). We derived estimates of excess relative rate per 100 working level months (ERR/100 WLM) for mortality from lung cancer.

RESULTS

The analysis included 1.9  million person-years of observation and 1,217 deaths due to lung cancer. The relative rate of lung cancer increased in a linear fashion with cumulative exposure to radon progeny (ERR/100 WLM = 1.33 ; 95% CI: 0.89, 1.88). The association was modified by attained age, age at exposure, and annual exposure rate; for attained ages < 55    y , the ERR/100 WLM was 8.38 (95% CI: 3.30, 18.99) among miners who were exposed at ≥ 35   years of age and at annual exposure rates of < 0.5 working levels. This association decreased with older attained ages, younger ages at exposure, and higher exposure rates.

DISCUSSION

Estimates of association between radon progeny exposure and lung cancer mortality among relatively contemporary miners are coherent with estimates used to inform current protection guidelines. https://doi.org/10.1289/EHP10669.

Trends in fatal occupational injuries in Latino/a workers relative to other groups, North Carolina 2000-2017

BACKGROUND

Latino/a workers may experience higher fatal occupational injury rates than non-Latino/a workers. In North Carolina, the Latino/a population more than doubled between 2000 and 2017. We examined fatal occupational injuries among Latino/a and non-Latino/a workers in North Carolina over this period.

METHODS

Information on fatal occupational injuries was abstracted from records of the North Carolina Office of the Chief Medical Examiner and the death certificate records held by the North Carolina Office of Vital Records. Estimates of the working population were derived from the decennial census and American Community Survey. Estimates of annual rates of fatal occupational injury for the period January 1, 2000 to December 31, 2017 were derived for Latino/a workers and compared to Black and White workers not identified as Latino/a.

RESULTS

Over the study period, 1,783 fatal occupational injuries were identified among non-Latino/a workers and 259 fatal occupational injuries among Latino/a workers in North Carolina. The majority of fatal occupational injuries among Latino/a workers occurred among males employed in construction and agriculture. While the fatal occupational injury rate among Latino/a workers declined over the study period, the rate among Latino/a workers was higher than among non-Latino/a White and Black workers; moreover, fatal occupational injury rates for Latino/a workers trended upwards during the most recent years of the study period.

CONCLUSIONS

Latino/a workers in North Carolina have the highest fatal occupational injury rate of any race/ethnicity group.

Bespoke Instruments: A new tool for addressing unmeasured confounders

Suppose that an investigator is interested in quantifying an exposure-disease causal association in a setting where the exposure, disease, and some potential confounders of the association of interest have been measured. However, there remains concern about residual confounding of the association of interest by unmeasured confounders. We propose an approach to account for residual bias due to unmeasured confounders. The proposed approach uses a measured confounder to derive a "bespoke" instrumental variable that is tailored to the study population and is used to control for bias due to residual confounding. The approach may provide a useful tool for assessing and accounting for bias due to residual confounding. We provide a formal description of the conditions for identification of causal effects, illustrate the method using simulations, and provide an empirical example concerning mortality among Japanese atomic bomb survivors.

Effects of short-term ambient PM2.5 exposure on cardiovascular disease incidence and mortality among U.S. hemodialysis patients: a retrospective cohort study

BACKGROUND

Ambient PM_2.5 is a ubiquitous air pollutant with demonstrated adverse health impacts in population. Hemodialysis patients are a highly vulnerable population and may be particularly susceptible to the effects of PM_2.5 exposure. This study examines associations between short-term PM_2.5 exposure and cardiovascular disease (CVD) and mortality among patients receiving maintenance in-center hemodialysis.

METHODS

Using the United State Renal Data System (USRDS) registry, we enumerated a cohort of all US adult kidney failure patients who initiated in-center hemodialysis between 1/1/2011 and 12/31/2016. Daily ambient PM_2.5 exposure estimates were assigned to cohort members based on the ZIP code of the dialysis clinic. CVD incidence and mortality were ascertained through 2016 based on USRDS records. Discrete time hazards regression was used to estimate the association between lagged PM_2.5 exposure and CVD incidence, CVD-specific mortality, and all-cause mortality 1 t adjusting for temperature, humidity, day of the week, season, age at baseline, race, employment status, and geographic region. Effect measure modification was assessed for age, sex, race, and comorbidities.

RESULTS

Among 314,079 hemodialysis patients, a 10 µg/m³ increase in the average lag 0-1 daily PM_2.5 exposure was associated with CVD incidence (HR: 1.03 (95% CI: 1.02, 1.04)), CVD mortality (1.05 (95% CI: 1.03, 1.08)), and all-cause mortality (1.04 (95% CI: 1.03, 1.06)). The association was larger for people who initiated dialysis at an older age, while minimal evidence of effect modification was observed across levels of sex, race, or baseline comorbidities.

CONCLUSIONS

Short-term ambient PM_2.5 exposure was positively associated with incident CVD events and mortality among patients receiving in-center hemodialysis. Older patients appeared to be more susceptible to PM_2.5-associated CVD events than younger hemodialysis patients.

Estimating Associations Between Annual Concentrations of Particulate Matter and Mortality in the United States, Using Data Linkage and Bayesian Maximum Entropy

BACKGROUND

Exposure to fine particulate matter (PM2.5) is an established risk factor for human mortality. However, previous US studies have been limited to select cities or regions or to population subsets (e.g., older adults).

METHODS

Here, we demonstrate how to use the novel geostatistical method Bayesian maximum entropy to obtain estimates of PM2.5 concentrations in all contiguous US counties, 2000-2016. We then demonstrate how one could use these estimates in a traditional epidemiologic analysis examining the association between PM2.5 and rates of all-cause, cardiovascular, respiratory, and (as a negative control outcome) accidental mortality.

RESULTS

We estimated that, for a 1 log(μg/m3) increase in PM2.5 concentration, the conditional all-cause mortality incidence rate ratio (IRR) was 1.029 (95% confidence interval [CI]: 1.006, 1.053). This implies that the rate of all-cause mortality at 10 µg/m3 would be 1.020 times the rate at 5 µg/m3. IRRs were larger for cardiovascular mortality than for all-cause mortality in all gender and race-ethnicity groups. We observed larger IRRs for all-cause, nonaccidental, and respiratory mortality in Black non-Hispanic Americans than White non-Hispanic Americans. However, our negative control analysis indicated the possibility for unmeasured confounding.

CONCLUSION

We used a novel method that allowed us to estimate PM2.5 concentrations in all contiguous US counties and obtained estimates of the association between PM2.5 and mortality comparable to previous studies. Our analysis provides one example of how Bayesian maximum entropy could be used in epidemiologic analyses; future work could explore other ways to use this approach to inform important public health questions.

Hurricane flooding and acute gastrointestinal illness in North Carolina

Hurricanes often flood homes and industries, spreading pathogens. Contact with pathogen-contaminated water can result in diarrhea, vomiting, and/or nausea, known collectively as acute gastrointestinal illness (AGI). Hurricanes Matthew and Florence caused record-breaking flooding in North Carolina (NC) in October 2016 and September 2018, respectively. To examine the relationship between hurricane flooding and AGI in NC, we first calculated the percent of each ZIP code flooded after Hurricanes Matthew and Florence. Rates of all-cause AGI emergency department (ED) visits were calculated from NC's ED surveillance system data. Using controlled interrupted time series, we compared AGI ED visit rates during the three weeks after each hurricane in ZIP codes with a third or more of their area flooded to the predicted rates had these hurricanes not occurred, based on AGI 2016-2019 ED trends, and controlling for AGI ED visit rates in unflooded areas. We examined alternative case definitions (bacterial AGI) and effect measure modification by race and age. We observed an 11% increase (rate ratio (RR): 1.11, 95% CI: 1.00, 1.23) in AGI ED visit rates after Hurricanes Matthew and Florence. This effect was particularly strong among American Indian patients and patients aged 65 years and older after Florence and elevated among Black patients for both hurricanes. Florence's effect was more consistent than Matthew's effect, possibly because little rain preceded Florence and heavy rain preceded Matthew. When restricted to bacterial AGI, we found an 85% (RR: 1.85, 95% CI: 1.37, 2.34) increase in AGI ED visit rate after Florence, but no increase after Matthew. Hurricane flooding is associated with an increase in AGI ED visit rate, although the strength of effect may depend on total storm rainfall or antecedent rainfall. American Indians and Black people-historically pushed to less desirable, flood-prone land-may be at higher risk for AGI after storms.

Lung and extrathoracic cancer incidence among underground uranium miners exposed to radon progeny in the Příbram region of the Czech Republic: a case-cohort study

OBJECTIVES

Radon is carcinogenic, but more studies are needed to understand relationships with lung cancer and extrathoracic cancers at low exposures. There are few studies evaluating associations with cancer incidence or assessing the modifying effects of smoking.

METHODS

We conducted a case-cohort study with 16 434 underground uranium miners in the Czech Republic with cancer incidence follow-up 1977-1996. Associations between radon exposure and lung cancer, and extrathoracic cancer, were estimated with linear excess relative rate (ERR) models. We examined potential modifying effects of smoking, time since exposure and exposure rate.

RESULTS

Under a simple ERR model, assuming a 5-year exposure lag, the estimated ERR of lung cancer per 100 working level months (WLM) was 0.54 (95% CI 0.33 to 0.83) and the estimated ERR of extrathoracic cancer per 100 WLM was 0.07 (95% CI -0.17 to 0.72). Most lung cancer cases were observed among smokers (82%), and the estimated ERR of lung cancer per 100 WLM was larger among smokers (ERR/100 WLM=1.35; 95% CI 0.84 to 2.15) than among never smokers (ERR/100 WLM=0.12; 95% CI -0.05 to 0.49). Among smokers, the estimated ERR of lung cancer per 100 WLM decreased with time since exposure from 3.07 (95% CI -0.04 to 10.32) in the period 5-14 years after exposure to 1.05 (95% CI 0.49 to 1.87) in the period 25+ years after exposure.

CONCLUSIONS

We observed positive associations between cumulative radon exposure and lung cancer, consistent with prior studies. We observed a positive association between cumulative radon exposure and extrathoracic cancers, although the estimates were small. There was evidence that the association between radon and lung cancer was modified by smoking in a multiplicative or super-multiplicative fashion.

Amplification of Bias Due to Exposure Measurement Error

Observational epidemiologic studies typically face challenges of exposure measurement error and confounding. Consider an observational study of the association between a continuous exposure and an outcome, where the exposure variable of primary interest suffers from classical measurement error (i.e., the measured exposures are distributed around the true exposure with independent error). In the absence of exposure measurement error, it is widely recognized that one should control for confounders of the association of interest to obtain an unbiased estimate of the effect of that exposure on the outcome of interest. However, here we show that, in the presence of classical exposure measurement error, the net bias in an estimate of the association of interest may increase upon adjustment for confounders. We offer an analytical expression for calculating the change in net bias in an estimate of the association of interest upon adjustment for a confounder in the presence of classical exposure measurement error, and we illustrate this problem using simulations.

Pregnancy exposure to organophosphate esters and the risk of attention-deficit hyperactivity disorder in the Norwegian mother, father and child cohort study

BACKGROUND

Organophosphate esters (OPEs) are a class of flame retardants in common use. OPEs can easily leach from materials, resulting in human exposure. Increasing concentrations have been reported in human populations over the past decade. Recent studies have linked prenatal OPE exposure to hyperactivity and attention problems in children. Such behaviors are often found among children with attention-deficit hyperactivity disorder (ADHD), however, no study has investigated OPEs in relation to clinically assessed ADHD.

OBJECTIVE

To evaluate prenatal exposure to OPEs as risk factors for clinically assessed ADHD using a case-cohort study nested within the Norwegian Mother, Father, and Child Cohort Study (MoBa).

METHODS

We included in the case group 295 ADHD cases obtained via linkage with the Norwegian Patient Registry, and the sub-cohort group 555 children sampled at baseline, irrespective of their ADHD case status. Prenatal concentrations of OPE metabolites were measured in maternal urine collected at 17 weeks of gestation, and included diphenyl phosphate (DPHP), di-n-butyl phosphate (DNBP), bis(2-butoxyethyl) hydrogen phosphate (BBOEP), and bis(1,3-dichloro-2-propyl) phosphate (BDCIPP). We estimated risk ratios and the corresponding 95% confidence intervals [95% CI] using logistic regression, adjusting for season of urine collection, child sex, birth year, and maternal depression, education, and sum of urinary di(2-ethylhexyl) phthalate metabolites (∑DEHP) concentration during pregnancy. To assess the overall impact of simultaneously decreasing exposure to all chemical constituents of an OPE-phthalate mixture, quantile based g-computation was implemented. The mixture constituents included OPE and phthalate metabolites commonly detected in our study. In all models, we considered effect measure modification by child sex and polymorphisms in genes encoding paraoxonase 1 (PON1) and cytochrome P450 (P450) enzymes. Mediation analysis was conducted using thyroid function biomarkers estimated from maternal blood collected at 17 weeks of gestation.

RESULTS

DPHP was detected in nearly all samples (97.2%), with a higher geometric mean among the case group (0.70 µg/L) as compared to the sub-cohort (0.52 µg/L). DNBP was commonly detected as well (93.8%), while BBOEP (52.9%) and BDCIPP (22.9%) were detected less frequently. A higher risk of ADHD was observed in children with greater than median exposure to DPHP during pregnancy (risk ratio: 1.38 [95% CI: 0.96, 1.99]), which was slightly higher among girls (2.04 [1.03, 4.02]) and children of mothers with PON1 Q192R genotype QR (1.69 [0.89, 3.19]) or PON1 Q192R genotype RR (4.59 [1.38, 15.29]). The relationship between DPHP and ADHD (total risk ratio: 1.34 [0.90, 2.02]) was partially mediated through total triiodothyronine to total thyroxine ratio (natural direct effect: 1.29 [0.87, 1.94]; natural indirect effect: 1.04 [1.00, 1.10]; 12.48% mediated). We also observed an elevated risk of ADHD in relation to BDCIPP detection during pregnancy (1.50 [0.98, 2.28]). We did not observe notable differences in ADHD by DNBP (0.88 [0.62, 1.26]) or BBOEP (1.03 [0.73, 1.46]) during pregnancy. Simultaneously decreasing all constituents of common-detect OPE-phthalate mixture, specifically DPHP, DNBP, and 6 phthalate metabolites, by a quartile resulted in an ADHD risk ratio of 0.68 [0.64, 0.72].

CONCLUSION

Prenatal exposure to DPHP and BDCIPP may increase the risk of ADHD. For DPHP, we observed potential modification by child sex and maternal PON1 Q192R genotype and partial mediation through maternal thyroid hormone imbalance at 17 weeks gestation.

Decreased Susceptibility of Marginal Odds Ratios to Finite-sample Bias

Parameters representing adjusted treatment effects may be defined marginally or conditionally on covariates. The choice between a marginal or covariate-conditional parameter should be driven by the study question. However, an unappreciated benefit of marginal estimators is a reduction in susceptibility to finite-sample bias relative to the unpenalized maximum likelihood estimator of the covariate-conditional odds ratio (OR). Using simulation, we compare the finite-sample bias of different marginal and conditional estimators of the OR. We simulated a logistic model to have 15 events per parameter and two events per parameter. We estimated the covariate-conditional OR by maximum likelihood with and without Firth's penalization. We used three estimators of the marginal OR: g-computation, inverse probability of treatment weighting, and augmented inverse probability of treatment weighting. At 15 events per parameter, as expected, all estimators were effectively unbiased. At two events per parameter, the unpenalized covariate-conditional estimator was notably biased but penalized covariate-conditional and marginal estimators exhibited minimal bias.

Pregnancy exposure to common-detect organophosphate esters and phthalates and maternal thyroid function

BACKGROUND

Contemporary human populations are exposed to elevated concentrations of organophosphate esters (OPEs) and phthalates. Some metabolites have been linked with altered thyroid function, however, inconsistencies exist across thyroid function biomarkers. Research on OPEs is sparse, particularly during pregnancy, when maintaining normal thyroid function is critical to maternal and fetal health. In this paper, we aimed to characterize relationships between OPEs and phthalates exposure and maternal thyroid function during pregnancy, using a cross-sectional investigation of pregnant women nested within the Norwegian Mother, Father, and Child Cohort (MoBa).

METHODS

We included 473 pregnant women, who were euthyroid and provided bio-samples at 17 weeks' gestation (2004-2008). Four OPE and six phthalate metabolites were measured from urine; six thyroid function biomarkers were estimated from blood. Relationships between thyroid function biomarkers and log-transformed concentrations of OPE and phthalate metabolites were characterized using two approaches that both accounted for confounding by co-exposures: co-pollutant adjusted general linear model (GLM) and Bayesian Kernal Machine Regression (BKMR).

RESULTS

We restricted our analysis to common-detect OPE and phthalate metabolites (>94%): diphenyl phosphate (DPHP), di-n-butyl phosphate (DNBP), and all phthalate metabolites. In GLM, pregnant women with summed di-isononyl phthalate metabolites (∑DiNP) concentrations in the 75th percentile had a 0.37 ng/μg lower total triiodothyronine (TT3): total thyroxine (TT4) ratio (95% credible interval: [-0.59, -0.15]) as compared to those in the 25th percentile, possibly due to small but diverging influences on TT3 (-1.99 ng/dL [-4.52, 0.53]) and TT4 (0.13 μg/dL [-0.01, 0.26]). Similar trends were observed for DNBP and inverse associations were observed for DPHP, monoethyl phthalate, mono-isobutyl phthalate, and mono-n-butyl phthalate. Most associations observed in co-pollutants adjusted GLMs were attenuated towards the null in BKMR, except for the case of ∑DiNP and TT3:TT4 ratio (-0.48 [-0.96, 0.003]).

CONCLUSIONS

Maternal thyroid function varied modestly with ∑DiNP, whereas results for DPHP varied by the type of statistical models.

Mortality among uranium miners in North America and Europe: the Pooled Uranium Miners Analysis (PUMA)

BACKGROUND

The Pooled Uranium Miners Analysis (PUMA) study draws together information from cohorts of uranium miners from Canada, the Czech Republic, France, Germany and the USA.

METHODS

Vital status and cause of death were ascertained and compared with expectations based upon national mortality rates by computing standardized mortality ratios (SMRs) overall and by categories of time since first hire, calendar period of first employment and duration of employment as a miner.

RESULTS

There were 51 787 deaths observed among 118 329 male miners [SMR = 1.05; 95% confidence interval (CI): 1.04, 1.06]. The SMR was elevated for all cancers (n = 16 633, SMR = 1.23; 95% CI: 1.21, 1.25), due primarily to excess mortality from cancers of the lung (n = 7756, SMR = 1.90; 95% CI: 1.86, 1.94), liver and gallbladder (n = 549, SMR = 1.15; 95% CI: 1.06, 1.25), larynx (n = 229, SMR = 1.10; 95% CI: 0.97, 1.26), stomach (n = 1058, SMR = 1.08; 95% CI: 1.02, 1.15) and pleura (n = 39, SMR = 1.06; 95% CI: 0.75, 1.44). Lung-cancer SMRs increased with duration of employment, decreased with calendar period and persisted with time since first hire. Among non-malignant causes, the SMR was elevated for external causes (n = 3362, SMR = 1.41; 95% CI: 1.36, 1.46) and respiratory diseases (n = 4508, SMR = 1.32; 95% CI: 1.28, 1.36), most notably silicosis (n = 814, SMR = 13.56; 95% CI: 12.64, 14.52), but not chronic obstructive pulmonary disease (n = 1729, SMR = 0.98; 95% CI: 0.93, 1.02).

CONCLUSIONS

Whereas there are important obstacles to the ability to detect adverse effects of occupational exposures via SMR analyses, PUMA provides evidence of excess mortality among uranium miners due to a range of categories of cause of death. The persistent elevation of SMRs with time since first hire as a uranium miner underscores the importance of long-term follow-up of these workers.

Reducing Bias Due to Exposure Measurement Error Using Disease Risk Scores

Suppose that an investigator wants to estimate an association between a continuous exposure variable and an outcome, adjusting for a set of confounders. If the exposure variable suffers classical measurement error, in which the measured exposures are distributed with independent error around the true exposure, then an estimate of the covariate-adjusted exposure-outcome association may be biased. We propose an approach to estimate a marginal exposure-outcome association in the setting of classical exposure measurement error using a disease score-based approach to standardization to the exposed sample. First, we show that the proposed marginal estimate of the exposure-outcome association will suffer less bias due to classical measurement error than the covariate-conditional estimate of association when the covariates are predictors of exposure. Second, we show that if an exposure validation study is available with which to assess exposure measurement error, then the proposed marginal estimate of the exposure-outcome association can be corrected for measurement error more efficiently than the covariate-conditional estimate of association. We illustrate both of these points using simulations and an empirical example using data from the Orinda Longitudinal Study of Myopia (California, 1989-2001).