The epidemiology and direct medical costs of diseases associated with human papillomavirus infection among men in Manitoba, Canada
==================================================================================================================================

* Christiaan H. Righolt
* Gurpreet Pabla
* Salaheddin M. Mahmud

## Abstract

**Background** There is little information on the economic burden of human papillomavirus-related diseases (HPV-RDs) among men. We used province-wide clinical, administrative and accounting databases to measure the direct medical costs of HPV infections in men in Manitoba (Canada).

**Methods** We included all males aged 9 years and older with health insurance coverage in Manitoba between January 1997 and December 2016. We identified HPV-RD patient cohorts and matched each patient to HPV-RD-free men. We estimated the net direct medical cost (excess cost of hospitalizations, outpatient visits, and prescription drugs) of patients compared to their matches for anogenital warts (AGWs) and HPV-caused cancers. We adjusted costs to 2017 Canadian dollars. For each condition, we attributed costs to HPV based on the etiological fraction caused by HPV infection.

**Results** We found that the median net direct medical cost was about $250 for AGW patients and $16,000 for invasive cancer patients. The total cost was about $49 million or $2.6 million per year. Overall, 54%-67% ($26-$33 million) was attributable to HPV infection according different estimates of the attributable fraction. The net annual attributable cost was $2.37-$2.95 per male resident and $161-$200 per male newborn. The estimated potential savings was 30% for the bivalent vaccine and 56%-60% for the quadrivalent and nonavalent vaccines.

**Conclusions** Overall, HPV’s economic burden on males remains significant, the average cost of treating all conditions attributable to HPV was about $180 per male newborn. Invasive cancer accounted for the majority of these costs.

Keywords
*   Human papillomavirus
*   direct medical cost
*   HPV-related diseases
*   HPV-related cost

## 1 Introduction

Anogenital infection with the human papillomavirus (HPV) is the most common sexually transmitted infection in Canada and worldwide [1–3]. Although HPV infection is as prevalent among men [4], the focus of most research in this field has been on women due to its strong association with cervical cancer, a more common and devastating cancer. However, persistent infection with 13 or so high-risk oncogenic HPV genotypes is known to causes most (90%) anal cancers and a significant proportion of oropharyngeal cancers (12%) in both men and women [5–7]. About 40% of penile cancers are also caused by HPV infection. In addition, “low-risk” types 6 and 11 are responsible for most external anogenital warts (AGWs) in both genders [8]. With the availability of prophylactic vaccines against HPV infection in males, it is important to understand both the epidemiologic and economic male burden associated with these human papillomavirus-related diseases (HPV-RDs).

Little is known about the epidemiology and long-term trends of anal and penile cancer among Canadian men. Because of their lower incidence, they are not regularly reported in national and provincial cancer surveillance reports. In addition, patterns of treatment of these cancers have not been characterized at the population level in Canada. Most studies that estimated HPV-related medical costs were not based on nationally representative data as information was often obtained for individual organizations (hospitals, health maintenance organizations, etc.) or from patient surveys [9–11]. Some studies used diseases incidence rates derived from mathematical modeling, these rates were then combined with estimates of average costs per case/service as derived from the literature to arrive at total expected costs [12,13]. Most studies were limited to costs associated with the prevention and treatment of one or two conditions and were focused on the cost of all HPV-RDs. We previously reported on the direct medical cost for both men and women in Manitoba and found that two-thirds of costs were for cervical screening; because of this, we reported limited data specific to males [14].

Information on the economic burden of disease is essential for planning and evaluating prevention and control programs [15]. Up-to-date cost estimates are essential to inform decisions on priorities, evaluate existing prevention and treatment programs, forecast expenditures, and estimate the impact of new interventions [16]. Estimating cost-effectiveness requires detailed knowledge of costs to facilitate decision-making on whether to expand HPV vaccine programs or to switch to using the nonavalent vaccine [17,18].

We aimed to assess the direct medical costs associated with the prevention and treatment of conditions attributable to HPV infection among men in Manitoba, Canada.

## 2 Methods

