Next Article in Journal
Nociception Level Index-Directed Erector Spinae Plane Block in Open Heart Surgery: A Randomized Controlled Clinical Trial
Previous Article in Journal
Could an Increased Percentage of Immature Granulocytes Accompanying Dyspepsia Predict COVID-19?
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Medicina
Volume 58
Issue 10
10.3390/medicina58101461
medicina-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Human Papillomavirus and Risk of Colorectal Cancer: An Analysis of Nationwide Claims Data

by
Chih-Hsiung Hsu
1,
Yu-Jyun Lin
1,
Yong-Chen Chen
2,3,
I-Lan Liu
1,
San-Lin You
2,3,
Je-Ming Hu
4,5,
Tzu-Chiao Lin
5,
Pi-Kai Chang
4,5,
Chao-Yang Chen
5,
Yu-Ching Chou
1,* and
Chien-An Sun
3,6,*
1
School of Public Health, National Defense Medical Center, Taipei City 114, Taiwan
2
Department of Medicine, College of Medicine, Fu-Jen Catholic University, New Taipei City 242, Taiwan
3
Data Science Center, College of Medicine, Fu-Jen Catholic University, New Taipei City 242, Taiwan
4
Graduate Institute of Medical Sciences, National Defense Medical Center, Taipei City 114, Taiwan
5
Division of Colorectal Surgery, Department of Surgery, Tri-Service General Hospital, National Defense Medical Center, Taipei City 114, Taiwan
6
Department of Public Health, College of Medicine, Fu-Jen Catholic University, New Taipei City 242, Taiwan
*
Authors to whom correspondence should be addressed.
Medicina 2022, 58(10), 1461; https://doi.org/10.3390/medicina58101461
Submission received: 25 August 2022 / Revised: 4 October 2022 / Accepted: 10 October 2022 / Published: 15 October 2022
Browse Figures
Versions Notes

Abstract

:
Background and Objectives: Although human papillomavirus (HPV) is a major etiology of cervical and anogenital cancers, whether it is associated with colorectal carcinogenesis is yet undetermined. Materials and Methods: The longitudinal association of HPV infection with colorectal cancer (CRC) was evaluated using 2000–2013 data from a nationwide Taiwanese claims database. In this retrospective cohort study, 358 patients with primary HPV diagnoses (HPV-infected cohort) and 1432 patients without such a diagnosis (HPV-uninfected cohort) were recruited between 2000 and 2006. Both cohorts were followed up to identify CRC incidences from 2006 to 2013. Hazard ratios (HRs) and their 95% confidence intervals (CIs) derived from Cox proportional hazards models were used to estimate the association between HPV and CRC risk. Results: The HPV-infected cohort had a significantly higher cumulative incidence of CRC than the HPV-uninfected cohort. The presence of HPV was associated with an increased risk of CRC (adjusted HR, 1.63; 95% CI, 1.02–3.62). Furthermore, the significant HPV–CRC risk association was evident in both sexes. Conclusions: This population-based cohort study reveals longitudinal evidence that HPV is associated with an increased risk of CRC. Further studies are required to verify the role of HPV in colorectal carcinogenesis.

1. Introduction

Colorectal cancer (CRC) is the third and second leading cause of cancer in men and women, respectively, based on the GLOBOCAN estimates of worldwide cancer incidence rates [1]. The incidence rates of CRC are rising rapidly in Asia, and the incidence rates of CRC in many Asian countries are similar to those in the western ones [2]. CRC is associated with lifestyle and inherited genetic factors [3]. In 2018, about 2.2 million cancer cases, or 13% of global cancer incidence rates, were attributable to infection [4]. The literature also suggests a link between CRC risk and infectious agents, especially human papillomavirus (HPV) [5,6,7].
A non-enveloped circular, double-stranded DNA virus, HPV, is the major cause of cervical cancer [8]. Similarly, the HPV–CRC cancer association was first reported by Kirgan et al. in 1990 [9], and HPV infection in relation to colorectal carcinogenesis has since been evaluated by numerous research groups [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32]. The possible association between CRC and HPV, however, remains controversial. Numerous studies have suggested an association of HPV infection with the development of CRC [10,11,12,13,14,15,16,18,22,24,27,28,30,31], but several investigations have indicated no evidence of this link [17,20,21,23,25,26,29,32]. The appropriate care for CRC and prevention of this disease would be strongly influenced by such a relationship. Therefore, we employed claims data from Taiwan’s nationwide National Health Insurance Research Database (NHIRD) to assess the association between HPV and CRC.

2. Methods

2.1. Data Source

