ORIGINAL RESEARCH article

Front. Med., 01 September 2021

Sec. Infectious Diseases: Pathogenesis and Therapy

Volume 8 - 2021 | https://doi.org/10.3389/fmed.2021.616873

The Potential Role of Human Papillomavirus Infection in Bell's Palsy: A Hypothesis-Generating Study Based on a Nationwide Cohort

  • 1. Department of Neurology, Kaohsiung Medical University Hospital, Kaohsiung, Taiwan

  • 2. Department of Neurology, School of Medicine, College of Medicine, Kaohsiung Medical University, Kaohsiung, Taiwan

  • 3. Neuroscience Research Center, Kaohsiung Medical University, Kaohsiung, Taiwan

  • 4. Shu-Zen Junior College of Medicine and Management, Kaohsiung, Taiwan

  • 5. Department of Recreation and Sports Management, Tajen University, Pingtung, Taiwan

  • 6. Department of Physical Medicine and Rehabilitation, Kaohsiung Veterans General Hospital, Pingtung, Taiwan

  • 7. Graduate Institute of Bioresources, National Pingtung University of Science and Technology, Pingtung, Taiwan

  • 8. Department of Emergency Medicine, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan

  • 9. Institute of Medicine, Chung Shan Medical University, Taichung, Taiwan

  • 10. Management Office for Health Data, China Medical University Hospital, Taichung, Taiwan

  • 11. College of Medicine, China Medical University, Taichung, Taiwan

  • 12. Institute of Public Health (Biostatistics), National Yangming University, Taipei, Taiwan

  • 13. Department of Internal Medicine, Kaohsiung Municipal United Hospital, Kaohsiung, Taiwan

  • 14. College of Health and Nursing, Meiho University, Pingung, Taiwan

  • 15. National Yang-Ming University, School of Medicine, Taipei, Taiwan

  • 16. Division of Allergy, Immunology and Rheumatology, Chung Shan Medical University Hospital, Taichung, Taiwan

  • 17. Graduate Institute of Integrated Medicine, China Medical University, Taichung, Taiwan

Abstract

Objective: Our purpose was to investigate whether people with a previous human papillomavirus (HPV) infection were associated with an increased risk of Bell's palsy (BP).

Methods: By using Taiwan population-based data, patients aged > 18 years with HPV infection (n = 22,260) from 2000 to 2012 were enrolled and compared with control subjects who had never been diagnosed with an HPV infection at a 1:4 ratio matched by sex, age, index date, and co-morbidities (n = 89,040). The index date was the first date of HPV diagnosis. All the patients were tracked until the occurrence of BP. Cox proportional hazards regression was applied to estimate the hazard ratios (HRs) for the development of BP in both groups.

Results: The HPV group had 1.25 [95% confidence interval (CI) = 1.03–1.51] times higher risk of BP compared with the non-HPV group after adjusting for sex, age, and co-morbidities. The association of HPV and BP was significant in the sensitivity analyses. In the subgroup analysis, the impact of HPV infection on the risk of BP was more pronounced in the elderly > 50 years [adjusted hazard ratio (aHR) =1.86; 95% CI = 1.37–2.52], hypertension (aHR = 1.65; 95% CI = 1.17–2.31), and chronic obstructive pulmonary disease (aHR = 2.14, 95% CI 1.333.43) subgroups.

Conclusions: Patients with HPV infection have a higher risk of subsequent BP compared with non-HPV patients. More rigorous studies are needed to confirm if and how specific HPV genotypes are associated with BP and the possible role of vaccines in disease prevention.

Introduction

Bell's palsy (BP), defined as acute peripheral facial nerve palsy of unknown cause, is the most common acute mono-neuropathy (1). This disease affects 11–40 people per 100,000 annually (2). Typically, BP results in a partial or complete inability to move the affected side of the facial muscles. There is no race, geographic, or gender preference (3). The risk factors include diabetes, hypertension, and pregnancy (3, 4). Although the facial paresis or paralysis usually resolves within several months, permanent symptoms may occur and affect the quality of life of the patient (5). Indeed facial paralysis was significantly associated with increased depression, anxiety, and a worse quality of life (6). The cause of BP remains unknown, and many events, such as ischemia, viral infection, acute cold exposure, and anatomic causes, have been proposed as causes of BP (7). Investigators proposed that immune alterations related to a previous viral infection have an important role in the etiopathogenesis of BP (8). The reasons for this include alterations in lymphocyte subsets in the peripheral blood during the acute stage of the disease (9, 10). Additionally, lymphocytes from patients with BP were stimulated in the presence of neuritogenic basic protein from the peripheral nerve myelin sheath (11). A similar pattern of response was also demonstrated in patients with Guillain–Barre syndrome (12, 13), an autoimmune demyelination neuritis triggered by a preceding infection (14). However, serological evidence and histopathology of facial nerves support the infectious nature of BP. Herpes simplex viruses 1 and 2 and varicella-zoster viruses are the most common viruses associated with BP (15). Other infections, such as Epstein–Barr virus, cytomegalovirus, and human immunodeficiency virus, have also been reported to be associated with BP (1517). Human papillomaviruses (HPV) are small double-stranded DNA viruses known as risk factors for multiple cancers (1820) and are also associated with autoimmune diseases, such as systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) (2123). The shared pathogenesis between these autoimmune diseases and BP is based on a viral infection and a subsequent autoimmune reaction (24, 25). These findings suggest a possible association between HPV infection and BP. However, there is no epidemiologic evidence to support this hypothesis. Thus, we conducted this original nationwide population-based cohort study to assess the incidence rate and risk of BP in people with a history of HPV infection when compared with those without HPV infection.

