Some novel insights on HPV16 related cervical cancer pathogenesis based on analyses of LCR methylation, viral load, E7 and E2/E4 expressions

The Hallmarks of Cervical Cancer: Molecular Mechanisms Induced by Human Papillomavirus

Human papillomaviruses (HPVs) and, specifically, high-risk HPVs (HR-HPVs) are identified as necessary factors in the development of cancer of the lower genital tract, with CaCU standing out as the most prevalent tumor. This review summarizes ten mechanisms activated by HR-HPVs during cervical carcinogenesis, which are broadly associated with at least seven of the fourteen distinctive physiological capacities of cancer in the newly established model by Hanahan in 2022. These mechanisms involve infection by human papillomavirus, cellular tropism, genetic predisposition to uterine cervical cancer (CaCU), viral load, viral physical state, regulation of epigenetic mechanisms, loss of function of the E2 protein, deregulated expression of E6/E7 oncogenes, regulation of host cell protein function, and acquisition of the mesenchymal phenotype.

Identification of HPV16 positive cervical cancer subsets characterized by divergent immune and oncogenic phenotypes with potential implications for immunotherapy

BACKGROUND

Cervical cancers (CaCx), like many other cancer types, portray high molecular heterogeneity that affects response to therapy, including immunotherapy. In India and other developing countries, CaCx mortality rates are very high because women report to the clinics with advanced cancers in absence of organized screening programs. This calls for implementation of newer therapeutic regimens for CaCx, like immunotherapy, which is again not used commonly in such countries.

OBJECTIVE

Therefore, we focused on dissecting tumour immune heterogeneity, if any, identify immune gene-based biomarkers of heterogeneity and subsets of such cancers with the potential for immunotherapy. We also attempted to characterize the cancer-associated phenotypes of such subsets, including viral load, to decipher the relationship of tumour immunogenicity with oncogenicity.

METHODS

Employing RNA-seq analysis of 44 HPV16 positive CaCx patients, immune subtypes were identified by unsupervised hierarchical clustering of global immune-gene expression profiles. Proportions of tumor infiltrating immune cells in the tumor milieu were estimated, employing Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT), using gene expression data from RNA-seq. The oncogenic phenotypes of the immune subtypes of CaCx were deciphered through differential gene expression (DEGs) and pathway enrichment analysis. Viral load was estimated through TaqMan-based qRT-PCR analysis.

RESULTS

Analysis revealed the presence of two immune subtypes of CaCx, A (26/44; 59.09%) and B (18/44; 40.90%). Compared to Subtype-A, Subtype-B portrayed overexpression of immune genes and high infiltration of immune cells, specifically CD8+ T cells (p < 0.0001). Besides, a significant correlation between PD-1 and PD-L1 co-expression among Subtype-B, as opposed to Subtype-A, confirmed the interactive roles of these immune checkpoint molecules in Subtype B. Stepwise discriminant analysis pin-pointed ten immune-genes that could classify 100% of the patients significantly (p < 0.0001) into the two immune subtypes and serve as potential biomarkers of CaCx immunity. Differential gene expression analysis between the subtypes unveiled that Subtype-B was more biologically aggressive than Subtype-A, reflecting loss of structural integrity and promotion of cancer progression. The viral load was significantly lower in Subtype-B (average viral load = 10.74/100 ng of genomic DNA) compared to Subtype-A (average viral load = 14.29/100 ng of genomic DNA). Thus viral load and the ten-gene panel underscore their association with immunogenicity and oncogenicity.

CONCLUSION

Our study provides strong evidence that only a subset, about 41% of HPV16 positive CaCx patients in India, portray immune enrichment of the tumor milieu coupled with aggressive phenotypes. Such subtypes are therefore likely to benefit through checkpoint molecule-based or tumor infiltrating lymphocyte-based immunotherapy, which could be a leap forward in tackling aggressive forms of such CaCx in India and other developing countries.

Predominance of genomically defined A lineage of HPV16 over D lineage in Indian patients from eastern India with squamous cell carcinoma of the cervix in association with distinct oncogenic phenotypes

•

The HPV16 lineage A viruses are predominantly associated with cervical squamous cell carcinoma (SCC) in the eastern region of India.

•

Among lineage A viruses, the E5(Y44L, I65V), E6(L83V) and LCR: C7577T variants are associated with SCC.

•

The risk alleles are distributed among 10 clades comprising of 64 HPV16 genomes.

•

