Next Article in Journal
Nucleolin Overexpression Predicts Patient Prognosis While Providing a Framework for Targeted Therapeutic Intervention in Lung Cancer
Next Article in Special Issue
Comparative Study of the Role of Interepithelial Mucosal Mast Cells in the Context of Intestinal Adenoma-Carcinoma Progression
Previous Article in Journal
An Invasive Disease Event-Free Survival Analysis to Investigate Ki67 Role with Respect to Breast Cancer Patients’ Age: A Retrospective Cohort Study
Previous Article in Special Issue
The Cellular Tumor Immune Microenvironment of Childhood Solid Cancers: Informing More Effective Immunotherapies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cancers
Volume 14
Issue 9
10.3390/cancers14092216
cancers-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

Characterizing the Inflammatory Microenvironment in K14-HPV16 Transgenic Mice: Mast Cell Infiltration and MicroRNA Expression

by
Alexandra C. Costa
1,2,3,†,
Joana M. O. Santos
1,2,†,
Beatriz Medeiros-Fonseca
4,
Paula A. Oliveira
4,
Margarida M. S. M. Bastos
5,6,
Haissa O. Brito
7,
Rui M. Gil da Costa
1,4,5,6,7 and
Rui Medeiros
1,2,3,8,9,*
1
Molecular Oncology and Viral Pathology Group, Research Center of IPO Porto (CI-IPOP)/RISE@CI-IPOP (Health Research Network), Portuguese Oncology Institute of Porto (IPO Porto), Porto Comprehensive Cancer Center (Porto.CCC), 4200-072 Porto, Portugal
2
Faculty of Medicine of the University of Porto (FMUP), 4200-319 Porto, Portugal
3
Research Department of the Portuguese League against Cancer—Regional Nucleus of the North (Liga Portuguesa Contra o Cancro—Núcleo Regional do Norte), 4200-177 Porto, Portugal
4
Centre for the Research and Technology of Agro-Environmental and Biological Sciences (CITAB), Inov4Agro, University of Trás-os-Montes e Alto Douro (UTAD), Quinta de Prados, 5000-801 Vila Real, Portugal
5
LEPABE—Laboratory for Process Engineering, Environment, Biotechnology and Energy, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, 4200-465 Porto, Portugal
6
ALiCE—Associate Laboratory in Chemical Engineering, Faculty of Engineering, University of Porto, Rua Dr. Roberto Frias, 4200-465 Porto, Portugal
7
Postgraduate Programme in Adult Health (PPGSAD), Department of Morphology, Federal University of Maranhão (UFMA), and UFMA University Hospital (HUUFMA), São Luís 65080-805, Brazil
8
Virology Service, Portuguese Oncology Institute of Porto (IPO Porto), 4200-072 Porto, Portugal
9
Biomedical Research Center (CEBIMED), Faculty of Health Sciences of the Fernando Pessoa University, 4249-004 Porto, Portugal
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to this work and are considered joint first authors.
Cancers 2022, 14(9), 2216; https://doi.org/10.3390/cancers14092216
Submission received: 18 March 2022 / Revised: 12 April 2022 / Accepted: 26 April 2022 / Published: 28 April 2022
(This article belongs to the Special Issue The Role of Immune Cells in the Tumor Microenvironment)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

K14-HPV16 transgenic mice have proved to be a useful model to study the carcinogenic cascade induced by HPV16, the tumor microenvironment and also the epigenetic and genetic factors associated with this type of malignancy. The aim of our study was to evaluate the infiltration of mast cells in two cutaneous regions with different severity of the lesions and to identify potential microRNAs that may regulate mast cell infiltration in this model. We were able to confirm that increased mast cell infiltration is associated with progression of HPV-induced lesions, and that miR-223-3p and miR-125b-5p might be assisting this process via the regulation of mast cell chemotactic proteins.

Abstract

High-risk human papillomavirus (HPV) is the etiologic agent of several types of cancer. Mast cells’ role as either a driving or opposing force for cancer progression remains controversial. MicroRNAs are dysregulated in several HPV-induced cancers, and can influence mast cell biology. The aim of this study was to evaluate mast cell infiltration and to identify microRNAs potentially regulating this process. Transgenic male mice (K14-HPV16; HPV+) and matched wild-type mice (HPV−) received 7,12-Dimethylbenz[a]anthracene (DMBA) (or vehicle) over 17 weeks. Following euthanasia, chest skin and ear tissue samples were collected. Mast cell infiltration was evaluated by immunohistochemistry. MicroRNAs associated with mast cell infiltration were identified using bioinformatic tools. MicroRNA and mRNA relative expression was evaluated by RT-qPCR. Immunohistochemistry showed increased mast cell infiltration in HPV+ mice (p < 0.001). DMBA did not have any statistically significant influence on this distribution. Ear tissue of HPV+ mice showed increased mast cell infiltration (p < 0.01) when compared with chest skin samples. Additionally, reduced relative expression of miR-125b-5p (p = 0.008, 2−ΔΔCt = 2.09) and miR-223-3p (p = 0.013, 2−ΔΔCt = 4.42) seems to be associated with mast cell infiltration and increased expression of target gene Cxcl10. These results indicate that HPV16 may increase mast cell infiltration by down-regulating miR-223-3p and miR-125b-5p.

1. Introduction

High-risk human papillomavirus (HPV) is responsible for cancers in several anatomic regions, including the cervix, vulva, vagina, anus, penis, and oropharynx [1]. HPV16 and HPV18 are the high-risk type most commonly found in most of those cancers. In high-risk HPV, the E6, E7 and E5 proteins play an oncogenic role due to their capacity to form specific complexes with tumor suppressor proteins [2]. Importantly, these oncoproteins have the capacity to induce pro-inflammatory signaling pathways [3,4]. Such pro-inflammatory stimuli, along with chronic inflammation resulting from host antiviral response to persistent HPV infection, are essential for tumor development [3,4]. The inflammatory tumor microenvironment is composed by pre-malignant and malignant cells, stromal cells, and innate and adaptative immune cells [5], which interact with each other, modulating cancer development [5]. Mast cells influence other cells, both through soluble mediators and cell-to-cell interactions [6]. The role of mast cells in cancer progression is ambiguous, since they can promote an anti-tumor response or contribute to tumor growth [7,8]. The pro-tumoral functions of mast cells include the promotion of metastasis, the release of pro-angiogenic factors and the participation in immunosuppression and cancer-related inflammation [9]. Additionally, mast cell accumulation has been associated with poor prognosis and reduced survival in several solid tumors [10]. The role of mast cells in high-risk HPV-induced carcinogenesis and the mechanisms that regulate their accumulation in the tumor microenvironment are poorly defined.
MicroRNAs (miRNAs/miRs) are small and conserved RNA molecules that play an important role in the regulation of gene expression [11]. The oncoproteins of high-risk HPV are able to alter the expression of cellular miRNAs, contributing to malignant transformation [11]. Further, miRNAs are key players in modulating the tumor microenvironment, including the infiltration and activity of immune cells [12,13]. However, there is not much information regarding the effects of miRNAs in mast cells in the context of cancer development. The HPV oncoproteins may deregulate miRNAs that control mast cell migration, modulating their infiltration into the tumor microenvironment. Therefore, in this study, we sought to evaluate the infiltration of mast cells into HPV16-driven lesions in transgenic mice carrying all of the HPV16 early genes [14,15]. Additionally, we aimed to identify potential microRNAs that may regulate mast cell infiltration in this model.

2. Materials and Methods

2.1. Animals

Generation of K14-HPV16 mice on an FVB/n background has been previously reported [16]. These transgenic mice were kindly donated by Dr. Jeffrey Arbeit and Dr. Douglas Hanahan (University of California) through the USA National Cancer Institute Mouse Repository. The animal experiments were approved by the University of Trás-os-Montes and Alto Douro Ethics Committee (approval no. 10/2013) and the Portuguese General Veterinary Directorate (0421/000/000/2014). Animals were housed and bred according to Portuguese (Decreto-Lei 113, 7 August) and European (EU Directive 2010/63/EU) legislation, under controlled temperature (23 ± 2 °C), light-dark cycle (12 h light/12 h dark) and relative humidity (50 ± 10%). Food and water were provided ad libitum.

2.2. Experimental Design and Sample Collection

The chest and ear samples and experimental design were previously used by Peixoto da Silva et. al., 2020 [17]. Briefly, 34 male mice were randomly allocated in four groups: Group 1 (n = 9, HPV− mice, without 7,12-Dimethylbenz[a]anthracene (DMBA)), Group 2 (n = 5, HPV− mice, with DMBA), Group 3 (n = 10, HPV+ mice, without DMBA) and Group 4 (n = 10, HPV+ mice, with DMBA). The DMBA (D3254, Sigma-Aldrich, Merck kGAa, Darmstadt, Germany) was dissolved in dimethyl sulfoxide (DMSO) (CARLO ERBA Reagents S.A.S., Val de Reuil, France) and topically administered to the penile mucosa once a week (0.031 mg/animal/week in 4 μL of DMSO), starting at 9–11 weeks during 17 weeks of the experiment. All of the mice were sacrificed at 31–33 weeks of age under ketamine (CLORKETAM 1000, injectable solution, Vétoquinol, Barcarena, Portugal) and xylazine (Rompun® 2%, Bayer Healthcare S.A., Kiel, Germany) anesthesia, by intracardiac puncture and exsanguination, as indicated by the Federation for Laboratory Animal Science Associations (FELASA). The chest skin and ear tissue samples were weighed and collected into TripleXtractor reagent (Grisp®, Porto, Portugal), macerated and kept at −80 °C until further use. Matched chest skin and ear samples were also collected for histological analysis (histological results previously published by [17]).

