Intended for healthcare professionals

  1. Research
  2. Association between...
  3. Association between human papillomavirus vaccination and serious adverse events in South Korean adolescent girls: nationwide cohort study
CCBYNC Open access
Research

Association between human papillomavirus vaccination and serious adverse events in South Korean adolescent girls: nationwide cohort study

BMJ 2021; 372 doi: https://doi.org/10.1136/bmj.m4931 (Published 29 January 2021) Cite this as: BMJ 2021;372:m4931
  1. Dongwon Yoon, doctoral student1,
  2. Ji-Ho Lee, family medicine doctor12,
  3. Hyesung Lee, postdoctoral research fellow1,
  4. Ju-Young Shin, associate professor13
  1. 1School of Pharmacy, Sungkyunkwan University, 2066, Seobu-ro, Jangan-gu, Suwon, Gyeonggi-do, South Korea
  2. 2Chung-Ang University Hospital, 06973, 102 Heukseok-ro, Dongjak-gu, Seoul, South Korea
  3. 3Samsung Advanced Institute for Health Sciences and Technology (SAIHST), Sungkyunkwan University, Seoul, South Korea
  1. Correspondence to: J-Y Shin shin.jy{at}skku.edu
  • Accepted 2 December 2020

Abstract

Objective To evaluate the association between human papillomavirus (HPV) vaccination and serious adverse events in adolescent girls in South Korea.

Design Cohort study.

Setting A large linked database created by linking the Korea Immunization Registry Information System and the National Health Information Database, between January 2017 and December 2019.

Participants 441 399 girls aged 11-14 years who had been vaccinated in 2017: 382 020 had been vaccinated against HPV and 59 379 had not been vaccinated against HPV.

Main outcome measures Outcomes were 33 serious adverse events, including endocrine, gastrointestinal, cardiovascular, musculoskeletal, haematological, dermatological, and neurological diseases. A cohort design was used for the primary analysis and a self-controlled risk interval design for the secondary analysis; both analyses used a risk period of one year after HPV vaccination for each outcome. Incidence rate and adjusted rate ratios were estimated using Poisson regression in the primary analysis, comparing the HPV vaccinated group with the HPV unvaccinated group, and adjusted relative risks were estimated using conditional logistic regression in the secondary analysis.

Results Among the 33 predefined serious adverse events, no associations were found with HPV vaccination in the cohort analysis, including Hashimoto’s thyroiditis (incidence rate per 100 000 person years: 52.7 v 36.3 for the vaccinated and unvaccinated groups; adjusted rate ratio 1.24, 95% confidence interval 0.78 to 1.94) and rheumatoid arthritis (incidence rate per 100 000 person years: 168.1 v 145.4 for the vaccinated and unvaccinated groups; 0.99, 0.79 to 1.25), with the exception of an increased risk observed for migraine (incidence rate per 100 000 person years: 1235.0 v 920.9 for the vaccinated and unvaccinated groups; 1.11, 1.02 to 1.22). Secondary analysis using self-controlled risk intervals confirmed no associations between HPV vaccination and serious adverse events, including migraine (adjusted relative risk 0.67, 95% confidence interval 0.58 to 0.78). Results were robust to varying follow-up periods and for vaccine subtypes.

Conclusions In this nationwide cohort study, with more than 500 000 doses of HPV vaccines, no evidence was found to support an association between HPV vaccination and serious adverse events using both cohort analysis and self-controlled risk interval analysis. Inconsistent findings for migraine should be interpreted with caution considering its pathophysiology and the population of interest.

Introduction

In 2012, about 630 000 women received a diagnosis of human papillomavirus (HPV) related cancer, of whom 530 000 (84%) had cervical cancer.1 Around 266 000 women with cervical cancer died in 2012, accounting for 8% of all cancer related mortality in women worldwide.2 To reduce the HPV attributable burden of disease, vaccines against HPV have been introduced globally and are indicated for both adolescent boys and adolescent girls.3 The current World Health Organization recommendation for HPV immunisation prioritises girls aged 9-14 years before they become sexually active, and in 2019, 124 countries adopted the policy of HPV vaccination as their national immunisation programme for adolescent girls.45 Since June 2016, the Korea Disease Control and Prevention Agency has also introduced a two dose schedule of bivalent or quadrivalent HPV vaccines for adolescent girls aged 12-13 years as part of its national immunisation programme.67 Coverage of the first HPV vaccine dose in the first year of implementation was 61.5%, and by 2018 this had increased to 87.2%.8

Despite the extensive use of the HPV vaccine worldwide, safety concerns about HPV vaccination have not diminished, and concerns about vaccine related serious adverse events have undermined public confidence about immunisation.910 A survey among Korean parents who did not want their children to receive HPV vaccination indicated that 73.5% were reluctant because of concerns about serious adverse events.11 To alleviate these concerns, good quality evidence on the safety of the vaccine in the Korean population is needed. Because most studies have focused on Western populations,12131415161718 information on the safety profile of HPV vaccine in Asian populations is insufficient. Furthermore, previous studies compared participants who had been vaccinated with those who had not, which could lead to biased results because of fundamental differences in personal and clinical characteristics. Moreover, studies have focused on the bivalent or quadrivalent HPV vaccine and thus provide information on the safety of certain types of HPV vaccines.121319

Given the public concerns and limited information in the Korean population, we investigated the association between HPV vaccination and serious adverse events among adolescent girls in South Korea using a nationwide linked database.

Methods

Data sources

We combined data from the Korea Immunization Registry Information System and the National Health Information Database to create a large database between January 2017 and December 2019 of adolescent girls aged 11-14 years in 2017. Using deidentified personal codes, we retrieved data from the Korea Immunization Registry Information System on type of vaccine administered, date of vaccination, dosing schedule, coadministered vaccines, and age at vaccination, and data from the National Health Information Database on medical diagnosis and date of diagnosis, prescription, treatment, medical insurance type, residence type, and total income of the affiliated households. This large linked database enabled us to evaluate the association between HPV vaccination and serious adverse events using both immunisation and clinical information.

