Skip to main content
  • Research article
  • Open access
  • Published:

Impact of a population-based HPV vaccination program on cervical abnormalities: a data linkage study



Australia was one of the first countries to introduce a publicly funded national human papillomavirus (HPV) vaccination program that commenced in April 2007, using the quadrivalent HPV vaccine targeting 12- to 13-year-old girls on an ongoing basis. Two-year catch-up programs were offered to 14- to 17- year-old girls in schools and 18- to 26-year-old women in community-based settings. We present data from the school-based program on population-level vaccine effectiveness against cervical abnormalities in Victoria, Australia.


Data for women age-eligible for the HPV vaccination program were linked between the Victorian Cervical Cytology Registry and the National HPV Vaccination Program Register to create a cohort of screening women who were either vaccinated or unvaccinated. Entry into the cohort was 1 April 2007 or at first Pap test for women not already screening. Vaccine effectiveness (VE) and hazard ratios (HR) for cervical abnormalities by vaccination status between 1 April 2007 and 31 December 2011 were calculated using proportional hazards regression.


The study included 14,085 unvaccinated and 24,871 vaccinated women attending screening who were eligible for vaccination at school, 85.0% of whom had received three doses. Detection rates of histologically confirmed high-grade (HG) cervical abnormalities and high-grade cytology (HGC) were significantly lower for vaccinated women (any dose) (HG 4.8 per 1,000 person-years, HGC 11.9 per 1,000 person-years) compared with unvaccinated women (HG 6.4 per 1,000 person-years, HGC 15.3 per 1,000 person-years) HR 0.72 (95% CI 0.58 to 0.91) and HR 0.75 (95% CI 0.65 to 0.87), respectively. The HR for low-grade (LG) cytological abnormalities was 0.76 (95% CI 0.72 to 0.80). VE adjusted a priori for age at first screening, socioeconomic status and remoteness index, for women who were completely vaccinated, was greatest for CIN3+/AIS at 47.5% (95% CI 22.7 to 64.4) and 36.4% (95% CI 9.8 to 55.1) for women who received any dose of vaccine, and was negatively associated with age. For women who received only one or two doses of vaccine, HRs for HG histology were not significantly different from 1.0, although the number of outcomes was small.


A population-based HPV vaccination program in schools significantly reduced cervical abnormalities for vaccinated women within five years of implementation, with the greatest vaccine effectiveness observed for the youngest women.

Peer Review reports


Randomized clinical trials have shown that prophylactic human papillomavirus (HPV) vaccines are highly effective at preventing infection with HPV types 16 and 18, which cause about 70 to 80% of cervical cancers worldwide [1, 2] and in preventing related cervical abnormalities, which are the precursors to cervical cancer [3, 4]. HPV16 and 18 are detected in about half of all high grade cervical intraepithelial lesions (CIN2/3), and in an even higher percentage of CIN2/3 diagnosed in women under 30 years of age [2, 57]. Two vaccines are commercially available - a bivalent vaccine targeting HPV types 16 and 18 and a quadrivalent vaccine which also targets HPV types 6 and 11, responsible for 90% of genital warts. Types 16/18/6/11 have been detected in about 30% of low-grade cervical cytological abnormalities, reflecting that these low-grade lesions indicate acute HPV infection, which can be caused by any of the 40 genital HPV types [8]. The trials showed greatest efficacy for women naïve to HPV at the time of vaccination, as the vaccines do not have any therapeutic effect on targeted HPV types.

These vaccines have been widely adopted in high income countries since 2007 with most countries targeting girls prior to commencement of sexual activity, between 9 and 13 years of age, with varying catch-up programs for older adolescent girls [9]. Australia was one of the first countries to offer a national, publicly funded vaccination program, using a three-dose schedule of the quadrivalent vaccine Gardasil©. Commencing in April 2007, the Australian vaccination program targeted 12- and 13-year-old girls in schools and offered catch-up vaccination for 14- to 17-year-old girls in schools and 18- to 26-year-old women in community based settings.

Australia is in a unique position to evaluate the impact of the HPV vaccine at the population level. The recommended starting age for cervical screening is 18 years or 2 years after the onset of sexual activity, whichever is later, and thus screening already overlaps with the vaccination program. Screening is recommended every two years using conventional cytology until the age of 69 years. Australia’s State based Pap test registers systematically record cervical cytology and histology results and send reminder and follow-up letters to women and their practitioners according to national guidelines. Participation in both screening and vaccination programs is high: catch-up HPV vaccine coverage in Victorian females 12 to 17-years old was 86/82/75% for doses 1/2/3 (JMLB personal communication) and 77% of 20 to 24-year-olds attended at least one screening between 2007 and 2011 [10].

In this retrospective cohort study, we linked data from the Victorian Cervical Cytology Registry (VCCR) and the National HPV Vaccination Program Register (NHVPR) and evaluated the effectiveness of the HPV vaccine against cervical abnormalities in a screening population of women eligible for vaccination in the school-based cohorts (aged 17 or younger in 2007). We focus on this cohort for two reasons: (1) effectiveness in this age group will most closely approximate the eventual impact of vaccination in the ongoing target group of girls prior to sexual activity and (2) vaccination status in the school cohorts is more completely notified than in the adult women catch-up program, where notification was not compulsory and under-notification has been estimated at 10 to 15% [11].


Data collection

The VCCR captures cervical screening results for all women resident in Victoria, Australia, which has a population of more than 2.7 million women, and receives data on cervical cytology and histopathology directly from laboratories in a timely manner. Fewer than 1% of women request to opt-off the VCCR [10]. Details regarding the information held on the VCCR have been published previously [10]. The VCCR supports the national screening policy by sending reminders and contacting women and their practitioners when results of indicated follow-up tests are not received by the register. Australian guidelines recommend referral to colposcopy for all women with a high-grade abnormality suspected or diagnosed on cytology and for all glandular abnormalities. It is recommended that women with low grade cytology have a repeat Pap smear at 12 months, except in women over 30 years of age with an inadequate screening history or when two consecutive low-grade abnormalities have been reported [12].

The NHVPR was established in June 2008 and records HPV vaccination doses administered nationwide. Details on the establishment of the NHVPR, and initial coverage in both the schools program and the catch-up program have been published elsewhere [11, 13, 14]. In brief, doses administered through schools are received from school immunization providers, which are predominantly local councils and health departments, and notification from these providers is believed to be virtually complete. Doses administered in the school programs prior to register establishment were also systematically uploaded to the register.

Data linkage

Cross-referencing of data against cervical cytology registers for monitoring of the effectiveness of HPV vaccine is allowed for by legislation enabling operation of the NHVPR. However, because the VCCR was established far earlier, there is no similar provision for the VCCR and thus individual identifiers could not be used for the data linkage in this study.

A deterministic linkage procedure was undertaken to link records from the VCCR, which included Victorian women who are age-eligible for the HPV vaccine as per product indication (date of birth on or after 1 July 1962), to the NHVPR which included all Australian women who received the HPV vaccine prior to 1 January 2012 (Figure 1).

Figure 1
figure 1

Data linkage process and exclusions for analysis.

At the time of extraction, data from the VCCR and NHVPR were de-identified in the same manner. Randomly generated record ID numbers were assigned to women in each register. The Data Linkage Unit generated a series of 22 linkage keys using combinations of variables, such as selected letters from given name and surname, perturbed date of birth, postcode and parts of the Medicare number. The quality of these linkage keys was assessed by calculating the percentage of unique combinations each would be likely to produce. The first 16 linkage keys produced a set of unique combinations with the best balance between incorrectly accepting a given record pair as the same person versus failing to accept a given record pair as the same person, and these were used to link records in a deterministic manner.

Study design and inclusion criteria

A retrospective cohort was constructed of women aged 17 or younger in 2007 who had a Pap test recorded on the VCCR during the study period, 1 April 2007 (the date the HPV vaccination program commenced) to 31 December 2011. Women were counted as at risk of a diagnosis of a cervical abnormality from the time they commenced cervical screening, and were entered into the cohort at their first Pap test (or on 1 April 2007 if their first Pap test was prior to that time). We found that our analyses were robust when the entry criteria to the cohort for those screening prior to April 2007 were altered to include person-time from their first Pap test after 1 April 2007. Women were followed until the outcome of interest, date of death, hysterectomy or the end of the study period.

Assignment of vaccination status

Unvaccinated women were those who had no doses of HPV vaccine recorded on the NHVPR; vaccinated women were those who received any doses of HPV vaccine. Completely vaccinated women were those vaccinated with three doses compliant with the vaccine schedule (and, if doses were closer than recommended minimum intervals, where no further doses were recommended according to Chief Medical Officer guidelines. Women with further doses needed according to the guidelines were excluded from the analyses (n = 25) [15]. Person-time and events were allocated to each “woman-dose” according to her vaccination dates. For unvaccinated women, person-time commenced from the date of first Pap test as above. For vaccinated women, a new vaccination status (number of doses) was assigned on the date of that dose and person-time at that dose number was included until receipt of a further dose (if any). Where screening and vaccination occurred on the same day, if a screening-related outcome was diagnosed on that test, the woman’s vaccination status for that outcome was assigned to her vaccination status during the period immediately prior to that date (dose delivered -1); hence, prevalent abnormalities were assigned to the relevant vaccination status.

We chose the dose assignment method, with no time-lag until abnormalities are ‘counted’ after vaccination, because we wanted to reflect apparent incident abnormalities as they occurred in the real world (that is, what the patient and clinician would experience after vaccination as a measure of vaccine effectiveness). However, we recognize that abnormalities diagnosed during the vaccination course are predominantly caused by HPV infection acquired prior to the onset of vaccination. Therefore, we also evaluated the cohort using a measure closer to vaccine efficacy assessment, excluding time for vaccinated women where they were in the process of being vaccinated, that is, censoring time occurring between dose 1 and dose 3. This second method (presented in the Additional file 1: Table S1) counts person-time according to the woman’s final vaccination status and thus only includes person-time for partial vaccination where this was the woman’s final vaccination status. This will more closely reflect the actual efficacy of receipt of only one or two doses.

Definition of outcomes

The primary outcome was histologically confirmed high-grade (HG) cervical disease (CIN2+/AIS), defined as CIN2, CIN3, adenocarcinoma in situ or mixed CIN3/AIS. Histological and cytological outcomes were assigned according to categorization used by the AIHW [12] and Australian Standardised Modified Bethesda System, respectively [16]. For all outcomes, a woman’s first relevant abnormality or her first in two years with at least two negative cytology tests in between, were counted. Abnormalities diagnosed on the first test were included as incident abnormalities.

Statistical methods

Analyses of demographic and exposure characteristics of women in the cohort by vaccination status used the Mann-Whitney U test for ordinal variables and the Pearson chi-square test for nominal variables. Detection rates were calculated as the number of events per 1,000 person-years at risk and vaccine effectiveness was (1-rate ratio) × 100.

Cox proportional hazards regression was used to estimate hazard ratios (HR) and 95% confidence intervals (CI), and included a priori the other predictors: age at first screening (as a proxy for age at sexual debut), socioeconomic status and remoteness. Remoteness areas were assigned based on postcode of residence and as defined according to the Australian Standard Geographic Classification. Socioeconomic status was assigned by postcode of residence according to the Australian Bureau of Statistic’s Socio-Economic Indexes for Areas (SEIFA) Index of Relative Socio-Economic Disadvantage [17]. Lower socioeconomic status is associated with poorer screening participation and higher cervical cancer incidence in Australia [12]. Individuals were observed until detection of outcome or 31 December 2011, whichever occurred first. All data were analyzed using the Proportional hazards regression (PHREG) procedure using SAS Institute Inc software, Cary, NC, USA [18].

Ethics approval was obtained from the Department of Health and Ageing and the Australian Institute of Health and Welfare’s Human Research and Ethics Committees. Approval for use of NHVPR data was given by the Department of Health and Ageing, the data custodian and for the VCCR data by the Victorian Department of Health.


Between 1 April 2007 and 31 December 2011, 24,871 women aged between 12 and 17 years who were vaccinated against HPV had commenced cervical screening. Of these women, 21,151 (85.0%) were completely vaccinated and 3,690 women had received one or two doses of vaccine. There were 14,085 unvaccinated women of the same age who had commenced cervical screening. The follow-up period was a maximum of 4.8 years with an average of 1.5 years for both vaccinated women and unvaccinated women. There was no difference in mean age at entry into the cohort or mean age at first screen between vaccinated and unvaccinated women, with almost all first Pap tests occurring after April 2007, Table 1. Because the youngest girls at the time of vaccination were only just becoming eligible for cervical screening at the time of analysis, most women in the study were vaccinated between 14 and 17 years of age, with a mean age at vaccination of 15.7 years. Seventy five percent of women in this analysis who completed the vaccine course received their first dose in 2007, at the commencement of the Australian HPV vaccination program, with 20% commencing in 2008 (predominantly those aged 14 years and under in 2007). Unvaccinated women were more likely to be from a major city and of lower socioeconomic status (SES) than vaccinated women. Partially vaccinated women (one or two doses) were younger at first screen than either unvaccinated women or vaccinated women, much older at vaccine commencement and of lower socioeconomic status than completely vaccinated women (P <0.05).

Table 1 Summary of descriptive characteristics of cohort

The number and rate of cervical abnormalities detected by vaccination status as per our primary analysis (dose assignment method for assessing vaccine effectiveness) is shown in Table 2. A lower risk of any histologically confirmed HG cervical abnormality was observed for vaccinated women (any dose) compared with unvaccinated women with a hazard ratio of 0.72 (95% CI 0.58 to 0.91), after adjusting for age at first screening, SES and remoteness. This effect was strongest for completely vaccinated women; there was no significant reduction among those partially vaccinated (Table 2), but the number of outcomes was small. Similarly, detection rates of CIN3/AIS and CIN2+ were significantly lower for vaccinated women compared to unvaccinated women during the study period. There was a reduced risk of LG cytological abnormalities for women who received one or two doses of vaccine HR 0.66 (95% CI 0.60 to 0.72) compared with unvaccinated women. However, this result was sensitive to the censoring of time during the vaccination course for fully vaccinated women using our secondary analytic method (more similar to trial efficacy calculations) (see Additional file 1: Table S1). There were only three women in this cohort with outcomes of AIS, one occurred in a vaccinated woman and two in unvaccinated women.

Table 2 Number and rate of cervical abnormalities for completely vaccinated, partially vaccinated and unvaccinated women

Figure 2a-f shows the detection rates by age in 2007 and vaccination status for each outcome of interest. Detection rates for any histological outcomes increased with age and were lower for women who had received any dose of HPV vaccine compared with unvaccinated women at all ages. For cytological outcomes, rates were also lower for vaccinated women at all ages as compared with unvaccinated women for LG cytology. There was a suggestion of a reduced risk of HG histology for younger women who received one or two doses of vaccine rather than none as compared with older women but the numbers were small (data not presented).

Figure 2
figure 2

Cervical abnormalities rates for vaccinated and unvaccinated women, by age in 2007. High-grade histology per 1,000 women years, by age in 2007. (a) (b) CIN 3/AIS histology per 1,000 women years, by age in 2007. (c) CIN 2 histology per 1,000 women years, by age in 2007. (d) CIN I histology per 1,000 women years, by age in 2007. (e) Low-grade cytology per 1,000 women years, by age in 2007. (f) High grade cytology per 1,000 women years, by age in 2007. All high grade histology is defined as CIN2, CIN3, AIS and mixed CIN3/AIS. Age is presented in years.

Vaccine effectiveness (VE), adjusted for remoteness, SES and age at first Pap test, was highest for CIN3/AIS at 47.5% (95% CI 22.7 to 64.4) for women who were completely vaccinated (compared to no doses), (Figure 3a) and was slightly lower for women who received any dose of vaccine 36.4% (95% CI 9.8 to 55.1), Figure 3b). VE for histological outcomes was highest in younger women (that is, those who were less likely to be sexually active before vaccination) declining thereafter with increasing age. For most age groups, VE was lowest for CIN1 outcomes, although it is important to note that in Australia only a small proportion of LG cytology lesions lead to further investigation [16]. For HG cytology, adjusted VE for all ages combined for any dose of vaccine was 24.9% (95% CI 13.3 to 35.0), similar to that for LG cytology, adjusted VE 23.9% (95% CI 19.6 to 27.9).

Figure 3
figure 3

HPV vaccine effectiveness for cervical histological outcome, by age in 2007, for (a) completed vaccine course, for (b) any vaccine dose. All high-grade histology is defined as CIN2, CIN3, AIS and mixed CIN3/AIS. Vaccine effectiveness is defined as (1-adjusted hazard rate) x 100. Age in years, as of 2007.


Our results show that quadrivalent HPV vaccination substantially reduced cervical abnormalities for women receiving any HPV vaccine dose through a school-based program, within five years of implementation. Vaccine effectiveness against high-grade histological disease was higher than that observed in the intent-to-treat (ITT) analyses of the clinical trials of the quadrivalent HPV vaccine, which vaccinated women between the ages of 16 and 25, and was highest in younger women, who as a group are least likely to have been sexually active prior to vaccination. Vaccine effectiveness was greatest for outcomes most strongly associated with vaccine related HPV types 16 and 18, in particular CIN3/AIS, and lowest for low-grade cytological abnormalities for which the majority of disease is accounted for by non-vaccine HPV types. The reduction in risk of cervical abnormalities was restricted to women who received three doses of vaccine.

The follow-up time in our study was five years and further follow-up will likely impact upon the effectiveness of the vaccine. Our data confirm that VE diminishes with age of administration in catch-up programs, which is a proxy for the onset of sexual activity. The ongoing vaccination programs in Australia and elsewhere target girls prior to onset of sexual debut and at the population level higher reductions will be observed as this large cohort reaches screening age, as a greater proportion of vaccinated women will have been unexposed to HPV at the time of vaccination. The VE was higher in the youngest women in our study than observed in the ITT analyses of either the quadrivalent or bivalent vaccine trials [1921]. The VE from ITT analyses, at an average of 3.6 years follow-up, for the quadrivalent vaccine irrespective of HPV type was 16.4% 95% CI (0.4 to 30.0), rate reduction 0.2 for CIN3, compared with 45.3% (95% CI 29.8 to 57.6), rate reduction 0.3, for the analyses restricted to disease related to vaccine types [3, 22]. For the bivalent vaccine, the VE for CIN3 at four years end of study follow-up was similar for TVC analyses related to HPV 16 and 18 endpoints 45.7% (95% CI 22.9 to 62.2), rate reduction 0.13, as compared with endpoints irrespective of HPV type VE 45.6% (95% CI 28.8 to 58.7), rate reduction 0.22 [20].

Studies have shown that a greater proportion of CIN3+ lesions in younger women (under 30 years of age) is due to HPV type 16, with other less oncogenic types more prevalent in women over 30 years of age [57]. HPV16 is likely to overwhelmingly predominate as a cause of high-grade disease in women aged under 20, which is the age group at screening of the youngest women in our study (age 14 years in 2007). It is also possible that the intensive surveillance of the trial populations detected more non-HPV16/18 related CIN than occurs in real world screening programs, where many of these abnormalities will develop and resolve within the screening interval and are thus never observed. This could explain why a larger percentage of disease appears to be prevented in the population setting than in the trial setting. The predominance of HPV16 related disease in younger women has implications for the post-vaccination surveillance of high-grade lesions, as a greater impact of the vaccine may be expected initially. As vaccinated cohorts age and span a broader age range that includes women in their 30s, where other less oncogenic HPV types are more prevalent and may take longer to progress, the impact of vaccination on the overall CIN3+ burden may be relatively reduced [2, 12].

We observed a significant impact of vaccination on LG cytological abnormalities that was similar in magnitude to that observed for HG cytological lesions. In addition, the effect of partial vaccination was significant for LG cytology, but not for HG cytology, although this may be due to the larger number of outcomes in this group. Because a greater proportion of HG cytology is due to vaccine-related HPV types than LG cytology, we would have expected to see a greater impact on HG as compared with LG cytological lesions. However, this result was sensitive to whether we counted person-time in between dose one to three for those women who ended up being fully vaccinated (see Additional file 1: Table S1). Abnormalities occurring during this period are likely to be a result of pre-existing HPV infection due to the short time frame between doses, but we believe that this reflects the real world vaccine effectiveness in the population. Two meta-analyses showed that the contribution of HPV16 and 18 worldwide in LSIL lesions is approximately 34% [8, 23] while in CIN3, it is approximately 57% [23], with some regional variation. Recent Australian data show a relatively high proportion, approximately 41%, of LG cytological lesions contain HPV types 16/18/6/11 in women aged 25 years and under (personal communication JMLB) and this may also in part explain our observed reduction in LG cytological abnormalities.

There are several potential limitations of our analysis. We did not have a complete date of migration out of Victoria and there may be a small degree of loss to follow-up on the VCCR. The data linkage for this analysis was conducted according to accepted protocols; although some inaccuracies are an inherent part of the linkage methodology, the effect of this is likely to be small [24]. Any under-reporting of doses to the NHVPR given in the community to women who were within the school cohort, or doses given to women who opt-off the register, would have biased results towards the null if a proportion of unvaccinated women were truly vaccinated. Follow-up time was limited due to the recency of screening commencement among women vaccinated at school.

As this is an observational study, rather than a randomized study, women who chose to be vaccinated are likely to be different from those who are not vaccinated. Women of lower socioeconomic status were somewhat less likely to be vaccinated in our study, although results were unchanged in the adjusted analysis. Uncontrolled confounding by other factors may also have influenced the magnitude of the observed impact of the vaccine. However, as only women who had ever had a Pap test are included in this analysis, the effect of any sociodemographic factors that cause both lower vaccination coverage and screening uptake (for example, being a member of a culturally and linguistically diverse group) is minimized.

If fewer than three doses could provide sufficient protection against disease, vaccine programs could achieve higher coverage at a lower cost [2527]. There are several possible explanations as to why we did not observe a protective effect of one or two doses. Notably, women who only received one or two doses of vaccine rather than the full course were younger at first screening (suggesting earlier sexual debut), older at vaccination and of lower socioeconomic status. Adjusting for age at first screening and SES may mitigate somewhat but cannot completely control for the underlying differences in this group of women, who appear to be more likely to have been infected prior to vaccination than the fully vaccinated group. Analysis including more women, vaccinated prior to sexual debut and taking into account underlying sociodemographic differences among women who do not complete the course, is needed to refine the estimates of impact against disease for partially vaccinated women.

Our results contribute to the growing body of evidence in Australia of a reduced burden of HPV-related disease since implementation of the vaccination program in 2007. Data from HPV prevalence studies show a significant reduction in the prevalence of vaccine HPV genotypes (6, 11, 16, 18) in a post-vaccination sample compared with a pre-vaccination sample and lower prevalence in both vaccinated and unvaccinated women compared with the pre-vaccine reference population [28]. A substantial impact of the HPV vaccine on genital warts has already been observed in Australia, with a decline by the end of 2009 of 59% in warts diagnosed at sexual health services among women in the target age group (26 and under in 2007), and a smaller decline of 28% among heterosexual men of similar age [29]. A similar decline in the incidence of genital warts has been reported in New Zealand, the United States and Sweden [3033]. Ecologic data [34] and routine monitoring data [10, 12] show significant declines in histologically confirmed high-grade cervical abnormalities since the Australian vaccination program commenced among women under 20 years of age. In 2011 we demonstrated an absolute decrease in incidence of high-grade cervical abnormalities in Victoria of 0.38% (CI 95% 0.61 to 0.16) for women under 18 years of age [34]. The main limitation of these analyses was that individual level data on vaccination status were not available and thus it was not possible to attribute the decline in disease solely to the vaccination program.


In conclusion, we have demonstrated the first population level impact of a school-based HPV vaccination program on cervical abnormalities. The program included a two-year catch-up for girls aged 14- to 17-years and achieved high vaccination coverage rates. The impact of the vaccine was greatest for women completely vaccinated at younger ages and for the more serious cervical abnormalities CIN3/AIS, which are more likely to be due to the vaccine HPV types. The reduction in risk of cervical abnormalities was restricted to women who received three doses of vaccine; however, relatively few women received fewer than three doses and these women were older at vaccination than the other women. Over the next five years, the effect of the vaccine is expected to increase as cohorts of women vaccinated prior to sexual debut move through the screening program.





Human papillomavirus


Hazard ratio


High-grade squamous intraepithelial lesion




Low-grade squamous intraepithelial lesions


National HPV Vaccination Program Register


Socio-Economic Indexes for Areas


Socioeconomic status


Total Vaccine Cohort


Victorian Cervical Cytology Registry


Vaccine effectiveness.


  1. Li N, Franceschi S, Howell-Jones R, Snijders PJ, Clifford GM: Human papillomavirus type distribution in 30,848 invasive cervical cancers worldwide: variation by geographical region, histological type and year of publication. Int J Cancer. 2010, 128: 927-935.

    Article  Google Scholar 

  2. Smith JS, Lindsay L, Hoots B, Keys J, Franceschi S, Winer R, Clifford GM: Human papillomavirus type distribution in invasive cervical cancer and high-grade cervical lesions: a meta-analysis update. Int J Cancer. 2007, 121: 621-632. 10.1002/ijc.22527.

    Article  CAS  PubMed  Google Scholar 

  3. The Future II Study Group: Quadrivalent vaccine against human papillomavirus to prevent high-grade cervical lesions. New Engl J Med. 2007, 356: 1915-1927.

    Article  Google Scholar 

  4. Paavonen J, Naud P, Salmeron J, Wheeler CM, Chow SN, Apter D, Kitchener H, Castellsague X, Teixeira JC, Skinner SR, Hedrick J, Jaisamrarn U, Limson G, Garland S, Szarewski A, Romanowski B, Aoki FY, Schwarz TF, Poppe WA, Bosch FX, Jenkins D, Hardt K, Zahaf T, Descamps D, Struyf F, Lehtinen M, Dubin G, HPV PATRICIA Study Group: Efficacy of human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine against cervical infection and precancer caused by oncogenic HPV types (PATRICIA): final analysis of a double-blind, randomised study in young women. Lancet. 2009, 374: 301-314. 10.1016/S0140-6736(09)61248-4.

    Article  CAS  PubMed  Google Scholar 

  5. Baandrup L, Munk C, Andersen KK, Junge J, Iftner T, Kjaer SK: HPV16 is associated with younger age in women with cervical intraepithelial neoplasia grade 2 and 3. Gynecol Oncol. 2012, 124: 281-285. 10.1016/j.ygyno.2011.10.020.

    Article  CAS  PubMed  Google Scholar 

  6. Castle PE, Schiffman M, Wheeler CM, Wentzensen N, Gravitt PE: Human papillomavirus genotypes in cervical intraepithelial neoplasia grade 3. Cancer Epidemiol Biomarkers Prev. 2010, 19: 1675-1681. 10.1158/1055-9965.EPI-10-0251.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Brotherton JM, Tabrizi SN, Garland SM: Does HPV type 16 or 18 prevalence in cervical intraepithelial neoplasia grade 3 lesions vary by age? An important issue for postvaccination surveillance. Future Microbiol. 2012, 7: 193-199. 10.2217/fmb.11.161.

    Article  PubMed  Google Scholar 

  8. Clifford GM, Rana RK, Franceschi S, Smith JS, Gough G, Pimenta JM: Human papillomavirus genotype distribution in low-grade cervical lesions: comparison by geographic region and with cervical cancer. Cancer Epidemiol Biomarkers Prev. 2005, 14: 1157-1164. 10.1158/1055-9965.EPI-04-0812.

    Article  PubMed  Google Scholar 

  9. Markowitz LE, Tsu V, Deeks SL, Cubie H, Wang SA, Vicari AS, Brotherton JM: Human papillomavirus vaccine introduction - the first five years. Vaccine. 2012, 30: F139-F148.

    Article  PubMed  Google Scholar 

  10. Victorian Cervical Cytology Registry: Statistical Report 2011. VCCR Update on Cervical Screening in Victoria. [http://www/]

  11. Brotherton J, Gertig D, Chappell G, Rowlands L, Saville M: Catching up with the catch-up: HPV vaccination coverage data for Australian women aged 18–26 year from the National HPV Vaccination Program Register. Commun Dis Intell Q Rep. 2011, 35: 197-201.

    PubMed  Google Scholar 

  12. Australian Institute of Health and Welfare: Cancer Series Number 67. Cervical Screening in Australia. 2009, Canberra: Australian Institute of Health and Welfare: Australian Institute of Health and Welfare Cat. no.: CAN63, [], -2010

    Google Scholar 

  13. Gertig DM, Brotherton JM, Saville M: Measuring human papillomavirus (HPV) vaccination coverage and the role of the National HPV Vaccination Program Register, Australia. Sex Health. 2011, 8: 171-178. 10.1071/SH10001.

    Article  PubMed  Google Scholar 

  14. Brotherton J: The first cab off the rank: the early Australian experience of universally funded HPV vaccine for women aged 12–26 years. HPV Today. 2009, Number 17. []

    Google Scholar 

  15. Chief Medical Officer: Guidance on revaccination where HPV vaccine doses have been given at less than recommended minimum intervals. 2009, [$File/CMO-full-advice-hpv.pdf]

    Google Scholar 

  16. The NHMRC Screening to Prevent Cervical Cancer: Guidelines for the Management of Asymptomatic Women with Screen Detected Abnormalities. []

  17. ABS: Information Paper: an introduction to Socio-Economic Indexes for Areas (SEIFA) 2006. 2008, Canberra: ABS,

    Google Scholar 

  18. SAS® Enterprise Guide®. []

  19. Schiller JT, Castellsague X, Garland SM: A review of clinical trials of human papillomavirus prophylactic vaccines. Vaccine. 2012, 30: F123-F138.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Lehtinen M, Paavonen J, Wheeler CM, Jaisamrarn U, Garland SM, Castellsagué X, Skinner SR, Apter D, Naud P, Salmerón J, Chow SN, Kitchener H, Teixeira JC, Hedrick J, Limson G, Szarewski A, Romanowski B, Aoki FY, Schwarz TF, Poppe WA, De Carvalho NS, Germar MJ, Peters K, Mindel A, De Sutter P, Bosch FX, David MP, Descamps D, Struyf F, Dubin G, HPV PATRICIA Study Group: Overall efficacy of HPV-16/18 AS04-adjuvanted vaccine against grade 3 or greater cervical intraepithelial neoplasia: 4-year end-of-study analysis of the randomised, double-blind PATRICIA trial. Lancet Oncol. 2012, 13: 89-99. 10.1016/S1470-2045(11)70286-8. Erratum in: Lancet Oncol 2012, 13:e1

    Article  CAS  PubMed  Google Scholar 

  21. Paavonen J, Jenkins D, Bosch FX, Naud P, Salmerón J, Wheeler CM, Chow SN, Apter DL, Kitchener HC, Castellsague X, de Carvalho NS, Skinner SR, Harper DM, Hedrick JA, Jaisamrarn U, Limson GA, Dionne M, Quint W, Spiessens B, Peeters P, Struyf F, Wieting SL, Lehtinen MO, Dubin G, HPV PATRICIA Study Group: Efficacy of a prophylactic adjuvanted bivalent L1 virus-like-particle vaccine against infection with human papillomavirus types 16 and 18 in young women: an interim analysis of a phase III double-blind, randomised controlled trial. Lancet. 2007, 369: 2161-2170. 10.1016/S0140-6736(07)60946-5. Erratum in: Lancet 2007, 370:1414

    Article  CAS  PubMed  Google Scholar 

  22. Muñoz N, Kjaer SK, Sigurdsson K, Iversen OE, Hernandez-Avila M, Wheeler CM, Perez G, Brown DR, Koutsky LA, Tay EH, Garcia PJ, Ault KA, Garland SM, Leodolter S, Olsson SE, Tang GW, Ferris DG, Paavonen J, Steben M, Bosch FX, Dillner J, Huh WK, Joura EA, Kurman RJ, Majewski S, Myers ER, Villa LL, Taddeo FJ, Roberts C, Tadesse A, et al: Impact of human papillomavirus (HPV)-6/11/16/18 vaccine on all HPV-associated genital diseases in young women. J Natl Cancer Inst. 2010, 102: 325-339. 10.1093/jnci/djp534.

    Article  PubMed  Google Scholar 

  23. Guan P, Howell-Jones R, Li N, Bruni L, de Sanjose S, Franceschi S, Clifford GM: Human papillomavirus types in 115,789 HPV-positive women: a meta-analysis from cervical infection to cancer. Int J Cancer. 2012, 131: 2349-2359. 10.1002/ijc.27485.

    Article  CAS  PubMed  Google Scholar 

  24. Bohensky MA, Jolley D, Sundararajan V, Evans S, Ibrahim J, Brand C: Development and validation of reporting guidelines for studies involving data linkage. Aust NZ J Public Health. 2011, 35: 486-489. 10.1111/j.1753-6405.2011.00741.x.

    Article  Google Scholar 

  25. Kreimer AR, Gonzalez P, Katki HA, Porras C, Schiffman M, Rodriguez AC, Solomon D, Jimenez S, Schiller JT, Lowy DR, van Doorn LJ, Struijk L, Quint W, Chen S, Wacholder S, Hildesheim A, Herrero R, CVT Vaccine Group: Efficacy of a bivalent HPV 16/18 vaccine against anal HPV 16/18 infection among young women: a nested analysis within the Costa Rica Vaccine Trial. Lancet Oncol. 2011, 12: 862-870. 10.1016/S1470-2045(11)70213-3. Erratum in: Lancet Oncol 2011, 12:1096

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Dobson SR, McNeil S, Dionne M, Dawar M, Ogilvie G, Krajden M, Sauvageau C, Scheifele DW, Kollmann TR, Halperin SA, Langley JM, Bettinger JA, Singer J, Money D, Miller D, Naus M, Marra F, Young E: Immunogenicity of 2 doses of HPV vaccine in younger adolescents vs 3 doses in young women: a randomized clinical trial. JAMA. 2013, 309: 1793-1802. 10.1001/jama.2013.1625.

    Article  CAS  PubMed  Google Scholar 

  27. Romanowski B, de Borba PC, Naud PS, Roteli-Martins CM, De Carvalho NS, Teixeira JC, Aoki F, Ramjattan B, Shier RM, Somani R, Barbier S, Blatter MM, Chambers C, Ferris D, Gall SA, Guerra FA, Harper DM, Hedrick JA, Henry DC, Korn AP, Kroll R, Moscicki AB, Rosenfeld WD, Sullivan BJ, Thoming CS, Tyring SK, Wheeler CM, Dubin G, Schuind A, Zahaf T, GlaxoSmithKline Vaccine HPV-007 Study Group, et al: Sustained efficacy and immunogenicity of the human papillomavirus (HPV)-16/18 AS04-adjuvanted vaccine: analysis of a randomised placebo-controlled trial up to 6.4 years. Lancet. 2009, 374: 1975-1985.

    Article  CAS  PubMed  Google Scholar 

  28. Tabrizi SN, Brotherton JM, Kaldor JM, Skinner SR, Cummins E, Liu B, Bateson D, McNamee K, Garefalakis M, Garland SM: Fall in human papillomavirus prevalence following a national vaccination program. J Infect Dis. 2012, 206: 1645-1651. 10.1093/infdis/jis590.

    Article  PubMed  Google Scholar 

  29. Ali H, Donovan B, Wand H, Read TR, Regan DG, Grulich AE, Fairley CK, Guy RJ: Genital warts in young Australians five years into national human papillomavirus vaccination programme: national surveillance data. BMJ. 2013, 346: f2032-10.1136/bmj.f2032. Erratum in: BMJ 2013, 346:F2942

    Article  PubMed  Google Scholar 

  30. Bauer HM, Wright G, Chow J: Evidence of human papillomavirus vaccine effectiveness in reducing genital warts: an analysis of California public family planning administrative claims data, 2007-2010. Am J Public Health. 2012, 102: 833-835. 10.2105/AJPH.2011.300465.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Oliphant J, Perkins N: Impact of the human papillomavirus (HPV) vaccine on genital wart diagnoses at Auckland Sexual Health Services. N Z Med J. 2011, 124: 51-58.

    PubMed  Google Scholar 

  32. Leval A, Herweijer E, Arnheim-Dahlstrom L, Walum H, Frans E, Sparen P, Simard JF: Incidence of genital warts in Sweden before and after quadrivalent human papillomavirus vaccine availability. J Infect Dis. 2012, 206: 860-866. 10.1093/infdis/jis405.

    Article  PubMed  Google Scholar 

  33. Sexually transmitted infections in New Zealand. 2011, [], Surveillance Report

  34. Brotherton JM, Fridman M, May CL, Chappell G, Saville AM, Gertig DM: Early effect of the HPV vaccination programme on cervical abnormalities in Victoria, Australia: an ecological study. Lancet. 2011, 377: 2085-2092. 10.1016/S0140-6736(11)60551-5.

    Article  PubMed  Google Scholar 

Pre-publication history

Download references


We would like to thank Cathryn May, Karen Winch and the information technology team at VCS for assistance in data extraction and testing, as well as the Data Linkage Unit from the Australian Institute of Health and Welfare for conducting the data linkage. We would also like to thank Michael Malloy at VCS for statistical support. The VCCR is fully funded by the Victorian Department of Health and operated by the Victorian Cytology Service. The NHVPR is owned by the Department of Health and Ageing and operated by the Victorian Cytology Service.

Author information

Authors and Affiliations


Corresponding author

Correspondence to A Marion Saville.

Additional information

Competing interests

DMG, JMLB and MS were investigators on an Australian Research Council Linkage grant, for which CSL biotherapies was a partner organization. JMLB has been an investigator on investigator designed unrestricted epidemiological research grants partially funded through bioCSL but received no personal financial benefits.

Authors’ contributions

MS, DMG and JMLB were principal investigators of the study, with assistance from KD and GC. JMLB, DMG and MS designed the study. AB conducted the statistical analysis and interpretation of results was led by DMG, JMLB and MS. DMG led the writing of the manuscript, together with AB, JMLB and MS, and all authors contributed to the final report. All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1: Rate of histological and cytological cervical abnormalities for completely vaccinated, partially vaccinated and unvaccinated women using final dose status method (secondary analysis using ‘efficacy’ method). (DOCX 14 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Gertig, D.M., Brotherton, J.M., Budd, A.C. et al. Impact of a population-based HPV vaccination program on cervical abnormalities: a data linkage study. BMC Med 11, 227 (2013).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: