Next Article in Journal
Improving the Safety of DIEP Flap Transplantation: Detailed Perforator Anatomy Study Using Preoperative CTA
Previous Article in Journal
Development and Validation of an Artificial Intelligence Electrocardiogram Recommendation System in the Emergency Department
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Human Papillomavirus Infection and the Risk of Erectile Dysfunction: A Nationwide Population-Based Matched Cohort Study

1
Department of Anesthesiology, Kaohsiung Chang Gung Memorial Hospital, Kaohsiung 833401, Taiwan
2
Center for Global Health, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA
3
Department of Epidemiology, Harvard T.H. Chan School of Public Health, Boston, MA 02115-5810, USA
4
Graduate Institute of Biomedical Electronics and Bioinformatics, National Taiwan University, Taipei 11114, Taiwan
5
School of Medicine, Chung Shan Medical University, Taichung 40402, Taiwan
6
Institute of Medicine, Chung Shan Medical University, Taichung 40402, Taiwan
7
Graduate Institute of Integrated Medicine, China Medical University, Taichung 40402, Taiwan
8
Division of Allergy, Immunology and Rheumatology, Chung Shan Medical University, Taichung 40402, Taiwan
9
Department of Management office for Health Data, China Medical University Hospital, Taichung 40402, Taiwan
10
College of Medicine, China Medical University, Taichung 40402, Taiwan
11
Institute of Public Health (Biostatistics), National Yangming University, Taipei 112304, Taiwan
12
Department of Recreation and Sports Management, Tajen University, Pingtung County 907101, Taiwan
13
Department of Physical Medicine and Rehabilitation, Kaohsiung Veterans General Hospital, Pingtung Branch, Pingtung County 907101, Taiwan
14
Department of Internal Medicine, Kaohsiung Municipal United Hospital, Kaohsiung 813414, Taiwan
15
Shu-Zen Junior College of Medicine and Management, Kaohsiung 813414, Taiwan
16
Department of Orthopedics, Taichung Veterans General Hospital, Taichung 40402, Taiwan
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Pers. Med. 2022, 12(5), 699; https://doi.org/10.3390/jpm12050699
Submission received: 20 February 2022 / Revised: 30 March 2022 / Accepted: 18 April 2022 / Published: 27 April 2022

Abstract

:
Background: Male patients with genital warts are known for higher rates of sexual dysfunction. This study was conducted to investigate whether human papillomaviruses (HPV) infection is associated with an increased risk of erectile dysfunction (ED). Methods: Patients aged over 18 with HPV infection (n = 13,296) and propensity score-matched controls (n = 53,184) were recruited from the Longitudinal Health Insurance Database (LHID). The primary endpoint was the diagnosis of ED. Chi-square tests were used to analyze the distribution of demographic characteristics. The Cox proportional hazards regression was used to estimate the hazard ratios (HRs) and 95% confidence intervals (CIs) for the development of ED in both groups, after adjusting for sex, age, relevant comorbidities, co-medication, and surgery. Results: ED developed in 181 patients of the study group. The incidence density of ED was 2.53 per 1000 person-years for the HPV group and 1.51 per 1000 person-years for the non-HPV group, with an adjusted HR (95% CI) of 1.63 (1.37–1.94). In stratification analysis, adjusted HR of diabetes-, chronic obstructive pulmonary disease (COPD-), and stroke-subgroup were 2.39, 2.51, and 4.82, with significant p values for interaction, respectively. Sensitivity analysis yields consistent findings. Conclusions: The patients with HPV infection had a higher risk of subsequent ED in comparison to the non-HPV controls. The mechanism behind such association and its possible role in ED prevention deserves further study in the future.

1. Introduction

Human papillomaviruses (HPVs) are small, double-stranded DNA viruses that infect the cutaneous and mucosa epithelium and are the most common sexually transmitted disease in women [1,2]. There are more than 200 recognized types of HPV [3], of which 40 types infect the genital area and a dozen or so types of HPV are known to be high risk and lead to cervical malignancy [4,5,6]. With the administration of HPV vaccines, a significant decrease in the prevalence of HPV infections, cervical lesions, and cervical malignancies in young females has been observed [7,8,9]. Despite the enhanced efficacy, proven safety, and recommendation by the World Health Organization [10], the public’s acceptance of HPV vaccines could still be improved, especially in men [11]. In the U.S. population-based study, the prevalence of genital HPV infection was estimated at 41.3% to 49.3% in males but HPV vaccination coverage was only approximated 7.8% to 14.6% among vaccine-eligible men [12]. Thus, the adverse impact of HPV infections on public health remains. Apart from malignancies, complications of HPV infections include genital warts, bladder cancer, and male infertility due to impaired sperm function [13,14]. However, compared to the extensive number of studies on HPV infections in women, there is a lack of knowledge on HPV infections in male patients [15].
Evidence suggests that males with genital warts have higher rates of sexual dysfunction [16]. Among sexual dysfunction, erectile dysfunction (ED), defined as the inability to maintain or achieve an erection necessary for sexual intercourse, afflicts an estimated 10–20 million men in the United States alone and a survey in Taiwan showed that the prevalence of ED was 27% among all respondents aged not less than 30 years [17,18,19]. The pathogenesis of ED depends on the vascular integrity of the erectile tissue and the HPV infections may correlate with vascular complications [20]. For instance, a case of isolating necrotizing granulomatous vasculitis of the genital area associated with HPV infection has been reported [21]. There are several potential mechanisms underlying the relationship between HPV and ED. First of all, ED is a neurovascular procession that depends on the vascular health of erectile tissue and the nervous system [20]. HPV infection may induce inflammation and cause atherosclerosis progress, which has been shown to seriously affect the vascular health of the erectile tissue, thus contributing to ED pathogenesis [22,23,24,25,26]. Common risk factors for atherosclerosis have been frequently found in patients with ED, with the extent of ED being associated with the number and severity of CVDs [22,23,24,25,26]. For instance, ED has been reported to be associated with elevated high-sensitivity C-reactive protein values and to be correlated with the magnitude of flow-mediated dilation of the brachial artery [23], with the latter being an indicator of endothelial dysfunction [23]. In addition, it has been reported that ED is correlated with coronary calcification and atherosclerotic plaque [23,26] and usually precedes the development of clinically evident CVDs [22,23,26]. These findings indicate that ED may share common pathophysiological mechanisms with CVDs, particularly coronary artery disease (CAD). For instance, both ED and CAD have endothelial dysfunction as part of their pathophysiology [27,28,29], which is evidenced by studies indicating that ED is associated with increased circulating levels of systemic inflammatory [30,31,32,33,34] and endothelial–prothrombotic markers/mediators [35]. As the diameters of the cavernosal arteries are smaller than large coronary arteries, they may be more sensitive to atherosclerosis [36]. Accordingly, HPV-associated atherosclerosis and affected vascular function may be involved in ED development.
Previous studies propose that inflammation caused by viral [37,38] or bacterial [39,40,41,42,43,44] pathogens plays an important role in atherosclerosis pathogenesis. Specifically, previous studies have shown that herpes zoster virus, Chlamydia pneumoniae, cytomegalovirus, and hepatitis A virus infection can increase the risk of CAD, myocardial infarction, and even cardiovascular death [45,46,47,48,49]. For instance, previous studies using data from the National Health and Nutrition Examination Survey (NHANES) demonstrated that HPV infection is associated with CVD among women [50]. In addition, infectious agents, including viruses, bacteria, and parasites, can trigger autoimmunity and are associated with secondary vasculitis [51], which may lead to cases such as HPV-associated granulomatous vasculitis of the genital area [21]. The etiologies of ED have also been reported to include chronic systemic vasculitis and severe venous leak caused by venous thrombosis, leading to the fact that HPV infection is related to vascular dysfunction [52,53,54].
Given the fact that genital warts negatively affect male sexual function [16], the manner in which HPV infections fit into our knowledge of ED etiologies remains undetermined. Therefore, we conducted this nationwide population-based cohort study using a longitudinal database in Taiwan to provide epidemiological evidence about the impact of HPV infection on ED development.