We used retrospective data from the Taiwan NHIRD, which comprises medical claims data from the National Health Insurance (NHI) program. The NHI is a single-payer, universal, and compulsory healthcare program in which more than 98% of the 23 million residents of Taiwan are enrolled [33]. The NHIRD contains demographic data of insurants; the clinical data on drug prescriptions, procedures, surgeries, and diagnostic codes in the format of the International Classification of Disease, Ninth Revision, Clinical Modification (ICD-9-CM); and the use of inpatient and outpatient care facilities. The NHIRD has been used for high-quality epidemiological research [34,35], and the information on diagnoses, drug prescriptions, and hospitalizations has been shown to be of good validity [36,37,38]. We used data from the Longitudinal Health Insurance Database (LHID), which covers a cohort of 2 million NHI beneficiaries sampled at random from the NHIRD’s registry for the period of 1 January 2000 to 31 December 2013. The individuals in the NHIRD cohort and those in the LHID cohort do not significantly differ in terms of sex, age, or healthcare costs [33]. The data provided to researchers are entirely anonymous because of de-identification through the scrambling of the beneficiaries’ ID codes. Written or verbal patient consent was thus waived. The Institutional Review Board of Fu-Jen Catholic University (FJUIRB No: C104014, approval date: 2015-10-06) approved this study’s protocol.

2.2. Study Population

All patient data analyzed in the present study were abstracted from the LHID. Patients with an HPV diagnosis between 1 January 2000 and 31 December 2006 were identified on the basis of ICD-9-CM codes [39] 079.4, 078.11, 795.05, 795.09, 795.15, 795.19, 796.75, and 796.79. Initially, we enrolled 2024 individuals who had new HPV diagnoses (the HPV cohort). Given that people aged more than 40 years were considered at-risk for CRC in Taiwan [40], we excluded patients aged younger than 40 years (n = 1362). In addition, patients with diagnosed HPV infection prior to 1 January 2000 (n = 70); patients with a cancer diagnosis (ICD-9-CM codes 140–209, and 235–239; n = 65), a colon adenomas diagnosis (ICD-9-CM code 211, n = 9), or an inflammatory bowel disease diagnosis (ICD-9-CM codes 555 and 556, n = 12); and patients with a CRC diagnosis prior to HPV infection (n = 21) were also excluded. Finally, 465 patients with a new HPV diagnosis were included in this study as the HPV-infected cohort. To ensure the accurate identification of patients with HPV, we enrolled only those patients who had two or more outpatient visits with HPV diagnostic codes. The date of HPV diagnosis was assigned as the index date. To recruit an HPV-uninfected cohort, we randomly selected individuals aged 40 years or older without an HPV diagnosis between 1 January 2000 and 31 December 2006. Patients in the HPV-uninfected cohort were assigned a reference date corresponding to the index date of the matched patients in the HPV-infected cohort. To control for differences in underlying conditions between the two cohorts, we applied propensity score matching for the HPV-infected cohort to the matched HPV-uninfected cohort. The propensity score was calculated for each patient using a logistic regression model with covariates of age, sex, comorbidities, and the frequency of outpatient visits. The matching ratio was 1:4 for the HPV-infected to HPV-uninfected cohorts for the consideration of increasing the power and efficiency of statistical analyses [41]. Finally, 358 patients with HPV infection and 1432 patients without HPV infection were included in the analysis. Follow-up was initiated from the index date until 31 December 2013, death (based on withdrawal from the NHI), or a first primary diagnosis of CRC, whichever occurred first (Figure 1). For the HPV-infected and HPV-uninfected cohorts, the mean follow-up period was 6.8 (±2.1) and 7.4 (±1.6) years, respectively.

2.3. CRC Identification

The main outcome was the incidence of CRC, identified on the basis of a primary diagnosis with one of the following ICD-9-CM codes: 153, 153.0, 153.1, 153.2, 153.3, 153.4, 153.5, 153.6, 153.7, 153.8, 153.9, 154, 154.0, 154.1, 154.2, 154.3, and 154.8. We consulted the Registry for Catastrophic Illness Patient Database (RCIPD), which is an NHIRD subset, to ensure that the identifications were sufficiently accurate. We can thoroughly and accurately track the incident CRC cases using data in the RCIPD [38]. In Taiwan, a patient with cancer is not given a certificate of catastrophic illness unless they have an official hospital-issued certificate of the diagnosis that is supported by laboratory results, diagnostic imaging, or histological findings.

2.4. Covariate Assessment and Adjustment

The covariates of interest were patient demographics, comorbidities, general health, and healthcare utilization. The following comorbidities were identified from inpatient and outpatient files: diabetes mellitus (ICD-9-CM code 250), hypertension (ICD-9-CM codes 401, 402, 403, 404, and 405), chronic obstructive pulmonary disease (COPD) (ICD-9-CM codes 491, 492, and 496), hypercholesterolemia (ICD-9-CM codes 272.0, 272.1, 272.2, and 272.4), and alcohol-related conditions (ICD-9-CM codes 571.0, 571.1, 571.2, and 571.3 for alcoholic liver disease, 303 for alcohol dependence, 305.0 for alcohol abuse, and 291 for alcohol-induced mental disorders). Comorbidities were defined in a patient if he or she was diagnosed with any of the aforementioned diseases on at least two or more outpatient claims or one or more inpatient claims in the study period. In the present study, hypertension, diabetes mellitus, or hypercholesterolemia was used as a surrogate measure of obesity. COPD was considered a proxy indicator of smoking, and alcohol-related conditions were utilized as a substitute variable for alcohol intake. In addition, general health status was indicated by a patient’s CCI score. In this scoring system, the presence of major concomitant diseases is summed, with each disease appropriately weighted; previous studies have validated the use of the CCI with administrative databases containing ICD-9-CM codes [42]. Furthermore, healthcare utilization was reflected by the frequency of outpatient visits. These covariates were included in the regression models for adjustment.

