L1 Gene Sequence Research Articles

HPV18 L1 and long control region sequences variation and E6/E7 differential expression in nasopharyngeal and cervical cancers: a comparative study

BackgroundThe role of high-risk human papillomaviruses (hr-HPVs) in cervical cancer (CC) pathogenesis has long been established. Knowledge about the involvement of hr-HPVs in the etiology of nasopharyngeal cancers (NPC) was not well appreciated until the early 2000s when a clear link began to emerge. However, it is not clear whether HPV oncogenesis in the different epithelial cancers is associated with L1 gene and long-control region (LCR) sequences variation. This study aimed to investigate the HPV18 L1 gene and LCR sequences variation in cervical and nasopharyngeal biopsies, and assessed E6 and E7 genes expression level in both cancers.MethodFour-hundred and three (403) formalin-fixed paraffin-embedded tissues originating from nasopharyngeal (NPC) (279) and cervical (CC) (124) sites were collected from a pathology laboratory, Pathologist Without Borders, Accra, Ghana. Haematoxylin and eosin staining was carried out to confirm the presence of cancer on prepared biopsy sections. DNA was extracted from the confirmed cancer biopsies, followed by PCR using MY09/GP5+ /6+ primers to detect the presence of HPV and specific primers for the amplification of L1 gene and LCR. Sanger sequencing was carried out to determine HPV genotypes, and L1 and LCR sequences variant of HPV18s in CC and NPC biopsies. The HPV18 E6/E7 mRNA expression pattern in both cancers was determined using RT-qPCR.ResultsMost of the NPC (45%) and CC (55%) biopsies were HPV18 positive. Comparison of HPV18 L1 sequences obtained from cervical and nasopharyngeal cancer tissues, the L1 sequences from the NPC were highly dissimilar with a 59–100% variation among themselves, and in relation to the reference strains. However, the L1 sequences from the CC were more similar with a 91.0–100% variation among the amplified sequences. Also, the LCR sequences from CC were quite different relative to that of NPC. Results for the differential expression of E6/E7 in the two cancers showed a higher fold change in E6 expression in the CC tissues than the NPC tissues while a reverse expression pattern was found for E7 gene.ConclusionThe current study reports for the first-time variations in HPV18 L1 and LCR sequences, and differential expression of E6/E7 genes in NPC compared to CC, suggesting a possible adaptation mechanism of the virus at different cancer sites.

Genetic variability analysis of humanpapillomavirus 58: Novel sublineage identification and persistent infection association.

This study aims to characterize the genetic variability of HPV58, identify novel lineages and sublineages, and explore the association between persistent/multiple HPV58 infections and genetic variation. In this study, samples from 124 women with HPV58 infection in Eastern China were collected and 81 isolates of E6 and L1 full-length genes were successfully amplified from 55 samples. We evaluated the diversity of genetic variants and performed correlation analyses between genetic variability and pathology, vaccination, multiple infections, and persistent infections. Among the E6 and L1 gene sequences collected, the dominant prevailing sublineages were A1 (46.2%) and A2 (23.1%). In addition, we found two potential novel sublineages denoted as the A4 and A5 sublineage. A total of 50 nucleotide substitutions, including 28 synonymous substitutions and 22 nonsynonymous substitutions, were observed in the E6 and L1 genes. Among them, variants with A388C/K93N substitutions in the E6 gene correlated with persistent infection (≥1 and ≥2 years) (p < 0.005), and C307T/C66C was associated with persistent infection (≥2 years) (p < 0.005). Notably, two mutations above were detected in the isolate from the patient with breakthrough vaccine infection. Our study found two novel sublineages and sites of genetic variability in multiple and persistent infection variants. In addition, we identified two mutational sites associated with persistent infection. This study provides new insight into the clinical characteristics of HPV 58 genetic variations and offers new ideas for research on next-generation vaccines in Eastern China.

Abortion storm of Yili horses is associated with Equus caballus papillomavirus 2 variant infection.

