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

Zigui Chen, Masanori Terai, Leiping Fu, Rolando Herrero, Rob DeSalle, Robert D. Burk

Research output: Contribution to journalArticle

119 Citations (Scopus)

Abstract

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

Original languageEnglish (US)
Pages (from-to)7014-7023
Number of pages10
JournalJournal of Virology
Volume79
Issue number11
DOIs
StatePublished - Jun 2005

Fingerprint

Human papillomavirus 16
Papillomaviridae
Open Reading Frames
open reading frames
Genome
genome
Amino Acids
amino acids
codons
Codon
uterine cervical neoplasms
Molecular Evolution
Viral Genome
etiological agents
Base Pairing
Uterine Cervical Neoplasms
niches
Nucleotides
nucleotides
neoplasms

ASJC Scopus subject areas

  • Immunology

Cite this

Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses. / Chen, Zigui; Terai, Masanori; Fu, Leiping; Herrero, Rolando; DeSalle, Rob; Burk, Robert D.

In: Journal of Virology, Vol. 79, No. 11, 06.2005, p. 7014-7023.

Research output: Contribution to journalArticle

Chen, Zigui ; Terai, Masanori ; Fu, Leiping ; Herrero, Rolando ; DeSalle, Rob ; Burk, Robert D. / Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses. In: Journal of Virology. 2005 ; Vol. 79, No. 11. pp. 7014-7023.
@article{698d454f28a4429ca6d4dbc73539da9e,
title = "Diversifying selection in human papillomavirus type 16 lineages based on complete genome analyses",
abstract = "Human papillomavirus type 16 (HPV16) is the primary etiological agent of cervical cancer, the second most common cancer in women worldwide. Complete genomes of 12 isolates representing the major lineages of HPV16 were cloned and sequenced from cervicovaginal cells. The sequence variations within the open reading frames (ORFs) and noncoding regions were identified and compared with the HPV16R reference sequence (50). This whole-genome approach gives us unprecedented precision in detailing sequence-level changes that are under selection on a whole-viral-genome scale. Of 7,908 base pair nucleotide positions, 313 (4.0{\%}) were variable. Within the 2,452 amino acids (aa) comprising 8 ORFs, 243 (9.9{\%}) amino acid positions were variable. In order to investigate the molecular evolution of HPV16 variants, maximum likelihood models of codon substitution were used to identify lineages and amino acid sites under selective pressure. Five codon sites in the E5 (aa 48, 65) and E6 (aa 10, 14, 83) ORFs were demonstrated to be under diversifying selective pressure. The ES ORF had the overall highest nonsynonymous/synonymous substitution rate (ω) ratio (M3 = 0.7965). The E2 gene had the next-highest ω ratio (M3 = 0.5611); however, no specific codons were under positive selection. These data indicate that the E6 and E5 ORFs are evolving under positive Darwinian selection and have done so in a relatively short time period. Whether response to selective pressure upon the E5 and E6 ORFs contributes to the biological success of HPV16, its specific biological niche, and/or its oncogenic potential remains to be established.",
author = "Zigui Chen and Masanori Terai and Leiping Fu and Rolando Herrero and Rob DeSalle and Burk, {Robert D.}",
year = "2005",
month = "6",
doi = "10.1128/JVI.79.11.7014-7023.2005",
language = "English (US)",
volume = "79",
pages = "7014--7023",
journal = "Journal of Virology",
issn = "0022-538X",
publisher = "American Society for Microbiology",
number = "11",

}

TY - JOUR

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

AU - Chen, Zigui

AU - Terai, Masanori

AU - Fu, Leiping

AU - Herrero, Rolando

AU - DeSalle, Rob

AU - Burk, Robert D.

PY - 2005/6

Y1 - 2005/6

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

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

UR - http://www.scopus.com/inward/record.url?scp=18744391386&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=18744391386&partnerID=8YFLogxK

U2 - 10.1128/JVI.79.11.7014-7023.2005

DO - 10.1128/JVI.79.11.7014-7023.2005

M3 - Article

VL - 79

SP - 7014

EP - 7023

JO - Journal of Virology

JF - Journal of Virology

SN - 0022-538X

IS - 11

ER -