2.3. Immunohistochemistry

Formalin-fixed, paraffin-embedded (FFPE) tissue sections were fixed, deparaffinized in xylene, rehydrated in a graded series of alcohol, and subjected to antigen retrieval through citrate buffer and microwave heating for 15 min and then cooled for 20 min at room temperature. The slides were covered for 5 min with peroxidase blocking reagent (Novolink Polymer Detection System, Novocastra, Leica Biosystems, Germany), followed by washing with PBS 1×-Tween 0.05% (PBS-T), and protein block (Novolink Polymer Detection System, Novocastra) was used for 5 min. Individual slides were incubated for 1 h 15 min with anti-TPSAB-1 (SC68-07; ThermoFisher) antibody diluted in PBS 1× (1:200 dilution), at room temperature. The slides were then washed with PBS-T and incubated with Novolink Polymer (30 min) and washed with PBS-T again. The color reaction was developed in DAB Chromogen (Novolink Polymer Detection System, Novocastra) according to the manufacturer’s instructions. The sections were then counterstained with hematoxylin, dehydrated and mounted. Positive controls consisted of sections of lymph node tissues, while for negative controls the primary antibody was omitted (Figure S1). The staining was assessed and validated by an experienced pathologist (RMGC). To evaluate mast cell infiltration, the number of immunostained cells was determined in at least five high magnification fields (500×) for each animal, randomly chosen from the superficial dermis, representing the dermo-epidermal interface, and counting was performed in ImageJ. The results were presented as the mean and standard deviation of the animals in each group.

2.4. MicroRNA Selection

To retrieve the proteins associated with mast cell infiltration we created a protein–protein interaction (PPI) network using Cytoscape 3.7.2. The Search Tool for the Retrieval of Interacting Genes/Proteins (STRINGapp 1.5.1), installed on Cytoscape, was used to search and visualize the proteins associated with mast cell infiltration and mast cell chemotaxis through Pubmed Query. A confidence (score) cutoff of 0.4 and a number of proteins of 100 were selected. The PPI network was further clustered using Clustermaker 2 v.1.3.1, to identify densely connected regions and clusters of proteins that were crucial in the PPI network. The CytoHubba (version 0.1) plugin Cytoscape was also utilized to identify the key proteins within the higher cluster, according to maximal clique centrality (MCC). Finally, for the top 10 proteins, we identified microRNAs that regulate their respective messenger RNAs, using miRNA databases (TargetScan, miRDB, miRmap, miRwalk, miRTarBase and TarBase) [18,19,20,21,22,23]. The four microRNAs with higher scores in those databases and related to the majority of the genes in our network were selected.

2.5. Total RNA Isolation

Total RNA extraction from chest skin and ear samples was performed using TripleXtractor reagent (Grisp®, Portugal) followed by a chloroform solution (EMSURE®, Merck kGAa, Darmstadt, Germany). Total RNA was then purified using a commercial kit, GRS total RNA kit (Grisp®, Portugal), according to the manufacturer’s instructions. DNAse treatment was included for all of the samples. RNA concentration and purity were measured using the NanoDrop Lite spectrophotometer (Thermo Scientific®).

2.6. cDNA Synthesis

The microRNA samples (50 ng) were used as templates for cDNA synthesis using the TaqMan® MicroRNA Reverse Transcription Kit (Applied Biosystems®) and sequence-specific stem-loop reverse transcription primers for miR-223-3p, miR-466l-3p, miR-466k, miR-125b-5p and snoRNA-202. The conversion was performed in a Biometra® Personal Cycler (Biocompare, USA) with the following conditions: 30 min at 16 °C, 60 min at 42 °C and 10 min at 85 °C. Total RNA samples (300 ng) were also used as templates for cDNA synthesis using the High-Capacity cDNA Reverse Transcription Kit (Applied BiosystemTM). The conversion was performed in a Biometra® Personal Cycler (Biocompare, USA) with the following conditions: 10 min at 25 °C, 120 min at 37 °C, and 5 min at 85 °C. Negative controls lacking RNA were also included in all reactions.

2.7. Relative Quantification of MicroRNAs

Each reaction was performed with 5 µL of 2× TaqMan® Fast Advanced Master Mix (Applied Biosystems®) with 0.5 µL of 20× TaqMan® MicroRNA Expression Assays (miR-466l-3p: 002804; miR-466k: 240990_mat; miR-223-3p: 002295; miR-125b-5p: 000449; snoRNA-202: 001232; Applied Biosystems®), 3 µL of nuclease-free water and 1.5 µL of cDNA sample, making a total volume of 10 µL. The quantification was performed in duplicate, and CT standard deviation values superior to 0.5 were excluded. Negative controls lacking cDNA were also included in all reactions. All target and endogenous controls for each sample were amplified in the same plate. The thermal cycling conditions for all assays were the following: 20 s at 95 °C followed by 45 cycles of 1 s at 95 °C and 20 s at 60 °C. The same baseline and threshold were set for each plate using the analysis software for qPCR from the Thermo Fisher Connect platform (Thermo Fisher Scientific, Waltham, MA, USA), in order to generate threshold cycle (Ct) values for all of the miRs/SnoR in each sample. Small nucleolar RNA 202 (snoR-202) was previously tested by our group in this mouse model, was the one that showed the lowest standard deviation values in chest and ear, and therefore, was used as endogenous control [14,15].

2.8. Relative Quantification of MRNAs

Each reaction was performed with 5 µL of 2× TaqMan® Fast Advanced Master Mix (Applied Biosystems®) with 0.5 µL of 20× TaqMan® Gene Expression Assays (Ccr2: Mm00438270_m1; Cxcl10: Mm00445235_m1; β-actin: Mm01205647_g1; Applied Biosystems®), 3.5 µL of nuclease-free water and 1 µL of cDNA sample, making a total volume of 10 µL. The quantification was performed in duplicate, and CT standard deviation values superior to 0.5 were excluded. Negative controls lacking cDNA were also included in all reactions. All target and endogenous controls for each sample were amplified in the same plate. The thermal cycling conditions for all assays were the following: 20 s at 95 °C followed by 45 cycles of 1 s at 95 °C and 20 s at 60 °C. The same baseline and threshold were set for each plate using the analysis software for qPCR from the Thermo Fisher Connect platform (Thermo Fisher Scientific, Waltham, MA, USA), in order to generate Ct values for all of the genes in each sample. β-actin was previously tested by our group in chest skin and ear tissue samples of this mice model, was the one that showed the lowest standard deviation values, and therefore, was used as endogenous control.

2.9. Statistical Analysis

Statistical analysis was performed using IBM SPSS Statistics for Windows (Version 27.0). Immune cell infiltration was evaluated using the Mann–Whitney U test and the Kruskal–Wallis test. Additionally, to analyze mast cell distribution in the different histological lesions, a Fisher–Freeman–Halton test was performed. In order to understand which microRNAs were influencing the various degrees of mast cell infiltration, we used the Kruskal–Wallis test. Finally, the presence of statistical differences in microRNA expression was evaluated using the Livak method along with Mann–Whitney U test. All of the graphics were constructed using GraphPad Prism 8 (GraphPad Software). The results were considered statistically significant when the p values were <0.05.

3. Results

3.1. Mast Cell Infiltration

To determine if HPV was able to modulate mast cell infiltration in both chest skin and ear tissues of the K14-HPV16 transgenic mice, we evaluated the infiltration of mast cells in FFPE tissue specimens by immunohistochemistry using a specific anti-mast cell tryptase antibody (TPSAB-1). We observed that, in general, mast cell infiltration increased substantially in the transgenic mice (HPV+) when compared with the wild type (HPV) (p < 0.001) (Figure 1 and Figure 2). The presence or absence of DMBA did not have any statistically significant influence on this distribution (p = 0.650 in wild-type mice and p = 0.568 in transgenic mice) (Figure 2).
We also evaluated mast cell infiltration in the two organs separately. The mean number of mast cells was higher in ear tissue when compared with chest skin (Figure 3). This alteration was perceptible in both wild-type (HPV−) and transgenic mice (HPV+), although it was much more evident in the HPV+ groups, both without and with DMBA (p < 0.01) (Figure 3).
Next, we wanted to understand if mast cell infiltration increased as lesions progressed from normal epidermis to epidermal hyperplasia, dysplasia and squamous cell carcinoma. Analyzing both organs together, we observed important differences between the normal epithelium and the hyperplastic epithelium (p < 0.001), as well as between the normal epithelium and the dysplastic epithelium (p < 0.001) (Figure 4a). No statistically significant difference was observed between the hyperplastic epithelium and the dysplastic epithelium (p = 0.348) (Figure 4a). The evaluation of chest skin alone showed, once again, differences between the normal epithelium when compared with both hyperplastic and dysplastic epithelium (p = 0.04 and p = 0.03, respectively) (Figure 4b). The same occurred when we analyzed only the ear tissue (p < 0.001) (Figure 4c). The single carcinoma sample could not be analyzed statistically, but mast cell infiltration was detected and mast cell counts were similar to those from dysplastic lesions (Figure 4c).
In order to determine how mast cells were distributed in the different histological lesions, we ranked mast cell averages in terciles, using the IBM SPSS software platform. A rank of 1 was assigned to a mean of 0 to 2 mast cells per camp and represented a low mast cell count; a rank of 2 was assigned to a mean of 2 to 9 mast cells per camp and represented an intermediate mast cell count; and finally, a rank of 3 was assigned to a mean of 9 to 33 mast cells per camp and represented a high mast cell count. A low mast cell count was mainly observed in normal samples, an intermediate mast cell count was associated with 4 hyperplastic samples and 8 dysplastic samples, and the highest mast cell count was found homogeneously distributed between the hyperplastic and dysplastic lesions (p < 0.001) (Table 1a). By analyzing each organ separately, we were able to confirm that more advanced lesions had an increase in the mast cell count (p = 0.002 and p < 0.001 for Table 1b,c respectively).

