J Cancer 2019; 10(19):4588-4595. doi:10.7150/jca.30081

Research Paper

Searching HPV genome for methylation sites involved in molecular progression to cervical precancer

Christine Kottaridi1 Corresponding address, Danai Leventakou1, Abraham Pouliakis1, Vasileios Pergialiotis2, George Chrelias2, Eugenia Patsouri1, Andriani Zacharatou1, Eleni Panopoulou1, Vasileia Damaskou1, Vasileios Sioulas2, Charalambos Chrelias2, Sofia Kalantaridou2, Ioannis G. Panayiotides1

1. 2 nd Department of Pathology, University General Hospital “ATTIKON”, School of Medicine, National and Kapodistrian University of Athens, Athens 12464, Greece
2. 3 rd Department of Gynaecology and Obstetrics, University General Hospital “ATTIKON”, School of Medicine, National and Kapodistrian University of Athens, Athens 12464, Greece

This is an open access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/). See http://ivyspring.com/terms for full terms and conditions.
How to cite this article:
Kottaridi C, Leventakou D, Pouliakis A, Pergialiotis V, Chrelias G, Patsouri E, Zacharatou A, Panopoulou E, Damaskou V, Sioulas V, Chrelias C, Kalantaridou S, Panayiotides IG. Searching HPV genome for methylation sites involved in molecular progression to cervical precancer. J Cancer 2019; 10(19):4588-4595. doi:10.7150/jca.30081. Available from http://www.jcancer.org/v10p4588.htm

Abstract

Background: Human Papilloma Virus has been considered as the main cause for cervical cancer. In this study we investigated epigenetic changes and especially methylation of specific sites of HPV genome. The main goal was to correlate methylation status with histological grade as well as to determine its accuracy in predicting the disease severity by establishing optimum methylation cutoffs.

Methods: In total, sections from 145 cases genotyped as HPV16 were obtained from formalin- fixed, paraffin-embedded tissue of cervical biopsies, conization or hysterectomy specimens. Highly accurate pyrosequencing of bisulfite converted DNA, was used to quantify the methylation percentages of UTR promoter, enhancer and 5' UTR, E6 CpGs 494, 502, 506 and E7 CpGs 765, 780, 790. The samples were separated in different groupings based on the histological outcome. Statistical analysis was performed by SAS 9.4 for Windows and methylation cutoffs were identified by MATLAB programming language.

Results: The most important methylation sites were at the enhancer and especially UTR 7535 and 7553 sites. Specifically for CIN3+ (i.e. HSIL or SCC) discrimination, a balanced sensitivity vs. specificity (68.1%, 66.2% respectively) with positive predictive value (PPV) and negative predictive value (NPV) (66.2%, 68.2% respectively) was achieved for UTR 7535 methylation of 6.1% cutoff with overall accuracy 67.1%, while for UTR 7553 a sensitivity 60.9%, specificity 69.0%, PPV=65.6%, NPV=64.5% and overall accuracy=65.0% at threshold 10.1% was observed.

Conclusion: Viral HPV16 genome was found methylated in NF-1 binding sites of UTR in cases with high grade disease. Methylation percentages of E6 and E7 CpG sites were elevated at the cancer group.

Keywords: HPV16, methylation, UTR, E6, E7

Introduction

Cervical cancer is the fourth most common cancer in women, and the seventh overall, with an estimated 528,000 new cases in 2012. A large majority (around 85%) of the global burden occurs in the less developed regions, where it accounts for almost 12% of all female cancers. Since 1999, Human Papilloma Virus (HPV) has been considered as the main cause for cervical cancer [1]. Over 100 different types of the virus have been isolated from human samples, out of which 20 are considered oncogenic and have been epidemiologically linked to this cancer [2].

A persistent infection by HPV aided by other parameters results in HPV integration and progression to high grade lesions. Most HPV infections are cleared by the host's immune system before they progress to lesions. Prophylactic HPV vaccines and traditional Pap-smear screening are undoubtedly capable of decreasing the incidence and mortality of cervical cancer. However, a large number of females succumb to the disease each year due to late diagnosis and resistance to conventional treatments. Nowadays, the more extensive use of modern molecular biological methods, have added to cervical cancer screening approaches [3].

