Next Article in Journal
Mutation Enrichment and Transcriptomic Activation Signatures of 419 Molecular Pathways in Cancer
Next Article in Special Issue
Hypoxia-Induced Centrosome Amplification Underlies Aggressive Disease Course in HPV-Negative Oropharyngeal Squamous Cell Carcinomas
Previous Article in Journal
Molecular Mechanisms of the Efficacy of Cold Atmospheric Pressure Plasma (CAP) in Cancer Treatment
Previous Article in Special Issue
Survival-Associated Metabolic Genes in Human Papillomavirus-Positive Head and Neck Cancers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multistate Markov Model to Predict the Prognosis of High-Risk Human Papillomavirus-Related Cervical Lesions

1
Department of Obstetrics and Gynecology, Graduate School of Medicine, The University of Tokyo, Tokyo 113-8655, Japan
2
Gynecology Division, Tokyo Metropolitan Cancer and Infectious Diseases Center, Komagome Hospital, Tokyo 113-8677, Japan
3
Graduate School of Economics, The University of Tokyo, Tokyo 113-0033, Japan
4
Hematology Division, Tokyo Metropolitan Cancer and Infectious Diseases Center, Komagome Hospital, Tokyo 113-8677, Japan
5
Department of Public Health, Graduate School of Medicine, The University of Tokyo, Tokyo 113-0033, Japan
6
Department of Obstetrics and Gynecology, School of Medicine, Nihon University, Tokyo 173-8610, Japan
*
Author to whom correspondence should be addressed.
Cancers 2020, 12(2), 270; https://doi.org/10.3390/cancers12020270
Submission received: 19 November 2019 / Revised: 10 January 2020 / Accepted: 20 January 2020 / Published: 22 January 2020
(This article belongs to the Special Issue Human Papillomavirus and Cancers)

Abstract

:
Cervical intraepithelial neoplasia (CIN) has a natural history of bidirectional transition between different states. Therefore, conventional statistical models assuming a unidirectional disease progression may oversimplify CIN fate. We applied a continuous-time multistate Markov model to predict this CIN fate by addressing the probability of transitions between multiple states according to the genotypes of high-risk human papillomavirus (HPV). This retrospective cohort comprised 6022 observations in 737 patients (195 normal, 259 CIN1, and 283 CIN2 patients at the time of entry in the cohort). Patients were followed up or treated at the University of Tokyo Hospital between 2008 and 2015. Our model captured the prevalence trend satisfactory, particularly for up to two years. The estimated probabilities for 2-year transition to CIN3 or more were the highest in HPV 16-positive patients (13%, 30%, and 42% from normal, CIN1, and CIN2, respectively) compared with those in the other genotype-positive patients (3.1–9.6%, 7.6–16%, and 21–32% from normal, CIN1, and CIN2, respectively). Approximately 40% of HPV 52- or 58-related CINs remained at CIN1 and CIN2. The Markov model highlights the differences in transition and progression patterns between high-risk HPV-related CINs. HPV genotype-based management may be desirable for patients with cervical lesions.

1. Introduction

Cervical cancer is the second most commonly diagnosed type of cancer, and the third leading cause of cancer-related death among females in developed countries [1]. In 2012, an estimated 527,600 new cases of cervical cancer were diagnosed and 265,700 related deaths were reported worldwide [1]. In Japan, 33,114 women were newly diagnosed with cervical cancer in 2013 [2], and the mortality rate was 4.1 per 100,000 person-years [3]. Infection with human papillomaviruses (HPVs) is the main cause of cervical cancer development [4,5,6]. The types of HPV are categorized on the basis of their carcinogenesis. The International Agency for Research on Cancer divided the HPV genotypes into the following groups according to their carcinogenesis: the highly carcinogenic Group 1 (HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, and 59); the probably carcinogenic Group 2A (HPV 68); and the possibly carcinogenic Group 2B (HPVs 26, 30, 34, 53, 66, 67, 69, 70, 73, 82, 85, and 97) [7]. The high-risk HPV (hrHPV)-derived oncogenes E6 and E7 are necessary for malignant conversion. E6 and E7 inactivate the p53 tumor suppressor, and suppress the expression of retinoblastoma proteins, resulting in resistance to apoptosis and the promotion of cell proliferation [6]. The stabilized expression of E6 and E7 is the critical step in the progression of cervical cancer [8,9].
Worldwide, patients with a high-grade squamous intraepithelial lesion (HSIL) undergo surgical interventions (e.g., conization and loop electrosurgical excision or laser vaporization), regardless of the HPV genotype [10,11]. In Japan, surgical intervention to patients with CIN2 was not common, because approximately half of these patients regress within two years [12]. A prospective cohort study involving patients with cervical intraepithelial neoplasia (CIN) has demonstrated that, without treatment, approximately 20% of CIN2 patients with HPVs 16, 18, 31, 33, 35, 45, 52 or 58 progress to CIN3 within 5 years [12]. On the basis of the available evidence, excision strategies can be considered in patients with CIN2 as a risk-reducing treatment when the lesion is caused by these hrHPVs [13]. However, surgical excision in patients with HSIL occasionally leads to poor obstetric outcome, including preterm birth, low-birth-weight infants, and cesarean delivery, due to cervical incompetence after surgery [14,15]. Therefore, surgical treatment is unfavorable for patients who desire to be pregnant when their risk for progression to CIN3 or cancer is low.
Among the hrHPVs, HPVs 16 and 18 are the most frequently observed genotypes in patients with cervical cancer, and approximately 70% of cervical cancers are HPV 16 or 18 positive [4,16]. Several studies, including prospective cohort studies, demonstrated that the risk of HPV 16-positive patients to develop CIN2 or CIN3 lesions is higher than that reported in patients positive for other hrHPVs [12,17,18,19]. Matsumoto et al. revealed that seven hrHPV types, including HPV 16, show a high rate of progression of CIN1–2 to CIN3 compared with the other hrHPVs [12]. In addition, several reports demonstrated that patterns of persistent infection or persistent cervical lesions may differ according to the HPV types [20,21,22]. Collectively, we consider that the natural history of HPV-infected CIN lesions differs depending on the HPV genotypes; i.e., some HPV genotype-related lesions are likely to progress to cervical cancer, some are likely to regress, and others are likely to result in persistent disease.
The Cox proportional hazards model is most frequently used to predict the risk of transition from one state to another. However, the natural history of CIN shows a bidirectional transition between different states, and patients intricately move between a series of states. For example, CIN2 can progress to CIN3 or regress to CIN1 during the follow-up period. A model capable of estimating the risks of transitions from and to multiple states is more suitable for assessing the fate of CIN patients compared with the traditional Cox proportional hazards model, accounting only for the transition between two states. A solution is the use of multistate Markov models that enable investigators to estimate the transitions between multiple states [23,24,25]. Use of this approach has recently been adopted in various clinical settings [26,27,28,29,30,31,32,33,34].
In the present study, we retrospectively analyzed a cohort of patients with HPV genotype-confirmed CIN to investigate the relationship between HPV genotypes and clinicopathological features in CIN/cervical cancer. We applied a continuous-time multistate Markov model to independently estimate the prognosis of cervical lesions for designated HPV categories.

2. Results

2.1. Patients