Nine different species of Equus caballus papillomavirus (EcPV) and three bovine papillomaviruses (BPVs) have been reported to infect horses. However, there are few descriptions of such infections in China. In our pioneer study on Chinese horses, we identified EcPV-2 in the nasal swabs (4/230, 1.7%) of Yili horses, and the semen (3/18, 16.7%) of thoroughbred horses. This indicated that EcPV is indeed hosted by horses in China, and that EcPV-2 might be transmitted though breeding. Further detection of EcPVs in the lung tissues of aborted fetuses of Yili horses, which were originally negative for equid herpes viruses, demonstrated EcPV-2 positivity in 19 of 50 samples, thereby indicating that EcPV-2 may be a new pathogen responsible for causing abortion. Thereafter, sequence analyses of the L1 genes of 26 EcPV-2 in China were performed, indicating that EcPV-2, which primarily infects horses in China, shared 98.3-99.9% nt identity with the published sequences for EcPV-2. These observations indicated that EcPV-2 identified in the current study were highly similar variants of the previously identified strains of EcPV-2. Phylogenetic analysis based on L1 gene sequences from GenBank showed that the EcPV-2 found in Chinese horses was closely related to and clustered together with an already known EcPV-2a lineage. Our study provides the first evidence related to EcPV-2 infection in Chinese horses, which can serve as a causative agent for Yili horse abortions, and may thus lay the foundation for a systematic and detailed epidemiological study of this infection in Chinese horses.

Characterization of a novel papillomavirus identified from a whale (Delphinapterus leucas) pharyngeal metagenomic library

Here, using viral metagenomic method, a novel whale papillomavirus (temporarily named wPV, GenBank accession number OP856597) was discovered in a whale (Delphinapterus leucas) pharyngeal metagenomic library. The complete genome size of wPV is 7179 bp, with GC content of 54.4% and a nucleotide composition of 23.4% A, 22.3% T, 28.4% G, and 25.9% C. The viral genome has a typical papillomavirus organization pattern, and five ORFs were predicted, including two late genes encoding L1 and L2, and three early genes encoding E1, E2, and E6. Pairwise sequence comparison and phylogenetic analysis based on the L1 gene sequence indicated that wPV may be a novel species within genus Dyodeltapapillomavirus. In addition, the E2 region of wPV was predicted to have a potential recombination event. The discovery of this novel papillomavirus increases our understanding of the viral ecology of marine mammals, providing insights into possible future infectious diseases.

Identification of Genetic Variants of Human Papillomavirus in a Group of Mexican HIV/AIDS Patients and Their Possible Association with Cervical Cancer.

Infections caused by the human immunodeficiency virus (HIV) and human papillomavirus (HPV) cause thousands of deaths worldwide each year. So far, there has been no consensus on whether there is a direct relationship between the incidence of neoplasms and the immunosuppression caused by HIV that could help understand if coinfection increases the likelihood of cervical cancer. The objective of the study was to identify the presence of genetic variants of HPV in a group of HIV-positive women and their possible association with cervical cancer. Cervical samples were taken from HIV-positive patients for cytological analysis to identify the HPV genotype by polymerase chain reaction (PCR) and sequencing. The most prevalent L1 capsid protein mutations in the HPV genotype were analyzed in silico. Various types of HPV were identified, both high-risk (HR) and low-risk (LR). The most prevalent genotype was HPV51. Analysis of the L1 gene sequences of HPV51 isolates showed nucleotide variations. Of the samples analyzed in Puebla, Mexico, HPV51 had the highest incidence (17.5%, 7/40). Different mutations, which could be used as population markers, were detected in this area, and they have not been reported in the L1 databases for HPV51 in Mexico. Genotypes 6, 14, 86, 87, 89, and 91, not detected or reported in samples from patients with HPV in Mexico, were also identified. Data from the population analyzed suggest no direct relationship between HIV immunosuppression and cervical cancer, regardless of the high- or low-risk HPV genotype. Furthermore, it is possible to develop regional population markers for the detection of HPV based on the mutations that occur in the sequence of nucleotides analyzed.

Prevalence of Human Papillomavirus Genotypes in Pterygia from Thai Individuals