In order to understand the biology of this tumor, it is of utmost importance to analyze its molecular dynamics aiding to improve the clinical outcome. Epigenetics is a well-established phenomenon that plays a major role in virus-associated neoplasms [4, 5] and one of the most widely studied epigenetic changes is DNA methylation. During cervical carcinogenesis, substantial changes in methylation are observed in both the host cell and the viral genome. Methylation of the viral DNA has been recently proposed as a novel biomarker with encouraging results [6]. The quantification of the percentage of cytosines with a covalently added methyl-group at individual CpG (Cytosine—phosphate—Guanine) dinucleotides reflects the degree of epigenetic changes of the viral genome. Many studies mainly focused on L1 viral gene, have already shown that the methylation percentage of HPV 16 specific CpG sites along this viral gene, is increasing gradually and it is highest in women with high-grade cervical neoplasia [7-10].

This study aims to assess whether quantitative measurement of methylation of CpGs along the HPV 16 UTR, E6 and E7 genes, could predict the presence of high-grade disease at histology in women testing positive for the HPV 16 genotype. More specifically, we aim to correlate specific sites with the histological grade and to determine their accuracy in predicting the disease severity by establishing optimum methylation cutoffs.

Material and Methods

In total, sections from 145 non-pregnant women, 21-62 years of age genotyped as HPV16 were obtained from formalin-fixed, paraffin-embedded tissue of cervical biopsies, conization or hysterectomy specimens of females that visited the gynecology clinic of Attikon University General Hospital, Athens, Greece, between May 2014 and May 2016. Women were included irrespective of their ethnicity, smoking habits, phase in their cycle, menopausal status and contraceptive practices. Women who were HIV or hepatitis B/C positive, with autoimmune disorders, or had a previous history of cervical treatment were excluded. The histological diagnoses included the following groups: normal, CIN1, CIN2, CIN3, squamous cell carcinoma (SCC), and adenocarcinomas.

DNA was extracted from formalin-fixed, paraffin-embedded tissue sections, using QIAamp DNA FFPE Tissue Kit (Qiagen, Heidelberg Germany). All steps in the purification procedure were done using the automated QIAcube technology (Qiagen, Heidelberg Germany). DNA typing was performed with the HPV Genotypes 14 Real-TM Quant (Sacace Biotechnologies, Como Italy) for the quantitative detection and genotyping of Human Papillomavirus types (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66 and 68). The extracted DNA concentrations were measured with QIAexpert technology (Qiagen, Heidelberg Germany).

DNAs were then bisulfite converted using the EpiTect Bisulfite Kit (Qiagen, Heidelberg Germany), according to the manufacturer's instructions and stored at -80°C. Biotin-labeled primer sets, sequencing primers and polymerase chain reaction conditions were used as previously described [11] to amplify HPV16 UTR region, while biotin-labeled primer sets and sequencing primers for E6 and E7 genomic regions were designed in the present study with the aid of PyroMark Assay Design 2.0 (table 1). PCR conditions for E6 amplification were preheating at 95 °C for 5 min, 40 cycles at 95 °C for 30 s, 53 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 10 min, while for E7 amplification were 95 °C for 5 min, 40 cycles at 95 °C for 30 s, 58 °C for 30 s, 72 °C for 30 s, and a final extension at 72 °C for 10 min. The methylation quantification was performed by Pyrosequencing technology (PyroMark Q24, Qiagen, Heidelberg, Germany), which provides a site-specific quantification of methylation at individual CpG sites. Pyrosequencing protocols were at first applied to bisulfite converted SiHa cells. This cervical cell line contains a single genome of HPV-16 integrated into chromosomal DNA which is completely unmethylated at LCR, E6 and E7 genomic regions [12] a finding that was confirmed by our protocols. In order to establish a limit of blank for each specific methylation site, we performed series of ten measurements using SiHa DNA. The highest value that was reported in the series of measurement was 3.4%. DNA from SiHa cells was run in every experiment as an unmethylated control. As of lack of an artificially methylated control for the studied CpG dinoucleotides spanning along HPV16 genome, we randomly chose a highly methylated to all sites clinical sample and performed ten independent experiments to check the reproducibility of our protocols to detect methylation. The measurements had at most 3.2% standard deviation around the mean value. The UTR methylation analysis included CpGs sites which are located within the p97 HPV 16 promoter (31, 37, 43, 52, 58, 7862), the enhancer (7535, 7553) and the 5'UTR (7270, 7461, 7455, 7428). The E6 methylation analysis included CpGs 494, 502 and 506, while the E7 methylation analysis included CpGs 765, 780 and 790 (reference sequence NC_001526).

 Table 1 

Primer sequences and CpGs sites that were studied

