Human Papillomavirus Is Associated With Adenocarcinoma of Lung: A Population-Based Cohort Study

Background Recent trends in the incidence of lung cancer have been reported despite the decreasing rate of smoking. Lung cancer is ranked among the top causes of cancer-related deaths. The ratio of adenocarcinoma to squamous cell carcinoma, as well as the ratio of women to men, is still increasing. Human papillomavirus (HPV) has been discovered in lung cancer tissues and blood specimens, particularly in Eastern countries. However, the association between HPV infection and lung adenocarcinoma remains unclear. Methods This population-based cohort study was conducted using data from Taiwan's single-payer national health insurance and cancer registry databases. Data on HPV infection, cancer, sex, age, comorbidities, urbanization, and occupation were collected. The cumulative incidence rates were generated using Kaplan–Meier curves and log-rank tests. COX regression analysis was used to estimate the hazard ratios of factors associated with cancer occurrence. We used data from 2007 and 2015. The cases were matched with sex and age in a 1:2 manner with 939,874 HPV+ and 1,879,748 HPV– individuals, respectively. Results The adjusted hazard ratios [95% confidence interval (CI)] for HPV infection in all lung cancers were 1.539 (1.436–1.649), male lung cancer 1.434 (1.312–1.566), female lung cancer 1.742 (1.557–1.948), squamous cell carcinoma (SCC) 1.092 (0.903–1.320), male SCC 1.092 (0.903–1.320), female SCC 0.949 (0.773–1.164), adenocarcinoma 1.714 (1.572–1.870), male adenocarcinoma 1.646 (1.458–1.858), and female adenocarcinoma 1.646 (1.458–1.858). The highest adjusted hazard ratio for lung cancer was chronic obstructive pulmonary disease (COPD) 1.799 (1.613–2.007), followed by male sex 1.567 (1.451–6.863) and HPV infection. The highest adjusted hazard ratio for adenocarcinoma was HPV infection 1.714 (1.572–1.870), followed by COPD 1.300 (1.102–1.533), and for SCC, male sex 5.645 (4.43–3.37), followed by COPD 2.528 (2.002–3.192). Conclusion Our study showed that HPV infection was associated with the occurrence of adenocarcinoma of the lung in both men and women but was not associated with SCC of the lung.


INTRODUCTION
Cancer is one of the primary causes of death (1) and its etiology involves multiple factors. Prominent factors include infectious agents, such as tuberculosis, human immunodeficiency virus, Epstein-Barr virus, hepatitis viruses, Helicobacter pylori, and human papillomavirus (HPV), which contribute to 20-25% of all globally registered cancer cases (2). Viruses are responsible for 12% of all human cancers, and HPV is responsible for approximately 30% of all cancers related to infectious agents (2). In the United States, 14.1 million new incidences of HPV infections were reported in 2008 and increased to 24.4 million by 2018 (3,4). HPV accounts for the majority of sexually transmitted disease cases in the United States (3). In summary, approximately 79 million women between the ages of 14 and 59 years contracted HPV, with the highest prevalence being in women between the ages of 20 and 24 years (5). HPV infection has been reported to correlate with head and neck, esophageal, and anogenital cancer (6)(7)(8)(9)(10).
Among all cancers, lung cancer is the most lethal in Taiwan (1). In addition to smoking, passive smoking, indoor fume pollution, and genetics, HPV has been proposed to be one of the etiologies of lung cancer (11)(12)(13). HPV was detected in lung squamous cell carcinoma in the pathological form of condylomas in 1979 (14). The metaanalyses conducted by Tsyganov et al. and Karnosky et al. found a correlation between HPV and lung cancer, as people with HPV had 5.4 and 4.7 times greater likelihood of developing lung cancer compared to people without. In addition, detection variation was shown based on geographic areas and HPV subtypes, with people living in Asia having a greater likelihood of having HPV. Currently, studies on the molecular mechanisms underlying HPV-related lung cancer are in progress (15,16). However, to date, no large-scale population-based studies have been conducted to demonstrate the association between HPV and various histological subtypes of lung cancer.
Although an association between HPV infection and the occurrence of Taiwanese lung cancer has been reported in a small number of clinical tumor tissues, no epidemiological data has been found to support the findings (17)(18)(19). Previously, the 2001-2004 Nationwide Population-Based Cohort Study showed a correlation between HPV and lung cancer in female individuals, whereas a correlation between lung cancer and chronic obstructive pulmonary disease (COPD) was found in male individuals (20). Furthermore, data from the National Cancer Registry Database, which includes cancer subtypes, support the association between HPV and lung cancer in cancer diagnosis. Additionally, lung adenocarcinoma was found to differ from squamous cell carcinoma (SCC) in terms of etiology, distribution, genetics, therapeutic medications, and survival (21). In this study, we investigated whether there was a relationship between HPV-related lung cancer and sex and lung cancer histological type.

