Next Article in Journal
When Bacteria and Viruses Collide: A Tale of Chlamydia trachomatis and Sexually Transmitted Viruses
Next Article in Special Issue
Unraveling Immunological Dynamics: HPV Infection in Women—Insights from Pregnancy
Previous Article in Journal
Real-World Experience of the Comparative Effectiveness and Safety of Combination Therapy with Remdesivir and Monoclonal Antibodies versus Remdesivir Alone for Patients with Mild-to-Moderate COVID-19 and Immunosuppression: A Retrospective Single-Center Study in Aichi, Japan
Previous Article in Special Issue
Diversity of Anal HPV and Non-HPV Sexually Transmitted Infections and Concordance with Genital Infections in HIV-Infected and HIV-Uninfected Women in the Tapajós Region, Amazon, Brazil
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Viruses
Volume 15
Issue 9
10.3390/v15091953
viruses-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The High-Risk Human Papillomavirus Type Influences the Tissue Microenvironment in Cervical Intraepithelial Neoplasia Grade 2

by
Mayumi Saito
1,
Aarthi Rajesh
1,
Carrie Innes
2,
Rachael van der Griend
3,
Peter Fitzgerald
4,
Bryony Simcock
2,
Peter Sykes
2 and
Merilyn Hibma
1,*
1
Department of Pathology, Dunedin School of Medicine, University of Otago, Dunedin 9054, New Zealand
2
Department of Obstetrics and Gynaecology, University of Otago, Christchurch 8011, New Zealand
3
Canterbury Health Laboratories, Christchurch 8011, New Zealand
4
Awanui Labs, Dunedin 9054, New Zealand
*
Author to whom correspondence should be addressed.
Viruses 2023, 15(9), 1953; https://0-doi-org.brum.beds.ac.uk/10.3390/v15091953
Submission received: 23 August 2023 / Revised: 11 September 2023 / Accepted: 13 September 2023 / Published: 19 September 2023
(This article belongs to the Special Issue Immune Responses to Papillomavirus Infections)
Browse Figures
Review Reports Versions Notes

Abstract

:
High-risk, cancer-causing human papillomavirus (HPV) types are associated with cervical precancer and cancer. A high proportion of high-risk HPV precancer lesions undergo immune-mediated regression. The purpose of this study was to determine if the tissue microenvironment of HPV16 and 18 (HPV16/18) cervical intraepithelial neoplasia grade 2 lesions differed from other high-risk types (HPV ‘other’). Consistent with other studies, we found that progression to higher-grade disease was more frequent in HPV16/18 lesions when compared with HPV ‘other’ lesions. HPV16/18 lesions were significantly more likely to be indoleamine 2,3,-dioxygenase 1 (IDO1)-positive and were associated with reduced CD8 and FoxP3 T cells in the lesion. In the stroma, reduced Tbet- and CD32-positive cells and increased Blimp1-positive cells were significantly associated with HPV16/18 lesions when compared with HPV ‘other’ types. On analysis of the IDO1-positive tissues, lesional IDO1 was associated with significantly decreased numbers of CD4-, CD8-, and FoxP3-positive cells in the stroma compared with IDO1-negative tissues. These data suggest that IDO1 expression may impair infiltration of CD4, CD8, and FoxP3 cells into the stroma beneath the precancer lesion. Increased expression of IDO1 may contribute to immune avoidance and an increased frequency of disease progression in HPV16- and 18-positive lesions.

Graphical Abstract

1. Introduction

Human papillomaviruses (HPVs) are double-stranded DNA viruses that exclusively infect cutaneous and mucosal epithelium. There are over 170 different HPV types, divided into five genera. Types included in the alpha genera are of particular interest because of their association with cancer. Thirteen HPV types within the α5 (type 51), α6 (type 56 and 66), α7 (type 18, 45, 39, and 59), and α9 (type 16, 31, 33, 35, 52, and 58) species are classified by the International Approval and Registration Centre (IARC) as Group 1 carcinogens (sufficient evidence for carcinogenicity in humans) [1].
Neoplastic disease associated with HPV infection has well-defined grades of cervical intraepithelial neoplasia (CIN), from CIN1–3. The high-risk (hr) cancer-causing HPV types are more frequently associated with high-grade disease (CIN2 and 3), and progression is more rapid in HPV16/18-positive tissue, compared to tissue positive with other hr HPV types (HPV31/33/45/52/58) [2]. Nearly half of women with HPV16 that persists for over two years develop precancer in the following five years [3].
Hr HPV types mediate several changes to the keratinocytes that they infect, facilitating S-phase progression and proliferation. Cellular changes include degradation of p53 mediated by the viral E6 protein, and degradation of pRB family members by E7 [4]. In addition, viral immune evasion helps facilitate persistent infection. Several immune evasion mechanisms have been reported, including depletion of Langerhans cells, downregulation of major histocompatibility complex I (MHCI), and disruption of the interferon pathway, including via cGAS-STING [5,6,7]. These mechanisms have been identified in vitro or in mouse models and have generally been associated with HPV16 E6 and/or E7 proteins, and in the case of downregulation of STING, both HPV 16 E5 and E7 [8,9]. The contribution of regulation of the tissue and immune microenvironment to disease associated with other hr HPV types is less clear.
Indoleamine 2,3,-dioxygenase 1 (IDO1) is an immune-suppressive molecule that is upregulated by HPV16 E7 [10]. It is an intracellular cytosolic enzyme that mediates local suppression of T cell proliferation and activation by catalysing tryptophan to kynurenine. Tryptophan degradation and accumulation of kynurenine in the extracellular microenvironment suppresses effector T cell proliferation and activation [11,12]. Inhibition of IDO1 function has been shown to elevate tumour-infiltrating CD8 T cells in the mouse [13], and in vitro [14].
The purpose of this study is to determine if hr HPV type (HPV16/18 vs. hr HPV ‘other’) is associated with differences in the tissue and immune microenvironment in CIN2, including differences in expression of IDO1. We found that HPV16/18 was associated with a significantly higher probability of progression over a 24-month follow-up, whereas HPV ‘other’ disease more frequently regressed. HPV16/18 tissues were also significantly more likely to be positive for IDO1. Additionally, fewer CD8- and FoxP3-positive T cells were present in HPV16/18 lesional tissue and Tbet- and CD32-positive cells were less frequent in the stroma when compared with other hr HPV types. When we looked specifically at the IDO1-positive tissues, the numbers of CD4, CD8, and FoxP3+ cells were reduced compared to IDO1-negative tissues. The marked reduction in the numbers of T cells in the lesional and stromal tissue of HPV16/18 lesions is likely to contribute to less frequent immune-mediated regression, when compared with HPV ‘other’ tissues.