We used several Manitoba Health (MH) clinical and administrative databases to identify all males diagnosed with an HPV-RD in Manitoba between January 1997 and December 2016 and measure all costs incurred by MH in relation to the diagnosis and treatment of these conditions (*direct medical costs*).

### 2.1 Data sources

MH is the publicly funded health insurance agency providing comprehensive health insurance, including coverage for hospital and outpatient physician services, to the province’s 1.3 million residents. Coverage is universal, with no eligibility distinction based on age or income, and participation rates are very high (>99%) [19]. Insured services include hospital, physician and preventive services including vaccinations. MH maintains several centralized, administrative electronic databases that are linkable using a unique personal health identification number (PHIN). The completeness and accuracy of these databases are well established [20,21]. These databases have been used extensively in studies of disease surveillance and post-marketing evaluation of various vaccines and drugs [22–24].

The Manitoba Health Population Registry (MHPR) tracks addresses and dates of birth, death and insurance coverage for all insured persons. Since 1971, the Hospital Abstracts Database (HAD) recorded virtually all services provided by hospitals in the province, including admissions and day surgeries [20]. The data collected comprise demographic as well as diagnosis and treatment information including primary diagnosis and service or procedure codes, coded using the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) before April, 2004, and the ICD-10-CA (Canadian adaptation of the ICD-10) and the Canadian Classification of Health Interventions (CCI) afterwards. The Medical Services Database (MSD), also in operation since 1971, collects similar information, based on physician fee-for-service or shadow billing, on services provided by physicians in offices, hospitals and outpatient departments across the province [20]. The Drug Program Information Network (DPIN), in operation since 1995, records all prescription drugs dispensed to Manitoba residents, including most personal care home residents [25]. The DPIN database captures data from pharmacy claims for formulary drugs dispensed to all Manitobans even those without prescription drug insurance [25].

CancerCare Manitoba (CCMB) maintains one of the oldest population-based cancer registries in the world (Manitoba Cancer Registry [MCR], in operation since 1956). Reporting of cancer cases is mandated by provincial regulations and required for payments of physicians’ service claims. The MCR is regularly audited by the North American Association of Central Cancer Registries. The quality of cancer registration has been consistently very high [21]. Most cases are pathologically-confirmed (94% for cases registered between 2003 and 2007) and less than 2% of registrations originate from death certificates [21]. It also tracks all cancer-specific treatment using ICD-9 procedure codes prior to 2005 and ICD-10 Canadian Classification of Health Interventions codes since then. The Emergency Department Information System (EDIS) captures all visits to Winnipeg (Manitoba’s main urban center with over half the provincial population and the province’s only tertiary care hospitals) emergency departments since April 1999.

### 2.2 Eligibility and study cohort

All males aged 9 years and older who were registered with MH at any point between January 1997 and December 2016 (the *study period*) were eligible for inclusion in this study. Males joined the population cohort on the latest of their 9th birthdate, date of immigration to Manitoba or January 1, 1997. They exited the population cohort on the earliest of December 31, 2016 or the date of loss of MH coverage for any reason, including death or emigration.

To build the study (HPV-RD) cohort, we identified all eligible males who were diagnosed with an incident HPV-RD during the study period. We matched each HPV-RD patient on age (±1 year), length of MH coverage (±1 year) and region of residence to up to 3 eligible males who met all the inclusion criteria but did not have a HPV-RD at the diagnosis date of the case.

We identified males diagnosed with invasive cancer or carcinoma *in situ* of the penis, anus, oral cavity, or oropharynx by linking with the MCR. We used ICD-O-3 codes to identify these cases (Supplementary Table 1). For carcinoma *in situ* (CIS) and cancer patients, the index date was the date of diagnosis as reported in the MCR (defined as the earliest of symptom onset, laboratory diagnosis or cancer treatment date).

We obtained information on incident medically-attended AGWs from the HAD and MSD using previously validated algorithms [26]. In short, persons treated for AGWs were identified from the MSD using tariff codes specific to condyloma (Supplementary Table 2) or from the HAD using period-dependent diagnostic and procedure codes (Supplementary Table 3-4). Because patients can have multiple diagnoses of AGWs, each episode of care had to be separated by 12 months without any related claims to be considered an incident episode (Supplementary Table 6).

