Intron Definition and a Branch Site Adenosine at nt 385 Control RNA Splicing of HPV16 E6*I and E7 Expression

HPV16 E6 and E7, two viral oncogenes, are expressed from a single bicistronic pre-mRNA. In this report, we provide the evidence that the bicistronic pre-mRNA intron 1 contains three 5′ splice sites (5′ ss) and three 3′ splice sites (3′ ss) normally used in HPV16⁺ cervical cancer and its derived cell lines. The choice of two novel alternative 5′ ss (nt 221 5′ ss and nt 191 5′ ss) produces two novel isoforms of E6E7 mRNAs (E6*V and E6*VI). The nt 226 5′ ss and nt 409 3′ ss is preferentially selected over the other splice sites crossing over the intron to excise a minimal length of the intron in RNA splicing. We identified AACAAAC as the preferred branch point sequence (BPS) and an adenosine at nt 385 (underlined) in the BPS as a branch site to dictate the selection of the nt 409 3′ ss for E6*I splicing and E7 expression. Introduction of point mutations into the mapped BPS led to reduced U2 binding to the BPS and thereby inhibition of the second step of E6E7 splicing at the nt 409 3′ ss. Importantly, the E6E7 bicistronic RNA with a mutant BPS and inefficient splicing makes little or no E7 and the resulted E6 with mutations of ⁹¹QYNK⁹⁴ to ⁹¹PSFW⁹⁴ displays attenuate activity on p53 degradation. Together, our data provide structural basis of the E6E7 intron 1 for better understanding of how viral E6 and E7 expression is regulated by alternative RNA splicing. This study elucidates for the first time a mapped branch point in HPV16 genome involved in viral oncogene expression.     

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.     

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.     

Sequence duplication and internal deletion in the integrated human papillomavirus type 16 genome cloned from a cervical carcinoma.

Integrated human papillomavirus type 16 (HPV16) sequences were cloned from a cervical carcinoma and analyzed by restriction mapping and nucleotide sequencing. The viral integration sites were mapped within the E1 and E2 open reading frames (ORFs). The E4 and E5 ORFs were entirely deleted. An internal deletion of 376 base pairs (bp) was found disrupting the L1 and L2 ORFs. Sequencing analysis showed that an AGATGT/ACATCT inverted repeat marked the deletion junction with two flanking direct repeats 14 and 8 bp in length. A 1,330-bp sequence duplication containing the long control region (LCR) and the E6 and E7 ORFs was also found. The duplication junction was formed by two 24-bp direct repeats with 79% (19 of 24) homology located within the LCR and the E2 ORF of the prototype viral genome, respectively. This observation leads us to propose that the initial viral integration involved an HPV16 dimer in which the direct repeats in tandem units recombined, resulting in reiteration of only a portion of the original duplication. A guanosine insertion between nucleotides 1137 and 1138 created a continuous E1 ORF which was previously shown to be disrupted. Results from this study indicate that sequence reiteration and internal deletion in the integrated, and possibly in the episomal, HPV16 genome are influenced by specific nucleotide sequences in the viral genome. Moreover, reiteration of the LCR/E6/E7 sequences further supports the hypothesis that the E6/E7 ORFs may code for oncogenic proteins and that regulatory signals in the LCR may play a role in cellular transformation.     

Human Papillomavirus Type 6 and 11 Genetic Variants Found in 71 Oral and Anogenital Epithelial Samples from Australia

Genetic variation of 49 human papillomavirus (HPV) 6 and 22 HPV11 isolates from recurrent respiratory papillomatosis (RRP) (n = 17), genital warts (n = 43), anal cancer (n = 6) and cervical neoplasia cells (n = 5), was determined by sequencing the long control region (LCR) and the E6 and E7 genes. Comparative analysis of genetic variability was examined to determine whether different disease states resulting from HPV6 or HPV11 infection cluster into distinct variant groups. Sequence variation analysis of HPV6 revealed that isolates cluster into variants within previously described HPV6 lineages, with the majority (65%) clustering to HPV6 sublineage B1 across the three genomic regions examined. Overall 72 HPV6 and 25 HPV11 single nucleotide variations, insertions and deletions were observed within samples examined. In addition, missense alterations were observed in the E6/E7 genes for 6 HPV6 and 5 HPV11 variants. No nucleotide variations were identified in any isolates at the four E2 binding sites for HPV6 or HPV11, nor were any isolates found to be identical to the HPV6 lineage A or HPV11 sublineage A1 reference genomes. Overall, a high degree of sequence conservation was observed between isolates across each of the regions investigated for both HPV6 and HPV11. Genetic variants identified a slight association with HPV6 and anogenital lesions (p = 0.04). This study provides important information on the genetic diversity of circulating HPV 6 and HPV11 variants within the Australian population and supports the observation that the majority of HPV6 isolates cluster to the HPV6 sublineage B1 with anogenital lesions demonstrating an association with this sublineage (p = 0.02). Comparative analysis of Australian isolates for both HPV6 and HPV11 to those from other geographical regions based on the LCR revealed a high degree of sequence similarity throughout the world, confirming previous observations that there are no geographically specific variants for these HPV types.     

Identification of a Novel Human Papillomavirus,Type HPV199, Isolated from a Nasopharynx and Anal Canal,and Complete Genomic Characterization of Papillomavirus Species Gamma-12

The novel human papillomavirus type 199 (HPV199) was initially identified in a nasopharyngeal swab sample obtained from a 25 year-old immunocompetent male. The complete genome of HPV199 is 7,184 bp in length with a GC content of 36.5%. Comparative genomic characterization of HPV199 and its closest relatives showed the classical genomic organization of Gammapapillomaviruses (Gamma-PVs). HPV199 has seven major open reading frames (ORFs), encoding five early (E1, E2, E4, E6, and E7) and two late (L1 and L2) proteins, while lacking the E5 ORF. The long control region (LCR) of 513 bp is located between the L1 and E6 ORFs. Phylogenetic analysis additionally confirmed that HPV-199 clusters into the Gamma-PV genus, species Gamma-12, additionally containing HPV127, HV132, HPV148, HPV165, and three putative HPV types: KC5, CG2 and CG3. HPV199 is most closely related to HPV127 (nucleotide identity 77%). The complete viral genome sequence of additional HPV199 isolate was determined from anal canal swab sample. Two HPV199 complete viral sequences exhibit 99.4% nucleotide identity. To the best of our knowledge, this is the first member of Gamma-PV with complete nucleotide sequences determined from two independent clinical samples. To evaluate the tissue tropism of the novel HPV type, 916 clinical samples were tested using HPV199 type-specific real-time PCR: HPV199 was detected in 2/76 tissue samples of histologically confirmed common warts, 2/108 samples of eyebrow hair follicles, 2/137 anal canal swabs obtained from individuals with clinically evident anal pathology, 4/184 nasopharyngeal swabs and 3/411 cervical swabs obtained from women with normal cervical cytology. Although HPV199 was found in 1.4% of cutaneous and mucosal samples only, it exhibits dual tissue tropism. According to the results of our study and literature data, dual tropism of all Gamma-12 members is highly possible.     

Genetic Variation of Human Papillomavirus Type 16 in Individual Clinical Specimens Revealed by Deep Sequencing

Viral genetic diversity within infected cells or tissues, called viral quasispecies, has been mostly studied for RNA viruses, but has also been described among DNA viruses, including human papillomavirus type 16 (HPV16) present in cervical precancerous lesions. However, the extent of HPV genetic variation in cervical specimens, and its involvement in HPV-induced carcinogenesis, remains unclear. Here, we employ deep sequencing to comprehensively analyze genetic variation in the HPV16 genome isolated from individual clinical specimens. Through overlapping full-circle PCR, approximately 8-kb DNA fragments covering the whole HPV16 genome were amplified from HPV16-positive cervical exfoliated cells collected from patients with either low-grade squamous intraepithelial lesion (LSIL) or invasive cervical cancer (ICC). Deep sequencing of the amplified HPV16 DNA enabled de novo assembly of the full-length HPV16 genome sequence for each of 7 specimens (5 LSIL and 2 ICC samples). Subsequent alignment of read sequences to the assembled HPV16 sequence revealed that 2 LSILs and 1 ICC contained nucleotide variations within E6, E1 and the non-coding region between E5 and L2 with mutation frequencies of 0.60% to 5.42%. In transient replication assays, a novel E1 mutant found in ICC, E1 Q381E, showed reduced ability to support HPV16 origin-dependent replication. In addition, partially deleted E2 genes were detected in 1 LSIL sample in a mixed state with the intact E2 gene. Thus, the methods used in this study provide a fundamental framework for investigating the influence of HPV somatic genetic variation on cervical carcinogenesis.     

Evolution and Taxonomic Classification of Alphapapillomavirus 7 Complete Genomes: HPV18, HPV39, HPV45, HPV59, HPV68 and HPV70

Background

The species Alphapapillomavirus 7 (alpha-7) contains human papillomavirus genotypes that account for 15% of invasive cervical cancers and are disproportionately associated with adenocarcinoma of the cervix. Complete genome analyses enable identification and nomenclature of variant lineages and sublineages.

Methods

The URR/E6 region was sequenced to screen for novel variants of HPV18, 39, 45, 59, 68, 70, 85 and 97 from 1147 cervical samples obtained from multiple geographic regions that had previously been shown to contain an alpha-7 HPV isolate. To study viral heterogeneity, the complete 8 kb genome of 128 isolates, including 109 sequenced for this analysis, were annotated and analyzed. Viral evolution was characterized by constructing phylogenic trees using maximum-likelihood and Bayesian algorithms. Global and pairwise alignments were used to calculate total and ORF/region nucleotide differences; lineages and sublineages were assigned using an alphanumeric system. The prototype genome was assigned to the A lineage or A1 sublineage.

Results

The genomic diversity of alpha-7 HPV types ranged from 1.1% to 6.7% nucleotide sequence differences; the extent of genome-genome pairwise intratype heterogeneity was 1.1% for HPV39, 1.3% for HPV59, 1.5% for HPV45, 1.6% for HPV70, 2.1% for HPV18, and 6.7% for HPV68. ME180 (previously a subtype of HPV68) was designated as the representative genome for HPV68 sublineage C1. Each ORF/region differed in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) / noncoding region 2 (NCR2) > upstream regulatory region (URR) > E6 / E7 > E2 / L2 > E1 / L1.

Conclusions

These data provide estimates of the maximum viral genomic heterogeneity of alpha-7 HPV type variants. The proposed taxonomic system facilitates the comparison of variants across epidemiological and molecular studies. Sequence diversity, geographic distribution and phylogenetic topology of this clinically important group of HPVs suggest an independent evolutionary history for each type.     

Establishment of an efficient multiplex real‐time PCR assay for human papillomavirus genotyping in cervical cytology specimens: comparison with hybrid capture II

J.‐H. Lee, N.‐W. Lee, S.‐W. Hong, Y.‐S. Nam, J.‐W. Choi and Y.‐S. Kim Establishment of an efficient multiplex real‐time PCR assay for human papillomavirus genotyping in cervical cytology specimens: comparison with hybrid capture II Objective: To establish an efficient multiplex real‐time PCR assay for 15 human papillomavirus (HPV) genotypes, we designed multiplexing parameters and compared our PCR system with the hybrid capture (HC) II test using cervical cytology samples. Methods: For preventing cross‐reactive amplifications, variable HPV genes (E1, E2, E6, E7 and L1) were targeted. The melting temperatures of all primers and probes, and the size of the PCR product were optimized for the multiplex PCR. Our PCR system was compared with the HC II assays in the detection and genotyping of HPV infection using 173 cytology smears. Discordant cases between the two assays were verified by direct HPV DNA sequencing. Results: Of 173 women, 93 (53.8%) were HPV‐positive by the HC II assay and/or the multiplex real‐time PCR assay. The HPV genotypes were determined in 92 (98.9%) of 93 cases by the multiplex real‐time PCR and/or DNA sequencing. The agreement rate between multiplex PCR and HC II methods was 91.9% (kappa = 0.84). Although the sample size of this study needs to be increased to have epidemiological significance, multiple infections and HPV 16 were the predominant type. HPV 58, 52 and 18 accounted for 25% of HPV infections. HPV 52, 58 and 31 constituted 30% of CIN 2/3. Conclusion: The multiplex real‐time PCR system shows a good and reliable clinical performance. This in house PCR assay is fast and cost‐effective for HPV genotyping and the detection of HPV co‐infection in the post‐HPV vaccination era.     

Tumourigenesis Driven by the Human Papillomavirus Type 16 Asian-American E6 Variant in a Three-Dimensional Keratinocyte Model

Infection with a transforming human papillomavirus (HPV) such as type 16 (of species Alphapapillomavirus 9) causes ano-genital and oral tumours via viral persistence in human squamous cell epithelia. Epidemiological studies showed that the naturally occurring HPV16 Asian-American (AA) variant (sublineage D2/D3) is found more often than the European Prototype (EP) (sublineage A1) in high-grade cervical neoplasia and tumours compared to non-cancer controls. Just three amino acid changes within the early gene, E6, of HPV16 AA have been linked to this augmented tumourigenicity. The AAE6 variant''s greater immortalizing and transforming potential over EPE6 has recently been confirmed in retrovirally-transduced keratinocytes expressing the E6 gene only. However, the tumourigenic role of the full-length viral genome of HPV16 has not yet been addressed with regard to these E6 variants. To investigate this process in the context of these two HPV16 E6 genotypes, an organotypic tissue culture model was used to simulate the HPV infectious life cycle. The AAE6 variant demonstrated an enhanced ability over EPE6 to drive the viral life cycle toward tumourigenesis, as evidenced phenotypically—by a more severe grade of epithelial dysplasia with higher proliferation and deregulated differentiation, and molecularly—by high viral oncogene E6 and E7 expression, but lack of productive viral life cycle markers. In contrast, EPE6 had low E6 and E7 but high E1∧E4 expression, indicative of a productive life cycle. We suggest increased viral integration into the host genome for AAE6 as one possible mechanism for these observed differences from EPE6. Additionally, we found downstream effects on immortalization and host innate immune evasion. This study highlights how minor genomic variations in transforming viruses can have a significant affect on their tumourigenic ability.     

Oncogenic comparison of human papillomavirus type 58 E7 variants

Human papillomavirus 58 (HPV58) ranks the second or third in East Asian cervical cancers. Current studies on HPV58 are scarce and focus on the prototype. Previously, we identified the three most common circulating HPV58 E7 strains contained amino acid alterations: G41R/G63D (51%), T20I/G63S (22%) and T74A/D76E (14%) respectively. Among them, the T20I/G63S variant (V1) had a stronger epidemiological association with cervical cancer. We therefore suggested that V1 possessed stronger oncogenicity than the other two variants. Here, we performed phenotypic assays to characterize and compare their oncogenicities with HPV58 E7 prototype. Our results showed that overexpression of V1 conferred a higher colony‐forming ability to primary murine epithelial cells than prototype (P < 0.05) and other variants, implicating its higher immortalising potential. Further experiments showed that both V1 and prototype enhanced the anchorage‐independent growth of NIH/3T3 cells (P < 0.001), implicating their stronger transforming power than the two other variants. Moreover, they possessed an increased ability to degrade pRb (P < 0.001), which is a major effector pathway of E7‐driven oncogenesis. Our work represents the first study to compare the oncogenicities of HPV58 E7 prototype and variants. These findings deepened our understanding of HPV58 and might inform clinical screening and follow‐up strategy.     

Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses

Human papillomavirus type 16 (HPV16) is the primary etiological agent of cervical cancer, the second most common cancer in women worldwide. Complete genomes of 12 isolates representing the major lineages of HPV16 were cloned and sequenced from cervicovaginal cells. The sequence variations within the open reading frames (ORFs) and noncoding regions were identified and compared with the HPV16R reference sequence. This whole-genome approach gives us unprecedented precision in detailing sequence-level changes that are under selection on a whole-viral-genome scale. Of 7,908 base pair nucleotide positions, 313 (4.0%) were variable. Within the 2,452 amino acids (aa) comprising 8 ORFs, 243 (9.9%) amino acid positions were variable. In order to investigate the molecular evolution of HPV16 variants, maximum likelihood models of codon substitution were used to identify lineages and amino acid sites under selective pressure. Five codon sites in the E5 (aa 48, 65) and E6 (aa 10, 14, 83) ORFs were demonstrated to be under diversifying selective pressure. The E5 ORF had the overall highest nonsynonymous/synonymous substitution rate (omega) ratio (M3 = 0.7965). The E2 gene had the next-highest omega ratio (M3 = 0.5611); however, no specific codons were under positive selection. These data indicate that the E6 and E5 ORFs are evolving under positive Darwinian selection and have done so in a relatively short time period. Whether response to selective pressure upon the E5 and E6 ORFs contributes to the biological success of HPV16, its specific biological niche, and/or its oncogenic potential remains to be established.     

Evolution and taxonomic classification of human papillomavirus 16 (HPV16)-related variant genomes: HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67

Background

Human papillomavirus 16 (HPV16) species group (alpha-9) of the Alphapapillomavirus genus contains HPV16, HPV31, HPV33, HPV35, HPV52, HPV58 and HPV67. These HPVs account for 75% of invasive cervical cancers worldwide. Viral variants of these HPVs differ in evolutionary history and pathogenicity. Moreover, a comprehensive nomenclature system for HPV variants is lacking, limiting comparisons between studies.

Methods

DNA from cervical samples previously characterized for HPV type were obtained from multiple geographic regions to screen for novel variants. The complete 8 kb genomes of 120 variants representing the major and minor lineages of the HPV16-related alpha-9 HPV types were sequenced to capture maximum viral heterogeneity. Viral evolution was characterized by constructing phylogenic trees based on complete genomes using multiple algorithms. Maximal and viral region specific divergence was calculated by global and pairwise alignments. Variant lineages were classified and named using an alphanumeric system; the prototype genome was assigned to the A lineage for all types.

Results

The range of genome-genome sequence heterogeneity varied from 0.6% for HPV35 to 2.2% for HPV52 and included 1.4% for HPV31, 1.1% for HPV33, 1.7% for HPV58 and 1.1% for HPV67. Nucleotide differences of approximately 1.0% - 10.0% and 0.5%–1.0% of the complete genomes were used to define variant lineages and sublineages, respectively. Each gene/region differs in sequence diversity, from most variable to least variable: noncoding region 1 (NCR1) /noncoding region 2 (NCR2) >upstream regulatory region (URR)> E6/E7 > E2/L2 > E1/L1.

Conclusions

These data define maximum viral genomic heterogeneity of HPV16-related alpha-9 HPV variants. The proposed nomenclature system facilitates the comparison of variants across epidemiological studies. Sequence diversity and phylogenies of this clinically important group of HPVs provides the basis for further studies of discrete viral evolution, epidemiology, pathogenesis and preventative/therapeutic interventions.     

Function of a Bovine Papillomavirus Type 1 Exonic Splicing Suppressor Requires a Suboptimal Upstream 3′ Splice Site

Alternative splicing is an important mechanism for the regulation of bovine papillomavirus type 1 (BPV-1) gene expression during the virus life cycle. Previous studies in our laboratory have identified two purine-rich exonic splicing enhancers (ESEs), SE1 and SE2, located between two alternative 3′ splice sites at nucleotide (nt) 3225 and nt 3605. Further analysis of BPV-1 late-pre-mRNA splicing in vitro revealed a 48-nt pyrimidine-rich region immediately downstream of SE1 that inhibits utilization of the nt 3225 3′ splice site. This inhibitory element, which we named an exonic splicing suppressor (ESS), has a U-rich 5′ end, a C-rich central part, and an AG-rich 3′ end (Z. M. Zheng, P. He, and C. C. Baker, J. Virol. 70:4691–4699, 1996). The present study utilized in vitro splicing of both homologous and heterologous pre-mRNAs to further characterize the ESS. The BPV-1 ESS was inserted downstream of the 3′ splice site in the BPV-1 late pre-mRNA, Rous sarcoma virus src pre-mRNA, human immunodeficiency virus tat-rev pre-mRNA, and Drosophila dsx pre-mRNA, all containing a suboptimal 3′ splice site, and in the human β-globin pre-mRNA, which contains a constitutive 3′ splice site. These studies demonstrated that suppression of splicing by the BPV-1 ESS requires an upstream suboptimal 3′ splice site but not an upstream ESE. Furthermore, the ESS functions when located either upstream or downstream of BPV-1 SE1. Mutational analyses demonstrated that the function of the ESS is sequence dependent and that only the C-rich region of the ESS is essential for suppression of splicing in all the pre-mRNAs tested.