Complete nucleotide sequence, genomic organization and phylogenetic analysis of a novel genital human papillomavirus type, HLT7474-S

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.

Cross-roads in the classification of papillomaviruses

Acceptance of an official classification for the family Papillomaviridae based purely on DNA sequence relatedness, was achieved as late as 2003. The rate of isolation and characterization of new papillomavirus types has greatly depended on and subjected to the development of new laboratory techniques. Introduction of every new technique led to a temporarily burst in the number of new isolates. In the following, the bumpy road towards achieving a classification system combined with the controversies of implementing and accepting new techniques will be summarized. An update of the classification of the 170 human papillomavirus (HPV) types presently known is presented. Arguments towards the implementation of metagenomic sequencing for this rapidly growing family will be presented.

Human genomic diversity, viral genomics and proteomics, as exemplified by human papillomaviruses and H5N1 influenza viruses

The diversity of hosts, pathogens and host-pathogen relationships reflects the influence of selective pressures that fuel diversity through ongoing interactions with other rapidly evolving molecules in the environment. This paper discusses specific examples illustrating the phenomenon of diversity of hosts and pathogens, with special reference to human papillomaviruses and H5NI influenza viruses. We also review the influence of diverse host-pathogen interactions that determine the pathophysiology of infections, and their responses to drugs or vaccines.

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

Human papillomavirus type 18 (HPV18) and HPV45 account for approximately 20% of all cervix cancers. We show that HPV18, HPV45, and the recently discovered HPV97 comprise a clade sharing a most recent common ancestor within HPV alpha7 species. Variant lineages of these HPV types were classified by sequence analysis of the upstream regulatory region/E6 region among cervical samples from a population-based study in Costa Rica, and 27 representative genomes from each major variant lineage were sequenced. Nucleotide variation within HPV18 and HPV45 was 3.82% and 2.39%, respectively, and amino acid variation was 4.73% and 2.87%, respectively. Only 18 nucleotide variations, of which 10 were nonsynonymous, were identified among three HPV97 genomes. Full-genome comparisons revealed maximal diversity between HPV18 African and non-African variants (2.6% dissimilarity), whereas HPV18 Asian-American [E1 (AA)] and European (E2) variants were closely related (less than 0.5% dissimilarity); HPV45 genomes had a maximal difference of 1.6% nucleotides. Using a Bayesian Markov chain Monte Carlo (MCMC) method, the divergence times of HPV18, -45, and -97 from their most recent common ancestors indicated that HPV18 diverged approximately 7.7 million years (Myr) ago, whereas HPV45 and HPV97 split off around 5.7 Myr ago, in a period encompassing the divergence of the great ape species. Variants within the HPV18/45/97 lineages were estimated to have diverged from their common ancestors in the genus Homo within the last 1 Myr (<0.7 Myr). To investigate the molecular basis of HPV18, HPV45, and HPV97 evolution, regression models of codon substitution were used to identify lineages and amino acid sites under selective pressure. The E5 open reading frame (ORF) of HPV18 and the E4 ORFs of HPV18, HPV45, and HPV18/45/97 had nonsynonymous/synonymous substitution rate ratios (d(N)/d(S)) over 1 indicative of positive Darwinian selection. The L1 ORF of HPV18 genomes had an increased proportion of nonsynonymous substitutions (4.93%; average d(N)/d(S) ratio [M3] = 0.3356) compared to HPV45 (1.86%; M3 = 0.1268) and HPV16 (2.26%; M3 = 0.1330) L1 ORFs. In contrast, HPV18 and HPV16 genomes had similar amino acid substitution rates within the E1 ORF (2.89% and 3.24%, respectively), while HPV45 E1 was highly conserved (amino acid substitution rate was 0.77%). These data provide an evolutionary history of this medically important clade of HPVs and identify an unexpected divergence of the L1 gene of HPV18 that may have clinical implications for the long-term use of an L1-virus-like particle-based prophylactic vaccine.

Classification of papillomaviruses