According to the system of universal insurance coverage in South Korea, the National Health Information Database represents the entire population (50 million).20 This database includes information on personal details and clinical records on reimbursed inpatient and outpatient visits. Diagnosis and drug prescription information was coded using ICD-10 (international classification of diseases, 10th revision) and Anatomical Therapeutic Chemical classification codes, respectively. The positive predictive value for diagnosis codes of claims data was reported to be 82% in a validation study.21

Study population

To control for the unmasking phenomenon and fundamental differences between vaccinated and unvaccinated individuals,22 we selected Japanese encephalitis vaccine and tetanus, diphtheria, and acellular pertussis vaccine as comparators in the cohort analysis as they have well established safety profiles and are recommended in the same population of interest.232425 Hence, the study population comprised girls aged 11-14 years who were vaccinated against HPV, Japanese encephalitis, and tetanus, diphtheria, and acellular pertussis in 2017. We excluded girls who were immunised with both bivalent and quadrivalent HPV vaccines in 2017, as well as those with a diagnosis of malignancy or congenital or neonatal disease (congenital iodine deficiency syndrome, congenital hypothyroidism, transitory neonatal hyperthyroidism, other specified diseases of the pancreas, transient neonatal thrombocytopenia, neonatal diabetes mellitus, cerebral palsy, and type 1 diabetes mellitus) before the first date of vaccination.

Vaccination

In South Korea, healthcare providers report detailed information on immunisation to the Korea Immunization Registry Information System and the government reimburses immunisation fees for all vaccines covered by the national immunisation programme. We obtained data on HPV, Japanese encephalitis, and tetanus, diphtheria, and acellular pertussis vaccinations from the Korea Immunization Registry Information System. A total of 443 290 adolescents girls received the study vaccines (HPV, Japanese encephalitis, and tetanus, diphtheria, and acellular pertussis) in 2017. We classified girls vaccinated with HPV as the HPV vaccinated group, and those vaccinated with comparator vaccines without HPV vaccination as the HPV unvaccinated group. Cohort entry was defined as the date of vaccination. Each girl could have two records of cohort entry because HPV vaccination comprises a two dosing schedule with an interval of six months between doses.

Outcomes

We selected 33 serious adverse events as the main outcomes, predefined based on literature review, including case reports, pharmacovigilance studies, observational studies, meta-analysis studies, and a review of the Korean medical advisory committee on vaccine safety.1314161826272829303132 Outcomes were classified into eight disease categories: endocrine (Graves’ disease, Hashimoto’s thyroiditis, hyperthyroidism, hypothyroidism, and type 1 diabetes), gastrointestinal (Crohn’s disease, ulcerative colitis, peptic ulcer, and pancreatitis), cardiovascular (Raynaud’s disease, venous thromboembolism, vasculitis, and hypotension), musculoskeletal and systemic (ankylosing spondylitis, Behcet’s syndrome, juvenile arthritis, rheumatoid arthritis, and systemic lupus erythematosus), haematological (idiopathic thrombocytopenic purpura, and Henoch-Schönlein’s purpura), dermatological (erythema nodosum and psoriasis), neurological (Bell’s palsy, epilepsy, narcolepsy, paralysis, migraine, Guillain-Barré syndrome, optical neuritis, neuralgia and neuritis, intracerebral haemorrhage, and extrapyramidal and movement disorders), and tuberculosis. Incident serious adverse events were identified using diagnosis records coded by ICD-10 (supplementary table S1).

Potential confounders

We included sociodemographic characteristics and anaemia as potential confounders in our analyses based on a causal diagram using a directed acyclic graph (supplementary figure S1). Because sociodemographic characteristics are well known confounders in vaccine safety studies,1213 we assessed age, region of residence, type of health insurance, and income level at cohort entry. In addition, because anaemia is related to HPV vaccination and could affect the health outcome of interest (supplementary section A), we assessed anaemia within one year before cohort entry.

Statistical analyses

To assess the association between HPV vaccination and the predefined serious adverse events separately we used both a cohort and a self-controlled risk interval analysis (fig 1). Sociodemographic, clinical, and vaccination related characteristics are summarised in our study population. In addition, we present the mean days (interquartile range) to disease diagnosis after each dose of HPV vaccine.

Fig 1
Fig 1

Overview of cohort and self-controlled risk interval study designs. HPV=human papillomavirus; SAE=serious adverse event; JEV=Japanese encephalitis vaccine; Tdap=tetanus, diphtheria, and acellular pertussis vaccine

Cohort analyses

We used a cohort design for our primary analysis because the results are more generalisable than those of the self-controlled risk interval analysis when vaccinated groups are compared with unvaccinated groups. Moreover, use of a nationwide linked database enables the association between HPV vaccine and serious adverse events to be identified without loss of statistical power. We defined the follow-up period as one year (365 days) to detect a potential insidious outcome and a delayed onset to diagnosis. To include incident cases only, we excluded girls who experienced a serious adverse event within one year before cohort entry. In the cohort analysis, we used a Poisson regression model to estimate the incidence rate ratio with 95% confidence interval. The incidence rate was calculated by dividing the number of serious adverse events by the total person years at risk and then multiplying the estimates by 100 000. The model was adjusted for age, region of residence, type of health insurance, income level, and anaemia.

Self-controlled risk interval analyses

We used the self-controlled risk interval design for secondary analysis by including girls with a history of both HPV vaccination and predefined serious adverse events. A major benefit of this within person comparison is that time invariant confounders (sex, region of residence, type of health insurance, income level, lifestyle) are self-adjusted, and the design is less subject to unmeasured confounders compared with other study designs. Using the self-controlled risk interval design, we tried to minimise the impact of bias from exposure misclassification and unmeasured time invariant confounders, which are limitations of the cohort design.33 The risk interval was defined as one year (365 days) after HPV vaccination to allow for latent diseases, and the control interval was defined as 466-730 days after HPV vaccination. We adopted a wash-out period of 90 days between the risk interval and control interval considering the possibility of difference between the disease onset and date of disease diagnosis.34 Only incident cases were included, and girls with serious adverse events in the preceding one year were excluded. We used a conditional logistic regression model to estimate the relative risk with 95% confidence interval for each serious adverse event. The model was adjusted for age, which is a time varying variable between the risk interval and the control interval.