2.5. Statistical Analysis

We used the chi-squared or Student t-test to determine between cohort differences in the baseline comorbidities and demographic characteristics. Cumulative CRC incidence was evaluated using a Kaplan–Meier curve, and the cumulative CRC risks in the cohorts were compared through the log-rank test. Hazard ratios (HRs) and 95% confidence intervals (CIs) derived from the multivariable Cox proportional hazards models were used to assess HPV–CRC risk association. We evaluated the proportionality assumption of Cox models using the log minus log plot of survival [43]. Statistical significance was considered to be indicated by a two-tailed p < 0.05. We performed all analyses in SAS version 9.4 (SAS Institute, Cary, NC, USA).

3. Results

The characteristics of the study cohorts are presented in Table 1. The mean age (± standard deviation) of the HPV-infected and uninfected cohorts was 51.9 (±10.1) and 51.8 (±10.0) years, respectively. There were no significant differences in the distributions of age, sex, CCI score, frequency of outpatient visits, and comorbidities between the HPV-infected cohort and the HPV-uninfected cohort due to the propensity score matching schemes.
The mean follow-up period was 6.9 years. During the follow-up, 7 patients were identified as incident CRC in the HPV-infected cohort and 13 were noted in the HPV-uninfected cohort for incidence rates of 26.48 and 14.17 per 10,000 person-years, respectively. The HPV-infected cohort had a significantly greater risk of CRC in comparison with the HPV-uninfected cohort, with an adjusted HR of 1.63 (95% CI, 1.02–3.62) (Table 2). Of note, the significant HPV–CRC risk association was evident in both sexes (Table 2). The Kaplan–Meier curves for the cumulative incidence of CRC among the two cohorts are displayed in Figure 2. The cumulative incidence of CRC is significantly higher in the HPV-infected cohort than in the HPV-uninfected cohort (p = 0.036).
A differential pattern of association between HPV infection and CRC risk was demonstrated according to the frequency of clinical visits. Relative to patients without HPV infection, a higher HR was noted among patients having the highest frequency of clinical visits for a diagnosed HPV infection (adjusted HR, 2.03; 95% CI, 1.02–7.18). A trend, with a significant increase in HR with increasing frequency of clinical visits for a diagnosed HPV infection, was noted (p = 0.004) (Table 3).

4. Discussion

Hereditary and environmental factors contribute to the development of CRC and, thus, affect CRC risk. In particular, the greatest increase in risk is based on inherited susceptibility [44]. If this cancer were found to have an infectious etiology, the potential for diagnostic tools and therapeutic interventions would be substantial. In the present large-scale population-based cohort study, we documented a significant association of HPV with increased CRC risk in both sexes.
HPV is widely recognized to play an important role in the development of some carcinoma types, including cervical cancer, which is positive for HPV in 90% of cases, and head and neck tumors, for which the percentage is approximately 30% [45,46]. Although how HPV precisely affects CRC risk has been the topic of numerous studies, an association between HPV and CRC is controversial, with several publications suggesting that HPV infection of the colonic mucosa may contribute to the development of CRC [10,11,12,13,14,15,16,18,22,24,27,28,30,31]. However, there have been several investigations that failed to establish a link between HPV and CRC [17,20,21,23,25,26,29,32]. These discrepant findings could result from geographic differences in the study populations, regional differences in HPV prevalence, or variations in methods of HPV detection and choice of the material for analysis (fresh/paraffin). The findings of this study support the findings of investigations reporting that the risk of CRC is higher in individuals with HPV [10,11,12,13,14,15,16,18,22,24,27,28,30,31]. Zhang et al. presented evidence that, in the Chinese population, HPV infection is related to CRC risk [31]. In addition, in their HPV-type-specific polymerase chain reaction (PCR) analysis of tissues from Taiwanese patients, Cheng et al. determined that 29% of adenomas and 53% of colorectal carcinomas contained HPV DNA [47]. Indeed, several studies have pointed out that high-risk HPVs are present in human CRCs, specifically types 16, 18, 31, 33, and 35 [19,24,48,49]. Furthermore, several meta-analyses have indicated associations of high-risk HPV types 16, 18, 31, 33, and 35 with CRC risk [6,19,24]. It is well known that high-risk HPV early proteins, including E5, E6, and E7 oncoproteins, increase cellular alteration and probably lead to HPV-induced carcinogenesis [50,51,52]. The high prevalence of infection in colorectal tumors with a high-risk HPV type (HPV-16) and evidence of viral genome integration [19,24,27] support the possibility that HPV may have a role in colorectal carcinogenesis. Collectively, these data quantitatively support the notion that HPV infection is associated with an elevated risk of CRC.
We were able to analyze an unselected study population in a real-life setting because we employed data from large and representative population-based NHI cohorts; this is a strength of the study. A further strength is that we employed routine database records and could thus minimize selection and nonresponse biases. However, the limitations of this study must be considered in the interpretation of the results. Data on potential confounders are often lacking from insurance claims data [53], and this could have affected the present study. Information on the virus types of HPV or variables related to lifestyle, such as alcohol consumption and cigarette smoking, were not available in the studied claims dataset. Furthermore, data on sexual behaviors, HPV vaccination, and body mass index were also not factored in the claims database. Thus, we were unable to account for these factors in our models, and it is likely that our findings could be subject to residual confounding. Of note, HPV was significantly associated with cervical cancer and head and neck cancers [45,46]. Thus, ignoring competing events in the analysis of HPV–CRC association may lead to biased results. Overall, although we did not have data on the types of HPV, Cheng et al. found HPV DNA, especially type 16, in 29.7% and 52.9% of adenomas and carcinomas, respectively, after PCR and Southern blot hybridization of colorectal tumors from Taiwanese patients [47]. These findings implicate HPV DNA in the development of colorectal carcinomas in the context of Taiwanese patients.