We added direct medical costs and the number of healthcare encounters for disease-specific periods: from 180 days prior to the diagnosis date to 730 days after diagnosis for invasive cancer, and from the diagnosis date to 365 days after diagnosis for AGWs and CIS.

### 2.3 Costs and healthcare utilization

We summed the costs, obtained from the administrative databases, of all prescription drugs (including dispensing fees), in-hospital and outpatient diagnostic and therapeutic procedures, and physician visits for each HPV-RD case and their matches. We attributed all costs to the year of diagnosis of the case and indexed costs to 2017 Canadian dollars using the all-item and prescription drugs consumer price indices. We counted all hospitalizations, physician visits, emergency department visits, and dispensed prescriptions for each HPV-RD case and their matches. We attributed all healthcare utilization to the year of diagnosis of the case.

### 2.5 Analysis

We calculated crude and age-standardized annual incidence rates of each HPV-RD using the 2016 Canada census population as the standard population. We calculated incidence rates, age-standardized when appropriate, for different subsets of the population, defined by age group, region of residence, and year of diagnosis.

We compared costs for each HPV-RD patient with their matches. We reported individual net costs (the difference in cost between the case and their specific matches) for all HPV-RDs and for each HPV-RD separately for the entire population and for specific age groups, regions of residence and years of diagnosis. To estimate the net cost attributable to HPV infection, we repeated all analyses after multiplying the cost differential of each HPV-RD by the HPV attributable fraction of that condition (Supplementary Table 7) using two sets of estimates derived from two detailed systematic reviews [27,28]. To estimate potential savings from various vaccines (bivalent, quadrivalent and nonavalent), we multiplied the cost differential of each HPV-RD by the percentage of that HPV-RD attributable to the HPV types included in the specific vaccine using estimates derived from Saraiya et al. [28] (Supplementary Table 8). We repeated the above methods for each patient’s healthcare utilization.

We used SAS 9.4 (SAS Institute, Cary, North Carolina) and Stata 16 (StataCorp, College Station, Texas) for all analysis. This study was approved by the University of Manitoba Research Ethics Board and by MH’s Health Information Privacy Committee.

## 3 Results

During the study period, AGWs were the mostly commonly treated HPV-RD with an overall rate of 149/100,000 person-years (Table 1). AGW rates were highest among 18-30 year-olds (367/100,000). Invasive cancer of the oropharynx was the most common HPV-related cancer (5.2/100,000), peaking among 51-64 year-olds at 13.4/100,000. The HPV-RD cohorts consisted of 14,784 AGW and 1,236 CIS/invasive cancer patients matched to 43,330 and 2,960 HPV-RD-free men (Table 2). About two-thirds of HPV-RD patients lived in Winnipeg, where incidence was highest in the southern suburbs. Incidence rates varied by age group, birth cohort, region of residence, and year of diagnosis (Supplementary Tables 9-33).

View this table:
[Table 1](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T1)

Table 1 
Crude incidence rates (95% confidence interval) per 100,000 person -years of HPV-related diseases among Manitoba males by age group (1997 - 2016).

View this table:
[Table 2](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T2)

Table 2 
Number of HPV-related disease patients and their matches among Manitoba males by certain socio-economic and clinical characteristics

Invasive cancer patients have the highest average direct medical cost (about $42,000; Table 3). Among these cancers, anal ($51,000) and oral cavity ($53,000) cancer patients incurred the most costs. The average cost for AGW patients was about $2,100. The median net cost (compared to matched men) of invasive cancer was about $16,000 (Table 4). The highest median net cost was for the cancer of the oral cavity ($27,000). Hospitalization was the biggest contributor to costs for invasive cancer, but outpatient costs were similar for CIS and dominant for AGWs. Overall, 65% of the net cost was due to hospitalization, the cost increased during the study period (more than doubling over 2 decades after inflation-adjustment) and varied by age group and region of residence (Supplementary Tables 34-61).

