Evolutionary dynamics of variant genomes of human papillomavirus types 18, 45, and 97

Vaccine and Non-Vaccine HPV Types Presence in Adolescents with Vertically Acquired HIV Five Years Post Gardasil Quadrivalent Vaccination: The ZIMGARD Cohort

BACKGROUND

Human papillomavirus (HPV) vaccination programs are a key intervention in protecting individuals against HPV-related disease. HIV1-infected individuals are at increased risk of HPV-associated cancers. This study was conducted to evaluate the potential role of prophylactic HPV vaccines in preventing new HPV infections among participants with perinatally acquired HIV who received the quadrivalent HPV vaccine at least five years before this study.

METHODS

This cross-sectional study was conducted at Newlands Clinic, Harare, Zimbabwe. The clinic provided the Gardasil quadrivalent HPV vaccine (4vHPV) to 624 adolescents living with HIV starting in December 2015. Vaginal and penile swabs were collected and tested for HPV types from the study participants who had received the 4vHPV vaccine 5-6 years before enrolment.

RESULTS

We present the results of 98 participants (44.6% female) vaccinated at a median age of 15 years (IQR 12-16). The mean amount of time since vaccination was 6 years (SD: ±0.4). The HPV-positive rate amongst the analyzed swabs was 69% (68/98). Among 30/98 (31%) HPV-positive participants, 13/98 (13%) had low-risk HPV types, and 17/98 (17%) had high-risk HPV types. Twelve participants tested positive for HPV18, only one participant tested positive for HPV16, and an additional four (4.3%) tested positive for either type 6 or 11, with respect to vaccine-preventable low-risk HPV types.

CONCLUSION

The Gardasil quadrivalent HPV vaccine (4vHPV) was expected to protect against infection with HPV types 16, 18, 6, and 11. We demonstrated a possible waning of immunity to HPV18 in 17% of the participants, and an associated loss in cross-protection against HPV45. We observed a relatively high prevalence of 'opportunistic non-vaccine HPV types' or 'ecological niche occupiers' in this cohort, and suggest further research on the involvement of these types in cervical and other genital cancers. Our study is one of the few, if not the first, to report on HPV vaccine immunoprotection among people living with HIV (PLWH), thereby setting a baseline for further studies on HPV vaccine effectiveness among PLWH.

Type distribution of human papillomaviruses in ThinPrep cytology samples and HPV16/18 E6 gene variations in FFPE cervical cancer specimens in Fars province, Iran

BACKGROUND

There exists strong evidence that human papillomavirus (HPV) is associated with cervical cancer (CC). HPV E6 is a major oncogene whose sequence variations may be associated with the development of CC. There is not sufficient data on the distribution of HPV types in ThinPrep cytology specimens and HPV 16/18 E6 gene variations among CC patients in the southwest of Iran. This study was conducted to contribute to HPV screening and vaccination in Iran.

METHODS

A total of 648 women screened for cervicitis, intraepithelial neoplasia or CC were included in the study. All participants underwent ThinPrep cytology testing, single-step HPV DNA detection and allele-specific reverse hybridization assays. Moreover, a total of 96 specimens previously tested positive for single infection with HPV16 or 18 were included for variant analysis. HPV16/18 lineages and sublineages were determined by PCR assays followed by sequencing the E6 gene and the construction of neighbor-joining phylogenetic trees.

RESULTS

Overall, HPV DNA was detected in 62.19% of all the screened subjects. The detection rates of HPV DNA among individuals with normal, ASC-US, ASC-H, LSIL, and HSIL cervical cytology were 48.9%, 93.6%, 100%, 100%, and 100%, respectively. Low-risk HPVs were detected more frequently (46.9%) than high-risk (38.9%) and possible high-risk types (11.1%). Of 403 HPV-positive subjects, 172 (42.7%) had single HPV infections while the remaining 231 (57.3%) were infected with multiple types of HPV. Our results indicated a remarkable growth of high-risk HPV66 and 68 and low-risk HPV81 which have rarely been reported in Iran and HPV90 and 87 that are reported for the first time in the country. In addition, 3 lineages (A, D, and C) and 6 sublineages (A1, A2, A4, C1, D1, and D2) of HPV16, and one lineage and 4 sublineages (A1, A3, A4, and A5) of HPV18 were identified. The studied HPV16 and 18 variants mainly belonged to the D1 and A4 sublineages, respectively.

CONCLUSION

The present study suggests that the prevalence of HPV infection in women of all age groups with or without premalignant lesions in the southwestern Iran is high and the predominant HPV types in the southwest of Iran may differ from those detected in other parts of the country. This study also highlights the necessity of not only initiating HPV vaccination for the general population but also developing new vaccines that confer immunity against the prevalent HPV types in the area and national cervical screening programs using a combination of thinPrep cytology test and HPV detection assays in order to improve the accuracy of the screening.

Human papillomavirus genomics: Understanding carcinogenicity

Human papillomavirus (HPV) causes virtually all cervical cancers and many cancers at other anatomical sites in both men and women. However, only 12 of 448 known HPV types are currently classified as carcinogens, and even the most carcinogenic type - HPV16 - only rarely leads to cancer. HPV is therefore necessary but insufficient for cervical cancer, with other contributing factors including host and viral genetics. Over the last decade, HPV whole genome sequencing has established that even fine-scale within-type HPV variation influences precancer/cancer risks, and that these risks vary by histology and host race/ethnicity. In this review, we place these findings in the context of the HPV life cycle and evolution at various levels of viral diversity: between-type, within-type, and within-host. We also discuss key concepts necessary for interpreting HPV genomic data, including features of the viral genome; events leading to carcinogenesis; the role of APOBEC3 in HPV infection and evolution; and methodologies that use deep (high-coverage) sequencing to characterize within-host variation, as opposed to relying on a single representative (consensus) sequence. Given the continued high burden of HPV-associated cancers, understanding HPV carcinogenicity remains important for better understanding, preventing, and treating cancers attributable to infection.

Determination of Human Papillomavirus Type 18 Lineage of E6: A Population Study from Iran

The epidemiological studies in Iran on HPV18 nucleotide changes are rare. This type of virus is prevalent in the Iranian population. Therefore, in the present study, we aimed to identify the genetic variability in HPV18 in the E6 region to evaluate the prevalence of lineage distribution and sublineages in a sample population in Iran. Overall, 60 HPV18 confirmed cases were investigated between 2019 and 2021. The specimens were collected, and molecular genotyping was done using the Linear Array HPV Genotyping Test. DNA extraction was performed by a viral DNA/RNA kit. The HPV E6 gene was amplified by using type-specific primers designed according to the HPV18 genome prototype sequence. The sequencing of the E6 region was successfully done on 43 samples which were then compared to the reference sequence. The most frequent sublineage of HPV18 in this study was A4 (69.7%), followed by A1 (18.6%) and A3 (11.6%). Neither A2 nor A5 sublineage was not detected in this study. The related nucleotide acid changes according to the main references were as follows: A3: T104C/T232G/T485C/C549A, A4: T104C/T485C/C549A. The predominance of A lineage with the high frequency of A4 sublineage was found in the present research. The importance of sublineages in susceptibility to a progressive form of infection requires to be more investigated among the different population.

Genetic variability of the HPV16 early genes and LCR. Present and future perspectives

Human papillomavirus 16 (HPV16) infection is the aetiologic factor for the development of cervical dysplasia and is regarded as highly carcinogen, because it is implicated in more than 50% of cervical cancer cases, worldwide. The tumourigenic potential of HPV16 has triggered the extensive sequence analysis of viral genome in order to identify nucleotide variations and amino acid substitutions that influence viral oncogenicity and subsequently the initiation and progression of cervical cancer. Nowadays, specific mutations of HPV16 DNA have been associated with an increased risk of high-grade squamous intraepithelial lesions and invasive cervical cancer (ICC) development, including E6: Q14H, H78Y, L83V, Ε7: N29S, S63F, E2: H35Q, P219S, T310K, E5: I65V, whereas highly conserved regions of viral DNA have been extensively characterised. In addition, numerous novel HPV16 mutations are observed among the studied populations from various geographic regions, hence advocating that different HPV16 strains seem to emerge with different tumourigenic capacities. The present review focuses on the variability of the early genes and the long control region, emphasising on the association of specific mutations with the development of severe dysplasia. Finally, it evaluates whether specific regions of HPV16 DNA are able to serve as valuable biomarkers for cervical cancer risk.