ABSTRACT Purpose Pterygium, a common ocular growth, has an unknown pathogenesis and aetiology. Environmental factors such as ultraviolet light, genetic factors and viral infections may be implicated in the development of pterygia. Human papillomavirus (HPV), an oncogenic virus, has previous been detected in individuals with pterygia. The aim of this study was to assess the prevalence of HPV genotypes in pterygia from Thai individuals. Methods DNA was extracted from 389 pterygia. HPV was detected by nested PCR and HPV genotyping was conducted using reverse hybridization. The DNA sequences of HPV-L1 genes were analyzed. Results HPV was detected in only 6.8% (25/389) of pterygia from Thai individuals. The majority (16/25, 64%) of strains were genotyped as HPV-16 and the remainder (9/25, 36%) could not be typed. Four pterygia showed evidence of coinfection by HPV-16 and either HPV-18 (2/25, 8%) or HPV-58 (2/25, 8%). Nine of 11 samples showed the same HPV-16 L1 gene sequence that was identical to a HPV-16 reference sequence in GenBank. The remaining two samples each bore silent single nucleotide mutations (T1078G and T1081A) that did not result in amino acid changes. Conclusion HPV, especially HPV-16, may be one of the pathogens causing pterygia in Thai individuals. Genotyping data suggested that HPV-16 from pterygia may be similar in sequence to HPV-16 causing cervical cancer.

Worldwide genetic variations in high-risk human papillomaviruses capsid L1 gene and their impact on vaccine efficiency.

Human papillomavirus is the most common sexually transmitted infection. It is associated with different cancers, mainly cervical cancer, which remains the fourth most frequent cancer among women worldwide; it is also related to anogenital (anus, vulvar, vagina, and penis) and oropharyngeal cancers. Vaccination against HPV infection is the major way of prevention, and it has demonstrated impressive efficacy in reducing cervical cancer incidence. Nowadays, all the licensed HPV recombinant vaccines were designed based on HPV major capsid L1 protein. However, some variations in the HPV L1 gene sequence may induce structural changes within the L1 protein, which may alter the affinity and interaction of monoclonal antibodies (MAbs) with L1 protein epitopes, and influence host immune response and recognition. Hence, the importance of accuracy in delineating epitopes relevant to vaccine design and defining genetic variations within antigenic regions in the L1 gene to predict its impact on prophylactic vaccine efficiency. The present review reports the sequence variations in HR-HPV L1 gene isolates from different countries around the world, which may help to understand the effect of HPV L1 gene variations on vaccine efficiency. Research studies were retrieved from PubMed, Google Scholar, Science direct, and the National Center for Biotechnology Information (NCBI) database. A total of 31 articles describing genetic variations within the major capsid L1 gene and conducted in Africa, Europe, America and Asia were found. Only 26 studies conducted on HPV16, 18, 31, 33, 58, 45 and 52 which are the targets of HPV prophylactic vaccines, and which reported genetic variations within the L1 gene, were selected and evaluated in this review. We found a total of 87, 49, 11, 7, 22, 3, and 17 non-synonymous single nucleotide polymorphisms (SNPs) within HPV16, HPV18, HPV31, HPV58, HPV45, and HPV52 L1 gene, respectively. Four mutations were frequently observed in HPV16 L1 sequences: T353P in the HI loop, H228D in the EF loop, T266A in the FG loop, and T292A in the FG loop. Two mutations in HPV58 L1 sequences: T375N in the HI loop and L150F in the DE loop. Three mutations in HPV33 L1 sequences: T56N in the BC loop, G133S in the DE loop, T266K in the FG loop. Other mutations were found in HPV18, HPV45, and HPV52 L1 sequences. Some were found in different countries, and others were specific to a given population. Furthermore, some variations were located on peptide binding epitopes and lead to a modification of epitopes, which may influence MAbs interactions. Others need further investigations due to the lack of studies. This study investigated the major capsid L1 genetic diversity of HPV16, 18, 31, 33, 58, 45, and 52 circulating in different populations around the world. Further investigations should be conducted to confirm their effect on immunogenicity and prophylactic vaccine efficiency.