View this table:
[Table 3](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T3)

Table 3 
Average direct medical cost1 of HPV-related disease patients and their matches in 2017 Canadian dollars among Manitoba males (1997-2015)

View this table:
[Table 4](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T4)

Table 4 
Distribution of net direct medical cost 1 of HPV-related diseases (for patients compared to their matches) in 2017 Canadian dollars among Manitoba males (1997-2015)

Annual healthcare utilization (outpatient visits, hospitalizations, ED visits, and prescriptions) varied by age group, region of residence, and year of diagnosis (Supplementary Tables 62-153). Overall, HPV-RD patients had about 11 outpatient visits and 13 prescriptions per year. These patients had less than 1 annual hospitalization or ED visit. As expected, healthcare utilization was highest for invasive cancer patients.

During 1997–2015, the total net cost of HPV-RDs was about $49 million or about $2.6 million per year. Overall, about 54% (about $26 million) of the total net direct medical costs were attributed to treatment of HPV-caused diseases according to attributable fraction estimates by Volesky *et al*. [27] and about 67% (about $33 million) according to estimates by Saraiya *et al*. [28] (Table 5). The difference was mostly due to different attributable fractions of invasive cancer of the oral cavity. The net attributable cost of AGWs was about $14 million and invasive diseases accounted for $12-$19 million. About half the net attributable cost of HPV-caused cancer was for the treatment of invasive cancer of the oropharynx.

View this table:
[Table 5](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T5)

Table 5 
Total net direct medical cost1 of HPV-related diseases (for patients compared to their matches) in 2017 Canadian dollars among Manitoba males according to different estimates of the attributable fraction (1997-2015)

The net attributable annual cost was $2.37-$2.95 per male resident and $161-$200 per male newborn over the course of study period (Table 6), although these costs varied by year and were higher at the end of the study period compared to the start.

View this table:
[Table 6](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T6)

Table 6 
Total net direct medical cost1 of HPV-related diseases (for patients compared to their matches) in 2017 Canadian dollars among Manitoba males by calendar year according to different estimates of the attributable fraction

The estimated potential saving from various HPV vaccines was about 30% (about $15 million) for the bivalent vaccine, 56% (about $28 million) for the quadrivalent vaccine, and 60% (about $29 million) for the nonavalent vaccine, according to attributable fractions estimated by Saraiya *et al*. [28] (Table 7). The major difference in these amounts is the exclusion of AGW-causing HPV types in the bivalent vaccine.

