Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

The association between cigarette smoking and inflammation: The Genetic Epidemiology Network of Arteriopathy (GENOA) study

  • Martin Tibuakuu,

    Roles Conceptualization, Formal analysis, Methodology, Supervision, Writing – original draft, Writing – review & editing

    Affiliations Ciccarone Center for the Prevention of Heart Disease, Johns Hopkins School of Medicine, Baltimore, United States of America, St. Luke’s Hospital, Department of Medicine, Chesterfield, United States of America

  • Daisuke Kamimura,

    Roles Formal analysis, Writing – review & editing

    Affiliation University of Mississippi Medical Center, Division of Cardiology, Jackson, United States of America

  • Sina Kianoush,

    Roles Writing – original draft, Writing – review & editing

    Affiliation Ciccarone Center for the Prevention of Heart Disease, Johns Hopkins School of Medicine, Baltimore, United States of America

  • Andrew P. DeFilippis,

    Roles Writing – review & editing

    Affiliation Division of Cardiology, University of Louisville School of Medicine, Louisville, United States of America

  • Mahmoud Al Rifai,

    Roles Writing – original draft, Writing – review & editing

    Affiliations Ciccarone Center for the Prevention of Heart Disease, Johns Hopkins School of Medicine, Baltimore, United States of America, Department of Medicine, University of Kansas School of Medicine, Wichita, United States of America

  • Lindsay M. Reynolds,

    Roles Writing – review & editing

    Affiliation Department of Epidemiology and Prevention, Wake Forest School of Medicine, Winston-Salem, United States of America

  • Wendy B. White,

    Roles Writing – review & editing

    Affiliation Jackson Heart Study, Tougaloo College, Jackson, United States of America

  • Kenneth R. Butler,

    Roles Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – review & editing

    Affiliation University of Mississippi Medical Center, Department of Medicine, Jackson, United States of America

  • Thomas H. Mosley,

    Roles Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – review & editing

    Affiliation University of Mississippi Medical Center, Department of Medicine, Jackson, United States of America

  • Stephen T. Turner,

    Roles Data curation, Funding acquisition, Investigation, Methodology, Project administration, Resources, Writing – review & editing

    Affiliation Division of Nephrology and Hypertension, Mayo Clinic, Rochester, United States of America

  • Iftikhar J. Kullo,

    Roles Data curation, Funding acquisition, Investigation, Methodology, Project administration, Writing – review & editing

    Affiliation Department of Cardiovascular Diseases, Mayo Clinic, Rochester, United States of America

  • Michael E. Hall,

    Roles Supervision, Writing – review & editing

    Affiliation University of Mississippi Medical Center, Division of Cardiology, Jackson, United States of America

  • Michael J. Blaha

    Roles Conceptualization, Formal analysis, Funding acquisition, Methodology, Supervision, Writing – original draft, Writing – review & editing

    mblaha1@jhmi.edu

    Affiliation Ciccarone Center for the Prevention of Heart Disease, Johns Hopkins School of Medicine, Baltimore, United States of America

Abstract

To inform the study and regulation of emerging tobacco products, we sought to identify sensitive biomarkers of tobacco-induced subclinical cardiovascular damage by testing the cross-sectional associations of smoking with 17 biomarkers of inflammation in 2,702 GENOA study participants belonging to sibships ascertained on the basis of hypertension. Cigarette smoking was assessed by status, intensity (number of cigarettes per day), burden (pack-years of smoking), and time since quitting. We modeled biomarkers as geometric mean (GM) ratios using generalized estimating equations (GEE). The mean age of participants was 61 ±10 years; 64.5% were women and 54.4% African American. The prevalence of smoking was 12.2%. After adjusting for potential confounders, 6 of 17 biomarkers were significantly higher among current smokers at a Bonferroni adjusted p-value threshold (p<0.003). High sensitivity C-reactive protein was the most elevated biomarker among current smokers when compared to never smokers [GM ratio = 1.39 (95% CI: 1.23, 1.57); p <0.001]. Among former smokers, each pack-year of cigarettes smoked was associated with a 0.4% higher serum level of hsCRP [GM ratio = 1.004 (95% CI: 1.001, 1.006); p = 0.002] and each 5-year lapsed since quitting was associated with a 4% lower serum level of hsCRP [GM ratio = 0.96 (95% CI: 0.93, 0.99); p = 0.006]. However, we found no significant association of smoking intensity or burden with biomarkers of inflammation among current smokers. HsCRP appears to be the most sensitive biomarker of inflammation associated with cigarette smoking of those investigated, and could be a useful biomarker of smoking-related injury for the study and regulation of emerging tobacco products.

Introduction