Molecular Characterization and Developing a Point-of-Need Molecular Test for Diagnosis of Bovine Papillomavirus (BPV) Type 1 in Cattle from Egypt.

Simple SummaryBovine papillomatosis is a disease caused by bovine papillomavirus (BPV), which is a diverse group of oncogenic viruses that challenge cattle industry, resulting in significant economic losses. The present study investigated the occurrence of bovine papillomatosis among cattle (n = 308) with cutaneous warts on the head and neck from New valley Province, Egypt through molecular detection of BPV-1, -2, -4, -5, and -10. The work also involved a phylogenetic analysis of the positive samples for detection of the genetic relatedness of the virus. Interestingly, BPV-1 DNA was detected in 84.6% of the collected samples. Furthermore, the study included the development of an isothermal nucleic acid amplification test, which is a field test combining molecular and lateral flow immunoassays for point-of-need testing appropriate for veterinary use in resource-limited settings. Collectively, our study provided interesting data related to the combined use of molecular and immunoassays methods in the detection of the virus besides better understanding the genetic relatedness of the circulating genotypes of BPV-1 in Egypt. Our study suggested further research to explore more about the other genotypes of BPV in the Egyptian environment that could be helpful for the implementation of control strategies for combating this disease.Bovine papillomatosis is a viral disease of cattle causing cutaneous warts. A diagnosis of this viral infection is very mandatory for combating the resulting economic losses. Given the limited data available about bovine papillomavirus (BPV) in Egypt, the present study involved the molecular diagnosis of bovine papillomavirus type-1 (BPV-1), -2, -4, -5, and -10 in cattle presenting cutaneous warts on the head and neck from New Valley Province, Egypt. The phylogenetic analysis of the detected types of BPV was also performed, followed by developing a point-of-need molecular assay for the rapid identification of identified BPV types. In this regard, a total of 308 cattle from private farms in Egypt were clinically examined, of which 13 animals presented cutaneous warts due to suspected BPV infection. The symptomatic animals were treated surgically, and biopsies from skin lesions were collected for BPV-1, -2, -4, -5, and -10 molecular identification using polymerase chain reaction (PCR). The presence of BPV-1 DNA was confirmed in 11 collected samples (84.6%), while BPV-2, -4, -5, and -10 were not detected. Sequencing of the PCR products suggested the Egyptian virus is closely related to BPV found in India. An isothermal nucleic acid amplification test (NAAT) with labeled primers specific for the BPV-1 L1 gene sequence, and based on recombinase polymerase amplification (RPA), in combination with a lateral flow strip assay for the detection of RPA products, was developed and tested. The point-of-need molecular assay demonstrated a diagnostic utility comparable to PCR-based testing. Taken together, the present study provides interesting molecular data related to the occurrence of BPV-1 in Egypt and reveals the genetic relatedness of the Egyptian BPV-1 with BPV-1 found in buffalo in India. In addition, a simple, low-cost combined test was also validated for diagnosis of the infection. The present study suggests the necessity of future investigations about the circulating strains of the virus among the cattle in Egypt to assess their genetic relatedness and better understand the epidemiological pattern of the disease.

Genomic sequencing of a virus representing a novel type within the species Dyopipapillomavirus 1 in an Indian River Lagoon bottlenose dolphin.