3.2. Identification of Proteins and MicroRNAs Potentially Involved in Mast Cell Infiltration

To identify proteins related with mast cell infiltration, we used StringApp 1.6.0 in Cytoscape 3.7.2. The top 100 proteins were retrieved with a confidence cut-off of 0.4. Next, using clusterMaker2, we applied the Markov clustering (MCL) to promote the clustering. We obtained clusters of proteins that had more interaction among themselves than with the rest of the network. Then, for the major cluster obtained, the Cytohubba 0.1 app was used to identify the top 10 hub proteins (Table 2). Having found the essential proteins associated with mast cell infiltration (Table 2), we then searched in several miRNA online databases for miRNAs that target the messenger RNA of those top 10 proteins and retrieved the four microRNAs with higher scores, namely miR-466l-3p, miR-466k, miR-223-3p and miR-125b-5p (Table 3). Next, we tried to understand if these four microRNAs could be associated with mast cell infiltration in chest skin and ear tissue of K14-HPV16 transgenic mice. Performing a Kruskal–Wallis test between the relative expression of the four microRNAs and the terciles obtained previously (low, intermediate and high mast cell count), we found that only miR-223-3p and miR-125b-5p had statistically significant results (p = 0.035 and p = 0.024, respectively), meaning that only these two microRNAs were influencing mast cell infiltration (Table 4). For miR-466l-3p and miR-466k, no statistical differences were found, so no further evaluation was performed (Table 4 and Figure 5).

3.3. MiR-223-3p Is Associated with Mast Cell Infiltration

To confirm our hypothesis that miR-223-3p could be associated with mast cell infiltration in this transgenic mouse model, we first conducted a joint analysis of both organs, and observed that there were statistically significant differences between the low mast cell count and the high mast cell count (p = 0.013; 2−ΔΔCt = 4.42) and also between the low mast cell count and the intermediate mast cell count (p = 0.040; 2−ΔΔCt = 4.85) (Figure 6a).
Since Cxcl10 is one of the essential proteins that is targeted by miR-223-3p and mast cell infiltration, we evaluated its expression and observed that there were statistically significant differences between the low mast cell count and the intermediate and high mast cell counts (p = 0.002; 2−ΔΔCt = 5.37 and p < 0.001; 2−ΔΔCt = 9.30, respectively) and also between the intermediate mast cell count and the high mast cell count (p = 0.023; 2−ΔΔCt = 1.73) (Figure 6b).
When we performed a more specific analysis and centered our attention on chest skin, we did not find any statistically significant difference when analyzing the miR-223 relative expression (Figure 6c,e). Nonetheless, when analyzing Cxcl10 expression, we observed a statistically significant difference between the low mast cell count and the intermediate mast cell count (p = 0.008; 2−ΔΔCt = 6.38) (Figure 6d). The same happened when we analyzed the ear tissue separately, where miR-223 did not have any statistically significant differences (Figure 6e), but when we focused on Cxcl10 expression, we observed statistically significant differences between the low mast cell count and high mast cell count (p = 0.009; 2−ΔΔCt = 9.08) and also between the intermediate mast cell count and high mast cell count (p = 0.019; 2−ΔΔCt = 2.81) (Figure 6f).

3.4. MiR-125b-5p Is Associated with Mast Cell Infiltration

Next, we performed the same analysis for miR-125b-5p involving both organs and observed statistically significant results between the low mast cell count and the high mast cell count (p = 0.008; 2−ΔΔCt = 2.09) (Figure 7a). When focusing on each organ separately, there was a statistically significant difference between the low mast cell count and the intermediate mast cell count (p = 0.026; 2−ΔΔCt = 2.21) in chest skin, and between the intermediate mast cell count and high mast cell count (p = 0.031; 2−ΔΔCt = 2.27) in ear tissue, suggesting that this miRNA influences the increased mast cell infiltration in both organs (Figure 6e and Figure 7c). Additionally, we also evaluated Ccr2 relative expression, since it is one of the essential proteins related to miR-125b. There were statistically significant differences in Ccr2 expression between the low mast cell count and the intermediate mast cell count (p = 0.035; 2−ΔΔCt = 1.61) and between the intermediate mast cell count and the high mast cell count (p = 0.036; 2−ΔΔCt = 1.55) when we analyzed both organs together (Figure 7b). However, when analyzing chest skin (Figure 7d) and ear (Figure 7f) separately, no statistical differences were found.

4. Discussion

4.1. Mast Cell Infiltration

High-risk HPV is considered one of the main etiological factors associated with the development of cancer [28]. The carcinogenic process associated with a persistent HPV infection develops in several locations through a multistep process, that goes from hyperplastic lesions (alteration in the cell number, shape and size, but still maintaining some differentiation), to dysplastic lesions (with cytological aberrations and loss of tissue differentiation), carcinoma in situ (CIS), and subsequently to invasive squamous cell carcinoma (SCC) [29]. There is a critical connection between the inflammatory microenvironment and cancer [30]. Numerous studies have documented this link in several human carcinomas and the influence of local immunity over patient prognosis and response to therapy [31,32,33,34]. In the early stages of tumor development, immunogenic cancer cells are recognized and eliminated by cytotoxic immune cells such as natural killer (NK) cells and CD8+ T cells (process called immunosurveillance) [35]. As the neoplastic tissue evolves into a clinically detectable tumor, there are different inflammatory cells that will impact tumor fate. These so-called tumor-infiltrating immune cells (TIICs) can be trained by tumor cells to favor tumor growth and development. An example of such TIICs are mast cells. Mast cells can secrete a wide range of bioactive molecules (that are contained inside their cytoplasmic granules) that exert both pro- and anti-tumor effects [36]. Several in vivo and in vitro studies showed that mast cells’ pro-tumor activity promotes lymphatic and blood vessel formation, tumor growth and metastasis [37].
With regard to the model that we used in this study, the integration of the early region of HPV16 DNA in these transgenic mice allows the spontaneous development of lesions that are very similar to the ones observed in humans, with analogous multistep development [38]. We were able to observe pre-malignant lesions and squamous cell carcinomas on the epidermis of the ear and chest tissue of these transgenic mice [15,38]. Our team and others also used this mouse strain to develop models for tongue base cancer, cervical cancer and penile cancer [39,40,41]. It is thought that a progression of the lesions associated with a persistent HPV infection is accompanied by increased immune cell infiltration, and that was what we wanted to verify whether this happened with mast cells [15,42]. There are previous reports from HPV16-E7 mice showing that accumulation of mast cells in the ear lesions causes local immune suppression [42]. In K14-HPV16 mice, it was shown that cutaneous mast cells release proteases and reorganize stromal architecture and hyperactivate angiogenesis [43]. However, the comparison between cutaneous regions with different susceptibility to HPV-induced carcinogenesis and the mechanisms that promote mast cell accumulation in this mouse model are poorly defined.
The K14-HPV16 mice used in this study had intraepithelial hyperplastic lesions in the chest skin, and dysplastic lesions were restricted to the ear skin [17]. DMBA induced dysplastic lesions in the chest skin and a squamous cell carcinoma in the ear skin [17]. There is evidence that suggests that differences in the histological microenvironment of chest skin and ear tissue of these transgenic mice are correlated with a differential microRNA expression in both of these organs that ultimately affects their susceptibility to cancer [14,15].
Our results were consistent with these previous studies, as we identified not only that HPV presence influenced a higher infiltration of mast cells, but also that this accumulation was more prominent in the ear tissue when compared to the chest skin (Figure 1, Figure 2 and Figure 3). Indeed, this increase in mast cells in ear tissue compared to chest skin was observed in both HPV− and HPV+ mice. Since DMBA presence did not influence the results, mast cell infiltration seems to only be promoted by HPV presence.
Since our transgenic animals (HPV+) had different histological lesions, we were also able to evaluate if the infiltration of mast cells was consistent during the entire carcinogenic process or if different lesions possessed distinct levels of mast cell infiltration. We concluded that as lesions progressed, there was also an increase in mast cell infiltration (Figure 4). This supports the hypothesis that increased mast cell counts are associated with lesion progression and the consequent development of cancer, although additional mechanistic studies are required to establish the tumorigenic contributions of mast cells in this model.

4.2. MicroRNAs and Their Putative Roles

MiRNAs possess a dual role, being able to act as tumor suppressors or oncogenes, depending on their target gene functions [44]. Some miRNAs have been described as putative biomarkers for the occurrence and development of HPV-induced cancers [45,46]. Numerous miRNAs have been reported to have a crucial role in HPV-induced carcinogenesis [47,48], including in studies by our group using K14-HPV16 transgenic mice [14,15,49,50]. MiRNAs also modulate many aspects of mast cell behavior (mast cell cycle, proliferation and maturation), being able to act as either inducers and/or suppressors of mast cell biological functions by affecting different targets [51]. Despite the importance of this emerging field of research, the available data on the influence of miRNAs in mast cell function, particularly in cancer, remain very limited.