Sensitivity and subgroup analyses

To confirm the robustness of the main analysis and to examine the possibility of an insidious outcome that was not identified even with a one year follow-up period, we varied the follow-up period in sensitivity analyses: short term (90 days and 180 days) and long term (two years (730 days)). We conducted subgroup analyses according to the subtypes of HPV (bivalent and quadrivalent) and then compared the HPV vaccinated group with the HPV unvaccinated groups to confirm the difference in risk among those factors.

All statistical analyses were performed using SAS Enterprise Guide 7.1 for Windows (SAS Institute, Cary, NC). A two tailed P value of less than 0.05 was considered significant.

Patient and public involvement

As we used deidentified data in our analyses, the need for informed consent was waived, and the study was conducted in accordance with the Guidelines for Good Pharmacoepidemiology Practices. No patients were involved in establishing the research question or outcome measures or in developing plans for recruitment, design, or implementation of the study. No patients were asked for advice on interpretation or documentation of the results. We have no plans to disseminate the results of the research to study participants or the relevant patient community.

Results

Cohort

Of the total 676 677 girls aged 11-14 years in South Korea, 441 399 (65%) were included in this study (fig 2). Among these girls, 382 020 were inoculated with 429 377 doses of the HPV vaccine and 59 379 with 87 099 doses of Japanese encephalitis vaccine or tetanus, diphtheria, and acellular pertussis vaccine without HPV vaccination (table 1 and fig 2). Among the HPV vaccinated girls, 295 365 received the quadrivalent vaccine and 86 655 received the bivalent vaccine (supplementary table S9). A total of 516 476 HPV vaccine doses were administered to the study population. The mean follow-up time for each outcome in the HPV vaccinated and HPV unvaccinated groups was 407 400 and 60 500 person years, respectively (table 2).

Fig 2
Fig 2

Flowchart of enrolment using large database created by linking data from the Korea Immunization Registry Information System with data from the National Health Information Database. HPV=human papillomavirus; Tdap=tetanus, diphtheria, and acellular pertussis vaccine; JEV=Japanese encephalitis vaccine

Table 1

Characteristics of girls aged 11-14 years vaccinated in South Korea in 2017. Values are numbers (percentages) unless stated otherwise

View this table:
Table 2

Cohort analysis of one year follow-up of girls aged 11-14 years vaccinated in South Korea in 2017

View this table:

Risk of serious adverse events after HPV vaccination

The mean number of days to diagnosis of the 33 predefined outcomes after HPV vaccination was generally consistent with each other (supplementary table S2). Among the 33 predefined outcomes, the overall estimates in the primary analyses showed no association between HPV vaccines and serious adverse events (fig 3). Rate ratios could not be calculated for ankylosing spondylitis, systemic lupus erythematosus, and erythema nodosum because none occurred in the HPV unvaccinated group. The risk for type 1 diabetes was not significantly increased (adjusted rate ratio 0.66, 95% confidence interval 0.38 to 1.14), and the result for the association between HPV vaccination and Hashimoto’s thyroiditis was not significant (1.24, 0.78 to 1.94). Other conditions not associated with HPV vaccination were cardiovascular diseases, including Raynaud’s disease (1.34, 0.40 to 4.49), and venous thromboembolism (0.62, 0.22 to 1.70). In musculoskeletal and systemic diseases, the rate ratio for the association between HPV vaccination and Behcet’s syndrome was 0.61 (0.22 to 1.67). No significant increase in risk was found for haematological and dermatological diseases. Although neurological diseases, including epilepsy, narcolepsy, paralysis, and Guillain-Barré syndrome showed no significant increase in risk, a significant increase was observed for migraine. However, the lower limit of the 95% confidence interval was close to 1 (1.11, 1.02 to 1.22). In line with the results of the cohort analysis, the results of the self-controlled risk interval analysis showed no significant increase in risk for all individual outcomes, including migraine (adjusted relative risk 0.67, 95% confidence interval 0.58 to 0.78).

Fig 3
Fig 3

Association between human papillomavirus vaccination and serious adverse events in both the cohort and the self-controlled risk interval (SCRI) analyses among girls aged 11-14 years and vaccinated in South Korea in 2017

Sensitivity and subgroup analyses

Results of the sensitivity analyses with diverse follow-up periods (90 days, 180 days, and 2 years) were consistent with those of the main findings, which showed no association between HPV vaccination and serious adverse events, except for migraine. Although a significant increase in risk was observed for migraine at the two year follow-up (adjusted rate ratio 1.10, 95% confidence interval 1.03 to 1.16), no significant increase was observed at 90 days (1.04, 0.85 to 1.26) and 180 days (1.09, 0.94 to 1.25) (supplementary section C).

Results from subgroup analyses for vaccine subtype were also consistent with those of the main findings. Although no significant increase was found for migraine in the bivalent HPV vaccinated group compared with HPV unvaccinated group (1.07, 0.96 to 1.20), the quadrivalent HPV vaccinated group showed a statistically significant increase in migraine compared with the HPV unvaccinated group (1.13, 1.03 to 1.24) (fig 4).

Fig 4
Fig 4

Association between human papillomavirus vaccination and serious adverse events in subgroup analyses by vaccine subtypes among girls aged 11-14 years and vaccinated in South Korea in 2017

Discussion

In this nationwide study of girls aged 11-14 years who received HPV, Japanese encephalitis, and tetanus, diphtheria, and acellular pertussis vaccines in South Korea in 2017, we used two complementary study designs to investigate the association between HPV vaccination and serious adverse events. All 33 predefined outcomes except for migraine showed no association between HPV vaccination and serious adverse events based on more than 400 000 girls and around 500 000 vaccine doses, and these findings were consistent across the follow-up periods (90 days, 180 days, and 2 years). Subgroup analysis by HPV vaccine subtype also showed consistent results with our main findings.

