Abstract
While the association between smoking and human papillomavirus infection, cervical cancer, and anal cancer has been well studied, evidence on the association between cigarette smoking and anal warts is limited. The purpose of this study was to investigate if cigarette smoking status influences the size of anal warts over time in HIV-infected women in a sample of 976 HIV-infected women from the Women's Interagency HIV Study (WIHS). A linear mixed model was used to determine the effect of smoking on anal wart size. Even though women who were currently smokers had larger anal warts at baseline and slower growth rate of anal wart size after each visit than women who were not current smokers, there was no association between size of anal wart and current smoking status over time. Further studies on the role of smoking and interaction between smoking and other risk factors, however, should be explored.
INTRODUCTION
The prevalence of genital warts among sexually active adults in the USA is estimated to be approximately 1%. 1 HIV-infected persons are more likely to have genital warts than HIV-uninfected persons; they also have a greater risk for recurrence of warts. 2,3 Specifically, anal warts pose a major problem for HIV-infected individuals, and thus should receive special attention for several reasons. They have not been studied separately from genital warts, and yet there are indications that anal warts might be more common than cervical warts in women. 4 Furthermore, once infected with one type of human papillomavirus (HPV), patients are more likely to be infected with other HPV types (both low- and high-risk). In fact, some studies have recently shown that 20–50% of genital warts are co-infected with HPV high-risk types. 5,6 There is also evidence for a strong association between the presence of anal warts and the development of anal intraepithelial neoplasia (AIN), a precancerous lesion for anal cancer. 7 Accordingly, Carter et al. 7 reported that men with anal warts were 4.70 times (95% CI 1.81–12.20) more likely to develop AIN than men without anal warts. Recently, from a Danish study of approximately 50,000 patients with genital warts, Blomberg et al. 8 found that genital wart diagnosis is strongly associated with anal cancer (standardized incidence ratio: 12.5 and 7.8 for men and women, respectively). Finally, the serious economical 9,10 and psychological 11,12 burdens associated with having anal warts must also be considered. In 1997, the estimated cost of HPV burden was $3.8 billion (excluding HPV-related cervical cancer) or more than one-third of the $10 billion spent annually on common sexually transmitted infections (STIs) (excluding HIV) and related syndromes. 9 Also different studies have reported decreased self-esteem and increased psychological distress, embarrassment, anger, shame, negative self-perception, anxiety and relationship difficulties among patients with anogenital warts. 11,12 It is therefore important to determine risk factors for the development, progression or regression of anal warts.
The role of smoking in cervical cancer was first reported by Naguib et al. in 1966. 13 The concentration of nicotine was 45 times higher in cervical tissue than in the serum of smoking women. 14 Tobacco smoke is likely to contribute to carcinogenesis through its impact on immune function, thereby altering the natural history of HPV infection and acting as a co-carcinogen in cervical tissue. 15 Also, the relationships between cigarette smoking and HPV infection, 16–18 as well as cervical and anal cancer 19–30 have been explored and are well-documented. Evidence on the association between cigarette smoking and anogenital warts, in general, and anal warts, in particular, is limited and has only recently received attention in the literature. Accordingly, Feldman et al. 31 reported that the incidence of genital warts was almost three times higher in smokers than non-smokers, both in HIV-infected women (13.3 versus 5.0, respectively) and HIV-uninfected women (1.5 versus 0.5, respectively). However, to date, there has been no study published on the association between cigarette smoking and anal warts.
Given the burden of disease and lack of understanding of risk factors for development of anal warts (i.e. smoking), we investigated whether cigarette smoking status is associated with the size of the largest anal wart in HIV-infected women over time using the public data-set obtained from the Women's Interagency HIV Study (WIHS), an ongoing cohort study of HIV-infected and HIV-uninfected women in the USA.
METHODS
Study population
Data used for the current analysis were obtained from the public data-set (release 09) of WIHS. WIHS is an ongoing prospective study of HIV-infected and HIV-uninfected women from six locations in the USA: Bronx/Manhattan, NY; Brooklyn, NY; Washington, DC; Los Angeles/southern California/Hawaii; San Francisco/Bay Area, CA; and Chicago, IL. Details on the WIHS study were described previously. 32,33 Briefly, WIHS recruited study participants through two enrolment phases: 1994–1995 and 2001–2002. The first enrolment phase was between October 1994 and November 1995. Initially, 2059 HIV-infected women and 569 HIV-uninfected women were recruited from both clinic-based and population-based sources. Inclusion criteria in the first enrolment phase were: (1) being at least 13 years of age; (2) giving informed consent; (3) being tested for HIV; (4) ability to complete the interview in either English or Spanish; (5) ability to travel to and from the clinic site to participate in a baseline visit; and (6) giving blood for laboratory testing. During the second enrolment phase between October 2001 and September 2002, 1144 women were recruited. Besides the above criteria, participants who were recruited in the second enrolment also met the following criteria: (1) documented results of an HIV enzyme-linked immunosorbent assay and confirmatory Western blot for HIV infection or documented HIV-negative status (within 30 days before recruitment); (2) no history of clinical AIDS-related conditions (confirmed by medical record abstraction); (3) documented laboratory testing results of HIV RNA levels and CD4 counts surrounding the highly active antiretroviral therapy (HAART) period for those enroled as HAART exposed; and (4) consent to give specimens. 33 During the first enrolment phase, frequency matching (age, ethnicity, education level, injection drug use since 1978 and total number of sexual partners since 1980) was employed to ensure the comparability between HIV-infected and HIV-uninfected groups.
The WIHS study protocols included a baseline visit and follow-up visits every six months, conducted by trained interviewers and examiners. Information obtained during interview included general medical history, obstetric and gynaecological history, HAART use, alcohol and cigarette use, and sexual behaviours. Medical examination, gynaecological examination and medical record abstraction were conducted during baseline and follow-up visits. Medical examination included height/weight/vital signs, lymph nodes and abdomen. Gynaecological examination included external genitalia, internal vagina and cervix, cervical–vaginal lavage, bimanual and rectal exam, and colposcopy, biopsy, and dysplasia treatment if necessary. Medical record abstraction included development of cancer, infectious diseases or opportunistic infections, and any biopsies, surgeries or hospitalization as well as medications received. 32,33
Details on specimens and laboratory techniques were previously reported by Barkan et al. 32 and Bacon et al. 33 Key clinical information collected in follow-up visits was CD4 and CD8 cell count, HIV sero-status among HIV-uninfected women and Pap smear outcomes. 34
Variables of interest and measurement
Outcome variable
The outcome variable for the current analysis was size of the largest anal wart present at the given visit. For this purpose, only those who had at least one anal wart during the course of follow-up were included in the analysis. During the gynaecological examination, a trained examiner identified the presence of anal warts, and then measured the length and width (in millimeters) of the largest wart. The size (area) of the anal wart was calculated by multiplying the width and length of the reported largest anal wart. The wart was identified as an anal wart if it presented in one of the following locations: ‘anus upper left’, ‘anus lower left’, ‘anus upper right’, ‘anus lower right’, ‘perineum left’ and ‘perineum right’. We assumed that the largest wart is an anal wart if there were multiple warts reported.
Independent variable
The independent variable for this analysis was current smoking status. This was obtained as ‘Yes'/‘No’ answer to the following question at baseline and each subsequent visit: ‘Do you currently smoke cigarettes?’
Other variables
CD4+ cell count (<200, 200–500 and >500 cells/mm3) and HIV viral load (<4000, 4000–20,000, 20,001–100,000, and >100,000 copies/mL) were analysed descriptively, as they are important variables in the WIHS study. Potentially confounding variables included in the current analysis were: race/ethnicity (African-American, Caucasian and others), number of sexual partners in the past six months (0 and ≥1 sexual partners), education level (less than high school education, high school education or general education development [GED], some college, and college graduate or graduate school), annual household income (≤$6000, $6001–$12,000, $12,001–$24,000, and ≥24,001), marital status (married or living with partner; widowed, separated or divorced; and never married), enrolment phase (enrolment phase 1 and enrolment phase 2), and HAART use (‘Yes'/‘No’ if on HAART at that particular visit).
Blood was drawn at each visit for determination of HIV status, CD4+ cell count and HIV viral load. Laboratories certified by the AIDS Clinical Trial Groups measured CD4+ cell count level using an established flow cytometry technique. 32 Serum HIV viral load was measured using the nucleic acid sequence-based amplification assay (NASBA) by Organon Teknika (Oklahoma City, OK, USA). HIV viral load tests were conducted at the National Institute of Allergy and Infectious Diseases, AIDS Program, Virology Assurance HIV RNA Proficiency Program. 32 A person was considered on HAART if she met one of the following criteria: (1) two or more nucleoside reverse transcriptase inhibitors (NRTIs) in combination with at least one protease inhibitor (PI) or one non-nucleoside reverse transcriptase inhibitors (NNRTI); (2) one NRTI in combination with at least one PI and at least one NNRTIs; (3) regimen containing ritonavir and saquinavir in combination with one NRTI and no NNRTI; (4) an abacavir or tenofovir containing regimen of three or more NRTIs in the absence of both PI and NNRTIs, except for the three-NRTI regimens consisting of: abacavir + tenofovir + lamivudine or didanosine + tenofovir + lamivudine. Combination of zidovudine (AZT) and stavudine (d4T) with either a PIT or NNRTI were not considered HAART. Monotherapy is considered as taking one NRTI, or only PI, or only NNRTI. This definition of HAART use followed the guidelines of the US Department of Health and Human Services, 35 the International AIDS-Society Panel Antiretroviral Guidelines 36 and is consistent with previous analyses from WIHS. 37,38
Statistical analysis
The distributions of sociodemographic characteristics were examined. We calculated mean and standard deviation of continuous variables and counts and their respective frequencies of categorical variables. Initially, CD4+ cell count and HIV viral load were provided as continuous variables. We used 10 copies/mL for those whose HIV viral load was suppressed to undetectable levels as it was identified and validated by Norteman et al. 39 CD4+ cell count was categorized into three groups (<200, 200–500 and >500 cells/mm3) and HIV viral load was categorized into four groups (<4000, 4000–20,000, 20,001–100,000 and >100,000 copies/mL) as in previous WIHS studies. 37,40–43
A linear mixed model was employed to determine the relationship between size of the largest anal wart and current smoking status at each visit. We chose this model over other statistical methods because of the following advantages. First, it is able to deal with missing values which are common in longitudinal studies. Second, it adjusts for the highly correlated nature of repeated measurements within and between individuals in longitudinal studies. Third, it is able to deal with the problem of unbalanced measurements (i.e. number of visit in our study) of subjects and the time interval between measurements. 44 In WIHS study, the time interval between measurements was approximately equal (i.e. 6 months between visits).
An unadjusted model was first developed to determine the total variation of growth velocity. 44 Next, an adjusted model was built including the following covariates: number of sexual partners in the past six months, education level, marital status, enrolment phase, HAART use and annual household income. Current smoking status (i.e. independent variable) was treated as a time-dependent variable in both the unadjusted and adjusted models. In the adjusted models, number of sexual partners in the past six months, education level, marital status, annual household income and HAART use were also treated as time-dependent variables. In other words, they were entered into the adjusted model both as a main effect and as a product with time (represented by visit number). The time-independent variables included race/ethnicity and enrolment phase and were entered in to the adjusted model as a main effect only. We use those covariates because they have been identified as potential confounders and been used consistently throughout other analyses in WIHS study. 37,43,45–48 CD4+ cell count and HIV viral load were not included in the final model because we previously reported that there was no association between them and the size of anal warts. 49 The PROC MIXED command of SAS 9.2 statistical package (Cary, NC, USA) was used in the modelling process. 50 All tests were two-sided and P = 0.05 was used as the cut-off for significance.
RESULTS
In the current study (between October 1994 and March 2006), the follow-up has 23 possible visits with 3766 HIV-infected and HIV-uninfected women. Exclusion criteria were HIV-uninfected women (n = 958), women who sero-converted during the study (n = 16) or those with unknown HIV sero-status (n = 1), women without anal warts during the entire follow-up period (n = 1777), and women who underwent treatment for anal warts during the study (n = 38). Women who received treatment for anal warts were excluded because the various treatment modalities could greatly influence the size of anal warts in differing ways during follow-up, and there were not enough participants in this group to conduct a meaningful subgroup analysis. Finally, a sample of 976 women was available for this analysis (Figure 1).
Flowchart of inclusion and exclusion of participants in the current analysis
Baseline sociodemographic characteristics of the WIHS HIV-infected participants in the current study
GED = general education development; SD = standard deviation; WIHS = Women's Interagency HIV Study; HAART = highly active antiretroviral therapy
*Among those with anal warts at baseline (n = 417)
Linear mixed model of size of anal warts and current smoking status in the WIHS HIV-infected participants of the current study in unadjusted and adjusted models
HAART = highly active antiretroviral therapy; SE = standard error; WIHS = Women's Interagency HIV Study
†Model adjusted for number of sexual partner in the last six month, race/ethnicity, HAART use, enrolment, marital status, annual household income and education level
‡Statistically significant at P value <0.05
*Type 3 P = 0.41; **Type 3 P = 0.44; ***Type 3 P = 0.56; ****Type 3 P = 0.75
DISCUSSION
In the current analysis, we did not find an association between smoking status and the size of the largest anal wart over time in HIV-infected women from an ongoing prospective cohort study in the USA. To our knowledge, this is the first study using a linear mixed model to investigate whether smoking status is a predictor for the size of anal warts over time among HIV-infected women. We, therefore, cannot compare our findings directly with any other study. There are, however, a few studies reporting the relationship between smoking status and presence of genital warts. In a previous analysis of a subset of WIHS participants, Feldman et al. 31 reported that current smokers were 5.2 times (95% CI 1.02–26.0) more likely to develop genital warts than non-smokers. The major difference between our study and their study 31 is the outcome variable examined. In our study we examined anal warts only, while Feldman et al. 30 also investigated genital warts. We looked into the changes in the size of anal warts over time (i.e. where anal wart were already present); whereas Feldman et al. 31 examined the presence or absence of genital warts.
One interesting issue is that the largest wart in current smokers was larger than that in non-smokers at baseline, even though it was insignificant in the adjusted model. It is plausible that when a current smoker has large wart at baseline, the growth rate (or speed of development) will be slower in subsequent visits when compared with non-smokers.
Even though we did not find an association between the size of the largest anal wart and current smoking status among HIV-infected women in the current analysis, this relationship should be further explored for several reasons. The effects of cigarette smoke by-products on HPV infection, in general, and the risk of cancer, in particular, have been examined previously with the majority of support coming from the cervical cancer literature. Accordingly, McArdle et al. 51 reported that in the occurrence of cervical neoplasia, there was a reduction of Langerhans cells in smokers that leads to mitigation of the effect of host immunity against HPV. Tobacco smoke likely exerts its actions via two classes of compounds: nitrosamines and polycyclic aromatic hydrocarbons (PAHs). 52 In particular, nitrosamine 4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone (NNK), which is the most active carcinogen in animal models, 53 has been reported in high levels in the cervical mucus of women who smoke compared with non-smokers (mean ± SD: 46.9 ± 32.5 ng/g versus 13.0 ± 9.3 ng/g). 54 Melikian et al. 55 identified benzo[a]pyrene metabolites in cervical mucus and DNA adducts in cervical tissues and suggested that PAHs from tobacco smoke and other environmental sources can be transported to the cervix where they are then metabolized in the cervical epithelium. Recently, Alam et al. 56 found that exposure of cervical cells to benzo[a]pyrene also induced high levels of HPV synthesis, thus facilitating the HPV-associated disease process. Since most of the literature on smoking and HPV disease has focused on cervical neoplasia, we have little to draw from for anal HPV disease. However, we feel confident that the interaction between HPV and smoking could be similar for the development of anal warts.
Another reason is that the relationship between smoking and HPV infection, 15–18 cervical cancer 19–28 and anal cancer 29,30 has been well studied in numerous epidemiological studies. Accordingly, current smokers are 1.6–4.6 times more likely to have pre-cancerous and invasive cervical cancer than non-smokers and that the risk increases with the intensity or duration of smoking (odds ratio [OR] 5.9, 95% CI 1.0–35.6) for those who smoke more than 10 cigarettes per day. Furthermore, cigarette smoking influences not only the incidence or prevalence but also the natural history and pathogenesis of HPV infection. Giuliano et al. 57 found that ‘ever’ smokers maintained an HPV infection significantly longer than women who never smoked (mean duration: 10.7 versus 8.5 months). Smokers were also found to have a lower probability of clearing oncogenic infections than women who never smoked (hazard ratio [HR] 0.44, 95% CI 0.20–0.96, ≤8 cigarettes/day). Recently, Matsumoto et al. 58 reported that smokers have significantly lower regression probability of low-grade cervical abnormalities than non-smokers (55.0% versus 68.8%, P = 0.004).
Smoking has been shown to be both an independent risk factor and a co-factor that interacts with other risk factors, such as CD4+ cell count or prior history of external genital warts to enhance the development of genital warts. In a study of 5622 asymptomatic men, Wiley et al. 59 found that the risk of external genital wart development in a smoker who had a history of external genital warts and who had CD4+ cell count <200 cells/mm3 was 6.9 (95% CI 4.7–10.1), compared with those who did not have a history of external genital warts and who had CD4+ cell count less than 200 cells/mm3. They, however, did not report the difference in the incidence of genitals warts between smokers and non-smokers. Our study focuses on clinical outcomes, and the data-set obtained from the WIHS study did not allow us to investigate the molecular mechanisms of smoking on the size of anal warts. For this reason, further studies of this association at the molecular level are warranted.
Our study has two major strengths. First, using linear mixed modelling allows us to clearly address the association between the size of anal wart and current smoking status over time while other statistical methods cannot. The ability of the model to deal with the problems of high correlation of repeated measurements within and between individuals, missing values, non-linear covariates, and unbalanced measurements greatly enhances our ability to utilize rich data-set to examine these important health outcomes among HIV-infected individuals. Furthermore, the use of linear mixed modelling also allows us to appropriately model the size of anal warts as a continuous outcome variable. Had we categorized the size of anal warts and used different approaches (i.e. logistic regression or Cox-proportional hazard regression), we would not have been able to detect subtle changes in the size of the wart during follow-up. One limitation to the current study is the use of the size of the largest anal wart at each visit as the outcome variable. This does not allow us to follow the same wart over time because the largest wart measured at one visit might not be the same in subsequent visits, especially when there are multiple warts. Even though studying the progression or regression of the same wart over time is not feasible with our data, we felt that studying overall disease burden (as measured by the largest wart) was the best proxy measure available. While we were unable to validate the largest anal wart size as a proxy measure for overall disease burden, no other studies to our knowledge have examined this relationship.
In summary, we did not find evidence for the association between the size of anal warts and current smoking status over time in HIV-infected women. However, our results suggest that, at baseline, women who smoke had much larger warts than those who did not smoke. Further exploration of the role of smoking, the interaction between smoking status with other risk factors (e.g. CD4+ cell count or HIV viral load), and the molecular study of the mechanism of smoking on anal warts over time are warranted.
Footnotes
ACKNOWLEDGEMENTS
Data in this manuscript were collected by the Women's Interagency HIV Study (WIHS) Collaborative Study Group with centres (Principal Investigators) at New York City/Bronx Consortium (Kathryn Anastos); Brooklyn, NY (Howard Minkoff); Washington, DC, Metropolitan Consortium (Mary Young); The Connie Wofsy Study Consortium of Northern California (Ruth Greenblatt); Los Angeles County/Southern California Consortium (Alexandra Levine); Chicago Consortium (Mardge Cohen); Data Coordinating Center (Stephen Gange). The WIHS is funded by the National Institute of Allergy and Infectious Diseases (UO1-AI-35004, UO1-AI-31834, UO1-AI-34994, UO1-AI-34989, UO1-AI-34993 and UO1-AI-42590) and by the Eunice Kennedy Shriver National Institute of Child Health and Human Development (UO1-HD-32632). The study is co-funded by the National Cancer Institute, the National Institute on Drug Abuse and the National Institute on Deafness and Other Communication Disorders. Funding is also provided by the National Center for Research Resources (UCSF-CTSI Grant Number UL1 RR024131). The contents of this publication are solely the responsibility of the authors and do not necessarily represent the official views of the National Institutes of Health.
HNL received funding from the Vietnam Education Foundation (VEF) under VEF Fellowship, the NIH/Fogarty Training Fellowship – grant #D43-TW007669 through the Center for International Training and Research (CITAR), School of Public Health, the University of Texas Health Science Center at Houston, and UTHealth Innovation for Cancer Prevention Research Pre- and Post-doctoral Fellowship, The University of Texas School of Public Health-Cancer Prevention and Research Institute of Texas grant #RP101503. The content is solely responsibility of the authors and does not necessarily represent official views of the Cancer Prevention and Research Institute of Texas.