Molecular Characterization of Human Papillomavirus Type 159 (HPV159)

Human papillomavirus type 159 (HPV159) was identified in an anal swab sample and preliminarily genetically characterized by our group in 2012. Here we present a detailed molecular in silico analysis that showed that the HPV159 viral genome is 7443 bp in length and divided into five early and two late genes, with conserved functional domains and motifs, and a non-coding long control region (LCR) with significant regulatory sequences that allow the virus to complete its life cycle and infect novel host cells. HPV159, clustering into the cutaneotropic Betapapillomavirus (Beta-PV) genus, is phylogenetically most similar to HPV9, forming an individual phylogenetic group in the viral species Beta-2. After testing a large representative collection of clinical samples with HPV159 type-specific RT-PCR, in addition to the anal canal from which the first HPV159 isolate was obtained, HPV159 was further detected in other muco-cutaneous (4/181, 2.2%), mucosal (22/764, 2.9%), and cutaneous (14/554, 2.5%) clinical samples, suggesting its extensive tissue tropism. However, because very low HPV159 viral loads were estimated in the majority of positive samples, it seemed that HPV159 mainly caused clinically insignificant infections of the skin and mucosa. Using newly developed, highly sensitive HPV159-specific nested PCRs, two additional HPV159 LCR viral variants were identified. Nevertheless, all HPV159 mutations were demonstrated outside important functional domains of the LCR, suggesting that the HPV159 viral variants were most probably not pathogenically different. This complete molecular characterization of HPV159 enhances our knowledge of the genome characteristics, tissue tropism, and phylogenetic diversity of Beta-PVs that infect humans.

Phylogeny and In Silico Structure Analysis of Major Capsid Protein (L1) Human Papillomavirus 45 from Indonesian Isolates

BACKGROUND

Human papillomavirus (HPV)-45 genotype circulates in high percentage in Bandung area - Indonesia, after HPV-16 and HPV-18. The aim of this study was to analyse variations of major capsid (L1) HPV-45 and its phylogeny. Furthermore in silico protein structure and epitope prediction was explored.

METHODS

L1 gene of HPV-45 was amplified, sequenced and aligned. Phylogenetic tree had been built and compared with a complete L1 HPV-45 sequence. Structure and epitope prediction of L1 protein were then developed in silico.

RESULTS

Of 5 L1 HPV-45 sequences collected, we have detected one variant of sub lineage A2 which was considered as a new variant, and two variants of B2. Superimposition of structure of these two variants with reference showed very similar structure. Furthermore, seven amino acid substitutions were found within these L1 variants of which two substitutions might change the polarity of corresponding amino acid I329T and S383G. The S383G occurred in surface loop (HI-Loop) of new L1 HPV-45 variant.

CONCLUSION

Similar structure of Indonesian variants indicates that amino acids variations do not affect the L1 structure. However, one substitution with altered amino acid polarity found within the area of surface loop suggests a potential impact in antibody recognition and neutralization.

Fifteen new nucleotide substitutions in variants of human papillomavirus 18 in Korea : Korean HPV18 variants and clinical manifestation

High-risk human papillomavirus (HPV) infection is an essential factor for the development of cervical cancer. HPV18 is the second most common carcinogenic HPV type following HPV16, but the lineages of HPV18 have been less well studied than those of HPV 16. The purpose of this study was to analyze the nucleotide variants in the E6, E7, and L1 genes of HPV18, to assess the prevalence of HPV18 variants in Korea and to explore the relationship between HPV18 genetic variants and the risk for cervical cancer.

A total of 170 DNA samples from HPV18-positive cervical specimens were collected from women admitted to a secondary referral hospital located in Seoul. Among them, the lineages of the 97 samples could be successfully determined by historical nomenclature.

All the studied HPV 18 variants were lineage A. Sublineages A1 and A4 comprised 91.7% (89/97) and 1.0% (1/97), respectively. Sublineages other than A1 or A4 comprised 7.2% (7/97). We identified 15 new nucleotide substitutions among 44 nucleotide substitutions: C158T, T317G, T443G, A560G, A5467G, A5560C, A5678C, A6155G, G6462A, T6650G, G6701A, T6809C, A6823G, T6941C and T6953C. Among them, 6 substitutions at positions 317, 443, 5467, 5560, 6462, and 6823 resulted in amino acid changes (E6: F71L and N113K; L1: H13R, H44P, A345T, and N465S, respectively). The pathologic results were classified as normal in 25.8% (25/97) of the women, atypical squamous cells of undermined significance (ASCUS) in 7.2% (7/97), cervical intraepithelial neoplasia (CIN) 1 in 36.1% (35/97), CIN2/3 in 19.6% (18/97), and carcinoma in 12.4% (12/97). There was no significant association between the HPV18 sublineages and the severity of pathologic lesion or the disease progression.

This study is the first to analyze the distribution of HPV18 variants in Korean and to associate the results with pathologic findings. Although the HPV18 variants had no significant effect on the degree and progression of the disease, the newly discovered nonsynonymous mutation in L1 might serve as a database to determine vaccine efficacy in Korean women.

Recombination Between High-Risk Human Papillomaviruses and Non-Human Primate Papillomaviruses: Evidence of Ancient Host Switching Among Alphapapillomaviruses

We use all the currently known 405 Papillomavirus (PV) sequences, 343 curated PV sequences from both humans and animals from the PAVE data base, to analyse the recombination dynamics of these viruses at the whole genome levels. After showing some evidence of human and non-human primate PV recombination, we report a comprehensive recombination analysis of all currently known 82 Alphapapillomaviruses (Alpha-PVs). We carried out an exploratory study and found novel recombination events between High-Risk HPV Types and Macaca fascicularis PV1 (MfPV1), Macaca Fuscata PV2 (MfuPV2) and Pan Paniscus PV1 (PpPV1) Papillomaviruses. This is the first evidence of interactions between PVs from different hosts and hence postulates the likelihood of ancient host switching among Alpha-PVs. Notwithstanding these results should be interpreted with caution because the major and minor parents indicated by RDP4 program are simply the sequences in the alignment that most closely resemble the actual parents. We found statistically significant differences between the phylogenies of the PV sequences with recombination regions and PV sequences without recombination regions using the Shimodaira–Hasegawa phylogenetic incongruence testing. We show that not more than 76MYA Alpha-PVs were in the same biological niche, a pre-requisite for recombination, and as the hosts evolved and diversified, the viruses adapted to specific host niches which eventually led to coevolution with specific hosts before speciation of primate species. Thus providing evidence that in ancient times no earlier than the Cretaceous period of the Mesozoic age, Alpha-PVs recombined and switched hosts, but whether this host switching is occurring currently is unknown. However, a clearer picture of the PVs evolutionary landscape can only be achieved with the incremental discovery of PV sequences, especially from the animal kingdom.

Molecular epidemiology of human papillomavirus 18 infections in Japanese Women

Among the oncogenic genotypes of human papillomavirus (HPV), HPV18 is the second most common type detected in cervical cancer worldwide and is primarily involved in the generation of cervical adenocarcinoma. Although HPV intra-type variants confer different risks of cervical carcinogenesis, there is little information on the genetic diversity of HPV18 compared to the most prevalent type, HPV16. In this study, we investigated the genetic variation of HPV18 in cervical specimens obtained from Japanese women with normal cervices or cervical cancers and precancers. Of the 101 HPV18-positive samples analyzed, viral whole genome amplification followed by next-generation sequencing led to the determination of viral complete genome sequences of 18 samples. Phylogenetic analysis of these HPV18 whole genome sequences identified a distinct variant cluster consisting of only Japanese samples (n = 7) belonging to sublineage A1. Viral genome sequences were also analyzed for the E6/E7 (n = 66) and E2 (n = 27) genes by Sanger sequencing. Phylogenetic analyses of these regions showed that the variant distribution among Japanese women was strongly biased toward sublineage A1 (72 of 87; 82.8%). No significant differences were observed in the prevalence of specific sublineages between cervical cancer/precancer cases and controls, and between squamous cell carcinoma and adenocarcinoma cases. These data contribute to our understanding of the genetic diversity of HPV18 in Japanese women.