4.2.1. MiR-223-3p Expression Analysis

MiR-223 in humans is located within the q12 locus of the X chromosome and is mainly expressed by hematopoietic cells [52]. MiR-223 remained highly conserved during evolution, which suggests that it has an important role in several physiological processes, including monocyte-macrophage differentiation, neutrophil recruitment and pro-inflammatory responses [53]. Additionally, miR-223 can also be transferred to non-myeloid cells through lipoproteins or extracellular vesicles [54]. This microRNA’s expression is also deregulated in many types of cancer, from hematological to solid malignancies [55]. Depending on the clinical context, miR-223 can act either as an oncogene (breast cancer) or tumor suppressor (leukemia and hepatocellular carcinoma) [55]. Additionally, studies suggest that downregulation of miR-223 seems to promote mast cell degranulation, while high levels of miR-223 may promote mast cell apoptosis [56,57].
In our study, miR-223 was chosen because it was associated with two proteins related to mast cell infiltration and chemotaxis, namely C-X-C motif chemokine ligand 10 (Cxcl10) and integrin subunit alpha M (Itgam). From miR-223 targets, we chose to study the expression of Cxcl10, because from the two the targets of miR-223 it was the one on top of the rank (Table 3). Cxcl10 is a chemokine produced by cancer cells that can activate mast cells and promote their recruitment to solid tumors [58]. Furthermore, this chemokine can stimulate mast cells to secrete diverse soluble factors that support invasion, proliferation and survival of tumor cells [59,60].
Our results suggest that miR-223 may interfere with mast cell infiltration in both chest skin and ear tissue of our transgenic mice, as we observed a tendency toward decreased miR-223 expression while mast cell counts increased. However, this tendency was not statistically significant when analyzing chest skin (Figure 6c) and ear (Figure 6e) separately, possibly due to the lack of statistical power due to the small sample size in some of the groups. Future studies should increase the sample size. Regarding Cxcl10, we were able to verify that an increase in Cxcl10 levels was associated with an increase in mast cell infiltration in chest skin and ear tissue lesions. Thus, we hypothesize that a decrease in miR-223 relative expression may lead to an increase in Cxcl10, which will probably increase mast cell infiltration in the tissue.

4.2.2. MiR-125b-5p Expression Analysis

MiR-125b has gained special interest in the field of cancer research because it is found aberrantly expressed in a great variety of tumors [61]. This miRNA is the human orthologue of lin-4, one of the first miRNAs to be identified in C. elegans [62]. The function of miR-125b diverges in different cancers depending on the molecular contexts and surrounding tumor microenvironment [61]. It can either act as an oncogene (e.g., nasopharyngeal carcinoma and gastric cancer) or as a tumor suppressor (e.g., breast cancer and colorectal cancer) [63]. Mature miR-125b is generated from two genes (miR-125b-1 and miR-125b-2), both situated in fragile sites commonly deleted in these types of cancer and that imply its loss of function [64].
In general, miR-125b is involved in different cellular processes such as inflammation, cell proliferation and cell cycle regulation. MiR-125b is capable of directly targeting p53, which is known for maintaining genome stability and playing a central role in apoptosis regulation [65]. Furthermore, miR-125b has been shown to interact with multiple mRNAs, including apoptosis regulators such as Bak-1, Bcl-2 and Bcl-w [66,67].
MiR-125b has a central role in productive HPV infection through the regulation of papillomavirus ORF L2 [68]. L2 is associated with the import of HPV DNA to the nucleus, and its inactivation by miR-125b results in loss of viral infectivity [69,70]. Shortly after HPV infection, there is an increased expression of miR-125b that can be explained by interaction with the viral oncoproteins, and that promotes the immune response and inflammation process [71]. As lesions progress, there is a significant decrease in the relative expression of miR-125b [72]. Regarding our bioinformatic analysis, miR-125b was also associated with two proteins related to mast cell infiltration and mast cell chemotaxis, namely C-C motif chemokine receptor 2 (Ccr2) and tumor necrosis factor (Tnf) (Table 3).
The results showed a decrease in miR-125b relative expression as the number of mast cells infiltrated in the tumor increased, which we hypothesized could be mediated by increased Ccr2 expression levels. Ccr2 is a chemokine that is thought to promote mast cell degranulation and to stimulate an increase in carcinogenesis-associated inflammation mediated by mast cells [73,74]. However, Ccr2 was not significantly regulated in relation to mast cell infiltration, suggesting that miR-125b may be regulating mast cell infiltration via other downstream effectors. In this context, the expression of Tnf, which was also identified as a major regulated protein in our network, should be explored in future studies.

5. Conclusions

With this work and its findings, we conclude that HPV16 promotes the progressive accumulation of mast cells in intraepithelial lesions, potentially promoting carcinogenesis. Additionally, we also came to the conclusion that miR-223-3p and miR-125b-5p may be orchestrating pathways that regulate mast cell infiltration in chest skin and ear tissue. There was increased mast cell infiltration in the microenvironment of ear tissue compared with chest skin in both HPV− and HPV+ mice. This finding may be a possible explanation for the differences in the severity of the lesions in these two cutaneous sites, observed in HPV+ mice.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cancers14092216/s1, Figure S1: Histological analysis of lymph node samples, DAB-hematoxylin. (a). Lymph node with mast cell infiltration observed in mice, 100×. (b). Same lymph node from figure (a). with high magnification 500×. (c). Lymph node without primary antibody, 100×. (d). Same lymph node from figure c. with high magnification 500×.

Author Contributions

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

Funding

This study was funded by the Portuguese League Against Cancer—Regional Nucleus of the North (Liga Portuguesa Contra o Cancro—Núcleo Regional do Norte), by the Research Center of the Portuguese Oncology Institute of Porto (project no. PI86-CI-IPOP-66-2017), by LA/P/0045/2020 (ALiCE), UIDB/00511/2020 and UIDP/00511/2020 (LEPABE), funded by national funds through FCT/MCTES (PIDDAC), and by National Funds from FCT—Portuguese Foundation for Science and Technology, under the project UID/AGR/04033/2020. 2SMART (NORTE-01-0145-FEDER-000054), supported by Norte Portugal Regional Operational Programme (NORTE 2020), under the PORTUGAL 2020 Partnership Agreement, through the European Regional Development Fund (ERDF). Joana M.O. Santos is a Ph.D. scholarship holder (SFRH/BD/135871/2018) supported by Fundação para a Ciência e Tecnologia (FCT), co-financed by European Social Funds (FSE) and national funds of MCTES.

Institutional Review Board Statement