Fecal samples collected from free-ranging Atlantic bottlenose dolphins (BDs) in the Indian River Lagoon of Florida were processed for viral discovery using a next-generation sequencing (NGS) approach. A 693-bp contig identified in the NGS data was nearly identical to the partial L1 gene sequence of a papillomavirus (PV) previously found in a penile papilloma in a killer whale (Orcinus orca). Based on this partial bottlenose dolphin papillomavirus (BDPV) sequence, a nested inverse PCR and primer-walking strategy was employed to generate the complete genome sequence. The full BDPV genome consisted of 7299bp and displayed a typical PV genome organization. The BDPV E6 protein contained a PDZ-binding motif, which has been shown to be involved in carcinogenic transformation involving high-risk genital human PVs. Screening of 12 individual fecal samples using a specific endpoint PCR assay revealed that the feces from a single female BD displaying a genital papilloma was positive for the BDPV. Genetic analysis indicated that this BDPV (Tursiops truncatus papillomavirus 8; TtPV8) is a new type of Dyopipapillomavirus 1, previously sequenced from an isolate obtained from a penile papilloma in a harbor porpoise (Phocoena phocoena). Although only a partial L1 sequence has been determined for a PV detected in a killer whale genital papilloma, our finding of a nearly identical sequence in an Atlantic BD may indicate that members of this viral species are capable of host jumping. Future work is needed to determine if this virus is a high-risk PV that is capable of inducing carcinogenic transformation and whether it poses a significant health risk to wild delphinid populations.

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.

Phylogenetic Analysis of Human Papillomavirus 16 and 52 L1 Gene from Cervical Cancer in Bandung, Indonesia

BACKGROUND: Chronic infection with high-risk type of human papillomavirus (HPV) can cause cervical cancer. Previous studies showed that multiple infections of HPV are found in cervical cancer caused by multiple HPV infections and the most common are HPV-16 and HPV-52. The origin of HPV-16 circulating in Indonesia varies. The purpose of this study was to explore the origin of multiple infections of HPV-52 and HPV-16 in cervical cancer by using a phylogenetic tree.METHODS: During July-November 2010, 100 women were diagnosed with cervical cancer in the Department of Obstetrics and Gynecology, Dr. Hasan Sadikin General Hospital, Bandung, Indonesia. Only 96 patients were involved in this study. Ninety-six samples of HPV deoxyribonucleic acid (DNA) were isolated from biopsied tissue of cervical cancer. Multiple infections of HPV genotypes HPV-16 and HPV-52 were confirmed by using the linear assay for HPV genotyping test. Afterward,HPV-52L1 gene was amplified by using self-designed primer. L1 gene was also sequenced and analyzed using phylogenetic program (MEGA6.06).RESULTS: The result of phylogenetic tree construction showed that isolated HPV-52 originated from multiple infections of HPV-16 and HPV-52 from cervical cancer patients in Bandung were in a subgroup with isolates originating from EU077219 Canada (America) and KT799980 southwest China (Asia). Isolate HPV-16 in one subgroup with isolates originating from KU951191.1 (Southwest China).CONCLUSION: L1 gene sequence from multiple infections isolated from HPV-16 and HPV-52 from cervical cancer patients in Bandung refers to the variation of L1 gene reported from Canada and southwest China. This proves that Indonesia’s HPV clusters are located in the strains found in America and Asia.KEYWORDS: multiple infections, HPV-16, HPV-52, L1 gene, phylogenetic

Molecular Characterization of Bovine Papillomaviruses Obtained From Cattle Cutaneous Warts in Nigeria

Bovine papillomatosis is a benign, self-limiting tumor in cattle that is caused by a wide diverse of bovine papillomaviruses (BPVs 1-23). The disease is of high economic importance to the beef and hide industries in Nigeria. To date, there is no information on the molecular characterization of BPVs in Nigeria. Consequently, we collected cutaneous wart lesions from 31 cows in Ibadan, Nigeria and screened them for BPV using a degenerate primer that amplify partial L1 gene sequence of the virus in PCR. Six samples (19.4%) were positive for bovine papillomavirus. Thereafter, amplicons were purified and sequenced. BLAST search, pairwise sequence alignments, pairwise genetic distance calculations, multiple sequence alignment and phylogenetic analysis revealed the viruses to be bovine papillomavirus-1 belonging to the genus Deltapapillomavirus.

Colorimetric Detection of 23 Human Papillomavirus Genotypes by Loop-Mediated Isothermal Amplification.