View this table:
[Table 7](http://medrxiv.org/content/early/2020/12/04/2020.12.02.20242891/T7)

Table 7 
Total net direct medical cost1 of HPV-related diseases (for patients compared to their matches) in 2017 Canadian dollars among Manitoba males attributable to the HPV types included in different HPV vaccines (1997-2015)

## 4 Discussion

We found that the average cost of treating all conditions attributable to HPV was about $2.60 per male resident per year or about $180 per male newborn (in 2017 Canadian dollars). Invasive cancer accounted for about half of these costs and AGW for roughly the other half. The costs of treating HPV-related cancers was about $16,000; others have estimated the annual costs for common chronic conditions (in 2006 Canadian dollars) at $3,400 for coronary heart disease, $3,800 for diabetes, and $6,000 for stroke [29] (corresponding to $4,100, $4,500, and $7,200 respectively in 2017 Canadian dollars).

The literature on the total economic burden of HPV infections is quite limited, especially when it comes to males. The costs of HPV infections are much higher than any other sexually transmitted infection, except for HIV/AIDS [30]. Differences in HPV-RDs examined, study periods, populations, currencies, methods, and assumptions make it difficult to compare individual studies [10]. Costs are often estimated for men and women combined. For example, the cost per case for oropharyngeal and anal cancer for both genders was estimated as $43,200 and $36,200 (in 2010 US dollars) from a literature review [12]. The same review estimated a cost of $19,800 (in 2010 US dollars) for penile cancer and $810 for AGW (men and women combined) [12]. The costs per anal cancer patient in the US was estimated as $60,913 (in 2014 US dollars) for 66+ Medicare-covered men [31] and $127,531 (2015 US dollars) for <65 commercially insured patients (men and women) [32]. The annual purchasing power parity-adjusted per-capita health expenditure over the last 15 years averaged about $3,800 in Canada and $7,400 in the US in 2010 US dollars [33]. This means the differences between Canada and the US could be explained by the different healthcare systems.

European health systems differ from Canada, so cost estimates from European studies are hard to compare to our results. An Italian study pegged the medical cost attributable to HPV 6, 11, 16, and 18 infections in men at €112.9 million (in 2011 Euros) based on a literature review [34]; they estimated €10,000 in medical cost per penile cancer case and €470 per male AGW case [34]. A French study found that over 95% of costs of HPV-related cancers in men was attributable to head and neck cancer [35]; we found 85% of cancer costs were for cancer of the oral cavity and oropharynx. The French cost for specific HPV-related cancers were estimated between €3,800 and €6,800 (in 2006/07 Euros) for men and women combined [35]. In Spain the male cost of HPV-related cancers ranged between €7,000 and €9,000 (in 2017 Euros), depending on the specific site [36]. These European costs per case all seem lower than the costs we reported from Manitoba, although we used a different methodology.

### 4.1 Strengths and limitations

A major strength of this study is the availability of high quality, comprehensive, population-based healthcare utilization databases and disease registries maintained by the sole health insurer in Manitoba. Since the original purpose of these databases was tracking payment for healthcare providers, our cost estimates were based on direct measurement of actual paid amounts rather than on using average cost estimates derived from micro-costing studies or surveys.

We may not have identified all HPV-RD patients, because of the limitations of retrospective studies and administrative databases. Patients may have been misdiagnosed or the diagnosis date may not have been correct. Disease misclassification could also have occurred due to under-reporting and coding errors (e.g., using the wrong ICD code). Because of the high accuracy of the MCR [21] and the prior validation study of AGW codes [26], we believe disease misclassification was rare in this study.

Our study is limited to the costs that are recorded in the MH datasets, which excludes costs of treating cancer patients in the community, e.g., home intravenous chemotherapy. We do not believe that these excluded costs were significant, because these programs were only recently introduced. We also did not include indirect costs (e.g., loss of employment) of costs incurred by patients (e.g., unpaid time off work), the total costs to all payers is certainly higher than the direct medical costs we report.

The HPV-attributable cost estimates are limited by the accuracy of the estimates of HPV-attributable fractions of disease. We used different estimates of attributable fractions from the literature to account for this uncertainty. Although this leads to roughly similar estimates, analysis based on our results should account for the range of results we report.

## 5 Conclusions

Overall, HPV’s economic burden on males remains significant, the average cost of treating all conditions attributable to HPV was about $2.60 per male resident per year or about $180 per male newborn. Invasive cancer accounted for about half of these costs.

## Supporting information

Supplement [[supplements/242891_file02.pdf]](pending:yes)

## Data Availability

Data used in this article was derived from administrative health and social data as a secondary use. The data was provided under specific data sharing agreements only for approved use at the Manitoba Centre for Health Policy (MCHP). The original source data is not owned by the researchers or MCHP and as such cannot be provided to a public repository. The original data source and approval for use has been noted in the acknowledgments of the article. Where necessary, source data specific to this article or project may be reviewed at MCHP with the consent of the original data providers, along with the required privacy and ethical review bodies.

## Data Availability Statement

Data used in this article was derived from administrative health and social data as a secondary use. The data was provided under specific data sharing agreements only for approved use at the Manitoba Centre for Health Policy (MCHP). The original source data is not owned by the researchers or MCHP and as such cannot be provided to a public repository. The original data source and approval for use has been noted in the acknowledgments of the article. Where necessary, source data specific to this article or project may be reviewed at MCHP with the consent of the original data providers, along with the required privacy and ethical review bodies.

## Conflict of interest statement

CHR has received an unrestricted research grant from Pfizer for an unrelated study. SMM has received unrestricted research grants from Merck, GlaxoSmithKline, Sanofi Pasteur, Pfizer and Roche-Assurex for unrelated studies. SMM has received fees as an advisory board member for Sanofi Pasteur. None of the other authors has any conflicts of interest to disclose.

## Funding sources

This work was supported by the Merck Investigator Studies Program with a grant to the Biomedical Association of Manitoba (BAM). SMM’s work is supported, in part, by funding from the Canada Research Chair Program. The sponsors had no role in the design or conduct of the study, including but not limited to, data identification, collection, management, analysis and interpretation, or preparation, review, or approval of the final results. The opinions presented in the report do not necessarily reflect those of the sponsors.

## Author’s contributions

CHR and SMM designed the study, GP and CHR analyzed the data, CHR, GP and SM interpreted the data, CHR and GP drafted the manuscript. All authors revised the manuscript and approved the final version.

## Acknowledgements

The authors acknowledge the Manitoba Centre for Health Policy for use of data contained in the Population Health Research Data Repository under project # 2019-008 (HIPC # 2018/2019-23; REB # HS21956 (H2018:265); RRIC #2018-028). The results and conclusions are those of the authors and no official endorsement by the Manitoba Centre for Health Policy, Manitoba Health, or the Winnipeg Regional Health Authority is intended or should be inferred.

## List of abbreviations

AGW
:   anogenital wart
CCMB
:   CancerCare Manitoba
CIS
:   carcinoma *in situ*
DPIN
:   Drug Program Information Network
ED
:   emergency department
EDIS
:   Emergency Department Information System
HAD
:   Hospital Abstracts Database
HPV
:   human papillomavirus
HPV-RD
:   human papillomavirus-related disease
MCHP
:   Manitoba Centre for Health Policy
MCR
:   Manitoba Cancer Registry
MH
:   Manitoba Health
MHPR
:   Manitoba Health Population Registry
MSD
:   Medical Services Database
PHIN
:   Personal Health Identification Number
WRHA
:   Winnipeg regional health authority

*   Received December 2, 2020.
*   Revision received December 2, 2020.
*   Accepted December 4, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/)