Cigarette smoking is one of the leading causes of preventable death in the United States, accounting for almost one in three deaths from cardiovascular disease (CVD) and 20 percent of deaths from ischemic heart disease in adults older than 35 years of age [1, 2]. Despite the progress made over the past 50 years to reduce tobacco use [1, 3], there has been an emergence of novel tobacco products, such as electronic cigarettes, that are claimed to be safer alternatives to combustible cigarettes. However, evidence on the potential cardiovascular toxicity of these products remains unclear [35]. The US Food and Drug Administration's (FDA) authority for tobacco regulation has now been extended to include all tobacco products including electronic nicotine delivery systems (ENDS), cigars, and hookahs [6]. To assist in the study and subsequent regulation of these products, prior to the availability of long-term cardiovascular event data, there is a need to identify sensitive and specific biomarkers of tobacco-mediated cardiovascular injury that could be used to evaluate safety and toxicity.

Despite the increased likelihood of cigarette smoking coexisting with chronic inflammatory conditions such as COPD and cancers, prior studies have demonstrated that inflammation itself is on the causal pathway linking cigarette smoking to CVD outside of these chronic inflammatory states [7, 8]. Indeed, a high burden of inflammation has been shown to identify smokers at high risk for CVD [911]. In a previous study in the Multi-Ethnic Study of Atherosclerosis (MESA) cohort, we identified hsCRP as the most sensitive marker of subclinical cardiovascular injury compared to markers of thrombosis, subclinical myocardial injury, endothelial damage, and vascular function [12]. However, these studies were limited to a few inflammatory biomarkers [9, 10, 1215]. We therefore sought to comprehensively study the association of smoking with 17 biomarkers representing different domains in the inflammatory cascade leading to atherosclerotic CVD [16].

Inflammatory biomarkers included biomarkers of systemic inflammation, including high sensitivity C-reactive protein (hsCRP) and serum amyloid A (SAA) [17]; cell adhesion molecules involved in the early phase of plaque development, including intercellular adhesion molecule (ICAM), vascular cell adhesion molecule (VCAM), E-selectin and P-selectin [18, 19]; cytokines, including interleukin-6 (IL-6), interleukin-18 (IL-18) and tumor necrosis factor alpha (effects mediated through tumor necrosis factor receptors 1 and 2, TNFR1 & TNFR2) [2022]; molecules attracting monocytes to developing plaques, including monocyte chemoattractant protein (MCP) [23]; markers of oxidative stress, including myeloperoxidase (MPO) [24] and receptor for advanced glycated endproducts (RaGE) [25]; and finally enzymes involved in vascular remodeling, including matrix metalloproteinase 2 and 9 (MMP-2 and MMP-9) and their inhibitors, tissue inhibitors of metalloproteinase 1 and 2 (TIMP-1 and TIMP-2) [2628].

Materials and methods

This study was IRB approved at the University of Mississippi Medical Center, Jackson and Mayo Clinic, Rochester. All participants provided informed consent.

Study population

The Genetic Epidemiology Network of Arteriopathy (GENOA) study is one of four research networks that form the National Heart, Lung, and Blood Institute (NHLBI) Family Blood Pressure Program (FBPP) [29]. From its inception in 1995, GENOA's long-term objective was to elucidate the genetics of hypertension and its arteriosclerotic and atherosclerotic target-organ damage, involving the heart, brain, kidneys, and peripheral arteries [30]. Two GENOA cohorts were originally established (1995–2000) through sibships in which at least 2 siblings had essential hypertension diagnosed prior to age 60 years. These include non-Hispanic White Americans from Rochester, MN (N = 1583 at the first exam) and African Americans from Jackson, MS (N = 1854 at the first exam). A third Mexican American cohort was also ascertained through diabetic rather than hypertensive sibships. Because the inflammatory markers of interest were not assessed in this sample, the Mexican American cohort is not included in this report. All siblings in the sibship were invited to participate, both normotensive and hypertensive. During the second exam (2000–2005), approximately 80% of participants were re-recruited to undergo physical examination and provided blood samples for measurement of inflammatory biomarkers.

Assessment of tobacco exposure

Smoking data were collected using questionnaires administered at the first exam and at the time of blood draw during the second GENOA exam. The latter was used in our analysis. Cigarette smoking was self-reported and assessed by status (never, former, current), intensity (number of cigarettes per day) among current smokers, burden (pack-years) among former and current smokers, and time since quitting (in years) among former smokers. We defined current smokers as those who reported smoking at least one cigarette per day during the year prior to exam date. Former smokers were defined as those who reported not smoking for at least a year prior to exam date, but reported smoking at least one cigarette per day in the past. Never smokers were defined as those who denied ever smoking. Pack-years of smoking was calculated by multiplying the number of packs of cigarettes smoked per day by duration of smoking in years. Per study protocol, participants were asked to refrain from smoking at least four hours before clinic visit.

Measurement of inflammatory biomarkers

Blood collection.

Blood was collected by venipuncture after an overnight fast and processed using standardized protocols at each collection site. Blood was centrifuged for 10 min at 4°C, aliquoted in 0.5–1 mL volumes of sodium-citrate plasma, EDTA plasma, and serum and stored at -80°C within 2 h of venipuncture. Aliquots of participants from Jackson, MS samples frozen to -80°C were shipped to Rochester, MN overnight on dry ice. Samples were visually inspected for evidence of thawing and then stored at -80°C. For protein measurements, samples were thawed on ice and aliquoted into bar-coded Eppendorf tubes. The new sample aliquots were re-frozen to -80°C until time of testing at which point they were thawed on ice again. Thus, samples from each collection site were exposed to identical numbers of freeze-thaw cycles for a given assay, ensuring that differential sample handling would not contribute to any subsequently noted ethnic differences in protein levels.