UTR
ForwardGTGTGTTTGTATGTATGGTATA
ReverseATACAATAAATAACCACAACACAATTAA-BTN
SequencingGTGTGTTTGTATGTATGGTATA
Sequenced CpGs:7270
Amplicon size:133bp
ForwardTGTAGGTTAGGAAAATAGGGAT
ReverseACTATATTTACTACATCCTATTTTTATT-BTN
SequencingTTGTGTTAAAAAGTATGTAAT
TAGGAAAATAGGGATTTGG
Sequenced CpGs:7428, 7455, 7461
7535, 7553
Amplicon size:205bp
ForwardGTAAAATTGTATATGGGTGTGTGTAAA
ReverseAATCCTAAAACATTACAATTCTCTTT-BTN
SequencingAATAATTTATGTATAAAATTAAGGG
GTATATGGGTGTGTGTAAAT
Sequenced CpGs:31, 37, 43, 52, 58
7862
Amplicon size:189 bp
E6
ForwardATGGGAATTTATATGTTGTATGTGAT
ReverseCTCCTCCTCTAAACTATCATTTAATTACTC-BTN
SequencingAGATTTTATAATATAAGGGGT
Sequenced CpGs:494, 502, 506
Amplicon size:397 bp
E7
ForwardATGAAATAGATGGTTTAGTTGGATAAG-BTN
ReverseAACACACAATTCCTAATATACCCATTAACA
SequencingCCATTAACAAATCTTCCAAAAT
Sequenced CpGs:765, 780, 790
Amplicon size:168 bp

Statistical analysis

The statistical analysis was performed by programming in SAS 9.4 for Windows (SAS Institute Inc. NC, USA) [13, 14]. Microsoft Excel for Windows was used for data storage and preprocessing. We applied the student's t-test to examine if the various methylation percentages were statistically different between various histological groupings. Finally, the algorithms for the determination of the optimum threshold values were developed in-house within the MATLAB software environment and programming language (The MathWorks, Inc. Natick, Massachusetts, U.S.A.).

We calculated different accuracy parameters for the ability of the mean methylation to detect the presence of disease for the previous histological cut- offs. These included the sensitivity (Sens), specificity (Spec), positive (PPV) and negative predictive value (NPV), the false positive (FPR) and false negative (FNR) rate, the overall accuracy (OA) and the positive (PLR) and negative likelihood ratio (NLR). These parameters were calculated from methylation percentage 0% and increased up to 100% using an increment step of 0.1%, as described in our previous studies [7, 15-17]. Subsequently graphs depicting sensitivity, specificity and overall accuracy for all methylation positions were produced and the above measures were reported for two positions, i.e. the methylation position that maximizes OA and the position that produces a more balanced result, i.e. whereas the difference between sensitivity and specificity is minimal.

In order to identify any possible correlations of age or disease severity with the methylation levels and the mean methylation per region it was used the Spearman correlation coefficient (Rs). Specifically concerning disease severity the histological status was numerically coded as 1 to 5 for the histological categories of: normal, CIN-1, CIN2, CIN-3 and carcinomas (SCC and adenocarcinomas) respectively.

Results

Study group profile

Age at diagnosis ranged from 21 to 62 years, with a mean of 35.52 years (SD: ±8.35) and a median of 34 years. Out of the 145 studied cases 69 (47.59%) had single HPV16 infection and from the remaining cases with multiple infections, 50 cases (34.48%) had only high risk infections and 26 (17.93%) cases had high risk and low risk infections. Multiple infections were more frequent in negative and CIN1 cases, specifically 71.43% (n=5) from negative cases and 60.87% (n=14) from CIN1 had multiple infections; for CIN2 and CIN3 the percentages of multiple infections were 53.33% (n=24) and 52.46% (n=32) respectively (Table 2). Finally out of the 6 SCC cases, 5 (83.33%) had single HPV16 infection, and similarly all three cases of adenocarcinomas.

Methylation analysis

Highly accurate pyrosequencing technology which performed after bisulfite treatment of DNA was used to calculate methylation percentage of specific sites harbored at 5'UTR, enhancer and promoter of the viral genome (nt 7270, 7461, 7455, 7428, 31, 37, 43, 52, 58, 7862) as well as the E6 (494, 502, 506) and E7 (765, 780, 790) gene. At figure 1 are presented pyrograms depicting samples either methylated or unmethylated at the genomic regions that were studied.

The mean methylation of 5' UTR (7270, 7461, 7455, 7428) CpG sites between the different histological groups was (Table 3) 5.6±3.4, 5.3±4.0, 5.8±2.5, 6.6±3.3, 13.3±13.1 and 14.5±8.3 (for Normal, CIN-1, CIN-2, CIN-3, SCC and adenocarcinomas respectively). The two CpG sites (7535, 7553) mapped at the enhancer of the genome, had a mean methylation that ranged from 6.3±3.2 to 39.2±19.2 when the studied groups were analyzed (Table 3). Respectively, the p97 promoter CpG sites (7862, 31, 37, 43, 52, 58) depicted methylation levels, varying from 1.4±0.5 to 9.6±13.2 (Table 3).

