Papillomavirus Episteme Research Articles

Generating human papillomavirus (HPV) reference databases to maximize genomic mapping.

Genomic experiments analyzing human papillomaviruses (HPVs) require a carefully selected list of sequences as a reference database to map millions of reads. The available sources, such as the Papillomavirus Episteme (PaVE), are organized based on variations in the L1 gene rather than the whole HPV sequence. Moreover, the PaVE process uses complex multiple sequence alignments containing hundreds or thousands of sequences. These issues complicate the generation of a reference database for genomics, leading to the generation of per-analysis-defined databases. Here, we propose a de novo strategy considering all HPV sequences reported in the NCBI database to define a subset of highly representative HPV sequences. The strategy is based on oligonucleotide frequency profiling of the whole sequence followed by hierarchical clustering. Using data from HPV capture experiments, we demonstrate that this strategy selects suitable sequences as a reference database to map most mappable reads unambiguously. We provide some recommendations to improve HPV mapping. The generated .fasta files can be accessed at https://github.com/vtrevino/HPV-Ref-Genomes .

HPV DeepSeq: An Ultra-Fast Method of NGS Data Analysis and Visualization Using Automated Workflows and a Customized Papillomavirus Database in CLC Genomics Workbench.

Next-generation sequencing (NGS) has actualized the human papillomavirus (HPV) virome profiling for in-depth investigation of viral evolution and pathogenesis. However, viral computational analysis remains a bottleneck due to semantic discrepancies between computational tools and curated reference genomes. To address this, we developed and tested automated workflows for HPV taxonomic profiling and visualization using a customized papillomavirus database in the CLC Microbial Genomics Module. HPV genomes from Papilloma Virus Episteme were customized and incorporated into CLC “ready-to-use” workflows for stepwise data processing to include: (1) Taxonomic Analysis, (2) Estimate Alpha/Beta Diversities, and (3) Map Reads to Reference. Low-grade (n = 95) and high-grade (n = 60) Pap smears were tested with ensuing collective runtimes: Taxonomic Analysis (36 min); Alpha/Beta Diversities (5 s); Map Reads (45 min). Tabular output conversion to visualizations entailed 1–2 keystrokes. Biodiversity analysis between low- (LSIL) and high-grade squamous intraepithelial lesions (HSIL) revealed loss of species richness and gain of dominance by HPV-16 in HSIL. Integrating clinically relevant, taxonomized HPV reference genomes within automated workflows proved to be an ultra-fast method of virome profiling. The entire process named “HPV DeepSeq” provides a simple, accurate and practical means of NGS data analysis for a broad range of applications in viral research.

Abstract 2264: Molecular detection and genotyping of human papillomavirus in presence of some oncogenic viruses among Egyptian breast cancer women