2. Materials and Methods

2.1. Participants

The patient cohort in this study are a subset of the participants of a larger observational management study, known as the ‘PRINCess’ trial, and have previously been reported [15,16]. For inclusion in the study reported here, women had to meet the following additional criteria: both hrHPV- and p16INK4A-positive test results were required and CIN2 had been confirmed in the tissue at review by a second pathologist.

2.2. HPV Testing

Participants were separated into HPV16/18-positive or hr HPV ‘other’, based on the result of HPV testing (Abbott Real Time High-risk HPV Assay) of their cytological sample taken prior to biopsy. The Abbott Test detects HPV16, HPV18, HPV ‘other’, and negative. The HPV ‘other’ category includes positivity for one or more of the following 12 HPV types: 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68.

2.3. Immunohistochemical Staining

Biopsies of CIN2 lesions obtained from women that met the criteria outlined above were formalin-fixed, paraffin-embedded, and sectioned to 4 µm thick. The antigen retrieval, staining, and evaluation of the sections for the following markers has been previously described: CD4, Tbet, GATA3, interleukin (IL)-17, FoxP3, CD8, granzyme B, langerin, Fascin, CD11c, CD32, CD138, high mobility group box (HMGB)-1, B-lymphocyte-induced maturation protein (Blimp)-1, thymic stromal lymphopoietin (TSLP), IDO1, and programmed death-ligand-1 (PD-L1) [15].

2.4. Statistical Analysis

Pairwise comparisons were made using the Mann–Whitney U (M–W U)-test. The chi-square and Fisher’s exact tests were used for categorical measures. Survival curves were plotted using CIN3+ as the endpoint, and the distribution of time to the CIN3+ outcome between the two groups was determined using the log-rank test (Mantel–Cox). Hazard ratios were determined using the Mantel–Haenszel test. Univariate analysis was carried out using simple logistic regression to determine the association between markers and HPV type. A Spearman’s correlation analysis was carried out to test the relationship between specific markers that were tested here. All statistical analyses were carried out using GraphPad Prism (Version 9.0).

3. Results

3.1. HPV16-Positive Lesions Were Frequently Also Positive for ‘Other’ hr HPV

This study reports results on a subset of 69 HPV-positive participants in the PRINCess study. Thirty of the women in this study were categorised as HPV16/18-positive, and 39 were HPV ‘other’. Of the HPV16/18-positive women, 15 women had only HPV16 (14) or HPV18 (1) and one woman had both HPV16 and HPV18. The remaining 14 HPV16/18 (13 HPV16, 1 HPV18) women had coinfections with other hr HPV types.
The clinical characteristics of the participants separated by HPV16/18 or HPV ’other’ high-risk types are shown in Table 1. The entry criteria for the PRINCess study were women under the age of 25 years old who had diagnosed CIN2, which is reflected in the narrow age range of 17–24 years in both the HPV16/18 and HPV ‘other’ cohorts described here. While smoking rates were higher in the HPV16/18 group compared with the HPV ’other’ group, this difference was not statistically significant. The HPV type (HPV16/18 or ‘other’) that was detected in participants had no significant impact on the reported lesion size.
Expectedly, significantly few vaccinated women were in the HPV16/18 group compared with the HPV ‘other’ group. Vaccination for HPV16 and 18 was introduced in New Zealand in 2008 and the change in genotypes associated with vaccination in the PRINCess study, from which this subset of participants was obtained, has been described in detail elsewhere [17].
While the majority of hrHPV+-vaccinated women were positive for HPV ’other’ rather than HPV16/18 in these cohorts, two vaccinated women were HPV16/18-positive. Of those women, one was positive for HPV16 and another non-HPV18 high-risk type and the other was HPV16- and 18-positive. Both individuals had disease that subsequently became more severe (CIN3+) during follow-up.
There was little apparent evidence for any cross-protection between the HPV16/18 vaccine and the HPV ‘other’ types, as a significantly more substantial proportion of the HPV ‘other’ women were vaccinated and had high-grade disease (CIN2).

3.2. HPV16/18-Positive Women Are More Likely to Progress to More Severe Disease than HPV ‘Other’-Positive Women

Participants in this study were CIN2 on entry and were monitored for 24 months. Women were treated and withdrawn from the study if disease became more severe (CIN3+) during that time. The progression to the CIN3+ endpoint for HPV16/18 and HPV ‘other’ women is shown in Figure 1. HPV16/18 women were 4.5 times more likely to progress to the CIN3+ endpoint than women infected with HPV ‘other’ types. This result is consistent with other studies that have identified the significantly greater risk of disease progression if women are infected with HPV16/18 compared with other high-risk types.

3.3. E1^E4 Staining Was Infrequent in hrHPV-Positive CIN2 Lesions

The HPV E1^E4 protein is a product of E1 and E4 spliced mRNA and is an abundant protein in the viral lifecycle during the productive phase. HPV16 E1^E4 induces G2/M cell cycle arrest and induces apoptosis in keratinocytes, and is expressed in the upper stratified layers of the epidermis during HPV infection [18,19,20]. While p16 overexpression is a direct result of HPV E7 degradation of pRb, the nuclear expression of Ki67 occurs in the nucleus during cell division as a result of cell proliferation. Ki67 expression is confined to the basal lower third of normal epidermis; however, expression in the upper layers is associated with increasing grades of CIN [21].
We looked at expression of the E1^E4 protein in HPV16/18 and HPV ‘other’ tissues and its association with Ki67 staining (Table 2). Overall, E1^E4 staining (16.9% positivity) was infrequent in the CIN2 lesions. E1^E4 staining was more frequent in HPV ‘other’ compared with HPV16/18 lesions (albeit not significantly; p = 0.08; Fisher’s exact test). Of the lesions with Ki67 staining in the lower third only, none of the HPV16/18-positive lesions were positive for E1^E4, compared with 57% of the HPV ‘other’ lesions.