## References

1.  1.Baseman JG, Koutsky LA. The epidemiology of human papillomavirus infections. J Clin Virol. 2005;32 Suppl 1:S16–24.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jcv.2004.12.008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15753008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228272000003&link_type=ISI) 

2.  2.Weinstock H, Berman S, Cates W. Sexually transmitted diseases among American youth: incidence and prevalence estimates, 2000. Perspect Sex Reprod Health. 2004;36:6–10.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1363/3600604&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14982671&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000189252700001&link_type=ISI) 

3.  3.Franco EL, Duarte-Franco E, Ferenczy A. Cervical cancer: epidemiology, prevention and the role of human papillomavirus infection. CMAJ. 2001;164:1017–25.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY21haiI7czo1OiJyZXNpZCI7czoxMDoiMTY0LzcvMTAxNyI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIwLzEyLzA0LzIwMjAuMTIuMDIuMjAyNDI4OTEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

4.  4.Dunne EF, Nielson CM, Stone KM, Markowitz LE, Giuliano AR. Prevalence of HPV infection among men: a systematic review of the literature. J Infect Dis. 2006;194:1044–1057.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/507432&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16991079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000240548500005&link_type=ISI) 

5.  5.Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV, et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348:518–27.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa021641&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12571259&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180801400005&link_type=ISI) 

6.  6.Forman D, de Martel C, Lacey CJ, Soerjomataram I, Lortet-Tieulent J, Bruni L, et al. Global burden of human papillomavirus and related diseases. Vaccine. 2012;30 Suppl 5:F12–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2012.07.055&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23199955&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000312974900004&link_type=ISI) 

7.  7.Pinkbook | HPV | Epidemiology of Vaccine Preventable Diseases | CDC [Internet]. 2019 [cited 2019 Oct 30]. Available from: [https://www.cdc.gov/vaccines/pubs/pinkbook/hpv.html](https://www.cdc.gov/vaccines/pubs/pinkbook/hpv.html)
    
    

8.  8.Franco EL, Mahmud SM, Tota J, Ferenczy A, Coutlée F. The Expected Impact of HPV Vaccination on the Accuracy of Cervical Cancer Screening: The Need for a Paradigm Change. Arch Med Res. 2009;40:478–85.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.arcmed.2009.06.003&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19853188&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000271772400008&link_type=ISI) 