Low integration, high episomal copy numbers and high E7 mRNA expression in SCC contrasts lineage A viruses from lineage D.

Human papillomavirus type-16 (HPV16) is classified into lineages, A, B, C and D and 10 sub-lineages portraying variable infectivity, persistence, and cytological outcomes, however, with geographical variations. Our objective was to delineate the distinctive features of lineages among cervical squamous cell carcinoma (SCC) in the eastern region of India. A total of 145 SCC cases and 24 non-malignant specimens, harboring episomal HPV16, were included. The presence of higher proportion of lineage A over D was observed among SCC cases (86.89% A1, 8.97% D1 and 4.14% D2), while only A1 sub-lineage viruses were found among control specimens. Among the A1 viruses, an association of variants in the E5 (Y44L, I65V), E6 (L83V) genes and LCR: C7577T with SCC, with combined Odd's ratio (95% CI) of 20.5(4.61–91.25) was observed. Network analyses revealed the presence of 10 clades of lineage A viruses comprising of 64 HPV16 genomes harboring the risk alleles. High episomal HPV16 DNA copy numbers (adjusted p-value= 0.0271) and E7 mRNA expression (p-value=0.000017) predominated in SCC with lineage A, over D. Our study highlights the distinctive modalities of oncogenicity among different HPV16 lineages.

Cervical cancer subtypes harbouring integrated and/or episomal HPV16 portray distinct molecular phenotypes based on transcriptome profiling of mRNAs and miRNAs

Heterogeneity in cervical cancers (CaCx) in terms of HPV16 physical status prompted us to investigate the mRNA and miRNA signatures among the different categories of CaCx samples. We performed microarray-based mRNA expression profiling and quantitative real-time PCR-based expression analysis of some prioritised miRNAs implicated in cancer-related pathways among various categories of cervical samples. Such samples included HPV16-positive CaCx cases that harboured either purely integrated HPV16 genomes (integrated) and those that harboured episomal viral genomes, either pure or concomitant with integrated viral genomes (episomal), which were compared with normal cervical samples that were either HPV negative or positive for HPV16. The mRNA expression profile differed characteristically between integrated and episomal CaCx cases for enriched biological pathways. miRNA expression profiles also differed among CaCx cases compared with controls (upregulation—miR-21, miR-16, miR-205, miR-323; downregulation—miR-143, miR-196b, miR-203, miR-34a; progressive upregulation—miR-21 and progressive downregulation—miR-143, miR-34a, miR-196b and miR-203) in the order of HPV-negative controls, HPV16-positive non-malignant samples and HPV16-positive CaCx cases. miR-200a was upregulated in HPV16-positive cervical tissues irrespective of histopathological status. Expression of majority of the predicted target genes was negatively correlated with their corresponding miRNAs, irrespective of the CaCx subtypes. E7 mRNA expression correlated positively with miR-323 expression among episomal cases and miR-203, among integrated cases. miR-181c expression was downregulated only among the episomal CaCx cases and negatively correlated with protein coding transcript of the proliferative target gene, CKS1B of the significantly enriched “G2/M DNA Damage Checkpoint Regulation” pathway among CaCx cases. Thus, the two CaCx subtypes are distinct entities at the molecular level, which could be differentially targeted for therapy. In fact, availability of a small molecule inhibitor of CKS1B, suggests that drugging CKS1B could be a potential avenue of treating the large majority of CaCx cases harbouring episomal HPV16.

Viral DNA integration and methylation of human papillomavirus type 16 in high-grade oral epithelial dysplasia and head and neck squamous cell carcinoma

This study evaluated the integration and methlyation of human papillomavirus type 16 (HPV16) in head and neck squamous cell carcinoma (HNSCC) and its oral precursor, high-grade oral epithelial dysplasia (hgOED). Archival samples of HPV16-positive hgOED (N = 19) and HNSCC (N = 15) were evaluated, along with three HNSCC (UMSCC-1, -47 and -104) and two cervical cancer (SiHa and CaSki) cell lines. HgOED cases were stratified into three groups with increasing degrees of cytologic changes (mitosis, karyorrhexis and apoptosis). The viral load was higher and the E2/E6 ratio lower (indicating a greater tendency toward viral integration) in group 3 than in groups 1 or 2 (p = 0.002, 0.03). Methylation was not observed in hgOED cases and occurred variably in only three HNSCC cases (26.67%, 60.0% and 93.3%). In HNSCC cell lines, lower E7 expression correlated with higher levels of methylation. HgOED with increased cytologic change, now termed HPV-associated oral epithelial dysplasia (HPV-OED), exhibited an increased viral load and a tendency toward DNA integration, suggesting a potentially increased risk for malignant transformation. More detailed characterization and clinical follow-up of HPV-OED patients is needed to determine whether HPV-OED is a true precursor to HPV-associated HNSCC and to clarify the involvement of HPV in HNSCC carcinogenesis.