3.4. IDO1 Expression Is More Frequent in HPV16/18-Positive Lesions

IDO1 is an enzyme that is overexpressed in many cancer cells and is associated with immunosuppression [22]. PD-L1 is also overexpressed in some types of cancers and is a successfully targeted checkpoint molecule [23]. As both molecules play a role in cancer immune escape, and expression of either of these proteins is generally associated with a poor prognosis in cancer, we wished to determine if expression was increased in HPV16/18 versus HPV ‘other’ infected lesions.
In positive tissues, IDO1 expression in the epidermal cells of the lesion was generally strong, punctate, and perinuclear (Figure 2). There was also weaker diffuse cytoplasmic staining in the majority of the positive cells. Significantly (Fisher’s exact test; p < 0.0001), we found that approximately half of the HPV16/18 lesions were positive for IDO1, and that IDO1 expression on HPV ‘other’ tissues was rare (Table 3). In contrast, less than 20% of the tissues were PD-L1-positive, irrespective of hr HPV type.
We also looked at coexpression of PD-L1 and IDO1. While noting that the number of PD-L1-positive tissues was low, all the PD-L1-positive HPV16/18 cases (n = 5) were also positive for IDO1, whereas all of the PD-L1-positive cases in HPV ‘other’ (n = 7) were IDO1-negative.

3.5. Immune Cell Infiltrates Were Reduced in HPV16/18 Lesions Compared with HPV ‘Other’ Lesions

We carried out a detailed analysis of the immune microenvironment of the HPV16/18 and HPV ‘other’ tissues (Table 4). The markers tested were primarily to characterise T cells, B cells, and antigen presenting cells. Additionally, expression of the immune regulatory molecules HMGB1 and Blimp-1 were also assessed.
Significantly, around half the number of CD8 and CD4+FoxP3+ T cells were present in HPV16/18 compared with HPV ‘other’ lesional tissue. Although not statistically significant, for several of the other markers (CD4, granzyme B, Langerin, CD32, and CD138), there were, on average, more positive cells in the HPV ‘other’ compared with HPV16/18 lesional tissue.
When we examined the lamina propria below HPV ‘other’ and HPV16/18 lesions (Table 5), we found that numbers of Tbet- and CD32-positive cells were significantly lower in HPV16/18 compared with HPV ‘other’ lesions. In contrast, Blimp-1-positive cells were significantly more likely to be detected in the lamina propria beneath HPV16/18 lesions when compared with HPV ‘other’ lesions.

3.6. The Correlations between Cells Differ in HPV16/18 and HPV ‘Other’ Tissues

We examined the relationships between different cell types in the lesion and the lamina propria beneath the lesion, in particular looking for marked differences between HPV16/18 and HPV ‘other’ tissues (Figure 3). The presence of CD4+FoxP3+ T cells directly correlated most strongly with HMGB1-positive cells in HPV16/18 lesions but more strongly with Blimp-1-positive cells in HPV ‘other’ lesions. Additionally, there was a strong correlation between CD4+ and IL17+ cells, and CD4+ and CD8+ cells, in HPV ‘other’ lesions, which was not apparent in HPV16/18 lesions.
In the lamina propria, Blimp-1-positive cells were inversely correlated with CD32+ cells, irrespective of hrHPV type. CD4-, CD8-, and granzyme B-positive (but not double-positive) cells were correlated with CD11c-positive cells in the HPV ‘other’, but not HPV16/18 lamina propria. Conversely, IL17-positive cells were inversely related with FoxP3+ and CD4+FoxP3+ cells in the HPV16/18 and not HPV ‘other’ lamina propria.

3.7. The Lamina Propria Beneath IDO1-Positive Tissue Has Fewer CD4, CD8 T Cells and FoxP3-Positive Cells

We carried out a comparison of the cell populations between IDO1-positive and IDO1-negative lesions in the hr HPV-positive tissues (Table 5). While there was no significant difference in the frequency of particular cell subsets in the lesional tissue, there was a marked reduction in the numbers of cells beneath the lesion. Specifically, the CD4 and CD8 T cells were around half as frequent in the lamina propria beneath the IDO1-negative tissues compared to IDO1-positive tissues. Perhaps surprisingly, the frequency of FoxP3-positive cells in IDO1-positive tissue was around a quarter of that of IDO1-negative tissue. This reduction was also reflected in both the CD4+FoxP3+ and the CD4-FoxP3+ subsets.

4. Discussion