The methylation of E6 CpG sites ranged from 8.0±6.9 to 40.6±27.8 among the histological groups (table 3) and the methylation of E7 CpG sites ranged 3.6±1.8 to 29.6±27.2 among the various histological groups (table 3).

Methylation and age

Concerning the relation of the 5' UTR region sites analyzed in this study the analysis via Spearman correlation indicated that no specific site neither the mean methylation for all 5' UTR sites was related to the women age. Similarly, for the enhancer region it was found a rather low correlation (Rs=0.23, p=0.0053) of age with site 7535, while the correlation for site 7553 was Rs=0.15 (p=0.0802), and for the mean methylation of the enhancer region was Rs=0.20, p=0.0248. The correlation of the promoter sites with age showed a weak correlation with methylation levels at site 7862 (Rs=0.18504, p=0.0280) and there was no correlation of the remaining sites or of the mean methylation of the promoter region. Furthermore, no significant relation of age with the methylation of sites from E6 or E7 region, neither the mean methylation of the E6 or E7 region were identified (p>0.05 for all cases).

Methylation and disease severity

The statistical tests correlating disease severity (as formed in a numeric manner) with methylation percentages indicated generally weak correlations. Specifically, for the 5'UTR sites the correlation for sites 7270, 7461, 7428 and the mean methylation of all 5' UTR sites was Rs 0.22, 0.36, 0.24 and 0.29 respectively, p<0.05, while for sites 7535, 7553 and for the mean enhancer methylation a correlation of 0.33, 0.29 and 0.35 respectively, p<0.001 was found. A negative correlation for site 43 of the promoter (Rs -0.23, p=0.0052) was indicated. As far as E6 is concerned, a weak correlation existed for sites 494, 502 and the mean E6 methylation (Rs 0.18, 0.20 and 0.17 respectively, p<0.05) and finally for site 790 of the E7 region the correlation was 0.19, p=0.0265.

 Figure 1 

Pyrograms depicting unmethylated (left column) and methylated samples (right column). Grey blocs are presenting the under investigation CG sites and numbers at blue squares the methylation percentages. A & B, UTR genomic region sites 7535, 7553. C & D, E6 genomic region sites 494, 502, 506. E & F, E7 genomic region sites 765, 780, 790 (a reverse primer was used for pyrosequencing of this region).

J Cancer Image (Click on the image to enlarge.)
 Table 2 

Distribution of single (HPV16 only) and multiple HPV infections of the studied population, in relation to the histological groups and the age using as cut-off the age of 30 years.

Multiple InfectionsSingle Infection
Histology/age group<=30>30Subtotal<=30>30SubtotalTotal
Negative55 (71.4%)112 (28.5%)7 (4.8%)
CIN-17714 (60.8%)369 (39.1%)23 (15.8%)
CIN-2101424 (53.3%)81321 (46.6%)45 (31.0%)
CIN-3131932 (52.4%)72229 (47.5%)61 (42.0%)
SCC11 (16.6%)55 (83.3%)6 (4.1%)
Adenocarcinoma33 (100%)3 (2.0%)
Total30 (39.4%)46 (60.5%)76 (52.4%)19 (27.5%)50 (72.4%)69 (47.5%)145
 Table 3 

Mean methylation and standard deviation according to the histological status and methylation site, grouped according to HPV region.