Data Source
This study analyzed datasets managed by the Health and Welfare Sciences Center, Ministry of Health, Taiwan. Datasets included 2006-2015 ambulatory care expenditures by visits, inpatient expenditures by admissions, details of ambulatory care orders, details of inpatient orders, registry for beneficiaries, cause of death data, and Taiwan Cancer Registry-Long/Short Form. Each patient was identified with a Taiwanese national identification number that was coded to an encrypted identification number to protect personal data. The first database was the 2007-2015 National Health Insurance Research Dataset, which included 26 million unique identification numbers managed by the Health and Welfare Sciences Center. The second database was the 1979-2015 Taiwan Cancer Registry managed by the Health Promotion Administration, Ministry of Health and Welfare, Taiwan.
Stratifications only included cancer subtypes of squamous cell carcinoma and adenocarcinoma. Other subtypes, such as nonsmall cell carcinoma type unknown, large cell, and small cell carcinoma, were excluded because the disease incidence numbers were not large enough to be included in the statistics.

Study Design and Ethical Considerations
For the population-based cohort study design, the index date was identified as the starting point, and the incidence of cancer was from the index date to 31 December 2015. The male population was also excluded from the cervical cancer statistics. The study protocol was reviewed and approved by the institutional review board of Chung Shan University Hospital (IRB CS13168). All patient data from the National Health Insurance Research Dataset were de-identified; therefore, the requirement for signed informed consent was waived.
Statistical Analysis SAS 9.4 (SAS Institute Inc., Cary, North Carolina) was used for statistical analyses in this study. For the design of the cohort study, the index date was identified as the starting point and the incidences of cancer from the index date to 31 December 2015 were included. Patients with HPV before 2008 or those with cancer before the diagnosis of HPV were excluded. The noninfected cases were matched with HPV-infected cases according to age, sex, and index date with structured query language (SQL) in a 2:1 manner. Patients with cancer onset before the index date were excluded. In this study, comparisons were made between HPV-infected and non-HPV-infected individuals. Demographic data were compared using the chi-square test for categorical data and Student's t-test for numerical data. The relative risk of cancer was calculated using a multiple Cox regression model. The cumulative incidence rates of cancer were generated using Kaplan-Meier curves and tested by the long-tank test.