In this study, we report the effects of hr HPV type (HPV16/18 and HPV ‘other’) on the immune microenvironment of p16+ CIN2 lesions. We found that HPV16/18 CIN2 lesions, which were more likely to progress, are significantly more frequently IDO1-positive than other hr HPV-positive CIN2 lesions. Furthermore, we found an HPV16/18-associated decrease in lesional CD8+ and CD4+FoxP3+ cells and stromal Tbet+ and CD32+ cells, when compared with other hr HPV types. Additionally, Blimp1 expression in the stroma was increased in the HPV16/18 compared with HPV ‘other’ lesions.
While some HPV infections persist beyond two years, most infections undergo immune-mediated regression. In women with CIN2, others have reported that disease regresses in 50% of cases within two years [24]. In the patient cohort studied here, disease in at least 53% of the women under 25 with CIN2 spontaneously regressed [25]. Importantly, progression was more likely in HPV16/18 infection than with other hr HPV types. A contributing factor may be that the tissue microenvironment of HPV16/18 lesions is not conducive to an effective antiviral immune response.
We only included p16-positive, block-stained lesions in this study. Block staining for p16 (i.e., continuous staining of the cells of the basal and parabasal layers of cervical squamous epithelium, with or without staining in the intermediate or superficial layers) is strongly associated with high-grade hr HPV-positive disease [26]. Increased p16 expression is a result of a pRB-independent increase by E7 of the expression of the histone demethylase KDM6B. Consequential demethylation of the repressive histone mark H3K27me3 occurs, resulting in the induction of p16 [27].
While p16 would typically induce pRB-dependent senescence, E7 degrades pRB, and cells are able to proliferate. Increased proliferating cells, as evidenced by Ki67 staining, was seen in the majority (around 80%) of stained tissues irrespective of HPV type. Considering that only around 15% of the tissues were positive for E1^E4 and that p16 upregulation is a consequence of E7, the viral early promoter, which controls E7, is likely to be more active in these tissues than the late promoter, which drives E1^E4. This is consistent with the predicted cellular changes in viral gene expression that are associated with progression of cells to a more transformed phenotype [28].
We found that IDO1 was significantly more frequently expressed in HPV16/18 lesions compared with hr HPV ‘other’ lesions. In normal cervix IDO1 is expressed only in rare glandular cells but is frequently expressed in cervical carcinomas [29]. Heeren et al. (2018) reported IDO1 expression in primary squamous carcinoma cells and stromal immune cells, as well as in metastatic IDO1+ tumour cells in the lymph nodes [30]. In a study of liquid-based cytology samples, IDO1 was significantly increased in HSIL squamous cells compared with control samples. Consistent with what we observed here, IDO1 staining was detected in the upper layers of the epidermis in CIN [31]. In contrast, in head and neck squamous cell carcinoma, IDO1 expression was most commonly at the focal interface between the tumour and the stroma and was not associated with HPV status (HPV-positive versus -negative) [32].
How IDO1 expression is regulated by HPV16/18 is less clear. Increased IDO1 expression is typically associated with a loss of BIN1 or overexpression of COX2. COX2 expression is increased in cervical cancer [33]. HPV16 E6 and E7 expression induces COX-2 transcription by activating the epidermal growth factor receptor in vitro, and the induction of COX-2 expression may be the mechanism by which HPV16/18 increases IDO1 expression. The infrequent IDO1 expression in other hr HPV lesions suggests that viral regulation of IDO1 may be restricted to only selected hr HPV types.
We found that CD8+ T cells in the lesion and Tbet+ cells in the stroma were reduced in HPV16/18 lesions, compared with hr HPV ‘other’ types. Consistent with these observations, commonly reported high-grade cervical disease associated cellular changes when compared with normal tissue include decreased CD4 and CD8 cells and increased IDO1-positive cells [34]. Changes associated with subsequent disease progression from CIN2 to CIN3 include decreased CD4+, T bet-positive, CD8+, and CD11c+ cells [15]. Overall, the reduced numbers of T cells in HPV16/18 compared with hr HPV ‘other’ tissues in this study may reflect inhibition of T cell proliferation in the local tissue microenvironment. The frequent expression of IDO1 in HPV16/18 tissues may contribute to the reduced numbers of CD8 and CD4+FoxP3+ T cells that proliferated in these lesions. We also found significantly reduced numbers of CD4 and CD8 T cells in the stroma of IDO1-positive tissues, suggesting that the effects on T cell proliferation mediated through IDO1 may extend into the stroma beneath the lesion.
In the skin-expressing HPV16 E7 transgenic mouse model, inhibition of IDO1 actively promoted rejection of E7-skin grafted onto an immune competent mouse that would not otherwise be rejected [10]. This indicates a direct role for HPV16 E7 in the regulation of IDO1, and direct immunosuppression mediated by this molecule. We propose that the more frequent expression of IDO1 in HPV16/18 compared with HPV ‘other’ types in the CIN2 lesions may reduce the frequency of immune-mediated regression in those tissues.
Counterintuitively, we found that higher numbers of CD4+FoxP3+ cells were significantly associated with HPV ‘other’ lesions, whereas HPV16/18 lesions had reduced numbers of these cells. In HPV16/18 lesions, CD4+FoxP3+ cells were strongly correlated with HMGB1, whereas they were strongly correlated with Blimp1 expression in HPV ‘other’ lesions. Previously, we did not find any association between CD4+FoxP3+ cells and disease progression in CIN2+ lesions but did detect an increase in CD4-FoxP3+ cells in lesions that subsequently progressed to higher-grade disease [15]. Furthermore, we found that Blimp1 and HMGB1 expression in the lesion was associated with a higher likelihood of disease progression. Additionally, others have reported that Treg are increased in the stroma in high-grade disease [34]. CD4+FoxP3+ cells demonstrate functional plasticity [35]. We note that CD4+FoxP3+ cells were detected at a low frequency in the lesional tissue (around 10-fold less than the CD8 T cells), and it is unclear what functional impact this comparatively low number of cells may have.

5. Conclusions

Cell populations in the tissue microenvironment of CIN2 lesions differ depending on the hr HPV type(s) present (HPV16/18 vs. hr HPV ‘other’). Disease progression is more frequent in HPV16/18-positive lesions and is associated with IDO1 positivity, reduced lesional CD8 cells and lamina propria Tbet- and CD32-positive cells, and increased Blimp1 expression. Surprisingly, CD4+FoxP3-positive cells were also reduced in the lesion and stroma of HPV16/18 CIN2 tissue when compared with HPV ‘other’. We propose that the increased expression of IDO1 and reduction in T and B cells contributes to a reduced likelihood of immune-mediated regression in HPV16/18 tissues. The consequences of the reduction in CD4+FoxP3+-positive cells in HPV16/18 CIN2 stroma are yet to be elucidated.

Author Contributions

Conceptualization, M.H., M.S. and P.S.; methodology, M.S. and M.H.; validation, M.S., P.F., R.v.d.G., C.I. and M.H.; formal analysis, M.S. and M.H.; investigation, M.S., P.S., M.H., C.I., B.S., R.v.d.G., P.F. and A.R.; data curation, M.S. and C.I.; writing—original draft preparation, M.H. and M.S.; writing—review and editing, M.H., M.S., A.R., C.I., R.v.d.G., P.F., B.S. and P.S.; funding acquisition, M.H. and P.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Genesis Oncology Trust (2016).

Institutional Review Board Statement