5azadC treatment upregulates miR-375 level and represses HPV16 E6 expression

High-risk human papillomaviruses are the etiological agents of cervical cancer and HPV16 is the most oncogenic genotype. Immortalization and transformation of infected cells requires the overexpression of the two viral oncoproteins E6 and E7 following HPV DNA integration into the host cell genome. Integration often leads to the loss of the E2 open reading frame and the corresponding protein can no longer act as a transcriptional repressor on p97 promoter. Recently, it has been proposed that long control region methylation also contributes to the regulation of E6/E7 expression.

To determine which epigenetic mechanism is involved in HPV16 early gene regulation, 5-aza-2′-deoxycytidine was used to demethylate Ca Ski and SiHa cell DNA. Decreased expression of E6 mRNA and protein levels was observed in both cell lines in an E2-independent manner. E6 repression was accompanied by neither a modification of the main cellular transcription factor expression involved in long control region regulation, nor by a modification of the E6 mRNA splicing pattern. In contrast, a pronounced upregulation of miR-375, known to destabilize HPV16 early viral mRNA, was observed. Finally, the use of miR-375 inhibitor definitively proved the involvement of miR-375 in E6 repression. These results highlight that cellular DNA methylation modulates HPV16 early gene expression and support a role for epigenetic events in high-risk HPV associated-carcinogenesis.

Expression signatures of HOX cluster genes in cervical cancer pathogenesis: Impact of human papillomavirus type 16 oncoprotein E7

The Homeobox (HOX) genes encode important transcription factors showing deregulated expression in several cancers. However, their role in cervical cancer pathogenesis, remains largely unexplored. Herein, we studied their association with Human Papillomavirus type 16 (HPV16) mediated cervical cancers. Our previously published gene expression microarray data revealed a significant alteration of 12 out of 39 HOX cluster members among cervical cancer cases, in comparison to the histopathologically normal controls. Of these, we validated seven (HOXA10, HOXA13, HOXB13, HOXC8, HOXC9, HOXC11 and HOXD10) by quantitative real-time PCR. We identified decreased HOXA10 expression as opposed to the increased expression of the rest. Such decrease was independent of the integration status of HPV16 genome, but correlated negatively with E7 expression in clinical samples, that was confirmed in vitro. HOXA10 and HOXB13 revealed association with Epithelial-Mesenchymal Transition (EMT). While HOXA10 expression correlated positively with E-Cadherin and negatively with Vimentin expression, HOXB13 showed the reverse trend. Chromatin immunoprecipitation study in vitro revealed the ability of E7 to increase HOX gene expression by epigenetic regulation, affecting the H3K4me3 and H3K27me3 status of their promoters, resulting from a loss of PRC2-LSD1 complex activity. Thus, besides identifying the deregulated expression of HOX cluster members in HPV16 positive cervical cancer and their association with EMT, our study highlighted the mechanism of HPV16 E7-mediated epigenetic regulation of HOX genes in such cancers, potentially serving as bedrock for functional studies in the future.

Methylation of CpG sites in the upstream regulatory region, physical status and mRNA expression of HPV-6 in adult-onset laryngeal papilloma

The methylation status of HPV-6 upstream regulatory region (URR) in adult-onset laryngeal papillomatosis (AO-LP) remains unclear. The purpose of this study was to investigate the methylation status of URR and the physical status of HPV-6, as well as the dynamic variations of viral load and mRNA expression in AO-LP. We examined 18 specimens from 11 patients with AO-LP by real-time polymerase chain reaction (PCR), bisulfite-sequencing PCR, and amplification of papilloma oncogene transcripts. HPV-6 was identified in 9 of 11 patients (81.8%), and all the 15 specimens derived from 9 HPV-6-positive cases contained only episomal HPV-6 transcripts with intact E2. Three HPV-6-positive patients developed recurrent lesions, and HPV-6 copy numbers and mRNA expression decreased after surgical treatment. Among the 96 CpG sites (16/case), 67 (69.8%) were unmethylated, while 23 (30.2%) were heterogeneous (≥ 1 methylated CpG clone). High viral loads and episomal status of HPV-6 were frequently observed in AO-LP; thus, persistent E6/E7 mRNA expression of LR-HPV-6 may be associated with AO-LP recurrences. Hypomethylation and scattered patterns of methylated CpGs at the URR of HPV-6 were identified in AO-LP.

