Causal relationship between two occupational carcinogens and lung cancer
Description
This directed acyclic graph is used to inform the estimation of the causal effect of lifetime cumulative exposure to asbestos (main exposure) on lung cancer (outcome) in the context of individual participant data (IPD) meta-analysis of the SYNERGY pooled case-control study [1]. Co-exposure to other occupational lung carcinogens is represented with respiratory cristaline silica (RCS) as a motivating example. True exposure is unobserved (X_asb_true and X_rcs_true) and thus estimated with the job exposure matrix SYN-JEM (asbestos_cum0 and silica_cum0), resulting in loss of individual variation (Epsilon), and exposure error terms (U) that are independent of outcome status (Berkson-type error) [2]. The minimal set of adjustment for this model consists of age (agegroup), sex, smoking status and time since quitting (time_quit), lifetime smoking intensity (packyears), and study (study_name). Study is considered a proxy for other study-level covariates and differences in baseline lung cancer risk. Occupational lung carcinogens are assumed to be related through the common cause industry which serves as a proxy for shared working environment conditions. This backdoor pathway was intentionally not closed as co-exposure to occupational lung carcinogens is frequent [3]. Estimating conditional effects by setting such co-exposures to zero could violate positivity (if certain covariate combinations never occur for a given job title) and consistency (counterfactual scenarios such as high asbestos exposure with zero RCS are implausible in some industries). The model assumes monotonicity, implying that asbestos cannot prevent lung cancer through any causal mechanism, and no-interference (only individual exposure is relevant for the effect and not affected by others' exposures). Lastly, both age and sex are matching factors in the individual case-control studies, thus potentially causing selection bias (not shown in this DAG for simplicity). Adjustment for these matching variables through unconditional logistic regression results in unbiased estimates under the rare disease assumption [4].
Additional details and documentation on this analysis are available at: https://github.com/UtrechtUniversity/PoC-Asbestos
References:
1 International Agency for Research on Cancer (IARC). About SYNERGY. [cited 2024 Oct 17]. SYNERGY: Home. Available from: https://synergy.iarc.who.int/about
2 Armstrong BG. Effect of measurement error on epidemiological studies of environmental and occupational exposures. Occupational and Environmental Medicine 1998; 55: 651–656.
3 Olsson A, Bouaoun L, Schüz J, et al. Lung Cancer Risks Associated with Occupational Exposure to Pairs of Five Lung Carcinogens: Results from a Pooled Analysis of Case-Control Studies (SYNERGY). Environmental Health Perspectives 2024; 132: 17005.
4 Kamenetsky ME, Keil AP. (Re-)match: adjusting for matching factors in case-control studies may be unnecessary or insufficient. American Journal of Epidemiology 2025; kwaf116. doi: 10.1093/aje/kwaf116
Visualisation
dag {
Epsilon_asb [latent,pos="0.999,6.027"]
Epsilon_rcs [latent,pos="0.000,5.932"]
U_asb [latent,pos="0.623,3.993"]
U_rcs [latent,pos="-0.323,3.991"]
X_asbestos_true [latent,pos="1.000,5.000"]
X_rcs_true [latent,pos="0.000,5.000"]
agegroup [adjusted,pos="0.115,1.155"]
asbestos_cum0 [exposure,pos="1.000,4.000"]
industry [latent,pos="0.495,4.727"]
packyears [adjusted,pos="1.977,2.040"]
sex [adjusted,pos="0.618,0.804"]
silica_cum0 [exposure,pos="0.000,4.000"]
status [outcome,pos="2.000,5.000"]
study_name [adjusted,pos="-0.106,1.928"]
time_quit [adjusted,pos="1.496,0.909"]
Epsilon_asb -> X_asbestos_true
Epsilon_rcs -> X_rcs_true
U_asb -> asbestos_cum0
U_rcs -> silica_cum0
X_asbestos_true -> status
X_rcs_true -> status [pos="0.864,5.972"]
agegroup -> asbestos_cum0
agegroup -> packyears
agegroup -> silica_cum0
agegroup -> status
agegroup -> time_quit
asbestos_cum0 -> X_asbestos_true
industry -> asbestos_cum0
industry -> silica_cum0
packyears -> asbestos_cum0
packyears -> silica_cum0
packyears -> status
packyears -> time_quit
sex -> asbestos_cum0
sex -> packyears
sex -> silica_cum0
sex -> status
sex -> time_quit
silica_cum0 -> X_rcs_true
study_name -> asbestos_cum0
study_name -> silica_cum0
study_name -> status
time_quit -> asbestos_cum0
time_quit -> silica_cum0
time_quit -> status
}
Discussion
No comments yet. Be the first to comment!