9.  9.Insinga R, Dasbach E, Elbasha E. Assessing the annual economic burden of preventing and treating anogenital human papillomavirus-related disease in the US. Pharmacoeconomics. 2005;23:1107–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2165/00019053-200523110-00004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16277547&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233761100004&link_type=ISI) 

10. 10.Lairson DR, Fu S, Chan W, Xu L, Shelal Z, Ramondetta L. Mean direct medical care costs associated with cervical cancer for commercially insured patients in Texas. Gynecol Oncol. 2017;
    
    

11. 11.Insinga RP, Glass AG, Rush BB. The health care costs of cervical human papillomavirus--related disease. Am J Obstet Gynecol. 2004;191:114–20.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajog.2004.01.042&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15295351&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223139100017&link_type=ISI) 

12. 12.Chesson HW, Ekwueme DU, Saraiya M, Watson M, Lowy DR, Markowitz LE. Estimates of the annual direct medical costs of the prevention and treatment of disease associated with human papillomavirus in the United States. Vaccine. 2012;30:6016–9.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2012.07.056&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22867718&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000309025300003&link_type=ISI) 

13. 13.Helms LJ, Melnikow J. Determining costs of health care services for cost-effectiveness analyses: the case of cervical cancer prevention and treatment. Med Care. 1999;37:652–61.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00005650-199907000-00005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10424636&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000081494200005&link_type=ISI) 

14. 14.Righolt CH, Pabla G, Mahmud SM. The Direct Medical Costs of Diseases Associated with Human Papillomavirus Infection in Manitoba, Canada. Appl Health Econ Health Policy. 2018;16:195–205.
    
    

15. 15.Erickson LJ, De Wals P, Farand L. An analytical framework for immunization programs in Canada. Vaccine. 2005;23:2470–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2004.10.029&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15752833&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227898800008&link_type=ISI) 

16. 16.Insinga RP, Dasbach EJ, Elbasha EH. Assessing the annual economic burden of preventing and treating anogenital human papillomavirus-related disease in the US: analytic framework and review of the literature. Pharmacoeconomics. 2005;23:1107–22.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2165/00019053-200523110-00004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16277547&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233761100004&link_type=ISI) 

17. 17.Petrosky E, Bocchini Jr JA, Hariri S, Chesson H, Curtis CR, Saraiya M, et al. Use of 9-Valent Human Papillomavirus (HPV) Vaccine: Updated HPV Vaccination Recommendations of the Advisory Committee on Immunization Practices. MMWR Morb Mortal Wkly Rep. 2015;64:300–4.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25811679&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 

18. 18.Joura EA, Giuliano AR, Iversen O-E, Bouchard C, Mao C, Mehlsen J, et al. A 9-valent HPV vaccine against infection and intraepithelial neoplasia in women. N Engl J Med. 2015;372:711–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMoa1405044&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25693011&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 

19. 19.Martens P, Nickel N, Lix L, Turner D, Prior H, Walld R, et al. The cost of smoking: a Manitoba study. Winnipeg: Manitoba Centre for Health Policy; 2015.
    
    

20. 20.Roos LL, Mustard CA, Nicol JP, McLerran DF, Malenka DJ, Young TK, et al. Registries and administrative data: organization and accuracy. Med Care. 1993;31:201–12.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00005650-199303000-00002&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8450678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1993KU34500002&link_type=ISI) 

21. 21.Robinson JR, Young TK, Roos LL, Gelskey DE. Estimating the burden of disease: comparing administrative data and self-reports. Med Care. 1997;35:932–47.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/00005650-199709000-00006&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9298082&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XV88500006&link_type=ISI) 

