Skip to main content

HPV E6/E7 mRNA in situ hybridization in endocervical adenocarcinoma: implications for prognosis and diagnosis



Human papillomavirus (HPV) E6/E7 mRNA in situ hybridization (HPV E6/E7 RNAscope) appears to be a sensitive and specific technique in detecting transcriptionally active HPV. We aimed to examine the diagnostic utility of this technique in endocervical adenocarcinoma (ECA), to explore the prognostic factors for ECA patients and develop a clinically useful nomogram based on clinicopathological parameters to predict it.


We retrospectively analyzed 200 patients with ECA who had undergone surgery at Sun Yat-sen University Cancer Center from 2010 and 2014. The diagnostic performance of HPV E6/E7 RNAscope were evaluated by receiver operating characteristic (ROC) curve. A prognostic nomogram model including HPV E6/E7 RNAscope was generated based on multivariate Cox regression analysis, then we compared the predictive accuracy of the prognostic model with FIGO staging and treatment using concordance index (C-index), time-dependent ROC (tdROC), and decision curve analysis (DCA).


The sensitivity and specificity of HPV E6/E7 RNAscope for distinguishing HPV-associated adenocarcinoma (HPVA) from non-HPV-associated adenocarcinoma (NHPVA) in the whole cohort were 75.8% and 80%, respectively. According the univariate analysis and multivariate logistic regression analysis, age, lymphovascular invasion (LVI), lymph node involvement (LNI), and HPV E6/E7 RNAscope were valuable predictive factors for OS. These parameters were incorporated into the nomogram model (nomogram A) compared with FIGO stage and treatment. The C-index of nomogram A for predicting OS was 0.825, which was significantly higher than FIGO stage (C-index = 0.653, p = 0.002) and treatment (C-index = 0.578, p < 0.001).


 HPV E6/E7 RNAscope is highly specific for ECA, and the 4-variable nomogram showed more accurate prognostic outcomes in patients with ECA.


Endocervical adenocarcinoma (ECA) accounts for 15%–20% of all cervical carcinomas. Studies have shown that ECA is increasingly reported in young women [1], with a projected 5-year overall survival (OS) rate of only 20.3% (95% confidence interval [CI]: 14.2–27.1%) compared to the squamous type (31.3%, 95%CI: 29.2–33.3%) [2]. The International Federation of Gynecology and Obstetrics (FIGO) staging system is the gold standard for predicting outcomes in patients with ECA, and the mainstay treatment of advanced ECA is concurrent chemoradiotherapy. However, ECA patients with the same FIGO stage can have marked heterogeneities in their outcomes [3, 4]. The new FIGO 2018 amendments put forward by the gynecological oncology committee agree that staging is an ongoing process informed by data on prognosis and survival, and therefore, treatment should be individualized and not merely based on staging given the variability of the resource across regions. FIGO 2018 also recommends any imaging modality and/or pathological findings like lymphovascular invasion (LVI) for allocating staging [3].

A new classification—the International Endocervical Adenocarcinoma Criteria and Classification (IECC) has replaced the 2014 WHO classification and categorizes ECA based on the morphological features linked to the etiology (i.e., HPV infection), thus grouping them into HPV-associated (HPVA) and non-HPVA (NHPVA) types [4, 5]. Comparing HPVA with NHPVA types reveals essential differences in tumor biology and patient survival, with NHPVA-type tumors significantly associated with worse outcomes [6, 7]. It is apparent that a classification based on the pathogenesis and other clinicopathological factors may be more clinically useful and reproducible than the current FIGO scheme. The other shortcoming of FIGO is that it doesn't stratify patients based on histology and HPV status. Hence, novel prognostic determinants to complement the FIGO staging system is needed. The validity of this new classification is supported by HPV status but with limited clinical data[8, 9].

In the present studies, the following tests were used to detect HPV infection: HPV DNA, genotype, RNAscope, and immunochemistry (IHC) against p16 protein. Plasma HPV DNA levels can be determined using PCR and is a potential marker for cervical cancer [10]. The other 3 assays can be conducted on formalin-fixed paraffinembedded (FFPE) samples. Immunostaining for p16 is a cost-effective and highly sensitive marker but with a low specificity [11]. PCR for HPV genotyping is the gold-standard assay to diagnose active HPV infection [12]. Chromogenic in situ hybridization against RNA can identify HPV in situ via microscopic observation [13]. The technique is expensive but can be used on FFPE samples [14]. In situ detection of HPV E6/E7 mRNA using the RNAscope also appears to be a sensitive and a specific method; hence, we aimed to investigate the potential prognostic utility of this technique in ECAs.

