Wood smoke exposure affects lung aging, quality of life, and all-cause mortality in New Mexican smokers
Respiratory Research volume 23, Article number: 236 (2022)
The role of wood smoke (WS) exposure in the etiology of chronic obstructive pulmonary disease (COPD), lung cancer (LC), and mortality remains elusive in adults from countries with low ambient levels of combustion-emitted particulate matter. This study aims to delineate the impact of WS exposure on lung health and mortality in adults age 40 and older who ever smoked.
We assessed health impact of self-reported “ever WS exposure for over a year” in the Lovelace Smokers Cohort using both objective measures (i.e., lung function decline, LC incidence, and deaths) and two health related quality-of-life questionnaires (i.e., lung disease-specific St. George's Respiratory Questionnaire [SGRQ] and the generic 36-item short-form health survey).
Compared to subjects without WS exposure, subjects with WS exposure had a more rapid decline of FEV1 (− 4.3 ml/s, P = 0.025) and FEV1/FVC ratio (− 0.093%, P = 0.015), but not of FVC (− 2.4 ml, P = 0.30). Age modified the impacts of WS exposure on lung function decline. WS exposure impaired all health domains with the increase in SGRQ scores exceeding the minimal clinically important difference. WS exposure increased hazard for incidence of LC and death of all-cause, cardiopulmonary diseases, and cancers by > 50% and shortened the lifespan by 3.5 year. We found no evidence for differential misclassification or confounding from socioeconomic status for the health effects of WS exposure.
We identified epidemiological evidence supporting WS exposure as an independent etiological factor for the development of COPD through accelerating lung function decline in an obstructive pattern. Time-to-event analyses of LC incidence and cancer-specific mortality provide human evidence supporting the carcinogenicity of WS exposure.
Wood smoke (WS), as a major contributor to ambient and indoor combustion emitted particulate matter (PM), has emerged as a critical public health issue in the USA and other high-income countries. As a result of global climate change, wildland fire events have steadily increased over the past 15 years in the USA and other areas of the world and created catastrophically high levels of fine PM (PM2.5) for weeks to months in the affected areas, with toxic smoke travelling to areas hundreds miles away . In addition, WS emitted by wood stoves in millions of USA homes dramatically compromises indoor air quality and is also a major contributor (> 50%) to ambient PM pollution in many urban and rural communities during winter [2,3,4,5]. Studies of rural, wood stove heated homes in the USA found 30% of households having daily PM2.5 levels beyond the USA Environmental Protection Agency (EPA) limit (35 μg/m3) with peak levels resembling those seen in low-income countries [6, 7]. Recreational wood burning and prescribed fires also generate WS exposure, although with much less impact on the air quality.
Epidemiologic studies have linked WS exposure with increased prevalence of respiratory infections, symptoms, and diseases, cancer incidence, and deaths in children and adults in the settings of household heating and cooking, street cooking, fish or meat smoking, or charcoal production in low-income countries where PM2.5 levels from wood burning in those settings could be very high (> 500 µg/m3) [8, 9]. However, WS exposure is typically much lower in the US, which raises uncertainty for the direct extrapolation of the health findings from those high exposure scenarios . Associations with respiratory symptoms or lung function decrement among children living in homes heated by a wood stove were identified in several USA communities [11, 12]. Sustained effects on lung function indicative of airway obstruction were also observed for 2 years following a 45-day exposure to wildfire smoke in Seeley Lake, Montana with daily average PM2.5 levels at 220.9 µg/m3 . Our cross-sectional analyses using the Lovelace Smokers Cohort (LSC) identified associations between “self-reported ever WS exposure for over a year” and a lower percent predicted FEV1 and a higher prevalence of airflow obstruction and chronic mucous hypersecretion (CMH) . However, the impact of WS exposure on lung function decline and its sequelae (e.g., lung cancer [LC] incidence, and disease mortality) in middle-aged and older adults has not been adequately addressed [14,15,16,17,18], particularly in high-income countries.
Disease-specific or generic health-related quality-of-life (HRQoL) questionnaires provide a holistic and quantitative approach assessing the multi-dimensional impacts (e.g., physical, psychological, and social aspects) of diseases or environmental exposures on health in clinical and research settings [19,20,21]. The St. George's Respiratory Questionnaire (SGRQ) is a lung disease-specific instrument whereas the 36-item short-form health survey (SF-36) is a generic HRQoL tool . One study conducted in rural Bolivia found that ventilation remediation of biomass-burning cooking stoves greatly improved the SGRQ scores, supporting the responsiveness of SGRQ scores to PM driven indoor air quality .
Retrospective assessment of total WS exposure from all major sources (e.g., indoor and outdoor) in community-based studies is very challenging. Indoor WS exposure history was often assessed using self-report qualitative or semi-quantitative questionnaires developed by individual studies which are usually lack of standardization [24, 25]. Occupational WS exposure for firefighters, trappers/hunters, rangers, and cooks using wood fired appliances was barely considered in most if not all job exposure matrixes [26, 27]. Health effect studies of wildfire exposure mainly use exposure versus non-exposure comparison approach to identify health differences at aggregate level post exposure . A single question “Have you ever been exposed to WS for 12 months or longer” was implemented in the Lovelace Smokers cohort at baseline to presumably provide a qualitative assessment (yes or no) of cumulative WS exposure for at least 1 year from any exposure sources . With this binary exposure assessment, we identified a strong association between WS exposure and a lower percent predicted FEV1 and a higher prevalence of airflow obstruction and CMH, supporting the validity of this question to some degree . In this study, we further assessed the impact of WS exposure in the LSC using both objective measures (i.e., lung function decline, LC incidence, and death) and two self-assessed HRQoLs (e.g., SGRQ and SF-36). The goal was to find epidemiological evidence supporting (1) the role of WS exposure on the development of chronic obstructive pulmonary disease (COPD) and (2) carcinogenicity of WS exposure in humans.
The LSC was established in 2001 to study sputum and blood biomarkers for LC risk assessment and COPD development in current and former smokers enrolled from the greater Albuquerque area of New Mexico. Study design and inclusion/exclusion criteria were described elsewhere [28, 29]. In brief, subjects were included in the study if they were 40–75 years of age, former or current smokers with at least 10 pack-years of smoking history, free of prior LC history, and able to understand English. At study entry, cohort members completed a battery of questionnaires including demographics, smoking and medical history, Harvard food frequency questionnaire, SF-36, and SGRQ, underwent pre- and post- bronchodilator spirometry testing adhering to the 1994 American Thoracic Society guidelines , and provided blood (lymphocytes, plasma) and induction sputum samples. Cohort members returned every 18 months to update smoking status, general health status and respiratory symptoms, undergo spirometry, and provide biological samples. Active enrollment and on-site follow-up ended in the summer of 2017. Spirometry data were collected from 2511 unique subjects for a total of 11,328 person-visits. Since then, living cohort members are contacted annually by phone calls or mail to collect data for LC incidence or death. This study was approved by the Western Institutional Review Board and all participants signed consent forms. Association analyses in this study were conducted in 2372 LSC subjects with at least one post-bronchodilator spirometry and no missing data for WS exposure and covariates. One hundred and thirty nine LSC subjects were excluded from this study mainly due to missing WS exposure (n = 126) and lack of at least one complete spirometry test (n = 9). Compared with subjects (n = 2372) included in this study, excluded subjects had similar age, BMI, smoking history, and prevalence of airway obstruction (defined as FEV1/FVC ratio ≤ 70% ) and CMH (defined as self-reported cough productive of phlegm for at least 3 months per year for at least 2 consecutive years [i.e., the standard definition of chronic bronchitis] ) at baseline.
Definition of wood smoke exposure
WS exposure was self-reported in response to a question “Have you been exposed to wood smoke for 12 months or longer” as part of the general health survey at study entry.
Health-related quality-of-life (HRQoL) was assessed using the generic health SF-36 questionnaire and the lung disease-specific SGRQ with the recall period of past 4 weeks [32, 33]. The SF-36 encompasses eight domains including physical functioning, role physical, role emotional, social functioning, mental health, vitality, general health perceptions, and bodily pain. The SF-36 scores range from 0 to 100, with higher scores indicating better HRQoL . The SGRQ total score and its activity, symptom, and impact domain subscores range from 0 to 100, with higher score indicating a worse HRQoL . A minimal clinically important difference in SGRQ total score and domain subscores is 4 . SGRQ was collected for all cohort members at baseline and then predominantly for those with FEV1/FVC ratio < 75% at follow-up visits, whereas SF-36 was only collected at baseline visit.
Mortality and lung cancer incidence data
Two National Death Index (NDI) searches were completed in 2014 and mid 2020 and in total identified 380 deaths out of the 2372 subjects (Table 1). Primary causes of death were coded using International Classification of Diseases-10 (ICD-10). Incident lung cancer (n = 72) was identified through NDI searches, obituary data, or by self-report from study subjects or their next of kin. Pathology reports were collected to confirm diagnosis, cancer histology and stage.
Gene promoter methylation in sputum
Promoter methylation of a 12-gene panel originally optimized for LC risk stratification in smokers by our group was measured in eligible sputum samples using nested methylation-specific PCR [36,37,38]. These 12 genes (Additional file 1: Table S1) were selected based on their strong cancer relevance, diversity of function, and specificity of methylation in lung epithelial cells. Methylation status for each individual gene was scored as 0 (unmethylated) or 1 (methylated). A composite methylation index (MI) was created that summed the number of genes methylated.
Linear mixed effects (LME) models with a subject-specific random intercept and slope were used to assess whether WS exposure was associated with a more rapid decline of lung function (FEV1, FVC, and FEV1/FVC ratio) [38, 39]. An interaction term between WS exposure and time in cohort (TIC) and their main effects were included in the LME models. Fixed effects for baseline age, smoking status (current versus former), packyears of smoking, body mass index (BMI), and height, sex, and ethnicity (Hispanic and Others with NHW as the reference) were included in the LME models. A significant estimate for the interaction term indicates the slope of lung function change over time varies by WS exposure with the slope difference estimated by the coefficient of the interaction term. A significant main effect of WS exposure quantifies the difference in baseline lung function (i.e., intercept) between subjects with and without WS exposure. We further extended the LME models by including additional interaction terms between candidate variables (e.g., baseline current smoking, airway obstruction, and CMH) and TIC. Those analyses can further reveal whether baseline subject characteristics such as smoking status, airway obstruction, and CMH might further moderate or mediate the WS exposure and TIC interaction. Similar to the analyses of lung function, LME models with an additional adjustment for education levels (some college or above versus high school or lower) were used to assess whether WS exposure affects overall levels of SGRQ scores over time. Because SF-36 score was only available at baseline, we used linear models to assess the impact of WS exposure with adjustment for baseline age, smoking status, packyears, BMI, and education levels, sex, and ethnicity. Alternative models added age unadjusted Charlson comorbidity score (≥ 1 versus 0), airway obstruction, and CMH at baseline to assess the independent component of effects for WS exposure on SGRQ and SF-36 scores. Kaplan–Meier curve and Cox proportional hazards model assessed the impact of WS exposure on LC incidence, all-cause mortality, and cause-specific mortality. When assessing disease-specific mortality, deaths due to other or unknown causes were treated as competing risk. Baseline age, smoking status, packyears, BMI, and income, education, sex, and ethnicity were adjusted in the Cox model. Mediation analysis with permutation-based statistics was used to quantify the impact of subjective and objective health measurements on the associations between WS exposure and all-cause mortality. Addition analyses were conducted to address how likely the associations seen was due to differential misclassification or confounding by social-economic status and whether excluding cohort members with lung cancer incidence or who died shortly (< 4.86 year) after enrollment affect the associations seen in the entire cohort. All statistical analyses were conducted using SAS 9.4.
Demographics of cohort members with and without WS exposure
Approximately 29% of LSC subjects self-reported to be “ever exposed to WS for over a year” at baseline interview (Table 1). Those subjects were slightly younger and more likely to be Hispanic and current smokers, and had shorter abstinence period if quit smoking. Subjects with WS exposure had lower education and annual income reported at enrollment. Subjects with WS exposure also had shorter duration in cohort and fewer spirometry tests. Moreover, we did not find any associations between WS exposure and sputum methylation index (P = 0.90, Table 1), probably because moderate and heavy cigarette smoking is the driver for acquisition of sputum methylation in this smoker cohort.
Effects of WS exposure on lung function decline
WS exposure was associated with lower FEV1 (− 75.2 ml/s, P = 0.0016) and FEV1/FVC ratio (-1.792%, P < 0.0001) at baseline, but had no impact on baseline FVC (Table 2). A significant impact of WS exposure on annual decline of FEV1 and FEV1/FVC ratio was also identified, but not for the decline of FVC (Table 2; Fig. 1A, B). WS exposure at baseline was associated with steeper slopes indicative of more rapid decline for FEV1 and FEV1/FVC ratio by 4.3 ml/s (P = 0.025) and 0.093% (P = 0.015), respectively. The defense mechanisms against environmental insults compromises as people age. Three-way interaction among WS exposure (yes versus no), baseline age (by year), and TIC (by year) was assessed using LME models and identified significant interactions for FEV1 (estimate = − 8.0 ml/s, P = 0.038) and FEV1/FVC ratio (estimate = − 0.184%, P = 0.017), but not for FVC (estimate = − 2.4 ml/s, P = 0.61). Stratification analyses by the median age of 55.9 year at study entry identified significant impact of WS exposure on the decline of FEV1 and FEV1/FVC ratio in older subjects (age ≥ 55.9 year, Table 2; Fig. 1C, D) but not in the younger group. However, WS exposure affected baseline FEV1 and FEV/FVC ratio regardless of age groups.
Effects of WS exposure on lung function decline not confounded by smoking status, airway obstruction, and CMH
Our previous analyses identified associations of WS exposure with airway obstruction and CMH at baseline in 1861 LSC subjects enrolled prior to 2008 . We replicated these analyses in a larger sample size in this study (n = 2372, Table 1). Because airway obstruction and CMH at baseline may increase the predisposition for lung function decline  and there was a weak correlation between WS exposure and baseline current smoking (P = 0.015, Table 1), we further analyzed whether the associations between WS exposure and decline of FEV1 and FEV1/FVC ratio were independent of baseline current smoking, airway obstruction, and CMH. Our analyses did find baseline airway obstruction and CMH were associated with a more rapid decline of FEV1 and FEV1/FVC ratio (Table 3). Inclusion of these three factors as main effects and their interactions with TIC individually or in combination only slightly (< 15%) reduced the magnitude of impact of WS exposure on lung function decline (Table 3). We further added baseline FEV1 or FEV1/FVC ratio as continuous variables and their interactions with TIC in the linear mixed effect models. Addition of these variables and their interaction terms with TIC did not reduce the magnitude of associations between woodsmoke exposure and lung function decline (data not shown).
Multi-dimensional impact of WS exposure on health
WS exposure also had a strong impact on SGRQ and SF-36 scores and subscores even after adjusting for current smoking, comorbidity, airway obstruction, and CMH at baseline (Table 4). The difference of SGRQ scores and subscores between subjects with and without WS exposure all exceeded or approached the moderately clinically important difference for this instrument (≥ 8, [41, 42]). No interactions were identified between WS exposure and airway obstruction at baseline for affecting SGRQ scores (all Ps > 0.44), suggesting homogeneity of WS exposure effects across airway obstruction status. Promoter methylation of a 12-gene panel in sputum originally optimized for LC risk stratification in smokers by our group was recently shown to be capable of quantifying airway remodeling and predict lung function decline and all-cause mortality [37, 38]. Thus, to discover determinants for SGRQ scores, we further assessed whether sputum methylation affected SGRQ scores. Subjects with ≥ 3 genes methylated had higher symptom scores compared to subjects with 0–2 genes methylated (Additional file 1: Table S1). Additional adjustment of comorbidity, airway obstruction, and CMH at baseline in the model reduced the difference of symptom score by 40%, however the difference remained statistically significant. Individual gene methylation association analyses identified SULF2 and GATA4 as having the strongest associations with higher symptom scores if methylated. No significant associations were identified for activity, impact, or total scores.
Time to event analyses for lung cancer incidence and death
Cohort follow-up identified 72 LC incidences and 380 deaths through the middle of 2020. Time to event was calculated using date at birth (Table 5) or date at baseline visit (Additional file 1: Table S2) because the timing of ever WS exposure was unknown. Subjects with WS exposure had shorter time to event for LC incidence or death from all-causes or specific diseases, although the difference for LC incidence is of borderline (Fig. 2, data not shown). Subjects reporting WS exposure had a shortened lifespan by 3.5 year at death (Table 1). Cox regression analyses identified increased hazard ratio for lung cancer incidence or all-cause or disease-specific death, although the association for LC incidence is of borderline (Table 5, Supplemental Table 2). After the additional adjustment of baseline airway obstruction and comorbidity, the HR of all-cause mortality due to WS exposure remained significantly elevated (HR = 1.39, 95%CI = 1.10–1.75, P = 0.0054, not shown). These findings suggest that WS exposure increased hazard for LC incidence and all-cause and cardiopulmonary- and cancer-specific death by at least 50%.
Health measurements mediating the association between WS exposure and all-cause mortality
We further analyzed whether the associations between WS exposure and all-cause mortality were mediated by health measurements (Table 6). This analysis included subjective (i.e., SGRQ and SF-36 scores) and objective (i.e., FEV1 and FEV1/FVC ratio) measurements that were significantly associated with all-cause mortality (all P values < 0.0003, not shown). Mediational effect sizes ranged from 4.1% to 21.1% with SGRQ total score and SF-36 physical functioning score mediating 21% of magnitude of association. This is comparable to if not higher than the mediational effect size seen for FEV1. Moreover, with additional adjustment for spirometry, significant associations between WS exposure and all-cause mortality remained.
Restricting analyses in subjects with 2 or more visits for spirometry or SGRQ did not change the results observed in the entire study cohort. Stratification analyses by baseline Charlson comorbidity score (≥ 1 versus 0) or education level (college or above versus high school or lower) were conducted to address potential issues of differential classification or confounding of social economic status, respectively (Additional file 1: Table S3). We used all-cause mortality as the outcome to ensure sufficient power in individual subgroup analyses. Significant associations were identified in subgroups with and without baseline comorbidities, suggesting the impact of WS exposure was highly unlikely due to differential misclassification. Moreover, significant associations were also identified in subgroups with and without college education or above, suggesting the impact of WS exposure was unlikely due to confounding from socioeconomic status. Among 380 cohort members who died, 96 subjects died within 4.86 years after study entry. So, additional sensitivity analyses were conducted by excluding these subsets. The estimates for FEV1 decline, FEV1/FVC ratio decline, and SGRQ total scores associated with WS exposure were − 4.3 (1.9) ml/year (P = 0.026), − 0.094 (0.038) %/year (P = 0.015), and 6.8 (0.7) (P < 0.0001) in the 2276 LSC subjects. Furthermore, excluding 72 LC incidents did not change the associations between WS exposure and lung function decline observed in the entire cohort. The estimates for FEV1 decline and FEV1/FVC ratio decline associated with WS exposure were − 4.6 (2.0) ml/year (P = 0.019) and − 0.097 (0.039) %/year (P = 0.012) in the 2300 LSC subjects.
Our study employed a combination of objective health measures (i.e., lung function decline, LC incidence, and disease mortality) and psychometric HRQoL assessment to delineate the pulmonary and overall health effects of WS exposure in a Southwestern USA cohort of middle-aged and older smokers. Compared to subjects without WS exposure, subjects with WS exposure had a more rapid decline of FEV1 (− 4.3 ml/s, P = 0.025) and FEV1/FVC ratio (− 0.093%, P = 0.015) but not of FVC (− 2.4 ml, P = 0.30). Age modified the impacts of WS exposure on lung function decline with stronger impacts seen in older subjects. Because decline in FVC was not associated with WS exposure, we propose that WS exposure leads to an obstructive pattern in COPD development . It is important to note that the adverse effects of WS exposure on baseline FEV1 and FEV1/FVC ratio were independent of age at study entry. In addition to the explanation that WS exposure accelerates lung function decline, exposure to WS exposure at early life could also blunt lung development and impair the maximal attainable lung function, an established risk factor for development of COPD . In addition, strong psychometric evidence was provided that WS exposure impaired all health domains assessed using the generic health SF-36 questionnaire and the lung disease-specific SGRQ. The effect of WS exposure on SGRQ scores doubled the minimal clinically important difference. For both categories of health measurements, we found significant components of the impacts of WS exposure independent of current smoking, airway obstruction, CMH, and comorbidity (for HRQoL measurements) at baseline. Finally, WS exposure was associated with lower age at death and cardiopulmonary diseases and cancers drove this association. Collectively, these findings provide strong support for an accelerated aging of the lung due to chronic WS exposure in a USA population with low exposure to ambient combustion emitted PM.
Carcinogenicity of WS exposure is supported by the abundance of human carcinogens (e.g., polycyclic aromatic hydrocarbons) detected in gaseous and particulate phases of WS and the mutagenicity and carcinogenicity of WS extract in both in vitro and preclinical models. However, carcinogenicity of ambient WS exposure has not been assessed as an independent agent in humans; instead it is considered as a component of outdoor air pollution and PM2.5, which were classified as group 1 carcinogen in humans by the International Agency for Research on Cancer. In 2006, indoor emissions from household combustion of biomass were classified as probably carcinogenic to humans (Group 2A) due to limited evidence for risk of lung cancer from epidemiologic studies. Several review and original research articles published after 2006 provided additional evidence supporting the associations between indoor biomass burning and risk for lung cancer and upper aero-digestive tract cancers mostly using cross-sectional case–control study design in both high- and low-income countries [25, 45,46,47]. Our study provides further evidence based on time-to-event analyses that WS exposure was associated with a > 50% increase risk for LC incidence and cancer-specific mortality, supporting the etiological role of WS exposure in carcinogenesis including LC.
WS composition varies by the type of biomass burned, burning conditions and stages, and secondary atmospheric reactions . WS contains variable amounts of various toxic gases and particulate matters, with the latter containing carbonaceous cores covered by organic matters, metals, and inorganic salts. The majority of WS particles is in a nano-scale range and can deposit deep inside the lungs. Phagocytosis of WS particles by airway macrophages as a major clearance mechanism in acinar airway triggers persistent cytokine/chemokine secretion and generation of other mediators, such as exosomes for downstream toxicity. The development of COPD and lung cancer in susceptible people likely requires decades of repetitive exposure of airways to WS which may occur in the settings of household wood burning for winter heating and frequent and intense wild fires due to climate change in many regions (e.g., the Mountain West) in the United States. Similar to cigarette smoke or diesel exhaust, WS contains numerous established human carcinogens, e.g., polycyclic aromatic hydrocarbons. Moreover, WS contains much higher levels of carbonaceous cores than cigarette smoke [49, 50]. Our study of occupational exposure to black carbon identified a dose–response relationship between macrophage carbon load and genomic instability in peripheral blood, an established cancer biomarker which may be mediated by a mutagenic cytokine TNF-α in circulation . Recent studies from Tesfaigzi’s group identified in vivo and in vitro evidence supporting WS-induced mucous cell hyperplasia which is more prominent in p53 arginine carriers than proline carriers . Fractionation of WS extract identified oxalate that recapitulated the MUC5AC induction by whole WS extract in human primary airway epithelial cells with p53 arginine genotype . These studies provide initial evidence that WS exposure and its specific constituent could induce mucous cell hyperplasia which underlies a specific sub-phenotype of COPD, i.e., chronic bronchitis. Findings of these studies begin to delineate how WS exposure interacts with host factors in COPD development. More researches are needed in this field to explore the molecular mechanisms underlying WS-induced COPD and lung cancer and the involvement of secondary reactions in determining toxicity and health effects of WS exposure.
Concomitant gene methylation detected in exfoliated lung epithelial cells collected in sputum provides an assessment of the extent of field cancerization and is a validated biomarker for diagnosis of primary LC and its recurrence [53, 54]. Moreover, an optimized 12-gene methylation panel in sputum was strongly associated with higher prevalence of CMH, a more rapid decline of FEV1, and increased all-cause mortality in the LSC and Pittsburgh Lung Screening Study cohort [37, 38]. Those 12 genes confer essential functions in biological pathways of cell cycle, cellular senescence, DNA repair, apoptosis, RAS signaling, and invasion that are mechanistically implicated in chronic lung injury and COPD development [14, 15, 55]. The current study identified a strong association between sputum methylation and self-rated SGRQ symptom score with the majority of effect independent of comorbidity and lung diseases, further supporting sputum methylation as an epigenetic endo-type biomarker for lung health [38, 56]. Although only baseline sputum data is available for analyses in this study, the longitudinal collection of sputum samples in the LSC offers a unique opportunity to assess the trajectory of sputum methylation over time and the role of this trajectory in lung aging.
Differential misclassification occurs when misclassification of exposure is not equal between subjects with or without certain health outcomes. Although subjects with WS exposure have higher Charlson comorbidity score at baseline, it is highly unlikely that our results were biased by differential misclassification. First, the LSC is a community-based volunteer cohort enrolling moderate and heavy smokers who have concerns about their health due to past and/or ongoing smoking history. Previous WS exposure should be something that they are least concerned about. Second, our stratification analyses by baseline comorbidity did not identify any difference in magnitude of associations between WS exposure and all-cause mortality, further supporting the very unlikelihood of differential misclassification. It is also unlikely that our results were confounded by lower socioeconomic status associated with WS exposure. Because baseline annual income has high missing rate (21%) an does not capture the wealth information well for the elderly in this cohort, we used education levels as the proxy for social economic status . No difference was identified for the magnitude of associations of all-cause mortality with WS exposure between subjects with and without college education or above. Moreover, education levels and baseline income were included for adjustment in association analyses for mortality.
Although our question for assessing WS exposure presumably provides a qualitative readout of WS exposure for over a year from any exposure sources and has been initially validated by showing strong associations with lung function and pulmonary comorbidities , this binary exposure assessment does not allow us to detect dose–response relationship. Currently, we are developing an artificial intelligence based counting algorithm for macrophage carbon load assay which assesses the lung dose of combustion emitted particulate matter from all sources . Periodical assessment of macrophage carbon load will provide a valid dosimetry over time and allows precise assessment of dose–response relationship.
This study identified epidemiological evidence supporting WS exposure as an independent etiological factor for the development of COPD through accelerating lung function decline in an obstructive pattern. Our study also provides prospective analyses supporting the escalation of carcinogenicity for WS exposure. This study stresses the importance of further characterization of the dose–response relationship between WS exposure and health effects, a key component for risk assessment of residential WS exposure in the USA and other high-income countries. This field of research is very timely as millions of people are not only exposed to residential WS but also to ambient WS from yearly wildfire episodes that are becoming more frequent and prolonged .
Availability of data and materials
Data is available in de-identified form upon completion of Data Use Agreement with Lovelace Biomedical Research Institute who owns the data.
Chronic obstructive pulmonary disease
Lovelace smokers cohort
St. George’s respiratory questionnaire
Orr A, C ALM, Buford M, Ballou S, Migliaccio CT. Sustained effects on lung function in community members following exposure to hazardous PM2.5 levels from wildfire smoke. Toxics 2020; 8.
Wargo J, Alderman N, Wargo L, Bradley J, Addiss S. The harmful effects of wood smoke and the growth of recreational wood burning. 2018.
Ward T, Lange T. The impact of wood smoke on ambient PM2.5 in northern Rocky Mountain valley communities. Environ Pollut. 2010;158:723–9.
Ward T, Trost B, Conner J, Flanagan J, Jayanty RKM. Source apportionment of PM2.5 in a subarctic Airshed—Fairbanks, Alaska. Aerosol Air Qual Res. 2012;12:536–43.
Agency USEP. Strategies for reducing residential wood smoke. pp. 45: United States Environmental Protection Agency 2013:45.
Meek P, Petersen H, Irshad H, Vazquez Guillamet R, McDonald J, Sood A, Tesfaigzi Y. Age of wood stove does not correlate with household air pollution measures. In: ATS 2020 International Conference. 2020.
Semmens EO, Noonan CW, Allen RW, Weiler EC, Ward TJ. Indoor particulate matter in rural, wood stove heated homes. Environ Res. 2015;138:93–100.
Bede-Ojimadu O, Orisakwe OE. Exposure to wood smoke and associated health effects in Sub-Saharan Africa: a systematic review. Ann Glob Health. 2020;86:32.
Chan KH, Kurmi OP, Bennett DA, Yang L, Chen Y, Tan Y, Pei P, Zhong X, Chen J, Zhang J, et al. Solid fuel use and risks of respiratory diseases. A cohort study of 280,000 Chinese never-smokers. Am J Respir Crit Care Med. 2019;199:352–61.
Morandi MT, Ward TJ, Risk Assessment W. Wood smoke risk assessment: defining the questions. Inhal Toxicol. 2010;22:94–8.
Noonan CW, Semmens EO, Ware D, Smith P, Boyer BB, Erdei E, Hopkins SE, Lewis J, Ward TJ. Wood stove interventions and child respiratory infections in rural communities: KidsAir rationale and methods. Contemp Clin Trials. 2020;89: 105909.
Noonan CW, Ward TJ. Asthma randomized trial of indoor wood smoke (ARTIS): rationale and methods. Contemp Clin Trials. 2012;33:1080–7.
Sood A, Petersen H, Blanchette CM, Meek P, Picchi MA, Belinsky SA, Tesfaigzi Y. Wood smoke exposure and gene promoter methylation are associated with increased risk for COPD in smokers. Am J Respir Crit Care Med. 2010;182:1098–104.
Brandsma CA, de Vries M, Costa R, Woldhuis RR, Konigshoff M, Timens W. Lung ageing and COPD: is there a role for ageing in abnormal tissue repair? Eur Respir Rev. 2017; 26.
Thannickal VJ, Murthy M, Balch WE, Chandel NS, Meiners S, Eickelberg O, Selman M, Pardo A, White ES, Levy BD, et al. Blue journal conference. Aging and susceptibility to lung disease. Am J Respir Crit Care Med. 2015;191:261–9.
Meiners S, Eickelberg O, Konigshoff M. Hallmarks of the ageing lung. Eur Respir J. 2015;45:807–27.
Mannino DM, Reichert MM, Davis KJ. Lung function decline and outcomes in an adult population. Am J Respir Crit Care Med. 2006;173:985–90.
Papi A, Beghe B, Fabbri LM. Rate of decline of FEV1 as a biomarker of survival? Am J Respir Crit Care Med. 2021;203:663–5.
Chang KKP, Wong FKY, Chan KL, Wong F, Ho HC, Wong MS, Ho YS, Yuen JWM, Siu JY, Yang L. The impact of the environment on the quality of life and the mediating effects of sleep and stress. Int J Environ Res Public Health. 2020; 17.
Haraldstad K, Wahl A, Andenaes R, Andersen JR, Andersen MH, Beisland E, Borge CR, Engebretsen E, Eisemann M, Halvorsrud L, et al. A systematic review of quality of life research in medicine and health sciences. Qual Life Res. 2019;28:2641–50.
Hansel NN, Putcha N, Woo H, Peng R, Diette GB, Fawzy A, Wise RA, Romero K, Davis MF, Rule AM, et al. Randomized clinical trial of air cleaners to improve indoor air quality and COPD health: results of the Clean Air Study. Am J Respir Crit Care Med. 2021.
Ware JE Jr, Gandek B, Kosinski M, Aaronson NK, Apolone G, Brazier J, Bullinger M, Kaasa S, Leplege A, Prieto L, et al. The equivalence of SF-36 summary health scores estimated using standard and country-specific algorithms in 10 countries: results from the IQOLA Project. International Quality of Life Assessment. J Clin Epidemiol. 1998;51:1167–70.
Alexander D, Linnes JC, Bolton S, Larson T. Ventilated cookstoves associated with improvements in respiratory health-related quality of life in rural Bolivia. J Public Health (Oxf). 2014;36:460–6.
Sood A, Assad NA, Barnes PJ, Churg A, Gordon SB, Harrod KS, Irshad H, Kurmi OP, Martin WJ, 2nd, Meek P, et al. ERS/ATS workshop report on respiratory health effects of household air pollution. Eur Respir J. 2018; 51.
Hosgood HD 3rd, Boffetta P, Greenland S, Lee YC, McLaughlin J, Seow A, Duell EJ, Andrew AS, Zaridze D, Szeszenia-Dabrowska N, et al. In-home coal and wood use and lung cancer risk: a pooled analysis of the International Lung Cancer Consortium. Environ Health Perspect. 2010;118:1743–7.
Kromhout H, Heederik D, Dalderup LM, Kromhout D. Performance of two general job-exposure matrices in a study of lung cancer morbidity in the Zutphen cohort. Am J Epidemiol. 1992;136:698–711.
Sunyer J, Kogevinas M, Kromhout H, Anto JM, Roca J, Tobias A, Vermeulen R, Payo F, Maldonado JA, Martinez-Moratalla J, Muniozguren N. Pulmonary ventilatory defects and occupational exposures in a population-based study in Spain. Spanish Group of the European Community Respiratory Health Survey. Am J Respir Crit Care Med. 1998;157:512–7.
Stidley CA, Picchi MA, Leng S, Willink R, Crowell RE, Flores KG, Kang H, Byers T, Gilliland FD, Belinsky SA. Multivitamins, folate, and green vegetables protect against gene promoter methylation in the aerodigestive tract of smokers. Cancer Res. 2010;70:568–74.
Leng S, Stidley CA, Willink R, Bernauer A, Do K, Picchi MA, Sheng X, Frasco MA, Van Den Berg D, Gilliland FD, et al. Double-strand break damage and associated DNA repair genes predispose smokers to gene methylation. Cancer Res. 2008;68:3049–56.
Standardization of Spirometry, 1994 Update. American Thoracic Society. Am J Respir Crit Care Med. 1995; 152:1107–1136.
Pauwels RA, Buist AS, Calverley PM, Jenkins CR, Hurd SS, GS Committee. Global strategy for the diagnosis, management, and prevention of chronic obstructive pulmonary disease. NHLBI/WHO Global Initiative for Chronic Obstructive Lung Disease (GOLD) Workshop summary. Am J Respir Crit Care Med. 2001;163:1256–76.
Kelly A, Rush J, Shafonsky E, Hayashi A, Votova K, Hall C, Piccinin AM, Weber J, Rast P, Hofer SM. Detecting short-term change and variation in health-related quality of life: within- and between-person factor structure of the SF-36 health survey. Health Qual Life Outcomes. 2015;13:199.
Jones PW. Health status measurement in chronic obstructive pulmonary disease. Thorax. 2001;56:880–7.
Jones PW, Quirk FH, Baveystock CM. The St George’s respiratory questionnaire. Respir Med. 1991;85(Suppl B):25–31 (discussion 33–27).
Jones PW. St. George’s respiratory questionnaire: MCID. COPD. 2005;2:75–9.
Belinsky SA, Liechty KC, Gentry FD, Wolf HJ, Rogers J, Vu K, Haney J, Kennedy TC, Hirsch FR, Miller Y, et al. Promoter hypermethylation of multiple genes in sputum precedes lung cancer incidence in a high-risk cohort. Cancer Res. 2006;66:3338–44.
Bruse S, Petersen H, Weissfeld J, Picchi M, Willink R, Do K, Siegfried J, Belinsky SA, Tesfaigzi Y. Increased methylation of lung cancer-associated genes in sputum DNA of former smokers with chronic mucous hypersecretion. Respir Res. 2014;15:2.
Leng S, Diergaarde B, Picchi MA, Wilson DO, Gilliland FD, Yuan JM, Siegfried JM, Belinsky SA. Gene promoter hypermethylation detected in sputum predicts FEV1 decline and all-cause mortality in smokers. Am J Respir Crit Care Med. 2018;198:187–96.
Leng S, Picchi MA, Tesfaigzi Y, Wu G, Gauderman WJ, Xu F, Gilliland FD, Belinsky SA. Dietary nutrients associated with preservation of lung function in Hispanic and non-Hispanic white smokers from New Mexico. Int J Chron Obstruct Pulmon Dis. 2017;12:3171–81.
Petersen H, Vazquez Guillamet R, Meek P, Sood A, Tesfaigzi Y. Early endotyping: a chance for intervention in chronic obstructive pulmonary disease. Am J Respir Cell Mol Biol. 2018;59:13–7.
Jones PW. Interpreting thresholds for a clinically significant change in health status in asthma and COPD. Eur Respir J. 2002;19:398–404.
Jones PW. Quality of life, symptoms and pulmonary function in asthma: long-term treatment with nedocromil sodium examined in a controlled multicentre trial. Nedocromil Sodium Quality of Life Study Group. Eur Respir J. 1994;7:55–62.
Gonzalez-Garcia M, Torres-Duque C. Wood smoke COPD: a new description of a COPD phenotype? Eur Respir J. 2014;44:261–2.
Lange P, Celli B, Agusti A, Boje Jensen G, Divo M, Faner R, Guerra S, Marott JL, Martinez FD, Martinez-Camblor P, et al. Lung-function trajectories leading to chronic obstructive pulmonary disease. N Engl J Med. 2015;373:111–22.
Sapkota A, Zaridze D, Szeszenia-Dabrowska N, Mates D, Fabianova E, Rudnai P, Janout V, Holcatova I, Brennan P, Boffetta P, Hashibe M. Indoor air pollution from solid fuels and risk of upper aerodigestive tract cancers in central and eastern Europe. Environ Res. 2013;120:90–5.
Kurmi OP, Arya PH, Lam KB, Sorahan T, Ayres JG. Lung cancer risk and solid fuel smoke exposure: a systematic review and meta-analysis. Eur Respir J. 2012;40:1228–37.
Bruce N, Dherani M, Liu R, Hosgood HD 3rd, Sapkota A, Smith KR, Straif K, Lan Q, Pope D. Does household use of biomass fuel cause lung cancer? A systematic review and evaluation of the evidence for the GBD 2010 study. Thorax. 2015;70:433–41.
Black C, Tesfaigzi Y, Bassein JA, Miller LA. Wildfire smoke exposure and human health: significant gaps in research for a growing public health issue. Environ Toxicol Pharmacol. 2017;55:186–95.
Daisey JM, Mahanama KR, Hodgson AT. Toxic volatile organic compounds in simulated environmental tobacco smoke: emission factors for exposure assessment. J Expo Anal Environ Epidemiol. 1998;8:313–34.
Hildemann LM, Markowski GR, Cass GR. Chemical composition of emissions from urban sources of fine organic aerosol. Environ Sci Technol. 1991;25:744–59.
Cheng W, Liu Y, Tang J, Duan H, Wei X, Zhang X, Yu S, Campen MJ, Han W, Rothman N, et al. Carbon content in airway macrophages and genomic instability in Chinese carbon black packers. Arch Toxicol. 2020.
Tassew D, Fort S, Mebratu Y, McDonald J, Chu HW, Petersen H, Tesfaigzi Y. Effects of wood smoke constituents on mucin gene expression in mice and human airway epithelial cells and on nasal epithelia of subjects with a susceptibility gene variant in Tp53. Environ Health Perspect. 2022;130:17010.
Brock MV, Hooker CM, Ota-Machida E, Han Y, Guo M, Ames S, Glockner S, Piantadosi S, Gabrielson E, Pridham G, et al. DNA methylation markers and early recurrence in stage I lung cancer. N Engl J Med. 2008;358:1118–28.
Leng S, Do K, Yingling CM, Picchi MA, Wolf HJ, Kennedy TC, Feser WJ, Baron AE, Franklin WA, Brock MV, et al. Defining a gene promoter methylation signature in sputum for lung cancer risk assessment. Clin Cancer Res. 2012;18:3387–95.
Belinsky SA. Gene-promoter hypermethylation as a biomarker in lung cancer. Nat Rev Cancer. 2004;4:707–17.
Meek PM, Sood A, Petersen H, Belinsky SA, Tesfaigzi Y. Epigenetic change (GATA-4 gene methylation) is associated with health status in chronic obstructive pulmonary disease. Biol Res Nurs. 2015;17:191–8.
Daly MC, Duncan GJ, McDonough P, Williams DR. Optimal indicators of socioeconomic status for health research. Am J Public Health. 2002;92:1151–7.
Leng S, Nguyen K, Picchi M, Meek P, Tesfaigzi Y, Sood A, Belinsky S. Self-reported wood smoke exposure is associated with rapid decline of FEV1 and FEV1/FVC ratio in lovelace smokers cohort: a longitudinal analysis. In: ATS 2021 International Conference; San Diego. 2021.
Part of the findings in this study has been published as an ePoster in the 2021 and 2022 American Thoracic Society annual meeting.
Supported by National Cancer Institute Grants R01 CA097356 and P30 CA118100, National Institute of General Medical Sciences U54 GM104944 and by the State of New Mexico as a direct appropriation from the Tobacco Settlement Fund.
Ethics approval and consent to participate
This study was approved by the Western Institutional Review Board and all participants signed consent forms.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
: Table S1. Impact of sputum methylation on SGRQ scores (n=1798). Table S2. Ever WS exposure on lung cancer incidence and call-cause and disease-specific mortality. Table S3. Stratification analyses of the impacts of WS exposure on all-cause mortality by baseline comorbidity or education levels.
About this article
Cite this article
Leng, S., Picchi, M.A., Meek, P.M. et al. Wood smoke exposure affects lung aging, quality of life, and all-cause mortality in New Mexican smokers. Respir Res 23, 236 (2022). https://doi.org/10.1186/s12931-022-02162-y
- Wood smoke
- Lung function decline
- Health related quality-of-life
- Lung cancer