22. 22.Mahmud S, Hammond G, Elliott L, Hilderman T, Kurbis C, Caetano P, et al. Effectiveness of the pandemic H1N1 influenza vaccines against laboratory-confirmed H1N1 infections: Population-based case–control study. Vaccine. 2011;29:7975–81.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.vaccine.2011.08.068&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21884747&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000296683700015&link_type=ISI) 

23. 23.Singh H, Mahmud SM, Turner D, Xue L, Demers AA, Bernstein CN. Long-Term Use of Statins and Risk of Colorectal Cancer: A Population-Based Study. Am J Gastroenterol. 2009;104:3015–23.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ajg.2009.574&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19809413&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000272553100020&link_type=ISI) 

24. 24.Mahmud SM, Van Caeseele P, Hammond G, Kurbis C, Hilderman T, Elliott L. No Association between 2008–09 Influenza Vaccine and Influenza A (H1N1) pdm09 Virus Infection, Manitoba, Canada, 2009. Emerg Infect Dis. 2012;18:801.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid1805.111596&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22516189&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 

25. 25.Kozyrskyj AL, Mustard CA. Validation of an electronic, population-based prescription database. Ann Pharmacother. 1998;32:1152–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1345/aph.18117&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9825079&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000076959800003&link_type=ISI) 

26. 26.Kliewer EV, Demers AA, Elliott L, Lotocki R, Butler JR, Brisson M. Twenty-year trends in the incidence and prevalence of diagnosed anogenital warts in Canada. Sexually transmitted diseases. 2009;36:380–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/OLQ.0b013e318198de8c&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19556932&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000266361700011&link_type=ISI) 

27. 27.Volesky KD, El-Zein M, Franco EL, Brenner DR, Friedenreich CM, Ruan Y. Cancers attributable to infections in Canada. Prev Med. 2019;122:109–17.
    
    

28. 28.Saraiya M, Unger ER, Thompson TD, Lynch CF, Hernandez BY, Lyu CW, et al. US Assessment of HPV Types in Cancers: Implications for Current and 9-Valent HPV Vaccines. J Natl Cancer Inst [Internet]. 2015 [cited 2019 Oct 31];107. Available from: [https://academic.oup.com/jnci/article/107/6/djv086/872092](https://academic.oup.com/jnci/article/107/6/djv086/872092)
    
    

29. 29.Finlayson GS, Ekuma O, Yogendran M, Burland E, Forget E. The Additional Cost of Chronic Disease in Manitoba. Winnipeg; 2010.
    
    

30. 30.Owusu-Edusei K, Chesson HW, Gift TL, Tao G, Mahajan R, Ocfemia MC, et al. The estimated direct medical cost of selected sexually transmitted infections in the United States, 2008. Sex Transm Dis. 2013;40:197–201.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/OLQ.0b013e318285c6d2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23403600&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 

31. 31.Deshmukh AA, Zhao H, Franzini L, Lairson DR, Chiao EY, Das P, et al. Total Lifetime and Cancer-related Costs for Elderly Patients Diagnosed With Anal Cancer in the United States. American Journal of Clinical Oncology. 2018;41:121–127.
    
    

32. 32.Wu C-F, Xu L, Fu S, Peng H-L, Messick CA, Lairson DR. Health Care Costs of Anal Cancer in a Commercially Insured Population in the United States. J Manag Care Spec Pharm. Academy of Managed Care Pharmacy; 2018;24:1156–64.
    
    

33. 33.OECD. Health expenditure indicators.
    
    

34. 34.Baio G, Capone A, Marcellusi A, Mennini FS, Favato G. Economic Burden of Human Papillomavirus-Related Diseases in Italy. PLoS One. 2012;7:e49699.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0049699&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23185412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F12%2F04%2F2020.12.02.20242891.atom) 

35. 35.Borget I, Abramowitz L, Mathevet P. Economic burden of HPV-related cancers in France. Vaccine. 2011;29:5245–9.
    
    

36. 36.López N, Torné A, Franco A, San-Martin M, Viayna E, Barrull C, et al. Epidemiologic and economic burden of HPV diseases in Spain: implication of additional 5 types from the 9-valent vaccine. Infect Agents Cancer. 2018;13:15.