Technical assay performance.

Given the logistical difficulties of measuring a large number of markers in over 2,500 participants, a subset of markers was measured using multiplex assays; IL-6, IL-18, P-selectin, RAGE, TNFRI, E-selectin, MCP-1, MMP-2, MMP-9, TIMP-1, TIMP-2, TNFRII, and ICAM were measured using a contracted service with SearchLight Technologies (Boston, MA). The remaining proteins were either measured in the investigators’ laboratory or the Mayo Immunochemical Core Laboratory (Rochester, MN). We evaluated intra- and inter-assay imprecision at a minimum of one level for each analyte to assess technical assay performance (please see Table A in S1 File). For analytes measured at Mayo, we prospectively determined intra-assay variability, reported as coefficient of variation (CV), by measuring the candidate protein markers in blood samples from volunteers in 10 parallel measurements. We also measured inter-assay variability by measuring the same samples across 10 assay runs. We retrospectively determined precision for the assays performed by SearchLight based on data derived from a blinded, internal plasma control sample. Due to plate-to-plate variations in protein levels in the SearchLight data sets, we developed an algorithm to reduce inter-plate variability; normalized data were used for subsequent analyses [31].

Quality control.

Our quality control program included evaluation of intra-assay variability between duplicate sample measurements and inter-assay variability of quality control materials. Protein levels were measured in duplicates. Sample measurements with CVs >20% were either retested or excluded from the dataset. Acceptable imprecision of measurements from the SearchLight platform was set at <30% due to performance limitations; mean values of samples with CVs >30% were replaced with a singlet value closest to the plate median because retesting was not feasible. We monitored inter-assay variability by measuring 1–3 quality control materials as part of each assay run, and we interpreted the results using a multi-rule approach (13s and 22s Westgard rules) [32]. These rules reject all data included in an assay run if any level of quality control material was three standard deviations (SD) above or below the target value or if 2 or more levels were 2 SD beyond the target value in the same direction. Except for proteins measured on the SearchLight platform, acceptable CV between inter-plate measurements was <20% for all assays and analyses were performed in real-time. Two levels of SearchLight controls and one normal serum control were embedded randomly across study plates and evaluated retrospectively using a modified multi-rule approach as described elsewhere [31]. Sample measurements from failed plates were either repeated or excluded from the data set.

Measurement of covariates

Demographics and medical history were collected from standardized questionnaires. Height was measured by stadiometer and weight by electronic balance to calculate body mass index (BMI) (kg/m2). Resting systolic and diastolic blood pressure (BP) levels were measured with a random zero sphygmomanometer in the right arm. The diagnosis of hypertension was based on either BP measurements (systolic BP≥140 or diastolic BP ≥90 mmHg), previous diagnosis of hypertension, or current treatment with anti-hypertensive medications. Diabetes was considered present if a participant had fasting serum glucose levels ≥126 mg/dL or was receiving treatment with insulin, oral hypoglycemic agents. Information about physical activity, the use of alcohol, BP medications, statins and aspirin was obtained from questionnaires completed by the participants. Serum cholesterol, high-density lipoprotein (HDL) cholesterol, glucose and creatinine were measured by standard enzymatic methods. The Mayo Clinic Quadratic equation was used to estimate glomerular filtration rate (eGFR) [33].

Statistical analyses

A total of 2759 participants (1518 from Jackson, MS and 1241 from Rochester, MN) returned for the second GENOA visit. We excluded 57 participants for missing smoking and baseline confounding variables, resulting in 2,702 participants in our analytical sample. Baseline characteristics of the study population were summarized by categories of smoking status. Data were summarized using counts (proportions) for categorical variables, and means (standard deviations) or medians (interquartile ranges) for continuous variables. Chi-squared test, ANOVA, and Kruskal-Wallis testing were used for comparison across smoking categories where appropriate.

First, all biomarkers were ln-transformed to approximate a normal distribution and we imputed for missing biomarkers (range: 5.9% for hsCRP to 33.9% for MMP-9) (Table A in S1 File) using multiple imputations by chained equations (MICE). Our imputation model included all biomarkers and all our baseline confounding variables. By default, MICE imputed for missing biomarkers in order of increasing missingness using data available for other biomarkers and all confounding variables [34]. Second, to study the association between cigarette smoking and inflammatory biomarkers, we used generalized estimating equations (GEE) that account for intra-familial correlations within GENOA [35]. To compare biomarkers, beta coefficients from GEE models were exponentiated to express geometric mean (GM) ratios [36] of each biomarker with smoking status, intensity, burden, and time since quitting. Models were adjusted for demographics (age, sex, race, education); lifestyle variables (alcohol use, physical activity); other traditional CVD risk factors (body mass index, diabetes, systolic blood pressure, total cholesterol, high density lipoprotein cholesterol), estimated glomerular filtration rate, and medication use (antihypertensives, lipid lowering medications, aspirin). We tested for effect modification by race/ethnicity and sex using multiplicative interaction terms, because plasma levels of several of these biomarkers in this population have been previously shown to differ significantly by these factors. [16]