Very Low Prevalence of Vaccine Human Papillomavirus Types Among 18- to 35-Year Old Australian Women 9 Years Following Implementation of Vaccination

Introduction

A quadrivalent human papillomavirus vaccination program targeting females aged 12-13 years commenced in Australia in 2007, with catch-up vaccination of 14-26 year olds through 2009. We evaluated the program's impact on HPV prevalence among women aged 18-35 in 2015.

Methods

HPV prevalence among women aged 18-24 and 25-35 was compared with prevalence in these age groups in 2005-2007. For women aged 18-24, we also compared prevalence with that in a postvaccine study conducted in 2010-2012.

Results

For the 2015 sample, Vaccination Register-confirmed 3-dose coverage was 53.3% (65.0% and 40.3% aged 18-24 and 25-35, respectively). Prevalence of vaccine HPV types decreased from 22.7% (2005-2007) and 7.3% (2010-2012), to 1.5% (2015) (P trend < .001) among women aged 18-24, and from 11.8% (2005-2007) to 1.1% (2015) (P = .001) among those aged 25-35.

Conclusions

This study, reporting the longest surveillance follow-up to date, shows prevalence of vaccine-targeted HPV types has continued to decline among young women. A substantial fall also occurred in women aged 25-35, despite lower coverage. Strong herd protection and effectiveness of less than 3 vaccine doses likely contributed to these reductions.

Evolutionary dynamics of ten novel Gamma-PVs: insights from phylogenetic incongruence, recombination and phylodynamic analyses

Background

Human papillomaviruses (HPVs) are genetically diverse, belonging to five distinct genera: Alpha, Beta, Gamma, Mu and Nu. All papillomaviruses have double stranded DNA genomes that are thought to evolve slowly because they co-opt high-fidelity host cellular DNA polymerases for their replication. Despite extensive efforts to catalogue all the HPV species that infect humans, it is likely that many still remain undiscovered. Here we use the sequences of ten novel Gammapapillomaviruses (Gamma-PVs) characterized in previous studies and related HPVs to analyse the evolutionary dynamics of these viruses at the whole genome and individual gene scales.

Results

We found statistically significant incongruences between the phylogenetic trees of different genes which imply gene-to-gene variation in the evolutionary processes underlying the diversification of Gamma-PVs. We were, however, only able to detect convincing evidence of a single recombination event which, on its own, cannot explain the observed incongruences between gene phylogenies. The divergence times of the last common ancestor (LCA) of the Alpha, Beta, Mu, Nu and Gamma genera was predicted to have existed between 49.7–58.5 million years ago, before splitting into the five main lineages. The LCA of the Gamma-PVs at this time was predicted to have existed between 45.3 and 67.5 million years ago: approximately at the time when the simian and tarsier lineages of the primates diverged.

Conclusion

Consequently, we report here phylogenetic tree incongruence without strong evidence of recombination.

Non-Vaccine-Type Human Papillomavirus Prevalence After Vaccine Introduction: No Evidence for Type Replacement but Evidence for Cross-Protection

BACKGROUND

We examined non-vaccine-type human papillomavirus (HPV) prevalence in a community before and during the first 8 years after vaccine introduction, to assess for (1) type replacement with any non-vaccine-type HPV and (2) cross-protection with non-vaccine types genetically related to vaccine-type HPV.

METHODS

Sexually experienced 13- to- 26-year-old women were recruited for 3 cross-sectional studies from 2006 to 2014 (N = 1180). Outcome variables were as follows: (1) prevalence of at least 1 of 32 anogenital non-vaccine-type HPVs and (2) prevalence of at least 1 HPV type genetically related to HPV-16 and HPV-18. We determined changes in proportions of non-vaccine-type HPV prevalence across the study waves using logistic regression with propensity score inverse probability weighting.

RESULTS

Vaccine initiation rates increased from 0% to 71.3%. Logistic regression demonstrated that from 2006 to 2014, there was no increase in non-vaccine-type HPV among vaccinated women (adjusted odds ratio [AOR], 1.02; 95% confidence interval [CI], 0.73-1.42), but an increase among unvaccinated women (AOR, 1.88; 95% CI, 1.16-3.04). Conversely, there was a decrease in types genetically related to HPV-16 among vaccinated (AOR, 0.57; 95% CI, 0.38-0.88) but not unvaccinated women (AOR, 1.33; 95% CI, 0.81-2.17).

CONCLUSIONS

We did not find evidence of type replacement, but did find evidence of cross-protection against types genetically related to HPV-16. These findings have implications for cost-effectiveness analyses, which may impact vaccine-related policies, and provide information to assess the differential risk for cervical cancer in unvaccinated and vaccinated women, which may influence clinical screening recommendations. The findings also have implications for public health programs, such as health messaging for adolescents, parents, and clinicians about HPV vaccination.

Discovery, characterisation and genomic variation of six novel Gammapapillomavirus types from penile swabs in South Africa

Six novel human papillomaviruses from penile swabs were characterised. Multiple full genome clones for each novel type were generated, and complete genome sizes were: HPV211 (7253bp), HPV212 (7208bp), HPV213 (7096bp), HPV214 (7357), HPV215 (7186bp) and HPV216 (7233bp). Phylogenetically the novel papillomaviruses all clustered with Gammapapillomaviruses: HPV211 is most closely related to HPV168 (72% identity in the L1 nucleotide sequence) of the Gamma-8 species, HPV212 is most closely related to HPV144 (82.9%) of the Gamma-17 species, HPV213 is most closely related to HPV153 (71.8%) of the Gamma-13 species, HPV214 is most closely related to HPV103 (75.3%) of the Gamma-6 species, HPV215 and HPV216 are most closely related to HPV129 (76.8% and 79.2% respectively) of the Gamma-9 species. The novel HPV types demonstrated the classical genomic organisation of Gammapapillomavirusess, with seven open reading frames (ORFs) encoding five early (E1, E2, E4, E6 and E7) and two late (L1 and L2) proteins. Typical of Gammapapillomavirusess the novel types all lacked the E5 ORF and HPV214 also lacked the E6 ORF. HPV212 had nine unique variants, HPV213 had five and HPV215 had four variants. Conserved domains observed among the novel types are the Zinc finger Binding Domain and PDZ domains. A retinoblastoma binding domain (pRB) binding domain in E7 protein was additionally identified in HPV214. This study expands the knowledge of the rapidly growing Gammapapillomavirus genus.

Phylogeny and Classification of Human Papillomavirus (HPV)16 and HPV18 Variants Based on E6 and L1 genes in Tunisian Women with Cervical Lesions

Background:

High-risk human papillomavirus (HPV) types are the main etiological factors for cervical cancer. HPV16 and HPV18 are generally the most common forms associated with development of high-grade cervical lesions. This study was undertaken to identify intratypic variants of HPV16 and HPV18 among women with cervical lesions in Tunisia.

Materials and Methods:

DNA was extracted from cervical samples collected from 49 women. using a PureLinkTM Genomic DNA mini Kit (Invitrogen). E6 and L1 open reading frames (ORF) were amplified by PCR and viral DNA amplicons were subjected to automated sequencing using Big Dye Terminators technology (Applied Biosystems). The obtained sequences were analyzed using an appropriate software program to allow phylogenetic trees to be generated.

Results:

HPV16 and HPV18 were detected in 15 and 5 cases, respectively. HPV16 E6 sequences clustered with the European German lineage (A2) whereas one isolate diverged differently in the L1 region and clustered with the African sub-lineage (B1). HPV 18 E6 sequences clustered with the European sub-lineage (A1) but L1 sequences clustered as a new clade which diverged from A1-A5.

Conclusions:

Our results suggest that the distribution of HPV16 and HPV18 sequences in women with cervical lesions in Tunisia is mainly related to European epidemiological conditions and point to the presence of recombinant HPV forms.

Phylogenetic analysis of Alphapapillomavirus based on L1, E6 and E7 regions suggests that carcinogenicity and tissue tropism have appeared multiple times during viral evolution