Impact of viral and host DNA methylations on HPV16-related cervical cancer pathogenesis

Epigenetic alterations within human papillomavirus (HPV) and host cellular genomes are known to occur during cervical carcinogenesis. Our objective was to analyse the influence of (1) methylation within two immunostimulatory CpG motifs within HPV16 E6 and E7 genes around the viral late promoter and their correlation, if any, with expression deregulation of host receptor (TLR9) and DNA methyltransferases (DNMT1, DNMT3A and DNMT3B) and (2) global DNA methylation levels within CpGs of the repetitive Alu sequences, on cervical cancer (CaCx) pathogenesis. Significantly higher proportions of CaCx samples portrayed methylation in immunostimulatory CpG motifs, compared to HPV16-positive non-malignant samples, with cases harbouring episomal HPV16 showing decreased methylation compared to those with viral integration. A significant linear trend of TLR9 upregulation was recorded in the order of HPV-negative controls < HPV16-positive non-malignant samples < HPV16-positive CaCx cases. TLR9 upregulation in cases with episomal HPV16 was again higher among those with non-methylated immunostimulatory CpG motifs. Comparison of cases with HPV-negative controls revealed that DNMT3A was significantly downregulated only among integrated cases, DNMT3B was significantly overexpressed among both categories of cases, although at variable levels, while DNMT1 failed to show any deregulated expression among the cases. Global host DNA hypomethylation, also showed a significant linear increasing trend through the progressive CaCx development stages mentioned above and was most prominently higher among cases with episomal HPV16 as opposed to viral integration. Thus, HPV16 and host methylations appear to influence CaCx pathogenesis, with differential molecular signatures among CaCx cases with episomal and integrated HPV16.

Human papillomavirus (HPV) genome status & cervical cancer outcome--A retrospective study

Background & objectives:

Persistent infections with high-risk (HR) human papillomaviruses such as HPV 16, 18, 31, 33 and 45 have been identified as the major aetiological factor for cervical cancer. The clinical outcome of the disease is often determined by viral factors such as viral load, physical status and oncogene expression. The aim of the present study was to evaluate the impact of such factors on clinical outcome in HPV16 positive, locally advanced cervical cancer cases.

Methods:

One hundred and thirty two pretreatment cervical tumour biopsies were selected from patients undergoing radiotherapy alone (n=63) or concomitant chemo-radiation (n=69). All the samples were positive for HPV 16. Quantitative real time-PCR was carried out to determine viral load and oncogene expression. Physical status of the virus was determined for all the samples by the ratio of E2_copies/E7_copies; while in 73 cases, the status was reanalyzed by more sensitive APOT (amplification of papillomavirus oncogene transcripts) assay. Univariate analysis of recurrence free survival was carried out using Kaplan-Meier method and for multivariate analysis the Cox proportional hazard model was used.

Results:

The median viral load was 19.4 (IQR, 1.9- 69.3), with viral integration observed in 86 per cent cases by combination of the two methodologies. Both univariate and multivariate analyses identified viral physical status as a good predictor of clinical outcome following radiation treatment, with episomal form being associated with increased recurrence free survival.

Interpretation & conclusions:

The present study results showed that viral physical status might act as an important prognostic factor in cervical cancer.

The Role of Epigenomics in the Study of Cancer Biomarkers and in the Development of Diagnostic Tools

Epigenetics plays a key role in cancer development. Genetics alone cannot explain sporadic cancer and cancer development in individuals with no family history or a weak family history of cancer. Epigenetics provides a mechanism to explain the development of cancer in such situations. Alterations in epigenetic profiling may provide important insights into the etiology and natural history of cancer. Because several epigenetic changes occur before histopathological changes, they can serve as biomarkers for cancer diagnosis and risk assessment. Many cancers may remain asymptomatic until relatively late stages; in managing the disease, efforts should be focused on early detection, accurate prediction of disease progression, and frequent monitoring. This chapter describes epigenetic biomarkers as they are expressed during cancer development and their potential use in cancer diagnosis and prognosis. Based on epigenomic information, biomarkers have been identified that may serve as diagnostic tools; some such biomarkers also may be useful in identifying individuals who will respond to therapy and survive longer. The importance of analytical and clinical validation of biomarkers is discussed, along with challenges and opportunities in this field.

