Characterization of sequence variations within HPV16 isolates among Indian women: prediction of causal role of rare non-synonymous variations within intact isolates in cervical cancer pathogenesis

Molecular epidemiology of human papillomavirus variants in cervical cancer in India

Background & objectives Cervical cancer (CC) has been documented as the fourth most common cancer worldwide. Persistent infections with high-risk human papillomavirus (hr-HPV) have been suggested in the development of CC. Although prophylactic vaccines are available for the prevention of prevalent hr-HPV types, intra-type variations exist within a particular HPV type that has varying oncogenic potential as well as the mechanism of pathogenicity and varying neutralization by antibodies. Therefore, we carried out a systematic review to determine the distribution of HPV intra-typic variations in different geographical locations of India and their reported implications. Methods Preferred Reporting Items for Systematic Reviews and Meta-Analysis (PRISMA) guidelines were followed to retrieve relevant articles from the standard databases using appropriate keywords. Consequently, 17 articles were included in the current review after screening based on inclusion and exclusion criteria. Results The majority of articles included in this review reported variations within the HPV16 E6 gene, followed by the L1 and E7 genes. Analysis of available data indicated the differential regional distribution of some variations. These variations have also been reported to impact the biological functions of various viral proteins. Interpretation & conclusions The distribution of lineages varied with the different genomic regions sequenced. Additionally, there were certain unique and common variations in the HPV genome with respect to geographical regions. Hence, we suggest the identification of region-specific variations for the development of diagnostic and prognostic interventions.

E7 Oncogene HPV58 Variants Detected in Northeast Brazil: Genetic and Functional Analysis

Cervical cancer is associated with persistent infections by high-risk Human Papillomavirus (HPV) types that may have nucleotide polymorphisms and, consequently, different oncogenic potentials. Therefore, this study aimed to evaluate the genetic variability and structural effects of the E7 oncogene of HPV58 in cervical scraping samples from Brazilian women. The study was developed with patients from hospitals in the metropolitan area of Recife, PE, Brazil. The most frequent HPV types were, in descending order of abundance, HPV16, 31, and 58. Phylogenetic analysis demonstrated that the isolates were classified into sublineages A2, C1, and D2. Two positively selected mutations were found in E7: 63G and 64T. The mutations G41R, G63D, and T64A in the E7 protein reduced the stability of the protein structure. Utilizing an NF-kB reporter assay, we observed a decrease in the NK-kB pathway activity with the HPV58-E7 variant 54S compared to the WT E7. The other detected E7 HPV58 variants presented similar NF-kB pathway activity compared to the WT E7. In this study, it was possible to identify mutations that may interfere with the molecular interaction between the viral oncoproteins and host proteins.

Mutation Profile of HPV16 L1 and L2 Genes in Different Geographic Areas

The causal relationship between HPV and cervical cancer in association with the high prevalence of high risk HPV genotypes led to the design of HPV vaccines based on the major capsid L1 protein. In recent years, capsid protein L2 has also become a focal point in the field of vaccine research. The present review focuses on the variability of HPV16 L1 and L2 genes, emphasizing the distribution of specific amino acid changes in the epitopes of capsid proteins. Moreover, a substantial bioinformatics analysis was conducted to describe the worldwide distribution of amino acid substitutions throughout HPV16 L1, L2 proteins. Five amino acid changes (T176N, N181T; EF loop), (T266A; FG loop), (T353P, T389S; HI loop) are frequently observed in the L1 hypervariable surface loops, while two amino acid substitutions (D43E, S122P) are adjacent to L2 specific epitopes. These changes have a high prevalence in certain geographic regions. The present review suggests that the extensive analysis of the amino acid substitutions in the HPV16 L1 immunodominant loops may provide insights concerning the ability of the virus in evading host immune response in certain populations. The genetic variability of the HPV16 L1 and L2 epitopes should be extensively analyzed in a given population.

Genetic variability in minor capsid protein (L2 gene) of human papillomavirus type 16 among Indian women

