- Research article
- Open Access
- Open Peer Review
Increasing incidence of anogenital warts with an urban–rural divide among males in Manitoba, Canada, 1990–2011
BMC Public Health volume 16, Article number: 219 (2016)
Anogenital warts (AGW) are caused by the most common sexually transmitted infection, human papillomavirus. The objective of this study was to examine AGW incidence from 1990 to 2011 by sex, age, income quintile, and residential area category (urban/rural). The study period included the initiation of school-based HPV vaccination for girls in the sixth grade, which began in 2008. The data presented in this paper may also be useful for establishing baseline rates of AGW incidence which may be used to evaluate the success of the school-based HPV immunization program.
Cases of anogenital warts were identified using Manitoba’s administrative databases of Physician Claims and Hospital Discharge Abstracts. Annual age-standardized incidence in Manitoba from 1990 to 2011 was calculated. Incident AGW rates were compared by sex, age group, residential area category (urban/rural), and household income quintile using logistic regression. Joinpoint regression analyses were used to evaluate the time trends of AGW.
Prior to 2000, AGW incidence was higher among females than males. However, from 2000 to 2011 the incidence was higher among males and increased steadily over time. AGW incidence tended to peak in younger age groups among females compared to males. Females and males living in urban areas had nearly twice the odds of AGW occurrence compared to those in rural areas.
There is a need for education about AGW in male population. The upcoming initiation of HPV vaccination among boys may reduce the incidence and should be evaluated.
Anogenital warts incidence was higher among females prior to 2000, but since then has been higher among males than females and has increased steadily.
AGW incidence was higher among males at all age groups except 15–19; by age 25 the incidence was 1.5 to 2 times higher than females.
Both females and males living in urban areas had twice the odds of AGW occurrence compared to those in rural areas.
Anogenital warts (AGW) are caused by the most common sexually transmitted infection, human papillomavirus (HPV) . Asymptomatic cervical HPV infection may be detected among 5 % to 40 % of females of reproductive age . Over time, an individual may be infected with different HPV subtypes and may clear some HPV infections naturally [1, 2]. It is estimated that 90 % of AGW are caused by HPV type 6 and HPV type 11 [3, 4].
AGW may be visible as single or multiple flesh-coloured bumps in the anogenital region. Whether treated or not, AGW are highly infectious and recurrent  and are associated with relatively high treatment costs. A study of the costs associated with AGW in British Columbia from 1998 to 2006 found that on average the treatment of one episode of AGW cost $176 for males and $207 for females . A recent review based on 32 articles published globally on the topic of AGW incidence from 2001 to 2012 found that the incidence of AGW ranged from 160 to 289 per 100,000, and peaked before the age of 24 among females and between 25 and 29 years of age among males . Only one other article described AGW incidence in Manitoba, and it included data up to 2004 . In the year that AGW incidence was last reported, 2004, the incidence of AGW among females was 120 cases per 100,000 persons and among males it was 154 cases per 100,000 persons. From 1985 to 1999, the incidence of AGW in Manitoba was higher among females than males, and showed a downward trend for both from 1992 to 1999 . The incidence in almost all age groups continued to increase among males but decreased among females from 2000 to 2004, and beginning in the year 2000 the incidence has been higher among males than females . It is not known whether the increasing trend of AGW in men will continue. It is also lack information on the characteristics of this increasing trend of AGW in men.
An HPV vaccine (Gardasil), which provides protection against HPV types 6, 11, 16, and 18, has been available free-of-charge to all females in Grade 6 in Manitoba since September 2008. Since November 2012, all females between 9 and 26 years of age with an increased risk of HPV infection (determined by a health care provider) are also eligible to receive the HPV vaccine . Starting in September 2016, the Manitoba government will further expand the HPV vaccination program to include males in Grade 6 or born on or after January 1st 2005 . Population-based baseline data of the incidence of AGWs in females and males is useful for the evaluation of the success of these vaccination initiatives.
The objective of this study was to examine AGW incidence from 1990 to 2011 by sex, age, income quintile, and residential area category (urban/rural), determining whether the increase in AGW rate among males detected in 2004 continued. The study period included the initiation of school-based HPV vaccination for girls (which began in 2008) and but not boys (which will begin in 2016) in the sixth grade. The data presented in this paper will be useful for establishing baseline rates of AGW incidence which may be used to evaluate the success of the school-based HPV immunization program.
Annual age-standardized incidence of AGW in Manitoba from 1990 to 2011 was calculated by sex, age group, income quintiles, and geographic residential area (urban or rural). The two major cities in Manitoba are Winnipeg, Manitoba’s biggest urban centre (estimated population 672,000) and Brandon (estimated population 46,000 population) . These two cities are grouped as the “urban” residential category and include 60 % of Manitoba’s population. The rest of Manitoba is the “rural” residential category.
Data was obtained from the centralized Manitoba Health, Healthy Living and Seniors databases, including: the Manitoba Population Registry, Physician Claims, and Hospital Discharge Abstracts. The Manitoba Population Registry was used to identify all Manitobans who were eligible for provincial health insurance benefits. Physician Claims (including shadow billing) and Hospital Discharge Abstracts were used to identify cases of AGW. Specifically, the de-identified, individual, line-level data from physician and hospital claims (which captured AGW cases diagnosed and treated in Manitoba) were utilized. Shadow billings are claims submitted to the provincial government, by physicians who are not fee-for-service, for administrative purposes only (i.e. as a record of services provided).
Manitoba population data for June 30th of each year were provided by the Health Information Management Branch of Manitoba Health, Healthy Living and Seniors.
Standard protocol approvals, registrations, and patient consents
The University of Manitoba Health Research Ethics Board approved the study. As a public health surveillance activity, a notification by the Manitoba Health Information Privacy Committee for analyzing administrative data was obtained. This study involved the use of anonymous provincial health insurance administrative data and did not involve any direct patient contact. Consent from patients was not required.
AGW cases were identified based on the procedure codes specific to condyloma (Additional file 1: Table S1 and Table S4, S5) from the medical claims data. The condyloma procedure codes are analogous to the Current Procedure Terminology codes used in the United States. The diagnosis and procedure codes, International Classification of Diseases [(ICD), Ninth Revision and 10th Revision (ICD 9/10)] recorded in the hospital records were used to identify the AGW cases in hospital settings (Additional file 1: Table S2 and Table S3). Each incident AGW occurrence (case) was defined as one “episode of care” to distinguish it from subsequent occurrences (see below). AGW-specific condyloma tariff (physician billing) and procedure (Hospital Discharge Abstracts) codes listed in two administrative databases from 1990 to 2011 were used to identify AGW cases.
For those AGW that were treated during a hospital stay and did not generate a record in the physician claims database, a combination of ICD (9/10) diagnosis codes and procedure codes were used to identify the AGW cases.
The case definition used in this study is consistent with the case definition used for previous studies of AGW .
Episode of care of anogenital warts and incident anogenital warts cases
An “episode of care” includes all AGW diagnoses and AGW-related billing which occurred within a twelve month period of time. If another AGW claim occurred more than twelve months after a previous AGW episode of care, it was considered a new AGW episode (i.e. a new AGW case).
Episodes of care were identified using the hospital and physician claims databases. If any of the below occurred more than twelve months after a previous AGW episode of care, it is considered a new AGW episode:
an AGW-related tariff code, OR
an AGW ICD-9/10 code, OR
an AGW-related procedure code followed by another claim within two weeks that had an AGW-related ICD code (ICD9 = 078, Other disease due to viruses and Chlamydia).
If any of the above occurred within twelve months of a previous episode of care, then it is considered as part of the same AGW episode of care.
The diagnostic date for a new AGW case was defined as the earliest date listed in hospital or physician claims for that particular episode of care. If a person died or moved out of province within the twelve month period of the AGW episode of care, the coverage cancellation date, recorded in the Manitoba Population Registry, was used as the end of that episode of care.
Age-standardized incidence was calculated using incident cases as the numerator and mid-year age-specific populations as the denominator. All incidence rates are reported as cases per 100,000 persons in the Manitoba population. Ninety-five percent confidence intervals (95 % CI) of rates were calculated, and formal statistical testing for incidence rate ratios was undertaken using Poisson regression, with the natural log of the population as an offset variable. As personal income data is not available in the administrative database, we linked the median household income at the neighbourhood-level (forward sortation area; ie. first 3 digits of postal code) from Canadian Census data in 2001 to the residential area of patients recorded in the databases as a proxy for socioeconomic status [12–14]. Median household income was grouped into quintiles, less than CAD $17,414 as the lowest quintile (Q1), CAD$17,415–20,411 as second quintile (Q2), CAD$20,412–24,422 as the third quintile (Q3), CAD$24,423–27,848 as the fourth quintile (Q4), and CAD$27,845–35,122 as the highest quintile (Q5). Based on the residential postal codes, data were grouped as rural or urban category. The two major urban centers (Winnipeg and Brandon) were categorized as urban and the rest of the province was categorized as rural. Cases with a missing or unknown postal code or originating from out of province were excluded from the region-specific analyses. Household income quintile categories were obtained by linking the forward sortation area code (first three digits of the postal code) to Canadian census data .
Incident AGW rates were compared by sex, age group (0–14, 15–19, 20–24, 25–29, 30–39, 40–49, 50–59, and 60+ years), residential area category (urban versus rural), five-year period, household income quintile using logistic regression. Results of regression models are summarized as adjusted odds ratios with 95 % confidence intervals (CI). These analyses were conducted using SAS v9.4 (SAS Institute, Cary, NC).
Historical cohort analyses were performed to assess the rates and time trends of AGW in females and males. Joinpoint regression analyses [16–18], developed by SEER (Surveillance Epidemiology and End Results, National Cancer Institute, United States), were used to evaluate the changes in time trends of AGW in females and males, and further by urban and rural residential areas.
This study was approved by the Health Research Ethics Board at University of Manitoba.
Overall time trends of anogenital warts incidence in Manitoba
There were 31,510 AGW cases identifiable in the administrative databases in Manitoba from 1990 to 2011; 15,208 among females and 16,302 among males. The age-standardized incidence of AGW was higher among females compared to males until 1999 (Fig. 1). However, this trend has reversed since 2000 with AGW incidence higher among males than females. Among males, the age-standardized AGW incidence increased 31 % from 124 (95 % CI 114–134) cases per 100,000 persons in 2001 to 175 (95 % CI 164–187) cases per 100,000 persons in 2010 (Fig. 1). Consequently, the male to female incidence rate ratio steadily increased from 1.05 (95 % CI 0.94–1.18) in 2001 to 1.50 (95 % CI 1.35–1.65) in 2010 and 1.39 (95 % CI 1.26–1.54) in 2011. Time trend analysis revealed two distinctive time periods and changing rates of AGW in males and females, respectively. The average annual change in AGW incidence was–3.35 % among females from 1990 to 2000 (p < 0.001), and incidence did not change from 2000 to 2011 (annual percent change 0.79, p = 0.09). The average annual change of AGW incidence in males was −1.86 % from 1990 to 1998 (p = 0.016), and 2.80 % from 1998 to 2011 (p < 0.001).
Compared to the earliest time period (1990–1994), the odds of AGW occurrence among females were consistently low (between 0.75 and 0.83) after adjusting for age group, geographic residential area category, and income quintile (Table 1). Starting in 2005, males had a greater odds of AGW compared to the time period 1990–1994.
Incidence of anogenital warts by Age group and Sex in Manitoba
Among females, the incidence of AGW has been consistently higher among those aged 15–19, 20–24, and 25–29 years while among males it has been consistently highest in the 20–24, 25–29, and 30–39 age groups (Fig. 2). The highest incidence of AGW occurred among those aged 20–24 in both males and females (Table 1 and Fig. 2). As illustrated in Fig. 2, AGW incidence was highest among females aged 15–19 and 20–24 and males aged 20–24 during the period of 1990 to 1994 and among females in this age group the incidence has been substantially lower since then. Among males in this age group, the same high incidence was reached again in 2005–2011.
Urban effects on the incidence of anogenital warts
AGW incidence was higher in urban areas for both females and males (Fig. 3). The trend lines of AGW incidence were parallel for urban and rural for both females and males. After adjusting for age group, 5-year time periods, and income quintile, both urban females and males had an almost 2-fold greater odds of AGW compared to their counterparts living in rural areas (Table 1). Multivariable logistic regression analysis for the time periods of 1990–1999 and 2000–2011 separately revealed two trend changes for males and females. From 1990 to 1999, males had 20 % lower odds of AGW than females (OR = 0.80, 95 % CI 0.75–0.86). However, from 2000 to 2011, males had 7 % greater odds of AGW infection (OR = 1.07, 95 % CI 1.01–1.14) after adjusting for the interaction of residential area and age.
A substantial increasing trend of AGW incidence occurred among both urban and rural males after 1998 (Fig. 3). However, AGW rates in urban and rural females showed a steady downward trend since 1990, at an average downward change of 3.35 % (p < 0.001) per year, and stabilized from 2000 to 2011.
This study provided population-based descriptive epidemiologic information about AGW in Manitoba from 1990 to 2011. The incidence was found to be higher among females than males prior to 2000, which is consistent with previously published data . Since 2000, the incidence of AGW has been higher among males than females and has increased steadily, with the incidence among females remaining steady and therefore the gap in incidence between females and males widening over time. Consequently, the male to female incidence rate ratio also steadily increased as the incidence of AGW among males grew to 50 % higher than females by 2010. A similar pattern of increasing AGW incidence among males was observed in Quebec, where it increased from 83 per 100,000 in 1998 to 103 per 100,000 in 2007 . This is consistent with the sex-specific results of 32 studies conducted globally during the period of 2001 to 2012, whereby the median incidence among males (137 per 100,000 per year) was indeed higher than that among females (120.5 per 100,000 per year) . The changing sex ratio in sexually transmitted diseases has been suggested as a surrogate marker of changes in sexual behaviors among men who have sex with men (MSM) . The rising trends of AGW in males is concurrent with the increasing trends of gonorrhea and syphilis in males in Canada . The resurging syphilis in males is largely attributed to the high-risk sexual practices among some MSM, such as using club drugs, multiple sexual partners, and larger sexual networks formed using the internet [22–27]. Therefore, changes in the sexual practices among MSM may also be an underlying factor for the rising AGW incidence among males in Manitoba. However, because sexual behavior and sexual orientation data are not available in the health insurance administrative databases, this study cannot assess the possible relationship between rising AGW incidence and trends in sexual behaviours.
AGW incidence tended to peak in younger age groups among females compared to males. The male to female incidence ratio indicated an often much higher incidence among males at all age groups except 15–19, where the incidence among males was only 40 % that of females. However, by the age of 25 the incidence of AGW among males was one-and-a-half to two times the incidence among females. These results are consistent with a review of 32 articles published globally on the topic of AGW incidence from 2001 to 2012, which found that the incidence of AGW peaked before the age of 24 among females and between 25 and 29 years of age among males [7.] However, the incidence observed in Manitoba over this time period among those aged 20–24 (520 and 440 cases per 100,000 among females and males, respectively) was notably higher than that observed elsewhere in Canada for the same age group, which ranged from 338 to 342 per 100,000 among females and 270 to 300 among males in BC and Quebec, respectively, for a similar time period [6, 28].
Both females and males living in urban areas had almost twice the odds of AGW occurrence compared to those in rural areas. The same temporal trend in which the incidence among males exceeded that of females and increased over time beginning in 2000 was observed in both urban and rural areas. The higher incidence of AGW in urban areas compared to rural areas may relate to greater MSM populations in urban areas  and differences in sexual mixing patterns, which are known to be associated with HIV, gonorrhea, and chlamydia risk [30–32]. It has been reported that for HPV type 6 and HPV type 11, which are most commonly associated with AGW, the prevalence among men who have sex with men (MSM) was very high (39.2–53.2 %) . The higher AGW incidence and increasing trend that we observed among urban males could be partially related to high risk sexual behaviours in urban MSM populations . It is also possible that individuals in urban areas are more likely to seek healthcare for AGW, and therefore were identifiable as AGW cases using the provincial administrative databases. Similar to our results, a study of AGW incidence in Germany from 2005–2006 revealed a higher incidence in urban areas compared to other regions . To our knowledge these are the only reports of an observed urban–rural divide in AGW incidence. Further research is needed to elucidate possible explanations for the observed urban–rural differences in AGW incidence.
An online survey conducted across Canada in 2011 found that on average only half of the respondents believe AGW to be a serious STI, only one-quarter knew that it is associated with HPV, and 57 % were either not sure whether a condom would prevent transmission or believed that it would . Other studies have indicated not only a lack of perceived risk for AGW among men who have sex with men in North America but also a lack of knowledge about the etiologic cause of AGW and the link between HPV infection and anal cancer [35, 36]. However, once informed of the role of HPV, hypothetical acceptance of HPV vaccine was common [35, 36]. HPV has been framed as a female problem, with emphasis on its etiologic role in cervical cancer and the targeting of vaccination programs towards females. This may partially explain the lack of awareness and increasing incidence among males since 2000. Education about HPV, its prevention, and the diseases it may cause are needed among both female and male populations. It has been recommended that HPV vaccination for boys could be an effective and cost-saving strategy for oropharyngeal cancer  and AGW . An upcoming HPV immunization program among boys in Grade 6 in Manitoba will offer an important approach to reducing the high incidence of AGW among males. Research of the cost-effectiveness of the HPV vaccine in Manitoba and an evaluation of the acceptability and effectiveness of the Grade 6 vaccination program in reducing AGW incidence would provide useful information for programs.
The data presented in this report is subject to a number of limitations. It likely underestimates the actual occurrence of AGW in Manitoba as it relies on diagnosed cases and does not include those cases for which treatment was not sought. Further, incomplete shadow billing by physicians who are not fee-for-service and the exclusion of emergency department data may also result in an underestimate. The data included in this study were obtained from existing administrative databases compiled for other purposes and subject to the coding practices used over time. Although the accuracy of ICD codes and completeness of the diagnoses captured in hospital discharge records have been questioned [39, 40], the overall quality of the administrative databases has been demonstrated [41, 42]. Cases were based on clinical observations and not confirmed by laboratory results, though AGW clinical diagnoses have been found to correlate well with histological results .
Strengths of this population-based study include the fact that it is free from selection bias and the limitations of self-report data. Results of this study are generalizable to other similar populations with a similar medical system.
Since 2000 the incidence of AGW among males in Manitoba has been increasing and is particularly high among those residing in urban areas. This observation is consistent with the increasing incidence of syphilis and gonorrhea among males in Manitoba and may be related to increasing sexual risk behaviours among some urban MSM populations. Lower AGW rates in rural males may be related to lower risk behaviours or healthcare seeking behavior if fewer seek healthcare, resulting in a smaller proportion of AGW cases being diagnosed compared to their urban counterparts. Moreover, HPV has been framed largely as a female problem, which may partially explain the lack of awareness and increasing incidence among males. More sexual health education may help patients overcome the stigma of early diagnosis of sexually transmitted infections and improve sexual health in general. There is a need for education about HPV, its prevention, and the diseases it may cause among both female and male populations. Also, effectiveness studies of the HPV vaccination program among Grade 6 males in Manitoba, to begin in 2016, will generate useful evidence about the importance of targeting prevention activities towards males.
Gross G, Pfister H. Role of human papillomavirus in penile cancer, penile intraepithelial squamous cell neoplasias and in genital warts. Med Microbiol Immunol. 2004;193(1):35–44.
Franco EL, Villa LL, Sobrinho JP, et al. Epidemiology of acquisition and clearance of cervical human papillomavirus infection in women from a high-risk area for cervical cancer. J Infect Dis. 1999;180(5):1415–23.
Greer CE, Wheeler CM, Ladner MB, et al. Human papillomavirus (HPV) type distribution and serological response to HPV type 6 virus-like particles in patients with genital warts. J Clin Microbiol. 1995;33(8):2058–63.
Castellsague X, Ghaffari A, Daniel RW, Bosch FX, Munoz N, Shah KV. Prevalence of penile human papillomavirus DNA in husbands of women with and without cervical neoplasia: a study in Spain and Colombia. J Infect Dis. 1997;176(2):353–61.
Workowski KA, Berman SM. Sexually transmitted diseases treatment guidelines, 2006. MMWR Recomm Rep. 2006;55(Rr-11):1–94.
Marra F, Ogilvie G, Colley L, Kliewer E, Marra CA. Epidemiology and costs associated with genital warts in Canada. Sex Transm Infect. 2009;85(2):111–5.
Patel H, Wagner M, Singhal P, Kothari S. Systematic review of the incidence and prevalence of genital warts. BMC Infect Dis. 2013;13:39.
Kliewer EV, Demers AA, Elliott L, Lotocki R, Butler JRG, Brisson M. Twenty-year trends in the incidence and prevalence of diagnosed anogenital warts in canada. Sex Transm Dis. 2009;36(6):380–6.
Manitoba Health HLaS. Human Papillomavirus (HPV) Information for Health Professionals. 2013.
Manitoba Health HLaS. Manitoba’s Immunication Program: Vaccines Offered Free-of-Charge (Eligibility Criteria for Pbulicly-Funded Vaccines). http://www.gov.mb.ca/health/publichealth/cdc/vaccineeligibility.html. Online accessed Jan 12, 2016. 2015.
Statistics-Canada. Census Profile. 2014.
Booth GL, Hux JE. Relationship between avoidable hospitalizations for diabetes mellitus and income level. Arch Intern Med. 2003;163(1):101–6.
Lipscombe LL, Austin PC, Manuel DG, Shah BR, Hux JE, Booth GL. Income-related differences in mortality among people with diabetes mellitus. CMAJ. 2010;182(1):E1–e17.
Marrie RA, Elliott L, Marriott J, et al. Dramatically changing rates and reasons for hospitalization in multiple sclerosis. Neurology. 2014;83(10):929–37.
Booth GLHJ. Relationship between avoidable hospitalizatins for diabetes mellitus and income level. Arch Intern Med. 2003;163:6.
Chen HS, Portier K, Ghosh K, et al. Predicting US- and state-level cancer counts for the current calendar year: Part I: evaluation of temporal projection methods for mortality. Cancer. 2012;118(4):1091–9.
Kim HJ, Fay MP, Feuer EJ, Midthune DN. Permutation tests for joinpoint regression with applications to cancer rates. Stat Med. 2000;19(3):335–51.
Kim HJ, Luo J, Kim J, Chen HS, Feuer EJ. Clustering of trend data using joinpoint regression models. Stat Med. 2014;33(23):4087–103.
Steben M, Ouhoummane N, Rodier C, Brassard P. Temporal trends in genital warts among individuals covered by the public prescription drug insurance plan in the province of Quebec, Canada, from 1998 to 2007. Journal of lower genital tract disease. 2013;17(2):147–53.
Beltrami JF, Shouse RL, Blake PA. Trends in infectious diseases and the male to female ratio: possible clues to changes in behavior among men who have sex with men. AIDS Educ Prev. 2005;17(6 Suppl B):49–59.
Canada PHAo. Report on Sexually Transmitted Infections in Canada: 2011. In: Centre for Communicable Diseases and Infection Control, Infectious Disease Prevention and Control Branch, Public Health Agency of Canada. 2014.
Fenton KA, Imrie J. Increasing rates of sexually transmitted diseases in homosexual men in Western europe and the United States: why? Infect Dis Clin North Am Jun. 2005;19(2):311–31.
Jayaraman GC, Read RR, Singh A. Characteristics of individuals with male-to-male and heterosexually acquired infectious syphilis during an outbreak in Calgary, Alberta. Canada Sex Transm Dis. 2003;30(4):315–9.
Jiamton S, Leeyaphan C, Maneeprasopchoke P, Omcharoen V. Prevalence and clinical manifestations of male patients with anogenital warts attending a sexually transmitted disease clinic prior Hpv vaccine recommendation. Southeast Asian J Trop Med Public Health. 2014;45(6):1337–43.
Mansergh G, Flores S, Koblin B, Hudson S, McKirnan D, Colfax GN. Alcohol and drug use in the context of anal sex and other factors associated with sexually transmitted infections: results from a multi-city study of high-risk men who have sex with men in the USA. Sex Transm Infect. 2008;84(6):509–11.
Mayer KH. Sexually transmitted diseases in men who have sex with men. Clin Infect Dis. 2011;53 Suppl 3:S79–83.
Poon AF, Joy JB, Woods CK, et al. The impact of clinical, demographic and risk factors on rates of HIV transmission: a population-based phylogenetic analysis in British Columbia. Canada J Infect Dis. 2015;211(6):926–35.
Steben M, LaBelle D. Genital warts: Canadians’ perception, health-related behaviors, and treatment preferences. Journal of lower genital tract disease. 2012;16(4):409–15.
Oster AM, Sternberg M, Lansky A, Broz D, Wejnert C, Paz-Bailey G. Population size estimates for Men who have Sex with Men and persons who inject drugs. J Urban Health Aug. 2015;92(4):733–43.
Aral SO, Hughes JP, Stoner B, et al. Sexual mixing patterns in the spread of gonococcal and chlamydial infections. Am J Public Health. 1999;89(6):825–33.
Aral SO, Ward H. Behavioral convergence: implications for mathematical models of sexually transmitted infection transmission. The Journal of infectious diseases. 2014;210 Suppl 2:S600–604.
Fujimoto K, Williams ML. Racial/ethnic differences in sexual network mixing: a Log-linear analysis of HIV status by partnership and sexual behavior among most at-risk MSM. AIDS Behav. 2015;19(6):996–1004.
Poynten IM, Waterboer T, Jin F, et al. Human papillomavirus types 6 and 11 seropositivity: risk factors and association with ano-genital warts among homosexual men. J Infect Jun. 2013;66(6):503–11.
Kraut AA, Schink T, Schulze-Rath R, Mikolajczyk RT, Garbe E. Incidence of anogenital warts in Germany: a population-based cohort study. BMC Infect Dis. 2010;10:360.
Sanchez DM, Pathela P, Niccolai LM, Schillinger JA. Knowledge of human papillomavirus and anal cancer among men who have sex with men attending a New York City sexually transmitted diseases clinic. Int J STD AIDS Jan. 2012;23(1):41–3.
Nadarzynski T, Smith H, Richardson D, Jones CJ, Llewellyn CD. Human papillomavirus and vaccine-related perceptions among men who have sex with men: a systematic review. Sex Transm Infect. 2014;90(7):515–23.
Graham DM, Isaranuwatchai W, Habbous S, et al. A cost-effectiveness analysis of human papillomavirus vaccination of boys for the prevention of oropharyngeal cancer. Cancer. 2015;121(11):1785–92.
Smith LM, Strumpf EC, Kaufman JS, Lofters A, Schwandt M, Levesque LE. The early benefits of human papillomavirus vaccination on cervical dysplasia and anogenital warts. Pediatrics. 2015;135(5):e1131–1140.
O’Malley KJ, Cook KF, Price MD, Wildes KR, Hurdle JF, Ashton CM. Measuring diagnoses: ICD code accuracy. Health Serv Res. 2005;40(5 Pt 2):1620–39.
Surjan G. Questions on validity of international classification of diseases-coded diagnoses. Int J Med Inform. 1999;54(2):77–95.
Roos LL, Gupta S, Soodeen RA, Jebamani L. Data quality in an information-rich environment: Canada as an example. Canadian journal on aging = La revue canadienne du vieillissement Spring . 2005;24 Suppl 1:153–70.
Fransoo R, Yogendran M, Olafson K, Ramsey C, McGowan KL, Garland A. Constructing episodes of inpatient care: data infrastructure for population-based research. BMC Med Res Methodol. 2012;12:133.
Beutner KR, Reitano MV, Richwald GA, Wiley DJ. External genital warts: report of the american medical association consensus conference. AMA expert panel on external genital warts. Clinical infectious diseases : an official publication of the Infectious Diseases Society of America. 1998;27(4):796–806.
No funding was received for this study.
“The Corresponding Author has the right to grant on behalf of all authors and does grant on behalf of all authors, an exclusive licence (or non-exclusive for government employees) on a worldwide basis to the BMJ and co-owners or contracting owning societies (where published by the BMJ on their behalf), and its Licensees to permit this article (if accepted) to be published in Sexually Transmitted Infections and any other BMJ products and to exploit all subsidiary rights, as set out in our licence.”
The author declares that they have no competing interests.
LHT and NY collaborated in the interpretation of results and the writing of the manuscript. ZN performed most of the analyses for this study and critically reviewed the manuscript. CE, JB, and NY were involved in the design of this study and critically reviewed the manuscript. All authors read and approved the final manuscript.
Tariff codes used to identify a person with AGW in the medical claims. Table S2. ICD-9/10 diagnosis codes for anogenital warts. Table S3. ICD 9 procedure codes used to assist in the identification of a person with anogenital warts. Table S4. ICD 10 procedure codes used to assist in the identification of a person with anogenital warts. Table S5. Tariff codes used to assist in the identification of a person with anogenital warts. (DOCX 46 kb)