5. Conclusions

In conclusion, despite the limitations described above, our study provides longitudinal evidence that HPV infection is associated with an increased risk of developing CRC in the Taiwanese population. It should be pointed out, however, that this observation needs to be substantiated by further studies to elucidate the role of HPV in colorectal carcinogenesis.

Author Contributions

Conceptualization, C.-H.H., Y.-C.C. (Yu-Ching Chou), and C.-A.S.; Data curation, Y.-C.C. (Yong-Chen Chen); Formal analysis, Y.-J.L.; Investigation, Y.-C.C. (Yong-Chen Chen), I.-L.L., S.-L.Y., J.-M.H., T.-C.L., P.-K.C., and C.-Y.C.; Methodology, C.-H.H., Y.-C.C. (Yu-Ching Chou), and C.-A.S.; Project administration, C.-A.S.; Validation, C.-H.H., Y.-J.L., Y.-C.C. (Yong-Chen Chen), I.-L.L., S.-L.Y., J.-M.H., T.-C.L., P.-K.C., C.-Y.C., Y.-C.C. (Yu-Ching Chou), and C.-A.S.; Writing—original draft, C.-H.H. and C.-A.S.; Writing—review and editing, Y.-C.C. (Yu-Ching Chou) and C.-A.S. All authors will be informed about each step of manuscript processing, including submission, revision, revision reminder, etc., via emails from our system or the assigned Assistant Editor. 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 Fu-Jen Catholic University (FJU-IRB No: C104014, approval date: 2015-10-06).

Informed Consent Statement

Not applicable.

Data Availability Statement

No additional data are available.

Acknowledgments

This manuscript was edited by Wallace Academic Editing. The authors thank the enrollees of the National Health Insurance Research Database for their important contributions.

Conflicts of Interest