RegionSiteNormal (N=7)CIN1 (N=23)CIN2 (N=45)CIN3 (N=61)SCC (N=6)AdenoCa (N=3)
5'UTR727011.1±11.812.4±14.914.0±10.014.7±10.220.0±20.628.0±17.0
74613.7±2.32.6±1.32.7±1.14.3±2.512.0±10.110.0±9.5
74554.2±0.93.1±1.43.4±1.14.1±4.214.2±16.011.3±9.0
74283.2±1.22.9±1.63.1±1.63.4±1.511.0±8.69.0±6.5
Mean5.6±3.45.3±4.05.8±2.56.6±3.313.3±13.114.5±8.3
Enhancer75359.4±5.58.1±6.16.3±3.211.1±8.427.6±14.121.3±10.0
755311.1±9.512.0±9.110.1±5.514.6±11.239.2±19.221.3±5.0
Mean10.2±7.210.10±7.48.2±4.112.9±9.033.4±17.321.3±7.5
Promoter78622.8±0.92.4±1.21.7±0.61.7±0.83.1±1.73.6±2.0
311.4±0.52.2±1.02.3±1.71.8±1.37.5±9.99.0±11.3
372.8±1.52.3±1.23.0±1.62.6±1.96.8±9.29.0±10.4
432.0±1.02.3±1.12.2±3.61.4±1.22.6±3.29.3±12.7
522.4±1.51.8±0.92.2±1.51.9±1.47.0±9.29.6±13.2
582.8±2.42.3±1.22.4±1.72.2±1.57.0±10.39.6±11.5
Mean2.4±1.02.2±0.72.3±1.31.9±1.05.6±6.88.3±10.2
E64948.4±7.58.3±6.98.8±7.810.2±8.325.3±16.112.6±12.4
50212.0±11.913.7±12.014.3±13.516.2±12.340.6±27.818.0±14.0
5068.0±6.98.1±6.78.2±7.99.2±8.225.3±17.111.0±11.3
Mean9.4±8.710.0±8.410.4±9.611.9±9.530.4±20.113.8±12.5
E776510.0±4.88.0±3.910.3±6.910.2±6.829.6±17.24.5±2.1
7805.0±2.45.3±3.15.3±2.95.5±4.123.6±11.34.5±2.1
7903.6±1.85.6±3.95.0±3.06.2±5.324.0±13.33.0±0.0
Mean6.2±2.76.3±3.46.9±3.97.3±5.025.7±13.74.0±0.0

Calculation of thresholds discriminating between histological groups and performance characteristics

Accuracy parameters were determined for different histological cut-offs. We searched for a methylation threshold optimizing overall accuracy as well as for a threshold that produces a balanced sensitivity vs. specificity (table 4). The most important identified methylation sites were on the enhancer region, sites 7535 and 7553. Specifically for CIN3+ discrimination, a balanced sensitivity vs. specificity (68.1%, 66.2% respectively) with positive predictive value (PPV) and negative predictive value (NPV) (66.2%,68.2% respectively) was achieved for UTR 7535 methylation of 6.1% cutoff with overall accuracy 67.1%, while for UTR 7553 a sensitivity 60.9%, specificity 69.0%, PPV 65.6%, NPV 64.5% and overall accuracy 65.0% at threshold 10.1% was observed. For the remaining methylation sites we also calculated the thresholds that maximized overall accuracy or produced a balanced sensitivity vs. specificity. However, results produced were not performing well (data not shown).

 Table 4 

Performance characteristics for discriminating CIN3+ cases using two cut-offs, one that maximizes overall accuracy (Optimal threshold) and one that balances sensitivity vs. specificity (Balanced threshold).

Site 7535Site 7553
Optimal thresholdBalanced thresholdOptimal thresholdBalanced threshold
Sensitivity68.168.175.360.8
Specificity66.266.254.969.0
PPV66.266.261.965.6
NPV68.168.169.664.4
FPR33.833.845.030.9
FNR31.831.824.639.1
OA67.167.165.065.0
OR4.14.13.73.4

PPV: Positive Predictive Value, NPV: Negative Predictive Value, FPR: False Positive Rate, FNR: False Negative Rate, OA: Overall Accuracy, OR: Odds Ratio. Note that Optimal threshold and balanced threshold was the same for site 7535 while for site 7553 the optimal threshold resulted in high sensitivity at the cost of specificity.

 Figure 2 

Histological photos from a CIN3 case (A, magnification 40x) with increased UTR 7535, 7553 sites HPV16 methylation and from a CIN1 case (B, magnification 10x) with absence of methylation at the same genomic viral sites.

J Cancer Image (Click on the image to enlarge.)

Discussion

Many epigenetic alterations are observed including DNA hypomethylation, hypermethylation of tumour suppressor genes and histone modifications during all the stages of cervical cancer [18]. One of the epigenetic mechanisms that are increasingly studied is HPV genes' methylation. At present, a consistent correlation of increased methylation of capsid viral genes with histology severity is referred by the researchers [7, 9, 11, 19, 20]. On the other hand, studies on methylation status of the UTR, E6 and E7 regions reveal heterogeneous and rather inconclusive results [21-25].