Human papillomavirus (HPV) infection is linked to cervical cancer. With the technological development of molecular biology and epidemiology, detection and treatment of HPV has become an important mean to prevent cervical cancer. A simple, rapid, and sensitive colorimetric loop-mediated isothermal amplification (LAMP) method was established herein to detect 23 HPV genotypes. The sequences of the primers for the LAMP reaction were located in the L1 gene of the HPV genome. As it is a fluorescent dye, calcein was added before the reaction. The reaction was run under isothermal conditions at 65°C for 40 minutes. A positive reaction was indicated by a color change from yellow to fluorescent green. The fluorescence curve diagram represents the monitoring of real time quantitative instrument. 450 cervical swab samples from patients with single infections of 23 different HPV genotypes were examined to evaluate the specificity. The results revealed no cross-reaction with other HPV genotypes. A serial dilution of a cloned plasmid containing 23 HPV L1 gene sequences was employed to evaluate the sensitivity. Different HPV subtypes have different detection capability. The sensitivity of different HPV subtypes tested by LAMP assay was in the range from 1.0 x10 to 4.0 x 103 copies per reaction. The LAMP assay and the RDB (reverse dot blot) were compared for detecting and genotyping HPV among the 450 clinical samples. There were 385 (85.6%) and 375 (83.3%) HPV positive specimens detected by LAMP and RDB, respectively, as well as 306 (68.0%) and 296 (65.8%) for HR-HPV positive specimens. The agreement between the LAMP and RDB assays was 93.3% (κ = 0.75) for HPV positivity and 94.7% (κ = 0.88) for HR-HPV positivity. It was concluded that this colorimetric LAMP assay had potential application for the rapid screening of the HPV infection in resource-limited hospitals or rural clinics.

Cloning, Expression and Bioinformatic Analysis of Human Papillomavirus Type 52 L1 Capsid Gene from Indonesian Patient

Human papillomavirus (HPV) type 52 is the most prevalent type for causing cervical cancer in Indonesian population. Cervical cancer becomes the most common cancer suffered by Indonesian women. Prevention of HPV infection can be achieved using HPV virus-like particle (VLP) vaccine derived from L1 major capsid protein. This study aimed to clone and analyze HPV-52 L1 gene. DNA obtained from biopsy of a cervical cancer patient was amplified using specific primers designed from Asian originated HPV-52 L1 gene available in the GenBank. The isolated HPV-52 L1 gene sequence was submitted to GenBank with accession number [KF225497]. Expression of HPV-52 L1 gene was performed using pRSET/EmGFPEscherichia coli expression vector. We analyzed and compared the HPV-52 L1 gene expressions from recombinant E.coli BL21 (DE3) that had been induced for 3 hours with 1 mM IPTG and without induction. The protein was expressed in insoluble form. We performed the following bioinformatic analyses: construction of phlyogenetic tree, T-cell epitopes prediction and 3D proteins structure modelling. We utilized the following softwares: MEGA5 for phylogenetic tree, IEDBann for MHC prediction, CLC DNA Workbench 6.5 for hydrophobicity analysis and PDB-Viewer Deep for 3D protein structure analysis. The phylogenetic tree which was developed based on [KF225497] sequence showed that it shared a branch with Asian countries (Philippines and Thailand). The deduced amino acid sequences of the predicted epitopes that were consistent in all of the programs were 259GTLGDPVPGDLYIQGS274 and 345KKESTYKNE353. This information may be useful to design diagnostic strategies and vaccine suitable for Indonesian population.

First report of a variant bovine papillomavirus type 2 (BPV-2) in cattle in the Iberian Peninsula.

Infections caused by bovine papillomavirus (BPV) have been described worldwide. Some types, like BPV-1 and BPV-2, have been reported in association with skin warts and fibropapillomas in cattle and sarcoids in equids. In this study we have investigated the presence of BPV in cutaneous warts isolated from a steer in Spain. Cutaneous fibropapillomatosis was confirmed by histopathological analysis. Complete genome was amplified by multiple-primed rolling circle and the L1, E5 and E6 genes were sequenced. The isolate was classified as a variant of BPV-2 on the basis of the L1 gene sequences. Genetic variability of L1, E5 and E6 genes was compared with BPV-2 isolates from different hosts in several continents. Some mutations involved non-synonymous substitutions when compared to the prototype strain. One of these non-conservative mutations was located in the jelly roll β-barrel of the EF loop of the capsid protein (encoded by L1). This study presents the first report of a variant of BPV-2 infection in the Iberian Peninsula and contributes to extend the knowledge of the spreading and circulation of BPV.

