Evolution of BK Virus Based on Complete Genome Data

The human polyomavirus BK virus (BKV) is ubiquitous in humans, infecting children asymptomatically. BKV is the only primate polyomavirus that has subtypes (I–IV) distinguishable by immunological reactivity. Nucleotide (nt) variations in a major capsid protein (VP1) gene region (designated the epitope region), probably responsible for antigenic diversity, have been used to classify BKV isolates into subtypes. Here, with all the protein-encoding gene sequences, we attempted to elucidate the evolutionary relationships among 28 BKV isolates belonging to subtypes I, III, and IV (no isolate belonging to subtype II, a minor one, was included). First, using the GTR + Γ + I model, maximum likelihood trees were reconstructed for individual viral genes as well as for concatenated viral genes. On the resultant trees, the 28 BKV isolates were consistently divided into three clades corresponding to subtypes I, III, and IV, although bootstrap probabilities are not always high. Then we used more sophisticated likelihood models, one of which takes account of codon structure, to elucidate the phylogenetic relationships among BKV subtypes, but the phylogeny of the deep branchings remained ambiguous. Furthermore, the possibility of positive selection in the evolution of BKV was examined using the nonsynonymous/synonymous rate ratio as a measure of selection. An analysis based on entire genes could not detect any strong evidence for positive selection, but that based on the epitope region identified a few sites potentially under positive selection (these sites were among those showing subtype linked polymorphisms).

[1]  K. Shah,et al.  Characterization of SA12 as a simian virus 40-related papovavirus of chacma baboons , 1977, Infection and immunity.

[2]  A. Mclachlan Tests for comparing related amino-acid sequences. Cytochrome c and cytochrome c 551 . , 1971, Journal of molecular biology.

[3]  J. Booth,et al.  Genomic typing of BK virus in clinical specimens by direct sequencing of polymerase chain reaction products , 1993, Journal of medical virology.

[4]  R. Sharan,et al.  Detecting excess radical replacements in phylogenetic trees. , 2003, Gene.

[5]  H. Kishino,et al.  Time scale of eutherian evolution estimated without assuming a constant rate of molecular evolution. , 2003, Genes & genetic systems.

[6]  Qin Chen,et al.  Subtypes of BK virus prevalent in Japan and variation in their transcriptional control region. , 2004, The Journal of general virology.

[7]  P. Sharp,et al.  Analysis of 15 novel full-length BK virus sequences from three individuals: evidence of a high intra-strain genetic diversity. , 2004, The Journal of general virology.

[8]  Ziheng Yang,et al.  Statistical methods for detecting molecular adaptation , 2000, Trends in Ecology & Evolution.

[9]  A. Limaye,et al.  BK virus‐associated nephropathy in kidney transplant recipients , 2004, Reviews in medical virology.

[10]  T. Kitamura,et al.  JC virus genotyping offers a new paradigm in the study of human populations , 2004, Reviews in medical virology.

[11]  D. Coleman,et al.  New human papovavirus (B.K.) isolated from urine after renal transplantation. , 1971, Lancet.

[12]  J. Heritage,et al.  The persistence of papovavirus BK DNA sequences in normal human renal tissue , 1981, Journal of medical virology.

[13]  J. Felsenstein CONFIDENCE LIMITS ON PHYLOGENIES: AN APPROACH USING THE BOOTSTRAP , 1985, Evolution; international journal of organic evolution.

[14]  L. Jin Rapid genomic typing of BK virus directly from clinical specimens. , 1993, Molecular and cellular probes.

[15]  T. Ohta THE NEARLY NEUTRAL THEORY OF MOLECULAR EVOLUTION , 1992 .

[16]  J. Heritage,et al.  Persistence of DNA sequences of BK virus and JC virus in normal human tissues and in diseased tissues. , 1983, The Journal of infectious diseases.

[17]  Z. Yang,et al.  Models of amino acid substitution and applications to mitochondrial protein evolution. , 1998, Molecular biology and evolution.

[18]  E. Harley,et al.  Genome analysis of BK (WW) viral DNA cloned directly from human urine. , 1984, Intervirology.

[19]  Dirk Husmeier,et al.  TOPALi: software for automatic identification of recombinant sequences within DNA multiple alignments , 2004, Bioinform..

[20]  D. Coleman,et al.  A prospective study of human polyomavirus infection in pregnancy. , 1981, The Journal of infectious diseases.

[21]  G. Stoner,et al.  The Human Polyomaviruses: Past, Present, and Future , 2002 .

[22]  W. Blattner,et al.  BK virus and a new type of JC virus excreted by HIV-1 positive patients in rural Tanzania , 2005, Archives of Virology.

[23]  J. Tavis,et al.  Nucleotide sequence of the human polyomavirus AS virus, an antigenic variant of BK virus , 1989, Journal of virology.

[24]  J. Pipas,et al.  Simian agent 12 is a BK virus-like papovavirus which replicates in monkey cells , 1985, Journal of virology.

[25]  M. Hasegawa,et al.  Evolution of Human Polyomavirus JC: Implications for the Population History of Humans , 2002, Journal of Molecular Evolution.

[26]  Masami Hasegawa,et al.  Mitochondrial phylogenetics and evolution of mysticete whales. , 2005, Systematic biology.

[27]  H. Kishino,et al.  Maximum likelihood inference of protein phylogeny and the origin of chloroplasts , 1990, Journal of Molecular Evolution.

[28]  Ziheng Yang Maximum Likelihood Estimation on Large Phylogenies and Analysis of Adaptive Evolution in Human Influenza Virus A , 2000, Journal of Molecular Evolution.

[29]  K. Khalili,et al.  Human Polyomaviruses: Molecular and Clinical Perspectives , 2002 .