Among the 1417 patients, 737 patients were enrolled with the following diagnoses at the time of entry: normal, 195 patients; CIN1, 259 patients; and CIN2, 283 patients. The median age was 37.6 years (interquartile range [IQR]: 31.8–44.5 years), the median number of visits was seven (IQR: 4–12), and the median duration of the follow-up was 3.03 years (IQR: 1.21–4.91 years). Table 1 displays the summary statistics for the combination of the defined HPV categories and diagnoses at the time of entry. In these patients, HPVs 16, 52, and 58 were the three most frequently observed HPV types. Because HPV 18 is the second most frequently observed HPV genotype in cervical cancer [35], in this study, we focused on HPVs 16, 52, 58, and 18.
Table 2 displays a summary of the transitions from each diagnosis of cervical epithelial lesions for the six HPV categories. Of the total 6022 transitions, 702, 270, 1008, and 852 transitions were observed for HPVs 16, 18, 52, and 58, respectively. The remaining 1501 and 2300 transitions were observed for other hrHPVs and no hrHPVs, respectively. For normal patients, 75–90% remained at the same state (e.g., 84.4% of HPV 16-positive normal patients at a visit were also diagnosed as normal at the subsequent visit). For CIN1 patients, 29–55% regressed to normal, 34–44% remained at the same state, and the rest progressed to CIN2 or CIN3. For CIN2 patients, >50% remained at the same state and the probability of progression to CIN3 was HPV type dependent: 15.6%, 10.3%, 11.1%, 7.7%, and 8.5% for HPVs 16, 18, 52, 58, and other hrHPVs, respectively. HPV 16-positive patients tended to shift to more severe states compared with those who were positive for the other HPV genotypes (e.g., 37% of those with CIN1 shifted to more severe states at the subsequent visit).

2.2. Prognosis of HPV-Infected Cervical Lesions Estimated Using the Markov Model

Table 3 shows the predicted probabilities of transitions from states to states within two years for the six HPV categories, as estimated by the Markov model. Regarding HPV 16-positive patients, 13%, 30%, and 42% of normal, CIN1, and CIN2 patients, respectively, progressed to CIN3/cancer. In contrast, 43% and 34% of CIN1 and CIN2 patients, respectively, regressed to the normal state. The fates of HPV 18-infected cervical lesions were similar to those of HPV 16-infected cervical lesions. However, exceptions to these were the probabilities for progression to CIN3/cancer: from the normal state, 7.6%; from CIN1, 15%; and from CIN2, 32%, respectively. The fates of HPV 52/58-positive cervical lesions differed considerably versus those of HPV 16/18-positive cervical lesions: one-third to two-fifths of HPV 52/58-positive patients were eventually diagnosed with CIN1 or CIN2 regardless of the initial state, and their probability of progression from CIN2 to CIN3/cancer was approximately 25%. Among the hrHPV-positive patients, those with other hrHPVs were most likely to regress to the normal state and least likely to progress to CIN3/cancer: probability of regression to normal state, 66% and 52% from CIN1 and CIN2, respectively, and probability of progression to CIN3/cancer, 3.4%, 8.4%, and 22% from normal, CIN1, and CIN2, respectively. The probability of regression to the normal state of no hrHPVs patients was higher than that of hrHPV-positive patients: 81% and 71% for CIN1 and CIN2, respectively. The Cox model yielded the progression probability from CIN1 to CIN2 or more severe lesions to be 50%, 9%, 53%, 50%, 25%, and 13% for HPVs 16, 18, 52, 58 other hrHPVs, and no hrHPVs, respectively (Supplementary Table S1).
Figure 1 demonstrates the observed and simulated prevalence transition of each state for the HPV categories. Overall, our model captured the trend of the prevalence, especially up to two years. However, as time progressed, the overestimation (underestimation) of the prevalence of the normal (CIN3/cancer) state expanded.
Supplementary Tables S2 and S3 display the summary statistics and the summary of the transitions from the dataset of the alternative model. Supplementary Table S4 and Supplementary Figure S1 show the results of all sensitivity analyses using the dataset and categories of the alternative model. These data indicate that the results of the main model were fairly robust when we altered the study population and categories.

3. Discussion

This study is the first to estimate the parameters of the continuous-time multistate Markov model from one data set to shed light on the prognosis of cervical lesions based on the infected HPV types. The results indicated that the prognosis of hrHPV-infected cervical lesions differs according to the types of HPV and the grades of CIN.
The Markov model revealed that the rate of progression to CIN3 from any condition was highest in HPV 16-positive patients, followed by HPV 18-positive patients. These results are consistent with those reported in previous studies. In addition, our Markov model also revealed that HPV 52/58-positive patients tend to remain between CIN1 and CIN2. The Markov model clarifies the characteristics of HPV 52- and HPV 58-related lesions, as it demonstrates the probabilities of transitioning from state to state. Contrary to HPV 16-positive lesions, which highly progress to CIN3/cancer, HPV 52- and 58-related lesions were characterized by their stability between CIN1 and CIN2. In this study, we estimated the 2-year prognosis of patients with CIN using two types of models. The main model reflects transitions based upon objective results, and it is not always in line with clinical practice, because some patients with CIN2 undergo treatment without a final diagnosis of CIN3. On the other hand, the alternative model well reflects the trajectory of clinical practice. However, treatment is occasionally affected by the physicians’ decision and patients’ choice; therefore, the alternative model may lack objectivity. Admitting that our models both have advantages and disadvantages, we consider the results for the 2-year prognosis of CIN patients to be reasonably robust, as both models essentially yielded identical results.
The continuous-time multistate Markov model is more capable of reflecting the features of HPV-infected cervical lesions than the frequently used Cox proportional hazards model. Most of the previous studies have focused on the risk of progression to CIN3 or cervical cancer according to the types of HPV or grades of CIN, implicitly assuming a unidirectional disease progression. Considering the bidirectional nature of CIN, the Markov model has the potential to more accurately estimate the fates of HPV-infected cervical lesions than the standard Cox proportional hazards model. Indeed, the Markov model and the Cox model led to quite different predicted 2-year transition probabilities from CIN1 to CIN2 or more severe lesions. Notably, we could not detect higher disease progression risk in HPV 16-positive patients than in HPV 52/58-positive patients in the Cox model. The Markov model can estimate the transition between multiple states besides the probability of the customarily used endpoint, e.g., CIN3 and cervical cancer. The use of this model highlighted the difference between the natural history of cervical lesions infected with the focused four HPV types. We conclude that the model may improve our understanding of the natural history of cervical lesions and aid in the management of hrHPV-related lesions.
The Markov model was suitable for HPV-infected cervical lesions because the virus–host interaction (virus vs. host immune system) leads to bidirectional transitions between states. Activated host immune responses contribute to CIN regression, whereas immunosuppressive conditions promote CIN progression [36,37,38]. Likewise, this model can be applied to virus-induced persistent conditions, such as hepatitis B virus- or hepatitis C virus-infected liver conditions, and human immunodeficiency virus-related conditions. Indeed, previous studies applied the model to hepatitis C virus-related liver diseases with three states (i.e., none or mild fibrosis, moderate fibrosis and cirrhosis) to quantify differences in the prognosis among cohorts [34].
In this study, we could not fully estimate the transitions of HPV 18-positive patients with the present Markov model due to the limited sample size. HPV 18 is the second most frequently observed HPV genotype in patients with cervical cancer, whereas it is not frequently detected in those with precancer lesions [35]. In addition, the detection rate of HPV 18 is higher in adenocarcinoma than that in squamous cell carcinoma [39]. Further accumulation of cases may enable us to model the complicated transitions of HPV 18-positive patients.
This study has several limitations. First, we performed the HPV genotyping test once for each patient. Therefore, the change to negative for the HPV test (latent infection or clearance) was not taken into account. Previous reports demonstrate that patients with HPV clearance tend to regress compared with those with persistent HPV infection [40,41]. Patients with the transient HPV infection and persistent infection might take different courses.
Second, in this study, we did not evaluate the effect of concurrent multiple HPV infections. In previous reports, coinfections with multiple α9 species were associated with an increased risk of CIN2 or more severe lesions compared with single infections [42,43]. Because there are many different possible combinations of concurrent infection, the numbers of patients for each possible combination was expected to be too small for an analysis of the effect of concurrent infection. Further analysis with a larger sample size is warranted to unveil the fates of cervical lesions with multiple HPV infections.
Third, the evaluation of model fitness uncovered that the simulated prevalence of the normal state (CIN3/cancer) was overstated (understated) compared with the observed prevalence. These discrepancies imply a violation of the Markov property, i.e., the distribution of the next states depends only upon the present state. Misclassification and omitted covariates, such as age, can be a potential source of the violation. The application of a hidden Markov model that can incorporate the possibility of misclassification through the information from the previous or more distant visits may be an intriguing direction for further research.
Fourth, we have not assessed the validity of our model with the use of a test set separate from the training set, which is the dataset used for parameter estimation. Instead, we prioritized improving the parameter estimation accuracy by using all data for parameter estimation. This is currently a common practice in the literature [26,27,28,29,30,31,32,33,34]. If a sufficient sample size can be acquired, the external validity can be assessed to some degree using a test set, which will be an issue for future research.
Lastly, this study was based on data obtained from a single institute with some censoring. Especially, there is a selection bias; for example, patients enrolled in this study previously had abnormal cytology. In addition, some patients with relatively severe CIN2 were censored because of the treatment according to the physician’s decision. However, by conducting the sensitivity analysis, we confirmed that the impact on the results was limited. Nevertheless, as we have not assessed the external validity of our model in the general population, care must be taken when extrapolating the present results to the entire population, unless similar results are yielded from other institutes.