Entropy-based approach for selecting informative regions in the L1 gene of bovine papillomavirus for phylogenetic inference and primer design

Bovine papillomaviruses (BPVs) cause many benign and malignant lesions in cattle and other animals. Twelve BPV types have been identified so far, and several putative novel BPV types have been detected based on the analysis of L1 gene fragments, generated by FAP59/64 and MY11/09 primers. Phylogenetic trees are important in studies that describe novel BPV types. However, topological mistakes could be a problem in such studies. Therefore, we made use of entropy to find phylogenetic informative regions in the BPV L1 gene sequences from all 12 BPVs. Six data sets were created and phylogenetically compared to each other using neighbor-joining and maximum likelihood methods of phylogenetic tree reconstruction. We found two major regions in the L1 gene, using an entropy-based approach, which selects regions with low information complexity. More robust phylogenetic trees were obtained with these regions, when compared to the ones obtained with FAP59/64 and MY11/09 primers. More robust phylogenetic trees are important to accurately position novel BPV types, subtypes and variants. We conclude that an entropy-based approach is a good methodology for selecting regions of the L1 gene of BPVs that could be used to design more specific and sensitive degenerate primers, for the development of improved diagnostic methods.

Novel animal papillomavirus sequences and accurate phylogenetic placement

All amniotes are probably infected by specific papillomaviruses (PVs), but knowledge about PV diversity remains sparse. An insufficient taxon sampling, and a focus on humans as hosts, may perturb phylogenetic analyses leading to wrong conclusions about PV evolution. We performed a systematic approach to explore the diversity of PVs combining rolling circle amplification with the use of “universal” primers to search for the presence of novel PV sequences in animal samples. We communicate 12 sequences putatively corresponding to novel PVs gained from 10 host species in eight mammal families: Bovidae, Canidae, Cervidae, Equidae, Hominidae, Phocoenidae, Procyonidae and Pteropodidae. The phylogenetic position of the new sequences was inferred with an evolutionary placement algorithm under a Maximum Likelihood framework using a pre-computed, well-resolved tree constructed with the E1–E2–L1 gene sequences as a backbone. The new sequences were phylogenetically diverse and could be respectively placed with confidence within all four PV crown groups. The prevailing presence of sequences from the crown groups Alpha+Omikron-PVs and Beta+Xi-PVs may correspond to an increased viral diversity in these taxa, or rather reflect a combination of anthropocentric bias and preferential amplification from commonly used “universal” primers. Our results combined with literature data support the view that the number and diversity of animal PVs is overwhelmingly large.

Identification of a Novel Bat Papillomavirus by Metagenomics

The discovery of novel viruses in animals expands our knowledge of viral diversity and potentially emerging zoonoses. High-throughput sequencing (HTS) technology gives millions or even billions of sequence reads per run, allowing a comprehensive survey of the genetic content within a sample without prior nucleic acid amplification. In this study, we screened 156 rectal swab samples from apparently healthy bats (n = 96), pigs (n = 9), cattles (n = 9), stray dogs (n = 11), stray cats (n = 11) and monkeys (n = 20) using a HTS metagenomics approach. The complete genome of a novel papillomavirus (PV), Miniopterus schreibersii papillomavirus type 1 (MscPV1), with L1 of 60% nucleotide identity to Canine papillomavirus (CPV6), was identified in a specimen from a Common Bent-wing Bat (M. schreibersii). It is about 7.5kb in length, with a G+C content of 45.8% and a genomic organization similar to that of other PVs. Despite the higher nucleotide identity between the genomes of MscPV1 and CPV6, maximum-likelihood phylogenetic analysis of the L1 gene sequence showed that MscPV1 and Erethizon dorsatum papillomavirus (EdPV1) are most closely related. Estimated divergence time of MscPV1 from the EdPV1/MscPV1 common ancestor was approximately 60.2–91.9 millions of years ago, inferred under strict clocks using the L1 and E1 genes. The estimates were limited by the lack of reliable calibration points from co-divergence because of possible host shifts. As the nucleotide sequence of this virus only showed limited similarity with that of related animal PVs, the conventional approach of PCR using consensus primers would be unlikely to have detected the novel virus in the sample. Unlike the first bat papillomavirus RaPV1, MscPV1 was found in an asymptomatic bat with no apparent mucosal or skin lesions whereas RaPV1 was detected in the basosquamous carcinoma of a fruit bat Rousettus aegyptiacus. We propose MscPV1 as the first member of the novel Dyolambda-papillomavirus genus.