All statistical analyses were performed using Stata 13 (StataCorp LP, College Station, TX) and statistical significance was generally considered at P value<0.05. To reduce the chances of obtaining false positive results associated with simultaneously testing 17 biomarkers, a statistical significant level was set at a Bonferroni adjusted P-value of <0.003 (0.05/17).

Results

Baseline characteristics of participants

Baseline characteristics of the study population stratified by categories of smoking status are presented in Table 1. The mean age of study participants was 61 ±10 years; 64.5% were women and 54.4% African American. Current smokers constituted 12.2% of the sample and they were likely to be younger (58±10 years), African American (61.2%), currently consume alcohol (65.5%), and use antihypertensive medications (57.0%).

thumbnail
Table 1. Baseline characteristics of study participants by categories of smoking status (2000–2005).

https://doi.org/10.1371/journal.pone.0184914.t001

Inflammatory biomarkers and smoking status

Median plasma concentrations of most inflammatory biomarkers differed significantly by categories of smoking status (Table 2). After adjusting for potential confounders, hsCRP, ICAM, P-selectin, IL-6, TNFR1, and MPO were significantly higher among current smokers at a Bonferroni adjusted p-value threshold (p<0.003), while the following biomarkers did not reach the multiple testing threshold, but were nominally significant (p<0.05): E-selectin, MMP9, and TIMP-2 (Fig 1; Table B in S1 File). Of note, TIMP-2 was lower among current smokers compared to never smokers (p = 0.01). Using a common scale for comparison (GM ratios or percent change), hsCRP was the most elevated biomarker among current smokers when compared to never smokers [GM ratio = 1.39 (95% CI: 1.23, 1.57); p <0.001] (Fig 1; Table B in S1 File).

thumbnail
Table 2. Unadjusted baseline distributions of inflammatory biomarkers by categories of smoking status.

https://doi.org/10.1371/journal.pone.0184914.t002

thumbnail
Fig 1. Adjusted geometric mean (GM) ratios with 95% CI of inflammatory biomarkers for current and former smokers vs. referent never smokers.

Model adjusted for age, sex, education, race/ethnicity, alcohol use, body mass index, physical activity, estimated glomerular filtration rate, systolic blood pressure, diabetes status, total cholesterol, high density cholesterol, family history of myocardial infarction, antihypertensive use, lipid lowering medication use, aspirin use.

https://doi.org/10.1371/journal.pone.0184914.g001

Inflammatory biomarkers and smoking intensity

Although results were not significant for the association between the number of cigarettes smoked per day and levels of inflammatory biomarkers among current smokers, the point estimate for hsCRP was the most elevated for each cigarette smoked per day after adjustment for multiple confounders (Fig 2; Table D in S1 File).

thumbnail
Fig 2. Adjusted geometric mean (GM) ratios with 95% CI of inflammatory biomarkers per unit increase in number of cigarettes smoked/day among current smokers only.

Model adjusted for age, sex, education, race/ethnicity, alcohol use, body mass index, physical activity, estimated glomerular filtration rate, systolic blood pressure, diabetes status, total cholesterol, high density cholesterol, family history of myocardial infarction, antihypertensive use, lipid lowering medication use, aspirin use and smoking duration.

https://doi.org/10.1371/journal.pone.0184914.g002

Inflammatory biomarkers and smoking burden

After multivariable adjustment, smoking burden was not found to be significantly associated with any of our inflammatory biomarkers among current smokers (Fig 3). However, each unit increase in pack-years of smoking was significantly associated with a 0.4% [GM ratio = 1.004 (95% CI: 1.001, 1.006); p = 0.002] higher serum level of hsCRP among former smokers. (Fig 3; Table E in S1 File).

thumbnail
Fig 3. Adjusted geometric mean (GM) ratios with 95% CI of inflammatory biomarkers per unit increase in pack-years of cigarettes smoked among former and current smokers.

Model adjusted for age, sex, education, race/ethnicity, alcohol use, body mass index, physical activity, estimated glomerular filtration rate, systolic blood pressure, diabetes status, total cholesterol, high density cholesterol, family history of myocardial infarction, antihypertensive use, lipid lowering medication use, aspirin use.

https://doi.org/10.1371/journal.pone.0184914.g003

Inflammatory biomarkers and time since quitting

Among former smokers, each 5-year increase in time since quitting was significantly associated with a 4% lower serum level of hsCRP [GM ratio = 0.96 (95% CI: 0.93, 0.99); p = 0.006] after adjusting for multiple variables (Fig 4; Table F in S1 File). Other biomarkers were not significantly associated with time since quitting.