2. Materials and Methods

2.1. Data Source

This study was constructed using longitudinal data from the National Health Insurance Research Database (NHIRD) in Taiwan. The NHIRD contains registry and original claimed data, such as the diagnoses coding, out-patient visits, hospitalization records, medication, and personal information of over 99% of Taiwan’s population. Quarterly expert reviews on randomly sampled claimed data, with a sampling rate of 1/50–100, were performed by the Bureau of National Health Insurance (BNHI). The BNHI also randomly reviews the medical charts of 1/100 ambulatory and 1/20 inpatient claims. This study was approved by the Institutional Review Board of the China Medical University Hospital Research Ethics Committee (CMUH104-REC2-115(AR-4)).

2.2. Study Population

Participants with medical records of HPV infection (ICD-9 codes 079.4, 078.10, 078.11, 078.12, 078.19, 795.05, 795.09, 795.15, 795.19, 796.75, and 796.79) between 2000 and 2012 constituted the HPV group (Figure 1). The index date was defined as the date of diagnosis and coding for HPV infection. The control group was selected from the same dataset and was matched by age, sex, index date, co-morbidities, and co-medication at a ratio of 1:4. Individuals diagnosed with ED before the index date or younger than 18 years old were excluded. To ensure the validity of the diagnoses of HPV, patients diagnosed with an HPV infection who never received HPV-associated treatment within three months after the index date were excluded. These treatments included excision, electrocauterization, carbon dioxide (CO2) laser operation, or chemosurgery for condyloma, as well as simple or complicated electrocauterization, liquid nitrogen cryosurgery, or cryotherapy using CO2 freezing or liquid nitrogen. All participants were followed up until the first medical record of ED or censored on the date of death, or the end of the study (31 December 2013). Finally, 13,296 subjects were included in the HPV group and 53,184 non-HPV subjects served as the control group.

2.3. Definitions of the Outcomes and Covariates

The primary endpoint of this study was set to be the diagnosis of ED (ICD-9-CM code 607.84 and 302.72) that recorded at least one inpatient or three outpatients [55], withdrawal, or end of follow-up. To ensure the reliability and accuracy of the diagnoses, the diagnoses had to be assigned by urologists for three or more ambulatory visits or at least one admission within a year. In a previous study, ICD-9 codes for ED have been validated, and corresponded to an International Index of Erectile Function score of less than 8 to define severe ED; using ICD-9 codes to enroll cases, the sensitivity was 80%, and the positive predictive value was 88.9% [56].
To eliminate potential bias, we adjusted for demographic variables, pre-existing co-morbidities and baseline co-medication, including hypertension (ICD-9-CM codes 401–405), diabetes mellitus (ICD-9-CM code 250), hyperlipidemia (ICD-9-CM code 272), stroke (ICD-9-CM code 430–438), coronary artery disease (ICD-9-CM code 410–414), chronic kidney disease (CKD) (ICD-9-CM code 585), chronic obstructive pulmonary disease (COPD) (ICD-9-CM codes 491,492, and 496), alcohol-related illness (ICD-9-CM codes 291, 303, 305, 571.0, 571.1, 571.2, 571.3, 790.3, A215, and V11.3), human immunodeficiency virus (HIV) infection (ICD-9-CM codes 042–044, 795.8, and V08), autoimmune disease (AID) including systemic lupus erythematosus (SLE) (ICD-9-CM code 710.0), rheumatoid arthritis (ICD-9-CM code 714.0), Sjögren’s syndrome (ICD-9-CM code 710.2), and multiple sclerosis (ICD-9-CM code 340). Other metrics that can affect ED were also adjusted in the analysis, including spinal cord injury pelvic fracture, burn of genitalia, foreign body entering through an orifice, injury to nerves and spinal cord, and poisoning. Information on comorbidities was obtained by tracing back 2 years before the index date about records of inpatients and ambulatory medical care in the database.
The medication confounders analyzed in this study consist of alpha-blockers, beta-blockers, diuretics, and calcium channel blockers (CCBs). Drug use was defined as the use of a drug for more than three prescriptions (with each prescription for 28 days) from outpatient visits prior to the index date.

2.4. Negative Control Exposure

A negative control exposure was designated to examine potential external confounders [57,58]. One can conduct an analysis with an alternative exposure that is not expected to be associated with the outcome of interest. In this study, Nontyphoidal Salmonella infection (NTS) was selected as the negative control exposure for its potential association with ischemic events [59,60]. However, there is no association between NTS and ED from a comprehensive review of the literature. We hypothesized that any association between NTS infection and ED implies possible unmeasured confounders. Accordingly, patients diagnosed with NTS infection (ICD-9-CM code 003.5) in at least three outpatient claims or one hospitalization from 2001 through 2012 and propensity score (PS)-matched non-salmonellosis controls were compared.

2.5. Statistical Analysis

Chi-squared tests were applied to compare the distribution of age, sex, and baseline co-morbidities between the HPV and non-HPV groups. The age mean was compared using Student’s t-test. The incidence density of ED per 1000 person-years was calculated in both groups. Multivariable Cox proportional hazards regression models were utilized to estimate the crude hazard ratios (cHRs), adjusted HRs (aHRs), and 95% confidence intervals (CIs) of ED for the HPV group versus the non-HPV group. The covariates adopted in the multivariate models encompassed sex, age, co-morbidities, and co-medications listed in Table 1. The Kaplan–Meier method was used to plot the cumulative incidence curves of ED for the HPV and non-HPV groups. Differences between the two groups were evaluated by the log-rank test. The incidence rates of ED were estimated by dividing the number of events by follow-up person-years for both groups. Two statistical models were fitted to evaluate the effect of HPV on the risk of ED (Table 2). In model 1, we examined the temporal relationship between HPV exposure and the risk of developing ED for both groups, adjusted for age, co-morbidities, and co-medications at baseline. Model 2 was adjusted for all covariates in model 1, and further adjusted for relevant operations associated with subsequent ED to minimize confounding, (e.g., foreign body entering through an orifice, low anterior resection, and abdominoperineal resection) [61,62,63]. Subgroup analyses were conducted to reveal the potential interaction effect of sex, age, and follow-up between HPV infection and subsequent ED development. All data analyses were performed with SAS (version 9.4; SAS Institute, Inc., Carey, NC, USA). The statistical significance level was set at a two-tailed p-value of <0.05.

2.6. Sensitivity Analyses

We performed a sensitivity analysis excluding patients with cancer and stroke at the index date to examine whether our findings remain constant under different assumptions.

3. Results