Effects of Methylation Status of CpG Sites within the HPV16 Long Control Region on HPV16-Positive Head and Neck Cancer Cells

Objective

To map comprehensively the methylation status of the CpG sites within the HPV16 long control region (LCR) in HPV-positive cancer cells, and to explore further the effects of methylation status of HPV16 LCR on cell bioactivity and E6 and E7 expression. In addition, to analyze the methylation status of the LCR in HPV-positive oropharyngeal squamous cell carcinoma (OPSCC) patients.

Methods and Materials

Methylation patterns of HPV16 LCR in UM-SCC47, CaSki, and SiHa cells and HPV16-positiive OPSCC specimens were detected by bisulfite-sequencing PCR and TA cloning. For cells treated with 5-aza-2′-deoxycytidine and E6 and E7 knockdown, MTS and trypan blue staining, annexin-V and 7-AAD staining, and prodidium iodide were used to evaluate cell growth and cell proliferation, cell apoptosis, and cell cycle arrest, respectively. E6 and E7 mRNA and protein expression were analyzed by quantitative real-time PCR and immunocytochemistry, respectively.

Results

Hypermethylation status of the LCR in UM-SCC47 (79.8%) and CaSki cells (90.0%) and unmethylation status of the LCR in SiHa cells (0%) were observed. Upon demethylation, the cells with different methylation levels responded differently during growth, apoptosis, and cell cycle arrest, as well as in terms of their E6 and E7 expression. In HPV16-positive OPSCC patients, the methylation rates were 9.5% in the entire LCR region, 13.9% in the 5′-LCR, 6.0% in the E6 enhancer, and 9.5% in the p97 promoter, and hypermethylation of p97 promoter was found in a subset of cases (20.0%, 2/10).

Conclusions

Our study revealed two different methylation levels of the LCR in HPV16-positive cancer cells and OPSCC patients, which may represent different carcinogenesis mechanisms of HPV-positive cancers cells. Demethylating the meCpGs in HPV16 LCR might be a potential target for a subgroup of HPV16-positive patients with head and neck squamous cell carcinoma.

Systematical analysis of cutaneous squamous cell carcinoma network of microRNAs, transcription factors, and target and host genes

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNA molecules found in multicellular eukaryotes which are implicated in development of cancer, including cutaneous squamous cell carcinoma (cSCC). Expression is controlled by transcription factors (TFs) that bind to specific DNA sequences, thereby controlling the flow (or transcription) of genetic information from DNA to messenger RNA. Interactions result in biological signal control networks.

MATERIALS AND METHODS

Molecular components involved in cSCC were here assembled at abnormally expressed, related and global levels. Networks at these three levels were constructed with corresponding biological factors in term of interactions between miRNAs and target genes, TFs and miRNAs, and host genes and miRNAs. Up/down regulation or mutation of the factors were considered in the context of the regulation and significant patterns were extracted.

RESULTS

Participants of the networks were evaluated based on their expression and regulation of other factors. Sub-networks with two core TFs, TP53 and EIF2C2, as the centers are identified. These share self-adapt feedback regulation in which a mutual restraint exists. Up or down regulation of certain genes and miRNAs are discussed. Some, for example the expression of MMP13, were in line with expectation while others, including FGFR3, need further investigation of their unexpected behavior.

CONCLUSIONS

The present research suggests that dozens of components, miRNAs, TFs, target genes and host genes included, unite as networks through their regulation to function systematically in human cSCC. Networks built under the currently available sources provide critical signal controlling pathways and frequent patterns. Inappropriate controlling signal flow from abnormal expression of key TFs may push the system into an incontrollable situation and therefore contributes to cSCC development.

Bridging Links between Long Noncoding RNA HOTAIR and HPV Oncoprotein E7 in Cervical Cancer Pathogenesis