Members of the Alphapapillomavirus genus are causative agents for cervix cancer and benign lesions in humans. These viruses are classified according to sequence similarities in their L1 region. Yet, viral carcinogenicity has been associated with variations in the proteins encoded by the E6 and E7 genes. In order to relate evolutionary history with origin of carcinogenicity, we performed phylogenetic reconstructions using both nucleotide and predicted amino acid sequences of the L1, E6 and E7 genes. Whilst phylogenetic analysis of L1 reconstructed genus evolutionary history, phylogenies based on E6 and E7 proteins support the idea that mutations at amino acids S/Tx [V/L] (E6) and LxCxE (E7) might be responsible for carcinogenic potential. These findings indicate that virulence within Alphapapillomavirus have appeared multiple times during evolution. Our results reveal that oncogenic potential is not a monophyletic clade-specific adaptation but might be the result of positive selection on random mutations occurring on proteins involved in host infection during viral diversification.

Genetic diversity and phylogenetic analysis of HPV 16 & 18 variants isolated from cervical specimens of women in Saudi Arabia

Human papillomavirus (HPV) are well known to be associated with the development of cervical cancer. HPV16 and HPV 18 are known as high-risk types and reported to be predominantly associated with cervical cancer. The prevalence and genetic diversity of HPV have been well documented globally but, in the Kingdom of Saudi Arabia, data on HPV genetic diversity are lacking. In this study, we have analyzed the genetic diversity of both HPV16 and HPV18 based on their L1 gene sequence because L1 gene is a major capsid protein gene and has been utilized to develop a prophylactic vaccine. In January 2011-2012, a total of forty samples from cervical specimens of women in Saudi Arabia were collected. The association of HPV16, HPV18 was detected by polymerase chain reaction, sequenced and submitted to GenBank. The sequences identity matrix and the phylogenetic relationship were analyzed with selected HPVs. The highest sequence identity (99.5%) for HPV16 and (99.3%) for HPV was observed with selected HPVs. The phylogenetic analysis results showed that HPVs from Saudi Arabia formed a closed cluster with African, Asian, East Asian as well as American HPVs distributed into multiple linages from various geographical locations. The results provided the valuable information about genetic diversity, but there is an urgent need to generate full genome sequence information which will provide a clearer picture of the genetic diversity and evolution of HPVs in Saudi Arabia. In conclusion, the generated data will be highly beneficial for developing molecular diagnostic tools, analyzing and correlating the epidemiological data to determine the risk of cervical cancer and finally to develop a vaccine for Saudi Arabian population.

Human papillomavirus genome variants and head and neck cancers: a perspective

Human papillomaviruses (HPV) cause infections that are responsible for diverse clinical manifestations from benign conditions to invasive cancer. As different HPV types are associated with variable pathogenic potential, minor genetic variations within a given high-risk HPV type might also be associated with distinct oncogenic capacities, through variable ability of persistence or risk of progression to precancer/cancer. Most recent HPV variant studies in the cervix using latest sequencing technology confirmed that minor changes in the HPV genome can have a major influence on carcinogenesis and have revealed key data that help better understand the carcinogenicity of HPV at a molecular level. Here we review the limited number of studies on HPV genome variants in head and neck cancers (HNC) and discuss their implications for cancer research in the light of accumulated knowledge for the cervix. Challenges in transposing HPV variant studies from the lower anogenital to the upper aerodigestive tract are also discussed, highlighting the main gaps of knowledge in the field of HPV-induced HNC. Specifically in the head and neck region, the lack of characterisation of precancerous lesions and the difficulty in sampling normal tissue will challenge the development of accurate studies. Although there is so far no indication that HPV variant research in HNC could directly translate into clinical application, such research is expected to be useful to disentangle unanswered questions in the pathogenesis of HNC. Yet, history of HPV variant research suggests that, to be successful, studies will require large international collaborative efforts.

Classification and evolution of human papillomavirus genome variants: Alpha-5 (HPV26, 51, 69, 82), Alpha-6 (HPV30, 53, 56, 66), Alpha-11 (HPV34, 73), Alpha-13 (HPV54) and Alpha-3 (HPV61)

HPV variants from the same type can be classified into lineages and sublineages based on the complete genome differences and the phylogenetic topologies. We examined nucleotide variations of twelve HPV types within the species Alpha-5 (HPV26, 51, 69, 82), Alpha-6 (HPV30, 53, 56, 66), Alpha-11 (HPV34, 73), Alpha-13 (HPV54) and Alpha-3 (HPV61) by analyzing 1432 partial sequences and 181 complete genomes from multiple geographic populations. The inter-lineage and inter-sublineage mean differences of HPV variants ranged between 0.9-7.3% and 0.3-0.9%, respectively. The heterogeneity and phylogenies of HPV isolates indicate an independent evolutionary history for each type. The noncoding regions were the most variable regions whereas the capsid proteins were relatively conserved. Certain variant lineages and/or sublineages were geographically-associated. These data provide the basis to further classify HPV variants and should foster future studies on the evolution of HPV genomes and the associations of HPV variants with cancer risk.

High Whole-Genome Sequence Diversity of Human Papillomavirus Type 18 Isolates

Background: The most commonly found human papillomavirus (HPV) types in cervical cancer are HPV16 and HPV18. Genome variants of these types have been associated with differential carcinogenic potential. To date, only a handful of studies have described HPV18 whole genome sequencing results. Here we describe HPV18 variant diversity and conservation of persistent infections in a longitudinal retrospective cohort study. Methods: Cervical self-samples were obtained annually over four years and genotyped on the SPF10-DEIA-LiPA₂₅ platform. Clearing and persistent HPV18 positive infections were selected, amplified in two overlapping fragments, and sequenced using 32 sequence primers. Results: Complete viral genomes were obtained from 25 participants with persistent and 26 participants with clearing HPV18 infections, resulting in 52 unique HPV18 genomes. Sublineage A3 was predominant in this population. The consensus viral genome was completely conserved over time in persistent infections, with one exception, where different HPV18 variants were identified in follow-up samples. Conclusions: This study identified a diverse set of HPV18 variants. In persistent infections, the consensus viral genome is conserved. The identification of only one HPV18 infection with different major variants in follow-up implies that this is a potentially rare event. This dataset adds 52 HPV18 genome variants to Genbank, more than doubling the currently available HPV18 information resource, and all but one variant are unique additions.

Human papillomavirus (HPV) 18 genetic variants and cervical cancer risk in Taizhou area, China

Human papillomavirus (HPV) type 18 is predominantly associated with the development of cervical adenocarcinomas, whereas data on HPV18 genetic variability in China are limited. HPV18 genetic variants were formed phylogenetic tree, including lineages A, B, and C. We aimed to evaluate the diversity of HPV18 genetic variants by sequencing the entire E6, E7 and L1 genes. Between 2012 and 2015, a total of 138 (0.8%, 138/17669) women with single HPV18 infection were selected in this study. Finally, we observed 122 HPV18 isolates of the complete E6-E7-L1 sequences, and obtained 36 distinct variation patterns which the accession GenBank numbers as KY457805-KY457840. Except KY457805, KY457813, KY457819, KY457827, KY457829, the rest of HPV18 isolates (81.1%, 31/36) are novel variants. All of HPV18 variants belong to lineage A, while no lineage B, and C was found in our population of Taizhou region, Southeast China. Sublineage A1 was the most common variants (85.2%, 104/122), followed by sublineage A4, A3 and A5, while no sublineage A2 was obtained. Based on the tree topologies, there were three newly identified candidates' sublineages A6-A8. Out of 122 women, 67 (54.9%) had diagnosed by biopsy, including 49 women who diagnosed with cervicitis, 12 with cervical intraepithelial neoplasia (CIN)1, 4 with CIN2/3, and 2 with adenocarcinomas, respectively. Nevertheless, there was no association between HPV18 (sub) lineages and CIN1 or worse (CIN1+) lesions comparing with normal biopsies (P = .469). In conclusion, knowledge of the distribution of geographic/ethnical HPV18 genetic diversity provides critical information for developing diagnostic probes, epidemiologic correlate of cervical cancer risk and design of HPV vaccines for targeted populations.

Molecular epidemiology and phylogenetic analysis of human papillomavirus infection in women with cervical lesions and cancer from the coastal region of Ecuador