Methods

Data Source

The data in this study were obtained from Taiwan's National Health Insurance Research Database (NHIRD), which contains registration files and original claim data, such as diagnoses coding, out-patient visits, hospitalization records, and medication and personal information, for almost 99% of the population of Taiwan. Moreover, ~93% of healthcare facilities in Taiwan have signed contracts with the Bureau of National Health Insurance. To ensure the accuracy of these data, all claims data were validated by the Bureau of National Health Insurance. The data used in the study were obtained from the Longitudinal Health Insurance Database (LHID), a subset of NHIRD. The LHID contains original claims data of 1,000,000 people randomly sampled from the 1996 to 2013 registries for beneficiaries of the NHIRD, and the dataset has been used for medical and epidemiological research (26). There were no statistically significant differences in the distribution of age, sex, or healthcare costs between the two databases. With a strict definition of the study group selection, the positive predictive value of the NHIRD claim data can be up to 84.6%. The authority replaced the original identification numbers with surrogate numbers before the data were released to protect the privacy of the people. The Institutional Review Board of the China Medical University in Taiwan approved this study [CMUH140-REC2-115 (CR-4)].

Data Availability Statement

Data are available from the NHIRD published by Taiwan National Health Insurance (NHI) Bureau. Due to legal restrictions imposed by the government of Taiwan in relation to the “Personal Information Protection Act,” data cannot be made publicly available. Requests for data can be sent as a formal proposal to the NHIRD (http://nhird.nhri.org.tw).

Study Population

To ensure the accuracy of patient information during the study period, patients with missing demographic data and those who died before 2003 were excluded from this analysis. HPV infection was defined according to the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM) codes. We identified patients who had been diagnosed with HPV infection (ICD-9-CM codes 079.4, 078.1, 078.10–078.12, 078.19, 795.05, 795.09, 795.15, 795.19, 796.75, and 796.79) between 2000 and 2012 as the HPV group. The index date was defined as the date of HPV infection diagnosis. For each patient in the study group, the corresponding controls were selected at a ratio of 1:4 based on the nearest propensity score. The propensity scores were calculated using the probability of HPV assignment by using a logistic regression model and included the following baseline variables: year of index date; age; sex; comorbidities of hypertension, diabetes, hyperlipidemia, chronic kidney disease, and chronic obstructive pulmonary disease (COPD); alcohol-related illness; herpes simplex virus type 1; herpes zoster; systemic lupus erythematosus; rheumatoid arthritis; Sjögren's syndrome; and multiple sclerosis. Figure 1 shows the exclusion criteria and process of assembling the HPV group and the non-HPV group step by step. To ameliorate potential confounders in the analysis, patients with well-known risk factors contributing to BP, such as stroke, otitis media, Guillain–Barre syndrome, and head and neck cancer, or pregnancy 1 year before BP were excluded (Figure 1). The Supplementary Figure 1 shows an overview of the study. Data from 1996 to 2000 in the dataset were used as a wash-out period to ensure that the patients in both groups did not have a record of BP or HPV infection before study entry. To strengthen the validity of coding for HPV-associated diseases, patients with an HPV infection who did not receive any HPV-related treatment procedures in the 3 months after the index date were excluded. Such treatment procedure codes included “electrocauterization for condyloma (50005),” “condyloma, excision, and electrocauterization (55008),” “CO2 laser operation (62020),” “chemosurgery, condyloma (50015),” “electro-cauterization, simple (51005),” “electro-cauterization, complicated (51006),” “liquid nitrogen cryosurgery (51017),” “cryotherapy, simple, including CO2 freezing and liquid nitrogen (51021C),” and “cryotherapy, complicated, including CO2 freezing and liquid nitrogen (51022).” All participants were followed until the presence of BP, dropping out from the insurance program, death, or the end of the study (December 31, 2013). Finally, 22,260 subjects were included in the HPV group and 89,040 subjects served as the control group.

Figure 1

Covariates and Outcome

The co-morbidities analyzed in this study were hypertension (ICD-9-CM codes 401–405), diabetes mellitus (ICD-9-CM code 250**), hyperlipidemia (ICD-9-CM code 272), chronic obstructive pulmonary disease (ICD-9-CM codes 490–492 and 493–496), and herpes zoster (ICD-9-CM code 053**). All subjects in both groups were tracked from the index date to the first BP event. The end-point of this study was the occurrence of BP (ICD-9-CM code 351.0) (27). To improve the reliability of ICD coding, a BP event was defined as a patient having a record of BP diagnosis in three or more ambulatory visits or at least one admission within a year. The co-morbidities were obtained in the LHID 2 years before the index date of HPV diagnosis.