Numerous controversies surround the FIGO staging system, wherefore ECAs are morphologically stratified into A, B and C subgroups using the Silva based on their LVI or lymph node involvement (LNI) status [15]. Tumors characterized by well-demarcated glands with rounded contours are designated as pattern A; tumors with early or limited, localized destructive invasion or inflamed stroma adjacent to an intact gland are designated pattern B, while tumors with more aggressive features, including diffusely infiltrative glands, usually associated with extensive, diffuse desmoplastic response are designated pattern C. Therefore, clinicians could provide precision treatment to patients with pattern C tumors. However, stratifying patients based on a single biomarker, such as LVI or LNI alone, can lead to misclassification of tumors and poor patients' prognostication, hence the two parameters should be combined [16]. There is currently no such available nomogram model for patients with ECA that combines HPV, LVI and LNI status, And other survival prediction models that follow the FIGO staging system have been insufficient.

In the present study, we aim to validate the utility of the novel RNA scope method for detecting HPV RNA in tissue samples for proper initial risk stratification of patients. We also built a nomogram that will compliment FIGO staging for prognostic classification of EAC. We hope the link between these two modalities will bring new insights into the clinical care of EAC patients.

Patients and methods


The present retrospective study included 200 patients with pathologically diagnosed ECA who had been treated at the Sun Yat-sen University Cancer Center from January 2010 to December 2014. We enrolled patients with confirmed pathological diagnoses of ECA with available complete clinical data and test results. We excluded patients with systemic metastasis or with synchronous primaries and any patient with a prior cancer history. The last follow-up period was June 2020. The histological diagnoses were based on our proposal for a new ECA classification: the IECC 2017 [9], and patients were grouped into HPVA and NHPVA ECA based on HE staining and microscopic morphology.

Immunohistochemistry (IHC)

A tissue microarray consisting of the ECA and adjacent non-tumorous tissue was constructed using a tissue array instrument (Minicoreexcilone; Minicore, UK). For each patient, hematoxylin and eosin (H&E) stained slides were examined, and at least two areas from different regions were marked for sampling. A 1.0 mm diameter tissue core was punched from the marked fields and re-embedded. FFPE ECA sections were dewaxed in xylene and graded alcohols, hydrated, and washed in PBS. After pretreatment in a microwave oven, endogenous peroxidase was blocked using 3% hydrogen peroxide in methanol for 20 min. This was followed by avidin–biotin blocking using a kit (DAKO, Germany). The slides were then incubated with antibodies for p16 (Roche, Germany), MLH1 (Roche [M1], Germany), PMS2 (Dako [EP51], Germany), MSH2 (ZA0622, Zhongshan, China), MSH6 (Roche [SP93], Germany), and Ki-67 (ZA0502, Zhongshan, China) overnight at 4 °C, washed in PBS, and incubated with biotinylated goat anti-rabbit/mouse antibodies (DAKO, Germany). The slides were developed using diaminobenzidine (DAB) and counterstained with hematoxylin. Two independent, experienced pathologists evaluated the staining; and p16 was interpreted as positive if diffuse, block-like staining was found in all cores. No staining or patchy staining was interpreted as unfavorable. MSH2/MSH6/PMS2/MLH1 was interpreted as positive if ≥ 1% of the tumor cell nuclei were positive. Representative images were show in Additional file 1: Fig. S1.


DNA load of High-risk-HPV (HR-HPV) in plasma sample was evaluated by Digene second generation hybrid capture (HC2) DNA test and in vitro acid hybridization assay with signal amplification using microplate chemiluminescence for HR-HPV DNA detection. A commercial kit detecting high risk was used. The alkaline solution was used to denature the double-stranded DNA, and the liberated single-strand DNA was combined with the RNA probe mix. RNA–DNA hybrids were then transferred to a capture microplate coated with goat polyclonal for immobilization. Subsequently, alkaline phosphatase-conjugated antibodies were added to the RNA–DNA hybrids. After washing, the chemiluminescent substrate was added for signal amplification. Results were analyzed by DIGENE DML 2000 Software. DNA load was expressed by the unit of Relative Light Units/Cutoff Value (RLU/CO), representing the ratio of the light emission of the sample to the average of three positive control samples. HR-HPV DNA load of 1.0 (pg/ml) or more was defined as HR-HPV DNA positive.

HPV genotype

Tissues were scraped by coverslip as per the labeled region on the H&E slide and transferred into a clean 1.5 ml centrifuge tube. DNA was extracted using a QIAGEN DNA FFPE Kit (Qiagen, Dusseldorf, Germany) according to the manufacturer's instructions and was quantified with a Nano-Drop2000 (NanoDrop Technologies, Wilmington, DE, USA). HPV was detected by the Human papillomavirus Genotyping Test Kit (Anhui Dajian, Guangzhou, China) using the BioRad CFX96 real-time quantitative PCR instrument. The Roche Cobas 4800 system (Roche, Pleasanton, CA) was used to detect HPV and evaluate the presence of the following 14 types of HPV DNA: 16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68.