The authors declare that they have no potential conflicts of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; Laversanne, M.; Soerjomataram, I.; Jemal, A.; Bray, F. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2021, 71, 209–249. [Google Scholar] [CrossRef] [PubMed]
  2. Onyoh, E.F.; Hsu, W.F.; Chang, L.C.; Lee, Y.C.; Wu, M.S.; Chiu, H.M. The Rise of Colorectal Cancer in Asia: Epidemiology, Screening, and Management. Curr. Gastroenterol. Rep. 2019, 21, 36. [Google Scholar] [CrossRef] [PubMed]
  3. Johnson, C.M.; Wei, C.; Ensor, J.E.; Smolenski, D.J.; Amos, C.I.; Levin, B.; Berry, D.A. Meta-analyses of colorectal cancer risk factors. Cancer Causes Control 2013, 24, 1207–1222. [Google Scholar] [CrossRef] [Green Version]
  4. de Martel, C.; Georges, D.; Bray, F.; Ferlay, J.; Clifford, G.M. Global burden of cancer attributable to infections in 2018: A worldwide incidence analysis. Lancet Glob. Health 2020, 8, e180–e190. [Google Scholar] [CrossRef] [Green Version]
  5. Antonic, V.; Stojadinovic, A.; Kester, K.E.; Weina, P.J.; Brücher, B.L.; Protic, M.; Avital, I.; Izadjoo, M. Significance of infectious agents in colorectal cancer development. J. Cancer 2013, 4, 227–240. [Google Scholar] [CrossRef]
  6. Chen, H.; Chen, X.Z.; Waterboer, T.; Castro, F.A.; Brenner, H. Viral infections and colorectal cancer: A systematic review of epidemiological studies. Int. J. Cancer 2015, 137, 12–24. [Google Scholar] [CrossRef]
  7. Damin, D.C.; Ziegelmann, P.K.; Damin, A.P. Human papillomavirus infection and colorectal cancer risk: A meta-analysis. Colorectal Dis. 2013, 15, e420–e428. [Google Scholar] [CrossRef]
  8. Crosbie, E.J.; Einstein, M.H.; Franceschi, S.; Kitchener, H.C. Human papillomavirus and cervical cancer. Lancet 2013, 382, 889–899. [Google Scholar] [CrossRef]
  9. Kirgan, D.; Manalo, P.; Hall, M.; McGregor, B. Association of human papillomavirus and colon neoplasms. Arch. Surg. 1990, 125, 862–865. [Google Scholar] [CrossRef]
  10. Ibragimova, M.K.; Tsyganov, M.M.; Litviakov, N.V. Human papillomavirus and colorectal cancer. Med. Oncol. 2018, 35, 140. [Google Scholar] [CrossRef]
  11. Libera, L.S.; de Siqueira, T.; Santos, I.L.; Porto Ramos, J.E.; Milhomen, A.X.; Alencar, R.C. Detection of human papillomavirus and the role of P16INK4a in colorectal carcinomas. PLoS ONE 2020, 15, e0235065. [Google Scholar]
  12. Bodaghi, S.; Yamanegi, K.; Xiao, S.Y.; Da Costa, M.; Palefsky, J.M.; Zheng, Z.M. Colorectal papillomavirus infection in patients with colorectal cancer. Clin. Cancer Res. 2005, 11, 2862–2867. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pérez, L.O.; Abba, M.C.; Laguens, R.M.; Golijow, C.D. Analysis of adenocarcinoma of the colon and rectum: Detection of human papillomavirus (HPV) DNA by polymerase chain reaction. Colorectal Dis. 2005, 7, 492–495. [Google Scholar] [CrossRef] [PubMed]
  14. Damin, D.C.; Caetano, M.B.; Rosito, M.A.; Schwartsmann, G.; Damin, A.S.; Frazzon, A.P.; Ruppenthal, R.D.; Alexandre, C.O. Evidence for an association of human papillomavirus infection and colorectal cancer. Eur. J. Surg. Oncol. 2007, 33, 569–574. [Google Scholar] [CrossRef]
  15. Militello, V.; Trevisan, M.; Squarzon, L.; Biasolo, M.A.; Rugge, M.; Militello, C.; Palù, G.; Barzon, L. Investigation on the presence of polyomavirus, herpesvirus, and papillomavirus sequences in colorectal neoplasms and their association with cancer. Int. J. Cancer 2009, 124, 2501–2503. [Google Scholar] [CrossRef]
  16. Pérez, L.O.; Barbisan, G.; Ottino, A.; Pianzola, H.; Golijow, C.D. Human papillomavirus DNA and oncogene alterations in colorectal tumors. Pathol. Oncol. Res. 2010, 16, 461–468. [Google Scholar] [CrossRef]
  17. Gornick, M.C.; Castellsague, X.; Sanchez, G.; Giordano, T.J.; Vinco, M.; Greenson, J.K.; Capella, G.; Raskin, L.; Rennert, G.; Gruber, S.B.; et al. Human papillomavirus is not associated with colorectal cancer in a large international study. Cancer Causes Control 2010, 21, 737–743. [Google Scholar] [CrossRef]
  18. Liu, F.; Mou, X.; Zhao, N.; Lin, J.; Teng, L.; Xiang, C. Prevalence of human papillomavirus in Chinese patients with colorectal cancer. Colorectal Dis. 2011, 13, 865–871. [Google Scholar] [CrossRef]
  19. Lorenzon, L.; Ferri, M.; Pilozzi, E.; Torrisi, M.R.; Ziparo, V.; French, D. Human papillomavirus and colorectal cancer: Evidences and pitfalls of published literature. Int. J. Colorectal Dis. 2011, 26, 135–142. [Google Scholar] [CrossRef]
  20. Snietura, M.; Waniczek, D.; Nowakowska-Zajdel, E.; Piglowski, W.; Kopec, A.; Muc-Wierzgon, M. Does human papilloma virus participate in colorectal carcinogenesis? J. Biol. Regul. Homeost. Agents 2012, 26, 757–762. [Google Scholar]
  21. Burnett-Hartman, A.N.; Feng, Q.; Popov, V.; Kalidindi, A.; Newcomb, P.A. Human papillomavirus DNA is rarely detected in colorectal carcinomas and not associated with microsatellite instability: The Seattle colon cancer family registry. Cancer Epidemiol. Biomark. Prev. 2013, 22, 317–319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ranjbar, R.; Saberfar, E.; Shamsaie, A.; Ghasemian, E. The aetiological role of human papillomavirus in colorectal carcinoma: An Iranian population- based case control study. Asian Pac. J. Cancer Prev. 2014, 15, 1521–1525. [Google Scholar] [CrossRef]
  23. Taherian, H.; Tafvizi, F.; Fard, Z.T.; Abdirad, A. Lack of association between human papillomavirus infection and colorectal cancer. Prz. Gastroenterol. 2014, 9, 280–284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Baandrup, L.; Thomsen, L.T.; Olesen, T.B.; Andersen, K.K.; Norrild, B.; Kjaer, S.K. The prevalence of human papillomavirus in colorectal adenomas and adenocarcinomas: A systematic review and meta-analysis. Eur. J. Cancer 2014, 50, 1446–1461. [Google Scholar] [CrossRef]
  25. Aghakhani, A.; Hamkar, R.; Ramezani, A.; Bidari-Zerehpoosh, F.; Sabeti, S.; Ghavami, N.; Banifazl, M.; Rashidi, N.; Eslamifar, A. Lack of human papillomavirus DNA in colon adenocarcinama and adenoma. J. Cancer Res. Ther. 2014, 10, 531–534. [Google Scholar] [PubMed]
  26. Lorenzon, L.; Mazzetta, F.; Pilozzi, E.; Uggeri, G.; Torrisi, M.R.; Ferri, M.; Ziparo, V.; French, D. Human papillomavirus does not have a causal role in colorectal carcinogenesis. World J. Gastroenterol. 2015, 21, 342–350. [Google Scholar] [CrossRef] [PubMed]
  27. Bernabe-Dones, R.D.; Gonzalez-Pons, M.; Villar-Prados, A.; Lacourt-Ventura, M.; Rodríguez-Arroyo, H.; Fonseca-Williams, S.; Velazquez, F.E.; Diaz-Algorri, Y.; Lopez-Diaz, S.M.; Rodríguez, N.; et al. High Prevalence of Human Papillomavirus in Colorectal Cancer in Hispanics: A Case-Control Study. Gastroenterol. Res. Pract. 2016, 2016, 7896716. [Google Scholar] [CrossRef] [Green Version]
  28. Khabaz, M.N. HPV and the Development of Colorectal Cancer. Glob. J. Health. Sci. 2017, 9, 251. [Google Scholar] [CrossRef]
  29. Vuitton, L.; Jaillet, C.; Jacquin, E.; Monnien, F.; Heberle, M.; Mihai, M.I.; Lassabe, C.; Raffoul, J.; Puyraveau, M.; Lakkis, Z.; et al. Human papillomaviruses in colorectal cancers: A case-control study in western patients. Dig. Liver Dis. 2017, 49, 446–450. [Google Scholar] [CrossRef]
  30. Malekpour Afshar, R.; Deldar, Z.; Mollaei, H.R.; Arabzadeh, S.A.; Iranpour, M. Evaluation of HPV DNA positivity in colorectal cancer patients in Kerman, Southeast Iran. Asian Pac. J. Cancer Prev. 2018, 19, 193–198. [Google Scholar]
  31. Zhang, X.H.; Wang, W.; Wang, Y.Q.; Jia, D.F.; Zhu, L. Human papillomavirus infection and colorectal cancer in the Chinese population: A meta-analysis. Colorectal Dis. 2018, 20, 961–969. [Google Scholar] [CrossRef] [PubMed]
  32. Biesaga, B.; Janecka-Widła, A.; Kołodziej-Rzepa, M.; Słonina, D.; Darasz, Z.; Gasińska, A. The prevalence of HPV infection in rectal cancer—Report from South—Central Poland (Cracow region). Pathol. Res. Pract. 2019, 215, 152513. [Google Scholar] [CrossRef] [PubMed]
  33. Lin, L.Y.; Warren-Gash, C.; Smeeth, L.; Chen, P.C. Data resource profile: The National Health Insurance Research Database (NHIRD). Epidemiol. Health 2018, 40, e2018062. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, C.Y.; Chen, Y.J.; Ho, H.J.; Hsu, Y.C.; Kuo, K.N.; Wu, M.S.; Lin, J.T. Association between nucleoside analogues and risk of hepatitis B virus–related hepatocellular carcinoma recurrence following liver resection. JAMA 2012, 308, 1906–1914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Chang, S.H.; Chou, I.J.; Yeh, Y.H.; Chiou, M.J.; Wen, M.S.; Kuo, C.T.; See, L.C.; Kuo, C.F. Association Between Use of Non-Vitamin K Oral Anticoagulants With and Without Concurrent Medications and Risk of Major Bleeding in Nonvalvular Atrial Fibrillation. JAMA 2017, 318, 1250–1259. [Google Scholar] [CrossRef]
  36. Gau, C.S.; Chang, I.S.; Lin Wu, F.L.; Yu, H.T.; Huang, Y.W.; Chi, C.L.; Chien, S.Y.; Lin, K.M.; Liu, M.Y.; Wang, H.P. Usage of the claim database of national health insurance programme for analysis of cisapride-erythromycin co-medication in Taiwan. Pharmacoepidemiol. Drug Saf. 2007, 16, 86–95. [Google Scholar] [CrossRef]
  37. Cheng, C.L.; Kao, Y.H.; Lin, S.J.; Lee, C.H.; Lai, M.L. Validation of the National Health Insurance Research Database with ischemic stroke cases in Taiwan. Pharmacoepidemiol. Drug Saf. 2011, 20, 236–242. [Google Scholar] [CrossRef]
  38. Kao, W.H.; Hong, J.H.; See, L.C.; Yu, H.P.; Hsu, J.T.; Chou, I.J.; Chou, W.C.; Chiou, M.J.; Wang, C.C.; Kuo, C.F. Validity of cancer diagnosis in the National Health Insurance database compared with the linked National Cancer Registry in Taiwan. Pharmacoepidemiol. Drug Saf. 2018, 27, 1060–1066. [Google Scholar] [CrossRef]
  39. Gavrielov-Yusim, N.; Friger, M. Use of administrative medical databases in population-based research. J Epidemiol. Community Health 2014, 68, 283–287. [Google Scholar] [CrossRef]
  40. Chen, C.J.; You, S.L.; Lin, L.H.; Hsu, W.L.; Yang, Y.W. Cancer epidemiology and control in Taiwan: A brief review. Jpn. J. Clin. Oncol. 2002, 32 (Suppl. S1), S66–S81. [Google Scholar] [CrossRef] [Green Version]
  41. Sturmer, T.; Brenner, H. Degree of matching and gain in power and efficiency in case-control studies. Epidemiology 2001, 12, 101–108. [Google Scholar] [CrossRef] [PubMed]
  42. Deyo, R.A.; Cherkin, D.C.; Ciol, M.A. Adapting a clinical comorbidity index for use with ICD-9-CM administrative databases. J. Clin. Epidemiol. 1992, 45, 613–619. [Google Scholar] [CrossRef]
  43. Therneau, T.M.; Grambsch, P.M. Modeling Survival Data: Extending the Cox Model; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2013. [Google Scholar]
  44. Rocha, J.L.U.; Ferreiro, R.S.; Suarez, R.F.; Vazquez, J.T. Primary prevention of CRC. In Fundations of Colorectal Cancer; Sierra, A.P., Ed.; Academic Press: Cambridge, MA, USA, 2021. [Google Scholar]
  45. Schiffman, M.H.; Bauer, H.M.; Hoover, R.N.; Glass, A.G.; Cadell, D.M.; Rush, B.B.; Scott, D.R.; Sherman, M.E.; Kurman, R.J.; Wacholder, S.; et al. Epidemiologic evidence showing that human papillomavirus infection causes most cervical intraepithelial neoplasia. J. Natl. Cancer Inst. 1993, 85, 958–964. [Google Scholar] [CrossRef] [PubMed]
  46. Bychkov, V.A.; Nikitina, E.G.; Ibragimova, M.K.; Kaigorodova, E.V.; Choinzonov, E.L.; Litviakov, N.V. Comprehensive meta-analytical summary on human papillomavirus association with head and neck cancer. Exp. Oncol. 2016, 38, 68–72. [Google Scholar] [CrossRef]
  47. Cheng, J.Y.; Sheu, L.F.; Meng, C.L.; Lee, W.H.; Lin, J.C. Detection of human papillomavirus DNA in colorectal carcinomas by polymerase chain reaction. Gut 1995, 37, 87–90. [Google Scholar] [CrossRef] [Green Version]
  48. Ghabreau, L.; Segal, E.D.; Yasmeen, A.; Kassab, A.; Akil, N.; AI Moustafa, A.E. High-risk human papillomavirus infections in colorectal cancer in the Syrian population and their association with Fascin, P-cadherin and Id-1 expressions: A tissue microarray study. Clin. Cancer Investig. 2012, 1, 26–30. [Google Scholar]
  49. Al Moustafa, A.E.; Ghabreau, L.; Akil, N.; Rastam, S.; Alachkar, A.; Yasmeen, A. High-risk HPVs and human carcinomas in the Syrian population. Front. Oncol. 2014, 4, 68. [Google Scholar] [CrossRef] [Green Version]
  50. Jeon, S.; Allen-Hoffmann, B.L.; Lambert, P.F. Integration of human papillomavirus type 16 into the human genome correlates with a selective growth advantage of cells. J. Virol. 1995, 69, 2989–2997. [Google Scholar] [CrossRef] [Green Version]
  51. Grm, H.S.; Massimi, P.; Gammoh, N.; Banks, L. Crosstalk between the human papillomavirus E2 transcriptional activator and the E6 oncoprotein. Oncogene 2005, 24, 5149–5164. [Google Scholar] [CrossRef]
  52. Al Moustafa, A.E. E5 and E6/E7 of high-risk HPVs cooperate to enhance cancer progression through EMT initiation. Cell Adhes. Migr. 2015, 9, 392–393. [Google Scholar] [CrossRef]
  53. Hyman, J. The limitations of using insurance data for research. J. Am. Dent. Assoc. 2015, 146, 283–285. [Google Scholar] [CrossRef] [PubMed]