In cohort analyses, of a total of 33 serious adverse events, a significant increase in risk was observed only for migraine (adjusted rate ratio 1.11, 1.02 to 1.22), but the point estimate and lower limit of the 95% confidence interval were close to 1. Contrary to the results of the cohort analysis, the result of the self-controlled risk interval analysis presented no evidence of an association between HPV vaccination and migraine (relative risk 0.67, 95% confidence interval 0.58 to 0.78), which did not support the risk of migraine after HPV vaccination. This also means that migraine occurs more frequently after the risk interval. Possible explanations could be outcome misclassification, validity of migraine diagnosis, and a relatively low number of girls with migraine compared with those in the general population. Above all, we suspected that this inconsistency would probably be attributed to age. From a literature review and clinical consultation, we identified that hormonal changes and menstruation are well known risk factors for migraine.35 Considering the overall mean age at menarche (12.7 years) in the Korean population,36 our study population might encounter menstrual related migraine during the follow-up period. Moreover, given the mean age differences between the HPV vaccinated (12.42 years) and HPV unvaccinated (11.84 years) groups, menstrual related migraine is likely to occur more often in the HPV vaccinated group, resulting in statistically significantly increased results. Along with the results of the cohort analysis, our sensitivity analyses support an increased effect of age on migraine over time. Although sensitivity analyses for the short term period did not show a significant increase in the risk of migraine (adjusted rate ratios: 90 days, 1.04, 0.85 to 1.26; 180 days 1.09, 0.94 to 1.25), significant results were found for the long term period (2 years, 1.10, 1.03 to 1.16). Therefore, the effect of age would possibly contribute to this result, but this finding should be interpreted with caution given the pathophysiology of migraine.

Comparison with other studies

Our findings of no significant increase in serious adverse events show the safety of HPV vaccination and corroborate the findings of previous studies. Several studies have found no association between HPV vaccination and autoimmune diseases, a major concern in vaccine safety studies.13141937 Although a recent meta-analysis highlighted no clear association between HPV vaccines and autoimmune diseases, the study also suggested a significant increase in the risk of some types of autoimmune diseases: Hashimoto’s thyroiditis (odds ratio 1.25, 95% confidence interval 1.09 to 1.44) and Raynaud’s disease (1.63, 1.21 to 2.20).38 We speculated that this increased risk resulted from the unmasking phenomenon, which can occur when vaccinated groups are compared with unvaccinated groups.22 Visiting a medical institution for vaccination offers an opportunity to assess symptoms that might not have been assessed otherwise. The fundamental difference between vaccinated and unvaccinated groups would bias the results towards an increased risk after vaccination. Thus we used comparator vaccines to account for the unmasking phenomenon, and thereby Hashimoto’s thyroiditis showed no increased risk in both cohort and self-controlled risk interval analyses (adjusted rate ratio 1.24, 95% confidence interval 0.78 to 1.94; adjusted relative risk 1.28, 95% confidence interval 0.53 to 3.08, respectively), and Raynaud’s disease also showed no increased risk in the self-controlled risk interval analysis (adjusted relative risk 1.32, 95% confidence interval 0.30 to 5.76; rate ratios could not be calculated in the cohort analysis because the disease was not observed in the HPV unvaccinated group).

Several studies have focused on the association between HPV vaccination and migraine, and one study based on a passive surveillance database suggested migraine is a safety signal for HPV vaccination.3039 The studies did not confirm a clear relation between migraine and HPV vaccination, but further research is required to clarify the causal association behind the disease. A cohort study conducted in the Netherlands identified no significant increase in the risk of migraine after immunisation with the bivalent HPV vaccine.40 Although the study was limited by insufficient sample size and lack of power, this result was consistent with our result from subgroup analysis of the bivalent HPV vaccine (adjusted rate ratio 1.07, 95% confidence interval 0.96 to 1.20).

Strength and limitations of this study

By linking two different nationwide databases—a vaccine registry and health insurance database—our study provided results for the entire population of South Korea.41 Using this nationwide linked database allowed us to detect associations between HPV vaccine and rare or serious diseases without losing statistical power. In addition, as the South Korean vaccine registry compiles all immunisation data under the national immunisation programme, including vaccines against HPV, Japanese encephalitis, and tetanus, diphtheria, and acellular pertussis, the omission of information on vaccine use would be rare. Unlike previous studies, we adopted two complementary analyses—cohort and self-controlled risk interval. The cohort study design allowed us to provide more generalisable evidence and has higher statistical power, and the self-controlled risk interval study design enabled us to adjust time invariant confounders, including socioeconomic status, underlying chronic or congenital diseases, use of medical resources, and other unmeasured confounders.42 Specifically, in the cohort analysis we assessed the risk of serious adverse events after HPV vaccination by comparing the HPV vaccinated group with a comparator vaccinated group rather than with an unvaccinated group. Since the unvaccinated group includes those who decline vaccination, their clinical and personal characteristics and health seeking behaviour are likely to be different from those who are vaccinated.43 Thus, by comparing relatively homogeneous populations our approach secures internal validity and can produce reliable results by controlling for the unmasking phenomenon and fundamental differences between vaccinated and unvaccinated groups. In addition, our study assessed several serious adverse events suspected from previous studies and contains acute and insidious diseases through consultation of the medical advisory committee on vaccine safety, which enabled us to provide comprehensive safety information on HPV vaccines. Through our subgroup analyses we were also able to supply more information on the safety of each HPV vaccine subtype because most previous research was on vaccination using the quadrivalent HPV vaccine.