The aim of the present study was to gather information regarding the molecular epidemiology of Human papillomavirus (HPV) and related risk factors in a group of women with low- and high-grade cervical lesions and cancer from the coastal region of Ecuador. In addition, we studied the evolution of HPV variants from the most prevalent types and provided a temporal framework for their emergence, which may help to trace the source of dissemination within the region. We analyzed 166 samples, including 57 CIN1, 95 CIN2/3 and 14 cancer cases. HPV detection and typing was done by PCR-sequencing (MY09/MY11). HPV variants and estimation of the time to most recent common ancestor (tMRCA) was assessed through phylogeny and coalescence analysis. HPV DNA was found in 54.4% of CIN1, 74.7% of CIN2/3 and 78.6% of cancer samples. HPV16 (38.9%) and HPV58 (19.5%) were the most prevalent types. Risk factors for the development of cervical lesions/cancer were the following: three or more pregnancies (OR=4.3), HPV infection (OR=3.7 for high-risk types; OR=3.5 for HPV16), among others. With regard to HPV evolution, HPV16 isolates belonged to lineages A (69%) and D (31%) whereas HPV58 isolates belonged only to lineage A. The period of emergence of HPV16 was in association with human populations (tMRCA=91052 years for HPV16A and 27000 years for HPV16D), whereas HPV58A preceded Homo sapiens evolution (322257 years). This study provides novel data on HPV epidemiology and evolution in Ecuador, which will be fundamental in the vaccine era.

Differences between vaccinated and unvaccinated women explain increase in non-vaccine-type human papillomavirus in unvaccinated women after vaccine introduction

The aim of this study was to determine whether an observed increase in non-vaccine-type human papillomavirus (HPV) in unvaccinated women during the first eight years after vaccine introduction may be explained by differences in demographics or sexual behaviors, instead of type replacement. We analyzed data from three cross-sectional surveillance studies of 13-26 year-old women (total N = 1180). For women recruited from a health department clinic, older age (OR = 1.4, 95% CI: 1.2-1.6) and consistent condom use with main partner in the past 3 months (OR = 11.6, 95% CI: 3.4-40) were associated with being unvaccinated. For women recruited from a teen health center African American race (OR = 0.2, 95% CI: 0.07-0.7) and having Medicaid health insurance (OR = 0.3, 95% CI: 0.1-0.7) were inversely associated with being unvaccinated. The observed increase in non-vaccine-type HPV prevalence in unvaccinated women may be explained by differences between unvaccinated and vaccinated women.

Ancient Evolution and Dispersion of Human Papillomavirus 58 Variants

Human papillomavirus 58 (HPV58) is found in 10 to 18% of cervical cancers in East Asia but is rather uncommon elsewhere. The distribution and oncogenic potential of HPV58 variants appear to be heterogeneous, since the E7 T20I/G63S variant is more prevalent in East Asia and confers a 7- to 9-fold-higher risk of cervical precancer and cancer. However, the underlying genomic mechanisms that explain the geographic and carcinogenic diversity of HPV58 variants are still poorly understood. In this study, we used a combination of phylogenetic analyses and bioinformatics to investigate the deep evolutionary history of HPV58 complete genome variants. The initial splitting of HPV58 variants was estimated to occur 478,600 years ago (95% highest posterior density [HPD], 391,000 to 569,600 years ago). This divergence time is well within the era of speciation between Homo sapiens and Neanderthals/Denisovans and around three times longer than the modern Homo sapiens divergence times. The expansion of present-day variants in Eurasia could be the consequence of viral transmission from Neanderthals/Denisovans to non-African modern human populations through gene flow. A whole-genome sequence signature analysis identified 3 amino acid changes, 16 synonymous nucleotide changes, and a 12-bp insertion strongly associated with the E7 T20I/G63S variant that represents the A3 sublineage and carries higher carcinogenetic potential. Compared with the capsid proteins, the oncogenes E7 and E6 had increased substitution rates indicative of higher selection pressure. These data provide a comprehensive evolutionary history and genomic basis of HPV58 variants to assist further investigation of carcinogenic association and the development of diagnostic and therapeutic strategies.IMPORTANCE Papillomaviruses (PVs) are an ancient and heterogeneous group of double-stranded DNA viruses that preferentially infect the cutaneous and mucocutaneous epithelia of vertebrates. Persistent infection by specific oncogenic human papillomaviruses (HPVs), including HPV58, has been established as the primary cause of cervical cancer. In this work, we reveal the complex evolutionary history of HPV58 variants that explains the heterogeneity of oncogenic potential and geographic distribution. Our data suggest that HPV58 variants may have coevolved with archaic hominins and dispersed across the planet through host interbreeding and gene flow. Certain genes and codons of HPV58 variants representing higher carcinogenic potential and/or that are under positive selection may have important implications for viral host specificity, pathogenesis, and disease prevention.

Association of codon 72 polymorphism of p53 with the severity of cervical dysplasia, E6-T350G and HPV16 variant lineages in HPV16-infected women

Purpose. Polymorphic variability in the tumour-suppressor protein p53 at codon 72 has a considerable impact on cervical cancer development. The present study clarified the association between p53 codon 72 genotypes and the risk of cervical disease in Greek patients. We also examined whether the presence of specific p53 genotypes in combination with HPV16 variants or E6 T350G sequence variation can modify an individual's susceptibility to cervical disease.Methodology. The analysis of p53 genotypes was performed through PCR-RFLP. Sequence and phylogenetic tree analyses of the HPV16 E6 gene were also performed in order to identify HPV16 variants and T350G sequence variation.Results/Key findings. The outcomes of the present analysis revealed that women who are homozygous for the arg genotype are at a 4.17-fold higher risk of developing HPV16-associated HSIL+ (OR=4.17, 95 % CI:1.48-4.9, P=0.0049). Moreover, p53 arg/arg patients infected by an HPV16 prototype strain were associated with an increased risk of more severe lesions, while a significant relationship between the p53 arg/arg genotype in patients with T350G sequence variation and the risk of high-grade squamous intraepithelial lesions (HSILs) was revealed.Conclusion. The oncogenic potential of the virus is increased by the presence of the p53 arg/arg genotype in the Greek population in such a way that the specific protein interaction E6 (L83V)-p53 (Arg-72) can modify an individual's susceptibility to cervical disease.

Unique LCR variations among lineages of HPV16, 18 and 45 isolates from women with normal cervical cytology in Ghana

Background

In addition to being useful for classification, sequence variations of human Papillomavirus (HPV) genotypes have been implicated in differential oncogenic potential and a differential association with the different histological forms of invasive cervical cancer. These associations have also been indicated for HPV genotype lineages and sub-lineages. In order to better understand the potential implications of lineage variation in the occurrence of cervical cancers in Ghana, we studied the lineages of the three most prevalent HPV genotypes among women with normal cytology as baseline to further studies.

Methods

Of previously collected self- and health personnel-collected cervical specimen, 54, which were positive for HPV16, 18 and 45, were selected and the long control region (LCR) of each HPV genotype was separately amplified by a nested PCR. DNA sequences of 41 isolates obtained with the forward and reverse primers by Sanger sequencing were analysed.

Results

Nucleotide sequence variations of the HPV16 genotypes were observed at 30 positions within the LCR (7460 – 7840). Of these, 19 were the known variations for the lineages B and C (African lineages), while the other 11 positions had variations unique to the HPV16 isolates of this study. For the HPV18 isolates, the variations were at 35 positions, 22 of which were known variations of Africa lineages and the other 13 were unique variations observed for the isolates obtained in this study (at positions 7799 and 7813). HPV45 isolates had variations at 35 positions and 2 (positions 7114 and 97) were unique to the isolates of this study.

Conclusion

This study provides the first data on the lineages of HPV 16, 18 and 45 isolates from Ghana. Although the study did not obtain full genome sequence data for a comprehensive comparison with known lineages, these genotypes were predominately of the Africa lineages and had some unique sequence variations at positions that suggest potential oncogenic implications. These data will be useful for comparison with lineages of these genotypes from women with cervical lesion and all the forms of invasive cervical cancers.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-017-0755-z) contains supplementary material, which is available to authorized users.

HPV Genotyping of Modified General Primer-Amplicons Is More Analytically Sensitive and Specific by Sequencing than by Hybridization

Sensitive and specific genotyping of human papillomaviruses (HPVs) is important for population-based surveillance of carcinogenic HPV types and for monitoring vaccine effectiveness. Here we compare HPV genotyping by Next Generation Sequencing (NGS) to an established DNA hybridization method. In DNA isolated from urine, the overall analytical sensitivity of NGS was found to be 22% higher than that of hybridization. NGS was also found to be the most specific method and expanded the detection repertoire beyond the 37 types of the DNA hybridization assay. Furthermore, NGS provided an increased resolution by identifying genetic variants of individual HPV types. The same Modified General Primers (MGP)-amplicon was used in both methods. The NGS method is described in detail to facilitate implementation in the clinical microbiology laboratory and includes suggestions for new standards for detection and calling of types and variants with improved resolution.