Human papillomavirus type 16 (HPV-16) is the predominant genotype worldwide associated with invasive cervical cancer and hence remains as the focus for diagnostic development and vaccine research. L2, the minor capsid protein forms the packaging unit for the HPV genome along with the L1 protein and is primarily associated with transport of genomic DNA to the nucleus. Unlike L1, L2 is known to elicit cross-neutralizing antibodies and thus becomes a suitable candidate for pan-HPV prophylactic vaccine development. In the present study, a total of 148 cervical HPV-16 isolates from Indian women were analyzed by PCR-directed sequencing, phylogenetic analysis and in silico immunoinformatics tools to determine the L2 variations that may impact the immune response and oncogenesis. Ninety-one SNPs translating to 35 non-synonymous amino acid substitutions were observed, of these 16 substitutions are reported in the Indian isolates for the first time. T245A, L266F, S378V and S384A substitutions were significantly associated with high-grade cervical neoplastic status. Multiple substitutions were observed in samples from high-grade cervical neoplastic status as compared to those from normal cervical status (p = 0.027), specifically from the D3 sub-lineage. It was observed that substitution T85A was part of both, B and T cell epitopes recognized by MHC-I molecules; T245A was common to B and T cell epitopes recognized by MHC-II molecules and S122P/A was common to the region recognized by both MHC-I and MHC-II molecules. These findings reporting L2 protein substitutions have implications on cervical oncogenesis and design of next-generation L2-based HPV vaccines.

HPV16 and HPV18 type-specific APOBEC3 and integration profiles in different diagnostic categories of cervical samples

Human papillomavirus (HPV) 16 and 18 are the most predominant types in cervical cancer. Only a small fraction of HPV infections progress to cancer, indicating that additional factors and genomic events contribute to the carcinogenesis, such as minor nucleotide variation caused by APOBEC3 and chromosomal integration. We analysed intra-host minor nucleotide variants (MNVs) and integration in HPV16 and HPV18 positive cervical samples with different morphology. Samples were sequenced using an HPV whole genome sequencing protocol TaME-seq. A total of 80 HPV16 and 51 HPV18 positive samples passed the sequencing depth criteria of 300× reads, showing the following distribution: non-progressive disease (HPV16 n = 21, HPV18 n = 12); cervical intraepithelial neoplasia (CIN) grade 2 (HPV16 n = 27, HPV18 n = 9); CIN3/adenocarcinoma in situ (AIS) (HPV16 n = 27, HPV18 n = 30); cervical cancer (HPV16 n = 5). Similar numbers of MNVs in HPV16 and HPV18 samples were observed for most viral genes, with the exception of HPV18 E4 with higher numbers across clinical categories. APOBEC3 signatures were observed in HPV16 lesions, while similar mutation patterns were not detected for HPV18. The proportion of samples with integration was 13% for HPV16 and 59% for HPV18 positive samples, with a noticeable portion located within or close to cancer-related genes.

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.

Structural and functional impacts of E5 genetic variants of human papillomavirus type 31

Persistent infections caused by high-risk human papillomavirus (HR-HPV) are important, for the development of cervical lesions, but environmental and genetic factors are also related in the process of carcinogenesis. Among the genetic factors, the genetic variants of HR-HPV appear to be related to the risk of persistent infections. Therefore, the present study investigates variants of HPV31 E5 oncogene in cervical scraping samples from Brazilian women to assess their functional and structural effects, in order to identify possible repercussions of these variants on the infectious and carcinogenic process. Our results detected nucleotide changes previously described in the HPV31 E5 oncogene, which may play a critical role in the development of cancer due to its ability to promote cell proliferation and signal transmission. In our study, the interaction percentage of the 31E5 sequence generated by the Immune Epitope Server database and the Analysis Resource (IEDB) allowed us to include possible immunogenic epitopes with the MHC-I and MHC-II molecules, which may represent a possible relationship between protein suppression of the immune system. In the structural analysis of the HPV31 E5 oncoprotein, the N5D, I48 V, P56A, F80I and V64I polymorphisms can be found inserted within transmembrane regions. The P56A mutation has been predicted to be highly stabilizing and, therefore, can cause a change in protein function. Regarding the interaction of the E5 protein from HPV31 with the signaling of NF-kB pathway, we observed that in all variants of the E5 gene from HPV-31, the activity of the NF-kB pathway was increased compared to the prototype. Our study contributes to a more refined design of studies with the E5 gene from HPV31 and provides important data for a better understanding of how variants can be distinguished under their clinical consequences.