Human Papillomavirus (HPV) type 16 oncoprotein E7 plays a major role in cervical carcinogenesis by interacting with and functionally inactivating various host regulatory molecules. Long noncoding RNA (lncRNA) HOTAIR is one such regulator that recruits chromatin remodelling complex PRC2, creating gene silencing H3K27 me3 marks. Hence, we hypothesized that HOTAIR could be a potential target of E7, in HPV16 related cervical cancers (CaCx). We identified significant linear trend of progressive HOTAIR down-regulation through HPV negative controls, HPV16 positive non-malignants and CaCx samples. Majority of CaCx cases portrayed HOTAIR down-regulation in comparison to HPV negative controls, with corresponding up-regulation of HOTAIR target, HOXD10, and enrichment of cancer related pathways. However, a small subset had significantly higher HOTAIR expression, concomitant with high E7 expression and enrichment of metastatic pathways. Expression of HOTAIR and PRC2-complex members (EZH2 and SUZ12), showed significant positive correlation with E7 expression in CaCx cases and E7 transfected C33A cell line, suggestive of interplay between E7 and HOTAIR. Functional inactivation of HOTAIR by direct interaction with E7 could also be predicted by in silico analysis and confirmed by RNA-Immunoprecipitation. Our study depicts one of the causal mechanisms of cervical carcinogenesis by HPV16 E7, through modulation of HOTAIR expression and function.

Viral load, integration and methylation of E2BS3 and 4 in human papilloma virus (HPV) 16-positive vaginal and vulvar carcinomas

Objective

To investigate if viral load, integration and methylation of E2BS3 and 4 represent different ways of tumor transformation in vaginal and vulvar carcinoma and to elucidate its clinical impact.

Methods

Fifty-seven samples, positive for HPV16, were selected for the study. Detection of viral load was made with realtime-PCR using copy numbers of E6 and integration was calculated from comparing E2 to E6-copies. Methylation of E2BS3 and 4 was analysed using bisulphite treatment of tumor DNA, followed by PCR and pyrosequencing.

Results

Vaginal tumors were found to have a higher viral load (p = 0.024) compared to vulvar tumors but a high copy number (> median value, 15 000) as well as high methylation (>50%) was significantly (p = 0.010 and p = 0.045) associated with a worse cancer-specific survival rate in vulvar carcinoma, but not in vaginal carcinoma. Four groups could be defined for the complete series using a Cluster Two step analysis; (1) tumors holding episomal viral DNA, viral load below 150 000 copies not highly methylated (n = 25, 46.3%); (2) tumors harboring episomal viral DNA and being highly methylated (>50%; n = 6, 11.1%); (3) tumors with viral DNA fully integrated (n = 11, 20.4%), and (4) tumors harboring episomal viral DNA and being medium- or unmethylated (<50%) and having a high viral load (> total mean value 150 000; n = 12, 22.2%). The completely integrated tumors were found to be distinct group, whilst some overlap between the groups with high methylation and high viral load was observed.

Conclusion

HPV16- related integration, methylation in E2BS3 and 4 and viral load may represent different viral characteristics driving vaginal and vulvar carcinogenesis. HPV16- related parameters were found to be of clinical importance in the vulvar series only.

The role of the globular heads of the C1q receptor in HPV-16 E2-induced human cervical squamous carcinoma cell apoptosis via a mitochondria-dependent pathway

Background

Human papillomavirus type-16 (HPV-16) E2 protein acts as a transcriptional modulator and plays a key role in regulating many biological responses. The purpose of this study was to investigate the relationship between HPV-16 E2, the receptor for the globular heads of human C1q (gC1qR) gene expression, mitochondrial dysfunction and apoptosis regulation in human cervical squamous carcinoma cells (C33a and SiHa).

Methods

HPV-16 E2 and gC1qR expression was examined using real-time PCR and western blot analysis. Apoptosis in C33a and SiHa cells was assessed by flow cytometry. Mitochondrial function was detected via ROS generation, the amount of cytosolic Ca²⁺, and changes in the mitochondrial membrane potential (Δψm).

Results

The expression of the HPV-16 E2 and gC1qR gene significantly decreased in human cervical squamous carcinoma samples relative to the non-cancerous cervix samples. C33a and SiHa cells that were transfected with a vector encoding HPV-16 E2 displayed significantly increased gC1qR gene expression and mitochondrial dysfunction as well as an up-regulation of cellular apoptosis, which was abrogated by the addition of gC1qR small-interfering RNA (siRNA).

Conclusions

These data support a mechanism whereby gC1qR plays an important role in HPV-16 E2-induced human cervical squamous carcinoma cell apoptosis via a mitochondria-dependent pathway.

Electronic supplementary material

The online version of this article (doi:10.1186/s12967-014-0286-y) contains supplementary material, which is available to authorized users.

Folate and vitamin B12 may play a critical role in lowering the HPV 16 methylation-associated risk of developing higher grades of CIN