HPV E6/E7 mRNA in situ hybridization

The stained slides of each specimen were examined using the RNAscope scoring system, as described in previous studies [17, 18]. High-risk HPV subtypes were evaluated on all ECAs in the tissue microarray that had sufficient tissue to allow scoring (n = 200). HPV fluorescent in situ hybridization (FISH) was performed using a chromogen and the RNAscope system (Advanced Cell Diagnostics; catalog No. 312598). The RNAscope probe “HPV HR18” contains probes that target E6 and E7 mRNA in the following high-risk subtypes: HPV16, 18, 26, 31, 33, 35, 39, 45, 51, 52, 53, 56, 58, 59, 66, 68, 73, and 82. RNAscope results were classified into five degrees based on the following scoring guidelines: score 0, no staining or less than one dot in every ten cells (visible at 40 × magnification); score 1, 1–3 dots per cell (visible at 20–40 × magnification); score 2, 4–10 dots per cell, with very few dot clusters (visible at 20–40 × magnification); score 3, > 10 dots per cell, with < 10% of positive cells having dot clusters (visible at 20 × magnification); score 4, > 10 dots per cell with > 10% of positive cells having dot clusters (visible at 20 × magnification). Cases with RNAscope score ≥ 1 were identified as positive.

Data collection

Among the patients enrolled, the following variables were recorded: age, the modality of treatment (Surgery with or without chemotherapy/radiotherapy), HPV DNA, HPV subtype, HPV E6/E7 RNAscope, histological type, tumor size, differentiation, LVI, invasion level of the uterine cervix, LNI, parametrium invasion, surgical margin, MMR status, p16, Ki-67, and FIGO stage.

Statistical analysis