Medicina 58 01461 g001 550
Figure 1. Study flowchart. LHID, Longitudinal Health Insurance Database; HPV, human papillomavirus; CRC, colorectal cancer; CCI, Charlson comorbidity index.
Figure 1. Study flowchart. LHID, Longitudinal Health Insurance Database; HPV, human papillomavirus; CRC, colorectal cancer; CCI, Charlson comorbidity index.
Medicina 58 01461 g001
Medicina 58 01461 g002 550
Figure 2. Kaplan–Meier analysis of the cumulative incidence of colorectal cancer in the HPV-infected cohort and the HPV-uninfected cohort.
Figure 2. Kaplan–Meier analysis of the cumulative incidence of colorectal cancer in the HPV-infected cohort and the HPV-uninfected cohort.
Medicina 58 01461 g002
Table
Table 1. Baseline demographics and comorbidities in the HPV-infected and HPV-uninfected cohorts.
Table 1. Baseline demographics and comorbidities in the HPV-infected and HPV-uninfected cohorts.
VariableHPV-Infected Cohort
n = 358 		HPV-Uninfected Cohort
n = 1432		p-Value
Age, years (mean (SD))51.9 (10.8)51.8 (10.0)0.817
CCI (mean (SD))0.3 (0.7)0.3 (1.1)0.914
No. of annual outpatient visits (mean (SD))15.3 (37.0)12.2 (57.6)0.513
Sex (No. (%)) 0.857
Female146 (40.8)573 (40.2)
Male212 (59.2)857 (59.8)
Comorbidities (No. (%))
COPD78 (21.8)275 (19.2)0.267
Hypertension150 (41.9)526 (36.7)0.092
Diabetes mellitus81 (22.6)316 (22.1)0.831
Hypercholesterolemia120 (33.5)426 (29.8)0.178
Alcoholic liver disease8 (2.2)25 (1.8)0.513
Alcohol dependence6 (1.7)23 (1.6)0.987
Alcohol abuse2 (0.6)10 (0.7)0.976
Alcohol-induced mental disorders3 (0.8)8 (0.6)0.768
CCI, Charlson comorbidity index; COPD, chronic obstructive pulmonary.
Table
Table 2. Incidence rates of colorectal cancer (CRC) in study cohorts and multivariable Cox proportional hazards regression model analysis for the association between HPV infection and the risk of CRC.
Table 2. Incidence rates of colorectal cancer (CRC) in study cohorts and multivariable Cox proportional hazards regression model analysis for the association between HPV infection and the risk of CRC.
VariableHPV-Infected CohortHPV-Uninfected CohortAdjusted HR (95% CI)
No. of CRC CasesPYsIncidence Rate
(per 10,000)		No. of CRC CasesPYsIncidence Rate
(per 10,000)
Total7264426.48 13917714.171.63 (1.02–3.62)
Sex
Female3109627.375400312.49 1.91 (1.09–6.36)
Male4154825.84 8570914.01 1.54 (1.08–4.01)
PYs: person-years; HR: hazard ratio; CI: confidence interval. Hazard ratios for the whole group analysis were adjusted for age, sex, the Charlson comorbidity index, comorbidities of chronic obstructive pulmonary disease, hypertension, diabetes mellitus, hypercholesterolemia, alcohol-related conditions, number of annual outpatient visits, and the index date. Hazard ratios for the sex-specific analyses were adjusted for age, the Charlson comorbidity index, comorbidities of chronic obstructive pulmonary disease, hypertension, diabetes mellitus, hypercholesterolemia, alcohol-related conditions, number of annual outpatient visits, and the index date.
Table
Table 3. Multivariable Cox proportional hazards regression model analysis for the association between the risk of colorectal cancer (CRC) associated with frequency of clinical visits for a diagnosed HPV infection.
Table 3. Multivariable Cox proportional hazards regression model analysis for the association between the risk of colorectal cancer (CRC) associated with frequency of clinical visits for a diagnosed HPV infection.
Frequency of Clinical Visits for Diagnosed HPV InfectionNo. of CRC CasesPYsIncidence Rate
(per 10,000)		Adjusted HR
(95% CI)
HPV-uninfected cohort13971110.17 1.00 (reference)
1–24159923.00 1.81 (0.59–5.57)
>23104549.34 2.03 (1.02–7.18)
p for trend 0.004
PYs: person-years; HR, hazard ratio; CI, confidence interval. Hazard ratios were adjusted for age, sex, the Charlson comorbidity index, comorbidities of chronic obstructive pulmonary disease, hypertension, diabetes mellitus, hypercholesterolemia, alcohol-related conditions, number of annual outpatient visits, and the index date.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Hsu, C.-H.; Lin, Y.-J.; Chen, Y.-C.; Liu, I.-L.; You, S.-L.; Hu, J.-M.; Lin, T.-C.; Chang, P.-K.; Chen, C.-Y.; Chou, Y.-C.; et al. Human Papillomavirus and Risk of Colorectal Cancer: An Analysis of Nationwide Claims Data. Medicina 2022, 58, 1461. https://doi.org/10.3390/medicina58101461