RESULTS
During the period of 2007-2017, 1,103,771 people tested positive for HPV, yielding an incidence rate of 4.2% (Figure 1). Among those with HPV, 141,001 individuals who contracted HPV before 2008 were excluded from the study. Similarly, 22,896 individuals diagnosed with lung cancer before HPV infection were excluded from this study. Subsequently, the remaining 939,873 individuals with HPV were included in the study cohort, and we included double the number of individuals without HPV who were matched on sex, age, and index date (control group). However, identifying HPV infection was challenging because most HPV types that developed viral warts were coded 078.10 and 078.19, and the latter accounted for 91.04% ( Table 1). The age and sex of HPV-positive and HPV-negative individuals were perfectly matched (p = 1.0,000) ( Table 2). A positive correlation was found between the degree of urbanization and incidence rate (p < 0.0001). Officers and laborers had more HPV infections than farmers or unemployed workers did (p < 0.0001). As for comorbidities, the incidence of HPV infection was positively correlated with hyperlipidemia, ischemic heart disease, hypertension, peptic ulcer, gastrointestinal bleeding, chronic kidney disease, and COPD but not with diabetes mellitus.
Human papillomavirus was associated with an increased risk of lung, cervical, and head and neck cancer ( Table 3). The hazard ratio (HR) of cervical cancer between individuals with and without HPV was 2.225 (2.008-2.464). The incidence rate of cervical cancer was 41.65/100,000 person-years for individuals with HPV compared to 18.89/100,000 person-years for individuals without HPV. Similarly, the HR of head and neck cancer was 1.527 (1.436, 1.649, p < 0.01), with an incidence rate of 21.86/100,000 person-years in the HPV positive individuals and an incidence rate of 14.32/100,000 person-years in HPVnegative individuals. The HR for indiscriminate lung cancer was 1.539 (1.436-1.649, p < 0.001), with an incidence rate of 39.44/100,000 person-years for individuals with HPV and an incidence rate of 25.54/100,000 person-years for individuals without HPV. Among all lung cancer individuals, men had an HR of 1.539 (1.439-1.649, p < 0.001) between HPVpositive and HPV-negative individuals, with an incidence rate of 47.47/100,000 person-years for HPV-positive individuals and an incidence rate of 33.40/100,000 person-years for HPV-negative individuals. In comparison, women had an HR of 1.742 (1.557-1.948), with an incidence of 31.57/100,000 person-years for HPV-positive individuals and an incidence rate of 17.83/100,000 person-years for HPV-negative individuals. The incidence rate was higher in HPV-positive than in HPV-negative cervical cancer (p < 0.001) (Supplementary Figure 1), head and neck cancer (Supplementary Figure 2), and lung cancer (Figure 2A).
The results of subgroup analysis indicated that the incidence rate of adenocarcinoma was 3.43 times higher than that of squamous cell carcinoma in the HPV-negative population and 5.87 times higher in the HPV-positive population ( Table 3). For squamous cell carcinoma, the incidence rate was lower in the HPV-positive population (7.54/100,000 person-years) than in the HPV-negative population (7.86/100,000 person-years) among men ( Table 3). In comparison, the incidence rate was higher in the HPV-positive population (1.55/100,000 person-years) than in the HPV-negative population (0.99/100,000 person-years) among women. The HR of squamous cell carcinoma of the lung was 0.949 (0.773, 1.164, p = 0.6,886) in the male population and 1.540 (0.945-2.509, p = 0.0,643) in the female population, which was almost significant (Figures 2B,C). Next, the HPV incidence rate of squamous cell carcinoma was stratified by age ≥ 60 years. The incidence rate in the HPV-positive population (1.09/100,000 person-years) was higher than that in the HPVnegative population (0.94/100,000 person-years) for the group aged < 60 years. A significant outcome was found among the group above the age of 60 years, comparing the incidence rate of the HPV-positive group (32.06/100,000 person-years) with the HPV-negative group (31.74/100,000 person-years). The HR for squamous cell carcinoma was 1.151 (0.762-1.739) for the population below the age of 60 and 0.984 (0.796-1.215) years for the population above the age of 60 years.
The incidence of adenocarcinoma was also compared between HPV-positive and HPV-negative individuals. Stratified One person could have multiple codes. The total number of persons > recruited. ICD-9 795.X and 796.X were rarely used in Taiwan.
by sex, the incidence rate in the HPV-positive population (26.71/100,000 person-years) was higher than that in the HPVnegative population ( The factors associated with lung cancer were determined by COX regression analysis, demonstrated with the adjusted hazard ratio (aHR), 95% confidence interval (CI), and p-value (p). The contributing factors of lung cancer included HPV infection, sex, age, occupation, and other comorbidities ( Table 4). Individuals with HPV had a higher risk of lung cancer compared to individuals without HPV (aHR = 1.539, 1.436-1.649), p < 0.0,001). Similarly, male individuals had a higher risk of lung cancer compared to their female counterparts (aHR = 1.567, 1.459-1.683, p < 0.0,001). The aHR of lung cancer was also proportional to the stratified age. With age 20-40 years as a reference, an aHR of 0.046 (0.017-0.124, p < 0.0,001) was observed among individuals <20 years of age, an aHR of 12.679 (10.415-15.434, p < 0.0,001) was observed among individuals aged 40 and 60 years, an aHR of 45.477 (37.294-55.457, p < 0.001) was observed among individuals between the ages of 60 and 80 years, and an aHR of 73.521 (59.192-91.319, p < 0.0,001) was observed among individuals above the age of 80 years. Furthermore, individuals covered by the insurance of the farmers' and fishermen's unions had greater risk (aHR = 1.242, 1.111-1.388, p = 0.0,001) compared to individuals with laborers' union insurance. Among the common comorbidities listed in Table 4, individuals with hypertension (aHR = 1.183, 1.090-1.285, p < 0.0,001), COPD (aHR = 1.799, 1.613-2.007, p < 0.0,001), and peptic ulcer (aHR = 1.178, 1.076-1.290, p = 0.0,004) were at a greater risk of lung cancer than those without the respective comorbidities. However, the results were dissimilar between cancer cell types.
In individuals with squamous cell carcinoma, sex, age, urbanization, insurance coverage, and comorbidities were contributing factors ( Table 4) Significant causes of adenocarcinoma included HPV infection, age, urbanization, insurance coverage, and comorbidities ( Table 4). HPV-infected adenocarcinoma individuals had an aHR of 1.714 (1.572-1.870, p < 0.0,001) compared to non-HPV individuals, which was the highest among all factors, except age. An aHR was also proportional to age. With age 20-40 years as a reference, an aHR of 0.016 was seen in individuals below the age of 20, an aHR of 26 Table S1). The aHR of HPV-positive cervical cancer individuals was 2.225 (2.008-2.464, p < 0.0,001) with HPV-negative counterparts as a reference. Compared to the age of 20-40 years, however, the maximum aHR was found at the age of 40-60 years (aHR = 1.124, 1.001-1.262, p = 0.0,481), unlike the case in all lung cancer. In terms of insurance coverage, except for those covered by the government that had a lower aHR of 0.720 (0.565-0.917, p = 0.0,077), cervical cancer individuals covered by all other insurance plans had a relatively higher aHR. Furthermore, renal failure was the main contributor to comorbidities, with an aHR of 2.725 (1.235-6.009, p = 0.0,130).