One hundred eighteen papillomavirus (PV) types have been completely described, and a yet higher number of presumed new types have been detected by preliminary data such as subgenomic amplicons. The classification of this diverse group of viruses, which include important human pathogens, has been debated for three decades. This article describes the higher-order PV taxonomy following the general criteria established by the International Committee on the Taxonomy of Viruses (ICTV), reviews the literature of the lower order taxa, lists all known "PV types", and interprets their phylogenetic relationship. PVs are a taxonomic family of their own, Papillomaviridae, unrelated to the polyomaviruses. Higher-order phylogenetic assemblages of PV types, such as the "genital human PVs", are considered a genus, the latter group, for example, the genus "Alpha-Papillomavirus". Lower-order assemblages of PV types within each genus are treated as species because they are phylogenetically closely related, but while they have distinct genomic sequences, they have identical or very similar biological or pathological properties. The taxonomic status of PV types, subtypes, and variants remains unchanged and is based on the traditional criteria that the sequence of their L1 genes should be at least 10%, 2-10%, and maximally 2% dissimilar from one another.

The novel human MOST-1 (C8orf17) gene exhibits tissue specific expression, maps to chromosome 8q24.2, and is overexpressed/amplified in high grade cancers of the breast and prostate

AIMS

To elucidate genes that participate in the process of oncogenesis, primers based on the E6 genes of genital human papillomaviruses (HPVs) were used to amplify potential expressed sequence tags (ESTs) from the MOLT-4 T lymphoblastic leukaemia cell line.

METHODS

Using the polymerase chain reaction (PCR) with human papillomavirus E6 gene primers, an EST from the MOLT-4 T lymphoblastic leukaemia cell line was amplified. Via rapid amplification of cDNA ends (RACE) and cycle sequencing from MOLT-4 and fetal lung cDNA libraries, overlapping cDNAs of 2786 bp and 2054 bp of the corresponding novel human intronless gene designated MOST-1 (for MOLT-4 sequence tag-1) were characterised and assigned the symbol C8orf17 by the HUGO Nomenclature Committee.

RESULTS

Both cDNAs contained a potential open reading frame (ORF) of 297 bp incorporating a methionine codon with an ideal Kozak consensus sequence for translation initiation, and encoding a putative hydrophilic polypeptide of 99 amino acids. Although reverse transcription PCR (RT-PCR) demonstrated MOST-1 expression in all 19 cancer and two normal cell lines tested, differential expression was seen in only nine of 16 normal tissues tested (heart, kidney, liver, pancreas, small intestine, ovary, testis, prostate, and thymus). A 388 bp fragment was amplified from the NS-1 mouse myeloma cell line, the sequence of which was identical to that within the MOST-1 ORF. The MOST-1 gene was mapped by fluorescent in situ hybridisation to chromosome 8q24.2, a region amplified in many breast cancers and prostate cancers, which is also the candidate site of potential oncogene(s) other than c-myc located at 8q24.1. Analysis of paired biopsies of invasive ductal breast cancer and adjacent normal tissue by semiquantitative and real time RT-PCR revealed average tumour to normal ratios of MOST-1 expression that were two times greater in grade 3 cancers than in grade 1 and 2 cancers. Quantitative real time PCR of archival prostatic biopsies displayed MOST-1 DNA values that were 9.9, 7.5, 4.2, and 1.4 times higher in high grade carcinomas, intermediate grade carcinomas, low grade carcinomas, and benign hyperplasias, respectively, than in normal samples.

CONCLUSIONS

These data suggest a role for MOST-1 in cellular differentiation, proliferation, and carcinogenesis.

Complete sequence analyses of enterovirus 71 strains from fatal and non-fatal cases of the hand, foot and mouth disease outbreak in Singapore (2000)

Enterovirus 71 (EV71) is a major aetiological agent of hand, foot and mouth disease (HFMD). In recent years, several outbreaks in East Asia were associated with neurological complications and numerous deaths. An outbreak in Singapore in October 2000 afflicted thousands of children, resulting in four fatal cases from three of whom EV71 was isolated. The genomes of two representative EV71 strains isolated from a fatal case and a surviving patient were completely sequenced, and their nucleotide and amino acid sequences compared with known EV71 strains. The two outbreak strains were classified under genogroup B, together with those previously isolated in Singapore, Malaysia and Japan. Comparative sequence analysis of the two Singapore strains revealed 99% nucleotide similarity, while their deduced amino acid sequences were almost identical except for residue 1506 in the 3A non-structural region. Given that the outbreak involved closely related genetic variants of EV71, the broad spectrum of disease severity may be attributed to critical factors such as varying viral inoculation doses or differing host immune responses following infection, but is less likely to be due to the emergence of EV71 strains with heightened virulence.