thumbnail
Fig 4. Adjusted geometric mean (GM) ratios with 95% CI of levels of inflammatory biomarkers per 5-year quitting interval.

Model adjusted for age, sex, education, race/ethnicity, alcohol use, body mass index, physical activity, estimated glomerular filtration rate, systolic blood pressure, diabetes status, total cholesterol, high density cholesterol, family history of myocardial infarction, antihypertensive use, lipid lowering medication use, aspirin use.

https://doi.org/10.1371/journal.pone.0184914.g004

Effect modification

We observed evidence of interaction by ethnicity (p for interaction = 0.017) but not sex with smoking status. HsCRP for current smokers compared to never smokers was higher among African Americans [GM ratio = 1.51 (95% CI: 1.29, 1.78); p<0.001] than non-Hispanic whites [GM ratio = 1.20 (95% CI: 1.00, 1.43); p = 0.047] (Table C in S1 File).

Discussion

Among a panel of 17 inflammatory biomarkers that are associated with CVD risk, we demonstrated that current smoking is significantly associated with altered serum levels of 6 biomarkers (hsCRP, ICAM, IL-6, TNFR1, P-selectin and MPO) compared to never smokers at a Bonferroni adjusted p-value threshold (p<0.003). These biomarkers are known to mediate different domains of the inflammatory process of atherosclerotic plaque development, from early endothelial dysfunction to unstable fibrous plaques [37]. For current smokers compared to never smokers, hsCRP was the most elevated for all participants and was even higher among African Americans compared to non-Hispanic whites. Levels of all biomarkers in former smokers were similar to never smokers (Fig 1), highlighting the beneficial effects of smoking cessation on cardiovascular risk factors and disease. While not reaching statistical significance, hsCRP demonstrated the greatest percent elevation with smoking intensity and burden among current smokers. Furthermore, hsCRP was positively associated with smoking burden among former smokers, and inversely associated with time since quitting cigarettes.

Our results are consistent with prior findings that show a relationship between cigarette smoking and vascular inflammation [911, 13, 38, 39]. Smoking has been shown to modulate inflammation by activating the NF-KB pathway thereby inducing the transcription of genes involved in the innate immune response [40, 41]. Thus, exposure to cigarette smoke produces a complicated systemic inflammatory response through the release of cytokines such as IL-6, which subsequently mediates the release of molecules such as hsCRP and SAA by the liver [4042]. HsCRP has been identified as a particularly sensitive inflammatory biomarker and a predictor of cardiovascular events [17, 43]. HsCRP has also been shown to be useful in identifying high-risk smokers who may be candidates for intensive smoking cessation programs [10]. The inverse association of time since quitting with hsCRP in this study adds to the growing evidence that vascular inflammation due to smoking may be reversed by long-term smoking cessation and once again highlights the importance of intensive smoking cessation initiatives for high-risk smokers [10, 38]. Conversely, King et al. reported an inverse association of cessation with F2 isoprostane:creatinine [F2:Cr] ratio and white blood cell (WBC) count but not hsCRP after a year of follow-up among participants who made aided attempt to quit smoking [39]. Although participants had biochemically confirmed 7-day point-prevalence abstinence at 1 year in their study, it is possible that the duration of cessation was not long enough to impact serum levels of hsCRP, which is more distal to F2:Cr in the inflammatory cascade.

Surprisingly, we found no significant association of smoking intensity or burden with any of our inflammatory biomarkers among current smokers. These results though conflicting with our previous findings in MESA [12], are consistent with other studies [14, 38, 44, 45]. Shiels et al. also found no association of smoking intensity or burden with serum concentrations of inflammatory markers [38]. Interestingly, similar results were reported by Ohsawa et al. in a multi-center community-based cohort in Japan [44]. The lack of association of smoking intensity and burden with inflammatory biomarkers among current smokers in the present study may be due to the small proportion of current smokers [328/2702 (12.1%)] (Table 1) included in our sample, as evidenced by the wide confidence intervals in Fig 2. Alternatively, the lack of a significant association between smoking intensity or burden with inflammatory biomarkers among current smokers may suggest a low threshold for a ceiling effects of cigarette smoking on these biomarkers [38]. Thus any amount of smoking may lead to alterations in serum levels of these biomarkers, thereby increasing CVD risk and mortality [46].

The results of the present study have important implications for tobacco regulatory science. The goal of our research, which is funded by the American Heart Association Tobacco Regulation and Addiction Center (A-TRAC), a member of the FDA Tobacco Centers of Regulatory Science (TCORS), is to inform the FDA about sensitive biomarkers of subclinical CVD that could be used to study the potential cardiovascular toxicity of novel tobacco products. From our data, we conclude that measuring hsCRP levels among users of electronic cigarettes and comparing these to current combustible cigarette users and non-users may help demonstrate whether novel tobacco products, such as ENDS, carry the same pro-inflammatory potential as traditional cigarette smoking prior to availability of long-term cardiovascular outcomes data.