We previously reported that a higher degree of methylation of CpG sites in the promoter (positions 31, 37, 43, 52, and 58) and enhancer site 7862 of human papillomavirus (HPV) 16 was associated with a lower likelihood of being diagnosed with HPV 16-associated CIN 2+. The purpose of this study was to replicate our previous findings and, in addition, to evaluate the influence of plasma concentrations of folate and vitamin B12 on the degree of HPV 16 methylation (HPV 16m). The study included 315 HPV 16-positive women diagnosed with either CIN 2+ or ≤CIN 1. Pyrosequencing technology was used to quantify the degree of HPV 16m. We reproduced the previously reported inverse association between HPV 16m and risk of being diagnosed with CIN 2+. In addition, we observed that women with higher plasma folate and HPV 16m or those with higher plasma vitamin B12 and HPV 16m were 75% (P < 0.01) and 60% (P = 0.02) less likely to be diagnosed with CIN 2+, respectively. With a tertile increase in the plasma folate or vitamin B12, there was a 50% (P = 0.03) and 40% (P = 0.07) increase in the odds of having a higher degree of HPV 16m, respectively. This study provides initial evidence that methyl donor micronutrients, folate and vitamin B12, may play an important role in maintaining a desirably high degree of methylation at specific CpG sites in the HPV E6 promoter and enhancer that are associated with the likelihood of being diagnosed with CIN 2+.

Mapping of HPV transcripts in four human cervical lesions using RNAseq suggests quantitative rearrangements during carcinogenic progression

Two classes of Human papillomaviruses (HPV) infect the anogenital track: high risk viruses that are associated with risk of cervical cancer and low risk types that drive development of benign lesions, such as condylomas. In the present study, we established quantitative transcriptional maps of the viral genome in clinical lesions associated with high risk HPV16 or low risk HPV6b. Marked qualitative and quantitative changes in the HPV16 transcriptome were associated with progression from low to high grade lesions. Specific transcripts encoding essential regulatory proteins such as E7, E2, E1^E4 and E5 were identified. We also identified intrinsic differences between the HPV6b-associated condyloma transcript map and that of the HPV16-associated low grade CIN specifically regarding promoter usage. Characterization and quantification of HPV transcripts in patient samples thus establish the impact of viral transcriptional regulation on the status of HPV-associated lesions and may therefore help in defining new biologically-relevant prognosis markers.

Human Papillomavirus 16, 18, 31 and 45 viral load, integration and methylation status stratified by cervical disease stage

Background

Persistent infection with oncogenic Human Papillomavirus (HPV) is associated with the development of cervical cancer with each genotype differing in their relative contribution to the prevalence of cervical disease. HPV DNA testing offers improved sensitivity over cytology testing alone but is accompanied by a generally low specificity. Potential molecular markers of cervical disease include type-specific viral load (VL), integration of HPV DNA into the host genome and methylation of the HPV genome. The aim of this study was to evaluate the relationship between HPV type-specific viral load, integration and methylation status and cervical disease stage in samples harboring HPV16, HPV18, HPV31 or HPV45.

Methods

Samples singly infected with HPV16 (n = 226), HPV18 (n = 32), HPV31 (n = 75) or HPV45 (n = 29) were selected from a cohort of 4,719 women attending cervical screening in England. Viral load and integration status were determined by real-time PCR while 3’L1-URR methylation status was determined by pyrosequencing or sequencing of multiple clones derived from each sample.

Results

Viral load could differentiate between normal and abnormal cytology with a sensitivity of 75% and a specificity of 80% (odds ratio [OR] 12.4, 95% CI 6.2–26.1; p < 0.001) with some variation between genotypes. Viral integration was poorly associated with cervical disease. Few samples had fully integrated genomes and these could be found throughout the course of disease. Overall, integration status could distinguish between normal and abnormal cytology with a sensitivity of 72% and a specificity of 50% (OR 2.6, 95% CI 1.0–6.8; p = 0.054). Methylation levels were able to differentiate normal and low grade cytology from high grade cytology with a sensitivity of 64% and a specificity of 82% (OR 8.2, 95% CI 3.8–18.0; p < 0.001). However, methylation varied widely between genotypes with HPV18 and HPV45 exhibiting a broader degree and higher magnitude of methylated CpG sites than HPV16 and HPV31.

Conclusions

This study lends support for HPV viral load and CpG methylation status, but not integration status, to be considered as potential biomarkers of cervical disease.