The Biological Impact of Genomic Diversity in Cervical Cancer Development

Human papillomaviruses (HPVs) are the etiologic agents of cervical cancer, the unique human neoplasia that has one single necessary cause. The diversity of HPVs is well described, with 200 HPV types existing as distinct taxonomic units and each receiving an Arabic number. On a clinical basis, they are usually grouped by their site of occurrence and disease associations. Those types inhabiting the anogenital mucosa are more intensively studied and further divided into cancer-associated HPVs, which are termed 'high risk', while those linked to benign proliferative lesions are assigned as 'low risk'. HPV16 is responsible for approximately 50% of all ICC cases, and paradoxically is one of the most prevalent types among healthy women. Longitudinal studies have shown that when an incidental HPV16 infection becomes persistent it will result in an enhanced risk for the development of high-grade lesions. However, it is unknown why some persistent, HPV16 infections (or infections by other HR-HPV types) progress to CIN3+ while most clear spontaneously. Several epidemiological investigations have focused on cofactors, from the most obvious such as cigarette and other carcinogenic exposures, to coinfections by other STDs such as chlamydia, with no significant findings. Thus, the current focus is on genomic variation from both virus and host. Such studies have been potentialized by the enormous technical advances in nucleic acid sequencing, allowing this relationship to be broadly interrogated. Corroborating subgenomic data from decades ago, an association between HPV16 lineages and carcinogenesis is being revealed. However, this effect does not seem to apply across female populations from different continents/ethnicities, again highlighting a role played by HPV16 adaptation and evasion from the host over time.

Global Genomic Diversity of Human Papillomavirus 11 Based on 433 Isolates and 78 Complete Genome Sequences

Human papillomavirus 11 (HPV11) is an etiological agent of anogenital warts and laryngeal papillomas and is included in the 4-valent and 9-valent prophylactic HPV vaccines. We established the largest collection of globally circulating HPV11 isolates to date and examined the genomic diversity of 433 isolates and 78 complete genomes (CGs) from six continents. The genomic variation within the 2,800-bp E5a-E5b-L1-upstream regulatory region was initially studied in 181/207 (87.4%) HPV11 isolates collected for this study. Of these, the CGs of 30 HPV11 variants containing unique single nucleotide polymorphisms (SNPs), indels (insertions or deletions), or amino acid changes were fully sequenced. A maximum likelihood tree based on the global alignment of 78 HPV11 CGs (30 CGs from our study and 48 CGs from GenBank) revealed two HPV11 lineages (lineages A and B) and four sublineages (sublineages A1, A2, A3, and A4). HPV11 (sub)lineage-specific SNPs within the CG were identified, as well as the 208-bp representative region for CG-based phylogenetic clustering within the partial E2 open reading frame and noncoding region 2. Globally, sublineage A2 was the most prevalent, followed by sublineages A1, A3, and A4 and lineage B.

IMPORTANCE

This collaborative international study defined the global heterogeneity of HPV11 and established the largest collection of globally circulating HPV11 genomic variants to date. Thirty novel complete HPV11 genomes were determined and submitted to the available sequence repositories. Global phylogenetic analysis revealed two HPV11 variant lineages and four sublineages. The HPV11 (sub)lineage-specific SNPs and the representative region identified within the partial genomic region E2/noncoding region 2 (NCR2) will enable the simpler identification and comparison of HPV11 variants worldwide. This study provides an important knowledge base for HPV11 for future studies in HPV epidemiology, evolution, pathogenicity, prevention, and molecular assay development.

Naturally Occurring Major and Minor Capsid Protein Variants of Human Papillomavirus 45 (HPV45): Differential Recognition by Cross-Neutralizing Antibodies Generated by HPV Vaccines

We investigated naturally occurring variation within the major (L1) and minor (L2) capsid proteins of human papillomavirus genotype 45 (HPV45). Pseudoviruses (PsVs) representing HPV45 sublineages A1, A2, A3, B1, and B2 exhibited comparable particle-to-infectivity ratios and morphologies but demonstrated both increased (A2, A3, and B1) and decreased (B2) sensitivities to cross-neutralization by HPV vaccine antibodies compared to that of the A1 sublineage. Mutant PsVs identified HI loop residue 357 as being critical for conferring this differential sensitivity.

Could the human papillomavirus vaccines drive virulence evolution?

The human papillomavirus (HPV) vaccines hold great promise for preventing several cancers caused by HPV infections. Yet little attention has been given to whether HPV could respond evolutionarily to the new selection pressures imposed on it by the novel immunity response created by the vaccine. Here, we present and theoretically validate a mechanism by which the vaccine alters the transmission-recovery trade-off that constrains HPV's virulence such that higher oncogene expression is favoured. With a high oncogene expression strategy, the virus is able to increase its viral load and infected cell population before clearance by the vaccine, thus improving its chances of transmission. This new rapid cell-proliferation strategy is able to circulate between hosts with medium to high turnover rates of sexual partners. We also discuss the importance of better quantifying the duration of challenge infections and the degree to which a vaccinated host can shed virus. The generality of the models presented here suggests a wider applicability of this mechanism, and thus highlights the need to investigate viral oncogenicity from an evolutionary perspective.

High-level of viral genomic diversity in cervical cancers: A Brazilian study on human papillomavirus type 16

Invasive cervical cancer (ICC) is the third most frequent cancer among women worldwide and is associated with persistent infection by carcinogenic human papillomaviruses (HPVs). The combination of large populations of viral progeny and decades of sustained infection may allow for the generation of intra-patient diversity, in spite of the assumedly low mutation rates of PVs. While the natural history of chronic HPVs infections has been comprehensively described, within-host viral diversity remains largely unexplored. In this study we have applied next generation sequencing to the analysis of intra-host genetic diversity in ten ICC and one condyloma cases associated to single HPV16 infection. We retrieved from all cases near full-length genomic sequences. All samples analyzed contained polymorphic sites, ranging from 3 to 125 polymorphic positions per genome, and the median probability of a viral genome picked at random to be identical to the consensus sequence in the lesion was only 40%. We have also identified two independent putative duplication events in two samples, spanning the L2 and the L1 gene, respectively. Finally, we have identified with good support a chimera of human and viral DNA. We propose that viral diversity generated during HPVs chronic infection may be fueled by innate and adaptive immune pressures. Further research will be needed to understand the dynamics of viral DNA variability, differentially in benign and malignant lesions, as well as in tissues with differential intensity of immune surveillance. Finally, the impact of intralesion viral diversity on the long-term oncogenic potential may deserve closer attention.

Evolution and classification of oncogenic human papillomavirus types and variants associated with cervical cancer

The nomenclature of human papillomavirus (HPV) is established by the International Committee on Taxonomy of Virus (ICTV). However, the ICTV does not set standards for HPV below species levels. This chapter describes detailed genotyping methods for determining and classifying HPV variants.

Genetic characterization and clinical implications of human papillomavirus type 16 (HPV16) variants from northeastern Argentina

BACKGROUND

Human papillomavirus type 16 (HPV16) plays a central role in the development of cervical cancer. Worldwide studies indicate the existence of HPV16 variants that show different geographic distributions and oncogenic potential.

OBJECTIVE

Our goal was to describe the genetic variation of HPV16 isolates identified in urban women with different grades of cervical lesions living in northeastern Argentina.

STUDY DESIGN

We analyzed 116 HPV16-positive cervical samples (16 NLIM, 62 L-SIL, 16 H-SIL and 22 cervical cancer) from patients attending health centers in Misiones (Argentina) during 2006-13. HPV16 isolates were genetically characterized through PCR amplification and direct sequencing of 364 bp within the long control region, and the resulting sequences classified into variants based on phylogenetic analysis (lineages A, B, C and D). A potential association between HPV16 variants and lesion grade was evaluated through an odds ratio (OR) test. A temporal framework for the origin of HPV16 variants was assessed through coalescence analysis (BEAST v 1.7.5).

RESULTS