The main strength of this study is the well-characterized cohort with a comprehensive phenotyping of 17 inflammatory biomarkers of CVD risk. This study expands upon previous studies by including a broad spectrum of inflammatory biomarkers to find the most sensitive biomarker of early cardiovascular injury due to cigarette smoking in a unique cohort that included mostly African Americans predisposed to CVD. Nonetheless, this study has some limitations. First, the cross-sectional design limits our ability to demonstrate causality and investigate changes in concentrations of biomarkers following tobacco exposure. While unlikely, other unaccounted for chronic inflammatory triggers could explain the higher concentrations of inflammatory biomarkers among current smokers. For instance, African-Americans and Non-Hispanic Whites were recruited from different geographic regions and may have had a difference in circulating levels of biomarkers because of difference in environmental exposures not related to smoking. This could explain the ethnic difference in the levels of some of these biomarkers. Second, the analytical precision of the assays varied, and this may have resulted in the lack of significant association between smoking and some biomarkers with poor test performance. Third, although multiple imputations were used to account for missing biomarkers, the high proportion of missingness for some biomarkers (33.9% for MMP-9) could affect the interpretation of these results. Fourth, smoking exposure was self-reported and did not include information on urine cotinine levels to corroborate self-reported smoking behaviors and avoid misclassification of smoking exposure [47]. Finally, we did not assess immediate inhalational exposure to cigarettes in our study and therefore our results can only be applied to epidemiologic studies of novel tobacco products where there is no immediate exposure.

In conclusion, our data support previous findings that biomarkers of inflammation may allow for early detection of cardiovascular injury due to cigarette smoking. HsCRP appears to be the most sensitive biomarker of inflammation associated with cigarette smoking of the 17 biomarkers investigated, although other markers including MPO, ICAM, IL-6, P-selectin, E-selectin and MMP-9 were also elevated in current smokers. In an era when novel smokeless tobacco products are rapidly emerging in the United States and in other parts of the world, hsCRP and potentially other inflammatory biomarkers may be useful for the study and regulation of these products.

Supporting information

S1 File. Supplementary tables (Table A to Table F).

https://doi.org/10.1371/journal.pone.0184914.s001

(PDF)

Acknowledgments

The authors thank the participants, staff, and investigators of the Genetic Epidemiology Network of Arteriopathy (GENOA) Study for their important contributions.

