- 1Department of Epidemiology and Health Statistics, West China School of Public Health and West China Fourth Hospital, Sichuan University, Chengdu, China
- 2Sichuan Cancer Hospital & Institute, Sichuan Cancer Center, School of Medicine, University of Electronic Science and Technology of China, Chengdu, China
- 3The Healthcare Department, Chengdu Women's and Children's Central Hospital, School of Medicine, University of Electronic Science and Technology of China, Chengdu, China
Introduction: Since 2006, the human papillomavirus (HPV) vaccine has been recommended for females aged 9–26 years in the United States. Aiming to evaluate the early effect of the HPV vaccine on cervical cancer, this study assessed the incidence of cervical cancer by age and histology before and after the introduction of HPV vaccination.
Methods: Data on cervical cancer incidence from 1975–2019 were extracted from the Surveillance, Epidemiology, and End Results Program. Joinpoint regression was used to determine temporal trends over time. Future cervical cancer incidence (2015–2039) was projected using Bayesian age-period-cohort analysis. Age-period-cohort (APC) models were created to evaluate age, period, and cohort effects.
Results: For overall cervical cancer and cervical squamous cell carcinoma (SCC), incidence rate showed decreasing trends (–0.7%, and –1.0% annually, respectively), whereas cervical adenocarcinoma (AC) incidence continuously increased (2.6% annually). The incidence trends for AC were stable in the 20–24 and 25–29-year age groups, whereas there was an increasing trend in older age groups. Similarly, the projected trend for AC in females aged 20–30 years exhibited a decline, whereas an increase was predicted in the 31–40–year age group, especially in the 35–44 year age group. The birth cohort and period effects in SCC and AC were extracted from APC models.
Discussion: During the period of 1975–2019, the incidence of cervical AC remained almost unchanged in the age groups receiving HPV vaccines while increased in the age groups not receiving HPV vaccines. The birth cohort effects of SCC and AC of the cervix provided evidence supporting the effectiveness of the HPV vaccine in preventing cervical cancer.
Introduction
Cancer of the uterine cervix ranks as the fourth most common cancer worldwide among females accounting for 604 127 cases and 341 831 deaths in 2020 (1). It is estimated that up to 14,100 and 4,280 new invasive cervical cancer cases and deaths, respectively, will be diagnosed in the United States (US) in 2022 (2). The primary etiologic factor for cervical cancer is persistent human papillomavirus (HPV) infection (3). Bivalent/quadrivalent/9-valent HPV vaccines have been recommended for females aged 9–26 years by the US Advisory Committee on Immunization Practices since 2009 (4), 2006 (5), and 2015 (6), respectively. The bivalent vaccine (2vHPV) and quadrivalent vaccine (4vHPV) are directed against two oncogenic types (HPV 16 and 18) (4), which are responsible for ~70% of cervical cancers (7). Quadrivalent vaccine (4vHPV) covers other two low-risk HPV genotypes (HPV 6 and 11), which account for 90% of anogenital warts (8). The 9-valent vaccine (9vHPV) covers the other five oncogenic types (HPV 31, 33, 45, 52, and 58) that are responsible for another 20% of cervical cancers (7). In the US, the uptake of ≥1 dose of HPV vaccines increased from 37.2 to 69.9% from 2008 to 2018 among adolescent females aged 13–17 years, compared to 20.7–52.8% from 2010 to 2018 among females aged 19–26 years (9–11). Meanwhile, the proportion of female teens who completed the full series of HPV vaccination increased from 17.9% in 2008 to 53.7% in 2018 (9, 10).
Some studies investigated the vaccine's impact on HPV-associated health outcomes (7, 12–22); of these, the majority focused on the effect in reducing the prevalence of the HPV infection and high-grade cervical lesions (12–22). To date, only one study assessed trends in cervical cancer incidence after the introduction of HPV vaccination among young females in the US (7). However, that study did not comprehensively examine the birth cohort effects before and after initiation of HPV vaccination decades ago, nor did it predict future trends of cervical cancer incidence. Among sexually active females aged 14–24 years who participated in the National Health and Nutrition Examination Survey (NHANES), the rate of 4vHPv (types 6/11/16/18) prevalence was reduced compared to that between 2003 and 2006 (15). Previous modeling suggested that the full effect of HPV vaccination in this age group would not be observed until 2025, but half of the benefit would be seen in England by 2019, where the HPV vaccine was first introduced in 2008 (23). The association between quadrivalent/bivalent HPV vaccination and the risk of invasive cervical cancer was assessed in 2020 (24) and 2021(25), respectively. Data used in the studies was retrieved from nationwide Swedish demographic and health registers from 2006 to 2017 (24) and population-based cancer registry in England from 2006 to 2019 (25). The HPV vaccine was introduced in the US 13 years ago, and although it may be premature to assess its full benefits, its early impact can be studied by considering the incidence trends of cervical cancer. To better evaluate the early impact of the HPV vaccine on cervical cancer, we described temporal trends of cervical cancer incidence by age and histology, and built age-period-cohort (APC) models to estimate period, and birth cohort effects, before and after the introduction of HPV vaccination among females aged 20–44 years, using data from the Surveillance, Epidemiology, and End Results (SEER) Program. The period effects referred to policy and environmental factors that impact the populations of all ages. The cohort effects were due to historical factors on certain populations born in the same era. Further, we also projected the cervical cancer incidence by age and histology to 2,039 based on the HPV vaccine introduction by performing Bayesian age-period-cohort (BAPC) models that incorporate age, calendar period, and birth cohort effects.
Materials and methods
Data sources
The data of cervical cancer incidence for 1975–2019 was extracted from the Surveillance, Epidemiology, and End Results (SEER) Program: SEER 9 registries for 1975–1991, which cover 9.4% of the total US population (26), SEER 13 registries for 1992–1999, which cover 13.4% of the US population (26), and SEER 22 registries for 2000–2019, which covered 47.9% of the US population (27). The SEER Program is an authoritative source for cancer statistics in the United States as it is supported by the Surveillance Research Program (SRP) from the NCI's Division of Cancer Control and Population Sciences (DCCPS). The SEER data used here was released in November 2020 and November 2021 [SEER*Stat Database: Incidence—SEER 9, Nov 2020 Sub (1975–2018), SEER 13, Nov 2020 Sub (1992–2018), SEER 22, Nov 2021 Sub (2000–2019), Single Ages to 85+]. Cervical cancer data was divided into 5 age groups (20–24, 25–29, 30–34, 35–39, and 40–44 years) and was identified as site record ‘Cervix Uteri’ according to the International Classification of Diseases for Oncology, Third Edition (ICD-O-3), with histology codes of 8000-9992. Cervical cancer cases are classified into four histologic subtypes: squamous cell carcinoma (SCC) (ICD-O-3 8050-8130), adenocarcinoma (AC) (ICD-O-3 8140-8147, 8160-8162, 8180-8221, 8250-8506, 8520-8550, 8570-8573, 8940-8941), adenosquamous carcinoma (ICD-O-3 8560-8563), and others (including cervical cancers with unknown histology or unspecified carcinomas). Due to sparse data, cervical adenosquamous carcinoma and other histologic subtypes were excluded from subgroup analysis of this study.
A total of 160 868 newly diagnosed cases of cervical cancer were reported to SEER 9, SEER 13, and SEER 22 programs from 1975 to 2019, with 108 364 (67.36%) cases of SCC and 32 523 (20.22%) of AC. The SCC and AC types accounted for 74.43% and 8.86% of cervical cancers in 1975, whereas they accounted for 64.32% and 24.92% of cases in 2019, respectively.
Description of temporal trends
The incidence rates, for SCC, AC of the cervix, and all histopathologic subtypes combined were calculated from 1975 to 2019 for the 20–44-year-old group. In order to describe the incidence trends by birth cohort group, the majority of females aged 15–64 years were included to analyze considering the few cohort groups of females aged 20–44. Herein, there were ten cohort groups (1948–1952, 1953–1957, 1958–1962, 1963–1967, 1968–1972, 1973–1977, 1978–1982, 1983–1987, 1988–1992, and 1993–1997) in this section. The logarithmic transformation of rates in cervical cancer incidence by histology were used to show the pace of changes of incidence trends.
Joinpoint regression model
The trends in the incidence rates of cervical cancer during 1975–2019 by age group and histological subtype were depicted using joinpoint regression analyses, which provided annual percent changes (APCs) and average annual percent changes (AAPCs) (28). Joinpoint regression models were fitted to identify the joinpoints (year) when annual percentage changes (APCs) changed significantly. Annual percentage changes (APCs) were used to characterize the magnitude and direction of trends. It was calculated as {exp(β)−1} × 100, where the regression coefficient (β) was estimated by fitting a least squares regression line to a logarithmic transformation of incidence rates, using the calendar year as an independent variable. Joinpoint regression uses least squares regression to fit line segments to the log-scale rates, joined at discrete points that represent statistically significant changes in direction of the trend. Average annual percent change (AAPC) is a summary measure of the trend over a pre-specified fixed interval. It allows us to use a single number to describe the average APCs over a period of multiple years (28). We specified a logarithmic transformation of incidence rates and a maximum number of three joinpoints to avoid capturing unstable trends due to relatively small numbers of cases in some age groups.
Age-period-cohort model
Age-period-cohort (APC) models were built to evaluate age, period, and birth cohort effects for cervical cancer based on five 5-year age groups (20–24, 25–29, 30–34, 35–39, and 40–44 years old) and nine 5-year time periods (1975–1979, 1980–1984, 1985–1989, 1990–1994, 1995–1999, 2000–2004, 2005–2009, 2010–2014, and 2015–2019). The age effects are the changes in disease rate with age, independent of birth cohort effects and period effects, on account of risk factors linked to maturation. The period effects are caused by factors impacting all subjects during a particular period, such as the implementation of a cervical cancer screening policy. The cohort effects refer to historical factors on certain populations born in the same era, such as the HPV vaccine initiation for younger age groups which was not available for older generations. Assuming to use a log-linear Poisson regression model, the APC model is:
where yij represents the number of cases, and nij the population. The parameter αi represents the age group i effect, βj the time period j effect, and γk the birth cohort k effect. Due to the significant linear correlation among age, period, and birth cohort [period(P) − age(A) = cohort(C)], the APC model presented the non-identifiability issue, i.e., the independent effects of age, period, and cohort cannot be evaluated simultaneously. However, age and birth cohort effects could be estimated by constraining the slope of period effect to be 0 (βj = 0) (29). To verify the cohort effects, we performed sensitivity analyses under different assumptions regarding the slope of the period effect, βp = – 0.01 or + 0.01, which indicated that the period slope decreased or increased, respectively. Moreover, goodness of fit of the APC models was examined using residual deviance statistics by Clayton and Schiffler approach (30–32). The linear time variations are referred to as a parameter “drift” which can be partitioned into any linear combination of period and cohort effects. Therefore, the overall sub-models of the APC model by Clayton and Schiffler are as follows: ① age-only model (A), ② age-drift model, ③ age-cohort model (AC), ④ age-period model (AP), and ⑤ age-period-cohort model. Each sub-model fit was assessed by comparing each iterative model with sequentially adding period and cohort effects to the primary model of age alone to determine whether these added parameters significantly improved model fit. Significant differences were tested in residual deviance of each pairwise comparison using chi-squared tests.
Bayesian age-period-cohort model
Future age-specific cervical cancer incidence and cases from 2015 to 2039 by histology were projected using the Bayesian age-period-cohort analysis (BAPC) with integrated nested Laplace approximations (INLA) (33), which was proved to be the only method to achieve reasonable projections (34). BAPC models involve no parametric assumptions, expecting that all parameters are from appropriate prior distributions. Commonly, smoothing priors, i.e., first-order random walk (RW1) and second-order random walk (RW2), are used for age, period, and birth cohort effects when the assumption that effects adjacent in time might be similar (35). The RW2 prior assumes independent mean-zero normal distributions on the second differences of all time effects (34) and is suitable for a linear time trend, whereas the RW1 prior grants using constant extrapolation (33). Based on the APC model, with the addition of unstructured heterogeneity zij, the Bayesian APC model can be expressed as
For zij, Gaussian random effects can be added to the linear predictor (33). The hyperprior of gamma distribution k ~ G (a, b) was assumed for all precision parameters with shape parameter a = 1 and rate parameter b = 0.005.
Calculations were restricted to individuals targeted by the US HPV vaccination program, including cases in 20–24 to 40–44 years age groups for each year between 1975 and 2019 and excluding females younger than 20 years with extremely low cervical cancer incidence. Adults ≥34 years old in 2017 had low vaccine uptake (2.7%) (36), so we conservatively assumed no association between vaccination and observed cervical cancer in 35–39 and 40–44 age groups in 2015–2019. We projected future cervical cancer incidence under the scenario of no HPV vaccination in 35–39 and 40–44 age groups, assuming that the observed contemporary trends of cervical cancer incidence continue (because the model that used current levels of HPV vaccination to estimate their future impact).
Information on cases was exported from SEER*Stat (version 8.4.0) to R (version 3.6.1) for further statistical analyses. Joinpoint Regression analyses were conducted using Joinpoint Desktop Software version 4.9.0.1 (National Cancer Institute, Bethesda, MD, USA). All APC and BAPC analyses were performed using the apc.fit function from the R-package Epi and the R-package BAPC, respectively. The BAPC package was built upon INLA-package (version 20.03.07), to predict future cancer rates and cases within a fully Bayesian inference setting. The national projected population was used for estimated population sizes in 2017 and beyond. All tests were 2-sided at a significance level of α = 0.05, with 95% confidence intervals.
Data used in this research are publicly available and are deidentified, hence this study was exempt from approval of the institutional review board.
Results
A total of 63 304 cervical cancer cases in the 20–44 age groups were newly diagnosed from 1975 to 2019. Among them, the numbers of cervical SCC, cervical AC, cervical adenosquamous carcinoma and other histologic subtypes cases were 42 220 (66.69%), 13 694 (21.63%), 2 456 (3.88%), and 4 934 (7.79%), respectively.
Cases and rates on cervical cancer incidence (1975–2019) for females aged 20 to 44 years by year of diagnosis (i.e., time period) and age group were shown in Tables 1, 2, respectively. The incidence rates for cervical SCC, AC, and all histologic subtypes combined for the 20–44 years age group are presented in Figure 1. Similar to the overall histological incidence rate, the incidence rate of cervical SCC confirmed that cervical cancer generally decreased (overall: −0.7% per year, SCC: −1.0% per year), except for 1980–1994 or 2010–2019. However, the incidence rate of cervical AC continuously increased throughout 1975–2019 (2.6% annually). The trends of age-specific incidence rates over time period for SCC, AC, and all histological types combined according to age group are shown in Figure 2. APCs and AAPCs in the five age groups by histological subtype are presented in Table 3. For overall cervical cancer and SCC of the cervix, the increase during 2010–2019 was largely confined to ages 30–34, 35–39, and 40–44 years (Figures 2A,B). The trends for AC of the cervix were stable in the 20–24 and 25–29 age groups, and there was a distinct increase in older age groups (Figure 2C). The birth cohort trends of age-specific incidence rates for overall cervical cancer, SCC, and cervical AC by age groups are shown in Figure 3, which showed the trends over birth cohort The incidence rates for overall cervical cancer and SCC of the cervix reflected generally decreasing trends by cohort (Figures 3A,B). In contrast, except for a stable trend in cervical AC incidence rates in the 20–24 age group, the incidence visibly increased in the 30–44 age group and slightly increased in the 25–29 age group. It should be noted that in the 20–24 and 25–29-year groups, the 1983–1997 cohorts showed decreasing trends of incidence rates (Figure 3C).
Table 1. Cases and rates on cervical cancer incidence (1975–2019) for females aged 20–44 years by histology and year of diagnosis, United States.
Table 2. Cases and rates on cervical cancer incidence (1975–2019) for females aged 20–44 years by histology and age group, United States.
Figure 1. Temporal trends in cervical cancer incidence rates (1975–2019) among females aged 20 to 44 years by histology, United States. All data are expressed as the rate per 100,000.
Figure 2. Temporal trends in cervical cancer incidence rates (1975–2019) among females aged 20 to 44 years by age and histology, United States. (A) Temporal trends of all histopathologic subtypes combined. (B) Temporal trends of squamous cell carcinoma (SCC). (C) Temporal trends of adenocarcinoma (AC). All data are expressed as the rate per 100,000.
Table 3. Trends on cervical cancer incidence rates (1975–2019) for females aged 20–44 years by histology and age, United States.
Figure 3. Cohort trends in cervical cancer incidence rates (1975–2019) among females aged 20 to 44 years by age and histology, United States. (A) Cohort trends of all histopathologic subtypes combined. (B) Cohort trends of squamous cell carcinoma (SCC). (C) Cohort trends of adenocarcinoma (AC). All data are expressed as the rate per 100,000.
Considering the few cohort groups of females aged 20–44, females aged 15–64 years were included to analyze the incidence trends by birth cohort group (Figure 4). For SCC, AC, and overall cervical cancer, the incidence rates in the 1988–1997 birth cohorts were lower than those in the 1978–1987 birth cohorts. In particular, for SCC, the 1978–1997 cohort groups had a lower incidence than the 1948–1977 cohorts. By contrast, for cervical AC the incidence rates in the 1978–1997 cohorts were comparable to those of birth cohorts of earlier years. However, the log-scale curves showed a significantly lower rate of incidence in the 1993–1997 birth cohort than in other cohorts (Figure 5).
Figure 4. Age trends in cervical cancer incidence rates (1975-2019) among females in 1948-1997 birth cohorts by cohort and histology, United States. (A) Age trends of all histopathologic subtypes combined. (B) Age trends of squamous cell carcinoma (SCC). (C) Age trends of adenocarcinoma (AC). All data are expressed as the rate per 100,000.
Figure 5. The log-scale age trends in cervical cancer incidence rates (1975–2019) among females in 1948–1997 birth cohorts by age and histology, United States. (A) Log-scale age trends of all histopathologic subtypes combined. (B) Log-scale age trends of squamous cell carcinoma (SCC). (C) Log-scale age trends of adenocarcinoma (AC). All data are expressed as the rate per 100,000.
APC models indicated the cohort effects and period effects in overall cervical cancer, SCC, and cervical AC incidences (Table 4). Significant improvements in model fitting were found when drift (i.e., the overall linear trend in cervical cancer incidence), period, and cohort effects were added (P < 0.05). The best-fitted model included all three effects (age, period, cohort effects), demonstrated by the lowest residual deviance, for overall, SCC, and AC (residual deviance: 57.50, 58.70, 38.82). Females born between 1970 and 1985 had slightly lower cervical cancer incidence than females born in earlier years, while females born after 1985 had visibly lower cervical cancer incidence than those born earlier, as shown in Figures 6A,B. For AC of the cervix, females born between 1960 and 1970 had notably higher cervical cancer incidence than those born earlier, as shown in Figure 6C. Finally, females born between 1970 and 1985 and those born after 1985 experienced a slight increase and a sharp decrease in cervical cancer incidence, respectively, compared to cohorts born earlier.
Table 4. Age-period-cohort models for cervical cancer incidence among 20–44-year group females in US, by histological subtype, 1975–2019.
Figure 6. Age, period, and cohort* effects in cervical cancer incidence rates among females aged 20 to 44 years by histology, United States. (A) Effects of all histopathologic subtypes combined. (B) Effects of squamous cell carcinoma (SCC). (C) Effects of adenocarcinoma (AC). * Each plot's horizontal axis is divided into two parts: age, ranging from 20–44 years (left), and calendar time, ranging from 1935–2020 (right). Each plot contains two vertical axes: cervical cancer incidence per 100,000 person-years (left) and rate ratios (right), and three sets of curves: age effects, with corresponding 95% confidence intervals (left); cohort effects, with corresponding 95% confidence intervals (middle); and period effects, with corresponding 95% confidence intervals (right). In cohort effects, the slope of the period was constrained to be 0 to analyze the cohort effect, and the empty red circle represents the reference cohort of 1970 (median date of birth among cervical cancer cases).
Moreover, sensitivity analyses for APC models using different slopes of the period showed that regardless of the period effects, the trends of cohort effects remained the same for both SCC and AC (Figure 7), similar to those in Figure 6.
Figure 7. Sensitivity analyses of age-period-cohort model in cervical cancer incidence rates among females aged 20 to 44 years by histology, United States. (A) Results of all histopathologic subtypes combined. (B) Results of squamous cell carcinoma (SCC). (C) Results of adenocarcinoma (AC).The slope of the period effect was assumed to be β = 0 (solid lines), −0.01 (dotted lines), and +0.01 (dashed lines) to analyze the cohort effect. The scale on the left axis is for age effect. The scale for the cohort and period effects is on the right axis. The empty circle represents the reference cohort of 1970 (median date of birth among cervical cancer cases).
The predicted new cases of cervical cancer for 2035–2039 among 20–44 years old were 19 801, 10 029, 11 068 for overall, SCC and AC, respectively. Predicted incidence cases and rates for SCC, AC, and all histologic subtypes combined by age group between the ages of 20 and 44 from 2015 to 2039 are presented in Tables 5, 6. According to our projection, incidence rates of overall cervical cancer and SCC of the cervix might continue to decrease throughout 2015–2039 in each age group. Of note, the declining trend in the 35–39 and 40–44 age groups gradually slowed down (Figures 8A,B). As for AC of the cervix, the incidence in females aged 20–29 years took on a declining trend. In contrast, the predicted incidence among those aged 30–44 years seems to increase, especially for those in 35–39 and 40–44 age groups (Figure 8C).
Table 5. Projected number of new cases of cervical cancer among females aged 20–44 years by age, 2015–2039, United States.
Table 6. Projected incidence rates (per 100,000 person-years) of cervical cancer among females aged 20–44 years by age, 2015–2039, United States.
Figure 8. Predicted trends in cervical cancer incidence rates (2015-2039) among females aged 20 to 44 years by age and histology, United States. (A) Predicted trends of all histopathologic subtypes combined. (B) Predicted trends of squamous cell carcinoma (SCC) (C) Predicted trends of adenocarcinoma (AC). All data are expressed as the rate per 100,000. The fan represents the predictive distribution between the 5% and 95% quantiles. Each shaded band represents an additional 10% quantile. The solid line shows the predictive mean, and the white dot represents the observed incidence. The vertical dashed line indicates 2014, when prediction started.
Discussion
In this study, we calculated the temporal and age-period-cohort effects and future trends of cervical cancer incidence in the US, the results of which provide supporting evidence for the effectiveness of the HPV vaccine in preventing cervical cancer, regardless of the histological subtype.
Reasons for the changing trends of cervical cancer incidence might be multi-faceted. The increasing overall pattern could be due to in-utero exposure to maternal diethylstilbestrol (DES) (37), obesity (38), and long-term use of hormonal contraceptives (39). However, the overall increased incidence of cervical cancer could be due to two factors. First, from WHO classification of female genital tumors (5th Edition) (40), the vast majority (>90–95%) of squamous cell carcinoma of the cervix is caused by high-risk HPV genotypes, while HPV-associated adenocarcinoma of the cervix accounts for ~75% of all cervical adenocarcinomas. HPV 16 and 18 together account for 95% of HPV-associated cervical adenocarcinoma (41, 42). As the majority of cervical cancer were caused by HPV infections (3), HPV vaccination might contribute to reducing the incidence. Second, it could be the impact of the current screening policy targeted at cervical cancer (43). Since the 1960s and 1990s, the Papanicolaou (Pap) test and Thinprep cytologic test (TCT) were widely applied, successively (43). Thereafter, Pap/HPV DNA co-testing and primary testing were included in screening guidelines in 2012 and 2018 (44). For the Pap and TCT tests, sensitivity is superior for cervical SCC than for AC (45, 46). Since less is known about the natural history of AC of the cervix, the impact of Pap tests on its incidence is unclear (47). Chronologically increasing trends in the incidence rate of AC in the US from 1975 to 2013 resemble those of other non-screen-detected HPV-associated cancers, such as oropharyngeal and anal cancer (48). Therefore, the incidence trends and cohort effects for cervical AC, rather than SCC, could more likely reflect the impact of HPV vaccines.
The quadrivalent human papillomavirus vaccine was first recommended by the ACIP in the US in 2006 (5). The vaccine is recommended for girls 11–12 years old and can also be administered to girls as young as nine. Catch-up vaccination is recommended for females aged 13–26 years who were not previously vaccinated (5). In the US, the coverage of HPV vaccine initiations (received ≥1 dose) was 25.1, 53.8, 69.9, and 77.1% in 2007, 2012, 2018, and 2020, respectively, among adolescent females aged 13–17. The coverage among females aged 19–26 years increased from 20.7% in 2010 to 52.8% in 2018 (9–11); only 2.7% of adults aged ≥34 in 2017 had been vaccinated (38). Our results revealed that the temporal trends of AC incidence rate have been increasing throughout 1975–2019, where trends in the 20–29 age groups remained relatively stable, whereas the incidence displayed an increasing trend in the 30–34, 35–39, and 40–44 age groups, particularly in the latter two groups. Compared to the 35–44 age groups, the 20–29-year age groups who received HPV vaccines exhibited lower incidences of cervical AC. Since females 35–39, 40–44 years of age in 2019 were 23–27, 28–32 years of age in 2007 and 28–32, 33–37 years of age in 2012, they would not have been within the recommended window for catch-up HPV vaccination when first approved, and either not in the age intervals when the coverage of HPV vaccination reached over 50%. Similarly, according to the projected trends from Bayesian APC analyses, AC incidence of females aged 20–30 will decline, but increase among those aged 31–40, especially in the 35–44 age group.
Specifically, from the results of incidence trends by birth cohorts, regardless of the histological subtype, cervical cancer incidence rates in the 1988–1997 birth cohorts when there was greater coverage of HPV vaccination, were lower than those of the 1978–1987 birth cohorts. For AC of the cervix, considering increasing incidence, the incidence in the 1978–1997 cohorts was comparable to that from the 1948–1977 cohorts. The incidence in the 1993–1997 birth cohort with the highest vaccination coverage displayed a markedly slower growth rate than others. The predicted rates by birth cohort groups from BAPC analyses (Figure 9) showed that regarding SCC of the cervix, the younger the cohort population, the lower the incidence. For AC, after 30 years old, the younger cohorts exhibited higher incidence trends than older ones. Although the reason might be the increasing incidence of AC, the differences of predicted trends among 1998–2012 birth cohort, 1988–1997 birth cohort, and 1948–1987 birth cohort might be misleading. Because, based on the current risk factors, projections obtained through BAPC cannot foresee the changes in risk factors or new interventions in the future (33); thus, BAPC did not account for the potentially increasing coverage of HPV vaccination for 1998–2012 and 1988–1997 birth cohorts from 2015 to 2039.
Figure 9. Predicted age trends in cervical cancer incidence rates among females in 1948–1997 birth cohorts by cohort and histology, United States. (A) Predicted age trends of all histopathologic subtypes combined. (B) Predicted age trends of squamous cell carcinoma (SCC). (C) Predicted age trends of adenocarcinoma (AC). All data are expressed as the rate per 100,000.
The cohort effects for cervical AC also suggested the effectiveness of HPV vaccination. AC incidence in females born after 1985 exhibited a sharp decline and was less than or equal to that of those 26 years old in 2011. These findings were consistent with the cohort effects of cervical intraepithelial neoplasia grades 2 or 3 and AC in situ (collectively, CIN2+) incidence among TennCare-enrolled females in the state of Tennessee, US (22). Furthermore, the cohort effects of cervical SCC displaying a decreasing trend in the 1985–1997 birth cohorts were likely due to generational differences in HPV vaccine eligibility and vaccination behaviors.
National guidelines currently recommend that women undergo Pap testing at 3-year intervals starting at 21 years old and plus HPV DNA testing from the age of 30 years, regardless of HPV vaccination status (49–51). The screening policy has changed several times during the study period: the interval of screening, not recommended for young females aged < 21 years, was prolonged (50–53). In the US, the screening coverage rate declined modestly in recent years (54, 55), especially among the 20–29-year age group, consistent with the trending of overall incidence of cervical cancer and SCC of the cervix. Nevertheless, the screening uptake in the 30–39 age group showed an increasing trend during 2013–2015 (54), which could have contributed to the increasing incidence of overall cervical cancer and SCC of the cervix among the same age groups during 2012–2019. The period effects of overall cervical cancer and SCC were significant, which implied that the influences were caused by changes in the screening policy. Meanwhile, the period trends in AC exhibited a mild increase, and the period effects for AC were significant (P = 0.013). The potential reason could be that the Pap/HPV DNA co-test was sensitive to AC of the cervix.
Although screening with subsequent treatment for precancerous cervical lesions has significantly reduced the incidence of cervical cancer (56), not all females receive screening or follow-up treatment, and cervical cancer remains a considerable cause of mortality (57). Besides, based on our BAPC analyses, the incidence of cervical AC is likely to continue to increase from 2015 to 2039. The future disease burden of AC may be 11 068 cases among females aged 20–44 years in 2035–2039, while it might be 6 617 cases in 2015–2019. The burden of overall cervical cancer in 2035–2039 (19 801 cases) appears to resemble that in 2015–2019 (21 602 cases). Therefore, expanded HPV vaccine coverage may be an effective approach to achieving general reductions in cervical cancer over the coming years. Our projection of incidence rates during the 2015–2019 period was comparable with that from the SEER program, indicating that our predictions could be considered reliable.
Our study has some limitations. The data in our study did not include all individuals in the US. However, by analyzing the Seer 9, Seer 13, and Seer 22 program combined data, which covered nearly 47.9% of the US population, our research has an advantage over other studies that solely used the Seer 9 program (which covered 9.4% of the US population). Besides, we used the latest data from the Seer 22 program, which was released on April 15, 2022. Moreover, as an ecologic study, we were unable to examine individual-level potential confounders and vaccination data, but instead were able to use the age-period-cohort model to explain most confounding factors such as age effects and period effects and provide a methodologically acceptable design for studying the intervention effects of HPV vaccination. Our study was the first to comprehensively examine the birth-cohort effects before and after the HPV vaccination had been introduced in recent decades.
Conclusions
The cervical AC incidence increased from 1975 to 2019; however, this increase involved only those who did not receive the HPV vaccine, while the incidence was stable in the vaccinated. The projected trends showed similar results whereby vaccinated females will exhibit a decline in cervical AC incidence rate, whereas the unvaccinated will display an increased rate. The birth cohort effects of SCC, AC, and all histologic subtypes combined also provide evidence for the effectiveness of the HPV vaccine in preventing cervical cancer.
Data availability statement
Publicly available datasets were analyzed in this study. This data can be found here: the National Cancer Institute's (NCI's) Surveillance, Epidemiology, and End Results (SEER) Program 9+13+22.
Author contributions
JP: conceptualization, methodology [software], data curation, formal analysis, writing—original draft, and writing—review and editing. TS: methodology [software], data curation and analysis, validation, and writing—review. CW: data analysis and writing—review and editing. ML and MX: methodology [software], data analysis, and writing—review. CZ: conceptualization, methodology, formal analysis, validation, project administration, and writing—review and editing. MJ: conceptualization, methodology, formal analysis, validation, supervision, and writing—review and editing. All authors contributed to the article and approved the submitted version.
Funding
This study was supported by funding from the National Key R&D Program of China (Grant no. 2016YFC1302505-2).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. (2021) 71:209–49. doi: 10.3322/caac.21660
2. American Cancer Society. Cancer Facts & Figures 2022. American Cancer Society (2022). Available online at: https://www.cancer.org/content/dam/cancer-org/research/cancer-facts-and-statistics/annual-cancer-facts-and-figures/2022/2022-cancer-facts-and-figures.pdf (accessed April 25, 2022).
3. Sharma M, Bruni L, Diaz M, Castellsagué X, de Sanjosé S, Bosch F. X, et al. Using HPV prevalence to predict cervical cancer incidence. Int J Cancer. (2013) 132:1895–900. doi: 10.1002/ijc.27835
4. CDC. FDA licensure of bivalent human papillomavirus vaccine (HPV2, Cervarix) for use in females and updated HPV vaccination recommendations from the advisory committee on immunization practices (ACIP). Morb Mortality Wkly Rep. (2010) 59:626–9.
5. Markowitz LE, Dunne EF, Saraiya M, Lawson HW, Chesson H, Unger ER. Quadrivalent human papillomavirus vaccine: recommendations of the advisory committee on immunization practices (ACIP). MMWR Recomm Rep. (2007) 56:1–24. doi: 10.1037/e601292007-001
6. Petrosky E, Bocchini JA Jr, Hariri S, Chesson H, Curtis C.R, Saraiya M, et al. Use of 9-valent human papillomavirus (HPV) vaccine: updated HPV vaccination recommendations of the advisory committee on immunization practices. Morb Mortality Wkly Rep. (2015) 64:300–4.
7. Guo FJ, Cofie LE, Berenson AB. Cervical cancer incidence in young u.s. females after human papillomavirus vaccine introduction. Am J Prev Med. (2018) 55:197–204. doi: 10.1016/j.amepre.2018.03.013
8. Garland SM, Kjaer SK, Muñoz N, Block SL, Brown DR, DiNubile MJ, et al. Impact and effectiveness of the quadrivalent human papillomavirus vaccine: a systematic review of 10 years of real-world experience. Clin Infect Dis. (2016) 63:519–27. doi: 10.1093/cid/ciw354
9. CDC. National, state, and local area vaccination coverage among adolescents aged 13-17 years – United States, 2008. Morb Mortality Wkly Rep. (2009) 58:997–1001.
10. Walker TY, Elam-Evans LD, Yankey D, Markowitz LE, Williams CL, Fredua B, et al. National, regional, state, and selected local area vaccination coverage among adolescents aged 13-17 years – United States, 2018. Morb Mortality Wkly Rep. (2019) 68:718–23. doi: 10.15585/mmwr.mm6833a2
11. Lu PJ, Hung MC, Srivastav A, Grohskopf LA, Kobayashi M, Harris AM, et al. Surveillance of vaccination coverage among adult populations – United States, 2018. Morb Mortality Wkly Rep Surveill Summ. (2021) 70:1–26. doi: 10.15585/mmwr.ss7003a1
12. Hariri S, Johnson ML, Bennett NM, Bauer HM, Park IU, Schafer S, et al. Population-based trends in high-grade cervical lesions in the early human papillomavirus vaccine era in the United States. Cancer. (2015) 121:2775–81. doi: 10.1002/cncr.29266
13. Markowitz LE, Hariri S, Lin C, Dunne EF, Steinau M, McQuillan G, et al. 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–93. doi: 10.1093/infdis/jit192
14. Markowitz LE, Liu G, Hariri S, Steinau M, Dunne EF, Unger ER. Prevalence of HPV after introduction of the vaccination program in the United States. Pediatrics. (2016) 137:e20151968. doi: 10.1542/peds.2015-1968
15. Oliver SE, Unger ER, Lewis R, McDaniel D, Gargano JW, Steinau M, et al. Prevalence of human papillomavirus among females after vaccine introduction – National Health and Nutrition Examination Survey, United States, 2003-2014. J Infect Dis. (2017) 216:594–603. doi: 10.1093/infdis/jix244
16. Dunne EF, Naleway A, Smith N, Crane B, Weinmann S, Braxton J, et al. Reduction in human papillomavirus vaccine type prevalence among young women screened for cervical cancer in an integrated US healthcare delivery system in 2007 and 2012-2013. J Infect Dis. (2015) 212:1970–5. doi: 10.1093/infdis/jiv342
17. Kahn JA, Brown DR, Ding L, Widdice LE, Shew ML, Glynn S, et al. Vaccine-type human papillomavirus and evidence of herd protection after vaccine introduction. Pediatrics. (2012) 130:e249–56. doi: 10.1542/peds.2011-3587
18. Kahn JA, Widdice LE, Ding L, Huang B, Brown DR, Franco EL, et al. Substantial decline in vaccine-type human papillomavirus (HPV) among vaccinated young women during the first 8 years after HPV vaccine introduction in a community. Clin Infect Dis. (2016) 63:1281–7. doi: 10.1093/cid/ciw533
19. Benard VB, Castle PE, Jenison SA, Hunt WC, Kim JJ, Cuzick J, et al. Population-based incidence rates of cervical intraepithelial neoplasia in the human papillomavirus vaccine era. JAMA Oncol. (2017) 3:833–7. doi: 10.1001/jamaoncol.2016.3609
20. Flagg EW, Torrone EA, Weinstock H. Ecological association of human papillomavirus vaccination with cervical dysplasia prevalence in the United States, 2007-2014. Am J Public Health. (2016) 106:2211–8. doi: 10.2105/AJPH.2016.303472
21. Berenson AB, Laz TH, Rahman M. Reduction in vaccine-type human papillomavirus prevalence among women in the United States, 2009-2012. J Infect Dis. (2016) 214:1961–4. doi: 10.1093/infdis/jiw515
22. Shing JZ, Beeghly-Fadiel A, Griffin MR, Chang RS, Sudenga SL, Slaughter JC, et al. The impact of the human papillomavirus vaccine on high-grade cervical lesions in urban and rural areas: an age-period-cohort analysis. Cancers. (2021) 13:4215. doi: 10.3390/cancers13164215
23. Cuzick J, Castañón A, Sasieni P. Predicted impact of vaccination against human papillomavirus 16/18 on cancer incidence and cervical abnormalities in women aged 20-29 in the UK. Br J Cancer. (2010) 102:933–9. doi: 10.1038/sj.bjc.6605528
24. Lei J, Ploner A, Elfström KM, Wang J, Roth A, Fang F, et al. HPV vaccination and the risk of invasive cervical cancer. N Engl J Med. (2020) 383:1340–48. doi: 10.1056/NEJMoa1917338
25. Falcaro M, Castañon Castañon A, Ndlela B, Checchi M, Soldan K, Lopez-Bernal J, et al. The effects of the national HPV vaccination programme in England, UK, on cervical cancer and grade 3 cervical intraepithelial neoplasia incidence: a register-based observational study. Lancet. (2021) 398:2084–92. doi: 10.1016/S0140-6736(21)02178-4
26. National Cancer Institute Surveillance Epidemiology, and End Results Program. SEER*Stat Databases: November 2020 Submission. Available online at: https://seer.cancer.gov/data-software/documentation/seerstat/nov2020/ (accessed April 25, 2022).
27. National Cancer Institute Surveillance Epidemiology, and End Results Program. SEER*Stat Databases: November 2021 Submission. Available online at: https://seer.cancer.gov/data-software/documentation/seerstat/nov2021/ (accessed April 25, 2022).
28. Division of Cancer Control and Population Sciences, National Cancer Institute (NIH). Joinpoint Help Manual 4.8.0.1. Available online at: https://surveillance.cancer.gov/joinpoint/Joinpoint_Help_4.8.0.1.pdf (accessed April 25, 2022).
29. Holford TR. Understanding the effects of age, period, and cohort on incidence and mortality rates. Annu Rev Public Health. (1991) 12:425–57. doi: 10.1146/annurev.pu.12.050191.002233
30. Robertson C. Chapter7: age-Period-Cohort Models. In:Zheng T, Boffetta P, Boyle P, editors. Epidemiology and Biostatistics. Lyon: International Prevention Research Institute (2011). p. 127–41.
31. Clayton D, Schifflers E. Models for temporal variation in cancer rates. I: age-period and age-cohort models. Stat Med. (1987) 6:449–67. doi: 10.1002/sim.4780060405
32. Clayton D, Schifflers E. Models for temporal variation in cancer rates. II: age-period-cohort models. Stat Med. (1987) 6: 469–81. doi: 10.1002/sim.4780060406
33. Riebler A, Held L. Projecting the future burden of cancer: Bayesian age-period-cohort analysis with integrated nested Laplace approximations. Biom J. (2017) 59:531–49. doi: 10.1002/bimj.201500263
34. Bray, I. Application of Markov chain Monte Carlo methods to projecting cancer incidence and mortality. J R Stat Soc Ser C Appl Stat. (2002) 51:151–64. doi: 10.1111/1467-9876.00260
35. Knorr-Held L, Rainer E. Projections of lung cancer mortality in West Germany: a case study in Bayesian prediction. Biostatistics. (2001) 2:109–29. doi: 10.1093/biostatistics/2.1.109
36. Zhang Y, Fakhry C, D'Souza G. Projected association of human papillomavirus vaccination with oropharynx cancer incidence in the US, 2020-2045. JAMA Oncol. (2021) 7:e212907. doi: 10.1001/jamaoncol.2021.2907
37. Zhu C, Bassig BA, Zhang Y, Shi K, Boyle P, Li N, et al. A birth cohort analysis of the incidence of adenocarcinoma of the uterine cervix in the USA. Eur J Cancer Prev. (2015) 24:240–44. doi: 10.1097/CEJ.0000000000000062
38. Lacey JV, Swanson CA, Brinton LA, Altekruse SF, Barnes WA, Gravitt PE, et al. Obesity as a potential risk factor for adenocarcinomas and squamous cell carcinomas of the uterine cervix. Cancer. (2003) 98:814–21. doi: 10.1002/cncr.11567
39. Castellsagué X, Díaz M, de Sanjosé S, Muñoz N, Herrero R, Franceschi S, et al. Worldwide human papillomavirus etiology of cervical adenocarcinoma and its cofactors: implications for screening and prevention. J Natl Cancer Inst. (2006) 98:303–15. doi: 10.1093/jnci/djj067
40. WHO Classification of Tumours Editorial Board. Female Genital Tumours. WHO Classification of Tumours. 5th ed. Lyon: International Agency for Research on Cancer (2020) 4:347–69.
41. Pirog EC, Kleter B, Olgac S, Bobkiewicz P, Lindeman J, Quint WG, et al. Prevalence of human papillomavirus DNA in different histological subtypes of cervical adenocarcinoma. Am J Pathol. (2000) 157:1055–62. doi: 10.1016/S0002-9440(10)64619-6
42. Quint KD, de Koning MN, Geraets DT, Quint WG, Pirog EC. Comprehensive analysis of Human papillomavirus and Chlamydia trachomatis in in-situ and invasive cervical adenocarcinoma. Gynecol Oncol. (2009) 114:390–4. doi: 10.1016/j.ygyno.2009.05.013
43. Jemal A, Simard EP, Dorell C, Noone AM, Markowitz LE, Kohler B, et al. Annual Report to the Nation on the Status of Cancer, 1975-2009, featuring the burden and trends in human papillomavirus (HPV)-associated cancers and HPV vaccination coverage levels. J Natl Cancer Inst. (2013) 105:175–201. doi: 10.1093/jnci/djs491
44. US Preventive Services Task Force, Curry SJ, Krist AH, Owens DK, Barry MJ, Caughey AB, et al. Screening for cervical cancer: US Preventive Services Task Force recommendation statement. JAMA. (2018) 320:674–86. doi: 10.1001/jama.2018.10897
45. Castle PE, Kinney WK, Cheung LC, Gage JC, Fetterman B, Poitras NE, et al. Why does cervical cancer occur in a state-of-the-art screening program? Gynecol Oncol. (2017) 146:546–53. doi: 10.1016/j.ygyno.2017.06.003
46. Katki HA, Kinney WK, Fetterman B, Lorey T, Lorey NE, Cheung L, et al. Cervical cancer risk for women undergoing concurrent testing for human papillomavirus and cervical cytology: a population-based study in routine clinical practice. Lancet Oncol. (2011) 12:663–72. doi: 10.1016/S1470-2045(11)70145-0
47. Gradíssimo A, Burk RD. Molecular tests potentially improving HPV screening and genotyping for cervical cancer prevention. Expert Rev Mol Diagn. (2017) 17:379–91. doi: 10.1080/14737159.2017.1293525
48. Ryser MD, Rositch A, Gravitt PE. Modeling of US human papillomavirus (HPV) seroprevalence by age and sexual behavior indicates an increasing trend of HPV infection following the sexual revolution. J Infect Dis. (2017) 216:604–11. doi: 10.1093/infdis/jix333
49. US Preventive Services Task Force. The Guide to Clinical Preventive Services 2014. Available online at: www.ahrq.gov/sites/default/files/wysiwyg/professionals/clinicians-providers/guidelines-recommendations/guide/cpsguide.pdf (accessed April 25, 2022).
50. Saslow D, Solomon D, Lawson HW, Killackey M, Kulasingam SL, Cain J, et al. American Cancer Society, American Society for Colposcopy and Cervical Pathology, and American Society for Clinical Pathology screening guidelines for the prevention and early detection of cervical cancer. CA Cancer J Clin. (2012) 62:147–72. doi: 10.3322/caac.21139
51. Committee on Practice Bulletins-Gynecology. ACOG Practice Bulletin Number 131: screening for cervical cancer. Obstet Gynecol. (2012) 120:1222–38. doi: 10.1097/AOG.0b013e318277c92a
52. Committee on Practice Bulletins-Gynecology. ACOG Practice Bulletin no. 109: cervical cytology screening. Obstet Gynecol. (2009) 114:1409–20. doi: 10.1097/AOG.0b013e3181c6f8a4
53. Whitlock EP, Vesco KK, Eder M, Lin JS, Senger CA, Burda BU. Liquid-based cytology and human papillomavirus testing to screen for cervical cancer: a systematic review for the U.S. Preventive Services Task Force. Ann Intern Med. (2011) 155:687–97. doi: 10.7326/0003-4819-155-10-201111150-00376
54. Watson M, Benard V, King J, Crawford A, Saraiya M. National assessment of HPV and Pap tests: changes in cervical cancer screening, National Health Interview Survey. Prev Med. (2017) 100:243–47. doi: 10.1016/j.ypmed.2017.05.004
55. White A, Thompson TD, White MC, Sabatino SA, de Moor J, Doria-Rose PV, et al. Cancer screening test use - United States, 2015. Morb Mortality Wkly Rep. (2017) 66:201–06. doi: 10.15585/mmwr.mm6608a1
56. Adegoke O, Kulasingam S, Virnig B. Cervical cancer trends in the United States: a 35-year population-based analysis. J Women's Health. (2012) 21:1031–37. doi: 10.1089/jwh.2011.3385
Keywords: HPV vaccine, cervical cancer, incidence, projection, vaccine impact
Citation: Pei J, Shu T, Wu C, Li M, Xu M, Jiang M and Zhu C (2023) Impact of human papillomavirus vaccine on cervical cancer epidemic: Evidence from the surveillance, epidemiology, and end results program. Front. Public Health 10:998174. doi: 10.3389/fpubh.2022.998174
Received: 19 July 2022; Accepted: 11 November 2022;
Published: 06 January 2023.
Edited by:
Monica Catarina Botelho, Instituto Nacional de Saúde Doutor Ricardo Jorge (INSA), PortugalReviewed by:
Marco Mariani, ASUR Marche, ItalyElectra Diane Paskett, The Ohio State University, United States
Copyright © 2023 Pei, Shu, Wu, Li, Xu, Jiang and Zhu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Min Jiang, 43761128@qq.com; Cairong Zhu, cairong.zhu@hotmail.com
†These authors have contributed equally to this work