DISCUSSION
In this study, we utilized a nation-based cohort study to investigate the important etiological role of HPV in lung Frontiers in Medicine | www.frontiersin.org adenocarcinoma. Our study found that HPV was associated with adenocarcinoma of the lung in both sexes in Taiwan, with a stronger association in the female population, whereas SCC occurred more frequently in men with COPD.
Smoking is a well-established leading cause of lung cancer. However, approximately 50% of the global male population and 10% of the global female population are reported to possess smoking behaviors (22). The smoking rate in the Asian female population was 2%. The smoking rate in Taiwan has been surveyed since 1990 with a continual downward trend. According to data from 1990 to 2018, the smoking rate in Taiwan decreased from 59.4% to 23.4% in men and 3.8% to 2.4% in women (22). The leading cause of death before 2011 in the male population was hepatoma, followed by lung cancer, and lung cancer was the leading cause of death after 2011 (1). The standard lung cancer incidence rates were 44.2/10 6  According to a study from the Taiwan Cancer Registry from 1996 to 2008, the male-to-female ratio of the lung cancer incidence rate decreased due to the progressive discussion surrounding sex issues. However, the annual incidence rate of adenocarcinoma has increased in both sexes, with a greater increase in the female population. In contrast, the annual incidence rate of squamous cell carcinoma decreases in the female population and remains constant in the male population (21). Despite the drop in the smoking rate, the incidence rate of lung cancer increased from 45.04/10 6 people in 2002 to 49.86/10 6 people in 2011 (22). Notably, the incidence rate increased more rapidly in the female population than in the male population, with an annual 0.2% change in the male population and a 2.0% change in the female population. Overall, the incidence of adenocarcinoma increased from 43.45% to 64.89% between 2002 and 2014 (22).
Our results were temporally different from those of previous studies, in which our study included data from 2006 to 2015 compared to previous studies from 2001 to 2004 (20). During this period, the incidence rate of lung cancer and death rate continued to increase despite the decrease in the smoking rate, suggesting an increase in the incidence rate of non-smokers. The incidence rate of lung cancer in men was lower in this study than in a previous study (20). The latter also supports an association between lung cancer and HPV. Nonetheless, this study indicated an increase in the lung cancer rate among the non-smoking (COPD) population for both men and women, with women having a more rapid increase than men. In terms of cancer cell type, the incidence rate of squamous cell carcinoma decreased in this study, but the incidence rates of adenocarcinoma and lung cancer in the nonsmoking population increased, suggesting an association between HPV infection and lung adenocarcinoma in the non-smoking population. Contrary to a previous study, HPV infection was also associated with adenocarcinoma in the male population, which was supported by an increase in the incidence rate of adenocarcinoma in the non-smoking male population.
Human papillomavirus is prevalent worldwide, accounting for approximately 10% of the global population (23). There are 170 documented types of HPV (7) with additional subtypes. HPV types in most oncologic cell studies, mainly from investigating cervical cancer instead of lung cancer, include 12 viral types closely related to malignant neoplasms (HPV- 16, 18, 31, 33, 35,  The most prevalent HPV types are low-risk (LR)-HPVs-6 and 11, which cause 90% of genital and oral pharyngeal papillomas (6). In most of the studies, samples are collected from semen or vaginal swabs. In our studies, verruca vulgaris 1, 2, 4, 28, and 30, which resulted in common wart formation, were the most frequently investigated types of viruses. The treatment of such viral strains is estimated to cost 1.7 billion a year in the United States. Among the 86 tissue studies of lung cancer and high-risk HPV, including types 16, 18, 31, and 33 and types 16 and 18, were present in 58 articles (15). However, when HPV-6 was detected, HPV-16 and 18 were also detected in 74% (66/89) of the lung cancer tissues. Similarly, they were detected in 49% (21/43) of the lung cancer tissues when HPV-11 was detected. Some negative testing results may be due to the fewer subtypes tested in the studies. However, HPV-6 and 11 are predominant in lung adenocarcinoma; however, subtype and 1, 2, 4, 28, and 30 have rarely been studied in lung cancer tissue. The prevalence of HPV has also been suggested to be associated with the prevalence of cervical cancer. The prevalence rates of HPV differ across continents, with Asia having the highest prevalence rate, followed by Latin America, Australia, Europe, and North America (15,16). The incidence of lung cancer varies geographically. EGFR-mutant lung cancer, for instance, is more prevalent in Asia than in North America. Hence, the study of HPV-contracted lung cancer may be more important in areas with dense populations. This also implies that adenocarcinomas in Taiwan may be associated with HPV infection.
Interestingly, the localization of the virus does not fixate on the epithelium. It may travel beyond superficial layers to the semen, amniotic fluid, placenta, umbilical cord, fetus, breast milk, and blood. Evidence from about 10 studies revealed multiple types of HPV, including 6,11,16,18, and other 10s in blood monocytes and circulatory tumor cells of normal, cervical, and lung cancer subjects (2,24). HPV invades through epidermal defects or potentially through blood circulation and anchors to the basal membrane. However, asymptomatic HPV due to insufficient HPV concentration may escape clinical detection and lead to the development of cancer.
We suggest that squamous cell carcinoma is infected via direct contact with HPV, whereas adenocarcinoma develops through transmission in the bloodstream. Squamous cell carcinoma differs from adenocarcinoma in various aspects, yet HPV has been detected in equal amount in tissue studies: 25.8% of squamous cell carcinoma and 21% of adenocarcinomas are located in tissues, whereas 15.7% of squamous cell carcinomas and 25.4% of adenocarcinomas are located in the blood (15). Due to its origin in the bloodstream, adenocarcinoma has multiple primaries, with EGFR mutations associated with HPV (25,26). The amount of HPV present in the bloodstream that develops into lung cancer carcinoma remains unknown. Squamous cell carcinoma is caused by recurrent respiratory papillomatosis (RRP), which is a direct HPV infection of the respiratory system consisting of type 11 malignant transformations and both types 6 and 11 in other cancer cell types (27,28). From cohort studies, the incidence of lung involvement in RRP was estimated to be 3.3%. The incidence of cancer involving the lung was 5-16%, mostly type 11. In contrast to most molecular mechanisms in other cancers, HPV types 16 and 18 are unnecessary for oncogenesis. Further research is required to determine the mechanism of oncogenesis other than HPV types 16 and 18. As the viral load is much lower in lung cancer than in cervical cancer, they must share different molecular mechanisms of oncogenesis and virus subtypes. In cervical cancer, different survival rates were associated with different physical statuses of HPV. For instance, individuals with the episomal form of HPV have a higher survival rate than those with the integrated form (29). This is due to the integration of HPV vectors into the host genome to disrupt gene function, notably the tumor suppressor gene, as well as the expression of E6 and E7. Oncoproteins E6 and E7 are frequently transcribed after the transformation of HPV DNA into the cellular genome in these cancers. The virus-borne E6 and E7 oncoproteins then bind to the tumor suppressor genes p53 and Rb in the host, thus manipulating the cell cycle by disrupting the normal function of the tumor suppressor genes. In our previous studies, we demonstrated that HPV E6 involved in lung tumorigenesis is partially mediated through (1) p53 inactivation (30); (2) upregulation of IL-6, Mcl-1, hTERT, and expression (31,32); (3) reduced p21 via alteration of the p53/DDX3 (dead-box helicase 3) pathway (33); and (4) TIMP-3 inactivation (34). The reduction of p21 facilitates the release of the cyclin A/CDK2 complex (cyclin A/cyclindependent kinase 2), which phosphorylates and destroys the active form of pRb (RB transcriptional corepressor 1). On the contrary, the interaction of E7 and pRb leads to the dissociation of the HDAC/pRb/E2F (histone deacetylase/RB transcriptional corepressor 1/E2F transcriptional factor) complex. As a result, HDAC participates in p16INK4 hypermethylation (cyclindependent kinase inhibitor 2A), which coordinates with the E2F factor to inactivate tumor suppressor genes. Therefore, cell proliferation genes are activated and tumor progression ensues (35). In addition, HDAC can trigger angiogenesis via HIF-1α to induce VEGF and IL-8. The expression of E7 and the inhibition of p53 through the expression of E6 are involved in the regulation of the anti-apoptotic protein Mcl-1, as MCL1 and BCL2 families are apoptosis regulators that increase cell survival by inhibiting apoptosis via the PI3K/Akt-(IL-6)-(IL-17) pathway. Further studies suggested that p53 is regulated not only by E6 and E7 but also by E1 and E2 may also be involved (36). E5, a less-discussed oncoprotein, serves as a viroporin at the host cell membrane. E5 serves as an inhibitor of the degradation of epithelial growth factor that facilitates cell growth, prohibiting MHC major histocompatibility class I from floating to the cell membrane, which may elude the immune defense. E5 also prevents hydroxyl peroxidase-induced apoptosis (36).
Even though oropharyngeal and cervical cancers are recognized as the same oncogenic sexual transmission diseases and both are squamous cell epithelium, the pathways of oncogenesis are still different (23). We propose that the oncogenic process of lung adenocarcinoma is also different, and further studies are required to determine the underlying mechanism.
The data in our study were gathered from the National Insurance Database and Cancer Registry Database. However, details of the small incidence subtype, such as small cell and large cell lung cancers, were not provided in the database due to privacy protection policies. One of the main limitations of this study was the unrecognizable subtype and location of HPV. The data on cervical cancer and head and neck cancer were included as the standard of comparison to augment the validity of this study since the study of the HPV association of the two cancers was well established. Similarly, data supporting the association between diabetes mellitus and HPV were also included as the standard of comparison (22). Notably, the HPV incidence rate in our study was 39.44/100,000 person-years for men and 25.54/100,000 person-years for women. This aligned with data from the Cancer Registry Database in 2014, in which the standardized HPV incidence rates were 43.3/100,000 personyears for men and 28.3/100,000 person-years for women.
Although the presence of HPV in lung cancer is a possibility, large-scale full virus subtype tissue screening is required to clarify the role of HPV in various subtypes of lung cancer, such as small cell lung cancer and neuroendocrine cell lung cancers. Furthermore, a study of the virology types of lung cancer to investigate the subtype of HPV in different lung cancers is needed, since most current studies only focus on types 16 and 18. Hence, more pathophysiological studies on other types of HPV are required.
To increase the validity of the study, further population studies for small cell lung cancer are required by collecting a larger database of small and large cell cancers owing to the limited dataset on small cell lung cancer. Another limitation is the distinction between high-and low-risk HPV, as this study included mostly viral warts, unspecific plantar warts, and condyloma accumulates. In addition, the HPV-infected population in this study may be underestimated because individuals with subclinical HPV infection are unlikely to visit clinics, which is a natural limitation of the study population.
Further study designs may include investigation of HPV distribution within the bodies of HPV-infected individuals. The study of the oncogenic mechanisms of adenocarcinoma and squamous cell carcinoma of the lung is also needed because of the different pathways across different HPV and cancer subtypes. For instance, the incidence of HPV-infected lung cancer continues to rise despite the decline of CIS in cervical precancer incidence from 2008 to 2012 for women aged 18-20 years in the United States due to effective vaccination (37). Vaccination effectively decreased HPV types 16 and 18 and genital warts in female individuals aged 13-19 years, while the prevalence trends for anogenital warts varied by age. Nevertheless, overall survival rates of HPV-infected cervical, oral, and lung cancers were better in the past, and more emphasis should be placed on the study of HPV-infected cervical, oral, and lung cancers (29,38,39). A high prevalence of EGFR and ALK mutations was found in the lung adenocarcinoma of non-smokers with better survivals, which may benefit from the invention of suitable target therapies, and suggested different oncogenic mechanisms of the cancer itself that call for further investigations (40).
This is the first epidemiological study demonstrating that HPV infection was associated with the occurrence of lung adenocarcinoma, whereas squamous cell carcinoma was mostly associated with COPD. Further validation is warranted to delineate whether there is a causal relationship between HPV infection and lung adenocarcinoma, or simply coexistence.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Chung Shan Medical University Hospital IRB. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.