Statistical Analysis

First, we applied the chi-square test to compare the distribution of age, sex, and baseline co-morbidities between the two groups, with or without previous HPV infection. We compared the mean age using Student's t-test. The incidence density of BP per 1,000 person-years was calculated for both groups.

Second, we utilized multivariable Cox proportional hazard regression models to estimate the crude hazard ratios and adjusted HRs (aHRs) of BP and related 95% confidence intervals (CIs) of the HPV group compared to the non-HPV group adjusting for sex, age, and the aforementioned co-morbidities. We calculated the receiver operating characteristic curve for the main regression model and its area under the curve to assess the fit of the model. Two sensitivity analyses were conducted, including data on concomitant medication with anti-inflammatory effects and medical conditions that might be associated with immune system functions. In the first sensitivity analysis, potential anti-inflammatory medications were included in the covariates (Supplementary Tables 1, 2); in the second sensitivity analysis, we changed different coexisting medical co-morbidities as covariates (Supplementary Tables 3, 4).

Third, the Kaplan–Meier method was used to plot the cumulative incidence curves of BP for both groups. Differences between the two groups were evaluated by log-rank test.

Fourth, to reveal the association of sex, age, and follow-up time on the incidence rate of BP, we applied multivariable Cox regression models adjusted for sex, age, and co-morbidities. The HRs adjusted for covariates were calculated for males and females; age (years) <30, 31–40, 41–50, and above 50; and follow-up time <1 year, 1–5 years, and >5 years. All data analyses were performed with SAS® (version 9.4; SAS Institute, Inc., Cary, NC, USA). The statistical significance level was set at P-value < 0.05 by two-tailed test.

Results

Cohort Characteristics

The study consisted of 89,040 non-HPV subjects and 22,260 HPV subjects. The baseline characteristics of the two groups are listed in Table 1. In both groups, there was a slight female predominance (52%), and subjects aged below 30 years old (35%) comprised the largest age group.

Table 1

Non-HPVHPV
N=89,040N=22,260
Variablesn%n%SMD
Sex<0.001
   Female46,64852%11,66252%
   Male42,39248%10,59848%
Age, year
    <3030,76035%7,69035%<0.001
   30–391,807620%4,51920%<0.001
   40–4916,49219%4,12319%<0.001
   50–5911,48813%2,87213%<0.001
   60–696,0207%1,5057%<0.001
   70–794,4445%1,1115%<0.001
   ≥801,7602%4402%<0.001
   Mean (SD)40.4(16.6)40.4(16.6)0.003
Comorbidities
   Hypertension14,80417%4,02918%0.04
   Diabetes7,3808%1,9919%0.02
   Hyperlipidemia12,85614%4,04818%0.10
   CKD9341.0%3251.5%0.04
   COPD7,3178.2%2,29110%0.07
   Alcohol-related illness1,8252.0%4842.2%0.01
   HSV14,7065.3%1,9618.8%0.14
   Herpes zoster2,4692.8%9054.1%0.07
   SLE2410.27%970.44%0.03
Rheumatoid arthritis1,4561.6%4552.0%0.03
Sjögren's syndrome6660.7%3021.4%0.06
multiple sclerosis270.03%30.01%0.01

Baseline patient characteristics.

CKD, chronic kidney disease; COPD, chronic obstruction pulmonary disease; HSV1, genital herpes simplex virus; SLE, systemic lupus erythematosus.

Risk of BP

Table 2 shows the incidence rate and HRs for BP. There were 569 cases of newly diagnosed BP. The crude HR of HPV was 1.29 (95% CI = 1.07–1.57); the aHR after adjusting for demographic variables, including sex, age, and co-morbidities, was 1.25 (95% CI = 1.03–1.51). The risk of BP was higher in subjects above the age of 30 years, especially those in the age groups 30–39 and 40–49 years, with 1.36 (95% CI = 1.06–1.73) and 1.42 (95% CI = 1.10–1.81) times higher risk compared with individuals under 30 years old. The association between HPV infection and BP remained significant after considering coexisting medical conditions and medications (Supplementary Tables 2, 4). There were no statistically significant differences in risk between males and females. In Figure 2, the cumulative incidence curves plotted by the Kaplan–Meier method shows that the HPV group had a higher cumulative incidence of BP than that of the non-HPV group (log-rank test P-value < 0.008).

Table 2

Bell's palsy
VariablesnPYIRcHR(95% CI)aHR(95% CI)
HPV
 No429488,7420.881.001.00
   Yes140123,1731.141.29(1.07, 1.57)**1.25(1.03, 1.51)*
Sex
   Female301322,8580.931.00
   Male268289,0570.930.99(0.84, 1.17)