DNA methylation has been recognized as a frequent event in cervical cancer and as such, is referred as of valuable tool in the early detection of precancerous disease. In the present study we attempted to elucidate the methylation profile of HPV 16 genome for each of 4 CpG sites of the 5' UTR, 2 of the enhancer, 6 of the promoter, 3 of the E6 and 3 of the E7 gene in clinical specimens of different severities in a Greek women population. The studied sites were located in the 5' UTR at nt 7270, 7461, 7455 and 7428, in the enhancer at 7535, 7553, in the promoter at 7862, 31, 37, 43, 52, 58, in E6 at 494, 502, 506 and in E7 at 765, 780, 790.

The HPV16 UTR plays an important role in regulation of viral gene expression. HPV 16 E6 and E7 oncogenes are transcribed from the P97 promoter which is located at 3' UTR and is regulated by products of the viral E2 gene through a feedback mechanism. The transcription of the viral oncogenes depends also on the enhancer's activity, which is located between the positions 7454 and 7854 and acts as a cis-acting element that drives the transcription of the early HPV genes [26]. Several host transcription factors, such as like AP-1, NF1, SP1, TFIID, TF1, Oct-1 are bound to specific sites of this viral gene triggering the over-production of E6 and E7 oncoproteins gradually leading to neoplastic progression [27, 28]. There are studies supporting that this region is highly or moderately methylated [29-32] and the methylation is associated with the severity of cervical neoplasia. It is assumed that if E2 viral protein does not manage to bind at specific sites due to inhibition by methylation of cytosines within its binding site, the repression of oncogenes' transcription will be diminished. On the other hand, several studies claim that this region has an overall low percent methylation and there is no correlation of methylation with different severities of cervical carcinogenesis [10, 21, 24, 33, 34]. According to the results of the present study, the mean methylation of HPV 16 UTR showed constantly low methylation percentages between the different histological groups. The only sites with remarkable results were 7535 and 7553 mapped at the enhancer, were a correlation with CIN 3+ can be proposed. These sites are part of the binding positions of NF-1 [22, 35, 36]. One could assume that such an increase of methylation could eliminate the NF-1 binding activity, affecting the enhancer's cis-acting efficiency for the transcription of the early HPV genes. But according to our results, although these specific enhancer sites are methylated along with the progression of the disease, we can probably suggest that the massive production of E6, E7 oncoproteins is not affected by the addition of methyl groups at these sites and they may be selectively used to discriminate CIN3+ cases.

In the present study we also investigated the methylation of CpG sites that are located at E6 and E7 viral genes and specifically in those regions that are considered to be immunostimulatory motifs [25]. According to Hacker et al [37], a sequence motif that contains CpGs has the capacity to stimulate certain immune cells. As far as HPV 16 is concerned it has been shown that TLR9 is capable of recognizing a CpG motif between nt 496 and 514 of E6 gene [38] and methylation of the sites located into this genomic part could possibly lead to an escape from immune surveillance, having thus a significant biological impact in cervical carcinogenesis. One of the objectives of Sen et al. [25], was the analysis of the influence of methylation within two immunostimulatory CpG motifs within HPV16 E6 and E7 genes to cervical carcinogenesis. Presence of elevated methylation was shown at cervical cancer samples, with higher proportions at samples that had an integrated HPV 16 infection. In the present study, although we do not have the information concerning the integration of viral genome to the host genome, the percentages of methylation at CpG sites of E6 and E7 genes were elevated at the cervical cancer histology group when compared to precancer samples a finding that is in accordance to published results. As of statistical analysis revealed, these specific sites originated from E6 and E7 gene, cannot be proposed as biomarkers that could distinguish precancerous HPV 16 infections that have a true oncogenic potential from those that will dissolve without leading to disease.

In conclusion, the knowledge of viral genome's alterations during the viral life cycle is adding valuable information on understanding the biology of cervical cancer and to the exploring of new biomarkers. Frequently, different methodological approaches to the methylation study of viral genes may lead to inconsistent results, so the scientific community should feed the literature with findings of this area of HPV research. According to our results, methylation of HPV 16 UTR is not highly associated with severity of cervical neoplasm. However, two specific sites mapped at the enhancer region of UTR, could probably act as biomarkers and molecular determinants that could distinguish the rare HPV 16 infections that have a true oncogenic malignant potential from those common infections from HPV 16 that resolve spontaneously without leading to disease. Although a modest number of samples were studied, the reproducibility of these results should be assessed in large validation sets. Future studies should also analyze serial samples from larger cohorts to further assess the value of methylation as a predictive and diagnostic molecular determinant.

Competing Interests

The authors have declared that no competing interest exists.

References

1. zur Hausen H. Papillomaviruses and cancer: from basic studies to clinical application. Nature reviews Cancer. 2002;2:342-50