References

  1. 1. The Health Consequences of Smoking-50 Years of Progress: A Report of the Surgeon General. Atlanta GA 2014. pmid:24455788.
  2. 2. GBD 2015 Risk Factors Collaborators. Global, regional, and national comparative risk assessment of 79 behavioural, environmental and occupational, and metabolic risks or clusters of risks, 1990–2015: a systematic analysis for the Global Burden of Disease Study 2015. Lancet 2016; 388:1659–724. pmid:27733284.
  3. 3. Connolly GN, Alpert HR. Trends in the use of cigarettes and other tobacco products, 2000–2007. JAMA 2008; 299(22):2629–30. pmid:18544722
  4. 4. Southwell BG, Kim AE, Tessman GK, MacMonegle AJ, Choiniere CJ, Evans SE, et al. The marketing of dissolvable tobacco: social science and public policy research needs. Am J Health Promot 2012;26(6):331–2. Epub 2012/07/04. pmid:22747313.
  5. 5. Preventing Tobacco Use Among Youth and Young Adults: A Report of the Surgeon General. Atlanta GA 2012. pmid:22876391
  6. 6. Food and Drugs Administration. FDA proposes to extend its tobacco authority to additional tobacco products, including e-cigarettes. Website: http://www.fda.gov/NewsEvents/Newsroom/PressAnnouncements/ucm394667.htm: Accessed February 04, 2017.
  7. 7. Ambrose JA, Barua RS. The pathophysiology of cigarette smoking and cardiovascular disease: an update. J Am Coll Cardiol 2004; 43:1731–7. pmid:15145091.
  8. 8. Centers for Disease Control and Prevention (US); National Center for Chronic Disease Prevention and Health Promotion (US); Office on Smoking and Health (US). How Tobacco Smoke Causes Disease: The Biology and Behavioral Basis for Smoking-Attributable Disease: A Report of the Surgeon General. Atlanta (GA): Centers for Disease Control and Prevention (US); 2010. 6, Cardiovascular Diseases. https://www.ncbi.nlm.nih.gov/books/NBK53012/.
  9. 9. McEvoy JW, Nasir K, DeFilippis AP, Lima JA, Bluemke DA, Hundley WG, et al. Relationship of cigarette smoking with inflammation and subclinical vascular disease: the Multi-Ethnic Study of Atherosclerosis. Arterioscler Thromb Vasc Biol 2015;35(4):1002–10. pmid:25745060.
  10. 10. McEvoy JW, Blaha MJ, DeFilippis AP, Lima JA, Bluemke DA, Hundley WG, et al. Cigarette smoking and cardiovascular events: role of inflammation and subclinical atherosclerosis from the MultiEthnic Study of Atherosclerosis. Arterioscler Thromb Vasc Biol 2015;35(3):700–9. pmid:25573855.
  11. 11. Levitzky YS, Guo CY, Rong J, Larson MG, Walter RE, Keaney JF Jr., et al. Relation of smoking status to a panel of inflammatory markers: the framingham offspring. Atherosclerosis 2008;201(1):217–24. pmid:18289552.
  12. 12. Al Rifai M, DeFillippis AP, McEvoy JW, Hall ME, Acien AN, Jones MR, et al. The relationship between smoking intensity and subclinical cardiovascular injury: The Multi-Ethnic Study of Atherosclerosis (MESA). Atherosclerosis. http://dx.doi.org/10.1016/j.atherosclerosis.2017.01.021.
  13. 13. Liu J, Liang Q, Frost-Pineda K, Muhammad-Kah R, Rimmer L, Roethig H, et al. Relationship between biomarkers of cigarette smoke exposure and biomarkers of inflammation, oxidative stress, and platelet activation in adult cigarette smokers. Cancer Epidemiol Biomarkers Prev 2011; 20(8):1760–9. pmid:21708936.
  14. 14. Koenig W, Sund M, Frohlich M, Fischer HG, Lowel H, Doring A, et al. C-Reactive protein, a sensitive marker of inflammation, predicts future risk of coronary heart disease in initially healthy middle-aged men: results from the MONICA (Monitoring Trends and Determinants in Cardiovascular Disease) Augsburg Cohort Study, 1984 to 1992. Circulation 1999;99(2):237–42. Epub 1999/01/20. pmid:9892589.
  15. 15. Helmersson J, Larsson A, Vessby B, Basu S. Active smoking and a history of smoking are associated with enhanced prostaglandin F(2alpha), interleukin-6 and F2-isoprostane formation in elderly men. Atherosclerosis 2005; 181(1):201–7. pmid:15939073
  16. 16. Kim CX, Bailey KR, Klee GG, Ellington AA, Liu G, Mosley TH Jr., et al. Sex and ethnic differences in 47 candidate proteomic markers of cardiovascular disease: the Mayo Clinic proteomic markers of arteriosclerosis study. PLoS One 2010;5(2):e9065. pmid:20140090.
  17. 17. Ridker PM, Hennekens CH, Buring JE, Rifai N. C-reactive protein and other markers of inflammation in the prediction of cardiovascular disease in women. N Engl J Med 2000;342(12):836–43. Epub 2000/03/25. pmid:10733371.
  18. 18. Blann AD, Nadar SK, Lip GY. The adhesion molecule P-selectin and cardiovascular disease. Eur Heart J 2003; 24(24): 2166–79. pmid:14659768
  19. 19. Hwang SJ, Ballantyne CM, Sharrett AR, Smith LC, Davis CE, Gotto AM Jr., et al. Circulating adhesion molecules VCAM-1, ICAM-1, and E-selectin in carotid atherosclerosis and incident coronary heart disease cases: the Atherosclerosis Risk In Communities (ARIC) study. Circulation 1997;96(12):4219–25. Epub 1998/01/07. pmid:9416885.
  20. 20. Schulz R, Heusch G. Tumor necrosis factor-alpha and its receptors 1 and 2: Yin and Yang in myocardial infarction? Circulation 2009;119(10):1355–7. pmid:19255338
  21. 21. Ridker PM, Rifai N, Stampfer MJ, Hennekens CH. Plasma concentration of interleukin-6 and the risk of future myocardial infarction among apparently healthy men. Circulation 2000;101(15):1767–72. Epub 2000/04/19. pmid:10769275.
  22. 22. Trøseid M, Seljeflot I, Arnesen H. The role of interleukin-18 in the metabolic syndrome. Cardiovasc Diabetol 2010;9:11. pmid:20331890.
  23. 23. McDermott DH, Yang Q, Kathiresan S, Cupples LA, Massaro JM, Keaney JF Jr., et al. CCL2 polymorphisms are associated with serum monocyte chemoattractant protein-1 levels and myocardial infarction in the Framingham Heart Study. Circulation 2005; 112(8): 1113–20. pmid:16116069
  24. 24. Karakas M, Koenig W, Zierer A, Herder C, Rottbauer W, Baumert J, et al. Myeloperoxidase is associated with incident coronary heart disease independently of traditional risk factors: results from the MONICA/KORA Augsburg study. J Intern Med 2012;271(1):43–50. pmid:21535251.
  25. 25. Lindsey JB, Cipollone F, Abdullah SM, McGuire DK. Receptor for advanced glycation end-products (RAGE) and soluble RAGE (sRAGE): cardiovascular implications. Diab Vasc Dis Res 2009;6(1):7–14. pmid:19156622
  26. 26. Dollery CM, Libby P. Atherosclerosis and proteinase activation. Cardiovasc Res 2006; 69(3): 625–35. pmid:16376322
  27. 27. Galis ZS, Khatri JJ. Matrix metalloproteinases in vascular remodeling and atherogenesis: the good, the bad, and the ugly. Circ Res 2002;90(3):251–62. Epub 2002/02/28. pmid:11861412.
  28. 28. Kai H, Ikeda H, Yasukawa H, Kai M, Seki Y, Kuwahara F, et al. Peripheral blood levels of matrix metalloproteases-2 and -9 are elevated in patients with acute coronary syndromes. J Am Coll Cardiol 1998;32(2) 368–72. pmid:9708462
  29. 29. Granger CB, Van Eyk JE, Mockrin SC, Anderson NL. National Heart, Lung, And Blood Institute Clinical Proteomics Working Group report. Circulation 2004;109(14):1697–703. pmid:15078805.
  30. 30. Multi-center genetic study of hypertension: The Family Blood Pressure Program (FBPP). Hypertension 2002;39(1):3–9. Epub 2002/01/19. pmid:11799070.
  31. 31. Ellington AA, Kullo IJ, Bailey KR, Klee GG. Measurement and quality control issues in multiplex protein assays: a case study. Clin Chem 2009;55(6):1092–9. pmid:19372187.
  32. 32. Westgard JO, Barry PL, Hunt MR, Groth T. A multi-rule Shewhart chart for quality control in clinical chemistry. Clin Chem 1981;27(3):493–501. Epub 1981/03/01. pmid:7471403.
  33. 33. Rule AD, Larson TS, Bergstralh EJ, Slezak JM, Jacobsen SJ, Cosio FG. Using serum creatinine to estimate glomerular filtration rate: accuracy in good health and in chronic kidney disease. Ann Intern Med 2004; 141(12) 929–37. pmid:15611490
  34. 34. Royston P, White IR. Multiple Imputation by Chained Equations (MICE): Implementation in Stata. 2011. 2011;45(4):20. Epub 2011-12-12.
  35. 35. Zeger SL, Liang KY. Longitudinal data analysis for discrete and continuous outcomes. Biometrics 1986;42(1):121–30. Epub 1986/03/01. pmid:3719049.
  36. 36. Newson R. Stata tip 1: The eform() option of regress. The Stata Journal 2003; 3(4): 445. Download from http://www.stata-journal.com/article.html?article=st0054.
  37. 37. Ross R. Atherosclerosis—an inflammatory disease. N Engl J Med 1999;340(2):115–26. Epub 1999/01/14. pmid:9887164.
  38. 38. Shiels MS, Katki HA, Freedman ND, Purdue MP, Wentzensen N, Trabert B, et al. Cigarette smoking and variations in systemic immune and inflammation markers. J Natl Cancer Inst 2014;106(11). pmid:25274579.
  39. 39. King CC, Piper ME, Gepner AD, Fiore MC, Baker TB, Stein JH. Longitudinal Impact of Smoking and Smoking Cessation on Inflammatory Markers of Cardiovascular Disease Risk. Arterioscler Thromb Vasc Biol 2017; 37(2):374–379. pmid:27932354.
  40. 40. Goncalves RB, Coletta RD, Silverio KG, Benevides L, Casati MZ, da Silva JS, et al. Impact of smoking on inflammation: overview of molecular mechanisms. Inflamm Res 2011;60(5):409–24. Epub 2011/02/08. pmid:21298317.
  41. 41. Rom O, Avezov K, Aizenbud D, Reznick AZ. Cigarette smoking and inflammation revisited. Respir Physiol Neurobiol 2013; 187(1) 5–10. pmid:23376061
  42. 42. Pepys MB, Hirschfield GM. C-reactive protein: a critical update. J Clin Invest 2003;111(12):1805–12. pmid:12813013.
  43. 43. Ridker PM, Stampfer MJ, Rifai N. Novel risk factors for systemic atherosclerosis: a comparison of C-reactive protein, fibrinogen, homocysteine, lipoprotein(a), and standard cholesterol screening as predictors of peripheral arterial disease. JAMA 2001;285(19) 2481–5. pmid:11368701
  44. 44. Ohsawa M, Okayama A, Nakamura M, Onoda T, Kato K, Itai K, et al. CRP levels are elevated in smokers but unrelated to the number of cigarettes and are decreased by long-term smoking cessation in male smokers. Prev Med 2005; 41(2):651–6. pmid:15917065.
  45. 45. Tracy RP, Psaty BM, Macy E, Bovill EG, Cushman M, Cornell ES, et al. Lifetime smoking exposure affects the association of C-reactive protein with cardiovascular disease risk factors and subclinical disease in healthy elderly subjects. Arterioscler Thromb Vasc Biol 1997;17(10):2167–76. Epub 1997/11/14. pmid:9351386.
  46. 46. Inoue-Choi M, Liao LM, Reyes-Guzman C, Hartge P, Caporaso N, Freedman ND. Association of Long-term, Low-Intensity Smoking With All-Cause and Cause-Specific Mortality in the National Institutes of Health—AARP Diet and Health Study. JAMA Intern Med 2017; 177(1):87–95. pmid:27918784
  47. 47. Zhang Y, Florath I, Saum KU, Brenner H. Self-reported smoking, serum cotinine, and blood DNA methylation. Environ Res 2016;146:395–403. Epub 2016/02/02. pmid:26826776.