|Home | About | Journals | Submit | Contact Us | Français|
Although cigarette smoking has been associated with increased human papillomavirus (HPV) detection, its impact on HPV DNA load is unknown.
Study subjects were women who were positive for HPV16 and/or HPV18 at enrollment into the ASCUS-LSIL Triage Study. Assessments of exposure to smoke and sexual behavior were based on self-report. Viral genome copies per nanogram of cellular DNA were measured by multiplex real-time polymerase chain reaction. Linear or logistic regression models were used to assess the relationship between cigarette smoking and baseline viral load.
Of 1,050 women (752 with HPV16, 258 with HPV18, and 40 with both HPV16 and HPV18), 452 (43.0%) were current smokers and 101 (9.6%) were former smokers at enrollment. Baseline viral load was statistically significantly greater for current, compared to never, smokers (P = 0.03 for HPV16; P = 0.02 for HPV18), but not for former smokers. Among current smokers, neither HPV16 nor HPV18 DNA load appeared to vary appreciably by age of smoking initiation, smoking intensity or smoking duration. Results remained similar, when the analysis of smoking-related HPV16 DNA load was restricted to women without detectable cervical abnormality.
Higher baseline HPV16 and HPV18 DNA load was associated with status as a current but not former smoker. A lack of dose-response relationship between cigarette smoking and viral load may indicate a low threshold for the effect of smoking on HPV DNA load.
Cigarette smoking has been found to increase the risk of cervical cancer and its immediate precursor, cervical intraepithelial neoplasia grade 3 (CIN3), among women infected with oncogenic human papillomavirus (HPV) compared to women who do not smoke (1–9). A population-based case-control study in Sweden has showed a synergistic effect between cigarette smoking and both HPV16 positivity and HPV16 DNA load on the development of cervical carcinoma in situ (10). It is generally thought that, in addition to carcinogenic effects of cigarette smoke on cervical tissue, such as genotoxic DNA adducts (11) and the presence of mutagens in cervical mucus (12), the excess risk may be mediated by its effects on immunologic control of HPV infections (13).
Studies of the impact of cigarette smoking on the early natural history of HPV infections have focused on humoral immune response to HPV (14, 15) and on the prevalence (16–21), incidence (22–26), and persistence of HPV infections (27–33). Most previous studies showed some degree of the association between exposure to cigarette smoke and these measurable events in the natural course of HPV infections, although it might be that women who smoke have sexual behaviors that increase the risk of infection. Population data on the effects of cigarette smoking on viral load, an important attribute of HPV infections, have been so far unavailable. A recent in vitro study by Alam et al. (34) demonstrated that exposure of cervical cells to Benzo[a]pyrene, a major carcinogen in cigarette smoke, stimulated higher levels of virion synthesis in HPV-infected cell lines. Since Benzo[a]pyrene can be detected in the cervical mucus of smokers (35), cigarette smoke could increase the viral load of HPV infections.
We therefore evaluated the relationship between cigarette smoking and HPV16 and HPV18 DNA load among women infected with at least one of these two important genotypes, using the baseline data from women who participated in the Atypical Squamous Cells of Undetermined Significance (ASC-US) and Low-Grade Squamous Intraepithelial Lesion (LSIL) Triage Study (ALTS). ALTS was a 2-year multi-center randomized clinical trial that was designed to evaluate strategies for triaging women with mildly abnormal Pap smears.
Subjects eligible for the present study were ALTS participants who had HPV16 and/or HPV18 DNA detected in their enrollment cervical samples (baseline infections) by polymerase chain reaction (PCR)-based reverse line blot assay (36, 37). All ALTS participants had a referral Pap smear interpretation of ASC-US or LSIL. At enrollment, women underwent the same pelvic examination with collection of cervical specimens for cytology and HPV testing. Information on demographic, sexual behavior, lifestyle, and medical history was collected by personal interview at enrollment and during a 2-year follow-up. A detailed description of the ALTS design and the characteristics of the participants are available elsewhere (38, 39). Data and cervical specimens used in this study were those collected at trial enrollment.
In total, 1,071 women met the eligibility criteria, including 759 with baseline HPV16 infections, 258 with baseline HPV18 infections, and 54 with both HPV16 and 18 types (813 with HPV16 and 312 with HPV18). We excluded the following women: 1) 19 (11 positive for HPV16 and 8 positive for HPV18) whose enrollment samples were unavailable for viral DNA quantification; 2) one woman whose enrollment sample was positive for HPV18 but negative for cellular DNA by quantitative assay; and 3) 15 (10 positive for HPV16 and 5 positive for HPV18) who self-reported as current (n = 13) or former smokers (n = 2) but had provided no information on number of cigarettes per day. As a result, there were 792 women (97.4%) with baseline HPV16 infections and 298 women (95.5%) with baseline HPV18 infections in analysis (40 with both types). The protocol for this study was approved by the institutional review boards at the University of Washington and the National Cancer Institute.
Smoking status at enrollment and sexual behavior were ascertained based on self-reported data from the enrollment questionnaire. Women who had ever smoked cigarettes on a regular basis for 6 months or more by the time of enrollment were categorized as ever smokers, and otherwise, were never smokers. Among ever smokers, status as a current or former smoker was self-defined by ALTS participants. No information was collected from former smokers regarding the time of smoking cessation. The intensity and duration of cigarette smoking were dichotomized according to self-reported age of smoking initiation (<16 years versus ≥16 years), number of years smoked (<6 versus ≥6), and number of cigarettes per day (<20 versus ≥20). The number of pack-years of exposure (<4 versus ≥4) that takes into account both intensity and duration of smoking was estimated by number of years smoked multiplied by number of cigarettes per day and then divided by 20. Additionally, we subcategorized the number of cigarettes per day, stratifying by the number of years smoked. The lifetime number of male sex partners and number of male sex partners in the year before study entry were categorized as (0–4, 5–9, or ≥10) and (0–1 versus ≥2), respectively.
Cervical swab samples preserved in Specimen Transport Medium (STM) (Digene Corporation, Silver Spring, MD) were obtained from the NIH biorepository through the ALTS committee. DNAs were extracted from an aliquot of 100 μl STM samples using a procedure described previously (40). The DNA extracts were suspended in 30 μl TE (10mM Tris, pH 8.0, 1mM EDTA) and stored at −20°C until the assay. These samples were initially tested for characterization of HPV16 and HPV18 variants. The remaining was used for viral load quantification.
The number of HPV16 and HPV18 E7 copies and the amount of cellular DNA in enrollment cervical swab samples were measured by multiplex real-time PCR, as described previously (41, 42). Each sample was assayed in triplicate on the same plate. The viral load in each of the triplicate measurements was individually normalized to the input amount of cellular DNA and expressed as E7 copy number per 1 nanogram of cellular DNA.
DNA extracts underwent a thaw-freeze step at least once prior to the measurement of viral load. Some underwent 3 or more times if they were positive for both HPV16 and HPV18 or repeatedly tested for quality control. To examine whether the estimate of viral load would be affected by the limited number of thaw-freeze steps, we compared the amount of cellular DNA between samples positive for both HPV16 and HPV18 and those positive for HPV16 or HPV18 alone and between initial tests and retests. No meaningful difference was observed (data not shown).
Viral E7 DNA was undetectable by real-time PCR in 61 samples previously positive for HPV16 by the reverse line blot assay (33 from never smokers, 25 from current smokers, and 3 from former smokers) and 21 samples previously positive for HPV18 (11 from never smokers, 8 from current smokers, and 2 from former smokers). It is possible that the amount of viral DNA in these samples might have been too small to be detected. We therefore assigned a value of one viral copy per 1 nanogram of cellular DNA to each of these samples for analyses. Considering that the negative result could also be caused by factors such as deletion or integration, we performed parallel analyses with these samples excluded. The results were similar; for simplicity, they were not presented.
The normalized baseline HPV16 and HPV18 DNA load was individually log10-transformed; the mean value of the triplicate measurements was used for analyses. We used student t test to compare viral load by race (white versus nonwhite), current use of oral contraceptives (yes or no), number of male sex partners in the past year, parity (nulliparous versus parous), number of Pap tests per year in the last 5 years (<1 versus ≥1), referral Pap (ASC-US versus LSIL), HPV variant (European versus non-European), and coinfection with other types (yes or no). Differences in baseline viral load by age at enrollment (18–19, 20–24, 25–29, or ≥30 years) and lifetime number of male sex partners were tested by one way ANOVA. The proportions of having ≥2 recent male sex partners by smoking status were compared by chi-square test. Statistical tests were two-sided at the 5% significance level.
Linear regression models (43) were used to examine relationship between baseline viral load and self-reported smoking status at enrollment (current versus noncurrent and former or current versus never). The models were adjusted for age at enrollment, race, parity, referral Pap, current use of oral contraceptives, and coinfection with other HPV types. The linear regression models were also used to compare baseline viral load by a variety of measures of smoking intensity and smoking duration. Each measure of cigarette smoking was modeled independently of one another because of their strong correlation with one another.
To explore whether or not the relationship of interest would be altered by other factors, we evaluated the association of baseline viral load with smoking status overall and separately by race, parity, referral Pap, current use of oral contraceptives, and coinfection with other HPV types using polynomial logistic regression (44). Odds ratios (OR) and 95% confidence intervals (CI) for the overall association were adjusted for age at enrollment, race, parity, referral Pap, current use of oral contraceptives, and coinfection with other HPV types. In this model, the viral load, a continuous variable, was treated as a covariate. Thus, the magnitude of OR indicated risk of being a status of current or former smoker as per 1 log10-unit change of baseline viral load. We noted that most study subjects had a cytologic diagnosis of ≥ASC-US and/or histologic diagnosis of ≥CIN1 at enrollment. To avoid potential disturbance of these concurrent cervical abnormalities to the association of baseline viral load with cigarette smoking, the analysis of baseline HPV16 DNA load was repeated, restricted to women without any detectable cervical lesion.
Of 1,050 ALTS participants (752 with HPV16, 258 with HPV18, and 40 with both HPV16 and HPV18) included in this study, 452 (43.0%) were current smokers by self-report, and 101 (9.6%) were former smokers at enrollment. Of 452 current smokers, 39.8% began smoking prior to 16 years of age, 55.8% smoked ≥6 years, 38.1% smoked ≥20 cigarettes per day, and 48.5% reported 4 or more pack-years. The mean values of the log10-transformed HPV16 E7 copy number per 1 nanogram of cellular DNA was 2.65 (95% CI, 2.50–2.79), 2.87 (95% CI, 2.62–3.11), and 2.92 (95% CI, 2.78–3.06) for never, former, and current smokers, respectively. The corresponding values for women with baseline HPV18 infections were 3.73 (95% CI, 3.49–3.98), 3.25 (95% CI, 2.38–4.13), and 4.13 (95% CI, 3.82–4.43), respectively. Among never smokers, baseline HPV16 DNA load was related to race, parity, and referral Pap; baseline HPV18 DNA load was marginally related to coinfection with other HPV types (Table 1).
As shown in Table 2, after adjustment for age at enrollment, race, parity, referral Pap, current use of oral contraceptives, and coinfection with other HPV types, baseline HPV16 DNA load was statistically significantly higher among current smokers (P = 0.03) but not former smokers (P = 0.35) compared to women who never smoked. Likewise, HPV18 DNA load was statistically significantly higher among current smokers (P = 0.02) but not former smokers (P = 0.23) compared to women who never smoked. The proportion of having ≥2 male sex partners in the past year was highest among current smokers and lowest among never smokers (P = 0.001). Additional adjustment for the number of recent male sex partners did not appreciably alter the relationship between baseline viral load and smoking status (data not shown).
To determine whether there was a dose-response relationship between viral load and cigarette smoking, the smoking status was sub-classified according to age of smoking initiation, number of years smoked, number of cigarettes per day, number of cigarette pack-years, and combination of the years smoked and cigarettes per day. Among current smokers, neither HPV16 nor HPV18 DNA load appeared to vary appreciably by age of smoking initiation, smoking intensity or smoking duration (Table 3). The results remained the same when similar analyses were performed among former smokers (data not shown).
Overall, status as a current smoker at enrollment was statistically significantly associated with baseline viral load (HPV16: OR adjusted = 1.15, 95% CI, 1.02–1.30; HPV18: OR adjusted = 1.23, 95% CI, 1.04–1.45). To examine whether other determinants played a role in the association of interest, data were stratified according to race, parity, referral Pap, current use of oral contraceptives, and coinfection with other HPV types. As shown in Table 4, although the associations with current smoking were not statistically significant in all subgroups, a trend of positive association was the same across the strata. We did not observe an association between viral load and status as a former smoker compared to never smokers overall or in similar stratified analyses (data not shown).
A total of 667 (84.2%) of 792 women with baseline HPV16 infections and 246 (82.6%) of 298 women with baseline HPV18 infections had a cytologic interpretation of ≥ASC-US and/or a histologic diagnosis of ≥CIN1 at enrollment. To eliminate a potential influence of these concurrent cervical abnormalities, the analysis of smoking-related HPV16 DNA load was restricted to women without any cytologic or histologic abnormality at enrollment. The association between status as a current smoker and baseline HPV16 DNA load remained statistically significant (OR adjusted = 1.57, 95% CI, 1.07–2.33). Consistent with the results derived from the overall dataset, no appreciable difference was seen when smoking status was stratified by age of smoking initiation, smoking intensity or smoking duration (data not shown). We did not perform a similar analysis for HPV18 because of the small number of infections.
In this study of women with prevalent HPV16 and/or HPV18 infections at the time of enrollment into ALTS, we found that baseline viral load was statistically significantly higher among current smokers than among never smokers. As shown by an analysis of women with no detectable cytologic and histologic abnormality at enrollment, the association between cigarette smoking and baseline viral load did not appear to be mediated by underlying cervical lesion. Also, the association was not explained by race, parity, referral Pap, current use of oral contraceptives, and coinfection with other HPV types. Concerns remain that the correlation between baseline viral load and cigarette smoking might exist only among certain subgroups and thus adjustment for these factors might not adequately rule out confounding effects. We therefore stratified data by these factors individually. Failure to detect effect modifiers suggests that the association observed was independent of these factors. Further, because the samples used for viral load measurement were collected prior to biopsy or therapy, the association was not distorted by any of these procedures. In analyses of effects of other determinants on baseline viral load in the entire study population, we found that none of those listed in Table 1 except for the referral Pap and coinfection with other HPV types was significantly associated with both HPV16 and HPV18 DNA load (data not shown). The coinfection-related reduction of baseline viral load was reported previously (45).
Before interpreting our results, certain limitations to the present study should be considered. Only samples that tested positive for HPV16 and/or HPV18 by PCR-based reverse line blot assay were used for viral DNA quantification. Thus, samples that contained viral DNA below the detectable threshold were not included in this analysis. If this had occurred more frequently among current smokers than among never smokers, the association between viral load and current smoker status could have been overestimated. It is noteworthy though that among women with HPV16 infections, the proportion of samples that tested positive by PCR-based reverse line blot assay but negative by real-time PCR (presuming presence of tiny amount of viral DNA) was 7.1% for current smokers and 9.3% for never smokers. Among women with HPV18 infections, the proportions were similar between current and never smokers (6.7% versus 6.8%). In this study, exposure to cigarette smoke was based on participants’ self-reported information. Thus, recall bias could be a concern. However, there was no reason to suspect that women with high, compared to low, levels of baseline viral DNA would be more likely to report inaccurate smoking histories. The interview was conducted prior to HPV testing; HPV infections are asymptomatic. Data from a meta-analysis (46) indicated that as compared with biochemical assessments of smoking, self-reports of smoking were accurate in most studies. Information on time of smoking cessation was not collected in ALTS. Thus, we were unable to evaluate its effects on HPV DNA load.
We are aware that although additional adjustment for the number of recent male sex partners did not appreciably alter the relationship between baseline viral load and smoking status, residual confounding effects may remain. Smoking is closely related to sexual behavior. It is possible that the smoking behavior of sex partners, in addition to the number of sex partners, may affect the assessment of risk association. For example, there is no doubt that the extent of passive exposure to smoke is related to the extent a partner smokes. Information on partners’ behavior was not collected in ALTS.
Published data about the relationship between cigarette smoking and HPV DNA load are not available for comparison. Nevertheless, our results are generally consistent with those of other studies showing that cigarette smoking is associated with HPV prevalence (16–19), incidence (22, 23), and persistence (27, 28). Previous studies have shown that cervical mucus of smokers contains measurable amounts of cigarette constituents and their metabolites such as Benzo[a]pyrene (35), nicotine, and nicotine derived nitrosamines 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (47). One possible mechanism for the association between cigarette smoking and HPV DNA load is that smoking may increase cell proliferation and turnover in the transformation zone of the cervix. For example, smoking has been associated with Ki-67 staining, a marker for proliferation and metaplasia (9). HPV uses the host cell DNA replication machinery for its own replication. It is possible that a smoking-related increase in host cell proliferation leads to an increase in replication of HPV16 and HPV18 DNA and/or production of infectious virus. Another possible mechanism is that smoking may increase viral load by weakening the cellular immune response. This is supported by studies showing that smoking has adverse effects on both systemic and local immunity (48–52). For example, smoking has been correlated with impaired NK cytotoxic activity in peripheral blood and unbalanced systemic production of pro- and anti-inflammatory cytokines (50). In addition, smoking has been associated with a substantial decrease in Langerhans cells (51, 52) and helper/inducer T lymphocytes (52) in cervical epithelium. It is possible that the insufficient availability of these cells may cause a decrease in T cell response to HPV infections; consequently, the infected cells are less likely to be eliminated, leading to an increase in HPV DNA load.
We found no evidence that baseline viral load was related to long duration, early initiation, and greater intensity of smoking. The credibility of the association between smoking and viral load might be called into question by this absence of a dose-response relationship. A non-causal interpretation is that the results might be confounded by some unmeasured factors such as passive smoking, tobacco composition, and individual differences in metabolism. It is also possible that smoking might be a marker of recent sexual behavior which is related to recency of infection and consequently the viral load. However, it is more likely that the absence of a dose-response relationship might be due to a low threshold for the effect of smoke on replication of viral genome copies. Supporting for this comes from recent molecular evidence (34) that links dosage-related effects of a cigarette carcinogen to different stages of the HPV life cycle in a raft culture system, i.e., treatment with low concentrations of Benzo[a]pyrene (0.001μM) resulted in the highest magnitude of viral genome copies, while viral titers (a measure of infectious viral particles) remained the same as that of the control; conversely, treatment with high concentrations of Benzo[a]pyrene (1μM) significantly promoted virion synthesis but not viral genome amplification. Whereas the increased genome amplification results in more templates from which the viral oncogene transcripts are produced, the increased virion synthesis favors establishment of a persistent infection by increasing the chances of infecting additional sites. Data on concentrations of Benzo[a]pyrene in the cervical mucus of smokers have not been available. If the relationship between smoking and viral load mimics the Benzo[a]pyrene-mediated regulation of the HPV life cycle, the finding of the absence of a dose-response relationship would not be surprising, because increasingly heavy exposures to smoke might favor production of viral progeny rather than an incremental increase in replication of genome copies.
We noted that baseline HPV16 and HPV18 DNA load did not differ meaningfully between the former and never smokers. This agrees with the findings of an association between HPV infection and current but not former smoking (16–18). The absence of an association between baseline viral load and former smoker status could be due to a limited number of former smokers. Alternatively, it may suggest that concurrent co-localization of the virus and a smoking-related carcinogen at the cervix is a prerequisite for the effect of cigarette smoking on HPV DNA load. If the latter hypothesis is correct, smoking cessation may help with reduction of viral load among women with HPV infections.
In summary, we observed that in a population of women referred for a minor cytologic abnormality, higher HPV16 and HPV18 DNA load was associated with status as a current, but not former smoker. Among current smokers, the viral load did not appear to vary appreciably by the intensity and duration of cigarette smoking. This lack of dose-response relationship may suggest a low threshold for the effect of smoke on HPV DNA load. Data from this study, with viral load as an outcome, and those from others (22, 23, 27, 28), with outcomes of incident and persistent HPV infections, indicate that cigarette smoking affects the early natural history of HPV infections.
Grant supports: Public Health Service grant CA 84396 (L. F. Xi).
L.A.K. receives research funds from Merck Research Laboratories. C.M. is supported by Philip Morris USA Inc., and Philip Morris International. Other authors have no commercial or other associations that might pose a conflict of interest.
This study was part of the project ancillary to the ALTS clinical trial but does not represent the ALTS Group. The authors would like to thank the ALTS Group for providing the biological specimens and HPV typing results and Dr Kathrin Jansen for providing the HPV16 and HPV18 DNA standard. Dr. Craig Meyers was supported by the Pennsylvania Department of Health, Health Research Formula Funding program as a part of the Pennsylvania Tobacco Settlement Legislation and by Philip Morris USA Inc., and Philip Morris International.