Statistical investigations were conducted using SPSS 19.0 (IBM Corp., Armonk, NY, USA) and R 3.4.0 ( Continuous data were expressed as means ± standard deviations and were analyzed using Student’s t-test or Mann–Whitney U test, as appropriate based on the results of the Kolmogorov–Smirnov test. Categorical variables were presented using numbers (percentages) and were analyzed with Fisher’s exact test. Survival analysis was conducted utilizing the Kaplan–Meier curves. Univariate and multivariate Cox proportional hazard regression were used to explore the risk factor for OS. The best cutoff values were determined using X-tile [19]. The best cutoff values were determined using X-tile with the following results: age, 37 years; RNAscope, 3.3 (Additional file 1: Fig. S2A); tumor size, 4.5 cm, and Ki-67 (12.5%). Receiver operating characteristic curves (ROC) were generated to compare the diagnostic performance of different detection method for ECA patients. The model's predictive efficiency and discriminative capability were defined by a concordance index (C-index) and a calibration curve. A two-sided p-value < 0.05 was considered statistically significant.


Patients’ characteristics and survival

We enrolled 200 patients with ECA for analysis, of which 185 were HPVA and 15 NHPVA (Table 1). Across both groups, the median follow-up period was 68 months (range: 1.5–125.1 months). The 1, 3, and 5-year OS rates were 99.5%, 90.5%, and 85.0%, respectively.

Table 1 Patient demographics and clinical charateristics

Diagnostic performance of HPV E6/E7 mRNA in situ hybridization in ECA compared to other assays.

The RNAscope scores of 0 to 4 are shown in Fig. 1A and B. The diagnostic implications of the HPV E6/E7 RNAscope scores were evaluated in 200 ECA samples that included both HPVA and NHPVA cases (Additional file 2: Table S1). The positive rates of HPV DNA, p16 IHC, HPV genotype, and HPV E6/E7 RNAscope across all ECA cases were 76.4%, 59.0%, 62.5%, and 72.0%, respectively (Fig. 1C and Additional file 1: Fig. S3). HPV RNAscope and other assays were closely related (Additional file 2: Table S2). ROC curves suggested that HPV DNA and HPV E6/E7 RNAscope showed similar results in terms of distinguishing HPVA from NHPVA (area under the curve [AUC] = 0.802, sensitivity = 80.5%, specificity = 80% vs. AUC = 0.799, sensitivity = 75.8%, specificity = 80%, respectively), outperforming both p16 IHC (AUC = 0.751, sensitivity = 60.2%, specificity = 90%) and HPV genotype (AUC = 0.566, sensitivity = 63.3%, specificity = 50%; Fig. 1D and Table 2).

Fig. 1
figure 1

Representative images of human papillomavirus (HPV) E6/E7 RNAscope in paraffin-embedded endocervical adenocarcinoma (ECA) samples. A HPV mRNA was detected using RNAscope in ECA (T) and non-tumorous tissues adjacent to ECA (N). Representative images of scores 0 (T0), 1 (T1), 2 (T2), 3 (T3) and 4 (T4) are shown. B The proportion of RNAscope scores in ECA tissues. C Positive rates of HPV by PCR, RNAscope, HPV genotype and p16 IHC. (D) Receiver operating characteristics curve for HPV, detected by PCR, immunohistochemistry, genotyping and RNAscope for HPV-associated vs. non-HPV-associated types

Table 2 Diagnostic performances of studied testing for ECA patients

The correlations among serum HPV DNA level, p16 IHC, HPV genotype, and HPV E6/E7 RNAscope were also studied (Additional file 2: Table S2). A significantly higher proportion of ECA cases were positive for p16 IHC, HPV genotype, and HPV E6/E7 RNAscope in subgroups with HPV DNA levels that differed by tenfold. In patients with serum HPV DNA levels of 1–100 pg/mL, the p16 IHC positivity was similar to that of HPV genotype (60%), while those with serum HPV DNA levels of 100–10,000 pg/mL, HPV genotype positivity was similar to that of HPV E6/E7 RNAscope (93.8%). The HPV E6/E7 RNAscope positivity for the diagnosis of ECA was much higher compared to p16 IHC and HPV genotype (Additional file 1: Fig. S4A). We also elicited that HPV E6/E7 RNAscope was positive in almost all cases that were p16 IHC positive (88.5% of all cases, 89.5% of HPVA cases), only that more cases were negative by IHC and serum levels (Additional file 1: Fig. S4B). In different HPV subtypes, the positivity of RNAscope for the diagnosis of ECA was much higher compared to that of p16 IHC, but lower than HPV DNA (Additional file 1: Fig. S4C). The HPV subtypes in patients with HPVA and NHPVA are listed in Additional file 2: Table S3. In cases with RNAscope scores of 3 and 4, the HPV DNA positivity was higher than both p16 IHC and HPV genotype (Additional file 1: Fig. S4D).

Multivariate analysis of overall survival

Univariate analysis indicated that the following variables were related to OS in patients with ECA: age (p = 0.002), HPV E6/E7 RNAscope (p = 0.006), tumor size (p = 0.005), FIGO stage (p < 0.001), histological type (p = 0.004), LVI (p < 0.001), invasion level of uterine cervix (p < 0.001), LNI (p < 0.001), parametrium invasion (p < 0.001), surgical margin (p = 0.006), chemotherapy (p = 0.009), radiotherapy (p < 0.001), p16 (p = 0.039), and Ki-67 (p = 0.025). In a multivariate analysis, the following items remained independently prognostic: age (HR = 0.250, 95% confidence interval [CI]: 0.099–0.632, p = 0.003), HPV E6/E7 RNAscope (HR = 0.240, 95% CI: 0.093–0.616, p = 0.003), LVI (HR = 1.924, 95% CI: 1.295–2.857, p = 0.001), and LNI (HR = 3.047, 95% CI: 1.183–7.849, p = 0.021). The results are shown in Table 3, with a forest plot in Fig. 2A. Kaplan–Meier analysis showed significant diversity (Fig. 2B–E).

Table 3 Univariate and multivariate Cox proportional hazards regression analysis for OS
Fig. 2
figure 2

Forest plot and Kaplan–Meier curves for overall survival (OS) in patients with endocervical adenocarcinoma (ECA). A Forest plot showed the hazard ratio and 95% confidence interval for OS according to the Cox proportional hazards regression analysis in patients with ECA. BE Age, human papillomavirus E6/E7 RNAscope, lymphovascular invasion, and lymph node involvement in patients with ECA in the whole cohort are plotted as a distribution

Construction of the prognostic nomogram model and comparison of predictive accuracy between Nomogram model, FIGO stage and treatment.

Two different nomograms that predict the OS of ECA patients were built and compared for efficiency. Nomogram A included age, HPV E6/E7 RNAscope, LVI, and LNI, while nomogram B included age, HPV E6/E7 RNAscope, LVI, LNI, FIGO stage, and treatment. The results of time-dependent ROC curve for OS showed that AUCs of nomogram A and FIGO stage were higher than treatment (Fig. 3A). The resulting variables from the Cox proportional analysis were used to build the prognostic nomograms for OS (Fig. 3B and C). Each prognostic factor within the nomogram was assigned a point. By sum of the total points from all variables combined with the location the total point scale allowed us to obtain the probabilities of the outcomes by drawing a vertical line towards the axis labeled “1-, 3-, 5-Year survival probability” (Fig. 3B and C). The results of the comparison of efficiency between our nomogram and the conventional systems are shown in Table 4. In the two nomograms, no significant difference was observed (C-index: 0.825, 95% CI = 0.754–0.896 vs. 0.836, 95% CI = 0.771–0.902). The C-index of nomogram A (0.825, 95% CI = 0.754–0.896) was better than those of the FIGO system (0.653, 95% CI = 0.567–0.740) and treatment (0.578, 95% CI = 0.506–0.651). In Fig. 4A and B, the calibration plot for the OS rates for 1, 3, and 5-years was in line with both the nomogram and the actual observation. Besides, our model seemed to have a higher prediction accuracy (Fig. 4C).

Fig. 3
figure 3

Time-dependent ROC curve analysis and nomogram model for overall survival (OS) in patients with endocervical adenocarcinoma ECA. A Time-dependent receiver operating characteristic curves showing the sensitivity and specificity of the nomogram model for predicting OS. B, C The nomogram models A and B were used to summate the points identified on the points scale for each variable, indicating the probability of 1-, 3-, and 5-year survival

Table 4 The C-index of our model, FIGO stage, Treatment for prediction of OS in the ECA
Fig. 4
figure 4

The calibration curves and decision curve analysis for predicting patients’ overall survival (OS). The calibration curves for predicting OS at 1, 3, and 5 years in all cases. A, B Nomograms A and B, with the model-predicted OS plotted on the x-axis and actual OS plotted on the y-axis. Closer alignment with the diagonal line represents better estimation. C Decision curve analysis for 1-, 3-, and 5-year survival predictions

Risk stratification of prognosis

All patients were divided into a low (≤ 15.75) and a high-risk group (> 15.75) for OS (Additional file 1: Fig. S2B). The 1,3, and 5-year OS rates were 94.2%, 84.3%, 70.3%, respectively, in the low-risk group and 71.4%, 46.4%, 10.7%, respectively, in the high-risk group (Table 5). Moreover, significant differences in OS were noted between patients with stage I/II ECA and those with III/IV ECA (Fig. 5). Each risk subgroup represented a distinct prognosis, and our system accurately separated the OS rates of the two subgroups.

Table 5 Point assignment and prognostic score of the nomogram model
Fig. 5
figure 5

Graphs showing the Kaplan–Meier curves of both groups based on the predictors from the nomogram model. A stage I. B stage II. C stages III and IV


Our study is the first of its kind to describe a prognostic nomogram model for invasive ECA in a large Chinese series by reporting HPV E6/E7 RNAscope using a new RNA ISH assay that recognizes 18 high-risk HPV types. The most important finding of our research are as follows: (1) the validity of the HPVA and NHPVA categories were supported by p16 immunophenotyping and HPV status; (2) HPV E6/E7 RNAscope is more sensitive and specific than p16 and HPV genotype in identifying HPVA; (3) we established a prognostic nomogram model based on ECA pathogenesis that adjusted for age, LVI, HPV E6/E7 RNAscope, and LNI, which had a greater predictive efficiency than the current conventional staging system; (4) in the future, we believe, based on our model, ECA patients could be divided into low and high-risk groups.

HPV E6/E7 RNAscope measurement is not routinely available in most hospitals; hence the methodology cannot be extrapolated to standard practice. About 94% of adenocarcinomas are associated with HPV, particularly high risk strains like 16 and 18 [5], which is consistent with our results. In the present study, 87% (160/184) of HPVA samples overexpressed p16 or were HPV positive by RNAscope, validating the IECC criteria and previous studies [9, 20]. In terms of the IECC criteria, the HPV RNAscope is more specific for clinical practice, and it is more robust than p16 IHC, HPV DNA, and HPV genotype in identifying HPVA. Despite this, occasional p16 and/or HPV negative HPVAs were also identified. When outlying cases were excluded from the statistical analysis, 25 of the 184 HPVA cases were p16- and HPV-positive. All HPVA cases in the present study were systematically reviewed. The specimens were more than 5 years old and may not have been optimal tissues for the performance of HPV RNAscope and p16. A recent study reported that the HPV associated ECA might represent unusual morphological variants of gastric-type carcinoma [21]. Furthermore, rare HPV genotypes not included in the RNA ISH probe set may be responsible for the negative HPV results. There are precedents in the literature for p16- and HPV-associated neoplasia. Notably, p16-positive invasive squamous carcinomas of the cervix have been reported [22] to be closely associated with methylation-induced inactivation of the p16 gene and allelic loss of p16 [23, 24]. Other studies have reported that the HPV genome can be differentially expressed with varying heterogeneity between the primary and the metastatic sites suggesting that some HPV-related usual-type adenocarcinomas could have detectable HPV in multiple sections of primary carcinoma; hence all metastatic sites should be tested [20, 25, 26].

Clinicopathological variables carry a prognostic significance and usually affect clinical management [27, 28]. Depending on variables such as LVI and desire for fertility, histological-based risk evaluation should be considered during tissue sampling at initial diagnosis [15, 29]. The National Comprehensive Cancer Network guidelines recommend chemotherapy or radiotherapy as the standard treatment for patients with advanced ECA, hence, the performance of nomograms must be examined separately in patients treated using these therapies. The C-indices of nomograms A and B for predicting OS in treated patients were significantly higher than those of conventional classification, indicating that nomograms still have significant clinical value in patients with ECA. With the addition of the FIGO stage and treatment into nomogram B, the added values of these parameters over nomogram A were 0.836 and 0.825, respectively, across the entire population. There were no significant differences between nomograms A and B, probably due to the small sample size of patients with advanced-stage ECA.

The present study had several limitations. Firstly, there may have been selection bias in our patient cohort. Secondly, we lacked data on the impact of our nomogram on other prognostic endpoints, such as disease-free survival (DFS) prediction. Thirdly, in situ measurement of HPV RNA using RNAscope is not easily available and not globally standardized. Fourthly, our sample size was relatively small and from only one center. In order to validate our results, a study with a larger sample with multicenter data is required. Lastly, it remains unclear whether our nomogram can be applied to ECA patients with advanced-stage that is, stage III and IV.


We generated new nomograms that incorporate HPV E6/E7 mRNA in situ hybridization and patients’ clinical characteristics in order to prognosticate the OS rate in ECA patients. Our simple and explicit nomograms have a good clinical application value, with a good diagnostic discernment and calibration ability. They may be a useful tool for assessing the prognosis and the management of ECA patients.

Availability of data and materials

The authenticity of this article has been validated by uploading the key raw data onto the Research Data Deposit public platform ( with the approval RDD number RDDA2020001638.



Human papillomavirus

HPV E6/E7 RNAscope:

Human papillomavirus E6/E7 mRNA in situ hybridization;


Endocervical adenocarcinoma


Concordance index


HPV-associated adenocarcinoma


Non-HPV-associated adenocarcinoma


Overall survival


Lymphovascular invasion


Lymph node involvement


  1. Ward KK, Shah NR, Saenz CC, McHale MT, Alvarez EA, Plaxe SC. Changing demographics of cervical cancer in the United States (1973–2008). Gynecol Oncol. 2012;126(3):330–3.

    Article  Google Scholar 

  2. Galic V, Herzog TJ, Lewin SN, Neugut AI, Burke WM, Lu YS, Hershman DL, Wright JD. Prognostic significance of adenocarcinoma histology in women with cervical cancer. Gynecol Oncol. 2012;125(2):287–91.

    Article  Google Scholar 

  3. Bhatla N, Berek JS, Cuello Fredes M, Denny LA, Grenman S, Karunaratne K, Kehoe ST, Konishi I, Olawaiye AB, Prat J, et al. Revised FIGO staging for carcinoma of the cervix uteri. Int J Gynaecol Obstet. 2019;145(1):129–35.

    Article  Google Scholar 

  4. Stolnicu S, Hoang L, Chiu D, Hanko-Bauer O, Terinte C, Pesci A, Aviel-Ronen S, Kiyokawa T, Alvarado-Cabrero I, Oliva E, et al. Clinical Outcomes of HPV-associated and Unassociated Endocervical Adenocarcinomas Categorized by the International Endocervical Adenocarcinoma Criteria and Classification (IECC). Am J Surg Pathol. 2019;43(4):466–74.

    Article  Google Scholar 

  5. Rj K, Ml C, Cs H. RH Y: WHO Classification of Tumors of Female Reproductive Organs. 4th ed. Lyon: IARC, WHO Press; 2014.

    Google Scholar 

  6. Byun JM, Cho HJ, Park HY, Kim YN, Lee KB, Sung MS, Jeong CH, Jeong DH. Clinical significance of the pattern-based classification in endocervical adenocarcinoma, usual and variants. Int J Clin Oncol. 2019;24(10):1264–72.

    Article  Google Scholar 

  7. Karamurzin YS, Kiyokawa T, Parkash V, Jotwani AR, Patel P, Pike MC, Soslow RA, Park KJ. Gastric-type endocervical adenocarcinoma: an aggressive tumor with unusual metastatic patterns and poor prognosis. Am J Surg Pathol. 2015;39(11):1449–57.

    Article  Google Scholar 

  8. Stolnicu S, Barsan I, Hoang L, Patel P, Terinte C, Pesci A, Aviel-Ronen S, Kiyokawa T, Alvarado-Cabrero I, Oliva E, et al. Stromal invasion pattern identifies patients at lowest risk of lymph node metastasis in HPV-associated endocervical adenocarcinomas, but is irrelevant in adenocarcinomas unassociated with HPV. Gynecol Oncol. 2018;150(1):56–60.

    Article  CAS  Google Scholar 

  9. Stolnicu S, Barsan I, Hoang L, Patel P, Terinte C, Pesci A, Aviel-Ronen S, Kiyokawa T, Alvarado-Cabrero I, Pike MC, et al. International endocervical adenocarcinoma criteria and classification (IECC): a new pathogenetic classification for invasive adenocarcinomas of the endocervix. Am J Surg Pathol. 2018;42(2):214–26.

    Article  Google Scholar 

  10. Cheung TH, Yim SF, Yu MY, Worley MJ Jr, Fiascone SJ, Chiu RWK, Lo KWK, Siu NSS, Wong MCS, Yeung ACM, et al. Liquid biopsy of HPV DNA in cervical cancer. J Clin Virol. 2019;114:32–6.

    Article  CAS  Google Scholar 

  11. McCluggage WG, Jenkins D. p16 immunoreactivity may assist in the distinction between endometrial and endocervical adenocarcinoma. Int J Gynecol Pathol. 2003;22(3):231–5.

    Article  CAS  Google Scholar 

  12. Kenemans P. HPV genotype as a prognostic factor for progression to cervical carcinoma in young women. Eur J Obstet Gynecol Reprod Biol. 1994;55(1):24–5.

    Article  CAS  Google Scholar 

  13. Schache AG, Liloglou T, Risk JM, Jones TM, Ma XJ, Wang H, Bui S, Luo Y, Sloan P, Shaw RJ, et al. Validation of a novel diagnostic standard in HPV-positive oropharyngeal squamous cell carcinoma. Br J Cancer. 2013;108(6):1332–9.

    Article  CAS  Google Scholar 

  14. Mirghani H, Casiraghi O, Amen F, He M, Ma XJ, Saulnier P, Lacroix L, Drusch F, Ben Lakdhar A, Saint Guily JL, et al. Diagnosis of HPV-driven head and neck cancer with a single test in routine clinical practice. Mod Pathol. 2015;28(12):1518–27.

    Article  CAS  Google Scholar 

  15. De Vivar A, Roma AA, Park KJ, Alvarado-Cabrero I, Rasty G, Chanona-Vilchis JG, Mikami Y, Hong SR, Arville B, Teramoto N, et al. Invasive endocervical adenocarcinoma: proposal for a new pattern-based classification system with significant clinical implications: a multi-institutional study. Int J Gynecol Pathol. 2013;32(6):592–601.

    Article  Google Scholar 

  16. Roma AA, Park KJ, Xie H, De Vivar AD, Alvarado-Cabrero I, Rutgers JKL, Barbuto D, Silva EG. Role of lymphovascular invasion in pattern C invasive endocervical adenocarcinoma. Am J Surg Pathol. 2017;41(9):1205–11.

    Article  Google Scholar 

  17. Anderson CM, Zhang B, Miller M, Butko E, Wu X, Laver T, Kernag C, Kim J, Luo Y, Lamparski H, et al. Fully automated RNAscope in situ hybridization assays for formalin-fixed paraffin-embedded cells and tissues. J Cell Biochem. 2016;117(10):2201–8.

    Article  CAS  Google Scholar 

  18. Bakheet AMH, Zhao C, Chen JN, Zhang JY, Huang JT, Du Y, Gong LP, Bi YH, Shao CK. Improving pathological early diagnosis and differential biomarker value for hepatocellular carcinoma via RNAscope technology. Hepatol Int. 2020;14(1):96–104.

    Article  Google Scholar 

  19. Camp RL, Dolled-Filhart M, Rimm DL. X-tile: a new bio-informatics tool for biomarker assessment and outcome-based cut-point optimization. Clin Cancer Res. 2004;10(21):7252–9.

    Article  CAS  Google Scholar 

  20. Molijn A, Jenkins D, Chen W, Zhang X, Pirog E, Enqi W, Liu B, Schmidt J, Cui J, Qiao Y, et al. The complex relationship between human papillomavirus and cervical adenocarcinoma. Int J Cancer. 2016;138(2):409–16.

    Article  CAS  Google Scholar 

  21. Wada T, Ohishi Y, Kaku T, Aman M, Imamura H, Yasutake N, Sonoda K, Kato K, Oda Y. Endocervical adenocarcinoma with morphologic features of both usual and gastric types: clinicopathologic and immunohistochemical analyses and high-risk HPV detection by in situ hybridization. Am J Surg Pathol. 2017;41(5):696–705.

    Article  Google Scholar 

  22. Masoudi H, Van Niekerk DJ, Gilks CB, Cheang M, Bilek K, Fischer U, Ehlen T, Miller D, Horn LC. Loss of p16 INK4 expression in invasive squamous cell carcinoma of the uterine cervix is an adverse prognostic marker. Histopathology. 2006;49(5):542–5.

    Article  CAS  Google Scholar 

  23. Nuovo GJ, Plaia TW, Belinsky SA, Baylin SB, Herman JG. In situ detection of the hypermethylation-induced inactivation of the p16 gene as an early event in oncogenesis. Proc Natl Acad Sci U S A. 1999;96(22):12754–9.

    Article  CAS  Google Scholar 

  24. Poetsch M, Hemmerich M, Kakies C, Kleist B, Wolf E, vom Dorp F, Hakenberg OW, Protzel C. Alterations in the tumor suppressor gene p16(INK4A) are associated with aggressive behavior of penile carcinomas. Virchows Arch. 2011;458(2):221–9.

    Article  CAS  Google Scholar 

  25. Beyer-Finkler E, Girardi F, Sillem M, Pfister H. Human papillomavirus DNA in genital cancers, metastases, and lymph nodes. Intervirology. 1995;38(3–4):173–80.

    Article  CAS  Google Scholar 

  26. Liang WS, Aldrich J, Nasser S, Kurdoglu A, Phillips L, Reiman R, McDonald J, Izatt T, Christoforides A, Baker A, et al. Simultaneous characterization of somatic events and HPV-18 integration in a metastatic cervical carcinoma patient using DNA and RNA sequencing. Int J Gynecol Cancer. 2014;24(2):329–38.

    Article  Google Scholar 

  27. Ronnett BM. Endocervical adenocarcinoma: selected diagnostic challenges. Mod Pathol. 2016;29(Suppl 1):S12-28.

    Article  CAS  Google Scholar 

  28. Baalbergen A, Ewing-Graham PC, Hop WC, Struijk P, Helmerhorst TJ. Prognostic factors in adenocarcinoma of the uterine cervix. Gynecol Oncol. 2004;92(1):262–7.

    Article  CAS  Google Scholar 

  29. Pecorelli S. Revised FIGO staging for carcinoma of the vulva, cervix, and endometrium. Int J Gynaecol Obstet. 2009;105(2):103–4.

    Article  Google Scholar 

Download references


The authors gratefully acknowledge the assistance from professor Jing-Ping Yun from Sun Yat-sen University Cancer Center.


This work was supported by the National Natural Science Foundation of China (No. 82072853). The Natural Science Foundation of Guangdong province (No. 2021A1515010688).

Author information

Authors and Affiliations



Conception and design: L-LL. Performing experiments: ML, YL and XY. Drafting of the article: R-ZL and S-LC. Acquisition and interpretation of data, review, editing and approval of the manuscript: all authors. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Li-Li Liu.

Ethics declarations

Ethics approval and consent to participate

The Sun Yat-sen University Cancer Center Ethics Committee approved this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1

. Representative images for clinicopathological features in paraffin-embedded endocervical adenocarcinoma (ECA) samples. Figure S2. The best cutoff values for all variables as determined by using X-tile in all cases of endocervical adenocarcinoma (ECA). Figure S3. Endocervical adenocarcinoma (ECA): proportions of combinations of positive human papillomavirus (HPV) tests and p16 immunohistochemistry (IHC). Figure S4. The positive rates of human papillomavirus (HPV) subtypes, mRNA, protein, and p16 protein in the subgroups of endocervical adenocarcinoma (ECA) cases.

Additional file 2: Table S1.

The positive proportion of factors examined in ECA samples. Table S2. Correlation of P16 IHC, HPV DNA, HPV genotype and HPV RNAscope. Table S3. HPV genotypes in different ECA histological types.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Luo, RZ., Chen, SL., Li, M. et al. HPV E6/E7 mRNA in situ hybridization in endocervical adenocarcinoma: implications for prognosis and diagnosis. Cancer Cell Int 21, 643 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Endocervical adenocarcinoma
  • Prognosis
  • HPV E6/E7 mRNA in situ hybridization
  • Lymphovascular invasion
  • Lymph node involvement