The animal experiments were approved by the University of Trás-os-Montes and Alto Douro Ethics Committee (approval no. 10/2013) and the Portuguese General Veterinary Directorate (0421/000/000/2014).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kombe, A.J.K.; Li, B.; Zahid, A.; Mengist, H.M.; Bounda, G.A.; Zhou, Y.; Jin, T. Epidemiology and Burden of Human Papillomavirus and Related Diseases, Molecular Pathogenesis, and Vaccine Evaluation. Front. Public Health 2020, 8, 552028. [Google Scholar] [CrossRef] [PubMed]
  2. Estevao, D.; Costa, N.R.; da Costa, R.M.G.; Medeiros, R. Hallmarks of HPV carcinogenesis: The role of E6, E7 and E5 oncoproteins in cellular malignancy. Biochim. Biophys. Acta Gene Regul. Mech. 2019, 1862, 153–162. [Google Scholar] [CrossRef] [PubMed]
  3. Boccardo, E.; Lepique, A.P.; Villa, L.L. The role of inflammation in HPV carcinogenesis. Carcinogenesis 2010, 31, 1905–1912. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Hemmat, N.; Baghi, H.B. Association of human papillomavirus infection and inflammation in cervical cancer. Pathog. Dis. 2019, 77, ftz048. [Google Scholar] [CrossRef]
  5. de Visser, K.E.; Coussens, L.M. The inflammatory tumor microenvironment and its impact on cancer development. Infect. Inflamm. Impacts Oncog. 2006, 13, 118–137. [Google Scholar] [CrossRef]
  6. Dyduch, G.; Kaczmarczyk, K.; Okon, K. Mast cells and cancer: Enemies or allies? Pol. J. Pathol. 2012, 63, 1–7. [Google Scholar]
  7. Marichal, T.; Tsai, M.; Galli, S.J. Mast cells: Potential positive and negative roles in tumor biology. Cancer Immunol. Res. 2013, 1, 269–279. [Google Scholar] [CrossRef] [Green Version]
  8. Derakhshani, A.; Vahidian, F.; Alihasanzadeh, M.; Mokhtarzadeh, A.; Nezhad, P.L.; Baradaran, B. Mast cells: A double-edged sword in cancer. Immunol. Lett. 2019, 209, 28–35. [Google Scholar] [CrossRef]
  9. Maciel, T.T.; Moura, I.C.; Hermine, O. The role of mast cells in cancers. F1000Prime Rep. 2015, 7, 9. [Google Scholar] [CrossRef]
  10. Costa, A.C.; Santos, J.M.O.; da Costa, R.M.G.; Medeiros, R. Impact of immune cells on the hallmarks of cancer: A literature review. Crit. Rev. Oncol./Hematol. 2021, 168, 103541. [Google Scholar] [CrossRef]
  11. Santos, J.M.O.; da Costa, R.M.G.; Medeiros, R. Dysregulation of cellular microRNAs by human oncogenic viruses—Implications for tumorigenesis. Biochim. Biophys. Acta Gene Regul. Mech. 2018, 1861, 95–105. [Google Scholar] [CrossRef] [PubMed]
  12. Kuninty, P.R.; Schnittert, J.; Storm, G.; Prakash, J. MicroRNA Targeting to Modulate Tumor Microenvironment. Front. Oncol. 2016, 6, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Xing, Y.; Wang, Z.; Lu, Z.; Xia, J.; Xie, Z.; Jiao, M.; Liu, R.; Chu, Y. MicroRNAs: Immune modulators in cancer immunotherapy. Immunother. Adv. 2021, 1, ltab006. [Google Scholar] [CrossRef]
  14. Paiva, I.; da Costa, R.M.G.; Ribeiro, J.; Sousa, H.; Bastos, M.; Faustino-Rocha, A.; Lopes, C.; Oliveira, P.A.; Medeiros, R. A role for microRNA-155 expression in microenvironment associated to HPV-induced carcinogenesis in K14-HPV16 transgenic mice. PLoS ONE 2015, 10, e0116868. [Google Scholar] [CrossRef]
  15. Paiva, I.; da Costa, R.M.G.; Ribeiro, J.; Sousa, H.; Bastos, M.M.; Faustino-Rocha, A.; Lopes, C.; Oliveira, P.A.; Medeiros, R. MicroRNA-21 expression and susceptibility to HPV-induced carcinogenesis—Role of microenvironment in K14-HPV16 mice model. Life Sci. 2015, 128, 8–14. [Google Scholar] [CrossRef] [Green Version]
  16. Arbeit, J.M.; Münger, K.; Howley, P.M.; Hanahan, D. Progressive squamous epithelial neoplasia in K14-human papillomavirus type 16 transgenic mice. J. Virol. 1994, 68, 4358–4368. [Google Scholar] [CrossRef] [Green Version]
  17. da Silva, S.P.; Santos, J.M.O.; Mestre, V.F.; Medeiros-Fonseca, B.; Oliveira, P.A.; Bastos, M.M.S.M.; da Costa, R.M.G.; Medeiros, R. Human Papillomavirus 16-Transgenic Mice as a Model to Study Cancer-Associated Cachexia. Int. J. Mol. Sci. 2020, 21, 5020. [Google Scholar] [CrossRef]
  18. McGeary, S.E.; Lin, K.S.; Shi, C.Y.; Pham, T.M.; Bisaria, N.; Kelley, G.M.; Bartel, D.P. The biochemical basis of microRNA targeting efficacy. Science 2019, 366, eaav1741. [Google Scholar] [CrossRef]
  19. Chen, Y.; Wang, X. miRDB: An online database for prediction of functional microRNA targets. Nucleic Acids Res. 2020, 48, D127–D131. [Google Scholar] [CrossRef] [Green Version]
  20. Vejnar, C.E.; Zdobnov, E.M. miRmap: Comprehensive prediction of microRNA target repression strength. Nucleic Acids Res. 2012, 40, 11673–11683. [Google Scholar] [CrossRef]
  21. Sticht, C.; De La Torre, C.; Parveen, A.; Gretz, N. miRWalk: An online resource for prediction of microRNA binding sites. PLoS ONE 2018, 13, e0206239. [Google Scholar] [CrossRef] [PubMed]
  22. Huang, H.Y.; Lin, Y.C.; Li, J.; Huang, K.Y.; Shrestha, S.; Hong, H.C.; Tang, Y.; Chen, Y.G.; Jin, C.N.; Yu, Y.; et al. miRTarBase 2020: Updates to the experimentally validated microRNA-target interaction database. Nucleic Acids Res. 2020, 48, D148–D154. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Karagkouni, D.; Paraskevopoulou, M.D.; Chatzopoulos, S.; Vlachos, I.S.; Tastsoglou, S.; Kanellos, I.; Papadimitriou, D.; Kavakiotis, I.; Maniou, S.; Skoufos, G.; et al. DIANA-TarBase v8: A decade-long collection of experimentally supported miRNA-gene interactions. Nucleic Acids Res. 2018, 46, D239–D245. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Friedman, R.C.; Farh, K.K.; Burge, C.B.; Bartel, D.P. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009, 19, 92–105. [Google Scholar] [CrossRef] [Green Version]
  25. Vejnar, C.E.; Blum, M.; Zdobnov, E.M. miRmap web: Comprehensive microRNA target prediction online. Nucleic Acids Res. 2013, 41, W165–W168. [Google Scholar] [CrossRef]
  26. Dweep, H.; Gretz, N.; Sticht, C. miRWalk database for miRNA-target interactions. Methods Mol. Biol. 2014, 1182, 289–305. [Google Scholar] [CrossRef]
  27. Hsu, S.D.; Lin, F.M.; Wu, W.Y.; Liang, C.; Huang, W.C.; Chan, W.L.; Tsai, W.T.; Chen, G.Z.; Lee, C.J.; Chiu, C.M.; et al. miRTarBase: A database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011, 39, D163–D169. [Google Scholar] [CrossRef] [Green Version]
  28. de Martel, C.; Ferlay, J.; Franceschi, S.; Vignat, J.; Bray, F.; Forman, D.; Plummer, M. Global burden of cancers attributable to infections in 2008: A review and synthetic analysis. Lancet Oncol. 2012, 13, 607–615. [Google Scholar] [CrossRef]
  29. Shanmugasundaram, S.; You, J. Targeting Persistent Human Papillomavirus Infection. Viruses 2017, 9, 229. [Google Scholar] [CrossRef] [Green Version]
  30. Schreiber, R.D.; Old, L.J.; Smyth, M.J. Cancer immunoediting: Integrating immunity’s roles in cancer suppression and promotion. Science 2011, 331, 1565–1570. [Google Scholar] [CrossRef] [Green Version]
  31. Jochems, C.; Schlom, J. Tumor-infiltrating immune cells and prognosis: The potential link between conventional cancer therapy and immunity. Exp. Biol. Med. 2011, 236, 567–579. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Angell, H.; Galon, J. From the immune contexture to the Immunoscore: The role of prognostic and predictive immune markers in cancer. Curr. Opin. Immunol. 2013, 25, 261–267. [Google Scholar] [CrossRef] [PubMed]
  33. Binnewies, M.; Roberts, E.W.; Kersten, K.; Chan, V.; Fearon, D.F.; Merad, M.; Coussens, L.M.; Gabrilovich, D.I.; Ostrand-Rosenberg, S.; Hedrick, C.C.; et al. Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat. Med. 2018, 24, 541–550. [Google Scholar] [CrossRef] [PubMed]
  34. Fridman, W.H.; Zitvogel, L.; Sautès-Fridman, C.; Kroemer, G. The immune contexture in cancer prognosis and treatment. Nat. Rev. Clin. Oncol. 2017, 14, 717–734. [Google Scholar] [CrossRef]
  35. Gonzalez, H.; Hagerling, C.; Werb, Z. Roles of the immune system in cancer: From tumor initiation to metastatic progression. Genes Dev. 2018, 32, 1267–1284. [Google Scholar] [CrossRef] [Green Version]
  36. Salemme, V.; Centonze, G.; Cavallo, F.; Defilippi, P.; Conti, L. The Crosstalk between Tumor Cells and the Immune Microenvironment in Breast Cancer: Implications for Immunotherapy. Front. Oncol. 2021, 11, 610303. [Google Scholar] [CrossRef]
  37. Aponte-López, A.; Fuentes-Pananá, E.M.; Cortes-Muñoz, D.; Muñoz-Cruz, S. Mast Cell, the Neglected Member of the Tumor Microenvironment: Role in Breast Cancer. J. Immunol. Res. 2018, 2018, 2584243. [Google Scholar] [CrossRef] [Green Version]
  38. Coussens, L.M.; Hanahan, D.; Arbeit, J.M. Genetic predisposition and parameters of malignant progression in K14-HPV16 transgenic mice. Am. J. Pathol. 1996, 149, 1899–1917. [Google Scholar]
  39. Mestre, V.F.; Medeiros-Fonseca, B.; Estêvão, D.; Casaca, F.; Silva, S.; Félix, A.; Silva, F.; Colaço, B.; Seixas, F.; Bastos, M.M.; et al. HPV16 is sufficient to induce squamous cell carcinoma specifically in the tongue base in transgenic mice. J. Pathol. 2020, 251, 4–11. [Google Scholar] [CrossRef]
  40. Medeiros-Fonseca, B.; Mestre, V.F.; Estêvão, D.; Sánchez, D.F.; Cañete-Portillo, S.; Fernández-Nestosa, M.J.; Casaca, F.; Silva, S.; Brito, H.; Félix, A.; et al. HPV16 induces penile intraepithelial neoplasia and squamous cell carcinoma in transgenic mice: First mouse model for HPV-related penile cancer. J. Pathol. 2020, 251, 411–419. [Google Scholar] [CrossRef]
  41. Arbeit, J.M.; Howley, P.M.; Hanahan, D. Chronic estrogen-induced cervical and vaginal squamous carcinogenesis in human papillomavirus type 16 transgenic mice. Proc. Natl. Acad. Sci. USA 1996, 93, 2930–2935. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Bergot, A.S.; Ford, N.; Leggatt, G.R.; Wells, J.W.; Frazer, I.H.; Grimbaldeston, M.A. HPV16-E7 expression in squamous epithelium creates a local immune suppressive environment via CCL2- and CCL5-mediated recruitment of mast cells. PLoS Pathog. 2014, 10, e1004466. [Google Scholar] [CrossRef]
  43. Coussens, L.M.; Raymond, W.W.; Bergers, G.; Laig-Webster, M.; Behrendtsen, O.; Werb, Z.; Caughey, G.H.; Hanahan, D. Inflammatory mast cells up-regulate angiogenesis during squamous epithelial carcinogenesis. Genes Dev. 1999, 13, 1382–1397. [Google Scholar] [CrossRef] [PubMed]
  44. Chen, P.-S.; Su, J.-L.; Hung, M.-C. Dysregulation of MicroRNAs in cancer. J. Biomed. Sci. 2012, 19, 90. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Vojtechova, Z.; Tachezy, R. The Role of miRNAs in Virus-Mediated Oncogenesis. Int. J. Mol. Sci. 2018, 19, 1217. [Google Scholar] [CrossRef] [Green Version]
  46. Lajer, C.B.; Garnæs, E.; Friis-Hansen, L.; Norrild, B.; Therkildsen, M.H.; Glud, M.; Rossing, M.; Lajer, H.; Svane, D.; Skotte, L.; et al. The role of miRNAs in human papilloma virus (HPV)-associated cancers: Bridging between HPV-related head and neck cancer and cervical cancer. Br. J. Cancer 2012, 106, 1526–1534. [Google Scholar] [CrossRef]
  47. Pardini, B.; De Maria, D.; Francavilla, A.; Di Gaetano, C.; Ronco, G.; Naccarati, A. MicroRNAs as markers of progression in cervical cancer: A systematic review. BMC Cancer 2018, 18, 696. [Google Scholar] [CrossRef]
  48. Nunvar, J.; Pagacova, L.; Vojtechova, Z.; Azevedo, N.T.D.; Smahelova, J.; Salakova, M.; Tachezy, R. Lack of Conserved miRNA Deregulation in HPV-Induced Squamous Cell Carcinomas. Biomolecules 2021, 11, 764. [Google Scholar] [CrossRef]
  49. Araújo, R.; Santos, J.M.O.; Fernandes, M.; Dias, F.; Sousa, H.; Ribeiro, J.; Bastos, M.; Oliveira, P.A.; Carmo, D.; Casaca, F.; et al. Expression profile of microRNA-146a along HPV-induced multistep carcinogenesis: A study in HPV16 transgenic mice. J. Cancer Res. Clin. Oncol. 2018, 144, 241–248. [Google Scholar] [CrossRef]
  50. Santos, J.M.O.; Fernandes, M.; Araujo, R.; Sousa, H.; Ribeiro, J.; Bastos, M.; Oliveira, P.A.; Carmo, D.; Casaca, F.; Silva, S.; et al. Dysregulated expression of microRNA-150 in human papillomavirus-induced lesions of K14-HPV16 transgenic mice. Life Sci. 2017, 175, 31–36. [Google Scholar] [CrossRef]
  51. Montagner, S.; Orlandi, E.M.; Merante, S.; Monticelli, S. The role of miRNAs in mast cells and other innate immune cells. Immunol. Rev. 2013, 253, 12–24. [Google Scholar] [CrossRef] [PubMed]
  52. Ye, D.; Zhang, T.; Lou, G.; Liu, Y. Role of miR-223 in the pathophysiology of liver diseases. Exp. Mol. Med. 2018, 50, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Wang, Y.; Dakhlallah, D.; Moldovan, L.; Anderson, T.; Ezzie, M.; Nana-Sinkam, S.P.; Mohler, P.J.; Trgovich, J.; Cook, C.; Piper, M.G.; et al. Chapter 42—Circulating MicroRNAs as Biomarkers. In MicroRNA in Regenerative Medicine; Sen, C.K., Ed.; Academic Press: Oxford, UK, 2015; pp. 1093–1125. [Google Scholar]
  54. Roffel, M.P.; Bracke, K.R.; Heijink, I.H.; Maes, T. miR-223: A Key Regulator in the Innate Immune Response in Asthma and COPD. Front. Med. 2020, 7, 196. [Google Scholar] [CrossRef] [PubMed]
  55. Favero, A.; Segatto, I.; Perin, T.; Belletti, B. The many facets of miR-223 in cancer: Oncosuppressor, oncogenic driver, therapeutic target, and biomarker of response. Wiley Interdiscip. Rev. RNA 2021, 12, e1659. [Google Scholar] [CrossRef]
  56. Wang, Q.; Zhao, D.Y.; Xu, H.; Zhou, H.; Yang, Q.Y.; Liu, F.; Zhou, G.P. Down-regulation of microRNA-223 promotes degranulation via the PI3K/Akt pathway by targeting IGF-1R in mast cells. PLoS ONE 2015, 10, e0123575. [Google Scholar] [CrossRef]
  57. Gao, H.; Deng, H.; Xu, H.; Yang, Q.; Zhou, Y.; Zhang, J.; Zhao, D.; Liu, F. MicroRNA-223 promotes mast cell apoptosis by targeting the insulin-like growth factor 1 receptor. Exp. Ther. Med. 2016, 11, 2171–2176. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Oldford, S.A.; Marshall, J.S. Mast cells as targets for immunotherapy of solid tumors. Mol. Immunol. 2015, 63, 113–124. [Google Scholar] [CrossRef]
  59. Visciano, C.; Prevete, N.; Liotti, F.; Marone, G. Tumor-Associated Mast Cells in Thyroid Cancer. Int. J. Endocrinol. 2015, 2015, 705169. [Google Scholar] [CrossRef] [Green Version]
  60. Varricchi, G.; Galdiero, M.R.; Loffredo, S.; Marone, G.; Iannone, R.; Marone, G.; Granata, F. Are Mast Cells MASTers in Cancer? Front. Immunol. 2017, 8, 424. [Google Scholar] [CrossRef] [Green Version]
  61. Banzhaf-Strathmann, J.; Edbauer, D. Good guy or bad guy: The opposing roles of microRNA 125b in cancer. Cell Commun. Signal. 2014, 12, 30. [Google Scholar] [CrossRef] [Green Version]
  62. Peng, B.; Theng, P.Y.; Le, M.T.N. Essential functions of miR-125b in cancer. Cell Prolif. 2021, 54, e12913. [Google Scholar] [CrossRef] [PubMed]
  63. Wang, Y.; Zeng, G.; Jiang, Y. The Emerging Roles of miR-125b in Cancers. Cancer Manag. Res. 2020, 12, 1079–1088. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Calin, G.A.; Sevignani, C.; Dumitru, C.D.; Hyslop, T.; Noch, E.; Yendamuri, S.; Shimizu, M.; Rattan, S.; Bullrich, F.; Negrini, M.; et al. Human microRNA genes are frequently located at fragile sites and genomic regions involved in cancers. Proc. Natl. Acad. Sci. USA 2004, 101, 2999–3004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Le, M.T.; Teh, C.; Shyh-Chang, N.; Xie, H.; Zhou, B.; Korzh, V.; Lodish, H.F.; Lim, B. MicroRNA-125b is a novel negative regulator of p53. Genes Dev. 2009, 23, 862–876. [Google Scholar] [CrossRef] [Green Version]
  66. Gong, J.; Zhang, J.P.; Li, B.; Zeng, C.; You, K.; Chen, M.X.; Yuan, Y.; Zhuang, S.M. MicroRNA-125b promotes apoptosis by regulating the expression of Mcl-1, Bcl-w and IL-6R. Oncogene 2013, 32, 3071–3079. [Google Scholar] [CrossRef] [Green Version]
  67. Le, M.T.; Shyh-Chang, N.; Khaw, S.L.; Chin, L.; Teh, C.; Tay, J.; O’Day, E.; Korzh, V.; Yang, H.; Lal, A.; et al. Conserved regulation of p53 network dosage by microRNA-125b occurs through evolving miRNA-target gene pairs. PLoS Genet. 2011, 7, e1002242. [Google Scholar] [CrossRef] [Green Version]
  68. Nuovo, G.J.; Wu, X.; Volinia, S.; Yan, F.; di Leva, G.; Chin, N.; Nicol, A.F.; Jiang, J.; Otterson, G.; Schmittgen, T.D.; et al. Strong inverse correlation between microRNA-125b and human papillomavirus DNA in productive infection. Diagn. Mol. Pathol. Am. J. Surg. Pathol. Part B 2010, 19, 135–143. [Google Scholar] [CrossRef]
  69. Florin, L.; Becker, K.A.; Lambert, C.; Nowak, T.; Sapp, C.; Strand, D.; Streeck, R.E.; Sapp, M. Identification of a dynein interacting domain in the papillomavirus minor capsid protein l2. J. Virol. 2006, 80, 6691–6696. [Google Scholar] [CrossRef] [Green Version]
  70. Zhou, J.; Sun, X.Y.; Louis, K.; Frazer, I.H. Interaction of human papillomavirus (HPV) type 16 capsid proteins with HPV DNA requires an intact L2 N-terminal sequence. J. Virol. 1994, 68, 619–625. [Google Scholar] [CrossRef] [Green Version]
  71. Nicolete, L.D.; Nicolete, R.; Haddad, R.; Azevedo, R.; Castro, F.A.; Tanaka, Y.; Takayanagui, O.M.; Covas, D.T.; Kashima, S. Upregulation of hsa-miR-125b in HTLV-1 asymptomatic carriers and HTLV-1-associated myelopathy/tropical spastic paraparesis patients. Mem. Inst. Oswaldo Cruz 2012, 107, 824–827. [Google Scholar] [CrossRef]
  72. Ribeiro, J.; Marinho-Dias, J.; Monteiro, P.; Loureiro, J.; Baldaque, I.; Medeiros, R.; Sousa, H. miR-34a and miR-125b Expression in HPV Infection and Cervical Cancer Development. BioMed Res. Int. 2015, 2015, 304584. [Google Scholar] [CrossRef] [PubMed]
  73. Collington, S.J.; Hallgren, J.; Pease, J.E.; Jones, T.G.; Rollins, B.J.; Westwick, J.; Austen, K.F.; Williams, T.J.; Gurish, M.F.; Weller, C.L. The role of the CCL2/CCR2 axis in mouse mast cell migration in vitro and in vivo. J. Immunol. 2010, 184, 6114–6123. [Google Scholar] [CrossRef] [PubMed]
  74. Huang, B.; Lei, Z.; Zhao, J.; Gong, W.; Liu, J.; Chen, Z.; Liu, Y.; Li, D.; Yuan, Y.; Zhang, G.-M.; et al. CCL2/CCR2 pathway mediates recruitment of myeloid suppressor cells to cancers. Cancer Lett. 2007, 252, 86–92. [Google Scholar] [CrossRef] [PubMed]