We enrolled 13,296 patients newly diagnosed with an HPV infection from 2000 to 2012 and 53,184 non-HPV controls in this study. The baseline characteristics of patients with and without an HPV infection were not significantly different, as all standard mean differences between the groups were <0.1 (Table 1). The mean ages of participants in the HPV and non-HPV groups were 39.4 (16.3) and 39.0 (15.7) years after PS matching, respectively.
In Table 2, model 1 and model 2 both show that patients with a history of HPV infection had a higher risk of developing ED (model 1, aHR, 1.64; 95% CI 1.38–1.95 and model 2, aHR, 1.63; 95% CI 1.37–1.94; p < 0.001). Table 2 also shows a sensitivity analysis, in which participants with severe comorbidities (cancer, stroke, and autoimmune diseases) that can cause (not just a risk factor for) ED were excluded. The results show consistent findings (aHR, 1.59; 95% CI 1.30–1.95; p < 0.001; details in Supplement Tables S1 and S2).
Negative control analysis showed that there was no association between NTS and subsequent ED with aHR, 0.69; 95% CI 0.31–1.53 (detail in Supplement Tables S3 and S4). This measurement revealed that there is minimal uncontrolled confounding from the point of view.
The finding that HPV infection is associated with higher ED risks was shown in the survival analysis, in which the cumulative incidence of ED for the HPV cohort was higher than that for the non-HPV cohort (log-rank test, p < 0.001) (Figure 1).
Table 3 shows that ED occurred significantly at around 40 years in this study and this finding was as consistent as the epidemiologic report [64]. Compared with individuals aged less than 30 years (this age group as reference), those aged 31–40 years had a higher risk of reporting ED (aHR, 1.89; 95% CI 1.40–2.56); for those aged 41–50 years, the aHR was 3.17 (95% CI 2.39–4.20). In Table 3, higher risks of ED were also identified among patients with hyperlipidemia (aHR, 1.63; 95% CI 1.34–1.97), alcohol-related illness (aHR, 1.43; 95% CI 1.01–2.05), AID (aHR, 1.66; 95% CI 1.09–2.53) and use of alpha- and beta-blockers (aHR, 4.07; 95% CI 3.33–4.97 and 1.30; 95% CI 1.03–1.64, respectively).
Table 4 shows the Cox proportional hazards regression models for comparing HPV and non-HPV patients in subgroup analyses. In the age subgroup of 31–40 years, compared with the age-matched non-HPV participants, those with HPV infection had a higher risk of developing ED (aHR, 1.65; 95% CI 1.05–2.60); in the age subgroup of 41–50 years, compared with the age-matched non-HPV participants, those with HPV infection had a higher risk of having ED (aHR, 1.63; 95% CI 1.17–2.27). The p-value for the interaction of age was not significant (p = 0.4). In comorbidities-subgroup analysis, among participants with diabetes, there was a positive association between HPV infection with ED (aHR, 2.39; 95% CI 1.63–3.51; p for interaction 0.05). Such positive association was also found among patients with stroke (aHR, 4.82; 95% CI 2.34–9.92), and COPD (aHR, 2.51; 95% CI 1.69–3.74) with p values for interaction 0.001 and 0.03, respectively.
In Table 5, subgroup analysis stratified by follow-up duration demonstrated that the risk of ED among HPV-infected patients was significantly higher in general after HPV infection. The risk of ED onset at <1, 1–5, and >5 years after HPV infection was 1.64 (95% CI = 1.38, 1.95), 1.48 (95% CI = 1.22, 1.80), and 1.7 (95% CI = 1.24, 2.34), respectively.

4. Discussion

In this large-scale population-based cohort study with up to 14 years of follow-up, we reported that individuals with HPV infection had a higher risk of developing ED. We selected symptomatic HPV patients as the study population because both inpatients and outpatients with symptoms are rarely studied. The subgroup analyses indicated that HPV-infected patients had a significantly higher risk of ED if they had a stroke or COPD. That is, COPD and stroke acted as modifiers for the relationship between HPV infection and ED. COPD is a consequence of tobacco use which is a well-known risk for ED. A time-to-event analysis indicated a positive association between HPV infection and ED in all the follow-up periods, which may be due to the characteristic of recurrent HPV infection [65]. However, the effect of persistent or recurrent HPV infection on the risk of ED has not been explored in this study. These findings are consistent in the sensitivity analysis and negative control analysis.
Previous studies on HPV infection have reported that women of reproductive age with positive HPV tests for cervical screening have poorer sexual function [66]. On the contrary, HPV infections among males have been investigated less, even though HPV DNA has been proven to lie not only in the perianal region and external genitalia, including the penis foreskin, scrotum, and glans penis, but also in the urethra, ductus deferens, epididymis, or testis [67]. For instance, HPV may infect the semen, contributing to male infertility [68]. This is due to decreased sperm progressive motility and the normal morphology rate caused by HPV infection [68,69]. Moreover, higher risks of anxiety and depression have been observed in patients with genital warts [16]. Having HPV causes fear, feeling of stigmatization, and worry about disclosing the disease to partners, family, or friends [16]. Furthermore, the psychological impact of HPV infection includes a negative association between genital warts and Arizona Sexual Experience Scale (ASEX) score, with lower ASEX scores being associated with a decline in the Beck Anxiety Inventory (BAI) and Beck Depression Inventory (BDI) indices [16]. Since the common etiologies of psychogenic ED include performance anxiety, a strained relationship, lack of sexual arousability, and overt psychiatric disorders such as depression and schizophrenia [70], it is likely that the psychological impact of HPV infection may also lead to ED.
In this cohort study, we employed a nationwide population-based registry to provide epidemiological evidence on the risks of new-onset ED in patients with HPV infection. The effect of measurable confounders was minimized by PS matching [71,72,73,74,75]. Sensitivity and negative control analyses were conducted to manage residual confounding and to ascertain the impact of HPV infection on ED development.
However, there are several limitations to this study. First, the evaluation of HPV infection is challenging, as many infections might not be clinically recognized [76,77,78]. The selection bias resulting from a retrospective study design means that only patients with a symptomatic HPV infection were included. That is, whether patients with asymptomatic or self-resolving infection also have a higher risk of developing ED remains unknown. Second, due to inherent database limitations, this study lacked penile ultrasound data, other indicators that may affect erectile function such as the use of phosphodiesterase-5 (PDE5) inhibitors, and an erectile performance questionnaire; however, in a previous study, the ICD-9 code for ED has been validated with the International Index of Erectile Function score [56]. Third, as our study only included participants from one country of Asian origin, our findings may not be generalizable to other ethnic groups.

5. Conclusions

Findings of this cohort study suggest that ED is an overlooked long-term complication of HPV infections in males. Future studies on HPV infections in males are warranted to facilitate the diagnosis and management of HPV infections.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jpm12050699/s1, Table S1: Baseline characteristic in sensitivity analysis (excluding participants with cancer and stroke on the index date); Table S2: Sensitivity analysis (excluding participants with cancer and stroke on the index date); Table S3: Baseline characteristic in negative control analysis; Table S4: Negative control analysis (NTS infection as the negative control exposure).

Author Contributions

Conceptualization, S.-E.J., K.S.-K.M., P.-E.K., J.C.-C.W., H.-T.Y., M.-C.C., Y.-M.H. and N.-C.C.; formal analysis, S.-E.J., K.S.-K.M., P.-E.K., J.C.-C.W., H.-T.Y., M.-C.C., Y.-M.H. and N.-C.C.; data curation, H.-T.Y.; writing—original draft preparation, S.-E.J., K.S.-K.M., M.-C.C. and Y.-M.H.; writing—review and editing, P.-E.K., Y.-M.H. and M.-C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was approved by the Institutional Review Board of the China Medical University Hospital Research Ethics Committee (CMUH104-REC2-115(AR-4)).

Informed Consent Statement

The Institutional Review Board of the China Medical University Hospital Research Ethics Committee approved this study, and patient consent was waived due to anonymized data.

Data Availability Statement

The data used in this study was held by the Taiwan Ministry of Health and Welfare, and thus was not publicly available. The data can be accessed by submitting an application to the Ministry of Health and Welfare for access. (Taiwan Ministry of Health and Welfare, Address: No. 488, Sec. 6, Zhongxiao E. Rd., Nangang Dist., Taipei City 115, Taiwan (R.O.C.). Phone: +886-2-8590-6848).

Acknowledgments