Our study has some limitations. Firstly, our case definition relied on ICD-10 codes recorded in health insurance data. As the diagnosis records using ICD-10 codes are not perfectly consistent with the medical records of patients, the discrepancy between health insurance data and medical records can overestimate or underestimate the real incidence of selected outcomes.4445 However, in a validation study regarding the National Health Information Database, the overall positive predictive value of the diagnosis was 82%,21 which enabled us to provide reliable information on vaccine safety. Secondly, because of limited access to medical records we used the date of diagnosis rather than the date of disease onset. The delay between disease onset and diagnosis would have contributed to the misclassification of cases if onset and diagnosis did not occur in the same risk period. This random misclassification bias would skew the results towards a null association. To resolve any delay from onset to diagnosis, we selected one year as the primary risk period and conducted sensitivity analyses with varying risk periods (90 days, 180 days, and 2 years). Thirdly, a possibility of misclassification about vaccine use remains. As both the bivalent and the quadrivalent HPV vaccines have been covered by funding from the South Korean government through the national immunisation programme, the registered vaccination records are complete, whereas vaccination records of a nonavalent HPV vaccine, which is not under the national immunisation programme, are not complete. A girl inoculated with a comparator vaccine and the nonavalent HPV vaccine is likely to be classified as being HPV unvaccinated. However, given the equivalent effect of the nonavalent and quadrivalent HPV vaccines,46 few girls are likely to be inoculated with the nonavalent HPV vaccine by paying additional costs. Fourthly, because many of the predefined outcomes are rare, statistical power can be limited for some outcomes in our analyses. Notwithstanding, as we used a nationwide linked database, a few events might indicate that the event is not likely to be relevant. Finally, as our study population comprised vaccinated adolescent girls inoculated with either the HPV vaccine or the comparator vaccine, the generalisability of the findings to the entire South Korean population could be limited. However, considering that the study population included more than 65% of all girls aged 11-14 years in South Korea, we believe that our findings are representative of the general population.

Conclusion and policy implications

From a public health perspective, HPV vaccination and screening were adapted worldwide as general strategies for the prevention of cervical cancer.47 Fear of possible adverse events is, however, one of the main determinants of HPV vaccine hesitancy and of undermining public trust in vaccines.94849 The inevitable consequence of HPV vaccine hesitancy in Japan was estimated to result in around 5000 deaths due to cervical cancer during 2013-19.50 In light of the impact of HPV vaccine hesitancy on cervical cancer mortality, establishing the safety of HPV vaccines and alleviating public concerns is imperative. Our study evaluated the association between HPV vaccination and serious adverse events (those related to autoimmunity or other organ system disorders) to provide safety information on HPV vaccination in South Korea. Of 33 predefined serious adverse events in our study, we found no significant increase in risk for 32, and these findings were consistent across the self-controlled risk interval, sensitivity, and subgroup analyses. Although a slightly increased risk was observed for migraine in the cohort analysis, this should be interpreted with caution considering the pathophysiology of migraine and the population of interest. The results for safety profiles of the HPV vaccine in adolescent girls in South Korea are consistent with the evidence in Western populations.

What is already known on this topic

  • Serious adverse events after immunisation with the human papillomavirus (HPV) vaccine are one of the major concerns and barriers to HPV vaccination

  • Although the real world evidence on the safety of the HPV vaccine is established for Western populations, evidence remains insufficient for Asian populations

What this study adds

  • This study found no associations between HPV vaccination and 33 predefined serious adverse events, using both cohort and self-controlled risk interval analyses

  • Our results are consistent with those of studies that have shown the safety of HPV vaccination in Western populations

Acknowledgments

We thank the Korea Disease Control and Prevention Agency for providing anonymous data from individual immunisation records, and the National Health Insurance Service for providing invaluable data.

Footnotes

  • Contributors: DY collected the data, performed statistical analyses, interpreted the data, and wrote the manuscript. JHL designed the overall research, interpreted the data, and contributed to the writing of the manuscript. HL collected the data and contributed to the statistical analyses and writing of the manuscript. JYS designed the study, supervised the statistical analyses and interpretation of data, and critically revised the manuscript. JYS, guarantor of the study, accepts full responsibility for the research, had access to the data, and controlled the decision to publish. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.

  • Funding: This research was supported by Government-wide R&D Fund project for infectious disease research (GFID), Republic of Korea (grant No HG18C0068). The funder was not involved in the study design; in the collection, analysis, and interpretation of data; in the writing of the report; and in the decision to submit the article for publication.

  • Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/coi_disclosure.pdf and declare: all authors had financial support from government-wide research and development fund project for infectious disease research; J-YS received grants from the Ministry of Food and Drug Safety; the Ministry of Health and Welfare; the National Research Foundation of Korea; the Government-wide R&D Fund project for infectious disease research; and pharmaceutical companies, including Amgen, Pfizer, Hoffmann-La Roche, Dong-A ST, and Yungjin; no other relationships or activities that could appear to have influenced the submitted work.

  • Ethical approval: This study was approved by the institutional review board of Sungkyunkwan University (No 2019-03-017), which waived the requirement for informed consent as only deidentified data were used in this study.

  • Data sharing: No additional data available.

  • The lead author (J-YS) affirms that this manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned (and, if relevant, registered) have been explained.

  • Dissemination to participants and related patient and public communities: We have no plans to disseminate the results of the research to study participants or the relevant patient community. As anonymised data were used in this study, the authors had no direct contact information of the study participants to disseminate the results of the research.

  • Provenance and peer review: Not commissioned; externally peer reviewed.

http://creativecommons.org/licenses/by-nc/4.0/

This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