Phylogenetic analysis of HPV16 sequences showed that 92.1% of the samples clustered with lineage A, and 6.9% to lineage D. HPV16 variants from lineage D were more frequently associated with high-grade lesions and cancer (HSIL+) than lineage A variants at an OR of 13.8 (1.6-117.0). The time to most common recent ancestor (tMCRA) of all variants was 119,103 years before present (HPD 95%=48,486-197,239), a date consistent with the time frame for modern human evolution.

CONCLUSION

Our results suggest that HPV16 variants from lineage D may represent an additional risk factor for the development of cervical cancer in women living in northeastern Argentina. This study provides new information about viral isolates present in Argentina that will contribute to the monitoring of HPV16 infection in the vaccine era.

Prevalence of HPV 16 and HPV 18 lineages in Galicia, Spain

Genetic variants of human papillomavirus types 16 and 18 (HPV16/18) could differ in their cancer risk. We studied the prevalence and association with high-grade cervical lesions of different HPV16/18 variant lineages in a case-control study including 217 cases (cervical intraepithelial neoplasia grade 2 or grade 3 or worse: CIN2 or CIN3+) and 116 controls (no CIN2 or CIN3+ in two-year follow-up). HPV lineages were determined by sequencing the long control region (LCR) and the E6 gene. Phylogenetic analysis of HPV16 confirmed that isolates clustered into previously described lineages: A (260, 87.5%), B (4, 1.3%), C (8, 2.7%), and D (25, 8.4%). Lineage D/lineage A strains were, respectively, detected in 4/82 control patients, 19/126 CIN3+ cases (OR = 3.1, 95%CI: 1.0–12.9, p = 0.04), 6/1 glandular high-grade lesions (OR = 123, 95%CI: 9.7–5713.6, p<0.0001), and 4/5 invasive lesions (OR = 16.4, 95%CI: 2.2–113.7, p = 0.002). HPV18 clustered in lineages A (32, 88.9%) and B (4, 11.1%). Lineage B/lineage A strains were respectively detected in 1/23 control patients and 2/5 CIN3+ cases (OR = 9.2, 95%CI: 0.4–565.4, p = 0.12). In conclusion, lineages A of HPV16/18 were predominant in Spain. Lineage D of HPV16 was associated with increased risk for CIN3+, glandular high-grade lesions, and invasive lesions compared with lineage A. Lineage B of HPV18 may be associated with increased risk for CIN3+ compared with lineage A, but the association was not significant. Large well-designed studies are needed before the application of HPV lineage detection in clinical settings.

Human papillomavirus type 16 long control region and E6 variants stratified by cervical disease stage

•

We sequenced HPV16 LCR–E6 variants in cervical disease samples from England.

•

98% of variants were of the EUR lineage.

•

Site-specific entropy identified several variable sites in the LCR and E6.

•

No single or combination of sites were associated with disease, including E6 T350G.

Objective

Certain intra-type variants of HPV16 have been shown to be associated with an increased risk of developing high grade cervical disease, but their potential association is confounded by apparent geographic and phylogenetic lineage dependency. The objective of this study was to evaluate the relationship between HPV16 sequence variants and cervical disease stage in monospecific infection samples from a single lineage (European, EUR) in England.

Methods

One hundred and twelve women singly infected with HPV16 and displaying normal and abnormal cytology grades were selected. An 1187 bp fragment encompassing the entire LCR and a portion of the E6 open reading frame was sequenced to identify intra-type variants. Intra-type diversity was estimated using Shannon entropy.

Results

Almost all samples (110/112; 98%) were assigned to the EUR lineage, one sample was classified as European-Asian (EAS) and another African (Afr1a). The mean pairwise distance of the EUR sequences in this study was low (0.29%; 95%CI 0.13–0.45%) but there were nevertheless several sites in the LCR (n = 5) and E6 (n = 2) that exhibited a high degree of entropy. None of these sites, however, including the T350G non-synonymous (L83V) substitution in E6, alone or in combination, were found to be associated with cervical disease stage.

Conclusions

Despite using single infection samples and samples from a single variant lineage, intra-type variants of HPV16 were not differentially associated with cervical disease. Monitoring intra-lineage, site-specific variants, such as T350G, is unlikely to be of diagnostic value.

Human papillomavirus 45 genetic variation and cervical cancer risk worldwide

Human papillomavirus 45 (HPV45) is a member of the HPV18-related alpha-7 species and accounts for approximately 5% of all cervical cancer cases worldwide. This study evaluated the genetic diversity of HPV45 and the association of HPV45 variants with the risk of cervical cancer by sequencing the entire E6 and E7 open reading frames of 300 HPV45-positive cervical samples from 36 countries. A total of 43 HPV45 sequence variants were identified that formed 5 phylogenetic sublineages, A1, A2, A3, B1, and B2, the distribution of which varied by geographical region. Among 192 cases of cervical cancer and 101 controls, the B2 sublineage was significantly overrepresented in cervical cancer, both overall and in Africa and Europe separately. We show that the sequence analysis of E6 and E7 allows the classification of HPV45 variants and that the risk of cervical cancer may differ by HPV45 variant sublineage.

IMPORTANCE

This work describes the largest study to date of human papillomavirus 45 (HPV45)-positive cervical samples and provides a comprehensive reference for phylogenetic classification for use in epidemiological studies of the carcinogenicity of HPV45 genetic variants, particularly as our findings suggest that the B2 sublineage of HPV45 is associated with a higher risk of cervical cancer.

Human papillomavirus genomics: past, present and future

Human papillomaviruses (HPV) are a group of divergent DNA viruses, of which a select few evolutionarily related HPVs have emerged to be highly oncogenic and of significant medical importance. Essentially all cases of cervical cancer, as well as a subset of other anogenital and oral cancers are caused by this limited set of HPV types. At present, over 150 HPV types have been identified and may be classified into genera, species and types based upon comparison of the viral genome. Established nucleotide phylogenies sort the highly pathogenic HPV types to the genus Alphapapillomavirus (α-PV). A species group includes viral types with 60-70% genomic nucleotide similarity that share a most-recent common ancestor; for example the species group's alpha-9 (HPV16-related) and alpha-7 (HPV18-related), contain the majority of known oncogenic HPV types. Genomes from the same HPV type with 1-10% nucleotide differences designate HPV variant lineages. The established nucleotide variations observed in extant HPV genomes have been fixed through evolutionary processes prior to human population expansion and global dissemination. To characterize viral types and variants associated with pathology for clinical applications (e.g. screening), molecular epidemiological studies have proven essential for identifying links between HPV natural history and carcinogenicity. This chapter presents a historical account of HPV genomics in the context of major discoveries and advances over the past 2 thousand years.

Genomic diversity of human papillomaviruses (HPV) and clinical implications: an overview in adulthood and childhood

During the last years, several researchers have highlighted the importance of characterizing more than one genomic region in order to detect recombination and classify variants of human papillomaviruses (HPVs) properly. HPVs variants differ in their biological, molecular and chemical properties. Therefore, this genomic diversity can present differences in the natural history and pathogenicity of HPVs. Different 'high-risk' HPVs variants of the genotypes HPV 16 and 18 can confer varied risks of viral persistence in the human cervix and influence HPVs progression to cervical cancer. Moreover, different 'low-risk' HPVs variants of the genotypes HPV 6 and 11 can play a unique role in the development of anogenital and cutaneous warts, recurrent respiratory papillomatosis (RRP) and ophthalmic pterygium. In future, the precise impact of genomic HPVs diversity to the clinical course of HPVs-associated diseases as well as to the efficacy of the current HPVs vaccines remains to be elucidated.

Molecular and evolutionary analysis of HPV16 E6 and E7 genes in Greek women

Human papillomavirus type 16 (HPV16) non-European variants have been associated with persistent infection and cervical cancer development, while the L83V variant of the E6 gene has been correlated with the progression of cervical malignancy. The present study investigated the presence of the HPV16 L83V variant in Greek women. Molecular evolutionary analysis of the HPV16 E6 and E7 oncogenes was conducted in order to estimate the evolution of the HPV16 genome in the Greek population. The E6 L83V variant was found in 78.2 % of high- and 64.28 % of low-grade specimens. Moreover, the prototype and E6 L83V variants were both prevalent in high- and low-grade malignancies in Greek women. Selective pressure analysis of the individual amino acid residues of HPV16 sequences from the Greek population indicates that codon 83 of the E6 protein, as well as codon 85 of the E7 protein, are undergoing positive selection. Novel sequence variations were recorded within the E6 and E7 genes in cervical samples, characterized as (T350G) European variants. However, no signal of intratypic recombination event was identified within the E6–E7 region. Molecular and evolutionary analyses of HPV16 genomes from distinct geographical locations might provide valuable information about viral evolution and oncogenecity.