The authors express their appreciation to the Department of Medical Education and Research and Research Center of Medical Informatics in Kaohsiung Veterans General Hospital, as well as the Department of Medical Education and Research and Research Center of Medical Informatics in Kaohsiung Municipal United Hospital, for their comments. Special thanks to Renin Chang for comments about study design.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gheit, T. Mucosal and Cutaneous Human Papillomavirus Infections and Cancer Biology. Front. Oncol. 2019, 9, 355. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Zur Hausen, H. Papillomavirus infections—A major cause of human cancers. Biochim. Biophys. Acta 1996, 1288, F55–F78. [Google Scholar] [CrossRef]
  3. Kocjan, B.; Bzhalava, D.; Forslund, O.; Dillner, J.; Poljak, M. Molecular methods for identification and characterization of novel papillomaviruses. Clin. Microbiol. Infect. 2015, 21, 808–816. [Google Scholar] [CrossRef] [Green Version]
  4. Muñoz, N.; Bosch, F.X.; De Sanjosé, S.; Herrero, R.; Castellsagué, X.; Shah, K.V.; Snijders, P.J.F.; Meijer, C.J.L.M. Epidemiologic Classification of Human Papillomavirus Types Associated with Cervical Cancer. N. Engl. J. Med. 2003, 348, 518–527. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Bouvard, V.; Baan, R.; Straif, K.; Grosse, Y.; Secretan, B.; El Ghissassi, F.; Benbrahim-Tallaa, L.; Guha, N.; Freeman, C.; Galichet, L.; et al. A review of human carcinogens—Part B: Biological agents. Lancet Oncol. 2009, 10, 321–322. [Google Scholar] [CrossRef]
  6. Bodily, J.; Laimins, L.A. Persistence of human papillomavirus infection: Keys to malignant progression. Trends Microbiol. 2011, 19, 33–39. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Markowitz, L.E.; Hariri, S.; Lin, C.; Dunne, E.F.; Steinau, M.; McQuillan, G.; Unger, E. Reduction in Human Papillomavirus (HPV) Prevalence among Young Women Following HPV Vaccine Introduction in the United States, National Health and Nutrition Examination Surveys, 2003–2010. J. Infect. Dis. 2013, 208, 385–393. [Google Scholar] [CrossRef] [Green Version]
  8. Hariri, S.; Mph, M.L.J.; Bennett, N.M.; Bauer, H.M.; Park, I.U.; Schafer, S.; Niccolai, L.M.; Unger, E.R.; Markowitz, L.E.; HPV-IMPACT Working Group. Population-based trends in high-grade cervical lesions in the early human papillomavirus vaccine era in the United States. Cancer 2015, 121, 2775–2781. [Google Scholar] [CrossRef] [Green Version]
  9. Guo, F.; Cofie, L.E.; Berenson, A.B. Cervical Cancer Incidence in Young U.S. Females After Human Papillomavirus Vaccine Introduction. Am. J. Prev. Med. 2018, 55, 197–204. [Google Scholar] [CrossRef]
  10. World Health Organization. Human papillomavirus vaccines: WHO position paper, May 2017-Recommendations. Vaccine 2017, 35, 5753–5755. [Google Scholar] [CrossRef]
  11. Ozawa, K.; Hineno, A.; Kinoshita, T.; Ishihara, S.; Ikeda, S.-I. Suspected Adverse Effects After Human Papillomavirus Vaccination: A Temporal Relationship Between Vaccine Administration and the Appearance of Symptoms in Japan. Drug Saf. 2017, 40, 1219–1229. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Han, J.J.; Beltran, T.H.; Song, J.W.; Klaric, J.; Choi, Y.S. Prevalence of Genital Human Papillomavirus Infection and Human Papillomavirus Vaccination Rates Among US Adult Men: National Health and Nutrition Examination Survey (NHANES) 2013–2014. JAMA Oncol. 2017, 3, 810–816. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, Y.; Jia, C.-W.; Ma, Y.-M.; Zhou, L.-Y.; Wang, S.-Y. Correlation between HPV sperm infection and male infertility. Asian J. Androl. 2013, 15, 529–532. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Foresta, C.; Noventa, M.; De Toni, L.; Gizzo, S.; Garolla, A. HPV-DNA sperm infection and infertility: From a systematic literature review to a possible clinical management proposal. Andrology 2015, 3, 163–173. [Google Scholar] [CrossRef] [PubMed]
  15. Cheng, Y.-P.; Chen, C.-W.; Sheen, Y.-S.; Tsai, T.-F. Genotype distribution of human papillomavirus in anogenital warts of male patients in Taiwan. Dermatol. Sin. 2012, 30, 85–89. [Google Scholar] [CrossRef]
  16. Kucukunal, A.; Altunay, I.K.; Mercan, S. Sexual Dysfunction in Men Suffering from Genital Warts. J. Sex. Med. 2013, 10, 1585–1591. [Google Scholar] [CrossRef]
  17. Benet, A.E.; Melman, A. The Epidemiology of Erectile Dysfunction. Urol. Clin. N. Am. 1995, 22, 699–709. [Google Scholar] [CrossRef]
  18. Feldman, H.A.; Goldstein, I.; Hatzichristou, D.G.; Krane, R.J.; McKinlay, J.B. Impotence and Its Medical and Psychosocial Correlates: Results of the Massachusetts Male Aging Study. J. Urol. 1994, 151, 54–61. [Google Scholar] [CrossRef]
  19. Hwang, T.I.; Tsai, T.-F.; Lin, Y.-C.; Chiang, H.-S.; Chang, L.S. A Survey of Erectile Dysfunction in Taiwan: Use of the Erection Hardness Score and Quality of Erection Questionnaire. J. Sex. Med. 2010, 7, 2817–2824. [Google Scholar] [CrossRef]
  20. Traish, A.M.; Goldstein, I.; Kim, N. Testosterone and Erectile Function: From Basic Research to a New Clinical Paradigm for Managing Men with Androgen Insufficiency and Erectile Dysfunction. Eur. Urol. 2007, 52, 54–70. [Google Scholar] [CrossRef] [Green Version]
  21. Winegar, A.L.; Shepherd, M.D.; Lawson, K.A.; Richards, K.M. Comparison of the claim percent gross margin earned by Texas community independent pharmacies for dual-eligible beneficiary claims before and after Medicare Part D. J. Am. Pharm. Assoc. 2009, 49, 617–622. [Google Scholar] [CrossRef] [PubMed]
  22. Montorsi, F.; Briganti, A.; Salonia, A.; Rigatti, P.; Margonato, A.; Macchi, A.; Galli, S.; Ravagnani, P.M.; Montorsi, P. Erectile dysfunction prevalence, time of onset and association with risk factors in 300 consecutive patients with acute chest pain and angiographically documented coronary artery disease. Eur. Urol. 2003, 44, 360–365. [Google Scholar] [CrossRef]
  23. Chiurlia, E.; D’Amico, R.; Ratti, C.; Granata, A.R.; Romagnoli, R.; Modena, M.G. Subclinical Coronary Artery Atherosclerosis in Patients with Erectile Dysfunction. J. Am. Coll. Cardiol. 2005, 46, 1503–1506. [Google Scholar] [CrossRef] [Green Version]
  24. Ponholzer, A.; Temml, C.; Obermayr, R.; Wehrberger, C.; Madersbacher, S. Is erectile dysfunction an indicator for increased risk of coronary heart disease and stroke? Eur. Urol. 2005, 48, 512–518. [Google Scholar] [CrossRef]
  25. Thompson, I.M. Erectile Dysfunction and Subsequent Cardiovascular Disease. JAMA 2005, 294, 2996–3002. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Vlachopoulos, C.; Rokkas, K.; Ioakeimidis, N.; Aggeli, C.; Michaelides, A.; Roussakis, G.; Fassoulakis, C.; Askitis, A.; Stefanadis, C. Prevalence of Asymptomatic Coronary Artery Disease in Men with Vasculogenic Erectile Dysfunction: A Prospective Angiographic Study. Eur. Urol. 2005, 48, 996–1003. [Google Scholar] [CrossRef]
  27. Kaiser, D.R.; Billups, K.; Mason, C.; Wetterling, R.; Lundberg, J.L.; Bank, A.J. Impaired brachial artery endothelium-dependent and -independent vasodilation in men with erectile dysfunction and no other clinical cardiovascular disease. J. Am. Coll. Cardiol. 2004, 43, 179–184. [Google Scholar] [CrossRef] [Green Version]
  28. Montorsi, P.; Montorsi, F.; Schulman, C.C. Is erectile dysfunction the “tip of the iceberg” of a systemic vascular disorder? Eur. Urol. 2003, 44, 352–354. [Google Scholar] [CrossRef]
  29. Solomon, H.; Man, J.W.; Jackson, G. Erectile dysfunction and the cardiovascular patient: Endothelial dysfunction is the common denominator. Heart 2003, 89, 251–253. [Google Scholar] [CrossRef] [Green Version]
  30. Ma, K.S.-K.; Chung, W.H.; Hsueh, Y.-J.; Chen, S.-Y.; Tokunaga, K.; Kinoshita, S.; Ma, D.H.K.; Ueta, M. Human leucocyte antigen association of patients with Stevens-Johnson syndrome/toxic epidermal necrolysis with severe ocular complications in Han Chinese. Br. J. Ophthalmol. 2021, 106, 610–615. [Google Scholar] [CrossRef]
  31. Ma, K.S.-K.; Saeed, H.N.; Chodosh, J.; Wang, C.-W.; Chung, Y.-C.; Wei, L.-C.; Kuo, M.-T.; Liang, C.-M.; Chang, J.W.-C.; Chung, W.-H.; et al. Ocular manifestations of anti-neoplastic immune checkpoint inhibitor-associated Stevens-Johnson syndrome/toxic epidermal necrolysis in cancer patients. Ocul. Surf. 2021, 22, 47–50. [Google Scholar] [CrossRef] [PubMed]
  32. Huang, J.-W.; Kuo, C.-L.; Wang, L.-T.; Ma, K.S.-K.; Huang, W.-Y.; Liu, F.-C.; Yang, K.D.; Yang, B.-H. Case Report: In Situ Vaccination by Autologous CD16+ Dendritic Cells and Anti-PD-L 1 Antibody Synergized With Radiotherapy To Boost T Cells-Mediated Antitumor Efficacy In A Psoriatic Patient With Cutaneous Squamous Cell Carcinoma. Front. Immunol. 2021, 12, 752563. [Google Scholar] [CrossRef] [PubMed]
  33. Ma, K.S.; Wei, J.C.; Chung, W.H. Correspondence to ‘Hypersensitivity reactions with allopurinol and febuxostat: A study using the Medicare claims data’. Ann. Rheum. Dis. 2020. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, L.T.; Ma, K.S. Correspondence to ‘Normal human enthesis harbours conventional CD4+ and CD8+ T cells with regulatory features and inducible IL-17A and TNF expression’. Ann. Rheum. Dis. 2020. [Google Scholar] [CrossRef] [PubMed]
  35. Vlachopoulos, C.; Aznaouridis, K.; Ioakeimidis, N.; Rokkas, K.; Vasiliadou, C.; Alexopoulos, N.; Stefanadi, E.; Askitis, A.; Stefanadis, C. Unfavourable endothelial and inflammatory state in erectile dysfunction patients with or without coronary artery disease. Eur. Heart. J. 2006, 27, 2640–2648. [Google Scholar] [CrossRef] [Green Version]
  36. Billups, K.L.; Bank, A.J.; Padma-Nathan, H.; Katz, S.; Williams, R. Erectile dysfunction is a marker for cardiovascular disease: Results of the minority health institute expert advisory panel. J. Sex. Med. 2005, 2, 40–50, discussion 50–52. [Google Scholar] [CrossRef]
  37. Ma, K.S.; Lee, C.C.; Liu, K.J.; Wei, J.C.; Lee, Y.T.; Wang, L.T. Safety and Seroconversion of Immunotherapies against SARS-CoV-2 Infection: A Systematic Review and Meta-Analysis of Clinical Trials. Pathogens 2021, 10, 1537. [Google Scholar] [CrossRef]
  38. Chiang, C.-H.; Tang, P.-U.; Lee, G.H.; Chiang, T.-H.; Ma, K.S.-K.; Fang, C.-T. Prevalence of Nontuberculous Mycobacterium Infections versus Tuberculosis among Autopsied HIV Patients in Sub-Saharan Africa: A Systematic Review and Meta-Analysis. Am. J. Trop. Med. Hyg. 2020, 104, 628–633. [Google Scholar] [CrossRef]
  39. Ma, K.S.-K.; Chiang, C.-H.; Chen, Y.-W.; Wang, L.-T. Correspondence to ‘Bacterial citrullinated epitopes generated by Porphyromonas gingivalis infection-a missing link for ACPA production’. Ann. Rheum. Dis. 2021, 79, 1194–1202. [Google Scholar] [CrossRef]
  40. Wu, M.-C.; Ma, K.S.-K.; Wang, Y.-H.; Wei, J.C.-C. Impact of tonsillectomy on irritable bowel syndrome: A nationwide population-based cohort study. PLoS ONE 2020, 15, e0238242. [Google Scholar] [CrossRef]
  41. Luan, Y.Z.; Chen, B.S.; Ma, K.S. 16S rDNA Gene Sequencing and Virulence of Oral Microbiome in Patients with Rheumatoid Arthritis. Arthritis Rheumatol. 2022. [Google Scholar] [CrossRef] [PubMed]
  42. Kuo, C.-F.; Lio, C.F.; Chen, H.-T.; Wang, Y.-T.T.; Ma, K.S.-K.; Chou, Y.T.; Chang, F.-C.; Tsai, S.-Y. Discordance of vancomycin minimum inhibitory concentration for methicillin-resistant Staphylococcus aureus at 2 μg/mL between Vitek II, E-test, and Broth Microdilution. PeerJ 2020, 8, e8963. [Google Scholar] [CrossRef] [PubMed]
  43. Wu, M.-C.; Ma, K.S.-K.; Chen, H.-H.; Huang, J.-Y.; Wei, J.C.-C. Relationship between Helicobacter pylori infection and psoriasis: A nationwide population-based longitudinal cohort study. Medicine 2020, 99, e20632. [Google Scholar] [CrossRef]
  44. Wang, C.-Y.; Lee, B.-S.; Jhang, Y.-T.; Ma, K.S.-K.; Huang, C.-P.; Fu, K.-L.; Lai, C.-H.; Tseng, W.-Y.; Kuo, M.Y.-P.; Chen, Y.-W. Er:YAG laser irradiation enhances bacterial and lipopolysaccharide clearance and human gingival fibroblast adhesion on titanium discs. Sci. Rep. 2021, 11, 23954. [Google Scholar] [CrossRef] [PubMed]
  45. Roivainen, M.; Viik-Kajander, M.; Palosuo, T.; Toivanen, P.; Leinonen, M.; Saikku, P.; Tenkanen, L.; Manninen, V.; Hovi, T.; MänttäriM. Infections, inflammation, and the risk of coronary heart disease. Circulation 2000, 101, 252–257. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Zhu, J.; Nieto, F.J.; Horne, B.D.; Anderson, J.L.; Muhlestein, J.B.; Epstein, S.E. Prospective study of pathogen burden and risk of myocardial infarction or death. Circulation 2001, 103, 45–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Chiang, C.H.; Chiang, C.H.; Pickering, J.W.; Stoyanov, K.M.; Chew, D.P.; Neumann, J.T.; Ojeda, F.; Sörensen, N.A.; Su, K.Y.; Kavsak, P.; et al. Performance of the European Society of Cardiology 0/1-Hour, 0/2-Hour, and 0/3-Hour Algorithms for Rapid Triage of Acute Myocardial Infarction: An International Collaborative Meta-analysis. Ann. Intern. Med. 2022, 175, 101–113. [Google Scholar] [CrossRef]
  48. Ma, K.S.-K.; Liou, Y.-J.; Huang, P.-H.; Lin, P.-S.; Chen, Y.-W.; Chang, R.-F. Identifying Medically-compromised Patients with Periodontitis-Associated Cardiovascular Diseases Using Convolutional Neural Network-facilitated Multilabel Classification of Panoramic Radiographs. In Proceedings of the 2021 International Conference on Applied Artificial Intelligence (ICAPAI), Halden, Norway, 19–21 May 2021. [Google Scholar]
  49. Ma, K.S.-K.; Chiang, C.-H.; Lopez, A.A.V.; Wang, L.-T. Cohort study of periodontitis-associated signaling pathways in myocardial infarction following atherosclerotic cardiovascular diseases. Metabolism 2021, 116, 154478. [Google Scholar] [CrossRef]
  50. Kuo, H.K.; Fujise, K. Human papillomavirus and cardiovascular disease among U.S. women in the National Health and Nutrition Examination Survey 2003 to 2006. J. Am. Coll. Cardiol. 2011, 58, 2001–2006. [Google Scholar] [CrossRef] [Green Version]
  51. Muñoz-Grajales, C.; Pineda, J.C. Pathophysiological Relationship between Infections and Systemic Vasculitis. Autoimmune Dis. 2015, 2015, 286783. [Google Scholar] [CrossRef] [Green Version]
  52. Aksu, K.; Keser, G.; Günaydin, G.; Ozbek, S.S.; Colakoğlu, Z.; Gümüşdiş, G.; Doganavşargil, E. Erectile dysfunction in Behçet‘s disease without neurological involvement: Two case reports. Rheumatology 2000, 39, 1429–1431. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Kaul, N.; Bhat, A.; Singh, R.; Singh, I. Erectile. Erectile Dysfunction in Behcet‘s Disease. Indian J. Dermatol. 2017, 62, 217–219. [Google Scholar] [CrossRef] [PubMed]
  54. Mohammadi, A.; Pfeifer, J.D.; Lewis, J.S., Jr. Association between human papillomavirus DNA and temporal arteritis. BMC Musculoskelet. Disord. 2012, 13, 132. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Chen, Y.-F.; Lin, C.-S.; Hong, C.-F.; Lee, D.-J.; Sun, C.; Lin, H.-H. Design of a Clinical Decision Support System for Predicting Erectile Dysfunction in Men Using NHIRD Dataset. IEEE J. Biomed. Health Inform. 2019, 23, 2127–2137. [Google Scholar] [CrossRef]
  56. Chang, C.-H.; Chueh, S.-C.J.; Wu, V.-C.; Chen, L.; Lin, Y.-H.; Hu, Y.-H.; Wu, K.-D.; Tsai, Y.-C. Risk of severe erectile dysfunction in primary hyperaldosteronism: A population-based propensity score matching cohort study. Surgery 2019, 165, 622–628. [Google Scholar] [CrossRef] [Green Version]
  57. Lipsitch, M.; Tchetgen Tchetgen, E.; Cohen, T. Negative controls: A tool for detecting confounding and bias in observational studies. Epidemiology 2010, 21, 383–388. [Google Scholar] [CrossRef] [PubMed]
  58. Nørgaard, M.; Ehrenstein, V.; Vandenbroucke, J.P. Confounding in observational studies based on large health care databases: Problems and potential solutions—A primer for the clinician. Clin. Epidemiol. 2017, 9, 185–193. [Google Scholar] [CrossRef] [Green Version]
  59. Hitchcock, J.R.; Cook, C.N.; Bobat, S.; Ross, E.; Flores-Langarica, A.; Lowe, K.L.; Khan, M.; Dominguez-Medina, C.C.; Lax, S.; Carvalho-Gaspar, M.; et al. Inflammation drives thrombosis after Salmonella infection via CLEC-2 on platelets. J. Clin. Investig. 2015, 125, 4429–4446. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Chang, R.; Wei, J.C.-C.; Lin, M.-C.; Hung, Y.-M.; Hung, C.-H. Risk of subsequent ischemic stroke in patients with nontyphoidal salmonellosis: A nationwide population-based cohort study. J. Infect. 2020, 81, 396–402. [Google Scholar] [CrossRef]
  61. Zippe, C.D.; Raina, R.; Massanyi, E.Z.; Agarwal, A.; Jones, J.S.; Ulchaker, J.; Klein, E.A. Sexual function after male radical cystectomy in a sexually active population. Urology 2004, 64, 682–685, discussion 685–686. [Google Scholar] [CrossRef]
  62. Modh, R.A.; Mulhall, J.P.; Gilbert, S.M. Sexual dysfunction after cystectomy and urinary diversion. Nat. Rev. Urol. 2014, 11, 445–453. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Shieh, S.-I.; Shu-Ling, H.; Huang, C.-Y.; Kao, C.-C.; Hung, S.-L.; Yang, H.-Y.; Tung, H.-Y. Sexual dysfunction in males following low anterior resection. J. Clin. Nurs. 2016, 25, 2348–2356. [Google Scholar] [CrossRef] [PubMed]
  64. Kessler, A.; Sollie, S.; Challacombe, B.; Briggs, K.; Van Hemelrijck, M. The global prevalence of erectile dysfunction: A review. BJU Int. 2019, 124, 587–599. [Google Scholar] [CrossRef] [PubMed]
  65. Ranjeva, S.L.; Baskerville, E.B.; Dukic, V.; Villa, L.L.; Lazcano-Ponce, E.; Giuliano, A.R.; Dwyer, G.; Cobey, S. Recurring infection with ecologically distinct HPV types can explain high prevalence and diversity. Proc. Natl. Acad. Sci. USA 2017, 114, 13573–13578. [Google Scholar] [CrossRef] [Green Version]
  66. Uysal, G.; Bas, S.; Gokulu, S.G.; Okcu, N.T.; Destegul, E. Sexual dysfunction in human papillomavirus positive females during reproductive age. Reprod. Syst. Sex. Disord. 2018, 7, 224. [Google Scholar]
  67. Rintala, M.A.M.; Pöllänen, P.P.; Nikkanen, V.P.; Grénman, S.E.; Syrjänen, S.M. Human papillomavirus DNA is found in the vas deferens. J. Infect. Dis. 2002, 185, 1664–1667. [Google Scholar] [CrossRef] [Green Version]
  68. Rintala, M.A.M.; E Greénman, S.; Pöllänen, P.P.; Suominen, J.J.; Syrjänen, S.M. Detection of high-risk HPV DNA in semen and its association with the quality of semen. Int. J. STD AIDS 2004, 15, 740–743. [Google Scholar] [CrossRef]
  69. Perino, A.; Giovannelli, L.; Schillaci, R.; Ruvolo, G.; Fiorentino, F.P.; Alimondi, P.; Cefalù, E.; Ammatuna, P. Human papillomavirus infection in couples undergoing in vitro fertilization procedures: Impact on reproductive outcomes. Fertil. Steril. 2011, 95, 1845–1848. [Google Scholar] [CrossRef]
  70. Lue, T.F. Erectile dysfunction. N. Engl. J. Med. 2000, 342, 1802–1813. [Google Scholar] [CrossRef]
  71. Ma, K.S.; Lai, J.; Bds, J.J.V.; Chiu, L.; Van Dyke, T.E.; Wei, J.C. Fibromyalgia and periodontitis: Bidirectional associations in population-based 15-year retrospective cohorts. J. Periodontol. 2021. [Google Scholar] [CrossRef]
  72. Ma, K.; Hasturk, H.; Carreras, I.; Dedeoglu, A.; Veeravalli, J.; Huang, J.; Kantarci, A.; Wei, J. Dementia and the Risk of Periodontitis: A Population-Based Cohort Study. J. Dent. Res. 2022, 101, 270–277. [Google Scholar] [CrossRef] [PubMed]
  73. Ma, K.S.; Wu, M.; Thota, E.; Wang, Y.; Alqaderi, H.E.; Wei, J.C. Tonsillectomy as a risk factor of periodontitis: A population-based cohort study. J. Periodontol. 2021. [Google Scholar] [CrossRef] [PubMed]
  74. Ma, K.S.; Thota, E.; Huang, J.; Huang, Y.; Wei, J.C. Onset of oral lichen planus following dental treatments: A nested case-control study. Oral. Dis. 2021. [Google Scholar] [CrossRef] [PubMed]
  75. Ma, K.S.-K.; Ralda, M.M.I.; Veeravalli, J.J.; Wang, L.-T.; Thota, E.; Huang, J.-Y.; Kao, C.-T.; Wei, J.C.-C.; Resnick, C.M. Patients with juvenile idiopathic arthritis are at increased risk for obstructive sleep apnoea: A population-based cohort study. Eur. J. Orthod. 2022, 44, 226–231. [Google Scholar] [CrossRef] [PubMed]
  76. Strand, A.; Rylander, E. Human papillomavirus. Subclinical and atypical manifestations. Dermatol. Clin. 1998, 16, 817–822. [Google Scholar] [CrossRef]
  77. Plummer, M.; Schiffman, M.; Castle, P.E.; Maucort-Boulch, D.; Wheeler, C.M. ALTS (Atypical Squamous Cells of Undetermined Significance/Low-Grade Squamous Intraepithelial Lesions Triage Study) Group A 2-Year Prospective Study of Human Papillomavirus Persistence among Women with a Cytological Diagnosis of Atypical Squamous Cells of Undetermined Significance or Low-Grade Squamous Intraepithelial Lesion. J. Infect. Dis. 2007, 195, 1582–1589. [Google Scholar] [CrossRef] [PubMed]
  78. Rodríguez, A.C.; Schiffman, M.; Herrero, R.; Wacholder, S.; Hildesheim, A.; Castle, P.E.; Solomon, D.; Burk, R.; On behalf of the Proyecto Epidemiologico Guanacaste Group. Rapid Clearance of Human Papillomavirus and Implications for Clinical Focus on Persistent Infections. JNCI J. Natl. Cancer Inst. 2008, 100, 513–517. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The cumulative incidence of erectile dysfunction in HPV-infected patients and non-HPV controls.
