Bioinformatics Advance Access originally published online on January 12, 2005
Bioinformatics 2005 21(9):2104-2105; doi:10.1093/bioinformatics/bti263
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ProtTest: selection of best-fit models of protein evolution
1Department of Biochemistry, Genetics and Immunology, Universidad de Vigo 36310 Vigo, Spain
2Department of Biodiversity and Evolutionary Biology, Museo Nacional de Ciencias Naturales 28006 Madrid, Spain
*To whom correspondence should be addressed.
| Abstract |
|---|
|
|
|---|
Summary: Using an appropriate model of amino acid replacement is very important for the study of protein evolution and phylogenetic inference. We have built a tool for the selection of the best-fit model of evolution, among a set of candidate models, for a given protein sequence alignment.
Availability: ProtTest is available under the GNU license from http://darwin.uvigo.es
Contact: fabascal{at}uvigo.es
| 1 INTRODUCTION |
|---|
|
|
|---|
1.1 Models of protein evolution
Models of protein evolution, or amino acid replacement, describe the probabilities of change from one amino acid to another, and therefore become indispensable tools for the characterization of the process of protein evolution (Thorne, 2000; Thorne and Goldman, 2003). Indeed, these models provide the foundation for the reconstruction of protein phylogenies under distance, maximum likelihood and Bayesian methods. Dayhoff et al. (1978) introduced the most influential model of amino acid replacement, a 20-state time reversible homogeneous Markov model. Because the large number of parameters in a 20-state replacement matrix, estimates of these parameters are usually obtained from large datasets prior to the analysis of the dataset of interest. In this way, different empirical matrices with fixed relative rates of amino acid replacement have been already proposed, like the Dayhoff matrix (Dayhoff et al., 1978), the JTT matrix (Jones et al., 1992), the mtREV matrix (Adachi and Hasegawa, 1996) or the WAG matrix (Whelan and Goldman, 2001). While these models generally assume that the process of amino acid replacement is very similar across all positions, conservation of protein function and structure imposes constraints on which positions can change. This evolutionary information can be inferred by considering a fraction of amino acids to be invariable (+I) (Reeves, 1992), or assigning each site a probability to belong to given rate categories (+G) (Yang, 1993). Additionally, observed amino acid frequencies can also be considered (+F) (Cao et al., 1994).
1.2 Model selection and inference
Model selection may be seen as a way of identifying the model that, among a set of candidates, is closest to reality. Looking for a balance between accuracy and simplicity, Akaike (1973) found a simple relationship between the likelihood (L) and the number of parameters (K):
![]() |
![]() |
![]() |
| 2 THE PROGRAM: PROTTEST |
|---|
|
|
|---|
Although widely-used software exists for the selection of the best-fit nucleotide models (Posada and Crandall, 1998), no program has been developed until now for protein models. ProtTest is a java program to find the best model of amino acid replacement for a given protein alignment. It is based on the Phyml program (Guindon and Gascuel, 2003) for the ML optimizations, modified to support +F and four extra substitution matrices and uses the PAL library (Drummond and Strimmer, 2001) for handling protein alignments and trees. ProtTest is available for Mac OSX, Linux and Windows, and it can be run in three ways: using a GUI, at the command-line and through the web. Its basic workflow is summarized in Figure 1.
|
Given a protein alignment and a tree topology the program calculates the likelihood under each candidate model, and estimates model parameters. The current version 1.2 implements 64 empirical models: the eight matrices WAG, mtREV, Dayhoff, JTT, VT, Blosum62, CpREV and RtREV under +F, +G, +I and their combinations. Other models exist, particularly mechanical models, that are not implemented in ProtTest. For each model, the tree topology can be fixed [provided by the user or calculated by BIONJ (Gascuel, 1997)] or optimized under ML. After this, the user can choose a model selection strategy (AIC, AICc, BIC), and obtain a rank of model fits, model-averaged parameter estimates or measures of parameter importance. For the AICc and the BIC, sample size is set by default to the number of positions in the alignment. Other options to define sample size attempt to take into account both the number of sequences and their redundancy. Other valuable features include the ability to restrict the set of candidate models (only in the GUI version) and the possibility to output the tree corresponding to the best model.
| Acknowledgments |
|---|
Special thanks to Stephane Guindon (Phyml) and Matthew Goode (PAL) for their help. This work was financially supported from a research grant in bioinformatics from the Fundación BBVA (Spain).
Received on October 7, 2004; revised on December 15, 2004; accepted on January 5, 2005
| REFERENCES |
|---|
|
|
|---|
Adachi, J. and Hasegawa, M. (1996) Model of amino acid substitution in proteins encoded by mitochondrial DNA. J. Mol. Evol., 42, 459468[ISI][Medline].
Akaike, H. (1973) Information theory and an extension of the maximum likelihood principle. Proceedings of 2nd International Symposium on Information Theory, , Budapest, Hungary , pp. 267281.
Burnham, K.P. and Anderson, D.R. Model Selection and Multi-Model Inference. A Practical InformationTheoretic Approach, (2002) 2nd edn , New York Springer-Verlag.
Cao, Y., Adachi, J., Janke, A., Paabo, S., Hasegawa, M. (1994) Phylogenetic relationships among eutherian orders estimated from inferred sequences of mitochondrial proteins: instability of a tree based on a single gene. J. Mol. Evol., 39, 519527[ISI][Medline].
Dayhoff, M.O., Schwartz, R.M., Orcutt, B.C. (1978) A model of evolutionary change in proteins. In Dayhoff, M.O. (Ed.). Atlas of Protein Sequence and Structure, , Washington, DC National Biomedical Research Foundation, pp. 345352.
Drummond, A. and Strimmer, K. (2001) PAL: an object-oriented programming library for molecular evolution and phylogenetics. Bioinformatics, 17, 662663
Gascuel, O. (1997) BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol. Biol. Evol., 14, 685695[Abstract].
Guindon, S. and Gascuel, O. (2003) A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst. Biol., 52, 696704[CrossRef][ISI][Medline].
Jones, D.T., Taylor, W.R., Thornton, J.M. (1992) The rapid generation of mutation data matrices from protein sequences. Comp. Appl. Biosci., 8, 275282.
Posada, D. and Buckley, T.R. (2004) Model selection and model averaging in phylogenetics: advantages of AIC and Bayesian approaches over likelihood ratio tests. Syst. Biol., 53, 793808[CrossRef][ISI][Medline].
Posada, D. and Crandall, K.A. (1998) MODELTEST: testing the model of DNA substitution. Bioinformatics, 14, 817818
Reeves, J.H. (1992) Heterogeneity in the substitution process of amino acid sites of proteins coded for by mitochondrial DNA. J. Mol. Evol., 35, 1731[CrossRef][ISI][Medline].
Schwarz, G. (1978) Estimating the dimension of a model. Ann. Statist., 6, 461464.
Sugiura, N. (1978) Further analysis of the data by Akaikes information criterion and the finite correction. Comm. Statist. A-Theory. Meth., 7, 1326.
Thorne, J.L. (2000) Models of protein sequence evolution and their applications. Curr. Opin. Genet. Dev., 10, 602605[CrossRef][ISI][Medline].
Thorne, J.L. and Goldman, N. (2003) Probabilistic models for the study of protein evolution. In Balding, D.J., Bishop, M., Cannings, C. (Eds.). Handbook of Statistical Genetics, , Chichester, England John Wiley & Sons, Ltd, pp. 209226.
Whelan, S. and Goldman, N. (2001) A general empirical model of protein evolution derived from multiple protein families using a maximum-likelihood approach. Mol. Biol. Evol., 18, 691699
Yang, Z. (1993) Maximum-likelihood estimation of phylogeny from DNA sequences when substitution rates differ over sites. Mol. Biol. Evol., 10, 13961401[Abstract].
This article has been cited by other articles:
![]() |
A. Moya, S. Tambutte, A. Bertucci, E. Tambutte, S. Lotto, D. Vullo, C. T. Supuran, D. Allemand, and D. Zoccola Carbonic Anhydrase in the Scleractinian Coral Stylophora pistillata: CHARACTERIZATION, LOCALIZATION, AND ROLE IN BIOMINERALIZATION J. Biol. Chem., September 12, 2008; 283(37): 25475 - 25484. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Simon, N. A.C. Clarke, B. A. McNeil, I. Johnson, D. Pantuso, L. Dai, D. Chai, and S. Zimmerly Group II introns in Eubacteria and Archaea: ORF-less introns and new varieties RNA, September 1, 2008; 14(9): 1704 - 1713. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. D'Aniello, M. Irimia, I. Maeso, J. Pascual-Anaya, S. Jimenez-Delgado, S. Bertrand, and J. Garcia-Fernandez Gene Expansion and Retention Leads to a Diverse Tyrosine Kinase Superfamily in Amphioxus Mol. Biol. Evol., September 1, 2008; 25(9): 1841 - 1854. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Rokas and S. B. Carroll Frequent and Widespread Parallel Evolution of Protein Sequences Mol. Biol. Evol., September 1, 2008; 25(9): 1943 - 1953. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. M. Codoner and S. F. Elena The promiscuous evolutionary history of the family Bromoviridae J. Gen. Virol., July 1, 2008; 89(7): 1739 - 1747. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Dereeper, V. Guignon, G. Blanc, S. Audic, S. Buffet, F. Chevenet, J.-F. Dufayard, S. Guindon, V. Lefort, M. Lescot, et al. Phylogeny.fr: robust phylogenetic analysis for the non-specialist Nucleic Acids Res., July 1, 2008; 36(suppl_2): W465 - W469. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. Kirkham, S. J. Nixon, M. T. Howes, L. Abi-Rached, D. E. Wakeham, M. Hanzal-Bayer, C. Ferguson, M. M. Hill, M. Fernandez-Rojo, D. A. Brown, et al. Evolutionary analysis and molecular dissection of caveola biogenesis J. Cell Sci., June 15, 2008; 121(12): 2075 - 2086. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Hernandez-Sanchez, A. Mansilla, F. de Pablo, and R. Zardoya Evolution of the Insulin Receptor Family and Receptor Isoform Expression in Vertebrates Mol. Biol. Evol., June 1, 2008; 25(6): 1043 - 1053. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Larroux, G. N. Luke, P. Koopman, D. S. Rokhsar, S. M. Shimeld, and B. M. Degnan Genesis and Expansion of Metazoan Transcription Factor Gene Classes Mol. Biol. Evol., May 1, 2008; 25(5): 980 - 996. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Pantalacci, A. Chaumot, G. Benoit, A. Sadier, F. Delsuc, E. J. P. Douzery, and V. Laudet Conserved Features and Evolutionary Shifts of the EDA Signaling Pathway Involved in Vertebrate Skin Appendage Development Mol. Biol. Evol., May 1, 2008; 25(5): 912 - 928. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. H. Struck and F. Fisse Phylogenetic Position of Nemertea Derived from Phylogenomic Data Mol. Biol. Evol., April 1, 2008; 25(4): 728 - 736. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. A. Cymerman, I. Chung, B. M. Beckmann, J. M. Bujnicki, and G. Meiss EXOG, a novel paralog of Endonuclease G in higher eukaryotes Nucleic Acids Res., March 27, 2008; 36(4): 1369 - 1379. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Zhang, D. Bhattacharya, L. Maranda, and S. Lin Mitochondrial cob and cox1 Genes and Editing of the Corresponding mRNAs in Dinophysis acuminata from Narragansett Bay, with Special Reference to the Phylogenetic Position of the Genus Dinophysis Appl. Envir. Microbiol., March 1, 2008; 74(5): 1546 - 1554. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. B. Dacks, P. P. Poon, and M. C. Field Phylogeny of endocytic components yields insight into the process of nonendosymbiotic organelle evolution PNAS, January 15, 2008; 105(2): 588 - 593. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Podsiadlowski, A. Braband, and G. Mayer The Complete Mitochondrial Genome of the Onychophoran Epiperipatus biolleyi Reveals a Unique Transfer RNA Set and Provides Further Support for the Ecdysozoa Hypothesis Mol. Biol. Evol., January 1, 2008; 25(1): 42 - 51. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Casola, D. Hucks, and C. Feschotte Convergent Domestication of pogo-like Transposases into Centromere-Binding Proteins in Fission Yeast and Mammals Mol. Biol. Evol., January 1, 2008; 25(1): 29 - 41. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Cherif-Zahar, A. Durand, I. Schmidt, N. Hamdaoui, I. Matic, M. Merrick, and G. Matassi Evolution and Functional Characterization of the RH50 Gene from the Ammonia-Oxidizing Bacterium Nitrosomonas europaea J. Bacteriol., December 15, 2007; 189(24): 9090 - 9100. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. R. Tabita, T. E. Hanson, H. Li, S. Satagopan, J. Singh, and S. Chan Function, Structure, and Evolution of the RubisCO-Like Proteins and Their RubisCO Homologs Microbiol. Mol. Biol. Rev., December 1, 2007; 71(4): 576 - 599. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. R. Bench, T. E. Hanson, K. E. Williamson, D. Ghosh, M. Radosovich, K. Wang, and K. E. Wommack Metagenomic Characterization of Chesapeake Bay Virioplankton Appl. Envir. Microbiol., December 1, 2007; 73(23): 7629 - 7641. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Q. Matus, K. M. Halanych, and M. Q. Martindale The Hox gene complement of a pelagic chaetognath, Flaccisagitta enflata Integr. Comp. Biol., December 1, 2007; 47(6): 854 - 864. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. V. Goldstone, H. M. H. Goldstone, A. M. Morrison, A. Tarrant, S. E. Kern, B. R. Woodin, and J. J. Stegeman Cytochrome P450 1 Genes in Early Deuterostomes (Tunicates and Sea Urchins) and Vertebrates (Chicken and Frog): Origin and Diversification of the CYP1 Gene Family Mol. Biol. Evol., December 1, 2007; 24(12): 2619 - 2631. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Hausdorf, M. Helmkampf, A. Meyer, A. Witek, H. Herlyn, I. Bruchhaus, T. Hankeln, T. H. Struck, and B. Lieb Spiralian Phylogenomics Supports the Resurrection of Bryozoa Comprising Ectoprocta and Entoprocta Mol. Biol. Evol., December 1, 2007; 24(12): 2723 - 2729. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. B. Rounge, T. Rohrlack, A. Tooming-Klunderud, T. Kristensen, and K. S. Jakobsen Comparison of Cyanopeptolin Genes in Planktothrix, Microcystis, and Anabaena Strains: Evidence for Independent Evolution within Each Genus Appl. Envir. Microbiol., November 15, 2007; 73(22): 7322 - 7330. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Blanc, H. Ogata, C. Robert, S. Audic, J.-M. Claverie, and D. Raoult Lateral gene transfer between obligate intracellular bacteria: Evidence from the Rickettsia massiliae genome Genome Res., November 1, 2007; 17(11): 1657 - 1664. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Robinson and L. A. Katz Non-Mendelian Inheritance of Paralogs of 2 Cytoskeletal Genes in the Ciliate Chilodonella uncinata Mol. Biol. Evol., November 1, 2007; 24(11): 2495 - 2503. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. Ahlenstiel, K. Roomp, M. Daumer, J. Nattermann, M. Vogel, J. K. Rockstroh, N. Beerenwinkel, R. Kaiser, H.-D. Nischalke, T. Sauerbruch, et al. Selective Pressures of HLA Genotypes and Antiviral Therapy on Human Immunodeficiency Virus Type 1 Sequence Mutation at a Population Level Clin. Vaccine Immunol., October 1, 2007; 14(10): 1266 - 1273. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Yasui, E.-Y. Kim, H. Iwata, D. G. Franks, S. I. Karchner, M. E. Hahn, and S. Tanabe Functional Characterization and Evolutionary History of Two Aryl Hydrocarbon Receptor Isoforms (AhR1 and AhR2) from Avian Species Toxicol. Sci., September 1, 2007; 99(1): 101 - 117. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. J Agulleiro, M. Andre, S. Morais, J. Cerda, and P. J Babin High Transcript Level of Fatty Acid-Binding Protein 11 but Not of Very Low-Density Lipoprotein Receptor Is Correlated to Ovarian Follicle Atresia in a Teleost Fish (Solea senegalensis) Biol Reprod, September 1, 2007; 77(3): 504 - 516. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. G. Hoffmann and J. F. Storz The {alpha}D-Globin Gene Originated via Duplication of an Embryonic {alpha}-Like Globin Gene in the Ancestor of Tetrapod Vertebrates Mol. Biol. Evol., September 1, 2007; 24(9): 1982 - 1990. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. M. Hallstrom, M. Kullberg, M. A. Nilsson, and A. Janke Phylogenomic Data Analyses Provide Evidence that Xenarthra and Afrotheria Are Sister Groups Mol. Biol. Evol., September 1, 2007; 24(9): 2059 - 2068. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Bleidorn, I. Eeckhaut, L. Podsiadlowski, N. Schult, D. McHugh, K. M. Halanych, M. C. Milinkovitch, and R. Tiedemann Mitochondrial Genome and Nuclear Sequence Data Support Myzostomida As Part of the Annelid Radiation Mol. Biol. Evol., August 1, 2007; 24(8): 1690 - 1701. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. D. Hackett, H. S. Yoon, S. Li, A. Reyes-Prieto, S. E. Rummele, and D. Bhattacharya Phylogenomic Analysis Supports the Monophyly of Cryptophytes and Haptophytes and the Association of Rhizaria with Chromalveolates Mol. Biol. Evol., August 1, 2007; 24(8): 1702 - 1713. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Tarraga, I. Medina, L. Arbiza, J. Huerta-Cepas, T. Gabaldon, J. Dopazo, and H. Dopazo Phylemon: a suite of web tools for molecular evolution, phylogenetics and phylogenomics Nucleic Acids Res., July 13, 2007; 35(suppl_2): W38 - W42. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Slominski, J. Calkiewicz, P. Golec, G. Wegrzyn, and B. Wrobel Plasmids derived from Gifsy-1/Gifsy-2, lambdoid prophages contributing to the virulence of Salmonella enterica serovar Typhimurium: implications for the evolution of replication initiation proteins of lambdoid phages and enterobacteria Microbiology, June 1, 2007; 153(6): 1884 - 1896. [Abstract] [Full Text] [PDF] |
||||
![]() |
V. Ruano-Rubio and M. A. Fares Testing the Neutral Fixation of Hetero-Oligomerism in the Archaeal Chaperonin CCT Mol. Biol. Evol., June 1, 2007; 24(6): 1384 - 1396. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. A. Gladyshev and I. R. Arkhipova From the Cover: Telomere-associated endonuclease-deficient Penelope-like retroelements in diverse eukaryotes PNAS, May 29, 2007; 104(22): 9352 - 9357. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Poggio, C. Abreu-Goodger, S. Fabela, A. Osorio, G. Dreyfus, P. Vinuesa, and L. Camarena A Complete Set of Flagellar Genes Acquired by Horizontal Transfer Coexists with the Endogenous Flagellar System in Rhodobacter sphaeroides J. Bacteriol., April 15, 2007; 189(8): 3208 - 3216. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. C. Doxey, M. W. F. Yaish, B. A. Moffatt, M. Griffith, and B. J. McConkey Functional Divergence in the Arabidopsis {beta}-1,3-Glucanase Gene Family Inferred by Phylogenetic Reconstruction of Expression States Mol. Biol. Evol., April 1, 2007; 24(4): 1045 - 1055. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. G. Koziol, T. Borza, K.-I. Ishida, P. Keeling, R. W. Lee, and D. G. Durnford Tracing the Evolution of the Light-Harvesting Antennae in Chlorophyll a/b-Containing Organisms Plant Physiology, April 1, 2007; 143(4): 1802 - 1816. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. R. de la Vega, R. G. Sevilla, A. Hermoso, J. Lorenzo, S. Tanco, A. Diez, L. D. Fricker, J. M. Bautista, and F. X. Aviles Nna1-like proteins are active metallocarboxypeptidases of a new and diverse M14 subfamily FASEB J, March 1, 2007; 21(3): 851 - 865. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Abascal, D. Posada, and R. Zardoya MtArt: A New Model of Amino Acid Replacement for Arthropoda Mol. Biol. Evol., January 1, 2007; 24(1): 1 - 5. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. L. Porter, T. W. Cronin, D. A. McClellan, and K. A. Crandall Molecular Characterization of Crustacean Visual Pigments and the Evolution of Pancrustacean Opsins Mol. Biol. Evol., January 1, 2007; 24(1): 253 - 268. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Burki and J. Pawlowski Monophyly of Rhizaria and Multigene Phylogeny of Unicellular Bikonts Mol. Biol. Evol., October 1, 2006; 23(10): 1922 - 1930. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. Lieb, K. Dimitrova, H.-S. Kang, S. Braun, W. Gebauer, A. Martin, B. Hanelt, S. A. Saenz, C. M. Adema, and J. Markl Red blood with blue-blood ancestry: Intriguing structure of a snail hemoglobin PNAS, August 8, 2006; 103(32): 12011 - 12016. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Shalchian-Tabrizi, M. Skanseng, F. Ronquist, D. Klaveness, T. R. Bachvaroff, C. F. Delwiche, A. Botnen, T. Tengs, and K. S. Jakobsen Heterotachy Processes in Rhodophyte-Derived Secondhand Plastid Genes: Implications for Addressing the Origin and Evolution of Dinoflagellate Plastids Mol. Biol. Evol., August 1, 2006; 23(8): 1504 - 1515. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. Posada ModelTest Server: a web-based tool for the statistical selection of models of nucleotide substitution online. Nucleic Acids Res., July 1, 2006; 34(Web Server issue): W700 - W703. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. L. Afonso, E. R. Tulman, G. Delhon, Z. Lu, G. J. Viljoen, D. B. Wallace, G. F. Kutish, and D. L. Rock Genome of crocodilepox virus. J. Virol., |






