The study was conducted in accordance with the Declaration of Helsinki and approved by the Multi-region Ethics Committee of the Health and Disability Ethics Committees, New Zealand Government (MEC/09/07/079, approved on 14 April 2010). Informed consent was obtained from all subjects involved in the study.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The datasets analysed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We gratefully acknowledge the kind gift of the pan-E1^E4 antibody from John Doorbar.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Iarc Working Group on the Evaluation of Carcinogenic Risks to Humans. IARC monographs on the evaluation of carcinogenic risks to humans. Biological Agents. IARC Monogr. Eval. Carcinog. Risks Hum. 2012, 100B, 255–314. [Google Scholar]
  2. Karube, A.; Saito, F.; Waga, M.; Yokoyama, S.; Kanamori, K. Progression of cervical intraepithelial neoplasia grade 2 lesions among Japanese women harboring different genotype categories of high-risk human papillomaviruses. J. Rural. Med. 2021, 16, 91–97. [Google Scholar] [CrossRef]
  3. Kjaer, S.K.; Frederiksen, K.; Munk, C.; Iftner, T. Long-term absolute risk of cervical intraepithelial neoplasia grade 3 or worse following human papillomavirus infection: Role of persistence. J. Natl. Cancer Inst. 2010, 102, 1478–1488. [Google Scholar] [CrossRef] [PubMed]
  4. Doorbar, J.; Egawa, N.; Griffin, H.; Kranjec, C.; Murakami, I. Human papillomavirus molecular biology and disease association. Rev. Med. Virol. 2015, 25 (Suppl. 1), 2–23. [Google Scholar] [CrossRef]
  5. Matthews, K.; Leong, C.M.; Baxter, L.; Inglis, E.; Yun, K.; Backstrom, B.T.; Doorbar, J.; Hibma, M. Depletion of Langerhans cells in human papillomavirus type 16-infected skin is associated with E6-mediated down regulation of E-cadherin. J. Virol. 2003, 77, 8378–8385. [Google Scholar] [CrossRef] [PubMed]
  6. Kim, D.H.; Kim, E.M.; Lee, E.H.; Ji, K.Y.; Yi, J.; Park, M.; Kim, K.D.; Cho, Y.Y.; Kang, H.S. Human papillomavirus 16E6 suppresses major histocompatibility complex class I by upregulating lymphotoxin expression in human cervical cancer cells. Biochem. Biophys. Res. Commun. 2011, 409, 792–798. [Google Scholar] [CrossRef] [PubMed]
  7. Scott, M.L.; Woodby, B.L.; Ulicny, J.; Raikhy, G.; Orr, A.W.; Songock, W.K.; Bodily, J.M. Human Papillomavirus 16 E5 Inhibits Interferon Signaling and Supports Episomal Viral Maintenance. J. Virol. 2020, 94, e01582-19. [Google Scholar] [CrossRef] [PubMed]
  8. Bortnik, V.; Wu, M.; Julcher, B.; Salinas, A.; Nikolic, I.; Simpson, K.J.; McMillan, N.A.; Idris, A. Loss of HPV type 16 E7 restores cGAS-STING responses in human papilloma virus-positive oropharyngeal squamous cell carcinomas cells. J. Microbiol. Immunol. Infect. 2021, 54, 733–739. [Google Scholar] [CrossRef]
  9. Miyauchi, S.; Kim, S.S.; Jones, R.N.; Zhang, L.; Guram, K.; Sharma, S.; Schoenberger, S.P.; Cohen, E.E.W.; Califano, J.A.; Sharabi, A.B. Human papillomavirus E5 suppresses immunity via inhibition of the immunoproteasome and STING pathway. Cell Rep. 2023, 42, 112508. [Google Scholar] [CrossRef]
  10. Mittal, D.; Kassianos, A.J.; Tran, L.S.; Bergot, A.S.; Gosmann, C.; Hofmann, J.; Blumenthal, A.; Leggatt, G.R.; Frazer, I.H. Indoleamine 2,3-dioxygenase activity contributes to local immune suppression in the skin expressing human papillomavirus oncoprotein e7. J. Investig. Dermatol. 2013, 133, 2686–2694. [Google Scholar] [CrossRef]
  11. Munn, D.H.; Shafizadeh, E.; Attwood, J.T.; Bondarev, I.; Pashine, A.; Mellor, A.L. Inhibition of T cell proliferation by macrophage tryptophan catabolism. J. Exp. Med. 1999, 189, 1363–1372. [Google Scholar] [CrossRef]
  12. Munn, D.H.; Zhou, M.; Attwood, J.T.; Bondarev, I.; Conway, S.J.; Marshall, B.; Brown, C.; Mellor, A.L. Prevention of allogeneic fetal rejection by tryptophan catabolism. Science 1998, 281, 1191–1193. [Google Scholar] [CrossRef]
  13. Schafer, C.C.; Wang, Y.; Hough, K.P.; Sawant, A.; Grant, S.C.; Thannickal, V.J.; Zmijewski, J.; Ponnazhagan, S.; Deshane, J.S. Indoleamine 2,3-dioxygenase regulates anti-tumor immunity in lung cancer by metabolic reprogramming of immune cells in the tumor microenvironment. Oncotarget 2016, 7, 75407–75424. [Google Scholar] [CrossRef]
  14. Liu, H.; Liu, L.; Liu, K.; Bizargity, P.; Hancock, W.W.; Visner, G.A. Reduced cytotoxic function of effector CD8+ T cells is responsible for indoleamine 2,3-dioxygenase-dependent immune suppression. J. Immunol. 2009, 183, 1022–1031. [Google Scholar] [CrossRef] [PubMed]
  15. Saito, M.; Rajesh, A.; Innes, C.; van der Griend, R.; Fitzgerald, P.; Simcock, B.; Sykes, P.; Hibma, M. Blimp-1 is a prognostic indicator for progression of cervical intraepithelial neoplasia grade 2. J. Cancer Res. Clin. Oncol. 2022, 148, 1991–2002. [Google Scholar] [CrossRef] [PubMed]
  16. Sykes, P.; Innes, C.; Harker, D.; Whitehead, M.; van der Griend, R.; Lawton, B.; Hibma, M.; Fitzgerald, P.; Dudley, N.; Petrich, S.; et al. Observational Management of CIN 2 in Young Women: A Prospective Multicenter Trial. J. Low. Genit. Tract. Dis. 2016, 20, 343–347. [Google Scholar] [CrossRef]
  17. Innes, C.R.; Sykes, P.H.; Harker, D.; Williman, J.A.; Van der Griend, R.A.; Whitehead, M.; Hibma, M.; Lawton, B.A.; Fitzgerald, P.; Dudley, N.M.; et al. Changes in human papillomavirus genotypes associated with cervical intraepithelial neoplasia grade 2 lesions in a cohort of young women (2013–2016). Papillomavirus Res. 2018, 6, 77–82. [Google Scholar] [CrossRef] [PubMed]
  18. Davy, C.E.; Jackson, D.J.; Raj, K.; Peh, W.L.; Southern, S.A.; Das, P.; Sorathia, R.; Laskey, P.; Middleton, K.; Nakahara, T.; et al. Human papillomavirus type 16 E1 E4-induced G2 arrest is associated with cytoplasmic retention of active Cdk1/cyclin B1 complexes. J. Virol. 2005, 79, 3998–4011. [Google Scholar] [CrossRef]
  19. Raj, K.; Berguerand, S.; Southern, S.; Doorbar, J.; Beard, P. E1 empty set E4 protein of human papillomavirus type 16 associates with mitochondria. J. Virol. 2004, 78, 7199–7207. [Google Scholar] [CrossRef] [PubMed]
  20. Doorbar, J.; Foo, C.; Coleman, N.; Medcalf, L.; Hartley, O.; Prospero, T.; Napthine, S.; Sterling, J.; Winter, G.; Griffin, H. Characterization of events during the late stages of HPV16 infection in vivo using high-affinity synthetic Fabs to E4. Virology 1997, 238, 40–52. [Google Scholar] [CrossRef]
  21. Keating, J.T.; Cviko, A.; Riethdorf, S.; Riethdorf, L.; Quade, B.J.; Sun, D.; Duensing, S.; Sheets, E.E.; Munger, K.; Crum, C.P. Ki-67, cyclin E, and p16INK4 are complimentary surrogate biomarkers for human papilloma virus-related cervical neoplasia. Am. J. Surg. Pathol. 2001, 25, 884–891. [Google Scholar] [CrossRef]
  22. Tang, K.; Wu, Y.H.; Song, Y.; Yu, B. Indoleamine 2,3-dioxygenase 1 (IDO1) inhibitors in clinical trials for cancer immunotherapy. J. Hematol. Oncol. 2021, 14, 68. [Google Scholar] [CrossRef]
  23. Akinleye, A.; Rasool, Z. Immune checkpoint inhibitors of PD-L1 as cancer therapeutics. J. Hematol. Oncol. 2019, 12, 92. [Google Scholar] [CrossRef] [PubMed]
  24. Tainio, K.; Athanasiou, A.; Tikkinen, K.A.O.; Aaltonen, R.; Hernandes, J.C.; Glazer-Livson, S.; Jakobsson, M.; Joronen, K.; Kiviharju, M.; Louvanto, K.; et al. Clinical course of untreated cervical intraepithelial neoplasia grade 2 under active surveillance: Systematic review and meta-analysis. BMJ 2018, 360, k499. [Google Scholar] [CrossRef] [PubMed]
  25. Sykes, P.H.; Simcock, B.J.; Innes, C.R.; Harker, D.; Williman, J.A.; Whitehead, M.; van der Griend, R.A.; Lawton, B.A.; Hibma, M.; Fitzgerald, P.; et al. Predicting regression of cervical intraepithelial neoplasia grade 2 in women under 25 years. Am. J. Obstet. Gynecol. 2022, 226, 222.e1–222.e13. [Google Scholar] [CrossRef]
  26. Darragh, T.M.; Colgan, T.J.; Thomas Cox, J.; Heller, D.S.; Henry, M.R.; Luff, R.D.; McCalmont, T.; Nayar, R.; Palefsky, J.M.; Stoler, M.H.; et al. The Lower Anogenital Squamous Terminology Standardization project for HPV-associated lesions: Background and consensus recommendations from the College of American Pathologists and the American Society for Colposcopy and Cervical Pathology. Int. J. Gynecol. Pathol. 2013, 32, 76–115. [Google Scholar] [CrossRef] [PubMed]
  27. McLaughlin-Drubin, M.E.; Crum, C.P.; Munger, K. Human papillomavirus E7 oncoprotein induces KDM6A and KDM6B histone demethylase expression and causes epigenetic reprogramming. Proc. Natl. Acad. Sci. USA 2011, 108, 2130–2135. [Google Scholar] [CrossRef]
  28. Doorbar, J. Host control of human papillomavirus infection and disease. Best. Pract. Res. Clin. Obstet. Gynaecol. 2018, 47, 27–41. [Google Scholar] [CrossRef]
  29. Theate, I.; van Baren, N.; Pilotte, L.; Moulin, P.; Larrieu, P.; Renauld, J.C.; Herve, C.; Gutierrez-Roelens, I.; Marbaix, E.; Sempoux, C.; et al. Extensive profiling of the expression of the indoleamine 2,3-dioxygenase 1 protein in normal and tumoral human tissues. Cancer Immunol. Res. 2015, 3, 161–172. [Google Scholar] [CrossRef]
  30. Heeren, A.M.; van Dijk, I.; Berry, D.; Khelil, M.; Ferns, D.; Kole, J.; Musters, R.J.P.; Thijssen, V.L.; Mom, C.H.; Kenter, G.G.; et al. Indoleamine 2,3-Dioxygenase Expression Pattern in the Tumor Microenvironment Predicts Clinical Outcome in Early Stage Cervical Cancer. Front. Immunol. 2018, 9, 1598. [Google Scholar] [CrossRef]
  31. 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, 586–597. [Google Scholar] [CrossRef] [PubMed]
  32. Succaria, F.; Kvistborg, P.; Stein, J.E.; Engle, E.L.; McMiller, T.L.; Rooper, L.M.; Thompson, E.; Berger, A.E.; van den Brekel, M.; Zuur, C.L.; et al. Characterization of the tumor immune microenvironment in human papillomavirus-positive and -negative head and neck squamous cell carcinomas. Cancer Immunol. Immunother. 2021, 70, 1227–1237. [Google Scholar] [CrossRef] [PubMed]
  33. Kulkarni, S.; Rader, J.S.; Zhang, F.; Liapis, H.; Koki, A.T.; Masferrer, J.L.; Subbaramaiah, K.; Dannenberg, A.J. Cyclooxygenase-2 is overexpressed in human cervical cancer. Clin. Cancer Res. 2001, 7, 429–434. [Google Scholar] [PubMed]
  34. Muntinga, C.L.P.; de Vos van Steenwijk, P.J.; Bekkers, R.L.M.; van Esch, E.M.G. Importance of the Immune Microenvironment in the Spontaneous Regression of Cervical Squamous Intraepithelial Lesions (cSIL) and Implications for Immunotherapy. J. Clin. Med. 2022, 11, 1432. [Google Scholar] [CrossRef]
  35. Malviya, V.; Yshii, L.; Junius, S.; Garg, A.D.; Humblet-Baron, S.; Schlenner, S.M. Regulatory T-cell stability and functional plasticity in health and disease. Immunol. Cell Biol. 2023, 101, 112–129. [Google Scholar] [CrossRef] [PubMed]