Multiplex PCR 기법을 이용한 보통사마귀 내 인유두종바이러스 검출 및 분류

A number of epidemiological studies have identified human papillomavirus (HPV) types 1, 2, 3, 4, 7, 10, 27, 57, and 65 in cutaneous common warts. However, identification of the HPV subtype by conventional polymerase chain reaction (PCR) is time consuming with its multi-step laboratory process. In this study, we aim to develop a specific one-step multiplex polymerase chain reaction method which capably identifies six different HPV genotypes related to common warts. By HPV DNA sequence analysis, 6 pairs of specific primers were designed from the intergenic regions of genes L1 to E6, and from genes E2 to L2. DNA sequence analysis with the L1 gene sequence of the sample was performed to measure the specificity of multiplex PCR. HPV-1, -2, -3, -4, -27, and -57 were identified without cross amplification in 109 out of 129 samples. The sensitivity and specificity of our set of primers in detecting HPV were 85% and 99.5%, respectively. For the 20 samples where HPV type was not identifiable by our batch of primer sets, multiplex PCR with an additional set of HPV primers was done, where 7 were found positive for HPV-7 or -65. Our results demonstrate that the newly designed multiplex PCR can rapidly detect the specific HPV subtype involved in common warts with high accuracy.

Genetic variability in the major capsid L1 protein of human papillomavirus type 16 (HPV-16) and 18 (HPV-18)

HPV-16 and HPV-18 infections result in nearly 73% of cervical cancers worldwide. The L1 protein comprising HPV vaccine formulations elicit high-titre neutralizing antibodies. The aim of this study was to detect L1 HPV-16 and HPV-18 gene polymorphisms and analyze intratypic variations. HPV-16 (n=29) and HPV-18 (n=5) L1 gene sequences were obtained from cervical samples harvested from Italian women. Phylogenetic trees were constructed using the Neighbor-Joining and the Kimura 2-parameters methods (MEGA software). To estimate selection pressures acting on the L1 gene, codon-specific non-synonymous (dN) and synonymous (dS) substitutions were inferred using the Nei–Gojobori method and Jukes–Cantor model (MEGA software) and integrated analyses carried out using SLAC, FEL and REL methodologies. All the HPV-16 L1 sequences analyzed fell into the European branch (99.4–99.7% similarity). Thirty-four single nucleotide changes were observed and 18 (52.9%) were non-synonymous mutations (7/18 were identified in sequences encoding an immunodominant loop and one occurred in the sequence encoding the α-4 domain associated with VLP conformation). There was no evidence of positive selection in the sequence alignment of L1 HPV-16 genes (P-value<0.1). One mutation was identified in a negatively selected codon. HPV-18 L1 analyzed sequences fell into two phylogenetic branches: the HPV-18 European branch (99.5–100% similarity) and the HPV-18 African branch (99.8% similarity). Nine single nucleotide changes were observed and 4/9 (44.5%) of these nucleotide mutations were non-synonymous and one was present in a sequence encoding the immunodominant FG loop. There was no evidence of positive selection in the sequence alignment of L1 HPV-18 genes (P-value<0.1). This study identified polymorphisms of undefined biological activity in HPV-16 and HPV-18 L1 sequences. Information regarding the genetic diversity of HPV-16 and HPV-18 L1 gene sequences may help define the oncogenic potential of respective strains and to better understand immune escape mechanisms.