AMA Style

Hsu C-H, Lin Y-J, Chen Y-C, Liu I-L, You S-L, Hu J-M, Lin T-C, Chang P-K, Chen C-Y, Chou Y-C, et al. Human Papillomavirus and Risk of Colorectal Cancer: An Analysis of Nationwide Claims Data. Medicina. 2022; 58(10):1461. https://doi.org/10.3390/medicina58101461

Chicago/Turabian Style

Hsu, Chih-Hsiung, Yu-Jyun Lin, Yong-Chen Chen, I-Lan Liu, San-Lin You, Je-Ming Hu, Tzu-Chiao Lin, Pi-Kai Chang, Chao-Yang Chen, Yu-Ching Chou, and et al. 2022. "Human Papillomavirus and Risk of Colorectal Cancer: An Analysis of Nationwide Claims Data" Medicina 58, no. 10: 1461. https://doi.org/10.3390/medicina58101461

APA Style

Hsu, C. -H., Lin, Y. -J., Chen, Y. -C., Liu, I. -L., You, S. -L., Hu, J. -M., Lin, T. -C., Chang, P. -K., Chen, C. -Y., Chou, Y. -C., & Sun, C. -A. (2022). Human Papillomavirus and Risk of Colorectal Cancer: An Analysis of Nationwide Claims Data. Medicina, 58(10), 1461. https://doi.org/10.3390/medicina58101461

Article Metrics

Back to TopTop