Background:

While cigarette smoking is the leading cause of lung cancer, the majority of smokers do not develop the disease over their lifetime. The inter-individual differences in risk among smokers may in part be due to variations in exposure to smoking-related toxicants.

Methods:

Using data from a subcohort of 2,309 current smokers at the time of urine collection from the Multiethnic Cohort Study, we prospectively evaluated the association of ten urinary biomarkers of smoking-related toxicants [total nicotine equivalents (TNE), a ratio of total trans-3′-hydroxycotinine (3-HCOT)/cotinine (a phenotypic measure of CYP2A6 enzymatic activity), 4-(methylnitrosamino)-1–3-(pyridyl)-1-butanol (NNAL), S-phenylmercapturic acid (SPMA), 3-hydroxypropyl mercapturic acid (3-HPMA), phenanthrene tetraol (PheT), 3-hydroxyphenanthrene (PheOH), the ratio of PheT/PheOH, cadmium (Cd), and (Z)-7-(1R,2R,3R,5S)-3,5-dihydroxy-2-[(E,3S)-3-hydroxyoct-1-enyl]cyclopenyl]hept-5-enoic acid (8-iso-PGF)] with lung cancer risk (n = 140 incident lung cancer cases over an average of 13.4 years of follow-up). Lung cancer risk was estimated using Cox proportional hazards models.

Results:

After adjusting for decade of birth, sex, race/ethnicity, body mass index, self-reported pack-years, creatinine, and urinary TNE (a biomarker of internal smoking dose), a one SD increase in log total 3-HCOT/cotinine (HR, 1.33; 95% CI, 1.06–1.66), 3-HPMA (HR, 1.41; 95% CI, 1.07–1.85), and Cd (HR, 1.45; 95% CI, 1.18–1.79) were each associated with increased lung cancer risk.

Conclusions:

Our study demonstrates that urinary total 3-HCOT/cotinine, 3-HPMA, and Cd are positively associated with lung cancer risk. These findings warrant replication and consideration as potential biomarkers for smoking-related lung cancer risk.

Impact:

These biomarkers may provide additional information on lung cancer risk that is not captured by self-reported smoking history or TNE.

See related commentary by Etemadi et al., p. 289

This article is featured in Highlights of This Issue, p. 287

Lung cancer is the second most common cancer and the leading cause of cancer-related death in both men and women in the United States (1). It is well established that cigarette smoking is the primary risk factor for lung cancer. However, disease risk is not equal among individuals, as it is estimated that 11% to 24% of smokers will develop the malignancy over their lifetime (2). Along with nicotine, each cigarette delivers a complex mixture of tobacco-related carcinogens and toxicants. There are 80 established carcinogenic constituents in cigarette smoke, including several lung carcinogens such as 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), multiple polycyclic aromatic hydrocarbons (PAH), and volatiles such as 1,2-butadiene and acrolein (3). Inter-individual differences in the risk of smoking-related lung cancer may be, in part, attributable to the variation in harmful constituent exposure. Moreover, the investigation of the association of these biomarkers with lung cancer risk may improve our understanding of the mechanisms by which tobacco exposure causes lung cancer (4, 5).

Biomarkers of smoking, such as metabolites of nicotine, and those of tobacco carcinogens and toxicants, have been shown to better reflect the internal dose of these harmful constituents and in some studies, their levels correlate with lung cancer risk. For example, prior epidemiological studies have shown that internal smoking dose, as determined by either cotinine (6–15) or total nicotine equivalents (TNE; the molar sum of nicotine and five or six metabolites in urine; refs. 12, 13, 16), was associated with lung cancer risk, independent of self-reported smoking history. In addition, studies have shown that 4-(methylnitrosamino)-1–3-pyridyl)-1-butanol (NNAL; a biomarker for NNK; refs. 8–13) and phenanthrene tetraol (PheT; a biomarker for PAH; refs. 9–11, 17) are associated with lung cancer, even after adjusting for cotinine or TNE. However, these prior studies were conducted in single populations (e.g., only Whites or Chinese from Asia), and do not address the variation in smoking behaviors and lung cancer risk across populations. In addition, cotinine, along with self-reported smoking history does not fully account for individual variability in nicotine metabolism, which influences smoking dose and exposure to tobacco toxicants (12, 18–21). Cytochrome P450 2A6 (CYP2A6) is the primary catalyst of nicotine metabolism and among current smokers in the Multiethnic Cohort (MEC) study, smoking behavior and dose were influenced by a smoker's urinary nicotine metabolite ratio [total trans-3′-hydroxycotinine (3-HCOT)/cotinine; ref. 22]. This ratio is a phenotypic measure of CYP2A6 enzymatic activity and in MEC current smokers we previously reported that it was associated with lung cancer risk even after adjustment for self-reported smoking history and TNE (16). Acrolein, a well-known lung irritant, was recently reclassified by the International Agency for Research on Cancer (IARC) as probably carcinogenic to humans (23). However, the contribution of acrolein exposure to lung cancer risk in smokers has not been well studied. Therefore, a systematic examination of common biomarkers of tobacco toxicant exposure and metabolism is needed and may provide information on lung cancer risk beyond that of self-reported smoking history.

To better understand the individual and joint effects of smoking dose and specific smoking-related toxicant exposures on lung cancer risk, we evaluated the associations of a number of tobacco exposure biomarkers with the risk of smoking-related lung cancer among smokers in the MEC Study (n = 140 incident lung cancer cases over an average of 13.4 years of follow-up). These biomarkers included NNAL, S-phenylmercapturic acid (SPMA, a biomarker of benzene uptake), 3-hydroxypropyl mercapturic acid (3-HPMA, a biomarker of acrolein uptake), PheT, 3-hydroxyphenanthrene (PheOH, a biomarker of PAH detoxification), the ratio of PheT/PheOH (proposed as a biomarker of metabolic activation of polycyclic aromatic hydrocarbons), cadmium (Cd), and (Z)-7-[1R,2R,3R,5S)-3,5-dihydroxy-2-[(E,3S)-3-hydroxyoct-1-enyl]cyclopentyl]hept-5-enoic acid (8-iso-PGF, a biomarker of oxidative stress).

Study population

Details of the MEC have been previously described (24). Briefly, the cohort consists of 215,251 men and women recruited from Hawaii and California (primarily Los Angeles) between 1993 and 1996. Participants were between the ages of 45 and 75 years old at recruitment and primarily belonged to five ethnic/racial groups: African American, Japanese American, Latino, Native Hawaiian, and White. Approximately ten years after cohort entry, a subset of participants was asked to participate in a biorepository project by providing a blood and overnight urine sample (Hawaii) or a first-morning urine sample (California).

The present study includes a subcohort of MEC participants (N = 2,309) who were lung cancer-free current smokers at the time of urine collection and who have complete urinary biomarkers of both TNE and ratio of total 3-HCOT/cotinine measured (16). Approval for this study, including the consent procedure, was obtained from the Institutional Review Boards of the University of Minnesota, University of Hawaii, and University of Southern California. All study participants provided written informed consent.

Epidemiologic data

All participants completed an epidemiologic questionnaire at the time of MEC study enrollment and at the time of urine collection. These questionnaires included detailed information on demographic characteristics, medical history, average daily cigarettes smoked, cigarette smoking during the past two weeks (questionnaire at urine collection only), smoking duration (years), and a medication record. Measures of body-mass index (BMI) and smoking duration that were missing at the time of urine collection for a few subjects were imputed from other MEC questionnaires, as described previously (16, 20).

Urinary metabolites of smoking-related toxicants

Details of the analytic methods used to determine the levels of urinary biomarkers of nicotine metabolism and tobacco smoking toxicants and data for these biomarkers have been previously published (20, 25–29). Briefly, overnight or first-morning urine was used to measure TNE (the molar sum of nicotine N-oxide, total nicotine, total cotinine, and total 3-HCOT; ref. 20), total NNAL (25), SPMA (26), 3-HPMA (27), PheT, and PheOH (28) using LC/MS-MS [where “total” refers to the compound and its glucuronide conjugate(s)]. Urinary Cd was measured using inductively-coupled plasma mass spectrometry (ICP-MS; ref. 29). Urinary 8-iso-PGF levels were measured as described previously (30). The detection limits were 13 ng/mL for nicotine, 20 ng/mL for cotinine, 18 ng/mL for 3-HCOT, 0.14 pmol/mL for total NNAL, 0.01 pmol/mL for SPMA, 4.5 pmol/mL for 3-HPMA, 0.005 pmol/mL for PheT, 0.05 pmol/mL for PheOH, 0.02 ng/mL for Cd, and 0.03 pmol/mL for 8-iso-PGF. The coefficient of variation (CV) for the assays across runs were 16.7% for nicotine, 10.1% for cotinine, 11.4% for 3-HCOT, 16.2% for total NNAL, 15.0% for SPMA, 9.1% for 3-HPMA, 12.2% for PheT, 19.7% for PheOH, 3.1% for Cd, and 7.0% for 8-iso-PGF. The proportion of 2,309 participants with missing biomarkers or those < LOD can be found in Supplementary Table S1. Biomarker levels below the limit of detection (LOD) were replaced with LOD/2. Creatinine was quantified using a colorimetric microplate assay (CRE34-K01) purchased from Eagle Bioscience (https://eaglebio.com/product/creatinine-microplate-assay-kit/; ref. 20). A urinary biomarker of CYP2A6 enzymatic activity was computed as the ratio of total 3-HCOT/cotinine. A urinary biomarker of metabolic activation of PAH was computed as the ratio of PheT/PheOH.

Follow-up and identification of lung cancer cases and deaths

Participants’ follow-up began at the age of urine collection, for which urinary biomarkers were measured, and ended when one of the following events occurred: (i) diagnosis of lung cancer; (ii) death; or (iii) end of follow-up, December 31, 2017. Incident lung cancer cases were identified through linkages to two state-wide NCI Surveillance, Epidemiology and End Results (SEER) Program registries: the Hawaii Tumor Registry and the California State Cancer Registry. The International Classification of Diseases for Oncology (ICD-O-3; ref. 31) and ICD-10 (32) code C34 for malignant neoplasm of bronchus and lung were used for this purpose. Deaths were identified by linkages to the state death certificate files in Hawaii and California and to the National Death Index (NDI) for deaths occurring in other states. By the end of follow-up (average 13.4 years), 140 incident primary lung cancer cases were identified. ICD-O-3 codes were obtained from the tumor registry and classified into four common lung cancer histologic cell types [adenocarcinoma (ADC: 8140, 8250, 8481, and 8490), squamous cell carcinoma (SCC: 8070 and 8071), small-cell lung cancer (SCLC: 8041 and 8045), and large-cell carcinoma (8012 and 8013)] and unspecified malignant neoplasm (8000, 8010, 8020, 8033, 8046, and 8246; ref. 33)

Statistical analysis

The covariate-adjusted geometric mean values with estimated 95% confidence intervals (95% CI) were computed for each metabolite, adjusted for age, sex (male/female), race/ethnicity (African American, Native Hawaiian, White, Latino, Japanese American), body mass index (BMI, kg/m2; log), creatinine (mg/dL; log). Values of all metabolites were transformed by taking the natural log to better meet model assumptions. Geometric mean values presented in the tables are back-transformed to their natural scale for ease of interpretation. The risk of lung cancer was estimated using HRs and 95% CIs from Cox proportional hazards models, where age was used as the time metric, and follow-up began at the age of urine collection. TNE was used to confirm current smoking status. Subjects with TNE < 1.27 nmol/mL (four times the limit of quantitation) were excluded, as this would indicate that these participants were not current smokers at the time of urine collection. To compare our results across the different biomarkers with variation in ranges, we standardized each log-transformed biomarker by dividing an individual's value by the SD of the respective log-transformed biomarker for the overall population (Supplementary Table S1). Therefore, we present the HR for incident lung cancer associated with a one-unit SD change of the log urinary biomarker levels. Distribution plots for all biomarkers can be found in Supplementary Fig. S1. All analyses were adjusted for decade of birth (categorical:1910 to <1920, 1920 to <1930, 1930 to <1940, 1940 to <1950, and 1950 to <1960), sex (male/female), self-reported race/ethnicity (African American, Japanese American, Latino, Native Hawaiian, and White), BMI (kg/m2, log), smoking history (pack-years of smoking), and urinary creatinine (mg/dL, log; Model 1). Model 1 was further adjusted for urinary TNE (Model 2) to assess the independent effects of internal smoking dose beyond that provided by self-reported measures of smoking. For TNE, Model 2 was adjusted for the ratio of urinary total 3-HCOT/cotinine. The proportional hazards assumption underlying Cox regression were checked using the “proportionality test” for the SAS procedure PHREG which tests proportionality among groups by creating interaction terms between log-time and covariates. Assumptions were met for all the variables of interest. Associations by tumor histologic cell-types were conducted when n > 20 and specific race/ethnicity were explored using Model 2 but we note the small sample size in these analyses (n's < 47). For associations by histologic cell-type, we performed a competing risk analysis using cause-specific models for time to lung cancer histologic cell-type outcomes, with censoring at diagnosis for any lung cancer case with a histologic cell-type other than that being considered. To compare the parameters by histologic cell-type, an augmented data approach as described in Lunn and McNeil was implemented that computes simultaneous models for lung cancer of each histologic cell-type (34). Heterogeneity across lung cancer histologic cell-type is assessed by a Wald test comparing the interaction by cell-type and smoking metabolites using robust variance estimates. Lung cancer incidence rates (overall and by race/ethnicity) were left truncated at age 45 and age-standardized using the United States 2000 standard population. All statistical analyses were performed using SAS version 9.4 (Statistical Analysis System, RRID:SCR_008567; http://www.sas.com). Supplementary Figure S2 was generated using Python (IPython, RRID:SCR_001658; http://ipython.org).

Data availability

Restrictions apply to the availability of these data. Data were obtained via an approved proposal by the MEC Study Research Committee. Data requests should be made to the MEC Study (see “Data Sharing” on the MEC website: https://www.uhcancercenter.org/mec). Investigators need to submit a formal application that will be evaluated internally by the MEC Research Committee before any data are released. Documentation of IRB approval is required for all projects requesting to use MEC data.

Baseline characteristics of this cohort of MEC smokers (N = 2,309) and the 140 incident lung cancer cases are presented in Table 1. Overall, the eligible population was comprised of men (46%) and women (54%) from five racial/ethnic groups: African Americans (16%), Native Hawaiians (14%), Whites (19%), Latinos (20%), and Japanese Americans (31%). Participants’ median age was 63 years, and their smoking history was a reported median of 23 pack-years. After an average of 13.4 years of follow-up from biospecimen collection, 140 incident lung cancer cases were identified. Cases were more likely to be male (55%) and there was a slightly greater population of African Americans (21%) and Native Hawaiians (22%) compared with the eligible population. Consistent with the analysis in the overall MEC study (35), age-standardized incidence rates (ASIR) were highest in Native Hawaiians and African Americans, followed by Whites, Japanese Americans, and Latinos (Supplementary Table S2). Here, ASIRs were estimated only among this smoking population (N = 2,309). Cases compared to the MEC eligible cohort were slightly leaner (median BMI:24.4 vs. 25.7 kg/m2) and reported a greater number of pack-years (median:34.3 vs. 23.3). The majority of lung cancer cases were diagnosed with adenocarcinoma (34%), followed by squamous cell carcinoma (29%), other/unspecified/nonspecific (NOS; 17%), small-cell lung cancer (16%), and large-cell carcinoma (4%).

Table 1.

Characteristics of the MEC smokers and the lung cancer cases identified during study follow-up.

MEC eligible populationIncident lung cancer cases
(N = 2,309)(N = 140)
n (%)n (%)
Sex 
 Males 1,068 (46) 77 (55) 
 Females 1,241 (54) 63 (45) 
Race/ethnicity 
 African Americans 368 (16) 29 (21) 
 Native Hawaiians 331 (14) 31 (22) 
 Whites 445 (19) 26 (19) 
 Latinos 456 (20) 15 (11) 
 Japanese Americans 709 (31) 39 (28) 
Age at urine collection 
 45–<55 years 105 (5) 2 (1) 
 55–<65 years 1,271 (55) 71 (51) 
 65–<75 years 706 (31) 48 (34) 
 ≥75 years 227 (10) 19 (14) 
MEC eligible populationIncident lung cancer cases
(N = 2,309)(N = 140)
n (%)n (%)
Sex 
 Males 1,068 (46) 77 (55) 
 Females 1,241 (54) 63 (45) 
Race/ethnicity 
 African Americans 368 (16) 29 (21) 
 Native Hawaiians 331 (14) 31 (22) 
 Whites 445 (19) 26 (19) 
 Latinos 456 (20) 15 (11) 
 Japanese Americans 709 (31) 39 (28) 
Age at urine collection 
 45–<55 years 105 (5) 2 (1) 
 55–<65 years 1,271 (55) 71 (51) 
 65–<75 years 706 (31) 48 (34) 
 ≥75 years 227 (10) 19 (14) 
Median (25–75% IQR)Median (25–75% IQR)
Follow-up time, years 13.4 (11.7–14.3) 8.4 (5.7–10.9) 
Body mass index (BMI) at time of urine collection (kg/m225.7 (22.7–28.8) 24.4 (22.1–27.3) 
Urinary creatinine, mg/dL 64 (40–103) 65 (33–107.5) 
Pack-years 23.3 (11.6–37.5) 34.3 (20.8–50.6) 
Median (25–75% IQR)Median (25–75% IQR)
Follow-up time, years 13.4 (11.7–14.3) 8.4 (5.7–10.9) 
Body mass index (BMI) at time of urine collection (kg/m225.7 (22.7–28.8) 24.4 (22.1–27.3) 
Urinary creatinine, mg/dL 64 (40–103) 65 (33–107.5) 
Pack-years 23.3 (11.6–37.5) 34.3 (20.8–50.6) 
n (%)
Histologic cell-type (case only) 
 Adenocarcinoma 47 (34) 
 Squamous cell carcinoma 40 (29) 
 Large cell or non-small cell, unspecified 6 (4) 
 Small-cell lung cancer 23 (16) 
 Other/unspecified/not specific 24 (17) 
n (%)
Histologic cell-type (case only) 
 Adenocarcinoma 47 (34) 
 Squamous cell carcinoma 40 (29) 
 Large cell or non-small cell, unspecified 6 (4) 
 Small-cell lung cancer 23 (16) 
 Other/unspecified/not specific 24 (17) 

Abbreviations: N, number of participants in each category; IQR, interquartile range.

The adjusted geometric means and 95% CIs of all urinary biomarkers evaluated in this analysis are presented in Table 2 (minimally adjusted for just race/ethnicity in Supplementary Table S3). Compared with the overall subcohort of current smokers at the time of urine collection, incident lung cancer cases had higher geometric mean levels of all urinary biomarkers, except for PheOH, which was the same across groups. All urinary biomarkers were positively correlated with TNE levels (Supplementary Fig. S2).

Table 2.

Geometric mean (GM) and 95% CI for urinary biomarkers of smoking-related toxicants among the MEC subcohort of current smokers and incident lung cancer cases.

MEC eligible populationIncident lung cancer cases
Urinary biomarkers of smokersNGM (95% CI)bNGM (95% CI)a
TNE (nmol/mL)b 2,309 30.8 (30.0–31.6) 140 38.6 (35.3–42.3) 
Total 3-HCOT/cotinine 2,307 3.20 (3.10–3.30) 140 3.85 (3.42–4.29) 
Total NNAL (pmol/mL) 2,251 1.15 (1.13–1.18) 140 1.42 (1.32–1.54) 
SPMA (pmol/mL) 2,169 2.52 (2.42–2.62) 132 2.97 (2.52–3.55) 
3-HPMA (nmol/mL) 2,281 3.00 (2.91–3.10) 140 3.93 (3.47–4.44) 
PheT (pmol/mL) 2,295 0.95 (0.93–0.98) 139 1.05 (0.95–1.16) 
PheOH (pmol/mL) 2,253 0.71 (0.69–0.73) 136 0.71 (0.65–0.77) 
PheT/PheOH 2,239 1.35 (1.32–1.39) 135 1.50 (1.33–1.68) 
Cd (ng/mL) 1,976 0.60 (0.58–0.61) 125 0.77 (0.70–0.84) 
8-iso-PGF (pmol/mL) 1,913 0.80 (0.78–0.82) 122 0.83 (0.71–0.92) 
MEC eligible populationIncident lung cancer cases
Urinary biomarkers of smokersNGM (95% CI)bNGM (95% CI)a
TNE (nmol/mL)b 2,309 30.8 (30.0–31.6) 140 38.6 (35.3–42.3) 
Total 3-HCOT/cotinine 2,307 3.20 (3.10–3.30) 140 3.85 (3.42–4.29) 
Total NNAL (pmol/mL) 2,251 1.15 (1.13–1.18) 140 1.42 (1.32–1.54) 
SPMA (pmol/mL) 2,169 2.52 (2.42–2.62) 132 2.97 (2.52–3.55) 
3-HPMA (nmol/mL) 2,281 3.00 (2.91–3.10) 140 3.93 (3.47–4.44) 
PheT (pmol/mL) 2,295 0.95 (0.93–0.98) 139 1.05 (0.95–1.16) 
PheOH (pmol/mL) 2,253 0.71 (0.69–0.73) 136 0.71 (0.65–0.77) 
PheT/PheOH 2,239 1.35 (1.32–1.39) 135 1.50 (1.33–1.68) 
Cd (ng/mL) 1,976 0.60 (0.58–0.61) 125 0.77 (0.70–0.84) 
8-iso-PGF (pmol/mL) 1,913 0.80 (0.78–0.82) 122 0.83 (0.71–0.92) 

Abbreviations: N, number of events; GM, geometric mean; CI, confidence intervals; TNE, total nicotine equivalents; 3-HCOT, trans 3′-hydroxycotinine; NNAL, 4-(methylnitrosamino)-1–3-pyridyl)-1-butanol; SPMA, S-phenylmercapturic acid; 3-HPMA, 3-hydroxypropyl mercapturic acid; PheT, phenanthrene tetraol; PheOH, 3-hydroxyphenanthrene; PheT/PheOH, a proposed biomarker of metabolic activation of polycyclic aromatic hydrocarbons (PAH); Cd, cadmium; 8-iso-PGF, (Z)-7-[1R,2R,3R,5S)-3,5-dihydroxy-2-[(E,3S)-3-hydroxyoct-1-enyl]cyclopentyl]hept-5-enoic acid.

aGeometric means and corresponding 95% CI are adjusted for age, sex (male/female), race/ethnicity (African American, Native Hawaiian, White, Latino, Japanese American), body mass index (BMI, kg/m2; log), creatinine (mg/dL; log), cigarettes per day (CPD), and TNE (nmol/mL; where appropriate).

bFor TNE, the model does not include TNE.

The associations of known risk factors with lung cancer risk were confirmed (Supplementary Table S4). The association between smoking-related toxicant and carcinogen biomarkers and lung cancer risk in this multiethnic population of current smokers is presented in Table 3. After adjusting for decade of birth, sex, race/ethnicity, BMI, and creatinine, a one-SD increase in log-pack-years was associated with a 90% increase in lung cancer risk in current smokers (Table 3, Model 1). After adjustment for birth, sex, race/ethnicity, BMI, and creatinine and pack-years of smoking, we found that both a one SD increase of log TNE and a one SD increase of log total 3-HCOT/cotinine ratio were associated with lung cancer risk (HR, 1.36; 95% CI, 1.00–1.84 and HR, 1.37; 95% CI:1.11–1.71, respectively; Table 3, Model 1). When TNE was further adjusted for total 3-HCOT/cotinine, the association was attenuated and no longer remained statistically significant (HR per SD increase in log-TNE = 1.22; 95% CI, 0.91–1.64; Table 3, Model 2). Whereas the ratio of total 3-HCOT/cotinine remained significantly associated with lung cancer when adjusted for TNE (HR per SD increase in log-total 3-HCOT/cotinine = 1.33; 95% CI, 1.06–1.66; Table 3, Model 2).

Table 3.

Association of pack-years and urinary biomarkers of smoking-related toxicants with lung cancer incidence in the MEC subcohort of current smokers at time of urine collection.

Model 1bModel 2c
Urinary biomarkers of smoking-related toxicantsaNHR (95% CI)PHR (95% CI)P
Pack-yearsd 140 1.90 (1.50–2.41) <0.0001 1.70 (1.31–2.21) <0.0001 
TNEe 140 1.36 (1.00–1.84) 0.047 1.22 (0.91–1.64) 0.192 
Total 3-HCOT/cotinine 140 1.37 (1.11–1.71) 0.004 1.33 (1.06–1.66) 0.014 
Total NNAL 140 1.20 (0.95–1.52) 0.120 1.10 (0.82–1.47) 0.541 
SPMA 132 1.20 (0.96–1.49) 0.104 1.12 (0.89–1.41) 0.325 
3-HPMA 138 1.46 (1.14–1.88) 0.003 1.41 (1.07–1.85) 0.016 
PheT 139 1.08 (0.86–1.35) 0.500 0.98 (0.76–1.25) 0.842 
PheOH 136 1.00 (0.80–1.27) 0.976 0.87 (0.67–1.13) 0.308 
PheT/PheOH 135 1.05 (0.88–1.26) 0.579 1.03 (0.86–1.23) 0.758 
Cd 125 1.48 (1.21–1.82) 0.0002 1.45 (1.18–1.79) 0.0004 
8-iso-PGF 122 1.17 (0.91–1.49) 0.217 1.14 (0.89–1.46) 0.305 
Model 1bModel 2c
Urinary biomarkers of smoking-related toxicantsaNHR (95% CI)PHR (95% CI)P
Pack-yearsd 140 1.90 (1.50–2.41) <0.0001 1.70 (1.31–2.21) <0.0001 
TNEe 140 1.36 (1.00–1.84) 0.047 1.22 (0.91–1.64) 0.192 
Total 3-HCOT/cotinine 140 1.37 (1.11–1.71) 0.004 1.33 (1.06–1.66) 0.014 
Total NNAL 140 1.20 (0.95–1.52) 0.120 1.10 (0.82–1.47) 0.541 
SPMA 132 1.20 (0.96–1.49) 0.104 1.12 (0.89–1.41) 0.325 
3-HPMA 138 1.46 (1.14–1.88) 0.003 1.41 (1.07–1.85) 0.016 
PheT 139 1.08 (0.86–1.35) 0.500 0.98 (0.76–1.25) 0.842 
PheOH 136 1.00 (0.80–1.27) 0.976 0.87 (0.67–1.13) 0.308 
PheT/PheOH 135 1.05 (0.88–1.26) 0.579 1.03 (0.86–1.23) 0.758 
Cd 125 1.48 (1.21–1.82) 0.0002 1.45 (1.18–1.79) 0.0004 
8-iso-PGF 122 1.17 (0.91–1.49) 0.217 1.14 (0.89–1.46) 0.305 

Abbreviations: N, number of events; pack-years, Number of packs of cigarettes smoked per day × number of years the person has smoked.

aPack-years and all urinary biomarkers were standardized using log transformation and dividing the individual value by the overall population SD of the log biomarker and therefore the HR corresponds to a per one-unit SD change in log biomarker level. The SD of the log-biomarker can be found in Supplementary Table S1.

bModel 1 (base model): adjusted for decade of birth, sex (male/female), race/ethnicity (African American, Native Hawaiian, White, Latino, Japanese American), body mass index (BMI, kg/m2; log), creatinine (mg/dL; log), and pack-years of smoking.

cModel 2: Model 1 + TNE.

dFor pack-years, Model 1 adjusted for decade of birth, sex, race/ethnicity, BMI (kg/m2; log), creatinine (mg/dL; log); Model 2 additionally adjusted for TNE.

eFor TNE, Model 1 adjusted for decade of birth, sex, race/ethnicity, BMI (kg/m2; log), creatinine (mg/dL; log), and smoking history (pack-years); Model 2 additional adjusts for total 3-HCOT/cotinine.

Smoking-related toxicant levels of urinary 3-HPMA and Cd were individually associated with lung cancer risk after adjusting for decade of birth, sex, race/ethnicity, BMI, urinary creatinine, and self-reported pack-years (HR per SD increase in log-3-HPMA, 1.46; 95% CI, 1.14–1.88 and HR per SD increase in log-Cd, 1.48; 95% CI, 1.21–1.82; Table 3, Model 1). These associations remained even after adjustment for TNE (HR per SD increase in log-3-HPMA = 1.41; 95% CI, 1.07–1.85 and HR per SD increase in log-Cd, 1.45; 95% CI, 1.18–1.79). 3-HPMA association was driven primarily by those who smoked less intensely (TNE < median 32.4 nmol/mL HR, 1.78; 95% CI, 1.17–2.69 vs. TNE ≥ median 32.4 nmol/mL HR, 1.16; 95% CI, 0.80–1.70; Supplementary Table S5) and in a population who smoked less on average (Japanese Americans, n = 29; HR, 1.96; 95% CI, 1.14–3.36; Pinteraction race/ethnicity = 0.029; Table 4). The association for Cd was primarily driven among those who smoked more intensely (TNE < median 32.4 nmol/mL HR, 1.42; 95% CI, 1.07–1.89 vs. TNE ≥ median 32.4 nmol/mL HR, 1.44; 95% CI, 1.01–2.04; Supplementary Table S5) and by the results in Whites (n = 23; HR Cd, 1.91; 95% CI, 1.12–3.24; Pinteraction = 0.383).

Table 4.

Association of pack-years and urinary biomarkers of smoking-related toxicants with lung cancer incidence by race/ethnicity—MEC.

African AmericansNative HawaiiansWhitesLatinosJapanese Americans
Urinary biomarkers of smoking-related toxicantsanHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPPinteractionc
Pack-yearsd 29 1.17 (0.75–1.84) 0.494 31 1.79 (1.06–3.03) 0.029 26 4.19 (2.07–8.48) <0.0001 15 1.54 (0.58–4.09) 0.387 39 1.60 (0.93–2.75) 0.087 0.011 
TNEe 29 0.86 (0.51–1.44) 0.564 31 0.98 (0.51–1.88) 0.955 26 1.81 (0.77–4.29) 0.177 15 1.52 (0.57–4.10) 0.407 39 1.71 (0.91–3.21) 0.096 0.072 
Total 3-HCOT/ cotinine 29 1.42 (0.84–2.40) 0.195 31 1.21 (0.75–1.95) 0.437 26 1.99 (1.03–3.81) 0.039 15 1.05 (0.47–2.34) 0.899 39 1.19 (0.81–1.75) 0.382 0.318 
Total NNAL 29 0.81 (0.44–1.51) 0.511 31 1.02 (0.52–2.00) 0.946 26 0.77 (0.33–1.78) 0.539 15 1.15 (0.39–3.38) 0.799 39 1.45 (0.84–2.51) 0.187 0.044 
SPMA 28 1.41 (0.90–2.23) 0.138 30 0.84 (0.52–1.36) 0.478 25 1.19 (0.66–2.14) 0.555 15 0.64 (0.35–1.17) 0.145 34 1.27 (0.83–1.95) 0.276 0.318 
3-HPMA 28 0.96 (0.57–1.62) 0.882 30 1.26 (0.62–2.57) 0.521 26 0.95 (0.47–1.91) 0.877 15 2.02 (0.85–4.80) 0.110 39 1.96 (1.14–3.36) 0.015 0.029 
PheT 29 1.59 (1.11–2.28) 0.011 30 0.72 (0.37–1.37) 0.315 26 0.75 (0.38–1.49) 0.407 15 0.49 (0.20–1.20) 0.117 39 0.78 (0.44–1.36) 0.376 0.056 
PheOH 27 1.03 (0.66–1.60) 0.912 29 0.87 (0.46–1.65) 0.675 26 0.68 (0.31–1.52) 0.346 15 0.29 (0.11–0.76) 0.012 39 0.99 (0.57–1.71) 0.963 0.461 
PheT/PheOH 27 0.29 (0.11–0.76) 0.012 28 0.87 (0.57–1.34) 0.529 26 0.29 (0.11–0.76) 0.012 15 1.12 (0.65–1.94) 0.688 39 0.84 (0.57–1.24) 0.375 0.261 
Cd 26 1.72 (0.96–3.07) 0.067 29 1.29 (0.88–1.89) 0.197 22 1.91 (1.12–3.24) 0.017 14 1.38 (0.56–3.43) 0.482 34 1.05 (0.67–1.66) 0.821 0.383 
8-iso-PGF 22 1.26 (0.70–2.28) 0.435 29 0.90 (0.58–1.42) 0.660 23 1.38 (0.68–2.80) 0.366 14 1.28 (0.52–3.16) 0.593 34 1.06 (0.64–1.76) 0.819 0.429 
African AmericansNative HawaiiansWhitesLatinosJapanese Americans
Urinary biomarkers of smoking-related toxicantsanHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPPinteractionc
Pack-yearsd 29 1.17 (0.75–1.84) 0.494 31 1.79 (1.06–3.03) 0.029 26 4.19 (2.07–8.48) <0.0001 15 1.54 (0.58–4.09) 0.387 39 1.60 (0.93–2.75) 0.087 0.011 
TNEe 29 0.86 (0.51–1.44) 0.564 31 0.98 (0.51–1.88) 0.955 26 1.81 (0.77–4.29) 0.177 15 1.52 (0.57–4.10) 0.407 39 1.71 (0.91–3.21) 0.096 0.072 
Total 3-HCOT/ cotinine 29 1.42 (0.84–2.40) 0.195 31 1.21 (0.75–1.95) 0.437 26 1.99 (1.03–3.81) 0.039 15 1.05 (0.47–2.34) 0.899 39 1.19 (0.81–1.75) 0.382 0.318 
Total NNAL 29 0.81 (0.44–1.51) 0.511 31 1.02 (0.52–2.00) 0.946 26 0.77 (0.33–1.78) 0.539 15 1.15 (0.39–3.38) 0.799 39 1.45 (0.84–2.51) 0.187 0.044 
SPMA 28 1.41 (0.90–2.23) 0.138 30 0.84 (0.52–1.36) 0.478 25 1.19 (0.66–2.14) 0.555 15 0.64 (0.35–1.17) 0.145 34 1.27 (0.83–1.95) 0.276 0.318 
3-HPMA 28 0.96 (0.57–1.62) 0.882 30 1.26 (0.62–2.57) 0.521 26 0.95 (0.47–1.91) 0.877 15 2.02 (0.85–4.80) 0.110 39 1.96 (1.14–3.36) 0.015 0.029 
PheT 29 1.59 (1.11–2.28) 0.011 30 0.72 (0.37–1.37) 0.315 26 0.75 (0.38–1.49) 0.407 15 0.49 (0.20–1.20) 0.117 39 0.78 (0.44–1.36) 0.376 0.056 
PheOH 27 1.03 (0.66–1.60) 0.912 29 0.87 (0.46–1.65) 0.675 26 0.68 (0.31–1.52) 0.346 15 0.29 (0.11–0.76) 0.012 39 0.99 (0.57–1.71) 0.963 0.461 
PheT/PheOH 27 0.29 (0.11–0.76) 0.012 28 0.87 (0.57–1.34) 0.529 26 0.29 (0.11–0.76) 0.012 15 1.12 (0.65–1.94) 0.688 39 0.84 (0.57–1.24) 0.375 0.261 
Cd 26 1.72 (0.96–3.07) 0.067 29 1.29 (0.88–1.89) 0.197 22 1.91 (1.12–3.24) 0.017 14 1.38 (0.56–3.43) 0.482 34 1.05 (0.67–1.66) 0.821 0.383 
8-iso-PGF 22 1.26 (0.70–2.28) 0.435 29 0.90 (0.58–1.42) 0.660 23 1.38 (0.68–2.80) 0.366 14 1.28 (0.52–3.16) 0.593 34 1.06 (0.64–1.76) 0.819 0.429 

Abbreviations: N, number of events; pack-years, number of packs of cigarettes smoked per day × number of years the person has smoked.

aPack-years and all urinary biomarkers were standardized using log transformation and dividing the individual value by the overall population SD of the log biomarker and therefore the HR corresponds to a per one-unit SD change in log biomarker level. The SD of the log-biomarker can be found in Supplementary Table S1.

bModel adjusted for decade of birth, sex (male/female), body mass index (BMI, kg/m2; log), creatinine (mg/dL, log), pack-years (where appropriate), and TNE (nmol/mL; where appropriate).

cP value across race/ethnicity.

dFor pack-years, the model does not adjust for pack-years.

eFor TNE, the model does not include TNE and is additional adjusted for total 3-HCOT/cotinine.

We also explored the association by histologic cell-type (Table 5) and found that the association between incident lung cancer risk and pack-years was strongest for adenocarcinoma (ADC; P < 0.0001; Table 5). The association with TNE was strongest for risk of SCC (P = 0.038). The association with 3-HPMA was strongest for risk of lung cancer of other/unspecified histologic cell types (P = 0.008). The association with urinary Cd was strongest for risk of ADC (P = 0.001). These associations were not found heterogeneous across histologic cell types. Associations of race/ethnicity and by histologic cell-type were not conducted due to small sample size (Supplementary Table S6).

Table 5.

Association of pack-years and urinary biomarkers of smoking-related toxicants with lung cancer incidence – the MEC subcohort of current smokers at time of urine collection, stratified by histologic cell-type.

ADCSCCSCLCUnspecified
Urinary biomarkers of smoking-related toxicantsanHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPPheterogeneityc
Pack-yearsd 47 2.15 (1.41–3.27) 0.0004 40 1.65 (1.09–2.51) 0.019 23 1.87 (1.04–3.36) 0.037 27 1.96 (1.13–3.37) 0.016 0.841 
TNEe 47 0.94 (0.58–1.53) 0.806 40 1.82 (1.03–3.20) 0.038 23 1.27 (0.62–2.61) 0.511 27 1.03 (0.54–1.97) 0.931 0.365 
Total 3-HCOT/cotinine 47 1.12 (0.77–1.64) 0.559 40 1.41 (0.91–2.16) 0.121 23 1.72 (0.95–3.09) 0.071 27 1.30 (0.80–2.13) 0.293 0.652 
Total NNAL 47 1.38 (0.92–2.08) 0.120 40 1.05 (0.61–1.80) 0.869 23 1.08 (0.50–2.35) 0.844 27 0.63 (0.30–1.34) 0.233 0.172 
SPMA 46 1.06 (0.72–1.56) 0.766 36 1.37 (0.88–2.14) 0.166 22 1.04 (0.60–1.81) 0.891 25 0.94 (0.55–1.60) 0.812 0.587 
3-HPMA 46 1.27 (0.80–2.04) 0.312 40 1.28 (0.78–2.12) 0.333 22 1.03 (0.50–2.14) 0.930 27 2.21 (1.21–4.04) 0.010 0.762 
PheT 46 0.86 (0.55–1.34) 0.510 40 1.18 (0.75–1.84) 0.475 23 0.81 (0.42–1.55) 0.522 27 0.90 (0.50–1.62) 0.718 0.501 
PheOH 45 0.77 (0.48–1.25) 0.293 39 0.85 (0.52–1.38) 0.516 22 1.05 (0.56–1.96) 0.883 27 0.97 (0.56–1.67) 0.910 0.418 
PheT/PheOH 44 0.97 (0.71–1.33) 0.859 39 1.22 (0.88–1.68) 0.235 22 0.81 (0.51–1.29) 0.377 27 0.95 (0.65–1.39) 0.791 0.544 
Cd 42 1.75 (1.25–2.46) 0.001 38 0.96 (0.62–1.49) 0.870 21 1.54 (0.92–2.57) 0.101 22 1.64 (1.05–2.56) 0.030 0.338 
8-iso-PGF 40 1.42 (0.92–2.21) 0.115 37 0.99 (0.63–1.54) 0.948 19 1.49 (0.76–2.92) 0.243 24 0.86 (0.51–1.43) 0.557 0.371 
ADCSCCSCLCUnspecified
Urinary biomarkers of smoking-related toxicantsanHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPnHR (95% CI)bPPheterogeneityc
Pack-yearsd 47 2.15 (1.41–3.27) 0.0004 40 1.65 (1.09–2.51) 0.019 23 1.87 (1.04–3.36) 0.037 27 1.96 (1.13–3.37) 0.016 0.841 
TNEe 47 0.94 (0.58–1.53) 0.806 40 1.82 (1.03–3.20) 0.038 23 1.27 (0.62–2.61) 0.511 27 1.03 (0.54–1.97) 0.931 0.365 
Total 3-HCOT/cotinine 47 1.12 (0.77–1.64) 0.559 40 1.41 (0.91–2.16) 0.121 23 1.72 (0.95–3.09) 0.071 27 1.30 (0.80–2.13) 0.293 0.652 
Total NNAL 47 1.38 (0.92–2.08) 0.120 40 1.05 (0.61–1.80) 0.869 23 1.08 (0.50–2.35) 0.844 27 0.63 (0.30–1.34) 0.233 0.172 
SPMA 46 1.06 (0.72–1.56) 0.766 36 1.37 (0.88–2.14) 0.166 22 1.04 (0.60–1.81) 0.891 25 0.94 (0.55–1.60) 0.812 0.587 
3-HPMA 46 1.27 (0.80–2.04) 0.312 40 1.28 (0.78–2.12) 0.333 22 1.03 (0.50–2.14) 0.930 27 2.21 (1.21–4.04) 0.010 0.762 
PheT 46 0.86 (0.55–1.34) 0.510 40 1.18 (0.75–1.84) 0.475 23 0.81 (0.42–1.55) 0.522 27 0.90 (0.50–1.62) 0.718 0.501 
PheOH 45 0.77 (0.48–1.25) 0.293 39 0.85 (0.52–1.38) 0.516 22 1.05 (0.56–1.96) 0.883 27 0.97 (0.56–1.67) 0.910 0.418 
PheT/PheOH 44 0.97 (0.71–1.33) 0.859 39 1.22 (0.88–1.68) 0.235 22 0.81 (0.51–1.29) 0.377 27 0.95 (0.65–1.39) 0.791 0.544 
Cd 42 1.75 (1.25–2.46) 0.001 38 0.96 (0.62–1.49) 0.870 21 1.54 (0.92–2.57) 0.101 22 1.64 (1.05–2.56) 0.030 0.338 
8-iso-PGF 40 1.42 (0.92–2.21) 0.115 37 0.99 (0.63–1.54) 0.948 19 1.49 (0.76–2.92) 0.243 24 0.86 (0.51–1.43) 0.557 0.371 

Abbreviations: ADC, adenocarcinoma; SCC, squamous cell carcinoma; SCLC, small-cell lung cancer; unspecified, unspecified malignant neoplasm; N, number of events; pack-years, number of packs of cigarettes smoked per day x number of years the person has smoked.

aPack-years and all urinary biomarkers were standardized using log transformation and dividing the individual value by the overall population SD of the log biomarker and therefore the HR corresponds to a per one-unit SD change in log biomarker level. The SD of the log-biomarker can be found in Supplementary Table S1.

bModel adjusted for decade of birth, sex (male/female), race/ethnicity (African American, Native Hawaiian, White, Latino, Japanese American), body mass index (BMI, kg/m2; log), creatinine (mg/dL, log), pack-years (where appropriate), and TNE (nmol/mL; where appropriate).

cPheterogeneity across histologic cell-type in a competing risk model.

dFor pack-years, the model does not adjust for pack-years.

eFor TNE, the model does not include TNE and is additional adjusted for total 3-HCOT/cotinine.

This is the first study to examine the individual and joint effects of urinary biomarkers of multiple classes of smoking-related toxicants with the risk of smoking-related lung cancer in a multiethnic population. We demonstrated that urinary biomarkers of total 3-HCOT/cotinine, 3-HPMA (a metabolite of acrolein), and Cd (a biomarker of long-term Cd exposure) were significantly associated with lung cancer risk, independent of self-reported smoking history (pack-years) and internal smoking dose (TNE).

CYP2A6 catalyzed 5′-oxidation, the primary pathway of nicotine metabolism in most smokers produces cotinine (21), the most commonly used biomarker of tobacco exposure (21, 36). However, cotinine levels are influenced by individual variability in nicotine metabolism, which is known to influence exposure to tobacco toxicants (12, 18–21, 37). TNE, which includes metabolites from three pathways of nicotine metabolism accounts for approximately 85% of the nicotine dose and better reflects nicotine uptake, internal smoking dose, and lung cancer risk in some populations (4, 21). Because CYP2A6 also catalyzes the oxidation of cotinine to 3-HCOT, the ratio of 3-HCOT/cotinine is a phenotypic measure of CYP2A6 enzymatic activity. Variation in CYP2A6 activity can influence smoking intensity (30, 37, 38) and the ability to quit smoking over time (39) and therefore the urinary ratio of 3HCOT/cotinine may reflect longer-term smoking behaviors that are not captured by the short-term biomarker of dose (e.g., TNE).

Previously, we reported, in the same population studied here that a one log unit increase in the ratio of 3-HCOT/cotinine was associated with a 52% increase in lung cancer risk (n = 92 cases), independent of CPD, self-reported smoking duration, and TNE (16). In the current analysis with an additional 48 cases, the ratio of total 3-HCOT/cotinine was also significantly associated with lung cancer risk after adjustment for similar covariates (this analysis adjusted for the measure pack-years instead of CPD and self-reported smoking duration). These data support our hypothesis that this biomarker provides additional information regarding smoking history that may not be captured by self-reported data (16). In the Singapore Chinese Health Study, a positive association between total 3-HCOT/cotinine and lung cancer was observed; however, after adjustment for TNE, the association was attenuated (13). It is not surprising that in different populations with different genetics of nicotine metabolism and tobacco use, the ability of TNE and/or the ratio of total 3-HCOT/cotinine to reflect long- term smoking behavior varies.

Acrolein, an α,β-unsaturated aldehyde, is a ubiquitous environmental and dietary constituent and is produced endogenously as a product of lipid peroxidation, amino acid metabolism, and polyamine metabolism (40). Among current smokers, active tobacco smoking is the major source of exposure to acrolein (23, 41). One recent analysis reported a range of acrolein levels from 30.8–82.6 μg per unit in the mainstream smoke [under International Organization for Standardization (ISO) conditions] of 35 commercial cigarette tobacco products sold in the United States (42). Acrolein is a known lung irritant that produces inflammation and various other effects involved in carcinogenesis. As such, in 2020, IARC reclassified acrolein into Group 2A as “probably carcinogenic to humans” (23). Aligned with this reclassification, we found that 3-HPMA (the major metabolite of acrolein) was associated with lung cancer risk in our MEC smoker population after adjusting for lung cancer risk factors, the short-term biomarker of dose (TNE), and common urinary biomarkers of smoking. This association was primarily driven by those that smoked less intensely and populations who smoked less on average (e.g., Japanese Americans and Latinos). The Shanghai Cohort Study of male current smokers (n = 343 lung cancer cases; ref. 43) and never-smokers only (n = 82 lung cancer cases; ref. 17) and the Golestan Cohort Study in Iran of male exclusive cigarette smokers (n = 31 lung cancer cases; ref. 44) found an association with acrolein (measured by 3-HPMA) and lung cancer risk. However, these studies found that the association was no longer statistically significant after adjusting for total cotinine or smoking history, respectively, suggesting that the observed effect between 3-HPMA and lung cancer risk partly reflected an added measure of smoking dose.

Cd is a known human carcinogen (IARC group 1), and the primary source of exposure in most cases (e.g., leaving aside occupational exposure) is cigarette smoke (45, 46). While toxicologic studies in rodents and occupational studies in humans have shown that Cd is a respiratory toxicant associated with lung cancer, occupational studies cannot be extrapolated to the general population, and few studies have properly adjusted for smoking in their risk estimates (45, 47–49). Three non-occupationally based prospective cohort studies in the U.S. (NHANES III and Strong Heart Study) and Belgium (CadmiBel Study) have shown a positive association between urinary Cd levels, a biomarker of long-term Cd exposure, and lung cancer risk after adjusting for smoking status and/or pack-years (50–52). In two of these studies, after removing current smokers at baseline (Strong Heart Study; ref. 52) or including only never-smokers (U.S. NHANES III; ref. 51) in the analysis, the association remained positive but weaker. In the current study, we further demonstrated a positive association between Cd and lung cancer incidence, independent of measures of self-reported smoking history. Interestingly, when we explored associations by histologic cell-type, we demonstrated a strong positive association between urinary Cd and adenocarcinoma in our population, which is in agreement with previous evidence that has shown chronic Cd inhalation in rodents causes pulmonary adenocarcinomas (53). These results, in combination with those of other studies, suggest that urinary Cd may provide additional information regarding smoking-related lung cancer risk.

Prior studies have also detected associations with NNAL (8–10), PheT (10, 11, 17), PheOH (17), SPMA (43), and 8-iso-PGF (54, 55) with lung cancer risk, even after accounting for self-reported smoking and/or smoking dose (cotinine and/or TNE). In this study, we did not detect an association between these biomarkers and the risk of lung cancer with adjustment for pack-years and TNE. These reported differences across studies may be a result of differences in sample size (MEC with only 140 lung cancer cases), smoking behaviors, and study population (e.g., the Shanghai population smoked mainly Chinese cigarettes in the 1980s vs. the multiethnic U.S. population smoked mainly domestic cigarettes in the 2000s). While all lung cancer histologic cell types are attributed to tobacco smoking, the cigarette quantity and composition have been shown to influence lung cancer subtypes. For instance, cigarette smoking dose was associated with higher relative risks for squamous cell carcinoma and small-cell lung cancer than adenocarcinoma (56). In addition, increased filtration ventilation of cigarettes has been suggested as a potential contributor to the relative rise of adenocarcinoma (57). Our exploratory analyses stratified by race/ethnicity and by histologic cell-type further supports the possibility that variation in exposure to smoke toxicants may contribute to etiologic differences in lung cancer risk.

The main strength of this study is the use of a well-characterized multiethnic population of current smokers with a range of smoking intensity and our ability to evaluate the effects of nine urinary biomarkers of smoking-related toxicants with lung cancer risk independent of self-reported smoking history and urinary TNE (a biomarker for internal smoking dose). However, our study has some limitations. First, this study included a modest number of lung cancer cases, limiting the power of our racial/ethnic and histologic cell-type–specific analyses. Second, the biomarker levels were only measured from a single urine collection; therefore, we did not account for variations in smoking behavior and exposure to tobacco carcinogens over time. However, adult smoking behavior is relatively stable over time (58). Third, we did not have information on quitting attempts or quitting during follow-up or information on cigarette brands, which could contribute to differences in biomarker levels of smoking-related toxicants.

In conclusion, our findings suggest that urinary total 3-HCOT/cotinine, 3-HPMA, and Cd levels may provide additional information on smoking-related lung cancer risk that is not captured by self-reported smoking history or the short-term biomarker of internal smoking dose (TNE). Our findings also suggest that these biomarkers may provide distinct information for lung cancer risk prediction by population or histologic cell-type. Replication in large prospective studies is warranted.

S.S. Cigan reports grants from National Institute for Occupational Safety and Health (NIOSH), Centers for Disease Control and Prevention and grants from the National Institute of Health (NIH)/National Cancer Institute (NCI) during the conduct of the study. D.O. Stram reports grants from NIH during the conduct of the study; personal fees from US Environmental Protection Agency outside the submitted work. S.S. Hecht reports grants from National Cancer Institute during the conduct of the study. L. Le Marchand reports grants from NCI during the conduct of the study. I. Stepanov reports grants from NIH/NCI during the conduct of the study. S.L. Park reports grants from NIH/NCI during the conduct of the study. No disclosures were reported by the other authors.

The funders had no role in the study design, data collection, analysis, decision to publish, or preparation of the manuscript.

S.S. Cigan: Conceptualization, formal analysis, methodology, writing–original draft, writing–review and editing. S.E. Murphy: Conceptualization, data curation, funding acquisition, investigation, writing–review and editing. D.O. Stram: Funding acquisition, methodology, writing–review and editing. S.S. Hecht: Conceptualization, resources, data curation, funding acquisition, investigation, project administration, writing–review and editing. L. Le Marchand: Conceptualization, resources, data curation, supervision, funding acquisition, investigation, project administration, writing–review and editing. I. Stepanov: Conceptualization, resources, funding acquisition, writing–review and editing. S.L. Park: Conceptualization, data curation, formal analysis, supervision, funding acquisition, investigation, methodology, writing–original draft, project administration, writing–review and editing.

The authors gratefully acknowledge the time and efforts of all the MEC study participants. The MEC Study is supported by the NIH/NCI, grant numbers: P01 CA138338 (University of Minnesota, to S. Hecht), P30 CA071789 (University of Hawaii at Manoa, to L. Le Marchand), and U01 CA164973 (University of Hawaii Cancer Center to MPIs: L. Le Marchand, C. Haiman, L.R. Wilkens). S.S. Cigan was additionally supported in part by grants from the National Institute for Occupational Safety and Health (NIOSH), Centers for Disease Control and Prevention (5T42 OH008434, University of Minnesota, to S. Gerberich) and NCI (R01 CA179246, University of Minnesota, to I. Stepanov).

The publication costs of this article were defrayed in part by the payment of publication fees. Therefore, and solely to indicate this fact, this article is hereby marked “advertisement” in accordance with 18 USC section 1734.

Note: Supplementary data for this article are available at Cancer Epidemiology, Biomarkers & Prevention Online (http://cebp.aacrjournals.org/).

1.
Siegel
RL
,
Miller
KD
,
Fuchs
HE
,
Jemal
A
.
Cancer statistics, 2022
.
CA Cancer J Clinicians
2022
;
72
:
7
33
.
2.
International Agency for Research on Cancer.
Tobacco smoke and involuntary smoking
.
IARC Monographs on the Evaluation of Carcinogenic Risks to Humans
.
Lyon, France: IARC
;
2004
.
3.
Li
Y
,
Hecht
SS
.
Carcinogenic components of tobacco and tobacco smoke: A 2022 update
.
Food Chem Toxicol
2022
;
165
:
113179
.
4.
Murphy
SE
.
Biochemistry of nicotine metabolism and its relevance to lung cancer
.
J Biol Chem
2021
;
296
:
100722
.
5.
Hecht
SS
,
Hatsukami
DK
.
Smokeless tobacco and cigarette smoking: chemical mechanisms and cancer prevention
.
Nat Rev Cancer
2022
;
22
:
143
55
.
6.
de Waard
F
,
Kemmeren
JM
,
van Ginkel
LA
,
Stolker
AAM
.
Urinary cotinine and lung cancer risk in a female cohort
.
Br J Cancer
1995
;
72
:
784
7
.
7.
Boffetta
P
,
Clark
S
,
Shen
M
,
Gislefoss
R
,
Peto
R
,
Andersen
A
.
Serum cotinine level as predictor of lung cancer risk
.
Cancer Epidemiol Biomarkers Prev
2006
;
15
:
1184
8
.
8.
Yuan
JM
,
Koh
WP
,
Murphy
SE
,
Fan
Y
,
Wang
R
,
Carmella
SG
, et al
.
Urinary levels of tobacco-specific nitrosamine metabolites in relation to lung cancer development in two prospective cohorts of cigarette smokers
.
Cancer Res
2009
;
69
:
2990
5
.
9.
Church
TR
,
Anderson
KE
,
Caporaso
NE
,
Geisser
MS
,
Le
CT
,
Zhang
Y
, et al
.
A prospectively measured serum biomarker for a tobacco-specific carcinogen and lung cancer in smokers
.
Cancer Epidemiol Biomarkers Prev
2009
;
18
:
260
6
.
10.
Yuan
JM
,
Gao
YT
,
Murphy
SE
,
Carmella
SG
,
Wang
R
,
Zhong
Y
, et al
.
Urinary levels of cigarette smoke constituent metabolites are prospectively associated with lung cancer development in smokers
.
Cancer Res
2011
;
71
:
6749
57
.
11.
Hecht
SS
,
Murphy
SE
,
Stepanov
I
,
Nelson
HH
,
Yuan
JMM
.
Tobacco smoke biomarkers and cancer risk among male smokers in the Shanghai cohort study
.
Cancer Lett
2013
;
334
:
34
8
.
12.
Yuan
JM
,
Nelson
HH
,
Butler
LM
,
Carmella
SG
,
Wang
R
,
Kuriger-Laber
JK
, et al
.
Genetic determinants of cytochrome P450 2A6 activity and biomarkers of tobacco smoke exposure in relation to risk of lung cancer development in the Shanghai cohort study
.
Int J Cancer
2016
;
138
:
2161
71
.
13.
Yuan
JM
,
Nelson
HH
,
Carmella
SG
,
Wang
R
,
Kuriger-Laber
J
,
Jin
A
, et al
.
CYP2A6 genetic polymorphisms and biomarkers of tobacco smoke constituents in relation to risk of lung cancer in the Singapore Chinese Health Study
.
Carcinogenesis
2017
;
38
:
411
8
.
14.
Larose
TL
,
Guida
F
,
Fanidi
A
,
Langhammer
A
,
Kveem
K
,
Stevens
VL
, et al
.
Circulating cotinine concentrations and lung cancer risk in the Lung Cancer Cohort Consortium (LC3)
.
Int J Epidemiol
2018
;
47
:
1760
71
.
15.
Thomas
CE
,
Wang
R
,
Adams-Haduch
J
,
Murphy
SE
,
Ueland
PM
,
Midttun
Ø
, et al
.
Urinary cotinine is as good a biomarker as serum cotinine for cigarette smoking exposure and lung cancer risk prediction
.
Cancer Epidemiol Biomarkers Prev
2020
;
29
:
127
32
.
16.
Park
SL
,
Murphy
SE
,
Wilkens
LR
,
Stram
DO
,
Hecht
SS
,
Le Marchand
L
.
Association of CYP2A6 activity with lung cancer incidence in smokers: the multiethnic cohort study
.
PLoS One
2017
;
12
:
e0178435
.
17.
Yuan
JM
,
Butler
LM
,
Gao
YT
,
Murphy
SE
,
Carmella
SG
,
Wang
R
, et al
.
Urinary metabolites of a polycyclic aromatic hydrocarbon and volatile organic compounds in relation to lung cancer development in lifelong never smokers in the Shanghai cohort study
.
Carcinogenesis
2014
;
35
:
339
45
.
18.
Benowitz
NL
,
Dains
KM
,
Dempsey
D
,
Wilson
M
,
Jacob
P
.
Racial differences in the relationship between number of cigarettes smoked and nicotine and carcinogen exposure
.
Nicotine and Tobacco Research
2011
;
13
:
772
83
.
19.
Benowitz
NL
,
Jacob
P
,
Fong
I
,
Gupta
S
.
Nicotine metabolic profile in man: comparison of cigarette smoking and transdermal nicotine
.
J Pharmacol Exp Ther
1994
;
268
:
296
303
.
20.
Murphy
SE
,
Park
SSL
,
Thompson
EF
,
Wilkens
LR
,
Patel
Y
,
Stram
DO
, et al
.
Nicotine N-glucuronidation relative to N-oxidation and C-oxidation and UGT2B10 genotype in five ethnic/racial groups
.
Carcinogenesis
2014
;
35
:
2526
33
.
21.
Hukkanen
J
,
Jacob
P
,
Benowitz
NL
.
Metabolism and disposition kinetics of nicotine
.
Pharmacol Rev
2005
;
57
:
79
115
.
22.
Park
SL
,
Tiirikainen
MI
,
Patel
YM
,
Wilkens
LR
,
Stram
DO
,
Le Marchand
L
, et al
.
Genetic determinants of CYP2A6 activity across racial/ethnic groups with different risks of lung cancer and effect on their smoking intensity
.
Carcinogenesis
2016
;
37
:
269
79
.
23.
Marques
MM
,
Beland
FA
,
Lachenmeier
DW
,
Phillips
DH
,
Chung
FL
,
Dorman
DC
, et al
.
Carcinogenicity of acrolein, crotonaldehyde, and arecoline
.
Lancet Oncol
2021
;
22
:
19
20
.
24.
Kolonel
LN
,
Henderson
BE
,
Hankin
JH
,
Nomura
AMY
,
Wilkens
LR
,
Pike
MC
, et al
.
A multiethnic cohort in Hawaii and Los Angeles: baseline characteristics
.
Am J Epidemiol
2000
;
151
:
346
57
.
25.
Park
SL
,
Carmella
SG
,
Ming
X
,
Vielguth
E
,
Stram
DO
,
Le Marchand
L
, et al
.
Variation in levels of the lung carcinogen NNAL and its glucuronides in the urine of cigarette smokers from five ethnic groups with differing risks for lung cancer
.
Cancer Epidemiol Biomarkers Prev
2015
;
24
:
561
9
.
26.
Haiman
CA
,
Patel
YM
,
Stram
DO
,
Carmella
SG
,
Chen
M
,
Wilkens
LR
, et al
.
Benzene uptake and glutathione S-transferase T1 status as determinants of S-phenylmercapturic acid in cigarette smokers in the multiethnic cohort
.
PLoS One
2016
;
11
:
e0150641
.
27.
Park
SL
,
Carmella
SG
,
Chen
M
,
Patel
Y
,
Stram
DO
,
Haiman
CA
, et al
.
Mercapturic acids derived from the toxicants acrolein and crotonaldehyde in the urine of cigarette smokers from five ethnic groups with differing risks for lung cancer
.
PLoS One
2015
;
10
:
124841
.
28.
Patel
YM
,
Park
SL
,
Carmella
SG
,
Paiano
V
,
Olvera
N
,
Stram
DO
, et al
.
Metabolites of the polycyclic aromatic hydrocarbon phenanthrene in the urine of cigarette smokers from five ethnic groups with differing risks for lung cancer
.
PLoS One
2016
;
11
:
e0156203
.
29.
Cigan
SS
,
Murphy
SE
,
Alexander
BH
,
Stram
DO
,
Hatsukami
DK
,
Le Marchand
L
, et al
.
Ethnic differences of urinary cadmium in cigarette smokers from the multiethnic cohort study
.
Int J Environ Res Public Health
2021
;
18
:
2669
.
30.
Carmella
SG
,
Heskin
AK
,
Tang
MK
,
Jensen
J
,
Luo
X
,
Le
CT
, et al
.
Longitudinal stability in cigarette smokers of urinary eicosanoid biomarkers of oxidative damage and inflammation
.
PLoS One
2019
;
14
:
1
15
.
31.
Percy
C
,
van Holten
V
,
Muir
C
.,
World Health Organization
.
International classification of diseases for oncology, 3rd Edition (ICD-O-3)
.
World Health Organization
;
1990
.
32.
World Health Organization
.
ICD-10 : international statistical classification of diseases and related health problems : 10th revision
.
2nd ed
;
2004
.
33.
Lewis
DR
,
Check
DP
,
Caporaso
NE
,
Travis
WD
,
Devesa
SS
.
US lung cancer trends by histologic type
.
Cancer
2014
;
120
:
2883
92
.
34.
Lunn
M
,
McNeil
D
.
Applying Cox regression to competing risks
.
Biometrics
1995
;
51
:
524
32
.
35.
Stram
DO
,
Park
SL
,
Haiman
CA
,
Murphy
SE
,
Patel
Y
,
Hecht
SS
, et al
.
Racial/ethnic differences in lung cancer incidence in the multiethnic cohort study: an update
.
J Natl Cancer Inst
2019
;
111
:
811
9
.
36.
Benowitz
NL
,
Bernert
JT
,
Foulds
J
,
Hecht
SS
,
Jacob
P
,
Jarvis
MJ
, et al
.
Biochemical verification of tobacco use and abstinence: 2019 update
.
Nicotin Tob Res
2020
;
22
:
1086
97
.
37.
Derby
KS
,
Cuthrell
K
,
Caberto
C
,
Carmella
SG
,
Franke
AA
,
Hecht
SS
, et al
.
Nicotine metabolism in three ethnic/racial groups with different risks of lung cancer
.
Cancer Epidemiol Biomarkers Prev
2008
;
17
:
3526
35
.
38.
Benowitz
NL
,
Pérez-Stable
EJ
,
Herrera
B
,
Jacob
P
.
Slower metabolism and reduced intake of nicotine from cigarette smoking in Chinese-Americans
.
J Natl Cancer Inst
2002
;
94
:
108
15
.
39.
Tanner
JA
,
Tyndale
RF
.
Variation in CYP2A6 activity and personalized medicine
.
J Pers Med
2017
;
7
:
18
.
40.
Stevens
JF
,
Maier
CS
.
Acrolein: Sources, metabolism, and biomolecular interactions relevant to human health and disease
.
Mol Nutr Food Res
2008
;
52
:
7
25
.
41.
Alwis
KU
,
deCastro
BR
,
Morrow
JC
,
Blount
BC
.
Acrolein Exposure in U.S. tobacco smokers and non-tobacco users: NHANES 2005–2006
.
Environ Health Perspect
2015
;
123
:
1302
8
.
42.
Cecil
TL
,
Brewer
TM
,
Young
M
,
Holman
MR
.
Acrolein yields in mainstream smoke from commercial cigarette and little cigar tobacco products
.
Nicotine Tobacco Res
2017
;
19
:
865
70
.
43.
Yuan
JM
,
Gao
YT
,
Wang
R
,
Chen
M
,
Carmella
SG
,
Hecht
SS
.
Urinary levels of volatile organic carcinogen and toxicant biomarkers in relation to lung cancer development in smokers
.
Carcinogenesis
2012
;
33
:
804
9
.
44.
Rostron
BL
,
Wang
J
,
Etemadi
A
,
Thakur
S
,
Chang
JT
,
Bhandari
D
, et al
.
Associations between biomarkers of exposure and lung cancer risk among exclusive cigarette smokers in the Golestan Cohort Study
.
IJERPH
2021
;
18
:
7349
.
45.
International Agency for Research on Cancer (IARC).
IARC Monographs:
cadmium and cadmium compounds
.
Lyon, France: IARC
.;
2012
;
100C
:
121
45
.
46.
International Agency for Research on Cancer (IARC).
Evaluation of carcinogenic risks to humans: beryllium, cadmium, mercury, and exposures in the glass manufacturing industry
.
Lyon, France
:
IARC
;
1993
.
47.
National Toxicology Program.
Report on Carcinogens, 14th Edition
.
Research Triangle Park, NC
:
National Toxicology Program
;
2016
.
48.
Agency for Toxic Substances and Disease Registry (ATSDR).
Toxicological profile: cadmium
;
1999
.
49.
Järup
L
,
Åkesson
A
.
Current status of cadmium as an environmental health problem
.
Toxicol Appl Pharmacol
2009
;
238
:
201
8
.
50.
Nawrot
T
,
Plusquin
M
,
Hogervorst
J
,
Roels
HA
,
Celis
H
,
Thijs
L
, et al
.
Environmental exposure to cadmium and risk of cancer: a prospective population-based study
.
Lancet Oncol
2006
;
7
:
119
26
.
51.
Adams
SV
,
Passarelli
MN
,
Newcomb
PA
.
Cadmium exposure and cancer mortality in the third national health and nutrition examination survey cohort
.
Occup Environ Med
2012
;
69
:
153
6
.
52.
García-Esquinas
E
,
Pollan
M
,
Tellez-Plaza
M
,
Francesconi
KA
,
Goessler
W
,
Guallar
E
, et al
.
Cadmium exposure and cancer mortality in a prospective cohort: The Strong Heart Study
.
Environ Health Perspect
2014
;
122
:
363
70
.
53.
Waalkes
MP
.
Cadmium carcinogenesis
.
Mutat Res
2003
;
533
:
107
20
.
54.
Yuan
JM
,
Carmella
SG
,
Wang
R
,
Tan
YT
,
Adams-Haduch
J
,
Gao
YT
, et al
.
Relationship of the oxidative damage biomarker 8-epi-prostaglandin F2α to risk of lung cancer development in the Shanghai Cohort Study
.
Carcinogenesis
2018
;
39
:
948
54
.
55.
Gào
X
,
Brenner
H
,
Holleczek
B
,
Cuk
K
,
Zhang
Y
,
Anusruti
A
, et al
.
Urinary 8-isoprostane levels and occurrence of lung, colorectal, prostate, breast and overall cancer: results from a large, population-based cohort study with 14 years of follow-up
.
Free Radical Biol Med
2018
;
123
:
20
6
.
56.
Pesch
B
,
Kendzia
B
,
Gustavsson
P
,
Jöckel
KH
,
Johnen
G
,
Pohlabeln
H
, et al
.
Cigarette smoking and lung cancer-relative risk estimates for the major histological types from a pooled analysis of case-control studies
.
Int J Cancer
2012
;
131
:
1210
9
.
57.
Song
MA
,
Benowitz
NL
,
Berman
M
,
Brasky
TM
,
Cummings
KM
,
Hatsukami
DK
, et al
.
Cigarette filter ventilation and its relationship to increasing rates of lung adenocarcinoma
.
J Natl Cancer Inst
2017
;
109
:
djx075
.
58.
National Center for Health Promotion and Education (US) Office on Smoking and Health.
The Health Consequences of Smoking: Nicotine Addiction
.
Rockville, MD
:
Center for Disease Control and Prevention (US)
;
1988
.

Supplementary data