4. Materials and Methods

4.1. Study Design and Patients

This study was performed in accordance with the Declaration of Helsinki. Approval was obtained from the internal institutional review board (Approval number: 1390-1, G10082-7, and G0637-6) for this study.
This study was a retrospective cohort analysis of data extracted from the electronic health records of a teaching hospital in Japan. Since 2008, human papillomavirus (HPV) genotyping has been performed in the University of Tokyo Hospital (Tokyo, Japan) for outpatients. The patients who had abnormal cytology in the population-based screening, or whose abnormal cytology was found in the outpatient visits at the Obstetrics and Gynecology Department of the University of Tokyo Hospital or other hospitals, were enrolled. At the first visit, we confirmed the diagnosis by performing punch biopsy under colposcopic examination (based on the histological diagnosis). As follow-up screening, we routinely performed combination examinations of cytology and colposcopy because these examinations are non-invasive. We performed histological diagnosis only when the disease progressed. We reviewed the electronic health records of 1417 patients for whom genotyping was performed between October 1, 2008, and March 31, 2015, to construct a retrospective cohort. Patients were included in the study if they were (i) diagnosed with normal cervical lesion, CIN1, or CIN2 and (ii) observed for at least two visits during the study period. Patients were excluded if they had HPV 6- or HPV 11-single-positive lesions with only the diagnosis of condyloma during their follow-up period, or had only glandular lesions. Patients were followed up until they were diagnosed with cervical cancer, underwent treatment, transferred to another hospital, or until the date March 31, 2018, whichever came first. In our practical management, most patients with CIN2 were followed up without any treatment. Patients diagnosed with CIN3 underwent surgical intervention, such as conization and loop electrosurgical excision or laser vaporization.

4.2. Variables

The date of birth was extracted for each patient. Patient age was defined as their age at the time of entry. Follow-up interval was defined as the length of time between two consecutive visits.
Patient cytological and histological results and the date were recorded for each visit. Cytological and histological results were combined to classify the results into any of the following four diagnoses: normal, CIN1, CIN2, and CIN3/cancer. The investigators (experts in gynecologic oncology) convened and determined the criteria for pathological diagnosis as follows: (1) CIN1–2 was classified into CIN2; (2) CIN2–3 and carcinoma in situ were classified into CIN3; (3) uncertain diagnoses (e.g., atypical squamous cells of uncertain significance, atypical squamous cells that cannot exclude HSIL, and dysplasia without grading) were excluded from the study due to concerns regarding diagnostic reliability; and (4) in the presence of histological and cytological examinations, the most severe classification was adopted as the final diagnosis.
The results of the HPV genotyping in cervical samples collected using swabs were recorded. Genotyping was performed once for each patient; thus, the HPV type assigned to a patient did not change over time. It was permitted to assign multiple genotypes to a single patient. On the basis of the classification of the International Agency for Research on Cancer, we defined HPVs classified in Group 1 or Group 2A (HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68) as “high-risk HPVs (hrHPVs)”. Of these, HPVs 16, 18, 52, and 58 were separately categorized. hrHPVs other than HPVs 16, 18, 52, and 58 were classified as “other hrHPVs.” Patients who were not infected with any hrHPVs were referred to as “no hrHPVs” patients. Patients without HPV infection were also categorized in this group.

4.3. DNA Extraction and HPV Genotyping

Cervical samples were tested for HPV DNA using PGMY line-blot hybridization, as previously described [44]. DNA was extracted from cervical samples using the DNeasy Blood Mini Kit (Qiagen, Crawley, UK). HPV genotyping was performed using the PGMY-CHUV assay method. Briefly, standard polymerase chain reaction (PCR) was performed using the PGMY09⁄11 L1 consensus primer set and the human leukocyte antigen-DQ primer set, as previously described [44]. Reverse blotting hybridization was performed. Heat-denatured PCR amplicons were hybridized to probes specific for 32 HPV genotypes and the human leukocyte antigen-DQ references.

4.4. Continuous-Time Multistate Markov Model

We used the continuous-time multistate Markov model to estimate the prognosis of each patient with HPV-infected cervical lesions [23,24,25]. We defined the following four states: normal (state 1), CIN1 (state 2), CIN2 (state 3), and CIN3/cancer (state 4) (Figure 2). The arrows in Figure 2 specify possible transitions between the states defined in our model; all transitions between adjacent states, except the backward transition from CIN3/cancer to CIN2, were allowed. CIN3/cancer was the absorbing state.
Each transition parameter λ in Figure 2 indicates the transition intensity; i.e., λij is interpreted as an “instantaneous risk” of transition from state i to j [24]. Identification of transition parameters in the continuous-time multistate Markov model requires observation of the corresponding transitions.

4.5. Dataset Construction

Initially, an unbalanced panel data with the unit of observation being one visit was constructed on the basis of the chart review. We truncated observations after the diagnosis of CIN3 or cancer to make the data used in the estimation compatible with the model in Figure 2. Additionally, we excluded samples that left only one observation after the truncation to estimate the model, as these samples do not contribute to the estimation. Figure 3 shows the transitions of typical samples from the baseline dataset, which illustrates that the patterns of disease progression or regression varied among patients.

4.6. Statistical Analysis

Summary statistics are reported for the combination of defined HPV categories (HPV 16, HPV 18, HPV 52, HPV 58, other hrHPVs, and no hrHPVs) and diagnoses at the time of entry (normal, CIN1, and CIN2). In cases of observed coinfection with different HPV genotypes, it was possible to include the same patient in the summary statistics of multiple HPV categories.
A maximum likelihood estimation was performed to estimate the parameters using the msm package in R [45]. Intuitively, the likelihood function consists of the sum of the probabilities of all possible paths given the observations, and each observed transition contributes to it. Parameters for the defined HPV categories were independently estimated. In cases of coinfection, the patients could contribute to the parameter estimation of several HPV categories. According to the estimated parameters, we simulated the probabilities of transitioning from state to state after two years. This 2-year transition probability was selected for two reasons. First, since the typical follow-up duration of this study was <5 years (upper quartile value of the follow-up duration was 4.9 years), a prediction beyond this period was unwarranted. Second, a 2-year prediction is considered a good benchmark for the prognosis of hrHPV-related cervical lesions. Notably, ≥50% of these lesions regress to the normal state, and approximately 10% of these lesions progress to CIN3 within two years [12]. The Cox proportional hazards model was also used to predict the 2-year transition probability from CIN1 to CIN2 or more severe lesions for each HPV category. To this end, we restricted the patients to those who were diagnosed with CIN1 at the time of entry and followed them until they were diagnosed as CIN2 or more severe lesions or at the end of their observation period.
Additionally, we simulated the change in the prevalence of each state over the five years for the HPV categories. To make it easier to compare the observed and the simulated prevalence, the initial prevalence of each state was set to be the same as the data used for parameter estimation. The fit of the model was assessed on the basis of the comparison of the observed and the simulated prevalence.

4.7. Sensitivity Analysis

Some patients with CIN2 underwent treatment without a final diagnosis of CIN3. Some patients with CIN3 were observed without any treatment if the lesion immediately regressed to CIN2 or a lower state. Therefore, we specified an alternative model that encompasses these situations with the following four states: normal (state 1), CIN1 (state 2), CIN2/CIN3 (state 3), and treatment (state 4) (Supplementary Figure S2). The arrows in Supplementary Figure S2 specify possible transitions between these states; all transitions between adjacent states, except the backward transition from treatment to CIN2/CIN3, were allowed. Treatment was the absorbing state. We truncated observations after the diagnosis of treatment intervention in this alternative model, and excluded samples that left only one observation after the truncation to estimate the model. Typical sampling situations from this alternative dataset are illustrated in Supplementary Figure S3.

5. Conclusions

In this study, we applied a continuous-time multistate Markov model to predict the prognosis of patients with HPV-infected cervical lesions. We demonstrated that the Markov model is a promising analytical method considering the bidirectional nature of these lesions. The study revealed that the natural history of hrHPV-related cervical lesions differed among the focused four HPV types: HPV 16-positive lesions were likely to be upgraded to CIN states in a step-by-step manner; HPV 52/58-positive lesions were likely to be maintained between CIN1 and CIN2; lesions positive for the other hrHPVs were most likely to regress to the normal state and least likely to progress to CIN3/cancer. On the basis of the present findings, HPV genotype-based management may be desirable for patients with cervical lesions.

Supplementary Materials

The following are available online at https://www.mdpi.com/2072-6694/12/2/270/s1, Supplementary Materials: Figure S1: Observed and simulated prevalence transition of each state for the human papillomavirus (HPV) categories in the alternative model. Figure S2: Markov model for the disease progression and regression of cervical epithelial lesions of the alternative model. Figure S3: Possible transition paths of the alternative model for the selected patients. Table S1: Predicted 2-year transition probabilities from CIN1 to CIN2 or more severe lesions and their 95% confidence intervals for the six HPV categories derived from the Cox regression analysis. Table S2: Sample size and summary statistics for the combination of patients in the six HPV categories and diagnoses at the time of entry in the dataset for the alternative model. Table S3: Summary of the transitions from each diagnosis of cervical epithelial lesions for the six HPV categories in the dataset for the alternative model. Table S4: Predicted 2-year transition probabilities from states to states and their 95% confidence intervals for the six HPV categories in the alternative model.

Author Contributions

Conceptualization, A.T., K.H., and K.K.; methodology, A.T. and K.H.; data collection and data curation, A.T., K.H., T.T. (Tomoki Tanaka), S.B., A.K., S.E., M.M., T.T. (Tetsushi Tsuruga), and K.A.; formal analysis, K.H.; writing—original draft preparation, A.T. and K.H.; writing—review and editing, J.T., T.N., K.O., and K.K.; supervision, T.N., K.O., T.Y., J.T., K.K, Y.O., and T.F.; funding acquisition, A.T. and K.K. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Grants-in-Aid for Practical Research for Innovative Cancer Control (KK, Grant Number: 15656298) and J-PRIDE (AT, Grant number: 19fm0208013h0003) from the Japan Agency for Medical Research and Development (AMED).

Acknowledgments

We thank Terufumi Yokoyama for the excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Torre, L.A.; Bray, F.; Siegel, R.L.; Ferlay, J.; Lortet-Tieulent, J.; Jemal, A. Global cancer statistics, 2012. Ca Cancer J. Clin. 2015, 65, 87–108. [Google Scholar] [CrossRef] [Green Version]
  2. National Cancer Center. National Estimates of Cancer Incidence based on Cancer Registries in Japan (1975–2013). Available online: https://ganjoho.jp/en/professional/statistics/table_download.html (accessed on 1 January 2019).
  3. National Cancer Center. Cancer mortality from Vital Statistics in Japan (1958–2016). Available online: https://ganjoho.jp/en/professional/statistics/table_download.html (accessed on 1 January 2019).
  4. Lowy, D.R.; Schiller, J.T. Reducing HPV-Associated Cancer Globally. Cancer Prev. Res. 2012, 5, 18–23. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Viens, L.J.; Henley, S.J.; Watson, M.; Markowitz, L.E.; Thomas, C.C.; Thompson, T.D.; Razzaghi, H.; Saraiya, M. Human Papillomavirus–Associated Cancers—United States, 2008–2012. Mmwr. Morb. Mortal. Wkly. Rep. 2016, 65, 661–666. [Google Scholar] [CrossRef] [PubMed]
  6. Yim, E.; Park, J. The Role of HPV E6 and E7 Oncoproteins in HPV-associated Cervical Carcinogenesis. Cancer Res. Treat. 2005, 37, 319. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Bzhalava, D.; Guan, P.; Franceschi, S.; Dillner, J.; Clifford, G. A systematic review of the prevalence of mucosal and cutaneous human papillomavirus types. Virology 2013, 445, 224–231. [Google Scholar] [CrossRef] [Green Version]
  8. 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]
  9. McBride, A.A.; Warburton, A. The role of integration in oncogenic progression of HPV-associated cancers. Plos Pathog. 2017, 13, e1006211. [Google Scholar] [CrossRef] [Green Version]
  10. Wright, T.C.; Massad, L.S.; Dunton, C.J.; Spitzer, M.; Wilkinson, E.J.; Solomon, D. 2006 consensus guidelines for the management of women with cervical intraepithelial neoplasia or adenocarcinoma in situ. Am. J. Obs. Gynecol. 2007, 197, 340–345. [Google Scholar] [CrossRef]
  11. Wright, T.C.; Massad, L.S.; Dunton, C.J.; Spitzer, M.; Wilkinson, E.J.; Solomon, D. 2006 consensus guidelines for the management of women with abnormal cervical cancer screening tests. Am. J. Obs. Gynecol. 2007, 197, 346–355. [Google Scholar] [CrossRef]
  12. Matsumoto, K.; Oki, A.; Furuta, R.; Maeda, H.; Yasugi, T.; Takatsuka, N.; Mitsuhashi, A.; Fujii, T.; Hirai, Y.; Iwasaka, T.; et al. Predicting the progression of cervical precursor lesions by human papillomavirus genotyping: A prospective cohort study. Int. J. Cancer 2011, 128, 2898–2910. [Google Scholar] [CrossRef]
  13. Japan Society of Obstetrics and Gynecology; Japan Association of Obstetricians and Gynecologists (Eds.) Guideline for Gynecological Practice 2017; Japan Society of Obstetrics and Gynecology: Tokyo, Japan, 2017; pp. 53–57. (in Japanese) [Google Scholar]
  14. Kyrgiou, M.; Koliopoulos, G.; Martin-Hirsch, P.; Arbyn, M.; Prendiville, W.; Paraskevaidis, E. Obstetric outcomes after conservative treatment for intraepithelial or early invasive cervical lesions: Systematic review and meta-analysis. Lancet 2006, 367, 489–498. [Google Scholar] [CrossRef]
  15. Kyrgiou, M.; Athanasiou, A.; Paraskevaidi, M.; Mitra, A.; Kalliala, I.; Martin-Hirsch, P.; Arbyn, M.; Bennett, P.; Paraskevaidis, E. Adverse obstetric outcomes after local treatment for cervical preinvasive and early invasive disease according to cone depth: Systematic review and meta-analysis. BMJ 2016, 71, i3633. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Miura, S.; Matsumoto, K.; Oki, A.; Satoh, T.; Tsunoda, H.; Yasugi, T.; Taketani, Y.; Yoshikawa, H. Do we need a different strategy for HPV screening and vaccination in East Asia? Int. J. Cancer 2006, 119, 2713–3715. [Google Scholar] [CrossRef] [PubMed]
  17. Adebamowo, S.N.; Olawande, O.; Famooto, A.; Dareng, E.O.; Offiong, R.; Adebamowo, C.A. Persistent Low-Risk and High-Risk Human Papillomavirus Infections of the Uterine Cervix in HIV-Negative and HIV-Positive Women. Front. Public Heal. 2017, 5, 1–11. [Google Scholar] [CrossRef]
  18. Wright, T.C.; Stoler, M.H.; Behrens, C.M.; Sharma, A.; Zhang, G.; Wright, T.L. Primary cervical cancer screening with human papillomavirus: End of study results from the ATHENA study using HPV as the first-line screening test. Gynecol. Oncol. 2015, 136, 189–197. [Google Scholar] [CrossRef] [Green Version]
  19. Naucler, P.; Ryd, W.; Törnberg, S.; Strand, A.; Wadell, G.; Hansson, B.G.; Rylander, E.; Dillner, J. HPV type-specific risks of high-grade CIN during 4 years of follow-up: A population-based prospective study. Br. J. Cancer 2007, 97, 129–132. [Google Scholar] [CrossRef] [Green Version]
  20. Rositch, A.F.; Koshiol, J.; Hudgens, M.G.; Razzaghi, H.; Backes, D.M.; Pimenta, J.M.; Franco, E.L.; Poole, C.; Smith, J.S. Patterns of persistent genital human papillomavirus infection among women worldwide: A literature review and meta-analysis. Int. J. Cancer 2013, 133, 1271–1285. [Google Scholar] [CrossRef] [Green Version]
  21. Rodriguez, A.C.; Schiffman, M.; Herrero, R.; Wacholder, S.; Hildesheim, A.; Castle, P.E.; Solomon, D.; Burk, R. Rapid Clearance of Human Papillomavirus and Implications for Clinical Focus on Persistent Infections. Jnci. J. Natl. Cancer Inst. 2008, 100, 513–517. [Google Scholar] [CrossRef] [Green Version]
  22. Schiffman, M.; Herrero, R.; DeSalle, R.; Hildesheim, A.; Wacholder, S.; Cecilia Rodriguez, A.; Bratti, M.C.; Sherman, M.E.; Morales, J.; Guillen, D.; et al. The carcinogenicity of human papillomavirus types reflects viral evolution. Virology 2005, 337, 76–84. [Google Scholar] [CrossRef] [Green Version]
  23. Cox, D.R.; Miller, H.D. The Theory of Stochastic Processes, 1st ed.; Routledge: London, NY, USA; Chapman & Hall/CRC: London, UK, 1965. [Google Scholar]
  24. Kalbfleisch, J.D.; Lawless, J.F. The Analysis of Panel Data under a Markov Assumption. J. Am. Stat. Assoc. 1985, 80, 863–871. [Google Scholar] [CrossRef]
  25. Kay, R. A Markov Model for Analysing Cancer Markers and Disease States in Survival Studies. Biometrics 1986, 42, 855–865. [Google Scholar] [CrossRef] [PubMed]
  26. Buchman, A.S.; Leurgans, S.E.; Yu, L.; Wilson, R.S.; Lim, A.S.; James, B.D.; Shulman, J.M.; Bennett, D.A. Incident parkinsonism in older adults without Parkinson disease. Neurology 2016, 87, 1036–1044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Buter, T.C.; van den Hout, A.; Matthews, F.E.; Larsen, J.P.; Brayne, C.; Aarsland, D. Dementia and survival in Parkinson disease: A 12-year population study. Neurology 2008, 70, 1017–1022. [Google Scholar] [CrossRef] [PubMed]
  28. Dancourt, V.; Quantin, C.; Abrahamowicz, M.; Binquet, C.; Alioum, A.; Faivre, J. Modeling recurrence in colorectal cancer. J. Clin. Epidemiol. 2004, 57, 243–251. [Google Scholar] [CrossRef]
  29. Jack, C.R.; Therneau, T.M.; Wiste, H.J.; Weigand, S.D.; Knopman, D.S.; Lowe, V.J.; Mielke, M.M.; Vemuri, P.; Roberts, R.O.; Machulda, M.M.; et al. Transition rates between amyloid and neurodegeneration biomarker states and to dementia: A population-based, longitudinal cohort study. Lancet Neurol. 2016, 15, 56–64. [Google Scholar] [CrossRef] [Green Version]
  30. Pan, S.-L.; Lien, I.-N.; Yen, M.-F.; Lee, T.-K.; Chen, T.H.-H. Dynamic Aspect of Functional Recovery after Stroke Using a Multistate Model. Arch. Phys. Med. Rehabil. 2008, 89, 1054–1060. [Google Scholar] [CrossRef]
  31. Price, M.J.; Ades, A.E.; De Angelis, D.; Welton, N.J.; Macleod, J.; Soldan, K.; Simms, I.; Turner, K.; Horner, P.J. Risk of Pelvic Inflammatory Disease Following Chlamydia trachomatis Infection: Analysis of Prospective Studies With a Multistate Model. Am. J. Epidemiol. 2013, 178, 484–492. [Google Scholar] [CrossRef] [Green Version]
  32. Sharples, L.D.; Jackson, C.H.; Parameshwar, J.; Wallwork, J.; Large, S.R. Diagnostic accuracy of coronary angiography and risk factors for post-heart-transplant cardiac allograft vasculopathy. Transplantation 2003, 76, 679–682. [Google Scholar] [CrossRef]
  33. Skogvoll, E.; Eftestøl, T.; Gundersen, K.; Kvaløy, J.T.; Kramer-Johansen, J.; Olasveengen, T.M.; Steen, P.A. Dynamics and state transitions during resuscitation in out-of-hospital cardiac arrest. Resuscitation 2008, 78, 30–37. [Google Scholar] [CrossRef]
  34. Sweeting, M.J.; De Angelis, D.; Neal, K.R.; Ramsay, M.E.; Irving, W.L.; Wright, M.; Brant, L.; Harris, H.E. Estimated progression rates in three United Kingdom hepatitis C cohorts differed according to method of recruitment. J. Clin. Epidemiol. 2006, 59, 144–152. [Google Scholar] [CrossRef]
  35. Bulk, S.; Berkhof, J.; Rozendaal, L.; Fransen Daalmeijer, N.C.; Gök, M.; de Schipper, F.A.; van Kemenade, F.J.; Snijders, P.J.; Meijer, C.J. The contribution of HPV18 to cervical cancer is underestimated using high-grade CIN as a measure of screening efficiency. Br. J. Cancer 2007, 96, 1234–1236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. De Gruijl, T.D.; Bontkes, H.J.; Walboomers, J.M.M.; Stukart, M.J.; Doekhie, F.S.; Remmink, A.J.; Helmerhorst, T.J.M.; Verheijen, R.H.M.; Duggan-Keen, M.F.; Stern, P.L.; et al. Differential T helper cell responses to human papillomavirus type 16 E7 related to viral clearance or persistence in patients with cervical neoplasia: A longitudinal study. Cancer Res. 1998, 58, 1700–1706. [Google Scholar] [PubMed]
  37. Ahdieh, L.; Munoz, A.; Vlahov, D.; Trimble, C.L.; Timpson, L.A.; Shah, K. Cervical Neoplasia and Repeated Positivity of Human Papillomavirus Infection In Human Immunodeficiency Virus-seropositive and -seronegative Women. Am. J. Epidemiol. 2000, 151, 1148–1157. [Google Scholar] [CrossRef] [Green Version]
  38. Brennan, D.C.; Aguado, J.M.; Potena, L.; Jardine, A.G.; Legendre, C.; Säemann, M.D.; Mueller, N.J.; Merville, P.; Emery, V.; Nashan, B. Effect of maintenance immunosuppressive drugs on virus pathobiology: Evidence and potential mechanisms. Rev. Med. Virol. 2013, 23, 97–125. [Google Scholar] [CrossRef]
  39. De Sanjose, S.; Quint, W.G.V.; Alemany, L.; Geraets, D.T.; Klaustermeier, J.E.; Lloveras, B.; Tous, S.; Felix, A.; Bravo, L.E.; Shin, H.-R.; et al. Human papillomavirus genotype attribution in invasive cervical cancer: A retrospective cross-sectional worldwide study. Lancet Oncol. 2010, 11, 1048–1056. [Google Scholar] [CrossRef]
  40. Brown, D.R.; Shew, M.L.; Qadadri, B.; Neptune, N.; Vargas, M.; Tu, W.; Juliar, B.E.; Breen, T.E.; Fortenberry, J.D. A longitudinal study of genital human papillomavirus infection in a cohort of closely followed adolescent women. J. Infect. Dis. 2005, 191, 182–192. [Google Scholar] [CrossRef] [Green Version]
  41. Cho, H.W.; So, K.A.; Lee, J.K.; Hong, J.H. Type-specific persistence or regression of human papillomavirus genotypes in women with cervical intraepithelial neoplasia 1: A prospective cohort study. Obstet. Gynecol. Sci. 2015, 58, 40–45. [Google Scholar] [CrossRef] [Green Version]
  42. Wentzensen, N.; Schiffman, M.; Dunn, T.; Zuna, R.E.; Gold, M.A.; Allen, R.A.; Zhang, R.; Sherman, M.E.; Wacholder, S.; Walker, J.; et al. Multiple human papillomavirus genotype infections in cervical cancer progression in the study to understand cervical cancer early endpoints and determinants. Int. J. Cancer 2009, 125, 2151–2158. [Google Scholar] [CrossRef] [Green Version]
  43. Chaturvedi, A.K.; Katki, H.A.; Hildesheim, A.; Rodríguez, A.C.; Quint, W.; Schiffman, M.; Van Doorn, L.-J.; Porras, C.; Wacholder, S.; Gonzalez, P.; et al. Human Papillomavirus Infection with Multiple Types: Pattern of Coinfection and Risk of Cervical Disease. J. Infect. Dis. 2011, 203, 910–920. [Google Scholar] [CrossRef] [Green Version]
  44. Kojima, S.; Kawana, K.; Tomio, K.; Yamashita, A.; Taguchi, A.; Miura, S.; Adachi, K.; Nagamatsu, T.; Nagasaka, K.; Matsumoto, Y.; et al. The Prevalence of Cervical Regulatory T Cells in HPV-Related Cervical Intraepithelial Neoplasia (CIN) Correlates Inversely with Spontaneous Regression of CIN. Am. J. Reprod Immunol. 2013, 69, 134–141. [Google Scholar] [CrossRef]
  45. Jackson, C.H. Multi-State Models for Panel Data: The msm Package for R. J. Stat. Softw. 2011, 38, 1–29. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Observed and simulated prevalence transition of each state for the human papillomavirus (HPV) categories. The figure demonstrates the observed (solid line) and simulated (dotted line) prevalence transition of each state for the six HPV categories. Cytological and histological results were combined to classify the results into the states of the model. HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 were classified as high-risk HPVs (hrHPVs). Of these, HPV 16, 18, 52, and 58 were separately categorized. hrHPVs other than HPVs 16, 18, 52, and 58 were classified as “other hrHPVs.” Patients who were not infected with any hrHPVs were referred to as “no hrHPVs” patients. In cases of observed coinfection with different HPV genotypes, it was possible to include the same patient in multiple HPV categories.
Figure 1. Observed and simulated prevalence transition of each state for the human papillomavirus (HPV) categories. The figure demonstrates the observed (solid line) and simulated (dotted line) prevalence transition of each state for the six HPV categories. Cytological and histological results were combined to classify the results into the states of the model. HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 were classified as high-risk HPVs (hrHPVs). Of these, HPV 16, 18, 52, and 58 were separately categorized. hrHPVs other than HPVs 16, 18, 52, and 58 were classified as “other hrHPVs.” Patients who were not infected with any hrHPVs were referred to as “no hrHPVs” patients. In cases of observed coinfection with different HPV genotypes, it was possible to include the same patient in multiple HPV categories.
Cancers 12 00270 g001
Figure 2. Markov model for the disease progression and regression of cervical epithelial lesions. The figure displays the schema of the Markov model for the model specified in this study. We defined four states: normal (state 1), cervical intraepithelial neoplasia 1 (CIN1, state 2), CIN2 (state 3), and CIN3/cancer (state 4). The arrows in the figure specify possible transitions between these states; all transitions between adjacent states, except the backward transition from CIN3/cancer to CIN2, were allowed. CIN3/cancer was the absorbing state. Each transition parameter λ indicates the transition intensity; i.e., λij is interpreted as an “instantaneous risk” of transition from state i to j.
Figure 2. Markov model for the disease progression and regression of cervical epithelial lesions. The figure displays the schema of the Markov model for the model specified in this study. We defined four states: normal (state 1), cervical intraepithelial neoplasia 1 (CIN1, state 2), CIN2 (state 3), and CIN3/cancer (state 4). The arrows in the figure specify possible transitions between these states; all transitions between adjacent states, except the backward transition from CIN3/cancer to CIN2, were allowed. CIN3/cancer was the absorbing state. Each transition parameter λ indicates the transition intensity; i.e., λij is interpreted as an “instantaneous risk” of transition from state i to j.
Cancers 12 00270 g002
Figure 3. Possible transition paths for the selected patients. The figure shows possible transition paths for four selected patients (patient 01–patient 04). Cytological and histological results were combined to classify the results into the states of the model. We truncated observations after the diagnosis of cervical intraepithelial neoplasia 3 or cancer. The filled circles indicate actual observations or visits. The solid line is a possible transition path during the follow-up period. The possible transition paths were randomly selected on the basis of the observed states.
Figure 3. Possible transition paths for the selected patients. The figure shows possible transition paths for four selected patients (patient 01–patient 04). Cytological and histological results were combined to classify the results into the states of the model. We truncated observations after the diagnosis of cervical intraepithelial neoplasia 3 or cancer. The filled circles indicate actual observations or visits. The solid line is a possible transition path during the follow-up period. The possible transition paths were randomly selected on the basis of the observed states.
Cancers 12 00270 g003
Table 1. Sample size and summary statistics for the combination of patients in the six HPV categories and diagnoses at the time of entry.
Table 1. Sample size and summary statistics for the combination of patients in the six HPV categories and diagnoses at the time of entry.
Diagnosis at the Time of Entry HPV 16HPV 18HPV 52HPV 58Other hrHPVsNo hrHPVsAll
NormalN1311171330122195
Age at the time of entry (years), mean (SD)42.4 (13.8)39.5 (15.0)36.7 (10.0)41.3 (16.8)42.5 (16.9)41.2 (10.5)41.3 (12.1)
Number of visits, mean (SD)5.1 (4.8)7.2 (3.4)7.5 (5.7)6.8 (4.6)7.0 (4.8)6.5 (3.9)6.8 (4.3)
Follow-up interval (years), mean (SD)0.49 (0.37)0.46 (0.27)0.50 (0.35)0.49 (0.42)0.47 (0.27)0.53 (0.40)0.51 (0.37)
Follow-up period (years), mean (SD)2.1 (2.5)2.9 (1.6)3.7 (3.1)2.9 (2.3)2.9 (2.2)3.1 (2.1)3.1 (2.2)
CIN1N2311382479111259
Age at the time of entry (years), mean (SD)34.6 (8.2)33.0 (10.0)36.6 (8.6)36.1 (7.7)34.5 (7.2)39.0 (10.3)36.9 (9.2)
Number of visits, mean (SD)7.6 (5.3)7.6 (4.7)10.7 (5.9)10.5 (6.4)9.0 (4.9)9.3 (5.2)9.3 (5.2)
Follow-up interval (years), mean (SD)0.42 (0.43)0.51 (0.53)0.38 (0.18)0.45 (0.42)0.39 (0.22)0.42 (0.31)0.41 (0.29)
Follow-up period (years), mean (SD)3.2 (2.6)3.2 (2.2)4.0 (2.4)3.9 (2.6)3.4 (2.2)3.9 (2.3)3.7 (2.3)
CIN2N671565576754283
Age at the time of entry (years), mean (SD)37.6 (7.9)42.0 (5.7)41.2 (8.2)39.7 (8.5)37.9 (8.5)36.6 (9.2)39.1 (8.4)
Number of visits, mean (SD)6.8 (5.4)7.0 (4.4)7.2 (5.4)8.9 (5.8)8.6 (5.6)8.6 (5.8)8.0 (5.5)
Follow-up interval (years), mean (SD)0.35 (0.15)0.34 (0.14)0.38 (0.24)0.38 (0.24)0.39 (0.33)0.38 (0.33)0.37 (0.26)
Follow-up period (years), mean (SD)2.2 (2.2)2.1 (1.6)2.5 (2.2)3.5 (2.6)3.2 (2.3)3.3 (2.5)2.9 (2.3)
CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; hrHPV, high-risk human papillomavirus; SD, standard deviation. Cytological and histological results were combined to classify the results into the following diagnoses: normal, CIN1, and CIN2. HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 were classified as hrHPVs. Of these, HPV 16, 18, 52, and 58 were categorized separately. hrHPVs other than HPVs 16, 18, 52, and 58 were classified as “other hrHPVs.” Patients who were not infected with any hrHPVs were referred to as “no hrHPVs” patients. In cases of observed coinfection with different HPV genotypes, it was possible to include the same patient in the summary statistics of multiple HPV categories.
Table 2. Summary of the transitions from each diagnosis of cervical epithelial lesions for the six human papillomavirus (HPV) categories.
Table 2. Summary of the transitions from each diagnosis of cervical epithelial lesions for the six human papillomavirus (HPV) categories.
Diagnosis at tth Visit
Diagnosis at (t-1)th VisitHPV CategoryNormalCIN1CIN2CIN3Cancer
NormalHPV 16206 (84.4)13 (5.3)21 (8.6)4 (1.6)0 (0.0)
HPV 1889 (81.6)12 (11.0)8 (7.3)0 (0.0)0 (0.0)
HPV 52277 (75.8)54 (14.7)32 (8.7)2 (0.5)0 (0.0)
HPV 58230 (80.4)33 (11.5)21 (7.3)2 (0.6)0 (0.0)
Other hrHPVs611 (86.1)72 (10.1)23 (3.2)3 (0.4)0 (0.0)
No hrHPVs1289 (90.2)109 (7.6)26 (1.8)3 (0.2)1 (0.0)
CIN1HPV 1629 (28.9)34 (34.0)35 (35.0)2 (2.0)0 (0.0)
HPV 1818 (38.2)19 (40.4)8 (17.0)2 (4.2)0 (0.0)
HPV 5280 (35.0)90 (39.4)53 (23.2)5 (2.1)0 (0.0)
HPV 5851 (31.6)68 (42.2)40 (24.8)2 (1.2)0 (0.0)
Other hrHPVs132 (40.7)143 (44.1)45 (13.8)4 (1.2)0 (0.0)
No hrHPVs203 (54.5)132 (35.4)34 (9.1)3 (0.8)0 (0.0)
CIN2HPV 1631 (12.1)37 (14.4)147 (57.4)40 (15.6)1 (0.3)
HPV 1810 (12.9)8 (10.3)51 (66.2)8 (10.3)0 (0.0)
HPV 5241 (13.8)53 (17.9)168 (56.9)33 (11.1)0 (0.0)
HPV 5832 (10.2)45 (14.4)210 (67.5)24 (7.7)0 (0.0)
Other hrHPVs49 (16.7)52 (17.8)166 (56.8)25 (8.5)0 (0.0)
No hrHPVs58 (27.2)31 (14.5)114 (53.5)10 (4.6)0 (0.0)
CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; hrHPV, high-risk human papillomavirus. Values are the number (percentage) of transitions observed from prior diagnosis to current diagnosis for each HPV category. Cytological and histological results were combined to classify the results into the following diagnoses: normal, CIN1, CIN2, CIN3, and cervical cancer. HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 were classified as hrHPVs. Of these, HPV 16, 18, 52, and 58 were categorized separately. hrHPVs other than HPVs 16, 18, 52, and 58 were classified as “other hrHPVs.” Patients who were not infected with any hrHPVs were referred to as “no hrHPVs” patients. In cases of observed coinfection with different HPV genotypes, it was possible to include the same patient in multiple HPV categories.
Table 3. Predicted 2-year transition probabilities from states to states and their 95% confidence intervals (CIs) for the six HPV categories.
Table 3. Predicted 2-year transition probabilities from states to states and their 95% confidence intervals (CIs) for the six HPV categories.
Current State State after 2 Years
HPV CategoryNormalCIN1CIN2CIN3/Cancer
NormalHPV 160.598 (0.506–0.684)0.099 (0.074–0.128)0.169 (0.127–0.215)0.132 (0.090–0.183)
HPV 180.610 (0.479–0.719)0.156 (0.109–0.215)0.156 (0.093–0.230)0.076 (0.033–0.148)
HPV 520.533 (0.474–0.593)0.189 (0.162–0.219)0.180 (0.146–0.216)0.096 (0.070–0.130)
HPV 580.559 (0.484–0.627)0.171 (0.140–0.205)0.206 (0.162–0.255)0.062 (0.041–0.089)
Other hrHPVs0.723 (0.676–0.766)0.155 (0.132–0.182)0.085 (0.066–0.108)0.034 (0.023–0.050)
No hrHPVs0.838 (0.814–0.861)0.105 (0.090–0.121)0.042 (0.032–0.054)0.012 (0.007–0.020)
CIN1HPV 160.434 (0.349–0.512)0.089 (0.067–0.115)0.175 (0.134–0.223)0.300 (0.225–0.378)
HPV 180.535 (0.396–0.652)0.146 (0.100–0.207)0.172 (0.102–0.257)0.146 (0.069–0.267)
HPV 520.473 (0.413–0.529)0.178 (0.152–0.208)0.183 (0.150–0.221)0.164 (0.122–0.219)
HPV 580.469 (0.399–0.535)0.165 (0.135–0.197)0.239 (0.192–0.291)0.126 (0.084–0.181)
Other hrHPVs0.656 (0.606–0.702)0.156 (0.133–0.181)0.102 (0.079–0.128)0.084 (0.058–0.119)
No hrHPVs0.808 (0.781–0.835)0.107 (0.091–0.123)0.049 (0.038–0.065)0.034 (0.021–0.054)
CIN2HPV 160.335 (0.266–0.404)0.079 (0.059–0.101)0.165 (0.121–0.218)0.418 (0.330–0.512)
HPV 180.373 (0.245–0.501)0.119 (0.074–0.178)0.186 (0.099–0.302)0.320 (0.178–0.507)
HPV 520.381 (0.324–0.434)0.156 (0.129–0.184)0.175 (0.138–0.216)0.286 (0.220–0.367)
HPV 580.356 (0.291–0.419)0.150 (0.122–0.181)0.260 (0.209–0.319)0.232 (0.167–0.307)
Other hrHPVs0.518 (0.453–0.571)0.146 (0.122–0.169)0.117 (0.089–0.148)0.218 (0.159–0.299)
No hrHPVs0.706 (0.643–0.749)0.106 (0.090–0.123)0.063 (0.045–0.089)0.124 (0.079–0.191)
CIN, cervical intraepithelial neoplasia; HPV, human papillomavirus; hrHPV, high-risk human papillomavirus. Values are the predicted probabilities (95% confidence intervals) of transitions from the current state to the state after two years for each HPV category. Cytological and histological results were combined to classify the results into the following diagnoses: normal, CIN1, CIN2, CIN3, and cervical cancer. HPVs 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, and 68 were classified as hrHPVs. Of these, HPV 16, 18, 52, and 58 were categorized separately. hrHPVs other than HPVs 16, 18, 52, and 58 were classified as “other hrHPVs.” Patients who were not infected with any hrHPVs were referred to as “no hrHPVs” patients. In cases of observed coinfection with different HPV genotypes, it was possible to include the same patient in multiple HPV categories. We used the continuous-time multistate Markov model to estimate the prognosis of each patient. We defined four states: normal (state 1), cervical intraepithelial neoplasia 1 (CIN1, state 2), CIN2 (state 3), and CIN3/cancer (state 4). Arrows in Figure 2 specify possible transitions between the states defined in our model; all transitions between adjacent states, except the backward transition from CIN3/cancer to CIN2, were allowed. CIN3/cancer was the absorbing state. We truncated observations after the diagnosis of CIN3 or cancer.