Prevalence of Colistin-Resistant, Carbapenem-Hydrolyzing Proteobacteria in Hospital Water Bodies and Out-Falls of West Bengal, India

Indiscriminate use of antibiotics has resulted in a catastrophic increase in the levels of antibiotic resistance in India. Hospitals treat critical bacterial infections and thus can serve as reservoirs of multidrug resistant (MDR) bacteria. Hence, this study was conducted to gauge the prevalence patterns of MDR bacteria in hospital wastewater. Water samples collected from 11 hospitals and 4 environmental sources belonging to 5 most-densely populated districts of West Bengal, India were grown on MacConkey and Eosin Methylene Blue agar. A total of 84 (hospital-associated = 70, environmental water sources = 14) isolates were characterized. The predominant species found in water from hospital-associated areas (HAA) were Acinetobacter baumannii (22.9%), Escherichia coli (28.6 %), and Klebsiella pneumoniae (25.7%). Greater than 75% of the HAA isolates were found to be mcr-1 gene negative and colistinresistant. Meropenem non-susceptibility was also high among the HAA isolates at 58.6%, with the presence of the carbapenemase gene and bla_NDM in 67.1% of the non-susceptible isolates. Among the three predominant species, significantly higher numbers of E. coli isolates were found to be non-susceptible to meropenem ((80%), p-value = 0.00432) and amikacin (AK (90%), p-value = 0.00037). This study provides evidence for the presence of high numbers of colistin-resistant and carbapenem-hydrolyzing Proteobacteriain hospital wastewater.

HPV16 variant analysis in primary and recurrent CIN2/3 lesions demonstrates presence of the same consensus variant

INTRODUCTION

Recurrent cervical intraepithelial lesions (rCIN2/3) after treatment of CIN2/3 occur in 5-15% of cases. rCIN2/3 can result from incomplete resection of CIN2/3, where the same HPV type and variant remains present. rCIN2/3 could also occur following a new infection with a different HPV variant of the same HPV type as the initial lesion. This study investigates HPV16 consensus variants in paired HPV16 positive scrapes from baseline CIN2/3 and rCIN2/3 lesions.

METHODS

Paired HPV16 positive cervical scrapes of women with CIN2/3 at baseline and rCIN2/3 6 or 12 months after treatment were selected for whole-genome amplification and Illumina sequencing. Sequences were compared and nucleotide changes over time were characterized.

RESULTS

From 14 paired samples, 10 had identical consensus variants in baseline CIN2/3 and rCIN2/3. Four paired samples showed one to three nucleotide variations at recurrent disease compared to baseline.

CONCLUSION

Identical or nearly identical HPV16 consensus variants were found in scrapes of paired HPV16 positive baseline CIN2/3 and rCIN2/3 lesions after treatment, suggesting no need for HPV variant analysis when the same HPV type is found in both lesions. These results argue for either incomplete excision of baseline CIN2/3 or inability of clearance of the original HPV infection.

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.

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.

Disparities in Cervical Cancer Incidence and Mortality: Can Epigenetics Contribute to Eliminating Disparities?

Screening for uterine cervical intraepithelial neoplasia (CIN) followed by aggressive treatment has reduced invasive cervical cancer (ICC) incidence and mortality. However, ICC cases and carcinoma in situ (CIS) continue to be diagnosed annually in the United States, with minorities bearing the brunt of this burden. Because ICC peak incidence and mortality are 10-15 years earlier than other solid cancers, the number of potential years of life lost to this cancer is substantial. Screening for early signs of CIN is still the mainstay of many cervical cancer control programs. However, the accuracy of existing screening tests remains suboptimal. Changes in epigenetic patterns that occur as a result of human papillomavirus infection contribute to CIN progression to cancer, and can be harnessed to improve existing screening tests. However, this requires a concerted effort to identify the epigenomic landscape that is reliably altered by HPV infection specific to ICC, distinct from transient changes.

Complete Genome Sequences of Eight Human Papillomavirus Type 16 Asian American and European Variant Isolates from Cervical Biopsies and Lesions in Indian Women

Human papillomavirus type 16 (HPV16), a member of the Papillomaviridae family, is the primary etiological agent of cervical cancer. Here, we report the complete genome sequences of four HPV16 Asian American variants and four European variants, isolated from cervical biopsies and scrapings in India.

Identification and validation of immunogenic potential of India specific HPV-16 variant constructs: In-silico &in-vivo insight to vaccine development

Cervical cancer is one of the most common gynecological cancers in the world but in India, it is the top most cancer among women. Persistent infection with high-risk human papillomaviruses (HR-HPVs) is the most important risk factor. The sequence variation(s) in the most common HR-HPV i.e. HPV type 16 leads to altered biological functions with possible clinical significance in the different geographical locations. Sixteen major variants (V1-V16) in full length L1 gene of HPV-16 were identified following analysis of 250 prospectively collected cervical cancer tissue biopsies and their effect on immunogenicity was studied. The effect of these major variations on the epitopes were predicted by in silico methods and the immunogenicity of variants and respective reference DNA vaccine constructs were evaluated by administration of prepared DNA vaccine constructs in female BALB/c mice to evaluate antibody titer. In the present study, L500F (V16) variation showed a significant ~2.7 fold (p < 0.002) increase in antibody titer, whereas T379P (V8) showed ~0.4 fold (p < 0.328) decrease after final injection. These results showed a promising roadmap for the development of DNA based vaccine and for the generation of effective response, though there is a need to study more prevalent variants of HPV in the Indian population.

Identification of human papillomavirus-16 E6 variation in cervical cancer and their impact on T and B cell epitopes

The infection with high-risk human papillomavirus (HR-HPV) is the most important risk factor for development of cervical cancer. The intra-type variations of HPV have different biological and pathological consequences with respect to disease progression. In the present study, six major Indian variants were experimentally identified in E6 gene of HPV-16 and showed their impact on immunogenicity by in silico methods. Four different phylogenetic lineages were observed in sequences including European (E) prototype, European variant, Asian and American Asian variant classes and complete absence of African phylogenetic lineages. On the prediction of B- and T-cell epitopes, 18 and 23 potent epitopes for MHC-II alleles, 10 potent MHC-I and 15 B-cell epitopes in each reference and variant sequence were identified. Interestingly, the presence of variation H78Y and L83V result in creation of four new epitopes for the HLA-DQA1*0101/DQB1*0501. Out of 15 B-cell predicted epitopes, three most potent epitopes were identified in both reference and variant sequence. Notably the amino acid stretch from amino acid 16-60 and 76-94 are very important for the immunological properties of E6 protein because these regions contain majority of the predicted epitopes. In future, this could control the cervical cancer by targeting these amino acid stretches for the development of HPV-16 vaccine.

Interleukin-12 gene adjuvant increases the immunogenicity of virus-like particles of human papillomavirus type 16 regional variant strain

Objectives

To analyze the immunogenicity of virus-like particles (VLP) of human papillomavirus type 16 (HPV16) isolated in East China and the adjuvant potential of interleukin-12 (IL-12).

Methods

The variant HPV16 L1VLP expressed in sf9 insect cells were purified with cesium chloride gradient centrifugation. BALB/c mice were vaccinated with VLP (L1N), VLP with Freund's adjuvant (L1A) or VLP with IL-12 recombinant plasmid (L1P). HPV16 VLP specific IgG and IFN-γ level in the serum were detected by ELISA, and the percentage of CD4⁺ and CD8⁺ in spleen cells was detected with flow cytometry.

Results

The titers of serum IgG antibodies in vaccinated groups were higher than in negative control and the serum antibodies mainly recognized conformation-dependent HPV16 VLP epitopes. Splenic CD4⁺ and CD8⁺ T cell subsets increased after vaccination in every experimental group, and CD8⁺ increased obviously in L1P group. The ratio of CD4⁺/CD8⁺ decreased in L1P group and increased in the other two groups, compared to control group. Vaccination induced specific secretion of IFN-γ in the serum of vaccinated group (p < 0.05), especially in the L1P group.

Conclusions

VLP of HPV16 variant strain isolated in East China could induce humoral immunity and cellular immunity in mice, and IL-12 recombinant plasmid can enhance cellular immunity.

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.

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

This study was undertaken to decipher the interdependent roles of (i) methylation within E2 binding site I and II (E2BS-I/II) and replication origin (nt 7862) in the long control region (LCR), (ii) expression of viral oncogene E7, (iii) expression of the transcript (E7-E1/E4) that encodes E2 repressor protein and (iv) viral load, in human papillomavirus 16 (HPV16) related cervical cancer (CaCx) pathogenesis. The results revealed over-representation (p<0.001) of methylation at nucleotide 58 of E2BS-I among E2-intact CaCx cases compared to E2-disrupted cases. Bisulphite sequencing of LCR revealed overrepresentation of methylation at nucleotide 58 or other CpGs in E2BS-I/II, among E2-intact cases than E2-disrupted cases and lack of methylation at replication origin in case of both. The viral transcript (E7-E1/E4) that produces the repressor E2 was analyzed by APOT (amplification of papillomavirus oncogenic transcript)-coupled-quantitative-RT-PCR (of E7 and E4 genes) to distinguish episomal (pure or concomitant with integrated) from purely integrated viral genomes based on the ratio, E7 C(T)/E4 C(T). Relative quantification based on comparative C(T) (threshold cycle) method revealed 75.087 folds higher E7 mRNA expression in episomal cases over purely integrated cases. Viral load and E2 gene copy numbers were negatively correlated with E7 C(T) (p = 0.007) and E2 C(T) (p<0.0001), respectively, each normalized with ACTB C(T), among episomal cases only. The k-means clustering analysis considering E7 C(T) from APOT-coupled-quantitative-RT-PCR assay, in conjunction with viral load, revealed immense heterogeneity among the HPV16 positive CaCx cases portraying integrated viral genomes. The findings provide novel insights into HPV16 related CaCx pathogenesis and highlight that CaCx cases that harbour episomal HPV16 genomes with intact E2 are likely to be distinct biologically, from the purely integrated viral genomes in terms of host genes and/or pathways involved in cervical carcinogenesis.

Characterization and whole genome analysis of human papillomavirus type 16 e1-1374^63nt variants

BACKGROUND

The variation of the most common Human papillomavirus (HPV) type found in cervical cancer, the HPV16, has been extensively investigated in almost all viral genes. The E1 gene variation, however, has been rarely studied. The main objective of the present investigation was to analyze the variability of the E6 and E1 genes, focusing on the recently identified E1-1374^63nt variant.

METHODOLOGY/PRINCIPAL FINDINGS

Variation within the E6 of 786 HPV16 positive cervical samples was analyzed using high-resolution melting, while the E1-1374^63nt duplication was assayed by PCR. Both techniques were supplemented with sequencing. The E1-1374^63nt duplication was linked with the E-G350 and the E-C109/G350 variants. In comparison to the referent HPV16, the E1-1374^63nt E-G350 variant was significantly associated with lower grade cervical lesions (p = 0.029), while the E1-1374^63nt E-C109/G350 variant was equally distributed between high and low grade lesions. The E1-1374^63nt variants were phylogenetically closest to E-G350 variant lineage (A2 sub-lineage based on full genome classification). The major differences between E1-1374^63nt variants were within the LCR and the E6 region. On the other hand, changes within the E1 region were the major differences from the A2 sub-lineage, which has been historically but inconclusively associated with high grade cervical disease. Thus, the shared variations cannot explain the particular association of the E1-1374^63nt variant with lower grade cervical lesions.

CONCLUSIONS/SIGNIFICANCE

The E1 region has been thus far considered to be well conserved among all HPVs and therefore uninteresting for variability studies. However, this study shows that the variations within the E1 region could possibly affect cervical disease, since the E1-1374^63nt E-G350 variant is significantly associated with lower grade cervical lesions, in comparison to the A1 and A2 sub-lineage variants. Furthermore, it appears that the silent variation 109T>C of the E-C109/G350 variant might have a significant role in the viral life cycle and warrants further study.

Human papillomavirus type 16 genetic variants: phylogeny and classification based on E6 and LCR

Naturally occurring genetic variants of human papillomavirus type 16 (HPV16) are common and have previously been classified into 4 major lineages; European-Asian (EAS), including the sublineages European (EUR) and Asian (As), African 1 (AFR1), African 2 (AFR2), and North-American/Asian-American (NA/AA). We aimed to improve the classification of HPV16 variant lineages by using a large resource of HPV16-positive cervical samples collected from geographically diverse populations in studies on HPV and/or cervical cancer undertaken by the International Agency for Research on Cancer. In total, we sequenced the entire E6 genes and long control regions (LCRs) of 953 HPV16 isolates from 27 different countries worldwide. Phylogenetic analyses confirmed previously described variant lineages and subclassifications. We characterized two new sublineages within each of the lineages AFR1 and AFR2 that are robustly classified using E6 and/or the LCR. We could differentiate previously identified AA1, AA2, and NA sublineages, although they could not be distinguished by E6 alone, requiring the LCR for correct phylogenetic classification. We thus provide a classification system for HPV16 genomes based on 13 and 32 phylogenetically distinguishing positions in E6 and the LCR, respectively, that distinguish nine HPV16 variant sublineages (EUR, As, AFR1a, AFR1b, AFR2a, AFR2b, NA, AA1, and AA2). Ninety-seven percent of all 953 samples fitted this classification perfectly. Other positions were frequently polymorphic within one or more lineages but did not define phylogenetic subgroups. Such a standardized classification of HPV16 variants is important for future epidemiological and biological studies of the carcinogenic potential of HPV16 variant lineages.

E6 molecular variants of human papillomavirus (HPV) type 16: an updated and unified criterion for clustering and nomenclature

Reports on taxonomic identification of E6/HPV 16 variants, don't have a worldwide, updated and unified criterion for clustering and nomenclature. Our aim was to update the existing criterion and propose a new one for clustering and nomenclature for E6/HPV 16 molecular variants based on the descriptive and comparative analyses of nucleotide sequences. A systematic search of the publications between 1991 and 2010 was carried out in PUBMED and manually. 240 E6/HPV 16 variants were identified. 157 were classified as European (E), 24 as Asian (As), 14 as Asian American (AA), 11 as North American 1 (NA 1), 19 as African 1 (Af 1) and 15 as African 2 (Af 2). Three classes were determined for the E, 3 each for the As, Af 2 and AA branches, 4 classes for the NA 1 and 6 for the Af 1 branch. This study reports our results and proposes an updated criterion for clustering and nomenclature that will be useful for E6 variant identification.

Association of viral load with HPV16 positive cervical cancer pathogenesis: causal relevance in isolates harboring intact viral E2 gene

We tested the hypothesis that cervical cancers (CaCx) harbor high HPV16 viral load compared to controls and this is influenced by E2 status and age of subjects. Viral load (natural log transformed values) per 100ng genomic DNA was estimated (152 cases and 87 controls) by Taqman assay. Median viral load was significantly higher (Mann-Whitney U test) among cases (17.21) compared to controls (9.86), irrespective of E2 status or upon considering E2 status as a covariate in logistic regression model (p<0.001). Viral load of E2 intact cases (17.80) was significantly higher (p<0.001) compared to those with disrupted E2 (9.78). At equivalent probability of being a case, viral load was higher among individuals (i) of lower age, irrespective of E2 status, and (ii) with intact E2 but of similar age as those with disrupted E2. Thus viral load in association with E2 status and/or age might be of causal relevance in CaCx pathogenesis.