2. Munoz N, Bosch FX, de Sanjose S, Herrero R, Castellsague X, Shah KV. et al. Epidemiologic classification of human papillomavirus types associated with cervical cancer. N Engl J Med. 2003;348:518-27

3. Leal SM Jr, Gulley ML. Current and Emerging Molecular Tests for Human Papillomavirus-Related Neoplasia in the Genomic Era. J Mol Diagn. 2017;19:366-77

4. Minarovits J, Demcsak A, Banati F, Niller HH. Epigenetic Dysregulation in Virus-Associated Neoplasms. Adv Exp Med Biol. 2016;879:71-90

5. Kuss-Duerkop SK, Westrich JA, Pyeon D. DNA Tumor Virus Regulation of Host DNA Methylation and Its Implications for Immune Evasion and Oncogenesis. Viruses. 2018:10 82 doi:10.3390/v10020082

6. Nedjai B, Reuter C, Ahmad A, Banwait R, Warman R, Carton J. et al. Molecular progression to cervical precancer, epigenetic switch or sequential model?. Int J Cancer. 2018;143:1720-30

7. Kottaridi C, Kyrgiou M, Pouliakis A, Magkana M, Aga E, Spathis A. et al. Quantitative Measurement of L1 Human Papillomavirus Type 16 Methylation for the Prediction of Preinvasive and Invasive Cervical Disease. J Infect Dis. 2017;215:764-71

8. Torres-Rojas FI, Alarcon-Romero LDC, Leyva-Vazquez MA, Ortiz-Ortiz J, Mendoza-Catalan MA, Hernandez-Sotelo D. et al. Methylation of the L1 gene and integration of human papillomavirus 16 and 18 in cervical carcinoma and premalignant lesions. Oncol Lett. 2018;15:2278-86

9. Chaiwongkot A, Niruthisard S, Kitkumthorn N, Bhattarakosol P. Quantitative methylation analysis of human papillomavirus 16L1 gene reveals potential biomarker for cervical cancer progression. Diagn Microbiol Infect Dis. 2017;89:265-70

10. Clarke MA, Wentzensen N, Mirabello L, Ghosh A, Wacholder S, Harari A. et al. Human papillomavirus DNA methylation as a potential biomarker for cervical cancer. Cancer Epidemiol Biomarkers Prev. 2012;21:2125-37

11. Mirabello L, Schiffman M, Ghosh A, Rodriguez AC, Vasiljevic N, Wentzensen N. et al. Elevated methylation of HPV16 DNA is associated with the development of high grade cervical intraepithelial neoplasia. Int J Cancer. 2013;132:1412-22

12. Badal V, Chuang LS, Tan EH, Badal S, Villa LL, Wheeler CM. et al. CpG methylation of human papillomavirus type 16 DNA in cervical cancer cell lines and in clinical specimens: genomic hypomethylation correlates with carcinogenic progression. J Virol. 2003;77:6227-34

13. DiMaggio C. SAS for epidemiologists: applications and methods. New York, USA: Springer. 2013

14. Harris M, Taylor J. Medical statistics made easy. London, UK: Martin Dunitz, an imprint of the Taylor & Francis Group. 2003

15. Varlatzidou A, Pouliakis A, Stamataki M, Meristoudis C, Margari N, Peros G. et al. Cascaded learning vector quantizer neural networks for the discrimination of thyroid lesions. Anal Quant Cytol Histol. 2011;33:323-34

16. Pouliakis A, Margari C, Margari N, Chrelias C, Zygouris D, Meristoudis C. et al. Using classification and regression trees, liquid-based cytology and nuclear morphometry for the discrimination of endometrial lesions. Diagnostic cytopathology. 2014;42:582-91

17. Zygouris D, Pouliakis A, Margari N, Chrelias C, Terzakis E, Koureas N. et al. Classification of endometrial lesions by nuclear morphometry features extracted from liquid-based cytology samples: a system based on logistic regression model. Anal Quant Cytopathol Histpathol. 2014;36:189-98

18. Lu Q, Ma D, Zhao S. DNA methylation changes in cervical cancers. Methods Mol Biol. 2012;863:155-76

19. Yang-Chun F, Yuan Z, Cheng-Ming L, Yan-Chun H, Xiu-Min M. Increased HPV L1 gene methylation and multiple infection status lead to the difference of cervical epithelial cell lesion in different ethnic women of Xinjiang, China. Medicine (Baltimore). 2017;96:e6409

20. Sun C, Reimers LL, Burk RD. Methylation of HPV16 genome CpG sites is associated with cervix precancer and cancer. Gynecol Oncol. 2011;121:59-63