Construction and analysis of regulatory genetic networks in cervical cancer based on involved microRNAs, target genes, transcription factors and host genes

Over recent years, genes and microRNA (miRNA/miR) have been considered as key biological factors in human carcinogenesis. During cancer development, genes may act as multiple identities, including target genes of miRNA, transcription factors and host genes. The present study concentrated on the regulatory networks consisting of the biological factors involved in cervical cancer in order to investigate their features and affect on this specific pathology. Numerous raw data was collected and organized into purposeful structures, and adaptive procedures were defined for application to the prepared data. The networks were therefore built with the factors as basic components according to their interacting associations. The networks were constructed at three levels of interdependency, including a differentially-expressed network, a related network and a global network. Comparisons and analyses were made at a systematic level rather than from an isolated gene or miRNA. Critical hubs were extracted in the core networks and notable features were discussed, including self-adaption feedback regulation. The present study expounds the pathogenesis from a novel point of view and is proposed to provide inspiration for further investigation and therapy.

Differential expression of HPV16 L2 gene in cervical cancers harboring episomal HPV16 genomes: influence of synonymous and non-coding region variations

We tested the hypothesis that (i) synonymous variations within the coding regions, and (ii) variations within the non-coding regions of HPV, influence cervical cancer (CaCx) pathogenesis under the impact of intact HPV16 genomes. Whole genome sequence analysis of HPV16 isolates within 70 CaCx cases and 25 non-malignant samples revealed that synonymous variations were significantly higher within the E6 (p = 0.014), E5 (p = 0.001) and L2 (p = 0.0002) genes of HPV16 isolates within cases, compared to isolates within non-malignant samples. All of the 25 (100%) humanized codons identified within L2 ORF of the samples analyzed, were harbored by CaCx cases, while 8 out of 25 (32%) were harbored by HPV16 positive non-malignant samples (p = 3.87105E-07). L2 (mRNA and protein) expression was evident only among cases with episomal viral genomes and L2 mRNA expression correlated significantly with E2 gene copy numbers suggesting expression from all episomal genomes. Among such cases, Asian American (AA) isolates portrayed all of the humanized codons (100%; 4–6/sample) recorded within L2, which was significantly higher (p = 2.02E-7) compared to the European (E) isolates (22.8%; none or 1–2/sample). Additionally, majority of E variant isolates within cases (54/57; 94.7%) portrayed a variation (T4228C) within the short non-coding region (NCR2) between E5 and L2 genes, which portrays a weak promoter activity specific for L2 mRNA expression. This resulted in loss of 9 out of 14 miRNA binding sites (hsa-miR-548 family), despite the significant overexpression of miR548a-5p and miR548d-5p among such cases (28.64 and 36.25 folds, respectively), in comparison to HPV negative control samples. The findings exemplify the biological relevance of sequence variations in HPV16 genomes and highlight that episomal HPV16 in CaCx cases employ multiple mechanisms to sustain L2 expression, thereby justifying the potential role of L2 in such cancers, as opposed to those harboring viral integration.

Productive infection of bovine papillomavirus type 2 in the urothelial cells of naturally occurring urinary bladder tumors in cattle and water buffaloes

Background

Papillomaviruses (PVs) are highly epitheliotropic as they usually establish productive infections within squamous epithelia of the skin, the anogenital tract and the oral cavity. In this study, early (E) and late (L) protein expression of bovine papillomavirus type 2 (BPV-2) in the urothelium of the urinary bladder is described in cows and water buffaloes suffering from naturally occurring papillomavirus-associated urothelial bladder tumors.

Methods and Findings

E5 protein, the major oncoprotein of the BPV-2, was detected in all tumors. L1 DNA was amplified by PCR, cloned and sequenced and confirmed to be L1 DNA. The major capsid protein, L1, believed to be only expressed in productive papillomavirus infection was detected by Western blot analysis. Immunohistochemical investigations confirmed the presence of L1 protein both in the cytoplasm and nuclei of cells of the neoplastic urothelium. Finally, the early protein E2, required for viral DNA replication and known to be a pivotal factor for both productive and persistent infection, was detected by Western blot and immunohistochemically. Electron microscopic investigations detected electron dense particles, the shape and size of which are consistent with submicroscopic features of viral particles, in nuclei of neoplastic urothelium.

Conclusion

This study shows that both active and productive infections by BPV-2 in the urothelium of the bovine and bubaline urinary bladder can occur in vivo.