Age, year
    <30145222,6480.651.001.00
   30–39116126,6370.921.41(1.10, 1.80)**1.36(1.06, 1.73)*
   40–49119113,4741.051.62(1.27, 2.06)***1.42(1.10, 1.81)**
   ≥50189149,1571.271.96(1.58, 2.44)***1.26(0.96, 1.65)
Comorbidities
   Hypertension
No417519,2550.801.001.00
Yes15292,6601.642.06(1.71, 2.48)***1.51(1.18, 1.93)**
   Diabetes
No489566,4740.861.001.00
Yes8045,4411.762.05(1.62, 2.6)***1.29(0.98, 1.71)
   Hyperlipidemia
No436527,7170.831.001.00
Yes13384,1981.581.92(1.58, 2.34)***1.26(0.99, 1.61)
   CKD
No563606,8200.931.00
Yes65,0951.181.28(0.57, 2.87)
   COPD
No498563,8320.881.001.00
Yes7148,0831.481.68(1.31, 2.15)***1.19(0.91, 1.55)
   Alcohol-related illness
No555601,5470.921.00
Yes1410,3681.351.48(0.87, 2.52)
   HSV1
No547580,0420.941.00
Yes2231,8730.690.74(0.48, 1.13)
   Herpes zoster
No543596,5480.911.001.00
Yes2615,3671.691.88(1.27, 2.78)**1.45(0.97, 2.16)
   SLE
No567610,1830.931.00
Yes21,7321.151.25(0.31, 5.01)
   Rheumatoid arthritis
No558602,5020.931.00
Yes119,4131.171.27(0.7, 2.31)
   Sjögren's syndrome
No564607,6050.931.00
Yes54,3101.161.27(0.52, 3.05)
   Multiple sclerosis
No569611,7850.931.00
Yes01300.000.00(0, Inf)

Bell's palsy incidence rate and risk factors.

*

p-value < 0.05;

**

p-value < 0.01;

***

p-value < 0.001.

PY, person-years; IR, incidence rate (per 1,000 person-years); cHR, crude hazard ratio; aHR, adjusted hazard ratio; CKD, chronic kidney disease; COPD, chronic obstruction pulmonary disease; HSV1, genital herpes simplex virus; SLE, systemic lupus erythematosus.

Adjusted by age, hypertension, diabetes, hyperlipidemia, COPD and herpes zoster.

Figure 2

Subgroup Analysis

Table 3 shows the incidence rate and aHRs for BP among different subgroups stratified by sex, age, and follow-up time. The interaction between sex and HPV infection on the risk of BP was not significant (P-value for interaction = 0.55). The proportion of individuals with BP was highest in patients aged > 50 years in both groups. When we compared the risk of BP in both groups according to age, those aged > 50 years had the highest statistically significant risk of developing BP (aHR = 1.86; 95% CI = 1.37–2.52). The interaction between age and HPV infection on the risk of BP was significant (P-value for interaction = 0.01). In the comorbidity subgroup analysis, a higher relative risk for BP in the HPV group compared with the control group was observed for the subgroup with hypertension (aHR 1.65; 95% CI = 1.17–2.31; P-value for interaction = 0.04) and COPD (aHR 2.14; 95% CI = 1.33–3.43; P-value for interaction = 0.01). When stratified by follow-up time (Table 4), the HPV group had a significantly higher aHR for BP > 5 year (aHR = 1.61; 95% CI = 1.15–2.26). No statistically significant risks were found in follow-up time <1 year and 1–5 years.

Table 3

Non-HPVHPV
VariablesEventPYIREventPYIRcHR(95% CI)aHR(95% CI)P-value for interaction
Sex0.55
   Female230257,9290.897164,9291.091.23(0.94, 1.6)1.19(0.91, 1.56)
   Male199230,8130.866958,2441.181.37(1.04, 1.81)*1.33(1.01, 1.75)*
Age, year0.01
    <30118178,0910.662744,5570.610.91(0.6, 1.39)0.91(0.60, 1.39)
   30–3992101,1900.912425,4460.941.04(0.66, 1.63)1.00(0.64, 1.57)
   40–499390,6191.032622,8551.141.11(0.72, 1.71)1.07(0.69, 1.66)
   ≥50126118,8421.066330,3152.081.95(1.44, 2.65)***1.86(1.37, 2.52)***
Co-morbidities
   Hypertension0.04
No327416,6030.7890102,6520.881.12(0.88, 1.41)1.10(0.87, 1.39)
Yes10272,1391.415020,5212.441.71(1.22, 2.4)**1.65(1.17, 2.31)**
   Diabetes0.43
No373453,2250.82116113,2491.021.24(1.01, 1.53)*1.21(0.98, 1.49)
Yes5635,5181.58249,9242.421.53(0.95, 2.47)1.50(0.93, 2.42)
   Hyperlipidemia0.47
No338425,3110.7998102,4060.961.20(0.96, 1.51)1.20(0.96, 1.51)
Yes9163,4311.434220,7672.021.41(0.97, 2.03)1.38(0.95, 1.99)
   CKD0.36