Viruses 15 01953 g001
Figure 1. Survival curve with CIN3+ as the endpoint for HPV16/18 and HPV ‘other’.
Figure 1. Survival curve with CIN3+ as the endpoint for HPV16/18 and HPV ‘other’.
Viruses 15 01953 g001
Viruses 15 01953 g002
Figure 2. Representative examples of positive (A–C) and negative (D) epidermal IDO1 staining in hr HPV-positive CIN2 tissues.
Figure 2. Representative examples of positive (A–C) and negative (D) epidermal IDO1 staining in hr HPV-positive CIN2 tissues.
Viruses 15 01953 g002
Viruses 15 01953 g003
Figure 3. Correlations between variables in hrHPV16/18 and HPV ‘other’ CIN2 lesions (a,b) and lamina propria beneath the lesions (c,d). Blue indicates a direct correlation and red indicates an inverse correlation. P values are indicated on each matrix: * p < 0.05; ** p < 0.01; *** p < 0.001.
Figure 3. Correlations between variables in hrHPV16/18 and HPV ‘other’ CIN2 lesions (a,b) and lamina propria beneath the lesions (c,d). Blue indicates a direct correlation and red indicates an inverse correlation. P values are indicated on each matrix: * p < 0.05; ** p < 0.01; *** p < 0.001.
Viruses 15 01953 g003
Table 1. Baseline characteristics of HPV16/18 and HPV ’other’ groups.
Table 1. Baseline characteristics of HPV16/18 and HPV ’other’ groups.
HPV16/18
(n = 30)		HPV ‘Other’
(n = 39)		Statistical
Significance
AgeMean (Range)21.50 (17–24)21.91 (19–24)p = 0.1043
(M–W U *)
SmokerYes13 (43.3%)11 (28.2%)p = 0.4031
(chi-square)
No17 (56.7%)26 (66.7%)
Unknown-2 (5.1%)
VaccineYes2 (6.7%)25 (64.1%)p < 0.0001
(chi-square)
No12 (40.0%)4 (10.3%)
Unknown16 (53.3%)10 (25.6%)
Lesion areaMean
(Range)		0.4334
(0.0093–2.148)		0.3812
(0.0112–1.672)		p = 0.9462
(M–W U *)
* Mann–Whitney U.
Table 2. E1^E4 and Ki67 expression in hr HPV lesions.
Table 2. E1^E4 and Ki67 expression in hr HPV lesions.
HPV16/18
(n = 30)		HPV ‘Other’
(n = 39)
Ki67E4Ki67E4
1/3 *5 (16.7%)E4+0 (0%)1/38 (21.1%)E4+4 (57.1%)
E4−5 (100%) E4−3 (42.9%)
ND **1 ND1
2/323 (76.7%)E4+3 (13.0%)2/324 (63.2%)E4+4 (17.4%)
E4−20 (87.0%) E4−19 (82.6%)
ND1 ND1
3/32 (6.7%)E4+0 (0%)3/36 (15.8%)E4+0 (0%)
E4−2 (100%) E4−6 (100%)
ND- ND-
ND ND1E4−1
* Ki67 staining in 1/3 (the lower third), 2/3 (the lower and middle thirds), or 3/3 (full thickness) of the epidermis; ** no data.
Table 3. Expression of PD-L1 and IDO1 in HPV16/18 and HPV ‘other’ lesions.
Table 3. Expression of PD-L1 and IDO1 in HPV16/18 and HPV ‘other’ lesions.
HPV16/18
(n = 30)		 HPV ‘Other’
(n = 39)		Fisher’s Exact
Test
PD-L1
Positive5 (16.7%)* <CIN31 (20.0%) 7 (18.4%)<CIN35 (71.4%)p = 0.8505
CIN34 (80.0%) CIN32 (28.6%)
Negative25 (83.3%)<CIN37 (28.0%) 31 (81.6%)<CIN321 (67.7%)
CIN318 (72.0%) CIN310 (32.3%)
IDO1
Positive14 (46.7%)<CIN34 (28.6%) 2 (5.4%)<CIN31 (50.0%)p < 0.0001
CIN310 (71.4%) CIN31 (50.0%)
Negative16 (53.3%)<CIN34 (25.0%) 35 (94.6%)<CIN324 (68.6%)
CIN3 12 (75.0%) CIN311 (31.4%)
No data were available for one HPV ‘other’ PD-L1 sample and two HPV ‘other’ IDO1 samples. * <CIN3 indicates persistence or regression at follow-up; CIN3 indicates progression to CIN3 at follow-up.
Table 4. Cells in the epidermis and lamina propria of HPV16/18- and HPV ‘other’-positive tissues.
Table 4. Cells in the epidermis and lamina propria of HPV16/18- and HPV ‘other’-positive tissues.
Epidermis Lamina Propria
Cells/mm2 (Mean ± SD) Cells/mm2 (Mean ± SD)
HPV16/18HPV ‘Other’p Value #HPV16/18HPV ‘Other’p Value #
CD4298.13 ± 237.48407.11 ± 380.940.20722008.57 ± 1986.881993.33 ± 1506.650.7521
Tbet88.74 ± 117.54151.09 ± 146.810.1013371.78 ± 276.78834.04 ± 945.680.0335 *
GATA31579.15 ± 1968.701223.91 ± 1346.170.71893009.15 ± 2279.402399.54 ± 1665.510.3533
IL17489.24 ± 664.42380.62 ± 481.350.62731280.58 ± 680.081613.62 ± 851.250.1329
FoxP327.23 ± 100.0924.61 ± 45.770.416456.27 ± 69.8379.59 ± 79.140.2534
FoxP3+CD4−21.35 ± 99.1612.49 ± 39.860.874217.21 ± 28.9417.50 ± 24.480.6946
FoxP3+CD4+5.15 ± 9.1512.12 ± 14.540.0151 *39.06 ± 50.1262.09 ± 65.940.1555
CD8125.09 ± 104.54246.16 ± 210.550.0086 **609.07 ± 454.41888.83 ± 804.860.1275
Granzyme B95.85 ± 127.42129.85 ± 144.740.2477411.28 ± 417.18605.03 ± 817.210.3997
CD8+ GranzymeB+0.54 ± 1.271.89 ± 4.390.194411.83 ± 23.9166.75 ± 227.140.1811
Langerin+21.70 ± 26.7838.41 ± 58.940.1942.92 ± 5.546.05 ± 9.650.2664
Langerin+ Fascin+12.48 ± 14.7626.91 ± 43.320.1931.10 ± 2.972.36 ± 5.650.2948
CD11c87.07 ± 66.39109.04 ± 126.700.8521545.86 ± 1621.87353.28 ± 341.680.9857
CD32195.57 ± 473.50408.62 ± 866.340.064480.04 ± 74.57235.32 ± 297.030.0052 **
CD138740.95 ± 642.071043.74 ± 1912.310.4895845.84 ± 715.601706.27 ± 2549.270.1554
HMGB11615.30 ± 1559.871534.52 ± 1364.220.74633790.87 ± 1989.523841.23 ± 1995.400.9285
Blimp1537.32 ± 996.82343.93 ± 1094.080.10971760.00 ± 1563.051033.87 ± 1392.860.0447 *
TSLPa NDNDND1205.08 ± 793.24941.04 ± 647.560.207
Area0.43 ± 0.490.38 ± 0.420.748NDNDND
a ND: not detectable. * p < 0.05, ** p < 0.01; # simple logistic regression analysis of square root transformed data.
Table 5. Immune cell densities in IDO1-positive and IDO1-negative hr HPV-positive lesions.
Table 5. Immune cell densities in IDO1-positive and IDO1-negative hr HPV-positive lesions.
Epidermis Lamina Propria
IDO+
Cells/mm2		IDO−
Cells/mm2		p Value #IDO+
Cells/mm2		IDO−
Cells/mm2		p Value #
CD4304.57377.51 0.4081075.66 2265.73 0.01 **
Tbet92.61134.41 0.6268344.55 740.23 0.1176
GATA31384.551405.66 0.91362999.24 2546.18 0.4573
IL17426.39425.44 0.70621744.54 1377.64 0.1172
FoxP310.5530.84 0.500621.80 85.69 0.0006 ***
FoxP3+CD4−3.5321.10 0.56853.95 22.40 0.0034 **
FoxP3+CD4+6.799.74 0.578117.85 63.29 0.0023 **
CD8163.56201.05 0.6991477.46 874.07 0.0304 *
Granzyme B93.63 117.03 0.6063357.19 585.70 0.4148
CD8+ GranzymeB+0.88 1.45 0.750810.73 58.24 0.1161
Langerin+22.00 34.30 0.39682.50 5.37 0.4949
Langerin+ Fascin+14.42 22.83 0.66461.29 1.83 0.6469
CD11c83.50 105.68 0.6951193.95 506.50 0.084
CD32195.57 367.30 0.411991.10 197.40 0.1652
CD138616.64 1014.23 0.3613691.26 1549.87 0.1375
HMGB11823.29 1448.84 0.16343796.03 3778.15 0.7355
Blimp1943.97 281.20 0.07961634.13 1236.52 0.2363
TSLPa NDNDND1105.43 1028.97 0.7995
Area0.41 0.40 0.5534NDNDND
a ND: not detectable. * p < 0.05, ** p < 0.01, *** p < 0.001. # Univariate analysis on square root transformed data.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saito, M.; Rajesh, A.; Innes, C.; van der Griend, R.; Fitzgerald, P.; Simcock, B.; Sykes, P.; Hibma, M. The High-Risk Human Papillomavirus Type Influences the Tissue Microenvironment in Cervical Intraepithelial Neoplasia Grade 2. Viruses 2023, 15, 1953. https://0-doi-org.brum.beds.ac.uk/10.3390/v15091953

AMA Style

Saito M, Rajesh A, Innes C, van der Griend R, Fitzgerald P, Simcock B, Sykes P, Hibma M. The High-Risk Human Papillomavirus Type Influences the Tissue Microenvironment in Cervical Intraepithelial Neoplasia Grade 2. Viruses. 2023; 15(9):1953. https://0-doi-org.brum.beds.ac.uk/10.3390/v15091953

Chicago/Turabian Style

Saito, Mayumi, Aarthi Rajesh, Carrie Innes, Rachael van der Griend, Peter Fitzgerald, Bryony Simcock, Peter Sykes, and Merilyn Hibma. 2023. "The High-Risk Human Papillomavirus Type Influences the Tissue Microenvironment in Cervical Intraepithelial Neoplasia Grade 2" Viruses 15, no. 9: 1953. https://0-doi-org.brum.beds.ac.uk/10.3390/v15091953

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

Article Metrics

Back to TopTop