Cancers 14 02216 g001 550
Figure 1. Histological analysis of mice chest skin and ear tissue samples, 500×, DAB-hematoxylin, scale bar = 20 µm. (a) Chest skin with low mast cell infiltration observed in wild-type (HPV−) mice. (b) Normal ear with low mast cell infiltration observed in HPV− mice. (c) Mast cell infiltration observed in chest skin of K14-HPV16 (HPV+) transgenic mice. (d) Mast cell infiltration observed in ear tissue of HPV+ mice. (e) Chest skin of an HPV− mouse with DMBA treatment. (f) Ear tissue of an HPV− mouse with DMBA treatment. (g) Chest skin of HPV+ transgenic mouse with DMBA presence. (h) Ear tissue of HPV+ transgenic mouse with DMBA presence.
Figure 1. Histological analysis of mice chest skin and ear tissue samples, 500×, DAB-hematoxylin, scale bar = 20 µm. (a) Chest skin with low mast cell infiltration observed in wild-type (HPV−) mice. (b) Normal ear with low mast cell infiltration observed in HPV− mice. (c) Mast cell infiltration observed in chest skin of K14-HPV16 (HPV+) transgenic mice. (d) Mast cell infiltration observed in ear tissue of HPV+ mice. (e) Chest skin of an HPV− mouse with DMBA treatment. (f) Ear tissue of an HPV− mouse with DMBA treatment. (g) Chest skin of HPV+ transgenic mouse with DMBA presence. (h) Ear tissue of HPV+ transgenic mouse with DMBA presence.
Cancers 14 02216 g001
Cancers 14 02216 g002 550
Figure 2. Mast cell infiltration in chest skin and ear tissue in the different groups. Transgenic mice (HPV+) had a higher infiltration of mast cells when compared to the wild-type (HPV−) (*** p < 0.001). DMBA presence did not alter the predisposition for mast cell infiltration in any of the groups.
Figure 2. Mast cell infiltration in chest skin and ear tissue in the different groups. Transgenic mice (HPV+) had a higher infiltration of mast cells when compared to the wild-type (HPV−) (*** p < 0.001). DMBA presence did not alter the predisposition for mast cell infiltration in any of the groups.
Cancers 14 02216 g002
Cancers 14 02216 g003 550
Figure 3. Mast cell differential infiltration in each group (chest skin and ear tissue together). There was an increase in mast cell infiltration in the ear tissue when compared to the chest skin, in both wild-type and transgenic mice. HPV− had a statistically lower difference between both organs (* p < 0.05), whereas in HPV+ without and with DMBA, there was a higher statistically significant difference (** p < 0.01).
Figure 3. Mast cell differential infiltration in each group (chest skin and ear tissue together). There was an increase in mast cell infiltration in the ear tissue when compared to the chest skin, in both wild-type and transgenic mice. HPV− had a statistically lower difference between both organs (* p < 0.05), whereas in HPV+ without and with DMBA, there was a higher statistically significant difference (** p < 0.01).
Cancers 14 02216 g003
Cancers 14 02216 g004 550
Figure 4. Mast cell infiltration in lesion progression in both organs together (a), chest skin (b) and ear tissue (c). (a) Statistically significant differences were found when comparing the normal epithelium versus the hyperplastic epithelium (*** p < 0.001) and between the normal epithelium versus the dysplastic epithelium (*** p < 0.001). (b) In chest skin, comparing the normal epithelium with the hyperplastic and dysplastic epithelium, we observed an increase in mast cell infiltration (** p < 0.01). (c) In ear tissue, we also found statistically significant differences between the normal and hyperplastic epithelium (** p < 0.01) and between the normal and dysplastic epithelium (*** p < 0.001).
Figure 4. Mast cell infiltration in lesion progression in both organs together (a), chest skin (b) and ear tissue (c). (a) Statistically significant differences were found when comparing the normal epithelium versus the hyperplastic epithelium (*** p < 0.001) and between the normal epithelium versus the dysplastic epithelium (*** p < 0.001). (b) In chest skin, comparing the normal epithelium with the hyperplastic and dysplastic epithelium, we observed an increase in mast cell infiltration (** p < 0.01). (c) In ear tissue, we also found statistically significant differences between the normal and hyperplastic epithelium (** p < 0.01) and between the normal and dysplastic epithelium (*** p < 0.001).
Cancers 14 02216 g004
Cancers 14 02216 g005 550
Figure 5. MiR-466l-3p and MiR-466k relative expression in K14-HPV16 transgenic mice. (a) MiR-466l-3p relative expression in chest skin and ear tissue together. (b) MiR-466k relative expression in chest skin and ear tissue together. L: low mast cell count; I: intermediate mast cell count; H: high mast cell count.
Figure 5. MiR-466l-3p and MiR-466k relative expression in K14-HPV16 transgenic mice. (a) MiR-466l-3p relative expression in chest skin and ear tissue together. (b) MiR-466k relative expression in chest skin and ear tissue together. L: low mast cell count; I: intermediate mast cell count; H: high mast cell count.
Cancers 14 02216 g005
Cancers 14 02216 g006 550
Figure 6. MiR-223 and Cxcl10 relative expression in K14-HPV16 transgenic mice. (a) MiR-223 relative expression in chest skin and ear tissue, together. There are statistically significant differences between the low and intermediate mast cell counts (* p < 0.05) and between the low and high mast cell counts (** p < 0.01). (b) Cxcl10 relative expression in chest skin and ear tissue, together. We observed statistically significant results between the low and intermediate mast cell counts (** p < 0.01), between the low and high mast cell counts (*** p < 0.001), and between the intermediate and high mast cell counts (* p < 0.05). (c) Normalized relative expression of miR-223 in chest skin. (d) Normalized relative expression of Cxcl10 in chest skin. Statistically significant differences were found between the low and intermediate mast cell counts (** p < 0.01). (e) MiR-223 relative expression in ear tissue. (f) Cxcl10 relative expression in ear tissue. Statistically significant results were observed between the low and high mast cell counts (** p < 0.01) and between the intermediate and high mast cell counts (* p < 0.05). L: low mast cell count; I: intermediate mast cell count; H: high mast cell count.
Figure 6. MiR-223 and Cxcl10 relative expression in K14-HPV16 transgenic mice. (a) MiR-223 relative expression in chest skin and ear tissue, together. There are statistically significant differences between the low and intermediate mast cell counts (* p < 0.05) and between the low and high mast cell counts (** p < 0.01). (b) Cxcl10 relative expression in chest skin and ear tissue, together. We observed statistically significant results between the low and intermediate mast cell counts (** p < 0.01), between the low and high mast cell counts (*** p < 0.001), and between the intermediate and high mast cell counts (* p < 0.05). (c) Normalized relative expression of miR-223 in chest skin. (d) Normalized relative expression of Cxcl10 in chest skin. Statistically significant differences were found between the low and intermediate mast cell counts (** p < 0.01). (e) MiR-223 relative expression in ear tissue. (f) Cxcl10 relative expression in ear tissue. Statistically significant results were observed between the low and high mast cell counts (** p < 0.01) and between the intermediate and high mast cell counts (* p < 0.05). L: low mast cell count; I: intermediate mast cell count; H: high mast cell count.
Cancers 14 02216 g006
Cancers 14 02216 g007 550
Figure 7. MiR-125b and Ccr2 relative expression in K14-HPV16 transgenic mice. (a) MiR-125b relative expression in chest skin and ear tissue, together. There are statistically significant differences between the low and high mast cell counts (** p < 0.01). (b) Ccr2 relative expression in chest skin and ear tissue, together. We observed statistically significant results between the low and intermediate mast cell counts (* p < 0.05) and between the intermediate and high mast cell counts (* p < 0.05). (c) Normalized relative expression of miR-125b in chest kin. Statistically significant differences were observed between the low and intermediate mast cell counts (* p < 0.05). (d) Normalized relative expression of Ccr2 in chest skin. (e) MiR-125b relative expression in ear tissue. Statistically significant differences were found between the intermediate and high mast cell counts (* p < 0.05). (f) Ccr2 relative expression in ear tissue. L: low mast cell count; I: intermediate mast cell count; H: high mast cell count.
Figure 7. MiR-125b and Ccr2 relative expression in K14-HPV16 transgenic mice. (a) MiR-125b relative expression in chest skin and ear tissue, together. There are statistically significant differences between the low and high mast cell counts (** p < 0.01). (b) Ccr2 relative expression in chest skin and ear tissue, together. We observed statistically significant results between the low and intermediate mast cell counts (* p < 0.05) and between the intermediate and high mast cell counts (* p < 0.05). (c) Normalized relative expression of miR-125b in chest kin. Statistically significant differences were observed between the low and intermediate mast cell counts (* p < 0.05). (d) Normalized relative expression of Ccr2 in chest skin. (e) MiR-125b relative expression in ear tissue. Statistically significant differences were found between the intermediate and high mast cell counts (* p < 0.05). (f) Ccr2 relative expression in ear tissue. L: low mast cell count; I: intermediate mast cell count; H: high mast cell count.
Cancers 14 02216 g007
Table
Table 1. Differential mast cell count in different histological lesions of both organs together (a), and chest skin (b) and ear tissue (c) separately. Low, intermediate and high mast cell count values were obtained in the IBM SPSS software platform by ranking the mast cell count in terciles; Low represents a mean mast cell per camp of 0 to 2 mast cells, intermediate refers to a mean of 2 to 9 mast cells per camp, and high represents a mean of 9 to 33 mast cells per camp.
Table 1. Differential mast cell count in different histological lesions of both organs together (a), and chest skin (b) and ear tissue (c) separately. Low, intermediate and high mast cell count values were obtained in the IBM SPSS software platform by ranking the mast cell count in terciles; Low represents a mean mast cell per camp of 0 to 2 mast cells, intermediate refers to a mean of 2 to 9 mast cells per camp, and high represents a mean of 9 to 33 mast cells per camp.
a
HistologyTotal
NormalHyperplasiaDysplasiaCarcinoma
Mast cells Average
(Chest skin + Ear)		Low mast cells1130014
Intermediate mast cells348015
High mast cells077115
Total141415144
b
Histology Total
NormalHyperplasiaDysplasia
Mast cells Average
(Chest skin)		Low mast cells630 9
Intermediate mast cells046 10
High mast cells011 2
Total687 21
c
HistologyTotal
NormalHyperplasiaDysplasiaCarcinoma
Mast cells Average
(Ear)		Low mast cells50005
Intermediate mast cells30205
High mast cells066113
Total868123
Table
Table 2. Top 10 hub proteins.
Table 2. Top 10 hub proteins.
RankNameScoreProtein
1ENSMUSP000001324533.03 × 1027CCR2
2ENSMUSP000000252633.03 × 1027IL4
2ENSMUSP000000008893.03 × 1027TNF
4ENSMUSP000000313223.03 × 1027CXCL15
5ENSMUSP000000313273.03 × 1027CXCL1
6ENSMUSP000000476463.03 × 1027CXCL10
7ENSMUSP000000270613.03 × 1027IL13
7ENSMUSP000000206503.03 × 1027IL17A
9ENSMUSP000000748853.03 × 1027CXCL2
10ENSMUSP000000684683.03 × 1027ITGAM
Table
Table 3. MicroRNAs found in online databases with respective targets and scores (data last accessed in 10 April 2021). TargetScan 7.2 (https://www.targetscan.org/mmu_72/) relies on conservation of binding sites and divides miRNA families into broadly conserved (conserved across most vertebrates, usually to zebrafish), highly conserved (meaning that they are conserved across most mammals, but usually not beyond placenta) and poorly conserved (all others) [24]. miRDB (http://mirdb.org/) prediction scores are between 50 and 100, assigned by new computational target prediction algorithms (>80 considered most likely to be real) [19]. miRmap (http://mirmap.ezlab.org), ranks potential targets according to repression strength (mediated by RISC), from 0 to 100%, with 100% representing the strongest repression [25]. miRwalk (http://mirwalk.umm.uni-heidelberg.de), incorporates putative targets of 13 prediction datasets, with results being expressed in p-values between 0 and 1, representing the probability that a candidate target site is a true target site (high p-values are better) [26]. miRTarBase 8.0 (http://miRTarBase.mbc.nctu.edu.tw/)/TarBase v.8 (http://www.microrna.gr/tarbase) are experimentally validated by low throughput techniques (reporter assays, qPCR, Western blot or enzyme linked immunosorbent assays) or high throughput techniques (microarray or proteomics experiments) [23,27]. HITS-CLIP, high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation; LRA/WB, luciferase reporter assay/Western blot.
Table 3. MicroRNAs found in online databases with respective targets and scores (data last accessed in 10 April 2021). TargetScan 7.2 (https://www.targetscan.org/mmu_72/) relies on conservation of binding sites and divides miRNA families into broadly conserved (conserved across most vertebrates, usually to zebrafish), highly conserved (meaning that they are conserved across most mammals, but usually not beyond placenta) and poorly conserved (all others) [24]. miRDB (http://mirdb.org/) prediction scores are between 50 and 100, assigned by new computational target prediction algorithms (>80 considered most likely to be real) [19]. miRmap (http://mirmap.ezlab.org), ranks potential targets according to repression strength (mediated by RISC), from 0 to 100%, with 100% representing the strongest repression [25]. miRwalk (http://mirwalk.umm.uni-heidelberg.de), incorporates putative targets of 13 prediction datasets, with results being expressed in p-values between 0 and 1, representing the probability that a candidate target site is a true target site (high p-values are better) [26]. miRTarBase 8.0 (http://miRTarBase.mbc.nctu.edu.tw/)/TarBase v.8 (http://www.microrna.gr/tarbase) are experimentally validated by low throughput techniques (reporter assays, qPCR, Western blot or enzyme linked immunosorbent assays) or high throughput techniques (microarray or proteomics experiments) [23,27]. HITS-CLIP, high-throughput sequencing of RNAs isolated by crosslinking immunoprecipitation; LRA/WB, luciferase reporter assay/Western blot.
MicroRNAsTargetDatabases (Score)
miR-466l-3pIL-4TargetScan (poorly), miRDB (100), miRmap (93,99)
TNFmiRDB (100)
CXCL1miRDB (95)
IL-13miRDB (80)
IL-17ATargetScan (poorly), miRDB (96)
CXCL2miRDB (98), miRmap (98,64)
miR-466kCXCL15miRTarBase (HITS-CLIP), miRDB (100), miRmap (99,98), miRwalk (0.92)
ITGAMTargetScan (poorly), miRDB (100), miRmap (99,95), miRwalk (0.87)
miR-223-3pCXCL10miRDB (90), miRmap (94,46), TarBase (-) *
ITGAMmiRwalk (1), TarBase (-) *
miR-125b-5pTNFmiRTarBase (LRA/WB)
CCR2miRwalk (1), Tarbase (-) *
* Present in the database but with no score associated.
Table
Table 4. Kruskal–Wallis test for the normalized relative expression of the four microRNAs along mast cell count.
Table 4. Kruskal–Wallis test for the normalized relative expression of the four microRNAs along mast cell count.
Normalized Relative Expression (−ΔCt)p Value
miR-466l-3p0.408
miR-466k0.999
miR-223-3p0.035
miR-125b-5p0.024
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Costa, A.C.; Santos, J.M.O.; Medeiros-Fonseca, B.; Oliveira, P.A.; Bastos, M.M.S.M.; Brito, H.O.; Gil da Costa, R.M.; Medeiros, R. Characterizing the Inflammatory Microenvironment in K14-HPV16 Transgenic Mice: Mast Cell Infiltration and MicroRNA Expression. Cancers 2022, 14, 2216. https://doi.org/10.3390/cancers14092216

AMA Style

Costa AC, Santos JMO, Medeiros-Fonseca B, Oliveira PA, Bastos MMSM, Brito HO, Gil da Costa RM, Medeiros R. Characterizing the Inflammatory Microenvironment in K14-HPV16 Transgenic Mice: Mast Cell Infiltration and MicroRNA Expression. Cancers. 2022; 14(9):2216. https://doi.org/10.3390/cancers14092216

Chicago/Turabian Style

Costa, Alexandra C., Joana M. O. Santos, Beatriz Medeiros-Fonseca, Paula A. Oliveira, Margarida M. S. M. Bastos, Haissa O. Brito, Rui M. Gil da Costa, and Rui Medeiros. 2022. "Characterizing the Inflammatory Microenvironment in K14-HPV16 Transgenic Mice: Mast Cell Infiltration and MicroRNA Expression" Cancers 14, no. 9: 2216. https://doi.org/10.3390/cancers14092216

APA Style

Costa, A. C., Santos, J. M. O., Medeiros-Fonseca, B., Oliveira, P. A., Bastos, M. M. S. M., Brito, H. O., Gil da Costa, R. M., & Medeiros, R. (2022). Characterizing the Inflammatory Microenvironment in K14-HPV16 Transgenic Mice: Mast Cell Infiltration and MicroRNA Expression. Cancers, 14(9), 2216. https://doi.org/10.3390/cancers14092216

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