21. Fernandez AF, Rosales C, Lopez-Nieva P, Grana O, Ballestar E, Ropero S. et al. The dynamic DNA methylomes of double-stranded DNA viruses associated with human cancer. Genome Res. 2009;19:438-51

22. Ding DC, Chiang MH, Lai HC, Hsiung CA, Hsieh CY, Chu TY. Methylation of the long control region of HPV16 is related to the severity of cervical neoplasia. Eur J Obstet Gynecol Reprod Biol. 2009;147:215-20

23. Milutin Gasperov N, Sabol I, Planinic P, Grubisic G, Fistonic I, Corusic A. et al. Methylated Host Cell Gene Promoters and Human Papillomavirus Type 16 and 18 Predicting Cervical Lesions and Cancer. PLoS One. 2015;10:e0129452

24. Kalantari M, Calleja-Macias IE, Tewari D, Hagmar B, Lie K, Barrera-Saldana HA. et al. Conserved methylation patterns of human papillomavirus type 16 DNA in asymptomatic infection and cervical neoplasia. J Virol. 2004;78:12762-72

25. Sen S, Mandal P, Bhattacharya A, Kundu S, Roy Chowdhury R, Mondal NR. et al. Impact of viral and host DNA methylations on HPV16-related cervical cancer pathogenesis. Tumour Biol. 2017;39:1010428317699799

26. Cripe TP, Haugen TH, Turk JP, Tabatabai F, Schmid PG 3rd, Durst M. et al. Transcriptional regulation of the human papillomavirus-16 E6-E7 promoter by a keratinocyte-dependent enhancer, and by viral E2 trans-activator and repressor gene products: implications for cervical carcinogenesis. EMBO J. 1987;6:3745-53

27. Desaintes C, Demeret C. Control of papillomavirus DNA replication and transcription. Semin Cancer Biol. 1996;7:339-47

28. Thierry F, Howley PM. Functional analysis of E2-mediated repression of the HPV18 P105 promoter. New Biol. 1991;3:90-100

29. Cheung JL, Cheung TH, Yu MY, Chan PK. Virological characteristics of cervical cancers carrying pure episomal form of HPV16 genome. Gynecol Oncol. 2013;131:374-9

30. Bhattacharjee B, Sengupta S. CpG methylation of HPV 16 LCR at E2 binding site proximal to P97 is associated with cervical cancer in presence of intact E2. Virology. 2006;354:280-5

31. Thain A, Jenkins O, Clarke AR, Gaston K. CpG methylation directly inhibits binding of the human papillomavirus type 16 E2 protein to specific DNA sequences. J Virol. 1996;70:7233-5

32. Kim K, Garner-Hamrick PA, Fisher C, Lee D, Lambert PF. Methylation patterns of papillomavirus DNA, its influence on E2 function, and implications in viral infection. J Virol. 2003;77:12450-9

33. Brandsma JL, Sun Y, Lizardi PM, Tuck DP, Zelterman D, Haines GK 3rd. et al. Distinct human papillomavirus type 16 methylomes in cervical cells at different stages of premalignancy. Virology. 2009;389:100-7

34. Kalantari M, Garcia-Carranca A, Morales-Vazquez CD, Zuna R, Montiel DP, Calleja-Macias IE. et al. Laser capture microdissection of cervical human papillomavirus infections: copy number of the virus in cancerous and normal tissue and heterogeneous DNA methylation. Virology. 2009;390:261-7

35. Apt D, Chong T, Liu Y, Bernard HU. Nuclear factor I and epithelial cell-specific transcription of human papillomavirus type 16. J Virol. 1993;67:4455-63

36. Chong T, Apt D, Gloss B, Isa M, Bernard HU. The enhancer of human papillomavirus type 16: binding sites for the ubiquitous transcription factors oct-1, NFA, TEF-2, NF1, and AP-1 participate in epithelial cell-specific transcription. J Virol. 1991;65:5933-43

37. Hacker G, Redecke V, Hacker H. Activation of the immune system by bacterial CpG-DNA. Immunology. 2002;105:245-51

38. Hasan UA, Bates E, Takeshita F, Biliato A, Accardi R, Bouvard V. et al. TLR9 expression and function is abolished by the cervical cancer-associated human papillomavirus type 16. J Immunol. 2007;178:3186-97

Author contact

Corresponding address Corresponding author: Christine Kottaridi, MSc, MPH, PhD; phone: +30-210-5831959; E-mail: ckottaridicom


Received 2018-9-20
Accepted 2019-5-8
Published 2019-8-7