References

  1. ↵
    1. de Martel C,
    2. Plummer M,
    3. Vignat J,
    4. Franceschi S
    . Worldwide burden of cancer attributable to HPV by site, country and HPV type. Int J Cancer2017;141:664-70. doi:10.1002/ijc.30716 pmid:28369882
    OpenUrlCrossRefPubMed
  2. ↵
    1. Ferlay J,
    2. Soerjomataram I,
    3. Dikshit R,
    4. et al
    . Cancer incidence and mortality worldwide: sources, methods and major patterns in GLOBOCAN 2012. Int J Cancer2015;136:E359-86. doi:10.1002/ijc.29210. pmid:25220842
    OpenUrlCrossRefPubMed
  3. ↵
    1. Markowitz LE,
    2. Tsu V,
    3. Deeks SL,
    4. et al
    . Human papillomavirus vaccine introduction--the first five years. Vaccine2012;30(Suppl 5):F139-48. doi:10.1016/j.vaccine.2012.05.039 pmid:23199957
    OpenUrlCrossRefPubMedWeb of Science
  4. ↵
    Global HPV vaccine introduction overview: projected and current national introductions, demonstration/pilot projects, gender-neutral vaccination programs, and global HPV vaccine introduction maps (2006-2022). 2019 https://path.azureedge.net/media/documents/GlobalHPVVaccineIntroOverview_Nov2019.pdf (accessed 26 September 2020).
  5. ↵
    Schedules and strategies for human papillomavirus (HPV) immunization: World Health Organization; 2016 [updated 27 September 2016. https://www.who.int/immunization/sage/meetings/2016/october/1_HPV_vaccine_background_document_27Sept2016.pdf?ua=1
  6. ↵
    1. Kim MA,
    2. Han GH,
    3. Kim JH,
    4. Seo K
    . Current Status of Human Papillomavirus Infection and Introduction of Vaccination to the National Immunization Program in Korea: an Overview. J Korean Med Sci2018;33:e331. doi:10.3346/jkms.2018.33.e331 pmid:30584412
    OpenUrlCrossRefPubMed
  7. ↵
    National Immunization Program for children 2019 [updated 16 July 2019]. http://www.cdc.go.kr/contents.es?mid=a30301140000 (accessed 1 March 2020).
  8. ↵
    Recommendation of HPV immunization for adolescent girls born between 2006 to 2007 Korea Centers for Disease Control and Prevention, 2019. https://nip.cdc.go.kr/irgd/support.do?service=getNewsView&strNum=252&PROSEQNUM=429&SEARCHTYPE=&SEARCHWARD= (accessed 20 March 2020).
  9. ↵
    1. Bonanni P,
    2. Zanella B,
    3. Santomauro F,
    4. Lorini C,
    5. Bechini A,
    6. Boccalini S
    . Safety and perception: What are the greatest enemies of HPV vaccination programmes?Vaccine2018;36:5424-9. doi:10.1016/j.vaccine.2017.05.071 pmid:28610824
    OpenUrlCrossRefPubMed
  10. ↵
    1. Jørgensen L,
    2. Gøtzsche PC,
    3. Jefferson T
    . Benefits and harms of the human papillomavirus (HPV) vaccines: systematic review with meta-analyses of trial data from clinical study reports. Syst Rev2020;9:43. doi:10.1186/s13643-019-0983-y pmid:32106879
    OpenUrlCrossRefPubMed
  11. ↵
    Survey of attitudes toward human papillomavirus vaccine vaccination 2017 https://www.cdc.go.kr/board/board.es?mid=a205010000000&bid=0015&act=view&list_no=75836&tag=&nPage=112 (accessed 4 June 2020).
  12. ↵
    1. Scheller NM,
    2. Svanström H,
    3. Pasternak B,
    4. et al
    . Quadrivalent HPV vaccination and risk of multiple sclerosis and other demyelinating diseases of the central nervous system. JAMA2015;313:54-61. doi:10.1001/jama.2014.16946 pmid:25562266
    OpenUrlCrossRefPubMedWeb of Science
  13. ↵
    1. Arnheim-Dahlström L,
    2. Pasternak B,
    3. Svanström H,
    4. Sparén P,
    5. Hviid A
    . Autoimmune, neurological, and venous thromboembolic adverse events after immunisation of adolescent girls with quadrivalent human papillomavirus vaccine in Denmark and Sweden: cohort study. BMJ2013;347:f5906. doi:10.1136/bmj.f5906 pmid:24108159
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Miranda S,
    2. Chaignot C,
    3. Collin C,
    4. Dray-Spira R,
    5. Weill A,
    6. Zureik M
    . Human papillomavirus vaccination and risk of autoimmune diseases: A large cohort study of over 2million young girls in France. Vaccine2017;35:4761-8. doi:10.1016/j.vaccine.2017.06.030 pmid:28750853
    OpenUrlCrossRefPubMed
  15. ↵
    1. Hviid A,
    2. Thorsen NM,
    3. Valentiner-Branth P,
    4. Frisch M,
    5. Mølbak K
    . Association between quadrivalent human papillomavirus vaccination and selected syndromes with autonomic dysfunction in Danish females: population based, self-controlled, case series analysis. BMJ2020;370:m2930. doi:10.1136/bmj.m2930 pmid:32878745
    OpenUrlAbstract/FREE Full Text
  16. ↵
    1. Hviid A,
    2. Svanström H,
    3. Scheller NM,
    4. Grönlund O,
    5. Pasternak B,
    6. Arnheim-Dahlström L
    . Human papillomavirus vaccination of adult women and risk of autoimmune and neurological diseases. J Intern Med2018;283:154-65. doi:10.1111/joim.12694 pmid:29044769
    OpenUrlCrossRefPubMed
  17. ↵
    1. Grimaldi-Bensouda L,
    2. Guillemot D,
    3. Godeau B,
    4. et al.,
    5. PGRx-AID Study Group
    . Autoimmune disorders and quadrivalent human papillomavirus vaccination of young female subjects. J Intern Med2014;275:398-408. doi:10.1111/joim.12155 pmid:24206418
    OpenUrlCrossRefPubMed
  18. ↵
    1. Yih WK,
    2. Greene SK,
    3. Zichittella L,
    4. et al
    . Evaluation of the risk of venous thromboembolism after quadrivalent human papillomavirus vaccination among US females. Vaccine2016;34:172-8. doi:10.1016/j.vaccine.2015.09.087 pmid:26549364
    OpenUrlCrossRefPubMed
  19. ↵
    1. Willame C,
    2. Rosillon D,
    3. Zima J,
    4. et al
    . Risk of new onset autoimmune disease in 9- to 25-year-old women exposed to human papillomavirus-16/18 AS04-adjuvanted vaccine in the United Kingdom. Hum Vaccin Immunother2016;12:2862-71. doi:10.1080/21645515.2016.1199308 pmid:27428517
    OpenUrlCrossRefPubMed
  20. ↵
    1. Cheol Seong S,
    2. Kim YY,
    3. Khang YH,
    4. et al
    . Data Resource Profile: The National Health Information Database of the National Health Insurance Service in South Korea. Int J Epidemiol2017;46:799-800. doi:10.1093/ije/dyw253. pmid:27794523
    OpenUrlCrossRefPubMed
  21. ↵
    1. Park EJS,
    2. Jeon S,
    3. Lee S,
    4. Lee J,
    5. Choi D
    . Report of the evaluation for validity of discharged diagnoses in Korean Health Insurance database.Health Insurance Review and Assessment Service, 2017.
  22. ↵
    1. Jacobsen SJ,
    2. Sy LS,
    3. Ackerson BK,
    4. et al
    . An unmasking phenomenon in an observational post-licensure safety study of adolescent girls and young women. Vaccine2012;30:4585-7. doi:10.1016/j.vaccine.2012.04.103 pmid:22580356
    OpenUrlCrossRefPubMed
  23. ↵
    1. Park HJ,
    2. Kim SJ,
    3. Song R,
    4. et al
    . A 6-year Prospective, Observational, Multi-Center Post-Marketing Surveillance of the Safety of Tetanus Toxoid, Reduced Diphtheria Toxoid, and Acellular Pertussis (Tdap) Vaccine in Korea. J Korean Med Sci2019;34:e105. doi:10.3346/jkms.2019.34.e105 pmid:30940999
    OpenUrlCrossRefPubMed
  24. ↵
    1. Schuller E,
    2. Klingler A,
    3. Dubischar-Kastner K,
    4. Dewasthaly S,
    5. Müller Z
    . Safety profile of the Vero cell-derived Japanese encephalitis virus (JEV) vaccine IXIARO(®). Vaccine2011;29:8669-76. doi:10.1016/j.vaccine.2011.08.117 pmid:21907747
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    1. Chokephaibulkit K,
    2. Houillon G,
    3. Feroldi E,
    4. Bouckenooghe A
    . Safety and immunogenicity of a live attenuated Japanese encephalitis chimeric virus vaccine (IMOJEV®) in children. Expert Rev Vaccines2016;15:153-66. doi:10.1586/14760584.2016.1123097 pmid:26588242
    OpenUrlCrossRefPubMed
  26. ↵
    1. Kim CJ,
    2. Song R,
    3. Chen J,
    4. et al
    . Six-year multi-centre, observational, post-marketing surveillance of the safety of the HPV-16/18 AS04-adjuvanted vaccine in women aged 10-25 years in Korea. Pharmacoepidemiol Drug Saf2017;26:837-42. doi:10.1002/pds.4175 pmid:28266092
    OpenUrlCrossRefPubMed
  27. ↵
    1. Markowitz LE,
    2. Hariri S,
    3. Unger ER,
    4. Saraiya M,
    5. Datta SD,
    6. Dunne EF
    . Post-licensure monitoring of HPV vaccine in the United States. Vaccine2010;28:4731-7. doi:10.1016/j.vaccine.2010.02.019 pmid:20188681
    OpenUrlCrossRefPubMedWeb of Science
  28. ↵
    1. Angelo MG,
    2. Zima J,
    3. Tavares Da Silva F,
    4. Baril L,
    5. Arellano F
    . Post-licensure safety surveillance for human papillomavirus-16/18-AS04-adjuvanted vaccine: more than 4 years of experience. Pharmacoepidemiol Drug Saf2014;23:456-65. doi:10.1002/pds.3593 pmid:24644078
    OpenUrlCrossRefPubMed
  29. ↵
    1. Arana JE,
    2. Harrington T,
    3. Cano M,
    4. et al
    . Post-licensure safety monitoring of quadrivalent human papillomavirus vaccine in the Vaccine Adverse Event Reporting System (VAERS), 2009-2015. Vaccine2018;36:1781-8. doi:10.1016/j.vaccine.2018.02.034 pmid:29477308
    OpenUrlCrossRefPubMed
  30. ↵
    1. Ran J,
    2. Yang JY,
    3. Lee JH,
    4. Kim HJ,
    5. Choi JY,
    6. Shin JY
    . Signal detection of human papillomavirus vaccines using the Korea Adverse Events Reporting System database, between 2005 and 2016. Int J Clin Pharm2019;41:1365-72. doi:10.1007/s11096-019-00881-9 pmid:31313003
    OpenUrlCrossRefPubMed
  31. ↵
    1. Ozawa K,
    2. Hineno A,
    3. Kinoshita T,
    4. Ishihara S,
    5. Ikeda SI
    . Suspected Adverse Effects After Human Papillomavirus Vaccination: A Temporal Relationship Between Vaccine Administration and the Appearance of Symptoms in Japan. Drug Saf2017;40:1219-29. doi:10.1007/s40264-017-0574-6 pmid:28744844
    OpenUrlCrossRefPubMed
  32. ↵
    1. Poddighe D,
    2. Castelli L,
    3. Marseglia GL,
    4. Bruni P
    . A sudden onset of a pseudo-neurological syndrome after HPV-16/18 AS04-adjuvated vaccine: might it be an autoimmune/inflammatory syndrome induced by adjuvants (ASIA) presenting as a somatoform disorder?Immunol Res2014;60:236-46. doi:10.1007/s12026-014-8575-3 pmid:25388965
    OpenUrlCrossRefPubMed
  33. ↵
    1. Baker MA,
    2. Lieu TA,
    3. Li L,
    4. et al
    . A vaccine study design selection framework for the postlicensure rapid immunization safety monitoring program. Am J Epidemiol2015;181:608-18. doi:10.1093/aje/kwu322 pmid:25769306
    OpenUrlCrossRefPubMed
  34. ↵
    1. Chao C,
    2. Jacobsen SJ
    . Evaluation of autoimmune safety signal in observational vaccine safety studies. Hum Vaccin Immunother2012;8:1302-4. doi:10.4161/hv.21268 pmid:22871958
    OpenUrlCrossRefPubMed
  35. ↵
    1. Pavlović JM,
    2. Stewart WF,
    3. Bruce CA,
    4. et al
    . Burden of migraine related to menses: results from the AMPP study. J Headache Pain2015;16:24. doi:10.1186/s10194-015-0503-y pmid:25902814
    OpenUrlCrossRefPubMed
  36. ↵
    1. Lee MH,
    2. Kim SH,
    3. Oh M,
    4. Lee KW,
    5. Park MJ
    . Age at menarche in Korean adolescents: trends and influencing factors. Reprod Health2016;13:121. doi:10.1186/s12978-016-0240-y pmid:27662834
    OpenUrlCrossRefPubMed
  37. ↵
    1. Chao C,
    2. Klein NP,
    3. Velicer CM,
    4. et al
    . Surveillance of autoimmune conditions following routine use of quadrivalent human papillomavirus vaccine. J Intern Med2012;271:193-203. doi:10.1111/j.1365-2796.2011.02467.x pmid:21973261
    OpenUrlCrossRefPubMed
  38. ↵
    1. Willame C,
    2. Gadroen K,
    3. Bramer W,
    4. Weibel D,
    5. Sturkenboom M
    . Systematic Review and Meta-analysis of Postlicensure Observational Studies on Human Papillomavirus Vaccination and Autoimmune and Other Rare Adverse Events. Pediatr Infect Dis J2020;39:287-93. doi:10.1097/INF.0000000000002569 pmid:31876615
    OpenUrlCrossRefPubMed
  39. ↵
    1. Brinth L,
    2. Theibel AC,
    3. Pors K,
    4. Mehlsen J
    . Suspected side effects to the quadrivalent human papilloma vaccine. Dan Med J2015;62:A5064.pmid:25872549
    OpenUrlPubMed
  40. ↵
    1. Schurink-Van’t Klooster TM,
    2. de Ridder MA,
    3. Kemmeren JM,
    4. et al
    . Examining a possible association between human papilloma virus (HPV) vaccination and migraine: results of a cohort study in the Netherlands. Eur J Pediatr2015;174:641-9. doi:10.1007/s00431-014-2444-x pmid:25367054
    OpenUrlCrossRefPubMed
  41. ↵
    1. Kim L,
    2. Kim J-A,
    3. Kim S
    . A guide for the utilization of Health Insurance Review and Assessment Service National Patient Samples. Epidemiol Health2014;36:e2014008-08. doi:10.4178/epih/e2014008 pmid:25078381
    OpenUrlCrossRefPubMed
  42. ↵
    1. Li R,
    2. Stewart B,
    3. Weintraub E
    . Evaluating efficiency and statistical power of self-controlled case series and self-controlled risk interval designs in vaccine safety. J Biopharm Stat2016;26:686-93. doi:10.1080/10543406.2015.1052819 pmid:26098696
    OpenUrlCrossRefPubMed
  43. ↵
    1. Smith PJ,
    2. Chu SY,
    3. Barker LE
    . Children who have received no vaccines: who are they and where do they live?Pediatrics2004;114:187-95. doi:10.1542/peds.114.1.187 pmid:15231927
    OpenUrlAbstract/FREE Full Text
  44. ↵
    1. Lim S-S,
    2. Lee W,
    3. Kim Y-K,
    4. et al
    . The cumulative incidence and trends of rare diseases in South Korea: a nationwide study of the administrative data from the National Health Insurance Service database from 2011-2015. Orphanet J Rare Dis2019;14:49. doi:10.1186/s13023-019-1032-6 pmid:30777110
    OpenUrlCrossRefPubMed
  45. ↵
    1. O’Malley KJ,
    2. Cook KF,
    3. Price MD,
    4. Wildes KR,
    5. Hurdle JF,
    6. Ashton CM
    . Measuring diagnoses: ICD code accuracy. Health Serv Res2005;40:1620-39. doi:10.1111/j.1475-6773.2005.00444.x pmid:16178999
    OpenUrlCrossRefPubMedWeb of Science
  46. ↵
    1. Bergman H,
    2. Buckley BS,
    3. Villanueva G,
    4. et al
    . Comparison of different human papillomavirus (HPV) vaccine types and dose schedules for prevention of HPV-related disease in females and males. Cochrane Database Syst Rev2019;2019:CD013479. doi:10.1002/14651858.CD013479 pmid:31755549
    OpenUrlCrossRefPubMed
  47. ↵
    1. Jit M,
    2. Brisson M,
    3. Portnoy A,
    4. Hutubessy R
    . Cost-effectiveness of female human papillomavirus vaccination in 179 countries: a PRIME modelling study. Lancet Glob Health2014;2:e406-14. doi:10.1016/S2214-109X(14)70237-2 pmid:25103394
    OpenUrlCrossRefPubMed
  48. ↵
    1. Stretch R,
    2. Roberts SA,
    3. McCann R,
    4. et al
    . Parental attitudes and information needs in an adolescent HPV vaccination programme. Br J Cancer2008;99:1908-11. doi:10.1038/sj.bjc.6604766 pmid:18985038
    OpenUrlCrossRefPubMedWeb of Science
  49. ↵
    1. Egawa-Takata T,
    2. Ueda Y,
    3. Morimoto A,
    4. et al
    . Survey of Japanese mothers of daughters eligible for human papillomavirus vaccination on attitudes about media reports of adverse events and the suspension of governmental recommendation for vaccination. J Obstet Gynaecol Res2015;41:1965-71. doi:10.1111/jog.12822 pmid:26420072
    OpenUrlCrossRefPubMed
  50. ↵
    1. Simms KT,
    2. Hanley SJB,
    3. Smith MA,
    4. Keane A,
    5. Canfell K
    . Impact of HPV vaccine hesitancy on cervical cancer in Japan: a modelling study. Lancet Public Health2020;5:e223-34. doi:10.1016/S2468-2667(20)30010-4 pmid:32057317
    OpenUrlCrossRefPubMed
Back to top