Abstract Introduction High risk human papillomaviruses (HR-HPVs) has a possible role in development and pathogenesis of breast cancer (BC) disease. Other oncogenic viruses showed to be linked to BC such as Epstein Barr virus (EBV), and Human mammary tumor virus (HMTV). Presence of high risk HPV genotypes, EBV or HMTV could affect the severity course of BC. Therefore, the aim of the current study is to explore distribution of HPV with regards to its genotypes and the presence of EBV and HMTV among Egyptian BC women with its correlation to age and hormonal status. Patients & Methods This study was conducted on 40 fresh tissues of BC patients and blood samples from 30 apparently healthy controls. All samples were investigated for detection of HPV, EBV and HMTV DNAs using qualitative PCR. Detection of HPV viral load was performed using quantitative real time PCR. HPV genotypes were determined on 10 purified positive PCR product using Sanger sequencing of HPV L1 region. The obtained sequences were compared to the known HPV DNA sequences at National center for biotechnology information (NCBI) and/or Papillomavirus Episteme (Pave) data bases using their online BLAST and taxonomy tools. Data were correlated with age and hormonal status. Results All the tested 3 oncogenic viruses were not detected into 30 controls. HPV, EBV and HMTV were detected in 20 (50%), 23 (57.5%), 28 (70%) respectively of the 40 BC patients by using PCR assay. HPV median load of the 20 positive samples was 4800 copies/μL. Seven different HPV genotypes were detected included high risk genotypes ((HPV 18 in 3 samples), 70 and 52), low risk genotype (11) detected into only one sample along with 3 other genotypes (102, 114 and 151). One sample showed no significant similarity with any sequence registered in the used databases. Presence of EBV and HMTV DNAs and high risk HPV genotypes was more associated with young age (lower than 50 years), P value = 0.015, 0.001 and 0.05, respectively but were not associated with hormonal status of BC patients. Conclusion HMTV was the most detected viral DNAs among BC patients, High risk HPVs were identified into 10% of studied Egyptian BC patients and were associated with young age. Further larger sample size study is needed to evaluate role of HPV genotypes on pathogenesis and outcome of BC disease. STDF Acknowledgment This project was supported financially by the Science and Technology Development Fund (STDF), Egypt, Grant No.22944. Citation Format: Nasra F Abdel Fattah, Shimaa A Metwally, Samah A Loutfy, Maha A Abo-Shadi, Ahmed B Barakat, Omar A Rabee, Ahmed M Osman, Amany M Helal, Tarek M Hashem, Manar M Moneer. Molecular detection and genotyping of human papillomavirus in presence of some oncogenic viruses among Egyptian breast cancer women [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2264.

PuMA: A papillomavirus genome annotation tool.

High-throughput sequencing technologies provide unprecedented power to identify novel viruses from a wide variety of (environmental) samples. The field of ‘viral metagenomics’ has dramatically expanded our understanding of viral diversity. Viral metagenomic approaches imply that many novel viruses will not be described by researchers who are experts on (the genomic organization of) that virus family. We have developed the papillomavirus annotation tool (PuMA) to provide researchers with a convenient and reproducible method to annotate and report novel papillomaviruses. PuMA currently correctly annotates 99% of the papillomavirus genes when benchmarked against the 655 reference genomes in the papillomavirus episteme. Compared to another viral annotation pipeline, PuMA annotates more viral features while being more accurate. To demonstrate its general applicability, we also developed a preliminary version of PuMA that can annotate polyomaviruses. PuMA is available on GitHub (https://github.com/KVD-lab/puma) and through the iMicrobe online environment (https://www.imicrobe.us/#/apps/puma).

Detection of Alpha, Beta, Gamma, and Unclassified Human Papillomaviruses in Cervical Cancer Samples From Mexican Women.

Background: Cervical cancer (CC) is associated to high-risk human papillomavirus (HPV) infections, for this reason it is crucial to have sensitive and accurate HPV diagnostic tests. To date, most research is focused on HPVs within the Alphapapillomavirus (α-PVs) genus and little attention has been paid to cervical infections with other HPV genotypes, like those of the Betapapillomavirus (β-PVs) and Gammapapillomavirus (γ-PVs) genera. The aim of this study was to determine the HPV genotypes from different genera in women with CC using Next-Generation Sequencing (NGS).Methods: The study comprised 48 HPV positive CC samples evaluated with the Linear Array HPV Genotyping test and individually sequenced by 454 NGS using PGMY09/11 and FAP primers. To determine the HPV genotypes present in each sample, the obtained sequences were compared with all HPV L1 gene reference sequences from the Papillomavirus Episteme database (PaVE). Moreover, 50 HPV positive low-grade cervical lesion samples individually genotyped with NGS were also included to determine the genotypes present preferentially in CC patients.Results: Among the 48 CC samples, 68.75% consisted of multiple HPV infections, 51 different genotypes were detected, of which 7 are still unclassified, 28 belong to α-PVs (6, 11, 16, 18, 26, 30, 33, 35, 39, 42, 43, 44, 45, 51, 52, 53, 54, 59, 62, 66, 68, 69, 70, 71, 74, 81, 102, 114), 10 to β-PVs (5, 12, 21, 37, 38b, 47, 80, 107, 118, 122), and 6 to γ-PVs (101, 103, 123, 135, 147, 214). Among them, HPV16 was the most prevalent genotype (54.2%), followed by HPV18 (16.7%), HPV38b (14.6%), and HPVs 52/62/80 (8.3%). Some genotypes were exclusively found in CC when compared with Cervical Intraepithelial Neoplasia grade 1 (CIN1) samples, such as HPVs 5, 18, 38b, 107, 122, FA39, FA116, mSK_120, and mSK_136.Conclusions: This work demonstrates the great diversity of HPV genotypes detected by combining PGMY and FAP primers with NGS in cervical swabs. The relatively high attribution of β- and γ- PVs in CC samples suggest their possible role as carcinogenic cofactors, but deeper studies need to be performed to determine if they have transforming properties and the significance of HPV-coinfections.

An economical Nanopore sequencing assay for human papillomavirus (HPV) genotyping

BackgroundHuman papillomavirus (HPV) testing has been employed by several European countries to augment cytology-based cervical screening programs. A number of research groups have demonstrated potential utility of next-generation sequencing (NGS) for HPV genotyping, with comparable performance and broader detection spectrum than current gold standards. Nevertheless, most of these NGS platforms may not be the best choice for medium sample throughput and laboratories with less resources and space. In light of this, we developed a Nanopore sequencing assay for HPV genotyping and compared its performance with cobas HPV Test and Roche Linear Array HPV Genotyping Test (LA).MethodsTwo hundred and one cervicovaginal swabs were routinely tested for Papanicolaou smear, cobas HPV Test and LA. Residual DNA was used for Nanopore protocol after routine testing. Briefly, HPV L1 region was amplified using PGMY and MGP primers, and PCR-positive specimens were sequenced on MinION flow cells (R9.4.1). Data generated in first 2 h were aligned with reference sequences from Papillomavirus Episteme database for genotyping.ResultsNanopore detected 96 HPV-positive (47.76%) and 95 HPV-negative (47.26%) specimens, with 10 lacking β-globin band and not further analyzed (4.98%). Substantial agreement was achieved with cobas HPV Test and LA (κ: 0.83–0.93). In particular, Nanopore appeared to be more sensitive than cobas HPV Test for HPV 52 (n = 7). For LA, Nanopore revealed higher concordance for high-risk (κ: 0.93) than non-high risk types (κ: 0.83), and with similar high-risk positivity in each cytology grading. Nanopore also provided better resolution for HPV 52 in 3 specimens co-infected with HPV 33 or 58, and for HPV 87 which was identified as HPV 84 by LA. Interestingly, Nanopore identified 5 additional HPV types, with an unexpected high incidence of HPV 90 (n = 12) which was reported in North America and Belgium but not in Hong Kong.ConclusionsWe developed a Nanopore workflow for HPV genotyping which was economical (about USD 50.77 per patient specimen for 24-plex runs), and with comparable or better performance than 2 reference methods in the market. Future prospective study with larger sample size is warranted to further evaluate test performance and streamline the protocol.

The Papillomavirus Episteme: a major update to the papillomavirus sequence database.

The Papillomavirus Episteme (PaVE) is a database of curated papillomavirus genomic sequences, accompanied by web-based sequence analysis tools. This update describes the addition of major new features. The papillomavirus genomes within PaVE have been further annotated, and now includes the major spliced mRNA transcripts. Viral genes and transcripts can be visualized on both linear and circular genome browsers. Evolutionary relationships among PaVE reference protein sequences can be analysed using multiple sequence alignments and phylogenetic trees. To assist in viral discovery, PaVE offers a typing tool; a simplified algorithm to determine whether a newly sequenced virus is novel. PaVE also now contains an image library containing gross clinical and histopathological images of papillomavirus infected lesions. Database URL: https://pave.niaid.nih.gov/.

Role of the host restriction factor APOBEC3 on papillomavirus evolution.

More than 270 different types of papillomaviruses have been discovered in a wide array of animal species. Despite the great diversity of papillomaviruses, little is known about the evolutionary processes that drive host tropism and the emergence of oncogenic genotypes. Although host defense mechanisms have evolved to interfere with various aspects of a virus life cycle, viruses have also coevolved copious strategies to avoid host antiviral restriction. Our and other studies have shown that the cytidine deaminase APOBEC3 family members edit HPV genomes and restrict virus infectivity. Thus, we hypothesized that host restriction by APOBEC3 served as selective pressure during papillomavirus evolution. To test this hypothesis, we analyzed the relative abundance of all dinucleotide sequences in full-length genomes of 274 papillomavirus types documented in the Papillomavirus Episteme database (PaVE). Here, we report that TC dinucleotides, the preferred target sequence of several human APOBEC3 proteins (hA3A, hA3B, hA3F, and hA3H), are highly depleted in papillomavirus genomes. Given that HPV infection is highly tissue-specific, the expression levels of APOBEC3 family members were analyzed. The basal expression levels of all APOBEC3 isoforms, excluding hA3B, are significantly higher in mucosal skin compared with cutaneous skin. Interestingly, we reveal that Alphapapillomaviruses (alpha-PVs), a majority of which infects anogenital mucosa, display the most dramatic reduction in TC dinucleotide content. Computer modeling and reconstruction of ancestral alpha-PV genomes suggest that TC depletion occurred after the alpha-PVs diverged from their most recent common ancestor. In addition, we found that TC depletion in alpha-PVs is greatly affected by protein coding potential. Taken together, our results suggest that PVs replicating in tissues with high APOBEC3 levels may have evolved to evade restriction by selecting for variants that contain reduced APOBEC3 target sites in their genomes.

A large spectrum of alpha and beta papillomaviruses are detected in human stool samples.

Human papillomaviruses (HPVs) have been detected in urban wastewaters, demonstrating that epitheliotropic viruses can find their way into sewage through the washing of skin and mucous membranes. Papillomavirus shedding through faeces is still an unexplored issue. The objective of the present study was to investigate the presence of HPVs in stool samples. We analysed 103 faecal specimens collected from hospitalized patients with diarrhoea using validated primers able to detect α, β and γ HPVs. PCR products underwent sequencing analysis and sequences were aligned to reference genomes from the Papillomavirus Episteme database. A total of 15 sequences were characterized from the faecal samples. Thirteen samples (12.6 %) were positive for nine genotypes belonging to the α and β genera: HPV32 (LR, α1), HPV39 (HR, α7), HPV44 (LR, α10), HPV8 (β1), HPV9, HPV23, HPV37, HPV38 and HPV120 (β2). Two putative novel genotypes of the β genus, species 1 and 2, were also detected. The tissue(s) of origin is unknown, since faeces can collect HPVs originating from or passing through the entire digestive system. To our knowledge, this is the first investigation on the occurrence and diversity of HPVs in faecal samples. Results from this study demonstrate that HPVs can find their way into sewage as a consequence of shedding in the faeces. This highlights the need for further studies aimed at understanding the prevalence of HPV in different water environments and the potential for waterborne transmission.

Identification of a Novel Human Papillomavirus by Metagenomic Analysis of Samples from Patients with Febrile Respiratory Illness

As part of a virus discovery investigation using a metagenomic approach, a highly divergent novel Human papillomavirus type was identified in pooled convenience nasal/oropharyngeal swab samples collected from patients with febrile respiratory illness. Phylogenetic analysis of the whole genome and the L1 gene reveals that the new HPV identified in this study clusters with previously described gamma papillomaviruses, sharing only 61.1% (whole genome) and 63.1% (L1) sequence identity with its closest relative in the Papillomavirus episteme (PAVE) database. This new virus was named HPV_SD2 pending official classification. The complete genome of HPV-SD2 is 7,299 bp long (36.3% G/C) and contains 7 open reading frames (L2, L1, E6, E7, E1, E2 and E4) and a non-coding long control region (LCR) between L1 and E6. The metagenomic procedures, coupled with the bioinformatic methods described herein are well suited to detect small circular genomes such as those of human papillomaviruses.

Mucosal and cutaneous human papillomaviruses detected in raw sewages.

Epitheliotropic viruses can find their way into sewage. The aim of the present study was to investigate the occurrence, distribution, and genetic diversity of Human Papillomaviruses (HPVs) in urban wastewaters. Sewage samples were collected from treatment plants distributed throughout Italy. The DNA extracted from these samples was analyzed by PCR using five PV-specific sets of primers targeting the L1 (GP5/GP6, MY09/MY11, FAP59/64, SKF/SKR) and E1 regions (PM-A/PM-B), according to the protocols previously validated for the detection of mucosal and cutaneous HPV genotypes. PCR products underwent sequencing analysis and the sequences were aligned to reference genomes from the Papillomavirus Episteme database. Phylogenetic analysis was then performed to assess the genetic relationships among the different sequences and between the sequences of the samples and those of the prototype strains. A broad spectrum of sequences related to mucosal and cutaneous HPV types was detected in 81% of the sewage samples analyzed. Surprisingly, sequences related to the anogenital HPV6 and 11 were detected in 19% of the samples, and sequences related to the “high risk” oncogenic HPV16 were identified in two samples. Sequences related to HPV9, HPV20, HPV25, HPV76, HPV80, HPV104, HPV110, HPV111, HPV120 and HPV145 beta Papillomaviruses were detected in 76% of the samples. In addition, similarity searches and phylogenetic analysis of some sequences suggest that they could belong to putative new genotypes of the beta genus. In this study, for the first time, the presence of HPV viruses strongly related to human cancer is reported in sewage samples. Our data increases the knowledge of HPV genomic diversity and suggests that virological analysis of urban sewage can provide key information useful in supporting epidemiological studies.

The Papillomavirus Episteme: a central resource for papillomavirus sequence data and analysis

The goal of the Papillomavirus Episteme (PaVE) is to provide an integrated resource for the analysis of papillomavirus (PV) genome sequences and related information. The PaVE is a freely accessible, web-based tool (http://pave.niaid.nih.gov) created around a relational database, which enables storage, analysis and exchange of sequence information. From a design perspective, the PaVE adopts an Open Source software approach and stresses the integration and reuse of existing tools. Reference PV genome sequences have been extracted from publicly available databases and reannotated using a custom-created tool. To date, the PaVE contains 241 annotated PV genomes, 2245 genes and regions, 2004 protein sequences and 47 protein structures, which users can explore, analyze or download. The PaVE provides scientists with the data and tools needed to accelerate scientific progress for the study and treatment of diseases caused by PVs.