Share and Cite

MDPI and ACS Style

Taguchi, A.; Hara, K.; Tomio, J.; Kawana, K.; Tanaka, T.; Baba, S.; Kawata, A.; Eguchi, S.; Tsuruga, T.; Mori, M.; et al. Multistate Markov Model to Predict the Prognosis of High-Risk Human Papillomavirus-Related Cervical Lesions. Cancers 2020, 12, 270. https://doi.org/10.3390/cancers12020270

AMA Style

Taguchi A, Hara K, Tomio J, Kawana K, Tanaka T, Baba S, Kawata A, Eguchi S, Tsuruga T, Mori M, et al. Multistate Markov Model to Predict the Prognosis of High-Risk Human Papillomavirus-Related Cervical Lesions. Cancers. 2020; 12(2):270. https://doi.org/10.3390/cancers12020270

Chicago/Turabian Style

Taguchi, Ayumi, Konan Hara, Jun Tomio, Kei Kawana, Tomoki Tanaka, Satoshi Baba, Akira Kawata, Satoko Eguchi, Tetsushi Tsuruga, Mayuyo Mori, and et al. 2020. "Multistate Markov Model to Predict the Prognosis of High-Risk Human Papillomavirus-Related Cervical Lesions" Cancers 12, no. 2: 270. https://doi.org/10.3390/cancers12020270

APA Style

Taguchi, A., Hara, K., Tomio, J., Kawana, K., Tanaka, T., Baba, S., Kawata, A., Eguchi, S., Tsuruga, T., Mori, M., Adachi, K., Nagamatsu, T., Oda, K., Yasugi, T., Osuga, Y., & Fujii, T. (2020). Multistate Markov Model to Predict the Prognosis of High-Risk Human Papillomavirus-Related Cervical Lesions. Cancers, 12(2), 270. https://doi.org/10.3390/cancers12020270

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

Article Metrics

Back to TopTop