Differential display RT-PCR analysis of enterovirus-71-infected rhabdomyosarcoma cells reveals mRNA expression responses of multiple human genes with known and novel functions

In order to better understand cellular responses to viral infection at the transcriptional level, we employed differential display RT-PCR to analyze mRNAs from RD rhabdomyosarcoma cells following infection with a neurovirulent enterovirus 71 (EV71) strain, compared with mRNAs from uninfected cells. Of 250 expressed sequence tags (ESTs) isolated, sequenced, and identified, all were of cellular origin except 1 that was of viral origin. Of these, 156 were individual distinctive clones, comprising 45 mRNAs showing unaltered expression and 111 mRNAs exhibiting upregulation or downregulation. Of the 45 uniformly expressed mRNAs, 14 represented unknown genes. Of the 111 differentially expressed mRNAs, 63 did not match any known genes. Forty-eight of the 111 mRNAs modified by EV71 infection matched known genes, including those encoding components of cell cycle, cytoskeleton, and cell death mediators; protein degradation mediators; mitochondrial-related proteins; components of protein translation and modification; and cellular transport proteins. The altered expression profiles of representative genes were authenticated by semiquantitative RT-PCR and real-time RT-PCR. We also identified a novel alternatively spliced transcript of TRIP7 thyroid receptor interactor protein; the putative human homolog of murine mc7 mRNA predominantly expressed in the brain; and a novel mRNA similar to that encoding vacuolar protein 8 involved in protein targeting. These results underscore the applicability of the mRNA differential display technique for elucidating the expression profiles of known and even novel genes in response to cellular infection with pathogenic viruses.

Identification of six putative novel human papillomaviruses (HPV) and characterization of candidate HPV type 87

Six putative novel human papillomavirus (HPV) types were detected by using general primers for a conserved L1 HPV region in patients examined in gynecologic centers. One of the isolates, detected in samples from 4 patients with koilocytic atypia at cervical cytology (3 of whom were also infected with human immunodeficiency virus type 1), was completely sequenced, identified as a new HPV genotype, and designated candidate HPV87 (candHPV87) by the Reference Center for Human Papillomavirus. candHPV87 shows the classic HPV genome organization and the absence of a functional E5 coding region. Phylogenetic analysis documented that the candHPV87 genome clusters within the A3 group of HPVs, together with HPV61, HPV72, HPV83, HPV84 and candHPV86, which have been completely sequenced, and a number of other putative novel genotypes (two of which are described in this work), which have been partially characterized. To address the growth-enhancing potential of candHPV87, the E6 and E7 putative coding regions were cloned and expressed in tissue cultures. The data indicate that both proteins stimulate cell division in tissue cultures more than those of low-risk HPVs, though not as much as those of HPV16. Taken together, the clinical, molecular, and biological data suggest that the novel papillomavirus characterized in the present study is a low- to intermediate-risk HPV.

Characterization of a novel genital human papillomavirus by overlapping PCR: candHPV86 identified in cervicovaginal cells of a woman with cervical neoplasia

A novel human papillomavirus (HPV), candHPV86, was cloned and characterized from cervicovaginal cells obtained from a 37-year-old Hispanic woman with cervical intraepithelial neoplasia grade 1 (CIN1) using an overlapping PCR technique. Primers were designed by phylogenetic alignment of closely related HPV genomes using the L1 fragment sequence amplified by GP5+/6+. The 7983 bp complete nucleotide sequence of the HPV genome was determined by sequence walking. A basic local alignment sequence tool (BLAST) homology search using the L1 open reading frame demonstrated that this HPV was most closely related to HPVHAN2294 (GenBank, AJ400628; 86% homology) and HPV84 (84% homology). candHPV86 was placed in the HPV genome homology group A3 by phylogenetic analyses. The overlapping PCR technique is applicable for characterizing the complete spectrum and variation of HPVs in a population.