No426485,0240.88137121,7951.121.28(1.06, 1.55)*1.24(1.02, 1.5)*
Yes33,7180.81313772.182.68(0.54, 13.29)2.58(0.51, 13.17)
   COPD0.01
No388452,6200.86110111,2120.991.15(0.93, 1.43)1.13(0.91, 1.39)
Yes4136,1231.143011,9612.512.21(1.38, 3.54)***2.14(1.33, 3.43)**
   Alcohol-related illness0.33
No417480,5690.87138120,9781.141.31(1.08, 1.59)**1.27(1.04, 1.54)*
Yes128,1731.4722,1950.910.62(0.14, 2.77)0.62(0.14, 2.79)
   HSV10.62
No414466,6560.89133113,3861.171.32(1.09, 1.61)**1.28(1.05, 1.56)*
Yes1522,0860.6879,7870.721.07(0.44, 2.62)1.01(0.41, 2.49)
   Herpes zoster0.62
No413477,7670.86130118,7811.091.27(1.04, 1.54)*1.23(1.01, 1.5)*
Yes1610,9751.46104,3922.281.56(0.71, 3.44)1.52(0.69, 3.37)
   SLE0.94
No427487,5250.88140122,6581.141.30(1.08, 1.58)**1.26(1.04, 1.52)*
Yes21,2181.6405140.000.00(0, Inf)0.00(0, Inf)
   Rheumatoid arthritis0.89
No421481,6040.87137120,8981.131.30(1.07, 1.57)**1.25(1.03, 1.52)*
Yes87,1381.1232,2751.321.20(0.32, 4.51)1.06(0.28, 4.01)
   Sjögren's syndrome0.92
No426485,8070.88138121,7981.131.29(1.07, 1.57)**1.25(1.03, 1.51)*
Yes32,9351.0221,3751.451.42(0.24, 8.47)1.25(0.21, 7.62)
   Multiple sclerosis1.00
No429488,6280.88140123,1571.141.29(1.07, 1.57)**1.25(1.03, 1.51)*
Yes01140.000160.00

Subgroup analysis.

*

p-value < 0.05;

**

p-value < 0.01;

***

p-value < 0.001.

PY, person-years; IR, incidence rate (per 1,000 person-years); cHR, crude hazard ratio; aHR, adjusted hazard ratio; CKD, chronic kidney disease; COPD, chronic obstruction pulmonary disease; HSV1, genital herpes simplex virus; SLE, systemic lupus erythematosus.

Adjusted by age, hypertension, diabetes, hyperlipidemia, COPD and herpes zoster.

Table 4

Non-HPVHPV
EventPYIREventPYIRcHR(95% CI)aHR(95% CI)P-value for interaction
Follow-up time, year0.28
    <18587,7690.982121,9940.951.00(0.62, 1.61)0.97(0.60, 1.57)
   1–5231360,5380.647190,8330.781.22(0.93, 1.59)1.17(0.90, 1.53)
   >5113369,3080.314893,3910.511.68(1.2, 2.36)**1.61(1.15, 2.26)**

Bell's palsy in different follow-up period.

**

p-value < 0.01.

PY, person-years; IR, incidence rate (per 1,000 person-years); cHR, crude hazard ratio; aHR, adjusted hazard ratio.

Adjusted by age, hypertension, diabetes, hyperlipidemia, COPD and herpes zoster.

Receiver Operating Characteristic Curve

Finally, we calculated the receiver operating characteristic curve for the main regression model in years 1, 5, and 10 and the area under the curve to assess the fit of the model. The area under the curve increased with the increase in time, indicating that the model could predict the risk of BP (Supplementary Figure 2).

Discussion

In this population study, people with HPV infection had a significantly increased risk of BP compared with the sex-, age-, and index date-matched non-HPV group. After adjusting for sex, age, and co-morbidities, the HPV group had a 25% increased risk of developing BP compared with the non-HPV group. The consistency of a positive association from different sensitivity analyses validated our findings. In our stratified analysis, we found a prominent interaction effect between age and HPV infection on the risk of new-onset BP. Among participants aged 50 years and older, having an HPV infection associated with a higher risk of developing BP than those without an HPV infection. We also observed significant interactions between HPV infection and some comorbidities (hypertension and COPD). The simultaneous presence of HPV infection and hypertension or COPD revealed a higher aHR than HPV infection alone.

Our study showed consistent findings in reporting the epidemiology of BP with previous studies. In Table 2, the incidence rate was similar in women and men, which was compatible with a previous study in the UK (28). Our study showed that patients aged above 30 years were more susceptible to BP. This result is in agreement with a previous study that showed that BP was more common in mid- and late-life, with the highest incidence in the fifth decade of life (29). We reported that hypertension increased the risk of developing BP by 50%, which is in accord with a previous case–control study that enrolled 201 patients with BP aged over 40 years or older and reported a 4.5-fold increase in the odds of BP being associated with hypertension (30).

The underlying mechanism by which an HPV infection increases the risk of developing BP remains unclear. However, an HPV infection was associated with autoimmune diseases, such as SLE, RA, and psoriasis (31, 32). A previous study showed that an HPV infection triggered SLE via molecular mimicry (22) because the HPV protein has a high and widespread similarity to several human proteins (21, 33). Similarly, molecular mimicry plays an important role in the immunopathogenic mechanisms of Guillain–Barré syndrome and related diseases (14). An antecedent infection evokes an immune response, which, in turn, cross-reacts with peripheral nerve components due to the sharing of cross-reactive epitopes. The most common identified triggering agents are Campylobacter jejuni, Mycoplasma pneumoniae, Cytomegalovirus, and Epstein–Barr virus. All of these pathogens have carbohydrate sequences in common with peripheral nerve tissues (34). Therefore, molecular mimicry might explain how an HPV infection increases the risk of BP based on the high similarity between Guillain–Barré syndrome and BP. There is currently a lack of evidence for HPV protein sequences shared with human proteins which are associated with the peripheral nervous system. Myelin basic protein, an important antigen for T cells and CNS demyelinating disease, was reported to have a region with a structural similarity to human papillomavirus peptides (35). Furthermore, inflammatory cerebrospinal fluid findings are present in BP, suggesting that it is a generalized central nervous affliction rather than a disease of the peripheral nervous system (36).

Several limitations should be considered when interpreting our findings. First, the epidemiological evaluation of a true HPV infection is challenging because many cases are not clinically recognized. Using retrospective ICD-9-based methods to select study groups of HPV infection may have a selection bias because we can only identify symptomatic patients who required medical attention. A large proportion of HPV infections may be asymptomatic; therefore, individuals who were “clinically unrecognized” could be ignored and may have been included in the comparison group. However, if HPV is associated causally with subsequent BP, misclassification would bias the estimated adjusted HRs toward the null. Therefore, for a better ascertainment of HPV exposure, we excluded patients who were diagnosed with HPV without associated treatment procedures for warts. The purpose of accruing only those patients who received associated treatment into the HPV group was to validate the diagnosis of the physician and to avoid a misclassification bias. However, some patients with an HPV infection with low-grade cervical lesions not requiring treatment might have been enrolled into the control group; similarly, such misclassification may bias our results to null, and thus our findings are acceptable. Second, the ICD-9 codes range from cutaneous benign infections to cervical and anal lesions. Grouping conditions included mucosotropic (Alpha-papillomavirus) non-carcinogenic HPVs (genital warts) to those caused by carcinogenic HPVs as well as to those in a different genus (Beta- or Gamma-papillomavirus) (plantar warts), and this is an inherent limitation of NHIRD because the ICD codes are broad but cannot provide detailed information for biological subgroup analysis. Similarly, methods using ICD-9 codes to identify cases lack information on HPV infection disease severity, which might also lead to an estimated bias because the baseline severity also affects the outcome. However, the methods applied in this study could be a proxy for HPV infection. This is because patients who request medical aid are symptomatic related to pathogenic HPV infections (either cutaneous type or mucosal type or both). These pathogenic effects of HPV infections should be prevented or treated. To overcome these limitations, further studies linking NHIRD to other disease-specified databases are needed. Third, this was a mono-country evaluation, and therefore our findings may not be applicable to non-Asian ethnic groups. Considering possible ethnic and geographical differences in the incidence and serotypes of HPV, further studies should be conducted in other ethnic groups. Fourth, HPV vaccination may also alter the immune system of the host. There are no data regarding HPV vaccination in NHIRD because HPV vaccination is a self-paid service, not covered by the NHI, and is therefore not included in the NHIRD.

In conclusion, this population-based cohort study demonstrated a higher risk of BP in patients with a previous symptomatic HPV infection, and the association was consistent in three different models. The risk was more prominent in the elderly, those with hypertension, or those with COPD. Further studies are required to clarify the underlying biological mechanisms of these associations. HPV vaccination could be another important issue for preventing BP in patients with relevant demographic characteristics and co-morbidities.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

Data are available from the NHIRD published by Taiwan National Health Insurance (NHI) Bureau. Due to legal restrictions imposed by the government of Taiwan in relation to the “Personal Information Protection Act,” data cannot be made publicly available. Requests for data can be sent as a formal proposal to the NHIRD (http://nhird.nhri.org.tw).

Ethics statement

The studies involving human participants were reviewed and approved by The Institutional Review Board of the China Medical University in Taiwan approved this study [CMUH-104-REC2-115-(AR4)]. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

Author contributions

Acquisition of data and statistical analysis were conducted by H-TY. Study conception and design were carried out by K-YL, M-CC, JW, Y-MH, RC, and H-TY. Analysis and interpretation of data were performed by K-YL, M-CC, Y-MH, RC, JW, and H-TY. K-YL, M-CC, and Y-MH contributed to writing (original draft preparation). Y-MH, M-CC, and JW contributed to writing (review and editing). All the authors were involved in drafting the article or revising it and approved the final version to be published.

Acknowledgments

The authors express their appreciation to the Department of Medical Education and Research and the Research Center of Medical Informatics in Kaohsiung Municipal United Hospital for the comments.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2021.616873/full#supplementary-material

Supplementary Figure 1

Overview of the study.

Supplementary Figure 2

The receiver operating characteristic curves of the main regression model.

Supplementary Table 1

The baseline characteristics in sensitivity analysis 1.

Supplementary Table 2

The risk of Bell's palsy in sensitivity analysis 1.

Supplementary Table 3

The baseline characteristics in sensitivity analysis 2.

Supplementary Table 4

The risk of Bell's palsy in sensitivity analysis 2.

References

  • 1.

    BaughRFBasuraGJIshiiLESchwartzSRDrumhellerCMBurkholderRet al. Clinical practice guideline: Bell's palsy. Otolaryngol Head Neck Surg. (2013) 149:S127. 10.1177/0194599813505967

  • 2.

    DeJD-SPrim-EspadaMFernández-GarcíaF. The epidemiology of Bell's palsy. Rev Neurol. (2005) 41:28790. 10.33588/rn.4105.2004593

  • 3.

    EvistonTJCroxsonGRKennedyPGHadlockTKrishnanAV. Bell's palsy: aetiology, clinical features and multidisciplinary care. J Neurol Neurosurg Psychiatry. (2015) 86:135661. 10.1136/jnnp-2014-309563

  • 4.

    Hilsinger JrRLAdourKKDotyH. Idiopathic facial paralysis, pregnancy, and the menstrual cycle. Ann Otol Rhinol Laryngol. (1975) 84:43342. 10.1177/000348947508400402

  • 5.

    PeitersenE. The natural history of Bell's palsy. Am J Otol. (1982) 4:10711.

  • 6.

    WalkerDHallamMNi MhurchadhaSMcCabePNdukaC. The psychosocial impact of facial palsy: our experience in one hundred and twenty six patients. Clin Otolaryngol. (2012) 37:4747. 10.1111/coa.12026

  • 7.

    ZhangWXuLLuoTWuFZhaoBLiX. The etiology of Bell's palsy: a review. J Neurol. (2019) 267:1896905. 10.1007/s00415-019-09282-4

  • 8.

    GrecoAGalloAFusconiMMarinelliCMacriGDe VincentiisM. Bell's palsy and autoimmunity. Autoimmun Rev. (2012) 12:3238. 10.1016/j.autrev.2012.05.008

  • 9.

    AvielAOstfeldEMarshakGBursteinRBentwichZ. Peripheral blood T and B lymphocyte subpopulations in Bell's palsy. Ann Otol Rhinol Laryngol. (1983) 92:18791. 10.1177/000348948309200218

  • 10.

    JonssonLLarssonAThomanderL. Immune complexes and complement components in Bell's palsy. ORL. (1987) 49:294301. 10.1159/000275955

  • 11.

    AbramskyOWebbCTeitelbaumDArnonR. Cellular immune response to peripheral nerve basic protein in idiopathic facial paralysis (Bell's palsy). J Neurol Sci. (1975) 26:1320. 10.1016/0022-510X(75)90109-4

  • 12.

    NylandHNaessA. Lymphocyte subpopulations in blood and cerebrospinal fluid from patients with acute Guillain-Barré syndrome. Eur Neurol. (1978) 17:24752. 10.1159/000114953

  • 13.

    VedelerCANylandHMatreR. Antibodies to peripheral nerve tissue in sera from patients with acute Guillain-Barré syndrome demonstrated by a mixed haemagglutination technique. J Neuroimmunol. (1982) 2:20914. 10.1016/0165-5728(82)90055-8

  • 14.

    YukiNHartungH-P. Guillain–Barré syndrome. N Engl J Med. (2012) 366:2294304. 10.1056/NEJMra1114525

  • 15.

    PeitersenE. Bell's palsy: the spontaneous course of 2,500 peripheral facial nerve palsies of different etiologies. Acta Oto Laryngol. (2002) 122:430. 10.1080/000164802760370736

  • 16.

    TraavikTStørvoldGSundsfjordALundSMairIW. Peripheral facial palsy and coincidental cytomegalovirus infection or reactivation. Scand J Infect Dis. (1983) 15:2338. 10.3109/inf.1983.15.issue-3.02

  • 17.

    BrownMThompsonAGohBForsterGSwashM. Bell's palsy and HIV infection. J Neurol Neurosurg Psychiatry. (1988) 51:4256. 10.1136/jnnp.51.3.425

  • 18.

    MuñozNBoschFXDeSanjosé SHerreroRCastellsaguéXShahKVet al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. (2003) 348:51827. 10.1056/NEJMoa021641

  • 19.

    BouvardVBaanRStraifKGrosseYSecretanBEl GhissassiFet al. A review of human carcinogens–Part B: biological agents. Lancet Oncol. (2009) 10:321. 10.1016/S1470-2045(09)70096-8

  • 20.

    GeorgievaSIordanovVSergievaS. Nature of cervical cancer and other HPV-associated cancers. J BU ON. (2009) 14:3918.

  • 21.

    SegalYDahanSCalabròMKanducDShoenfeldY. HPV and systemic lupus erythematosus: a mosaic of potential crossreactions. Immunol Res. (2017) 65:56471. 10.1007/s12026-016-8890-y

  • 22.

    ShiLHuangJLiuYChiouJWuRWeiJC. Risk of systemic lupus erythematosus in patients with human papillomavirus infection: a population-based retrospective cohort study. Lupus. (2018) 27:227983. 10.1177/0961203318809179

  • 23.

    LiuY-TTsouH-KChiouJ-YWangY-HChouM-CWeiJC-C. Association of human papillomavirus infection with risk for rheumatoid arthritis: a population-based cohort study. Ann Rheum Dis. (2019) 78:17346. 10.1136/annrheumdis-2019-215931

  • 24.

    LiSYuYYueYZhangZSuK. Microbial infection and rheumatoid arthritis. J Clin Cell Immunol. (2013) 4:174. 10.4172/2155-9899.1000174

  • 25.

    DoriaACanovaMTononMZenMRampuddaEBassiNet al. Infections as triggers and complications of systemic lupus erythematosus. Autoimmun Rev. (2008) 8:248. 10.1016/j.autrev.2008.07.019

  • 26.

    WangIKChenHJChengYKWuYYLinSYChouCYet al. Subdural hematoma in diabetic patients. Eur J Neurol. (2015) 22:99105. 10.1111/ene.12538

  • 27.

    LeeCDCarnahanRMMcPheetersML. A systematic review of validated methods for identifying Bell's palsy using administrative or claims data. Vaccine. (2013) 31 (Suppl. 10):K711. 10.1016/j.vaccine.2013.04.040

  • 28.

    MoralesDRDonnanPTDalyFStaaTVSullivanFM. Impact of clinical trial findings on Bell's palsy management in general practice in the UK 2001-2012: interrupted time series regression analysis. BMJ Open. (2013) 3:e003121. 10.1136/bmjopen-2013-003121

  • 29.

    ZhaoHZhangXZhuJWangX-hLiS-t. Bell's palsy: clinical analysis of 372 cases and review of related literature. Eur Neurol. (2017) 77:16872. 10.1159/000455073

  • 30.

    Savadi-OskoueiDAbediASadeghi-BazarganiH. Independent role of hypertension in Bell's palsy: a case-control study. Eur Neurol. (2008) 60:2537. 10.1159/000151701

  • 31.

    ChenMLKaoWMHuangJYHungYMWeiJC. Human papillomavirus infection associated with increased risk of new-onset psoriasis: a nationwide population-based cohort study. Int J Epidemiol. (2020) 49:78697. 10.1093/ije/dyaa027

  • 32.

    WeiC-YLinJ-YWangY-THuangJ-YWeiJC-CChiouJ-Y. Risk of ankylosing spondylitis following human papillomavirus infection: a nationwide, population-based, cohort study. J Autoimmun. (2020) 113:102482. 10.1016/j.jaut.2020.102482

  • 33.

    KanducDShoenfeldY. Human papillomavirus epitope mimicry and autoimmunity: the molecular truth of peptide sharing. Pathobiology. (2019) 86:28595. 10.1159/000502889

  • 34.

    HaddenRKarchHHartungH-PZielasekJWeissbrichBSchubertJet al. Preceding infections, immune factors, and outcome in Guillain–Barré syndrome. Neurology. (2001) 56:75865. 10.1212/WNL.56.6.758

  • 35.

    Gonzalez-GronowMPizzoSV. Relevance of catalytic autoantibodies to myelin basic protein (MBP) in autoimmune disorders. J Neurol Neuromed. (2018) 3:758 (2018).

  • 36.

    RobergMErnerudhJForsbergPFridellEFrydénAHydenDet al. Acute peripheral facial palsy: CSF findings and etiology. Acta Neurol Scand. (1991) 83:5560. 10.1111/j.1600-0404.1991.tb03959.x

Summary

Keywords

Bell's palsy, human papillomavirus, cohort study, infection, HPV

Citation

Li K-Y, Chou M-C, Chang R, Yip H-T, Hung Y-M and Wei JC-C (2021) The Potential Role of Human Papillomavirus Infection in Bell's Palsy: A Hypothesis-Generating Study Based on a Nationwide Cohort. Front. Med. 8:616873. doi: 10.3389/fmed.2021.616873

Received

13 October 2020

Accepted

26 July 2021

Published

01 September 2021

Volume

8 - 2021

Edited by

Ana Afonso, University of São Paulo, Brazil

Reviewed by

Jun-Jun Yeh, Ditmanson Medical Foundation Chia-Yi Christian Hospital, Taiwan; Chun-Yu Lin, National Cheng Kung University Hospital, Taiwan

Updates

Copyright

*Correspondence: Yao-Min Hung Renin Chang James Cheng-Chung Wei

This article was submitted to Infectious Diseases – Surveillance, Prevention and Treatment, a section of the journal Frontiers in Medicine

†These authors share first authorship

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics