1. Introduction
Almost all cervical intraepithelial neoplasia (CIN) and cervical cancers are caused by infection with high-risk human papilloma virus (hrHPV) [1]. Most HPV infections are transient and cleared within a year [2]. From initial HPV infection, it takes between 10–30 years to develop invasive cancer [3,4,5]. This slow progression, through several levels of cellular changes, makes cervical cancer a disease that could be marginalized with optimized screening. CIN is considered a premalignant disease, even though only 1% of low-grade dysplasia (CIN1) [6] and 30% of high-grade dysplasia (CIN2-3) will progress to cancer [4]. Approximately, 20–40% of CIN2-3 will regress [6,7,8,9]. Despite this knowledge, all women with CIN2-3 are encouraged to undergo surgical treatment, mainly by removing a cone-shaped piece from the cervix, with risk of complications such as bleeding, infection, perforation, and cervical stenosis leading to dysmenorrhea, risk of occult endometrial cancer, and difficulties in follow-up regarding adequate cytology by positive HPV test [10,11]. The most severe late complication, however, is the increased risk of preterm birth, potentially leading to significantly higher neonatal morbidity and mortality, in addition to a higher number of adverse obstetric sequelae [12,13,14]. Consequently, a significant number of women, who would never be diagnosed with cervical cancer, end up being treated as high-risk CIN leading to psychological distress and potential overtreatment.
Several countries are currently in the process of substituting primary cytology screening with a primary hrHPV-test [15]. Extensive scientific studies show that HPV-based screening yields 23–27% higher sensitivity for detecting high-grade precancerous lesions compared to cytology, consequently preventing more cervical cancer [16,17]. However, more sensitive screening methods also result in higher detection rates of CIN2-3 lesions with a low risk of cancer development [18]. In addition, changed sexual behavior during the last decades, with younger age at first sexual intercourse and higher number of sexual partners, leads to higher exposure to sexually transmitted infections including HPV [19,20,21]. Furthermore, a steadily increasing HPV vaccinated population reaching screening age–presumably with less aggressive HPV infections and lower disease prevalence, is already causing challenges for clinicians in lack of adequate algorithms. Facing these aspects, robust prognostic biomarkers are urgently needed both for better risk stratification and more personalized screening approaches.
Due to poor reproducibility of CIN2, the College of American Pathologists and American Society for Colposcopy and Cervical Pathology strongly recommend implementation of Lower Anogenital Squamous Terminology (LAST), replacing CIN1, 2, and 3 with the categorization of low-grade and high-grade squamous intraepithelial neoplasia (LSIL and HSIL) [22,23]. Dealing with only two categories, the observational procedure will be more clinically challenging, as there will only be one step between low-grade CIN and cancer. Thus, prognostic biomarkers and adequate algorithms will be beneficial for the gynecologists, currently having frequent follow-ups in women with harmless HPV infections and low risk CIN lesions, as well as for the pathologists, currently experiencing a large proportion of normal/low grade biopsies or CIN2-3 with high probability of regression. A better risk stratification will also impact the society, reducing unnecessary health care spending both for follow-up and treatment of low-risk lesions, and reducing the burden of potential complications by over-treatment.
Despite steadily improving screening strategies, the prevalence of cervical cancer has remained stable over the last decades (Cancer Research UK,
2. Materials and Methods
2.1. Patient Cohorts and Biospecimen Collection
The CIN population consisting of 657 women, was included prospectively at the gynecologic outpatient clinic, Stavanger University Hospital (SUS), between 2007–2009 (cohort 1) and 2015–18 (cohort 2) as previously described [7,26,27]. The study was approved by the Regional ethical committee (REK number 2016/805, 2019/264, 2012/1292 and 2020/10399), and written informed consent was obtained by inclusion. Stavanger University Hospital is the only referral hospital for CIN treatment for a population of 370,000. The patients were all referred for further diagnostic evaluation of atypical cervical cytology and/or persistent HPV, according to the former national guidelines from the Norwegian cervical cancer screening program. In total, 333 women diagnosed with CIN2-3 in cervical biopsies were treated by cone excision, median 101 (23–1212) days after the diagnostic biopsy. In total, 63/333 (19%) of the women with proven CIN2-3 in punch biopsies, were diagnosed with only CIN1 or normal tissue in the cone biopsy and defined as regression cases. Two expert pathologists in the field of gynecology, blinded for each other’s diagnoses, evaluated all hematoxylin and eosin (HE) stained slides supported by Ki-67 and p16 immunohistochemical (IHC) staining. Forty-nine formalin fixed paraffin embedded (FFPE) biopsies evaluated preoperatively to represent CIN3 were used in the current study. Persistent CIN3 was defined by diagnosis of CIN3 both in the diagnostic and in the cone biopsy. Following these criteria, 21 were regression cases and 28 persistent CIN3. All patients were extensively followed up after cone excision. Median follow-up time after cone excision was 1886 days (range 119–5173 days) (Table S1). No patients with confirmed regression in the cone excision specimen experienced recurrence with high-grade CIN (CIN2/3) during follow-up.
2.2. RNA/DNA Extraction and p16 Immunohistochemistry
FFPE tissue was used for isolation of RNA and DNA using the miRNeasy FFPE kit (Qiagen) and the E.Z.N.A® Tissue DNA Kit (Omega Bio-tek Inc., Norcross, GA, USA) for cohort 1 and “Recover all total nucleic acid isolation” kit (Thermo Fisher Scientific, Waltham, MA, USA) for cohort 2. For all three methods, isolation was performed following the manufacturer’s instructions. RNA was isolated from 5 µm thick sections comprising the most severe dysplastic area of the epithelium and the adjacent stroma, marked by an expert pathologist, supported by p16 and Ki67 staining. To ensure continuous presence of CIN3, areas adjacent to the sections used for DNA/RNA isolation were HE-stained and evaluated by the expert pathologist.
2.3. Functional RNA Quantification and RNA Reverse Transcription
To quantify the exact amplifiable RNA concentrations from the FFPE samples, a one-step RT-qPCR procedure was applied (LightCycler® 480 System, Roche Diagnostics, Rotkreuz, Switzerland) to measure the RNA concentration with TaqMan Fast Advanced Master mix together with a TaqMan probe specific for the housekeeping gene GUSB (both ThermoFisher Scientific, Waltham, MA USA). To generate target gene standard curves, fourfold dilution series (range of 0.05–50 ng/µL) of a commercially available standard (HL-60, 100 ng/µL total RNA) was used. The RNA concentration for each lesion was calculated, by comparing the mean Ct of triplicates measured for test samples to the Ct measured for the different HL-60 dilutions of the standard curve. The Vilo Mastermix cDNA synthesis kit (Thermo Fisher Scientific) was used for transcription of 10 ng total RNA, as calculated by functional RNA quantification. p16 IHC scoring was performed by using a scale ranging from 0.0 (negative) to 1.0 (strong positive) depending on the proportion of positivity in the epithelium (exclusive basal cells). P16 scoring was available for cohort 1 (n = 32). The scoring was dichotomized; staining index 0.0–0.9 was defined as weak to moderate, and 1.0 as strong positive (strong expression through all epithelial layers). The extent and degree of immunopositivity was scored blinded and independently by two pathologists. In case of discrepancy >10%, which occurred in 8% of all cases, a consensus recount was performed.
2.4. HPV Testing
Biopsies from cohort 1 were genotyped by using Linear Array (Roche Molecular Systems, Roche Diagnostics, Mannheim, Germany), while biopsies from cohort 2 were genotyped using the InnoLipa HPV detection system (Fujirebio Europe N.V., Gent, Belgium). Both are line probe assays based on reverse hybridization principle for detection of both high- and low-risk, in addition to possible high-risk HPV types.
2.5. Next Generation Sequencing (NGS)
The Oncomine™ Immune Response Research Assay targets and quantifies the expression levels of a panel of 398 immune response related genes, including 10 housekeeping genes, using Next Generation Sequencing (NGS), Ion Torrent. Automated library preparation of 2 × 32 samples, each containing 10 ng/µL cDNA, were performed in batches of four libraries of eight samples, according to the Ion Chef system protocol. The library concentrations were measured by the Ion Library TQMN quantitation kit. The Ion OneTouch™ 2 System was used to prepare the enriched, template-positive Ion PI™ Ion Sphere™ Particles (ISPs). The samples were diluted to 100 pM prior to PCR on the Ion OneTouch™ 2 Instrument, using the Ion PI™ Hi-Q ™ OT2 200 kit. A quality control of template-positive ISPs was performed on a Qubit™ 2.0 Fluorometer by using the Ion Sphere™ Quality Control Kit. Enrichment of ISPs was performed by using the Ion OneTouch™ ES instrument. Target sequencing was performed on an Ion Proton instrument by using the Ion PI™ Hi-Q ™ Sequencing 200 chemistry and an Ion PI ™ chip. Subsequently, the sequencing results were downloaded to the Affymetrix Transcriptome Analysis Console (TAC) for further analyses. Mean housekeeping scaled log2 count data from the 398 genes in the Oncomine Immune Response panel, were obtained from the Torrent Suite™ Software. All kits, reagents and software were provided by ThermoFisher Scientific (Waltham, MA, USA) if not otherwise specified.
2.6. Gene Expression Analyses
Differentially expressed genes were called by the Feature Subset Selection (FSS) method within the JExpress software (
2.7. Creating a Six-Gene Signature
FSS analyses were performed to identify differentially expressed genes (DEGs) between CIN3 regression (n = 21) and persistent CIN3 (n = 28) lesions. Six genes matched the criteria of p < 0.05 and fold change <−1.75 or >1.75, consisting of five upregulated genes and one downregulated gene in CIN3 regression lesions. A signature score for each patient was created by subtracting the expression value of the downregulated gene from the sum of expression values of the five upregulated genes. A high signature score should thus associate with CIN3 regression.
2.8. Cancer Cohort
2.8.1. Patient Characteristics and Biospecimen Collection
The cancer cohort study was approved by the Regional ethical committee (REK-number 2018/591 and 2014/1907) with written informed consent from all patients. All patients were diagnosed and treated at the Haukeland University Hospital (HUS). HUS is a referral hospital for patients in Hordaland County in Western Norway, representing approximately 10% of the Norwegian population with similar patterns of incidence and prognosis as from whole of Norway (Cancer Registry of Norway,
2.8.2. Gene Expression Profiling
mRNA expression profiles were generated by the L1000 approach [32] for 239 patients. The L1000 expression data were obtained through an algorithm that extrapolates the expression of 978 directly measured (landmark) genes by a method involving ligation-mediated amplification and fluorescent labelling to generate a transcriptional profile of 12,328 genes in the full L1000 dataset [32]. Replicate-collapsed z-scores (level-5 data) were applied for subsequent L1000 analysis. Of the six genes within the identified gene signature in the precursor lesions, all genes except CCL3 overlapped with the L1000 data creating a 5-gene signature within the cervical cancer cohort. Following the same procedure as for the CIN cohort, a signature score was created for each cervical cancer patient. Optimal gene-signature cut-off values for dichotomization used in GSEA and Kaplan–Meier survival analyses were identified from ROC curves by applying the Youden index [33]. GSEA analyses were performed as described for the CIN cohort. An immune infiltration score was calculated for each patient with available L1000 data by using R version 3.6.3 (Massachusetts, USA) with the ESTIMATE (Estimation of Stromal and Immune cells in MAlignant tumor tissue using Expression) package version 1.0.13 [34].
2.9. Statistical Analyses
Statistical data analyses were performed within the Software package SPSS Statistics (Statistical Package of Social Science) version 27.0 (IBM, Armonk, NY, USA). All probability values were two-sided and considered statistically significant if <0.05. Multiple testing correction was calculated by the Benjamini–Hochberg method. For categorical variables, correlation between groups was assessed using Pearson χ2 or Fisher’s exact test as appropriate. For continuous variables, the Mann–Whitney U or the Kruskal–Wallis test was applied as appropriate. Spearman correlation was applied for detection of non-parametric relationships between pairs of continuous variables. Patient survival analyses were performed by using the Kaplan–Meier (product-limit) method, and survival differences were calculated by the log-rank test (Mantel–Cox). Receiver operating characteristics (ROC) analyses were utilized on the gene-signatures to compare performance related risk groups. Optimal gene-signature cut-off values for prediction of CIN3 regression and cervical cancer survival were identified from ROC curves by applying the Youden index [33] with regression as outcome in the CIN cohort and disease-free survival as outcome in the cancer cohort.
Table 1Distribution of clinicopathological characteristics for all CIN patients included in this study. The number of cases in each group is given followed by percentage for each row in parenthesis.
| Cone Excision Diagnosis | |||
|---|---|---|---|
| CIN3 Regression | Persistent CIN3 | p-Value | |
| n = 21 | n = 28 | ||
| Last cytology before biopsy | 0.71 a | ||
| AGUS | 0 (0) | 1 (100) | |
| ASC-H | 4 (36) | 7 (64) | |
| ASCUS | 0 (0) | 1 (100) | |
| HSIL | 10 (42) | 14 (58) | |
| LSIL | 6 (60) | 4 (40) | |
| Normal | 1 (50) | 1 (50) | |
| HPV Type in Biopsy | 0.79 a | ||
| HPV 16 | 9 (39) | 14 (61) | |
| HPV18 | 2 (40) | 3 (60) | |
| HPV 31 | 1 (50) | 1 (50) | |
| HPV 33 | 4 (36) | 7 (64) | |
| HPV 35 | 2 (100) | 0 (0) | |
| HPV 39 | 1 (50) | 1 (50) | |
| HPV 52 | 2 (50) | 2 (50) | |
| Age at diagnosis | 0.19 b | ||
| ≤29 | 8 (33) | 16 (67) | |
| >29 | 13 (52) | 12 (48) | |
| Interval between cytology and biopsy | 0.32 b | ||
| ≤41 | 12 (50) | 12 (50) | |
| >41 | 9 (36) | 16 (64) | |
a Pearson’s χ2 test. b Mann–Whitney U-test.
Figure 1Identification of a CIN regression signature. (A) Distribution of differentially expressed genes as defined by the criteria of p < 0.05 and fold change <−1.75 or >1.75. (B) Distribution of log 2 expression levels (scaled by housekeeping genes) of the six signature genes and the signature score in lesions of confirmed CIN3 regression versus persistent CIN3. The Man–Whitney U test was applied when if the distribution of the genes were different in CIN3 Regression versus Persistent CIN3. Abbreviations: CIN: Cervical intraepithelial Neoplasia.
[Figure omitted. See PDF]
3. Results
3.1. A Six-Gene Signature Predicting CIN3 Regression
No statistical differences in cytology prior to biopsy, HPV type, age, interval between last cytology and biopsy as well as biopsy cone interval between patients with CIN3 persistence or regression were detected (Table 1). By FSS, six genes matched the criteria of being significantly differentially expressed between CIN3 persistent versus CIN3 regression lesions (Figure 1A, Table S2). TDO2 (p = 0.005), CCL5 (p < 0.001), CCL3 (p = 0.04), CD38 (p = 0.02), and PRF1 (p = 0.005) were upregulated, whilst LCK (p = 0.008) was downregulated in CIN3 regression lesions (Figure 1B, Table S2).
A signature score was calculated for each patient (Figure 2A, for details: See Methods section). Although none of the genes remained significantly differentially expressed after multiple testing (FDR < 0.05), they displayed strong predictive power when considered together as a signature. With an optimal cut-off of 73.1408, 22 patients were classified with low and 27 with high signature score yielding an area under the ROC curve (AUC) of 0.85, a sensitivity of 91% and a specificity of 74% for predicting CIN3 regression (Figure 2B). A strong negative correlation between the six-gene regression signature and p16 protein expression was detected (Figure 2C).
3.2. Persistent CIN3 Associates to Proliferation
In GSEA analyses, within the Hallmark gene sets, “E2F targets genes” and “G2M checkpoint genes” were significantly enriched in persistent CIN3 (Table S3A). Within the C2 curated gene sets, nine out of the top twenty gene sets were related to aggressive cancer or cancer proliferation (Table S3B). Gene sets related to microtubules were significantly enriched in persistent CIN3 within the GO collection (Table S3C).
Figure 2(A) Heatmap illustrating signature score for each gene (column) within each individual patient (row) with risk group defined by red or blue color as indicated. Increasing darker blue color represents increasing gene expression level. In the last column, increasing darker blue color illustrates increasing signature score for each patient sorted by highest to lowest signature score. (B) Receiver operating (ROC) curve reflecting the sensitivity, specificity, and area under the curve (AUC) for the gene signature to predict risk group within the 49 precursor lesions. (C) Immunohistochemistry p16 Staining Index compared to the signature score.
[Figure omitted. See PDF]
3.3. High Regression Signature Score Associates to Favorable Survival and Less Aggressive Features within Cervical Cancer Patients
Next, we sought to characterize the prognostic value of the identified six-gene signature within a cervical cancer cohort with available transcriptional, clinicopathological and follow-up data on 239 patients (Table 2). Within the cancer cohort validation gene set, expression data for CCL3, was not available. Consequently, a five-gene signature was created for each cervical cancer patient (For details: See Methods). The five-gene signature strongly correlated to the six-gene signature (Supplementary Figure S1A, Pearson correlation coefficient of 0.98, p < 0.001), with an AUC of 0.86 (p < 0.01) and the same sensitivity and specificity to predict CIN3 regression as the six-gene signature (Supplementary Figure S1B). We thus concluded that further analyses on this five-gene signature within the cancer cohort was justified.
High signature score, which associated to CIN3 regression in the CIN cohort, correlated to favorable survival (Figure 3A, p = 0.007), MRI derived maximum tumor diameter of less than 2 mm (Figure 3B, p = 0.01), depth of invasion 7 mm or less (Figure 3C, p = 0.03), and no vascular space invasion (Figure 3D, p = 0.04) in the cancer cohort.
Table 2Distribution of available clinicopathological characteristics for the cervical cancer patients according to regression signature score. The number of cases in each group is given followed by percentage for each row in parenthesis.
| Variables (n) a | Regression Signature Score | p-Value b | |
|---|---|---|---|
| Low Score (n = 91) | High Score (n = 148) | ||
| Median age (n = 239) | 0.98 | ||
| <44 years | 42 (38) | 68 (62) | |
| ≥44 years | 49 (38) | 80 (62) | |
| BMI (n = 236) | 0.17 | ||
| <25 | 42 (34) | 83 (66) | |
| ≥25 | 47 (42) | 64 (58) | |
| FIGO-09 stage (n = 239) | 0.27 | ||
| I-IB1 | 40 (34) | 76 (66) | |
| IB2-IV | 51 (42) | 72 (58) | |
| Histologic type (n = 239) | 0.15 | ||
| Squamous cell carcinoma | 59 (35) | 112 (65) | |
| Adenocarcinoma | 24 (50) | 24 (50) | |
| Other histologic type | 8 (40) | 12 (60) | |
| Histologic grade (n = 225) | 0.06 | ||
| Grade 1/2 | 58 (34) | 111 (66) | |
| Grade 3 | 27 (48) | 29 (52) | |
Statistically significant p-values (p < 0.05) are in bold. a number of cases with available data for each variable. b Pearson’s χ2 test.
3.4. High Regression Signature Score Is Significantly Associated with Immune Activation in Cervical Carcinomas
To characterize activated pathways according to signature score, tumors with high versus low signature score were compared by GSEA. Within the Hallmarks, C2 and C5 gene set collections, tumors with high signature score were enriched for gene sets associated with immune activation and response to inflammation (Table S4A–C, respectively). Correspondingly, by histopathological inspection, tumors with high signature score associated with elevated inflammatory reaction (p < 0.001, Figure 4A) and increased expression of the major histocompatibility complex (MHC) Class II receptor, HLA-DQB1 [30] (p = 0.004, Figure 4B). To further investigate the association between signature score and immune activation, we calculated an Immune Score based on the ESTIMATE algorithm (For details: See Methods). Increasing signature score was significantly associated with increasing Immune Score (Shearman correlation = 0.69, p < 0.001, Figure 4C).
Figure 3(A) Disease specific survival for cervical cancer patients relative to signature score represented by Kaplan Meier curve with probability values for Mantel Cox log-rank test that compares categories. The number of patients and events are given within parentheses (patients/events). Signature score in relation to MRI derived maximal tumors diameter (B) histopathological depth of invasion (C) and, vascular space invasion (D) in tumors with available data.
[Figure omitted. See PDF]
4. Discussion
Despite solid evidence that around 50% CIN2 [35] and 30% of CIN3 lesions [6] will regress spontaneously over time, no systematic test for identifying these cases exists in current clinical practice. Consequently, women with low risk of cancer are offered the same follow-up and treatment regime as those with high risk. Although punch biopsies are included in the diagnostic work-up, so far, its potential to reveal features predicting risk-groups by transcriptomic profiling have not been exploited. To investigate this potential, we extracted mRNA from 21 CIN3 lesions with confirmed regression and 28 CIN3 lesions with confirmed persistence. By differential gene expression analyses, we identified a six-gene signature that can predict CIN3 regression with high accuracy. To our knowledge, this is the first study proposing a transcriptomic-based prognostic test, that may serve as a tool for stratifying CIN3 patients into treatment groups without the need for additional sampling.
Figure 4Signature score in relation to (A) inflammatory reaction (B) HLA-DQB1 protein expression and (C) immune score (ESTIMATE). Abbreviations: ESTIMATE: Estimation of Stromal and Immune cells in MAlignant tumor tissue using Expression.
[Figure omitted. See PDF]
The biological mechanisms determining whether the CIN3 lesions will persist or regress may prove relevant for understanding carcinogenic processes — also in the cervical cancer setting. Genes that are highly expressed in persistent CIN3 could be early drivers of cervical cancer and could serve as prognostic markers and/or therapeutic targets in cervical cancer. Within the six-gene signature, LCK was the only gene upregulated in persistent CIN3. Lymphocyte-specific protein tyrosine kinase (LCK) is a member of the Src family of tyrosine kinases predominantly expressed on T cells to regulate T cell development and homeostasis [36]. However, LCK has also been found expressed in several solid tumors including brain [37], breast [38], prostate [39], and colorectal [40] cancer. In breast cancer, LCK promotes angiogenesis and tumor progression [41], indicating an oncogenic role in line with our findings. So far, LCK expression has not been described in cervical precursor nor cancer lesions. Thus, a thorough characterization of its potential as prognostic and/or predictive marker in the precancer and cancer setting is warranted.
Genes that are highly expressed in CIN3 regression as compared to persistent CIN3 could be potential prognostic markers for CIN3 regression and of less aggressive cancer. TDO2, CD38, CCL5, CCL3, and PFR1 were upregulated in CIN3 regression lesions. Tryptophan 2,3-dioxygenase (TDO) is a heme enzyme encoded by the TDO2 gene that catalyzes the first and rate-limiting step of tryptophan degradation along the kynurenine pathway [42]. Tryptophan depletion has been linked to aggressive cervical cancer [43] and impaired immune responses by restraining T cell proliferation and restricting tumors immune cell infiltration [44]. Within the TCGA cohort, TDO2 expression associates with poor survival (p = 0.003,
CD38 is a type II transmembrane glycoprotein found on the surface of many immune cells including macrophages, CD4+, CD8+, B lymphocytes, and natural killer cells and is known to activate B and T cells. In addition, CD38 has enzymatic activity as a NADase [47] and regulates diverse pathways such as signal transduction, cell adhesion, cyclic ADP-ribose synthesis, and cell differentiation and activation. In cervical cancer, CD38 has been linked to increased proliferation, inhibition of apoptosis [48] and PI3K/AKT/mTOR signaling activation [49]. Yet, we found that CD38 associates with CIN3 regression. In line with our findings, CD38 associates with favorable prognosis in the TCGA cervical cancer cohort (p = 0.008,
CCL3 and CCL5 are chemokines involved in recruitment of a variety of leukocytes to inflammatory sites, including CD8+ T cells, macrophages, eosinophils, and basophils. These chemokines have both anti- and pro-cancer properties [52]. In a study investigating 23 cervical carcinomas, CCL5 levels were found to increase with advancing disease state [53]. However, CCL5 mRNA levels associate with good prognosis in the cervical cancer TCGA dataset (p = 0.002,
Perforin-1 is encoded by the PRF1 gene and is a pore forming cytolytic protein found in granules of cytotoxic T lymphocytes and natural killer cells. Perforin is known as a key enabler of granzyme-induced apoptosis of target cancer cells [54]. In cervical cancer, down-regulation of perforin has been linked to immune escape [55]. In CIN lesions, decreased perforin granule release in CD8+ T cells resulted in more severe lesions [56]. These studies are all in line with our findings that PRF1 is upregulated in CIN3 regression. Whether perforin expression could be used solely as an IHC detectable marker for CIN2-3 regression warrants further investigation.
GSEA analyses revealed that the Hallmark gene sets “E2F targets” and “G2M checkpoint genes” were significantly enriched in persistent CIN3. While the “G2M checkpoint genes” defines genes associated with G2/M progression, the “E2F targets” include gene linked to G1/S transition. In a non-proliferative state, E2F is bound to the tumor suppressor Rb which inhibits E2F signaling. When the HPV oncoprotein E7 binds to Rb, E2F becomes activated which in turn enables G1/S transition [57]. Within the C2 curated gene sets, 10 out of the top 20 genes sets were related to proliferation of cancer. Microtubules are important for proliferation of cells, and within the GO collection, gene sets related to microtubules were enriched in persistent CIN3. Hence, the GSEA results are unambiguously indicating that proliferation is a central characteristic of persistent CIN3 and that HPV oncoproteins, such as E7, are important drivers for increased proliferation. When E2F is released from Rb, protein levels of the tumor suppressor p16 increases due to a positive feedback loop. This can be detected on IHC and p16 positivity is regarded as a diagnostic marker for HPV infection [58,59]. Interestingly, we found that low regression signature score (associating with persistent CIN3) linked to p16 positivity. This could thus indicate a stronger HPV signaling in the persistent CIN3 cells.
When we applied the regression signature on the cancer cohort, we found that a high score associated significantly to good survival, small tumors, and no vascular space invasion. Vascular space invasion is a prognostic marker of poor survival [31] and lymph node metastases [60] in cervical cancer. Large tumors are also strongly linked to poor prognosis. These results confirm the theory that the presumably less aggressive CIN3 regression lesions share common features with less aggressive cervical carcinomas, whereas persistent CIN3 lesions resembles cancer lesions associated with poor survival.
By GSEA, we discovered that the tumors with high regression signature score and good prognosis were significantly linked to immune activation, immune cell infiltration and increased ESTIMATE immune score. Moreover, these tumors also expressed higher levels of the major histocompatibility complex Class II receptor HLA-DQB1 known as a marker for immune activation, infiltrating immune cells and good prognosis [30]. In order to progress, cancers and precancerous lesions need to bypass immune control. Reduced HLA Class I and II expression in HPV infected tissue may be linked to the HPV oncoprotein E5, which prevents presentation of MHC molecules instigating an immune evading mechanism [25]. In the less aggressive cancer lesions, HLA molecules may still actively express antigens allowing CD8+ T cells to identify the pathological cells and thereby restraining the tumor from uncontrolled growth. This may indicate that tumors that resemble the CIN3 regression lesions are, in addition to being less aggressive, also associated with immune activation and immune cell infiltration. When regression from CIN3 occurs, the immune cells overcome the immune evading and proliferative mechanisms initiated by the HPV infection. Our results indicate that these mechanisms are also in play in less aggressive cervical carcinomas decelerating further invasion or metastatic spread.
A limitation of the study is the relatively short observational period between biopsy and cone excision as regression can take from months to years [6,61,62]. However, a longer observational period in the CIN2-3 patients would have been unethical, considering the risk of progression to cancer. Furthermore, the sample size was moderate, which limits the separate interpretation of prognosis in some of the subgroups (i.e., smoking, CIN grade 2 or 3, CIN or cervical cancer in 1st or 2nd degree relatives). An increased sample size would increase the statistical power to identify differentially expressed genes. We were not able to identify differentially expressed genes that remained significant after multiple testing. Still, when applied as a signature, the six genes displayed high predictive power. In addition, Oncomine only includes a subset of 398 immune response related genes, thus other transcriptomic differences between CIN3 regression and persistent CIN3 will go undetected. Another challenge in the interpretation of gene expression data is the lack of spatial and cellular differentiation. To capture the complexity of immune regulation, differentiation between the epithelial, immune related and stromal cells, and the structural composition of the lesions and surrounding tissue, need to be characterized. In the cancer cohort, the application of the gene signature was suboptimal as one gene within the six-gene signature was not included in the transcriptomic cancer cohort data. Still, as a group, the five remaining genes proved useful in predicting disease-specific survival.
Strengths of the current study include the prospective sample collection and close follow-up. Stavanger University Hospital is the only hospital in the region recruiting patients to this study, ensuring a uniform inclusion strategy, and minimizing selection biases. The cervical dysplasia cohort consisted of healthy women with a first-time onset histologically (p16/Ki67 supported) confirmed CIN3, evaluated by two independent pathologists. All patients were followed up after cone excision, and none of the women diagnosed with regression experienced recurrence with high-grade CIN (CIN2/3) during follow-up. Furthermore, although the precancer cohort size is small and the transcriptomic data are confined to 398 immune response related genes, the unprecedented nature of this study and the acquisition of high-value lesions with confirmed regression, makes it highly relevant as a starting point for development of prognostic tools for CIN3 regression. The cervical cancer samples applied for transcriptomic data analysis were collected prospectively in a population-based setting. Extensive clinicopathological characterization and follow-up enabled validation and exploration of results obtained in a precancerous setting on a cancer cohort. By this approach, we have identified previously undetected molecular similarities as well as differences between CIN3 regression, persistent CIN3 and cervical cancer. We consider this study as an important first step towards the development of a clinically applicable pre-operative test for identifying CIN3 regression and hence avoid overtreatment. However, these findings should be validated in larger cohorts.
The development of CIN and cervical cancer is driven by a complex interaction between host and viral factors. Clinicopathological factors including HPV genotype, tobacco smoking, hormonal contraceptives, socioeconomic aspects, reproductive, and sexual factors have been reported as significant for CIN development [63,64,65,66]. Yet, only a few studies have investigated potential prognostic markers for CIN regression, and predictive power have been identified for HPV genotype [67], microbiota [68], age [69], parity, and smoking status [70]. CIN2 is known to have a higher probability of regression, and whether conservative treatment should be recommended for these patients is controversial. In Denmark, national guidelines have recently been changed allowing more than half of women diagnosed with CIN2 conservative treatment. Interestingly, half of these women experienced lesion regression [71], underlining the fact that conservative treatment could be an adequate approach for many CIN2 and potentially CIN3 patients. We present a regression signature that could contribute to risk stratification. However, while our regression signature was identified in CIN3 patients, all the above-mentioned studies were performed on CIN2 lesions. Hence, this study generates novel insights into possible models predicting CIN3 regression. Furthermore, the presented results highlight mechanisms driving persistent CIN3 lesions and provide clues to develop new treatment strategies for patients with CIN2-3.
5. Conclusions
CIN2-3 is characterized by a high spontaneous regression rate. Yet, currently, histological evaluation is inadequate to distinguish CIN2-3 lesions that will regress from those that will persist or progress [72]. To our knowledge, this is the first study to date, investigating broad immunological differences between CIN3 regression and persistence. By performing the Oncomine assay, we propose a six-gene signature from 398 immune related genes, predicting regression with a sensitivity of 91% (AUC = 0.85). Furthermore, LCK, Perforin, and CD38 were identified as potential prognostic IHC biomarkers and should be evaluated in larger patient cohorts of both precancer and cancer patients. By unravelling the immune profile of CIN3 regression versus persistence and by comparing this to a large cancer cohort, we have contributed to the etiological understanding of cervical carcinogenesis and have unraveled potential vulnerabilities that may be exploited in future conservative treatment strategies for CIN2-3. Some immune modulating treatments including local use of probiotics, interferon, and biological agents such as Imiquimod and 5-Fluorouracil, have already proven clinical benefit for CIN2-3 regression [73,74,75]. This focus on enabling CIN regression as a possible endpoint instead of quick surgical intervention is a welcoming approach to avoid overtreatment and underpins the growing importance of identifying robust markers for better patient stratification. The potential for development of prognostic markers as well as therapeutics utilizing these differences needs to be further explored in larger and independent high-grade CIN and cervical cancer patient cohorts.
Conceptualization, M.K.H., A.C.M., B.E., C.K., and I.T.Ø.; Data curation, M.K.H., A.C.M., B.E., S.A., E.A.H., K.W., I.S.H., B.I.B., E.A.M.J., E.G., C.K., and I.T.Ø.; Formal analysis, M.K.H., A.C.M., S.A., E.A.H., and D.F., Kristine E Fasmer and I.T.Ø.; Funding acquisition, M.K.H., B.E., C.K., and I.T.Ø.; Investigation, M.K.H., A.C.M., and S.A., Astri Frafjord and I.T.Ø.; Methodology, M.K.H., A.C.M., S.A., A.F., D.F., K.E.F., and I.T.Ø.; Project administration, M.K.H., A.C.M., B.E., C.K., and I.T.Ø.; Resources, M.K.H., A.C.M., B.E., K.W., I.S.H., B.I.B., E.A.M.J., E.G., C.K., and I.T.Ø.; Software, M.K.H., D.F., K.E.F., and I.T.Ø.; Supervision, I.S.H., E.A.M.J., E.G., C.K., and I.T.Ø.; Validation, M.K.H., A.C.M., B.E., S.A., A.F., C.K. and I.T.Ø.; Visualization, M.K.H., B.E., and C.K.; Writing–original draft, M.K.H., A.C.M., and I.T.Ø.; Writing—review and editing, M.K.H., A.C.M., B.E., S.A., A.F., E.A.H., D.F., K.E.F., K.W., I.S.H., B.I.B., E.A.M.J., E.G., C.K., and I.T.Ø. All authors have read and agreed to the published version of the manuscript.
This research was funded by Folke Hermansen Fond (year 2017-2019), Helse Vest (Grant Number 27804 and 911634), Universitetet i Bergen, Norges Forskningsråd (Grant Number 273280), and Kreftforeningen (Grant Number 190202).
The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee, REC Heath West: 2012/1292, 2016/805, 2019/264, 2020/10399, 2014/1907 and 2018/591.
Informed consent was obtained from all subjects involved in the study.
All data of this study are available within
The authors would like to acknowledge Olivera Bozickovic, Kadri Madissoo and Bendik Nordanger for excellent technical assistance.
The authors declare no conflict of interest.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.
References
1. Schiffman, M.; Castle, P.E.; Jeronimo, J.; Rodriguez, A.C.; Wacholder, S. Human papillomavirus and cervical cancer. Lancet; 2007; 370, pp. 890-907. [DOI: https://dx.doi.org/10.1016/S0140-6736(07)61416-0]
2. Rodríguez, 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. J. Natl. Cancer Inst.; 2008; 100, pp. 513-517. [DOI: https://dx.doi.org/10.1093/jnci/djn044]
3. Khan, M.J.; Castle, P.E.; Lorincz, A.T.; Wacholder, S.; Sherman, M.; Scott, D.R.; Rush, B.B.; Glass, A.G.; Schiffman, M. The elevated 10-year risk of cervical precancer and cancer in women with human papillomavirus (HPV) type 16 or 18 and the possible utility of type-specific HPV testing in clinical practice. J. Natl. Cancer Inst.; 2005; 97, pp. 1072-1079. [DOI: https://dx.doi.org/10.1093/jnci/dji187] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16030305]
4. McCredie, M.R.; Sharples, K.J.; Paul, C.; Baranyai, J.; Medley, G.; Jones, R.W.; Skegg, D.C. Natural history of cervical neoplasia and risk of invasive cancer in women with cervical intraepithelial neoplasia 3: A retrospective cohort study. Lancet Oncol.; 2008; 9, pp. 425-434. [DOI: https://dx.doi.org/10.1016/S1470-2045(08)70103-7]
5. Tjalma, W.A.; Fiander, A.; Reich, O.; Powell, N.; Nowakowski, A.M.; Kirschner, B.; Koiss, R.; O’Leary, J.; Joura, E.A.; Rosenlund, M. et al. Differences in human papillomavirus type distribution in high-grade cervical intraepithelial neoplasia and invasive cervical cancer in Europe. Int. J. Cancer; 2013; 132, pp. 854-867. [DOI: https://dx.doi.org/10.1002/ijc.27713] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22752992]
6. Ostör, A.G. Natural history of cervical intraepithelial neoplasia: A critical review. Int. J. Gynecol. Pathol.; 1993; 12, pp. 186-192. [DOI: https://dx.doi.org/10.1097/00004347-199304000-00018]
7. Munk, A.C.; Gudlaugsson, E.; Ovestad, I.T.; Lovslett, K.; Fiane, B.; Hidle, B.; Kruse, A.J.; Skaland, I.; Janssen, E.A.; Baak, J.P. Interaction of epithelial biomarkers, local immune response and condom use in cervical intraepithelial neoplasia 2-3 regression. Gynecol Oncol.; 2012; 127, pp. 489-494. [DOI: https://dx.doi.org/10.1016/j.ygyno.2012.09.010] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23017821]
8. Ovestad, I.T.; Gudlaugsson, E.; Skaland, I.; Malpica, A.; Munk, A.C.; Janssen, E.A.; Baak, J.P. The impact of epithelial biomarkers, local immune response and human papillomavirus genotype in the regression of cervical intraepithelial neoplasia grades 2-3. J. Clin. Pathol.; 2011; 64, pp. 303-307. [DOI: https://dx.doi.org/10.1136/jcp.2010.083626]
9. Trimble, C.L.; Piantadosi, S.; Gravitt, P.; Ronnett, B.; Pizer, E.; Elko, A.; Wilgus, B.; Yutzy, W.; Daniel, R.; Shah, K. et al. Spontaneous regression of high-grade cervical dysplasia: Effects of human papillomavirus type and HLA phenotype. Clin. Cancer Res.; 2005; 11, pp. 4717-4723. [DOI: https://dx.doi.org/10.1158/1078-0432.CCR-04-2599]
10. Hagen, B.; Skjeldestad, F.E.; Bratt, H.; Tingulstad, S.; Lie, A.K. CO2 laser conization for cervical intraepithelial neoplasia grade II-III: Complications and efficacy. Acta Obstet. Gynecol. Scand.; 1998; 77, pp. 558-563. [DOI: https://dx.doi.org/10.1034/j.1600-0412.1998.770515.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/9654180]
11. Santesso, N.; Mustafa, R.A.; Wiercioch, W.; Kehar, R.; Gandhi, S.; Chen, Y.; Cheung, A.; Hopkins, J.; Khatib, R.; Ma, B. et al. Systematic reviews and meta-analyses of benefits and harms of cryotherapy, LEEP, and cold knife conization to treat cervical intraepithelial neoplasia. Int. J. Gynaecol. Obstet.; 2016; 132, pp. 266-271. [DOI: https://dx.doi.org/10.1016/j.ijgo.2015.07.026] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/26643302]
12. Albrechtsen, S.; Rasmussen, S.; Thoresen, S.; Irgens, L.M.; Iversen, O.E. Pregnancy outcome in women before and after cervical conisation: Population based cohort study. BMJ; 2008; 337, a1343. [DOI: https://dx.doi.org/10.1136/bmj.a1343]
13. Arbyn, M.; Kyrgiou, M.; Simoens, C.; Raifu, A.O.; Koliopoulos, G.; Martin-Hirsch, P.; Prendiville, W.; Paraskevaidis, E. Perinatal mortality and other severe adverse pregnancy outcomes associated with treatment of cervical intraepithelial neoplasia: Meta-analysis. BMJ; 2008; 337, a1284. [DOI: https://dx.doi.org/10.1136/bmj.a1284] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18801868]
14. 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; 354, i3633. [DOI: https://dx.doi.org/10.1136/bmj.i3633] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27469988]
15. Maver, P.J.; Poljak, M. Primary HPV-based cervical cancer screening in Europe: Implementation status, challenges, and future plans. Clin. Microbiol. Infect.; 2020; 26, pp. 579-583. [DOI: https://dx.doi.org/10.1016/j.cmi.2019.09.006]
16. Arbyn, M.; Ronco, G.; Anttila, A.; Meijer, C.J.; Poljak, M.; Ogilvie, G.; Koliopoulos, G.; Naucler, P.; Sankaranarayanan, R.; Peto, J. Evidence regarding human papillomavirus testing in secondary prevention of cervical cancer. Vaccine; 2012; 30, pp. F88-F99. [DOI: https://dx.doi.org/10.1016/j.vaccine.2012.06.095]
17. Ronco, G.; Dillner, J.; Elfström, K.M.; Tunesi, S.; Snijders, P.J.; Arbyn, M.; Kitchener, H.; Segnan, N.; Gilham, C.; Giorgi-Rossi, P. et al. Efficacy of HPV-based screening for prevention of invasive cervical cancer: Follow-up of four European randomised controlled trials. Lancet; 2014; 383, pp. 524-532. [DOI: https://dx.doi.org/10.1016/S0140-6736(13)62218-7]
18. Orumaa, M.; Leinonen, M.K.; Campbell, S.; Møller, B.; Myklebust, T.; Nygård, M. Recent increase in incidence of cervical precancerous lesions in Norway: Nationwide study from 1992 to 2016. Int. J. Cancer; 2019; 145, pp. 2629-2638. [DOI: https://dx.doi.org/10.1002/ijc.32195]
19. Guleria, S.; Faber, M.T.; Hansen, B.T.; Arnheim-Dahlström, L.; Liaw, K.L.; Munk, C.; Nygård, M.; Kjær, S.K. Self-perceived risk of STIs in a population-based study of Scandinavian women. Sex. Transm Infect.; 2018; 94, pp. 522-527. [DOI: https://dx.doi.org/10.1136/sextrans-2017-053397]
20. Liu, G.; Hariri, S.; Bradley, H.; Gottlieb, S.L.; Leichliter, J.S.; Markowitz, L.E. Trends and patterns of sexual behaviors among adolescents and adults aged 14 to 59 years, United States. Sex. Transm Dis.; 2015; 42, pp. 20-26. [DOI: https://dx.doi.org/10.1097/OLQ.0000000000000231] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25504296]
21. Stigum, H.; Samuelsen, S.O.; Traeen, B. Analysis of first coitus. Arch. Sex. Behav.; 2010; 39, pp. 907-914. [DOI: https://dx.doi.org/10.1007/s10508-009-9494-6] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19330437]
22. Waxman, A.G.; Chelmow, D.; Darragh, T.M.; Lawson, H.; Moscicki, A.B. Revised terminology for cervical histopathology and its implications for management of high-grade squamous intraepithelial lesions of the cervix. Obstet. Gynecol.; 2012; 120, pp. 1465-1471. [DOI: https://dx.doi.org/10.1097/AOG.0b013e31827001d5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23168774]
23. Castle, P.E.; Adcock, R.; Cuzick, J.; Wentzensen, N.; Torrez-Martinez, N.E.; Torres, S.M.; Stoler, M.H.; Ronnett, B.M.; Joste, N.E.; Darragh, T.M. et al. Relationships of p16 Immunohistochemistry and Other Biomarkers with Diagnoses of Cervical Abnormalities: Implications for LAST Terminology. Arch. Pathol. Lab. Med.; 2020; 144, pp. 725-734. [DOI: https://dx.doi.org/10.5858/arpa.2019-0241-OA] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31718233]
24. Cohen, P.A.; Jhingran, A.; Oaknin, A.; Denny, L. Cervical cancer. Lancet; 2019; 393, pp. 169-182. [DOI: https://dx.doi.org/10.1016/S0140-6736(18)32470-X]
25. Shamseddine, A.A.; Burman, B.; Lee, N.Y.; Zamarin, D.; Riaz, N. Tumor Immunity and Immunotherapy for HPV-Related Cancers. Cancer Discov.; 2021; 11, pp. 1896-1912. [DOI: https://dx.doi.org/10.1158/2159-8290.CD-20-1760]
26. Munk, A.C.; Gudlaugsson, E.; Malpica, A.; Fiane, B.; Løvslett, K.I.; Kruse, A.J.; Øvestad, I.T.; Voorhorst, F.; Janssen, E.A.; Baak, J.P. Consistent condom use increases the regression rate of cervical intraepithelial neoplasia 2-3. PLoS ONE; 2012; 7, e45114. [DOI: https://dx.doi.org/10.1371/journal.pone.0045114]
27. Munk, A.C.; Ovestad, I.T.; Gudlaugsson, E.; Løvslett, K.; Fiane, B.; van Diermen-Hidle, B.; Kruse, A.J.; Skaland, I.; Janssen, E.A.; Baak, J.P. Consistent condom use increases spontaneous regression in high-risk non-HPV16 but not in HPV16 CIN2-3 lesions, a prospective population-based cohort study. Infect. Agent Cancer; 2012; 7, 30. [DOI: https://dx.doi.org/10.1186/1750-9378-7-30]
28. Dysvik, B.; Jonassen, I. J-Express: Exploring gene expression data using Java. Bioinformatics; 2001; 17, pp. 369-370. [DOI: https://dx.doi.org/10.1093/bioinformatics/17.4.369]
29. Liberzon, A.; Birger, C.; Thorvaldsdóttir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst.; 2015; 1, pp. 417-425. [DOI: https://dx.doi.org/10.1016/j.cels.2015.12.004]
30. Halle, M.K.; Sodal, M.; Forsse, D.; Engerud, H.; Woie, K.; Lura, N.G.; Wagner-Larsen, K.S.; Trovik, J.; Bertelsen, B.I.; Haldorsen, I.S. et al. A 10-gene prognostic signature points to LIMCH1 and HLA-DQB1 as important players in aggressive cervical cancer disease. Br. J. Cancer; 2021; [DOI: https://dx.doi.org/10.1038/s41416-021-01305-0]
31. Halle, M.K.; Ojesina, A.I.; Engerud, H.; Woie, K.; Tangen, I.L.; Holst, F.; Hoivik, E.; Kusonmano, K.; Haldorsen, I.S.; Vintermyr, O.K. et al. Clinicopathologic and molecular markers in cervical carcinoma: A prospective cohort study. Am. J. Obstet. Gynecol.; 2017; 217, pp. 432.e1-432.e17. [DOI: https://dx.doi.org/10.1016/j.ajog.2017.05.068]
32. Subramanian, A.; Narayan, R.; Corsello, S.M.; Peck, D.D.; Natoli, T.E.; Lu, X.; Gould, J.; Davis, J.F.; Tubelli, A.A.; Asiedu, J.K. et al. A Next Generation Connectivity Map: L1000 Platform and the First 1,000,000 Profiles. Cell; 2017; 171, pp. 1437-1452.e17. [DOI: https://dx.doi.org/10.1016/j.cell.2017.10.049]
33. Youden, W.J. Index for rating diagnostic tests. Cancer; 1950; 3, pp. 32-35. [DOI: https://dx.doi.org/10.1002/1097-0142(1950)3:1<32::AID-CNCR2820030106>3.0.CO;2-3]
34. Yoshihara, K.; Shahmoradgoli, M.; Martínez, E.; Vegesna, R.; Kim, H.; Torres-Garcia, W.; Treviño, V.; Shen, H.; Laird, P.W.; Levine, D.A. et al. Inferring tumour purity and stromal and immune cell admixture from expression data. Nature Commun.; 2013; 4, 2612. [DOI: https://dx.doi.org/10.1038/ncomms3612] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/24113773]
35. Tainio, K.; Athanasiou, A.; Tikkinen, K.A.O.; Aaltonen, R.; Cárdenas, J.; Hernándes,; Glazer-Livson, S.; Jakobsson, M.; Joronen, K.; Kiviharju, M. et al. Clinical course of untreated cervical intraepithelial neoplasia grade 2 under active surveillance: Systematic review and meta-analysis. BMJ; 2018; 360, k499. [DOI: https://dx.doi.org/10.1136/bmj.k499] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/29487049]
36. Bommhardt, U.; Schraven, B.; Simeoni, L. Beyond TCR Signaling: Emerging Functions of Lck in Cancer and Immunotherapy. Int. J. Mol. Sci.; 2019; 20, 3500. [DOI: https://dx.doi.org/10.3390/ijms20143500] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31315298]
37. Kim, R.K.; Yoon, C.H.; Hyun, K.H.; Lee, H.; An, S.; Park, M.J.; Kim, M.J.; Lee, S.J. Role of lymphocyte-specific protein tyrosine kinase (LCK) in the expansion of glioma-initiating cells by fractionated radiation. Biochem. Biophys. Res. Commun.; 2010; 402, pp. 631-636. [DOI: https://dx.doi.org/10.1016/j.bbrc.2010.10.072] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20971076]
38. Elsberger, B.; Fullerton, R.; Zino, S.; Jordan, F.; Mitchell, T.J.; Brunton, V.G.; Mallon, E.A.; Shiels, P.G.; Edwards, J. Breast cancer patients’ clinical outcome measures are associated with Src kinase family member expression. Br. J. Cancer; 2010; 103, pp. 899-909. [DOI: https://dx.doi.org/10.1038/sj.bjc.6605829]
39. Robinson, D.; He, F.; Pretlow, T.; Kung, H.J. A tyrosine kinase profile of prostate carcinoma. Proc. Natl. Acad. Sci. USA; 1996; 93, pp. 5958-5962. [DOI: https://dx.doi.org/10.1073/pnas.93.12.5958]
40. Veillette, A.; Foss, F.M.; Sausville, E.A.; Bolen, J.B.; Rosen, N. Expression of the lck tyrosine kinase gene in human colon carcinoma and other non-lymphoid human tumor cell lines. Oncogene Res.; 1987; 1, pp. 357-374.
41. Chakraborty, G.; Rangaswami, H.; Jain, S.; Kundu, G.C. Hypoxia regulates cross-talk between Syk and Lck leading to breast cancer progression and angiogenesis. J. Biol. Chem.; 2006; 281, pp. 11322-11331. [DOI: https://dx.doi.org/10.1074/jbc.M512546200]
42. Pilotte, L.; Larrieu, P.; Stroobant, V.; Colau, D.; Dolusic, E.; Frédérick, R.; De Plaen, E.; Uyttenhove, C.; Wouters, J.; Masereel, B. et al. Reversal of tumoral immune resistance by inhibition of tryptophan 2,3-dioxygenase. Proc. Natl. Acad. Sci. USA; 2012; 109, pp. 2497-2502. [DOI: https://dx.doi.org/10.1073/pnas.1113873109]
43. Ferns, D.M.; Kema, I.P.; Buist, M.R.; Nijman, H.W.; Kenter, G.G.; Jordanova, E.S. Indoleamine-2,3-dioxygenase (IDO) metabolic activity is detrimental for cervical cancer patient survival. Oncoimmunology; 2015; 4, e981457. [DOI: https://dx.doi.org/10.4161/2162402X.2014.981457] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25949879]
44. Platten, M.; Nollen, E.A.A.; Röhrig, U.F.; Fallarino, F.; Opitz, C.A. Tryptophan metabolism as a common therapeutic target in cancer, neurodegeneration and beyond. Nat. Rev. Drug. Discov.; 2019; 18, pp. 379-401. [DOI: https://dx.doi.org/10.1038/s41573-019-0016-5] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30760888]
45. Venancio, P.A.; Consolaro, M.E.L.; Derchain, S.F.; Boccardo, E.; Villa, L.L.; Maria-Engler, S.S.; Campa, A.; Discacciati, M.G. Indoleamine 2,3-dioxygenase and tryptophan 2,3-dioxygenase expression in HPV infection, SILs, and cervical cancer. Cancer Cytopathol.; 2019; 127, pp. 586-597. [DOI: https://dx.doi.org/10.1002/cncy.22172] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31412167]
46. Zhou, C.; Tuong, Z.K.; Frazer, I.H. Papillomavirus Immune Evasion Strategies Target the Infected Cell and the Local Immune System. Front. Oncol.; 2019; 9, 682. [DOI: https://dx.doi.org/10.3389/fonc.2019.00682] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31428574]
47. Konen, J.M.; Fradette, J.J.; Gibbons, D.L. The Good, the Bad and the Unknown of CD38 in the Metabolic Microenvironment and Immune Cell Functionality of Solid Tumors. Cells; 2019; 9, 52. [DOI: https://dx.doi.org/10.3390/cells9010052]
48. Liao, S.; Xiao, S.; Chen, H.; Zhang, M.; Chen, Z.; Long, Y.; Gao, L.; Zhu, G.; He, J.; Peng, S. et al. CD38 enhances the proliferation and inhibits the apoptosis of cervical cancer cells by affecting the mitochondria functions. Mol. Carcinog.; 2017; 56, pp. 2245-2257. [DOI: https://dx.doi.org/10.1002/mc.22677]
49. Liao, S.; Liang, L.; Yue, C.; He, J.; He, Z.; Jin, X.; Luo, G.; Zhou, Y. CD38 is involved in cell energy metabolism via activating the PI3K/AKT/mTOR signaling pathway in cervical cancer cells. Int. J. Oncol.; 2020; 57, pp. 338-354. [DOI: https://dx.doi.org/10.3892/ijo.2020.5040] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32319590]
50. Lam, J.H.; Ng, H.H.M.; Lim, C.J.; Sim, X.N.; Malavasi, F.; Li, H.; Loh, J.J.H.; Sabai, K.; Kim, J.K.; Ong, C.C.H. et al. Expression of CD38 on Macrophages Predicts Improved Prognosis in Hepatocellular Carcinoma. Front. Immunol.; 2019; 10, 2093. [DOI: https://dx.doi.org/10.3389/fimmu.2019.02093]
51. Ng, H.H.M.; Lee, R.Y.; Goh, S.; Tay, I.S.Y.; Lim, X.; Lee, B.; Chew, V.; Li, H.; Tan, B.; Lim, S. et al. Immunohistochemical scoring of CD38 in the tumor microenvironment predicts responsiveness to anti-PD-1/PD-L1 immunotherapy in hepatocellular carcinoma. J. Immunother. Cancer; 2020; 8, [DOI: https://dx.doi.org/10.1136/jitc-2020-000987] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/32847986]
52. Korbecki, J.; Grochans, S.; Gutowska, I.; Barczak, K.; Baranowska-Bosiacka, I. CC Chemokines in a Tumor: A Review of Pro-Cancer and Anti-Cancer Properties of Receptors CCR5, CCR6, CCR7, CCR8, CCR9, and CCR10 Ligands. Int. J. Mol. Sci.; 2020; 21, 7619. [DOI: https://dx.doi.org/10.3390/ijms21207619] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33076281]
53. Niwa, Y.; Akamatsu, H.; Niwa, H.; Sumi, H.; Ozaki, Y.; Abe, A. Correlation of tissue and plasma RANTES levels with disease course in patients with breast or cervical cancer. Clin. Cancer Res.; 2001; 7, pp. 285-289. [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/11234881]
54. Voskoboinik, I.; Whisstock, J.C.; Trapani, J.A. Perforin and granzymes: Function, dysfunction and human pathology. Nat. Rev. Immunol.; 2015; 15, pp. 388-400. [DOI: https://dx.doi.org/10.1038/nri3839] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/25998963]
55. Sheu, B.C.; Chiou, S.H.; Lin, H.H.; Chow, S.N.; Huang, S.C.; Ho, H.N.; Hsu, S.M. Up-regulation of inhibitory natural killer receptors CD94/NKG2A with suppressed intracellular perforin expression of tumor-infiltrating CD8+ T lymphocytes in human cervical carcinoma. Cancer Res.; 2005; 65, pp. 2921-2929. [DOI: https://dx.doi.org/10.1158/0008-5472.CAN-04-2108] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15805295]
56. Fernandes, A.T.; da Rocha, N.P.; Avvad, E.; Grinsztejn, B.J.; Russomano, F.; Tristão, A.; Quintana Mde, S.; Perez, M.A.; Conceição-Silva, F.; Bonecini-Almeida Mda, G. Balance of apoptotic and anti-apoptotic marker and perforin granule release in squamous intraepithelial lesions. HIV infection leads to a decrease in perforin degranulation. Exp. Mol. Pathol.; 2013; 95, pp. 166-173. [DOI: https://dx.doi.org/10.1016/j.yexmp.2013.06.006]
57. zur Hausen, H. Papillomaviruses causing cancer: Evasion from host-cell control in early events in carcinogenesis. J. Natl. Cancer Inst.; 2000; 92, pp. 690-698. [DOI: https://dx.doi.org/10.1093/jnci/92.9.690] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10793105]
58. Mahajan, A. Practical issues in the application of p16 immunohistochemistry in diagnostic pathology. Hum. Pathol.; 2016; 51, pp. 64-74. [DOI: https://dx.doi.org/10.1016/j.humpath.2015.12.021]
59. Wang, S.S.; Trunk, M.; Schiffman, M.; Herrero, R.; Sherman, M.E.; Burk, R.D.; Hildesheim, A.; Bratti, M.C.; Wright, T.; Rodriguez, A.C. et al. Validation of p16INK4a as a marker of oncogenic human papillomavirus infection in cervical biopsies from a population-based cohort in Costa Rica. Cancer Epidemiol. Biomarkers Prev.; 2004; 13, pp. 1355-1360.
60. Ronsini, C.; Anchora, L.P.; Restaino, S.; Fedele, C.; Arciuolo, D.; Teodorico, E.; Bizzarri, N.; Zannoni, G.F.; Ferrandina, G.; Scambia, G. et al. The role of semiquantitative evaluation of lympho-vascular space invasion in early stage cervical cancer patients. Gynecol. Oncol.; 2021; 162, pp. 299-307. [DOI: https://dx.doi.org/10.1016/j.ygyno.2021.06.002]
61. Holowaty, P.; Miller, A.B.; Rohan, T.; To, T. Natural history of dysplasia of the uterine cervix. J. Natl. Cancer Inst.; 1999; 91, pp. 252-258. [DOI: https://dx.doi.org/10.1093/jnci/91.3.252] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10037103]
62. van Oortmarssen, G.J.; Habbema, J.D. Epidemiological evidence for age-dependent regression of pre-invasive cervical cancer. Br. J. Cancer; 1991; 64, pp. 559-565. [DOI: https://dx.doi.org/10.1038/bjc.1991.350]
63. Appleby, P.; Beral, V.; Berrington de González, A.; Colin, D.; Franceschi, S.; Goodill, A.; Green, J.; Peto, J.; Plummer, M.; Sweetland, S. Carcinoma of the cervix and tobacco smoking: Collaborative reanalysis of individual data on 13,541 women with carcinoma of the cervix and 23,017 women without carcinoma of the cervix from 23 epidemiological studies. Int. J. Cancer; 2006; 118, pp. 1481-1495. [DOI: https://dx.doi.org/10.1002/ijc.21493]
64. Sand, F.L.; Munk, C.; Frederiksen, K.; Junge, J.; Iftner, T.; Dehlendorff, C.; Kjaer, S.K. Risk of CIN3 or worse with persistence of 13 individual oncogenic HPV types. Int. J. Cancer; 2019; 144, pp. 1975-1982. [DOI: https://dx.doi.org/10.1002/ijc.31883]
65. International Collaboration of Epidemiological Studies of Cervical Cancer Appleby, P.; Beral, V.; Berrington de Gonzalez, A.; Colin, D.; Franceschi, S.; Goodhill, A.; Green, J.; Peto, J.; Plummer, M. et al. Cervical cancer and hormonal contraceptives: Collaborative reanalysis of individual data for 16,573 women with cervical cancer and 35,509 women without cervical cancer from 24 epidemiological studies. Lancet; 2007; 370, pp. 1609-1621. [DOI: https://dx.doi.org/10.1016/S0140-6736(07)61684-5]
66. International Collaboration of Epidemiological Studies of Cervical Cancer. Cervical carcinoma and reproductive factors: Collaborative reanalysis of individual data on 16,563 women with cervical carcinoma and 33,542 women without cervical carcinoma from 25 epidemiological studies. Int. J. Cancer; 2006; 119, pp. 1108-1124. [DOI: https://dx.doi.org/10.1002/ijc.21953]
67. Salvadó, A.; Miralpeix, E.; Solé-Sedeno, J.M.; Kanjou, N.; Lloveras, B.; Duran, X.; Mancebo, G. Predictor factors for conservative management of cervical intraepithelial neoplasia grade 2: Cytology and HPV genotyping. Gynecol. Oncol.; 2021; 162, pp. 569-574. [DOI: https://dx.doi.org/10.1016/j.ygyno.2021.06.019]
68. Mitra, A.; MacIntyre, D.A.; Ntritsos, G.; Smith, A.; Tsilidis, K.K.; Marchesi, J.R.; Bennett, P.R.; Moscicki, A.B.; Kyrgiou, M. The vaginal microbiota associates with the regression of untreated cervical intraepithelial neoplasia 2 lesions. Nat. Commun.; 2020; 11, 1999. [DOI: https://dx.doi.org/10.1038/s41467-020-15856-y]
69. Zhang, J.; Lu, C.X. Spontaneous Regression of Cervical Intraepithelial Neoplasia 2: A Meta-analysis. Gynecol. Obstet. Invest.; 2019; 84, pp. 562-567. [DOI: https://dx.doi.org/10.1159/000497286] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31055567]
70. Koeneman, M.M.; Hendriks, N.; Kooreman, L.F.; Winkens, B.; Kruitwagen, R.F.; Kruse, A.J. Prognostic factors for spontaneous regression of high-risk human papillomavirus-positive cervical intra-epithelial neoplasia grade 2. Int. J. Gynecol. Cancer; 2019; 29, pp. 1003-1009. [DOI: https://dx.doi.org/10.1136/ijgc-2019-000343] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31079058]
71. Skorstengaard, M.; Lynge, E.; Suhr, J.; Napolitano, G. Conservative management of women with cervical intraepithelial neoplasia grade 2 in Denmark: A cohort study. BJOG; 2020; 127, pp. 729-736. [DOI: https://dx.doi.org/10.1111/1471-0528.16081] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31880054]
72. Andersen, E.S.; Nielsen, K.; Pedersen, B. The reliability of preconization diagnostic evaluation in patients with cervical intraepithelial neoplasia and microinvasive carcinoma. Gynecol. Oncol.; 1995; 59, pp. 143-147. [DOI: https://dx.doi.org/10.1006/gyno.1995.1281]
73. Mutombo, A.B.; Simoens, C.; Tozin, R.; Bogers, J.; Van Geertruyden, J.P.; Jacquemyn, Y. Efficacy of commercially available biological agents for the topical treatment of cervical intraepithelial neoplasia: A systematic review. Syst. Rev.; 2019; 8, 132. [DOI: https://dx.doi.org/10.1186/s13643-019-1050-4] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31174598]
74. Fonseca, B.O.; Possati-Resende, J.C.; Salcedo, M.P.; Schmeler, K.M.; Accorsi, G.S.; Fregnani, J.; Antoniazzi, M.; Pantano, N.P.; Santana, I.V.V.; Matsushita, G.M. et al. Topical Imiquimod for the Treatment of High-Grade Squamous Intraepithelial Lesions of the Cervix: A Randomized Controlled Trial. Obstet. Gynecol.; 2021; 137, pp. 1043-1053. [DOI: https://dx.doi.org/10.1097/AOG.0000000000004384]
75. Xiong, Y.; Cui, L.; Bian, C.; Zhao, X.; Wang, X. Clearance of human papillomavirus infection in patients with cervical intraepithelial neoplasia: A systemic review and meta-analysis. Medicine; 2020; 99, e23155. [DOI: https://dx.doi.org/10.1097/MD.0000000000023155] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33181688]
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2021 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/). Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
Simple Summary
Through implementation of HPV testing as standard primary screening method, the number of women diagnosed with high-grade cervical intraepithelial neoplasia (CIN2-3) has increased. Although only one third of CIN3 will progress to cancer, conization is standard treatment in high-income countries. The aim of this study was to identify tools with which to predict CIN regression relevant for individualizing treatment within this patient group. We compared the transcriptomic immune-profile from 21 lesions with confirmed regression and 28 lesions with confirmed persistent CIN3. A gene signature with high sensitivity to identify CIN3 lesions that regressed during follow-up was identified. When tested in a cervical cancer cohort (n = 239) with available transcriptomic data, a high regression signature score was associated with favorable survival, small tumors, and immune infiltration. This study presents a gene signature with the capacity to predict CIN regression, that may potentially guide treatment, and identifies common disease drivers in CIN and cervical cancer.
AbstractThe purpose of this study was to establish a gene signature that may predict CIN3 regression and that may aid in selecting patients who may safely refrain from conization. Oncomine mRNA data including 398 immune-related genes from 21 lesions with confirmed regression and 28 with persistent CIN3 were compared. L1000 mRNA data from a cervical cancer cohort was available for validation (n = 239). Transcriptomic analyses identified TDO2 (p = 0.004), CCL5 (p < 0.001), CCL3 (p = 0.04), CD38 (p = 0.02), and PRF1 (p = 0.005) as upregulated, and LCK downregulated (p = 0.01) in CIN3 regression as compared to persistent CIN3 lesions. From these, a gene signature predicting CIN3 regression with a sensitivity of 91% (AUC = 0.85) was established. Transcriptomic analyses revealed proliferation as significantly linked to persistent CIN3. Within the cancer cohort, high regression signature score associated with immune activation by Gene Set enrichment Analyses (GSEA) and immune cell infiltration by histopathological evaluation (p < 0.001). Low signature score was associated with poor survival (p = 0.007) and large tumors (p = 0.01). In conclusion, the proposed six-gene signature predicts CIN regression and favorable cervical cancer prognosis and points to common drivers in precursors and cervical cancer lesions.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Details
; Munk, Ane Cecilie 2 ; Engesæter, Birgit 3
; Akbari, Saleha 4 ; Frafjord, Astri 4
; Hoivik, Erling A 1 ; Forsse, David 1 ; Fasmer, Kristine E 5 ; Woie, Kathrine 6 ; Haldorsen, Ingfrid S 5
; Bertelsen, Bjørn I 7 ; Janssen, Emiel A M 8
; Gudslaugsson, Einar 4 ; Krakstad, Camilla 1
; Øvestad, Irene T 4 1 Centre for Cancer Biomarkers, Department of Clinical Science, University of Bergen, 5053 Bergen, Norway;
2 Department of Obstetrics and Gynaecology, Sørlandet Hospital Kristiansand, 4604 Kristiansand, Norway;
3 Section for Cervical Cancer Screening, Cancer Registry of Norway, 0304 Oslo, Norway;
4 Department of Pathology, Stavanger University Hospital, 4068 Stavanger, Norway;
5 Section for Radiology, Department of Clinical Medicine, University of Bergen, 5021 Bergen, Norway;
6 Department of Obstetrics and Gynaecology, Haukeland University Hospital, 5053 Bergen, Norway;
7 Department of Pathology, Haukeland University Hospital, 5021 Bergen, Norway;
8 Department of Pathology, Stavanger University Hospital, 4068 Stavanger, Norway;