[30]  W. Bruno,et al.  Performance of a divergence time estimation method under a probabilistic model of rate evolution. , 2001, Molecular biology and evolution.

[31]  W. Pauw,et al.  Isolation of a variant of BK virus with altered restriction endonuclease pattern , 2005, Archives of Virology.

[32]  Z. Yang,et al.  Preponderance of slightly deleterious polymorphism in mitochondrial DNA: nonsynonymous/synonymous rate ratio is much higher within species than between species. , 1998, Molecular biology and evolution.

[33]  W. Knowles,et al.  BK virus antigenic variants: Sequence analysis within the capsid VP1 epitope , 1993, Journal of medical virology.

[34]  M. O. Dayhoff,et al.  Atlas of protein sequence and structure , 1965 .

[35]  M. Kimura The Neutral Theory of Molecular Evolution: Introduction , 1983 .

[36]  H. Akaike A new look at the statistical model identification , 1974 .

[37]  J. Felsenstein Evolutionary trees from DNA sequences: A maximum likelihood approach , 2005, Journal of Molecular Evolution.

[38]  H. Kishino,et al.  Evaluation of the maximum likelihood estimate of the evolutionary tree topologies from DNA sequence data, and the branching order in hominoidea , 1989, Journal of Molecular Evolution.

[39]  J. Goudsmit,et al.  Human papovavirus isolated from urine of a child with acute tonsillitis. , 1981, British medical journal.

[40]  P. Howley,et al.  Physical map of the BK virus genome , 1975, Journal of virology.

[41]  C. Cole Polyomaviridae : the viruses and their replication , 2001 .

[42]  W. Wong,et al.  Bayes empirical bayes inference of amino acid sites under positive selection. , 2005, Molecular biology and evolution.

[43]  Takashi Miyata,et al.  Molecular evolution of mRNA: A method for estimating evolutionary rates of synonymous and amino acid substitutions from homologous nucleotide sequences and its application , 1980, Journal of Molecular Evolution.

[44]  Masami Hasegawa,et al.  Accuracies of the simple methods for estimating the bootstrap probability of a maximum-likelihood tree , 1994 .

[45]  R. Grantham Amino Acid Difference Formula to Help Explain Protein Evolution , 1974, Science.

[46]  Joseph P. Bielawski,et al.  Widespread Adaptive Evolution in the Human Immunodeficiency Virus Type 1 Genome , 2003, Journal of Molecular Evolution.

[47]  Y. Yogo,et al.  Growth efficiency of naturally occurring BK virus variants in vivo and in vitro , 1989, Journal of virology.

[48]  R. DeSalle,et al.  Adaptive Evolution of Genes and Genomes , 2000, Heredity.

[49]  C. Ryschkewitsch,et al.  Molecular genotyping of BK and JC viruses in human polyomavirus-associated interstitial nephritis after renal transplantation. , 2001, American journal of kidney diseases : the official journal of the National Kidney Foundation.

[50]  R. Nielsen,et al.  Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene. , 1998, Genetics.

[51]  D. Swofford PAUP*: Phylogenetic analysis using parsimony (*and other methods), Version 4.0b10 , 2002 .

[52]  K. Holsinger The neutral theory of molecular evolution , 2004 .

[53]  H Kishino,et al.  Converting distance to time: application to human evolution. , 1990, Methods in enzymology.

[54]  R. Yang,et al.  BK virus DNA: complete nucleotide sequence of a human tumor virus. , 1979, Science.

[55]  J. Oliver,et al.  The general stochastic model of nucleotide substitution. , 1990, Journal of theoretical biology.

[56]  J. Booth,et al.  Prevalence and distribution of BK virus subtypes in healthy people and immunocompromised patients detected by PCR-restriction enzyme analysis. , 1995, Clinical and diagnostic virology.

[57]  T. Gojobori,et al.  Reevaluation of Amino Acid Variability of the Human Immunodeficiency Virus Type 1 gp120 Envelope Glycoprotein and Prediction of New Discontinuous Epitopes , 2000, Journal of Virology.

[58]  M. O. Dayhoff,et al.  22 A Model of Evolutionary Change in Proteins , 1978 .

[59]  Z. Yang,et al.  Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution. , 1998, Molecular biology and evolution.

[60]  T. Iida,et al.  Typing of human polyomavirus JC virus on the basis of restriction fragment length polymorphisms , 1991, Journal of clinical microbiology.

[61]  N. Goldman,et al.  A codon-based model of nucleotide substitution for protein-coding DNA sequences. , 1994, Molecular biology and evolution.

[62]  S. Miyazawa,et al.  Two types of amino acid substitutions in protein evolution , 1979, Journal of Molecular Evolution.

[63]  G. B. Orsi,et al.  Detection of BK polyomavirus genotypes in healthy and HIV-positive children , 1997, European Journal of Epidemiology.

[64]  Ziheng Yang,et al.  PAML: a program package for phylogenetic analysis by maximum likelihood , 1997, Comput. Appl. Biosci..

[65]  David Posada,et al.  MODELTEST: testing the model of DNA substitution , 1998, Bioinform..

[66]  Z. Yang,et al.  Among-site rate variation and its impact on phylogenetic analyses. , 1996, Trends in ecology & evolution.

[67]  W. Knowles The Epidemiology of BK Virus and the Occurrence of Antigenic and Genomic Subtypes , 2002 .

[68]  H. Kishino,et al.  Estimating the rate of evolution of the rate of molecular evolution. , 1998, Molecular biology and evolution.

[69]  G. Khoury,et al.  The genome of human papovavirus BKV , 1979, Cell.

[70]  J. Pipas,et al.  Complete Nucleotide Sequence of Polyomavirus SA12 , 2005, Journal of Virology.