Figure 1. The cumulative incidence of erectile dysfunction in HPV-infected patients and non-HPV controls.
Jpm 12 00699 g001
Table 1. The baseline characteristics in patients with and without HPV. infection.
Table 1. The baseline characteristics in patients with and without HPV. infection.
VariablesNon-HPVHPVSMD
(N = 53,184)(N = 13,296)
n%n%
Age, year
18–3019,76837%489537%0.007
31–4010,96721%271020%0.006
41–50987319%239918%0.01
≥5112,57624%329225%0.03
mean, (SD)39.0(15.7)39.4(16.3)0.03
Comorbidities
hypertension982118%249219%0.007
diabetes mellitus 43988.0%11909.0%0.02
hyperlipidemia941918%245818%0.02
Stroke22434.0%6545.0%0.03
CAD6401.0%2172.0%0.04
CKD5481.0%1871.0%0.03
COPD49879%135110%0.03
Alcohol-related illness18383.0%4814.0%0.009
HIV800.0%320.0%0.02
AID7541.0%2382.0%0.03
Medication
α-blocker728314%190914%0.02
β-blocker37767%10178%0.02
CCB984219%247219%0.002
diuretics647012%167313%0.01
SMD: standard mean difference; CAD: cardiovascular disease; CKD: chronic kidney disease; COPD: chronic obstructive pulmonary disease; HIV: human immunodeficiency virus; AID: Autoimmune disease; CCB: calcium channel blocker.
Table 2. Hazard ratios for erectile dysfunction after diagnosis of HPV Infection.
Table 2. Hazard ratios for erectile dysfunction after diagnosis of HPV Infection.
Hazard Ratio (95% CI)
Primary analysis
Model 1 (adjusted comorbidities + comedications)1.64 (1.38, 1.95) ***
Model 2 (adjusted comorbidities + comedications + relevant operation)1.63 (1.37, 1.94) ***
Sensitivity analyses
HPV excluding stroke, cancer, AID1.59 (1.30, 1.95) ***
Alternative exposure (NTS)0.69 (0.31, 1.53)
***: p-value < 0.001. Model 1: adjusted by age and hypertension, diabetes mellitus, hyperlipidemia, stroke, COPD, alcohol-related illness, AID (SLE, RA, SS), and all medication (a-blocker, b-blocker, calcium channel blocker, diuretics); Use of medication was defined as the prescription for at least three outpatients visit before the index date. Model 2: adjusted covariates in model 1 + relevant operation (foreign body entering through orifice and LAR: Low anterior resection; APR: Abdominoperineal resection).
Table 3. The association of explanatory variables and the risk of erectile dysfunction.
Table 3. The association of explanatory variables and the risk of erectile dysfunction.
VariablesErectile Dysfunction (95% CI)
nPYIRcHR(95% CI)aHR
non-HPV432286,1371.511.00-1.00-
HPV18171,4782.531.69(1.42, 2.00) ***1.63(1.37, 1.94) ***
Age, year
18–3078141,3700.551.00-1.00-
31–409375,1741.242.25(1.67, 3.04) ***1.89(1.40, 2.56) ***
41–5017665,3172.694.93(3.78, 6.44) ***3.17(2.39, 4.2) ***
≥5126675,7543.516.49(5.04, 8.36) ***2.49(1.83, 3.39) ***
Comorbidities
hypertension
No416298,6321.391.00-1.00-
Yes19758,9833.342.43(2.05, 2.88) ***0.78(0.60, 1.00) *
diabetes mellitus
No501331,8831.511.00-1.00-
Yes11225,7324.352.92(2.38, 3.59) ***1.21(0.96, 1.52)
hyperlipidemia
No375298,6951.261.00-1.00-
Yes23858,9194.043.25(2.76, 3.82) ***1.63(1.34, 1.97) ***
Stroke
No581344,8851.681.00-1.00-
Yes3212,7292.511.52(1.06, 2.17) *0.60(0.41, 0.86) **
CAD
No602353,7901.71.00-
Yes1138252.881.70(0.94, 3.09)
CKD
No606354,7131.711.00-
Yes729022.411.45(0.69, 3.05)
COPD
No509325,8211.561.00-1.00-
Yes10431,7943.272.10(1.70, 2.59) ***1.03(0.82, 1.29)
Alcohol-related illness
No580347,6851.671.00-1.00-
Yes3399303.322.05(1.44, 2.91) ***1.43(1.01, 2.05) *
HIV
No613357,1091.72
Yes05060
AID
No590352,9181.671.00-1.00-
Yes2346964.92.95(1.94, 4.47) ***1.66(1.09, 2.53) *
Spinal cord injury
No613357,4451.71
Yes01700
Pelvic fracture
No613357,5991.71
Yes0160
Burn of unspecified degree of genitalia
No613357,6051.71
Yes090
Foreign body entering through orifice
No523321,8221.631.00-1.00-
Yes9035,7932.511.57(1.25, 1.96) ***1.20(0.96, 1.50)
Injury to nerves and spinal cord
No604354,3411.71.00-
Yes932732.751.64(0.85, 3.17)
Poisoning
No610354,6121.721.00-
Yes3300310.58(0.19, 1.82)
Medication
a-blocker
No316307,9841.031.00-1.00-
Yes29749,6315.985.81(4.96, 6.81) ***4.07(3.33, 4.97) ***
b-blocker
No502330,7561.521.00-1.00-
Yes11126,8594.132.70(2.19, 3.31) ***1.30(1.03, 1.64) *
CCB
No388291,2361.331.00-1.00-
Yes22566,3783.392.53(2.15, 2.98) ***1.02(0.80, 1.30)
diuretics
No474314,5051.511.00-1.00-
Yes13943,1103.222.12(1.76, 2.57) ***0.74(0.59, 0.93) *
Treatment
Radical cystectomy
No613357,5551.71
Yes0600
LAR + APR
No610357,3501.711.00-1.00-
Yes326511.36.56(2.11, 20.41) **2.70(0.86, 8.43)
*: p-value < 0.05; **: p-value < 0.01; ***: p-value < 0.001; PY: person-years; IR: incidence rate per 1000 person-years; cHR: crude hazard ratio; aHR: adjusted hazard ratio. CAD: cardiovascular disease; CKD: chronic kidney disease; COPD: chronic obstructive pulmonary disease; HIV: human immunodeficiency virus; AID: Autoimmune disease; CCB: calcium channel blocker; LAR: Low anterior resection; APR: Abdominoperineal resection. : adjusted by age and hypertension, diabetes mellitus, hyperlipidemia, stroke, COPD, alcohol-related illness, AID, and all medication.
Table 4. Subgroup analyses for the effect of HPV infection on ED risks under different stratification.
Table 4. Subgroup analyses for the effect of HPV infection on ED risks under different stratification.
VariablesNon-HPVHPV p for Interaction
nPYIRnPYIRcHR(95% CI)aHR (95% CI)
Age, year 0.40
18–3061113,4150.541727,9550.611.14(0.66, 1.94)1.14(0.67, 1.95)
31–406760,3351.112614,8391.751.61(1.03, 2.54) *1.65(1.05, 2.6) *
41–5012752,6462.414912,6713.871.61(1.16, 2.24) **1.63(1.17, 2.27) **
≥5117759,7402.968916,0135.561.87(1.45, 2.42) ***1.85(1.43, 2.39) ***
Comorbidities
hypertension 0.10
No303239,3031.2711359,3291.91.51(1.22, 1.87) ***1.5(1.21, 1.87) ***
Yes12946,8332.756812,1505.62.05(1.53, 2.76) ***1.99(1.48, 2.68) ***
diabetes mellitus 0.05
No364266,0141.3713765,8692.081.52(1.25, 1.86) ***1.5(1.23, 1.82) ***
Yes6820,1223.384456107.842.38(1.63, 3.49) ***2.39(1.63, 3.51) ***
hyperlipidemia 0.89
No265239,4301.1111059,2651.861.68(1.34, 2.1) ***1.69(1.36, 2.12) ***
Yes16746,7063.587112,2135.811.65(1.25, 2.18) ***1.58(1.19, 2.09) **
Stroke 0.001
No420276,4531.5216168,4332.351.56(1.3, 1.87) ***1.53(1.28, 1.84) ***
Yes1296841.242030466.575.18(2.53, 10.6) ***4.82(2.34, 9.92) ***
CAD 0.46
No424283,3491.517870,4412.531.7(1.42, 2.02) ***1.67(1.4, 1.99) ***
Yes827872.87310372.891.07(0.28, 4.07)1.3(0.32, 5.23)
CKD 0.27
No429283,9901.5117770,7232.51.66(1.4, 1.98) ***1.63(1.36, 1.94) ***
Yes321471.447555.33.64(0.81, 16.27)4.99(0.94, 26.6)
COPD 0.03
No370261,1781.4213964,6422.151.52(1.25, 1.85) ***1.5(1.23, 1.82) ***
Yes6224,9582.484268366.142.56(1.73, 3.81) ***2.51(1.69, 3.74) ***
Alcohol-related illness 0.17
No413278,3401.4816769,3452.411.63(1.36, 1.95) ***1.6(1.34, 1.92) ***
Yes1977962.441421346.562.68(1.34, 5.34) **2.52(1.25, 5.07) **
HIV 1.00
No432285,7901.5118171,3192.541.69(1.42, 2.01) ***1.64(1.38, 1.96) ***
Yes0347001590
AID 0.17
No420282,5711.4917070,3472.421.63(1.37, 1.95) ***1.62(1.35, 1.93) ***
Yes1235663.371111319.732.85(1.26, 6.46) *2.44(1.07, 5.6) *
Medication
a-blocker 0.47
No228246,7310.928861,2531.441.56(1.22, 1.99) ***1.55(1.21, 1.98) ***
Yes20439,4055.189310,2269.091.76(1.37, 2.25) ***1.75(1.37, 2.24) ***
b-blocker 0.82
No355264,7421.3414766,0142.231.66(1.37, 2.02) ***1.63(1.35, 1.98) ***
Yes7721,3943.63454646.221.75(1.16, 2.62) **1.7(1.13, 2.56) *
CCB 0.07
No283232,9961.2110558,2401.81.49(1.19, 1.86) ***1.45(1.16, 1.82) **
Yes14953,1402.87613,2385.742.06(1.56, 2.72) ***2.02(1.53, 2.67) ***
Diuretics 0.25
No340251,8481.3513462,6572.141.58(1.3, 1.93) ***1.54(1.26, 1.88) ***
Yes9234,2892.684788215.332.06(1.45, 2.93) ***2.07(1.45, 2.96) ***
*: p-value < 0.05; **: p-value < 0.01; ***: p-value < 0.001; PY: person-years; IR: incidence rate per 1000 person-years; cHR: crude hazard ratio; aHR: adjusted hazard ratio. CAD: cardiovascular disease; CKD: chronic kidney disease; COPD: chronic obstruction pulmonary disease; HIV: human immunodeficiency virus; AID: Autoimmune disease; CCB: calcium channel blocker; : adjusted by age and hypertension, diabetes mellitus, hyperlipidemia, stroke, COPD, alcohol-related illness, AID and all medication.
Table 5. Effect of HPV infection on ED development, stratified by follow-up duration to ED onset.
Table 5. Effect of HPV infection on ED development, stratified by follow-up duration to ED onset.
Non-HPVHPV
Follow-Up, YearnPYIR NPYIR cHR(95% CI)aHR (95% CI)
<15852,6910.114013,2100.301.67(1.41, 1.99) ***1.64(1.38, 1.95) ***
1–5238262,7050.098665,7570.131.51(1.24, 1.83) ***1.48(1.22, 1.80) ***
>5136212,5300.065553,0490.101.67(1.22, 2.29) **1.7(1.24, 2.34) ***
**: p-value < 0.01; ***: p-value < 0.001; PY: person-years; IR: incidence rate (per 1000 person-years): cHR: crude hazard ratio; aHR: adjusted hazard ratio. : adjusted by age and hypertension, diabetes mellitus, hyperlipidemia, stroke, COPD, alcohol-related illness, AID and all medication.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Juang, S.-E.; Ma, K.S.-K.; Kao, P.-E.; Wei, J.C.-C.; Yip, H.-T.; Chou, M.-C.; Hung, Y.-M.; Chin, N.-C. Human Papillomavirus Infection and the Risk of Erectile Dysfunction: A Nationwide Population-Based Matched Cohort Study. J. Pers. Med. 2022, 12, 699. https://doi.org/10.3390/jpm12050699