Human papillomavirus genome variants

Amongst the human papillomaviruses (HPVs), the genus Alphapapillomavirus contains HPV types that are uniquely pathogenic. They can be classified into species and types based on genetic distances between viral genomes. Current circulating infectious HPVs constitute a set of viral genomes that have evolved with the rapid expansion of the human population. Viral variants were initially identified through restriction enzyme polymorphisms and more recently through sequence determination of viral fragments. Using partial sequence information, the history of variants, and the association of HPV variants with disease will be discussed with the main focus on the recent utilization of full genome sequence information for variant analyses. The use of multiple sequence alignments of complete viral genomes and phylogenetic analyses have begun to define variant lineages and sublineages using empirically defined differences of 1.0-10.0% and 0.5-1.0%, respectively. These studies provide the basis to define the genetics of HPV pathogenesis.

The Merkel cell polyomavirus minor capsid protein

The surface of polyomavirus virions is composed of pentameric knobs of the major capsid protein, VP1. In previously studied polyomavirus species, such as SV40, two interior capsid proteins, VP2 and VP3, emerge from the virion to play important roles during the infectious entry process. Translation of the VP3 protein initiates at a highly conserved Met-Ala-Leu motif within the VP2 open reading frame. Phylogenetic analyses indicate that Merkel cell polyomavirus (MCV or MCPyV) is a member of a divergent clade of polyomaviruses that lack the conserved VP3 N-terminal motif. Consistent with this observation, we show that VP3 is not detectable in MCV-infected cells, VP3 is not found in native MCV virions, and mutation of possible alternative VP3-initiating methionine codons did not significantly affect MCV infectivity in culture. In contrast, VP2 knockout resulted in a >100-fold decrease in native MCV infectivity, despite normal virion assembly, viral DNA packaging, and cell attachment. Although pseudovirus-based experiments confirmed that VP2 plays an essential role for infection of some cell lines, other cell lines were readily transduced by pseudovirions lacking VP2. In cell lines where VP2 was needed for efficient infectious entry, the presence of a conserved myristoyl modification on the N-terminus of VP2 was important for its function. The results show that a single minor capsid protein, VP2, facilitates a post-attachment stage of MCV infectious entry into some, but not all, cell types.

Merkel cell polyomavirus (MCV or MCPyV) is a recently discovered member of the viral family Polyomaviridae. The virus plays a causal role in Merkel cell carcinoma, a highly lethal form of skin cancer. MCV encodes a major capsid protein, VP1, which forms the non-enveloped surface of the virion. Other polyomavirus species encode two minor capsid proteins, VP2 and VP3, which associate with the inner surface of the capsid and facilitate infectious entry. In this report we show that MCV does not have a VP3 minor capsid protein. Sequence analyses suggest that more than a quarter of known polyomavirus species share MCV's lack of a VP3 protein. In contrast to VP3, VP2-knockout MCV mutants displayed dramatically reduced infectivity. Consistent with native virion findings, MCV pseudovirions lacking VP2 or carrying mutations in the VP2 myristoylation motif displayed reduced infectivity on several cell lines. Puzzlingly, MCV pseudoviruses lacking VP2 successfully transduced other cell lines with high efficiency. Taken together, the data show that the lone MCV minor capsid protein, VP2, plays an important role during infectious entry into some cell types, but is dispensable for entry into other cell types.

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.

Amino acid sequence diversity of the major human papillomavirus capsid protein: implications for current and next generation vaccines

•

We evaluated amino acid diversity of the major capsid protein of HPV.

•

Residues displaying high entropy were found within surface-exposed domains.

•

We discuss the implications of this diversity on the current and next generation HPV vaccines.

Despite the fidelity of host cell polymerases, the human papillomavirus (HPV) displays a degree of genomic polymorphism resulting in distinct genotypes and intra-type variants. The current HPV vaccines target the most prevalent genotypes associated with cervical cancer (HPV16/18) and genital warts (HPV6/11). Although these vaccines confer some measure of cross-protection, a multivalent HPV vaccine is in the pipeline that aims to broaden vaccine protection against other cervical cancer-associated genotypes including HPV31, HPV33, HPV45, HPV52 and HPV58. Both current and next generation vaccines comprise virus-like particles, based upon the major capsid protein, L1, and vaccine-induced, type-specific protection is likely mediated by neutralizing antibodies targeting L1 surface-exposed domains.

The aim of this study was to perform an in silico analysis of existing full length L1 sequences representing vaccine-relevant HPV genotypes in order to address the degree of naturally-occurring, intra-type polymorphisms. In total, 1281 sequences from the Americas, Africa, Asia and Europe were assembled. Intra-type entropy was low and/or limited to non-surface-exposed residues for HPV6, HPV11 and HPV52 suggesting a minimal effect on vaccine antibodies for these genotypes. For HPV16, intra-type entropy was high but the present analysis did not reveal any significant polymorphisms not previously identified. For HPV31, HPV33, HPV58, however, intra-type entropy was high, mostly mapped to surface-exposed domains and in some cases within known neutralizing antibody epitopes. For HPV18 and HPV45 there were too few sequences for a definitive analysis, but HPV45 displayed some degree of surface-exposed residue diversity. In most cases, the reference sequence for each genotype represented a minority variant and the consensus L1 sequences for HPV18, HPV31, HPV45 and HPV58 did not reflect the L1 sequence of the currently available HPV pseudoviruses. These data highlight a number of variant amino acid residues that warrant further investigation for vaccine and natural history studies of HPV.

Human papillomavirus infection among sexual partners attending a Sexually Transmitted Disease Clinic in Rio de Janeiro, Brazil

Cervical cancer is a major source of illness and death among women worldwide and genital infection with oncogenic human papillomavirus (HPV) its principal cause. There is evidence of the influence of the male factor in the development of cervical neoplasia. Nevertheless, the pathogenic processes of HPV in men are still poorly understood. It has been observed that different HPV types can be found among couples. The objective of the present study was to investigate HPV infections in female patients (n = 60 females/group) as well as in their sexual partners and to identify the concordance of HPV genotypes among them. By using the polymerase chain reaction, we detected a 95% prevalence of HPV DNA in women with cervical intraepithelial neoplasia (CIN) compared to 18.3% in women with normal cervical epithelium, with a statistically significant difference (P < 0.001). The HPV DNA prevalence was 50% in male partners of women with CIN and 16.6% in partners of healthy women. In the control group (healthy women), only 9 couples were simultaneously infected with HPV, and only 22.2% of them had the same virus type, showing a weak agreement rate (kappa index = 0.2). Finally, we observed that HPV DNA was present in both partners in 30 couples if the women had CIN, and among them, 53.3% shared the same HPV type, showing moderate agreement, with a kappa index of 0.5. This finding supports the idea of circulation and recirculation of HPV among couples, perpetuating HPV in the sexually active population, rather than true recurrences of latent infections.

Context-dependent substitution models for circular DNA

The most general context-dependent Markov substitution process, where each substitution event involves only one site and substitution rates depend on the whole sequence, is presented for the first time. The focus is on circular DNA sequences, where the problem of specifying the behaviour of the first and last sites in a linear sequence does not arise. Important special cases include (1) the established models where each site behaves independently, (2) models which are increasingly applied to non-coding DNA, where each site depends on only the immediate neighbouring sites, and (3) models where each site depends on two closest neighbours on both sides, such as the codon models. These special cases are classified and illustrated by published models. It is shown that the existing codon substitution models mix up the mutation and selection processes, rendering the substitution rates challenging to interpret. The classification suggests the study of a more interpretable codon model, where the mutation and selection processes are clearly delineated. Furthermore, this model allows a natural accommodation of possibly different selection pressures in overlapping reading frames, which may contribute to furthering the understanding of viral diseases. Also included are brief discussions on the stationary distribution of a context-dependent substitution process and a simple recipe for simulating it on a computer.