Adverse Events After Immunization With HPV Vaccine
Adverse Events After Immunization With HPV Vaccine
This register based cohort study of serious adverse events associated with the qHPV vaccine was based on individual level data from all 10 to 17 year old adolescent girls in Denmark and Sweden, between 1 October 2006 and 31 December 2010. Every resident in both Denmark and Sweden has a unique personal identification number enabling individual level linkage between multiple registers. To define the study cohort in Sweden, we obtained information on birth date and date of death from the total population register and death register, respectively, from Statistics Sweden. To define the study cohort in Denmark, we used the Civil Registration System, which contains daily updated information on vital and demographic variables, such as birth date and place, and loss to follow-up due to emigration or disappearance from the registers, and death.
The qHPV vaccine (Gardasil; Sanofi Pasteur MSD SNC, Lyon, France (in the United States: Merck, Whitehouse Station, NJ)) was marketed in Europe on 20 September 2006. In Sweden, the qHPV vaccine has been subsidised for 13-17 year old adolescent girls since May 2007 (inclusion of HPV vaccination in the national vaccination programme for 10-12 year old girls was initially planned for January 2010 but was postponed to January 2012 and then coupled with catch-up vaccination of 13-17 year old adolescent girls). In Denmark, the qHPV vaccine has been included in the national vaccination programme since January 2009 for 12 year old girls, with catch-up vaccination of 13-15 year olds from October 2008.
We obtained information on exposure to HPV vaccine in Denmark from the childhood vaccination database at Statens Serum Institut. This database is continually updated from National Health Insurance data obtained from the National Board of Health. General practitioners in Denmark carry out all immunisations in the childhood vaccination programme and are reimbursed for reporting each instance to the National Health Insurance. Reimbursement takes place only after the general practitioner has submitted a form that details the vaccine, the date it was administered, and the personal identification number of the recipient. Therefore the database is thought to be close to complete for vaccines administered through the national programme. Because HPV vaccination was also available outside the national programme in the study period, we retrieved additional vaccination data from the national prescription register, which contains individual level information on all prescriptions filled at all Danish pharmacies. This includes the personal identification number of the recipient, the date the prescription was dispensed, and the Anatomic Therapeutic Chemical (ATC) code. The ATC code used to identify the qHPV vaccine was J07BM01. In Sweden, correspondingly, we obtained information on vaccination status from Svevac and the drug prescription register. Svevac is a national HPV vaccination register, established in 2006 and held by the Swedish Institute for Communicable Disease. Healthcare staff who administer the vaccines report to Svevac on a voluntary basis, and the register has an estimated completeness of about 80%. The drug prescription register contains all prescription drugs dispensed at pharmacies in Sweden since 1 July 2005. Adolescent girls aged between 13 and 17 years received subsidised HPV vaccination and vaccination had to be prescribed and expedited at a pharmacy, thereby generating a register entry in the drug prescription register. The register is held by the National Board of Health and Welfare. For adolescent girls aged 13-17 years it is assumed that almost 100% of administered HPV vaccine doses are registered in the drug prescription register.
We identified data on predefined adverse events from the national patient registers in both countries using ICD-10 codes (international classification of diseases, 10th revision). The patient registers include nationwide individual level information on dates of hospital contact and doctor assigned diagnoses according to the international classification of diseases. We did not have information on outcomes from primary healthcare. The Danish patient register was established in 1977, has included both inpatients and outpatients since 1994, and has been using ICD-10 codes since 1995, whereas the Swedish patient register has had nationwide coverage since 1987, has included both inpatients and outpatients since 2001, and has used ICD-10 codes since 1997. We predefined several serious adverse outcome events, as identified from records of inpatient admissions and hospital outpatient and emergency department visits, based on our earlier study of autoimmune events, and we added several neurological events. We also included venous thromboembolism because it represents a potential adverse event. The included outcomes are all well defined diseases. In total, we assessed 53 outcomes (see Supplementary Table 1 for all included outcome events, with ICD-10 codes).
From Statistics Denmark and Statistics Sweden we obtained data on parental educational level, country of birth, and socioeconomic status. We identified the parents from the Danish Civil Registration System and Swedish multigeneration register.
Adolescent girls were followed from age 10 years or 1 October 2006, whichever came latest, and until either the occurrence of an adverse event, receipt of bivalent HPV vaccine (Cervarix, GlaxoSmithKline Biologicals; Rixensart, Belgium; ATC code J07BM02), death, disappearance from the registers, emigration (this information was only available for Denmark), 18th birthday, or end of follow-up (31 December 2010). We aggregated the resulting person years of follow-up with counts of outcome events according to qHPV vaccine exposure status and analysed these using Poisson regression (log-linear regression of the counts using the logarithm of follow-up time as offset). This produced incidence rate ratios according to qHPV exposure status. Exposure to the qHPV vaccine was a time varying variable; thus adolescent girls could contribute person time to the study first as unvaccinated and later as vaccinated, but once vaccinated the girls could not be put into the unvaccinated category again. All individual outcomes were treated as separate analyses, and for each specific analysis girls were eligible only if free from the outcome event before entry to the cohort. Estimates were adjusted for country, age in two year categories, calendar year, parental educational level (highest attained level of either parent classified as: primary school (nine years) or shorter; secondary school (12 years); short tertiary education; and medium or long tertiary education), parental country of birth (categories: both parents, one parent, or no parent born in Scandinavia), and paternal socioeconomic status (categories: employment with basic, unknown, or no qualification; employment with medium level or high level qualifications; self employed; and not in labour market).
For all autoimmune and neurological outcomes, we defined the period at risk as 180 days after exposure to vaccine. This period was chosen to allow for the insidious onset of the diseases studied and because diagnostic investigations may take time; in a recent study of autoimmune outcomes after qHPV vaccination, the median time between first symptoms and diagnosis was 23 days (interquartile range 2-59 days). For venous thromboembolism, the period at risk was 90 days after vaccination. This was regarded as the maximal period where an acute event could be plausibly related to vaccination; furthermore, the mean time between vaccination and diagnosis of thromboembolism among cases reported to the Vaccine Adverse Event Reporting System was 42 days.
Because we acquired data on vaccine exposure from two sources that were partly overlapping—that is, some of the girls had both filled prescriptions and were registered in one of the vaccination databases—we applied an algorithm to harmonise the data and define dates of vaccination. Essentially this algorithm removed double data entries and doses appearing beyond the three dose schedule. Furthermore, on the basis of data from girls with both filled prescriptions and register entries in vaccination databases, it was established that the median lag between the dispensing of a prescription and the date of vaccination as registered in the vaccination databases was two days. Consequently, for vaccine doses that were defined by prescriptions alone, the date of exposure was defined as the date of filling the prescription plus two days.
As the recommended qHPV vaccine schedule includes three doses given at 0, 2, and 6 months, any girl could contribute up to three doses in the analyses; we counted exposed person time from the date the vaccine was administered, and each dose contributed up to 180 days (autoimmune and neurological events) of follow-up or up to 90 days (venous thromboembolism) of follow-up (Figure 1). We used SAS statistical software version 9.3 (SAS Institute, Cary, NC).
(Enlarge Image)
Figure 1.
Periods at risk for autoimmune and neurological events in adolescent girls after exposure to quadrivalent human papillomavirus (qHPV) vaccine. For venous thromboembolism, each period at risk was up to 90 days.
Given that a large number of serious adverse events were assessed and consequently there was a high probability of chance findings, we used the following predefined criteria for the analysis of data. As the first criterion, and in the interest of obtaining relatively reliable rate ratios, for any further assessment to take place we considered only outcomes with at least five vaccine exposed cases during the predefined period at risk. To be regarded as a safety signal, the rate ratio for an outcome with at least five vaccine exposed cases had to be significantly increased (lower bound of 95% confidence interval >1.0). We regarded three criteria as signal strengthening: analysis based on 20 or more vaccine exposed cases (reliability of analysis); a rate ratio of 3.0 or more (strength of association); and significantly increased rate ratios in both countries when analysed separately (consistency). For outcomes where the rate ratio was significantly increased, we additionally assessed clustering of events in time by plotting cases according to time since exposure to the vaccine and estimated rate ratios for a risk period that started on day 181.
Methods
Study Population
This register based cohort study of serious adverse events associated with the qHPV vaccine was based on individual level data from all 10 to 17 year old adolescent girls in Denmark and Sweden, between 1 October 2006 and 31 December 2010. Every resident in both Denmark and Sweden has a unique personal identification number enabling individual level linkage between multiple registers. To define the study cohort in Sweden, we obtained information on birth date and date of death from the total population register and death register, respectively, from Statistics Sweden. To define the study cohort in Denmark, we used the Civil Registration System, which contains daily updated information on vital and demographic variables, such as birth date and place, and loss to follow-up due to emigration or disappearance from the registers, and death.
Vaccination
The qHPV vaccine (Gardasil; Sanofi Pasteur MSD SNC, Lyon, France (in the United States: Merck, Whitehouse Station, NJ)) was marketed in Europe on 20 September 2006. In Sweden, the qHPV vaccine has been subsidised for 13-17 year old adolescent girls since May 2007 (inclusion of HPV vaccination in the national vaccination programme for 10-12 year old girls was initially planned for January 2010 but was postponed to January 2012 and then coupled with catch-up vaccination of 13-17 year old adolescent girls). In Denmark, the qHPV vaccine has been included in the national vaccination programme since January 2009 for 12 year old girls, with catch-up vaccination of 13-15 year olds from October 2008.
We obtained information on exposure to HPV vaccine in Denmark from the childhood vaccination database at Statens Serum Institut. This database is continually updated from National Health Insurance data obtained from the National Board of Health. General practitioners in Denmark carry out all immunisations in the childhood vaccination programme and are reimbursed for reporting each instance to the National Health Insurance. Reimbursement takes place only after the general practitioner has submitted a form that details the vaccine, the date it was administered, and the personal identification number of the recipient. Therefore the database is thought to be close to complete for vaccines administered through the national programme. Because HPV vaccination was also available outside the national programme in the study period, we retrieved additional vaccination data from the national prescription register, which contains individual level information on all prescriptions filled at all Danish pharmacies. This includes the personal identification number of the recipient, the date the prescription was dispensed, and the Anatomic Therapeutic Chemical (ATC) code. The ATC code used to identify the qHPV vaccine was J07BM01. In Sweden, correspondingly, we obtained information on vaccination status from Svevac and the drug prescription register. Svevac is a national HPV vaccination register, established in 2006 and held by the Swedish Institute for Communicable Disease. Healthcare staff who administer the vaccines report to Svevac on a voluntary basis, and the register has an estimated completeness of about 80%. The drug prescription register contains all prescription drugs dispensed at pharmacies in Sweden since 1 July 2005. Adolescent girls aged between 13 and 17 years received subsidised HPV vaccination and vaccination had to be prescribed and expedited at a pharmacy, thereby generating a register entry in the drug prescription register. The register is held by the National Board of Health and Welfare. For adolescent girls aged 13-17 years it is assumed that almost 100% of administered HPV vaccine doses are registered in the drug prescription register.
Outcomes
We identified data on predefined adverse events from the national patient registers in both countries using ICD-10 codes (international classification of diseases, 10th revision). The patient registers include nationwide individual level information on dates of hospital contact and doctor assigned diagnoses according to the international classification of diseases. We did not have information on outcomes from primary healthcare. The Danish patient register was established in 1977, has included both inpatients and outpatients since 1994, and has been using ICD-10 codes since 1995, whereas the Swedish patient register has had nationwide coverage since 1987, has included both inpatients and outpatients since 2001, and has used ICD-10 codes since 1997. We predefined several serious adverse outcome events, as identified from records of inpatient admissions and hospital outpatient and emergency department visits, based on our earlier study of autoimmune events, and we added several neurological events. We also included venous thromboembolism because it represents a potential adverse event. The included outcomes are all well defined diseases. In total, we assessed 53 outcomes (see Supplementary Table 1 for all included outcome events, with ICD-10 codes).
Covariate Information
From Statistics Denmark and Statistics Sweden we obtained data on parental educational level, country of birth, and socioeconomic status. We identified the parents from the Danish Civil Registration System and Swedish multigeneration register.
Statistical Analyses
Adolescent girls were followed from age 10 years or 1 October 2006, whichever came latest, and until either the occurrence of an adverse event, receipt of bivalent HPV vaccine (Cervarix, GlaxoSmithKline Biologicals; Rixensart, Belgium; ATC code J07BM02), death, disappearance from the registers, emigration (this information was only available for Denmark), 18th birthday, or end of follow-up (31 December 2010). We aggregated the resulting person years of follow-up with counts of outcome events according to qHPV vaccine exposure status and analysed these using Poisson regression (log-linear regression of the counts using the logarithm of follow-up time as offset). This produced incidence rate ratios according to qHPV exposure status. Exposure to the qHPV vaccine was a time varying variable; thus adolescent girls could contribute person time to the study first as unvaccinated and later as vaccinated, but once vaccinated the girls could not be put into the unvaccinated category again. All individual outcomes were treated as separate analyses, and for each specific analysis girls were eligible only if free from the outcome event before entry to the cohort. Estimates were adjusted for country, age in two year categories, calendar year, parental educational level (highest attained level of either parent classified as: primary school (nine years) or shorter; secondary school (12 years); short tertiary education; and medium or long tertiary education), parental country of birth (categories: both parents, one parent, or no parent born in Scandinavia), and paternal socioeconomic status (categories: employment with basic, unknown, or no qualification; employment with medium level or high level qualifications; self employed; and not in labour market).
For all autoimmune and neurological outcomes, we defined the period at risk as 180 days after exposure to vaccine. This period was chosen to allow for the insidious onset of the diseases studied and because diagnostic investigations may take time; in a recent study of autoimmune outcomes after qHPV vaccination, the median time between first symptoms and diagnosis was 23 days (interquartile range 2-59 days). For venous thromboembolism, the period at risk was 90 days after vaccination. This was regarded as the maximal period where an acute event could be plausibly related to vaccination; furthermore, the mean time between vaccination and diagnosis of thromboembolism among cases reported to the Vaccine Adverse Event Reporting System was 42 days.
Because we acquired data on vaccine exposure from two sources that were partly overlapping—that is, some of the girls had both filled prescriptions and were registered in one of the vaccination databases—we applied an algorithm to harmonise the data and define dates of vaccination. Essentially this algorithm removed double data entries and doses appearing beyond the three dose schedule. Furthermore, on the basis of data from girls with both filled prescriptions and register entries in vaccination databases, it was established that the median lag between the dispensing of a prescription and the date of vaccination as registered in the vaccination databases was two days. Consequently, for vaccine doses that were defined by prescriptions alone, the date of exposure was defined as the date of filling the prescription plus two days.
As the recommended qHPV vaccine schedule includes three doses given at 0, 2, and 6 months, any girl could contribute up to three doses in the analyses; we counted exposed person time from the date the vaccine was administered, and each dose contributed up to 180 days (autoimmune and neurological events) of follow-up or up to 90 days (venous thromboembolism) of follow-up (Figure 1). We used SAS statistical software version 9.3 (SAS Institute, Cary, NC).
(Enlarge Image)
Figure 1.
Periods at risk for autoimmune and neurological events in adolescent girls after exposure to quadrivalent human papillomavirus (qHPV) vaccine. For venous thromboembolism, each period at risk was up to 90 days.
Analytical Strategy
Given that a large number of serious adverse events were assessed and consequently there was a high probability of chance findings, we used the following predefined criteria for the analysis of data. As the first criterion, and in the interest of obtaining relatively reliable rate ratios, for any further assessment to take place we considered only outcomes with at least five vaccine exposed cases during the predefined period at risk. To be regarded as a safety signal, the rate ratio for an outcome with at least five vaccine exposed cases had to be significantly increased (lower bound of 95% confidence interval >1.0). We regarded three criteria as signal strengthening: analysis based on 20 or more vaccine exposed cases (reliability of analysis); a rate ratio of 3.0 or more (strength of association); and significantly increased rate ratios in both countries when analysed separately (consistency). For outcomes where the rate ratio was significantly increased, we additionally assessed clustering of events in time by plotting cases according to time since exposure to the vaccine and estimated rate ratios for a risk period that started on day 181.
Source...