AMA Style

Juang S-E, Ma KS-K, Kao P-E, Wei JC-C, Yip H-T, Chou M-C, Hung Y-M, Chin N-C. Human Papillomavirus Infection and the Risk of Erectile Dysfunction: A Nationwide Population-Based Matched Cohort Study. Journal of Personalized Medicine. 2022; 12(5):699. https://doi.org/10.3390/jpm12050699

Chicago/Turabian Style

Juang, Sin-Ei, Kevin Sheng-Kai Ma, Pei-En Kao, James Cheng-Chung Wei, Hei-Tung Yip, Mei-Chia Chou, Yao-Min Hung, and Ning-Chien Chin. 2022. "Human Papillomavirus Infection and the Risk of Erectile Dysfunction: A Nationwide Population-Based Matched Cohort Study" Journal of Personalized Medicine 12, no. 5: 699. https://doi.org/10.3390/jpm12050699

APA Style

Juang, S. -E., Ma, K. S. -K., Kao, P. -E., Wei, J. C. -C., Yip, H. -T., Chou, M. -C., Hung, Y. -M., & Chin, N. -C. (2022). Human Papillomavirus Infection and the Risk of Erectile Dysfunction: A Nationwide Population-Based Matched Cohort Study. Journal of Personalized Medicine, 12(5), 699. https://doi.org/10.3390/jpm12050699

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop