- Short report
- Open Access
Isolation of a novel monkey adenovirus reveals a new phylogenetic clade in the evolutionary history of simian adenoviruses
Virology Journal volume 8, Article number: 125 (2011)
Adenoviruses of primates include human (HAdV) and simian (SAdV) isolates classified into 8 species (Human Adenovirus A to G, and Simian Adenovirus A). In this study, a novel adenovirus was isolated from a colony of cynomolgus macaques (Macaca fascicularis) and subcultured in VERO cells. Its complete genome was purified and a region encompassing the hexon gene, the protease gene, the DNA binding protein (DBP) and the 100 kDa protein was amplified by PCR and sequenced by primer walking. Sequence analysis of these four genes showed that the new isolate had 80% identity to other primate adenoviruses and lacked recombination events. The study of the evolutionary relationships of this new monkey AdV based on the combined sequences of the four genes supported a close relationship to SAdV-3 and SAdV-6, lineages isolated from Rhesus monkeys. The clade formed by these three types is separated from the remaining clades and establishes a novel branch that is related to species HAdV-A, F and G. However, the genetic distance corresponding to the newly isolated monkey AdV considerably differs from these as to belong to a new, not yet established species. Results presented here widen our knowledge on SAdV and represents an important contribution to the understanding of the evolutionary history of primate adenoviruses.
Adenoviruses are non-enveloped, icosahedral, double-stranded DNA viruses known to cause gastroenteritis, keratoconjunctivitis, and acute respiratory disorders . Adenoviruses infecting primates belong to the genus Mastadenovirus. To date, at least 52 distinct human adenovirus (HAdV) types have been described. Phylogenetically, HAdVs cluster into 6 species (named Human Adenovirus A to F, HAdV-A to F), although recently a new species (HAdV-G) has been described to include the types HAdV-52 and SAdV-1 . Moreover, at least 25 simian adenovirus (SAdV) types have been recognized: 20 types (SAdV-1 to 20) were isolated from Old World monkeys and are more related to species HAdV-A and F [3, 4], whereas 5 types (SAdV-21 to 25) were isolated from chimpanzees and are closely related to HAdV-4 (HAdV-E) and to HAdV-B [5–7]. Recently, new isolates from chimpanzee, bonobo, orangutan, gorillas, and macaques have been described and characterized, largely expanding the SAdV taxonomy [8, 9].
HAdVs have traditionally been classified according to their immunochemical and biological properties. Nowadays, however, microbiologists favour to classify HAdVs based on their evolutionary relationships as inferred in phylogenetic analyses of DNA sequences of viral protein genes. This approach allows for a classification based on the evolutionary history of the viruses, besides avoiding misleading classification due to cross-reaction in neutralization and hemagglutination tests. In addition to the γ determinant of the fiber protein, the main type-specific epitope in adenovirus is the ε determinant present on the hexon capsid protein. DNA amplification and sequencing of these informative regions in the adenoviral genome have been used to provide molecular data for typing strains and identification of new prototypic isolates [10, 11]. Moreover, the hexon and the protease genes have been characterized in most adenovirus types providing a large amount of data for taxonomic classification of new isolates.
In this study we describe the isolation and identification of an adenoviral type strain from cages inhabited by a colony of cynomolgus monkeys (Macaca fascicularis). Viral concentrates were obtained from a mix of excreta and waste collected from bedding areas of a colony of M. fascicularis showing no signs of gastroenteritis or pathogenic disorders. Viruses were eluted and concentrated from the samples as described previously . African green monkey kidney epithelial cells (Vero) were infected blind with the viral concentrates and cytopathic effect was rapidly observed. Supernatant from infected Vero cells was prepared for analysis by transmission electron microscopy using copper grids coated with formvar film and carbon that were negatively stained with 2% phosphotungstic acid, pH 7.0. Analysis of the grids using a JEOL JEM-1010 identified adenoviral particles characterized by typical 70 nm icosahedral-shaped capsids (Figure 1). The virus (hereafter referred to as CynAdV) was purified by infection of end-point limiting dilutions in Vero cells for 3 passages. A purified stock was generated and the virus titer was determined as 107.5 TCID50 ml-1 by calculation using the Reed-Muench cumulative method.
The entire viral DNA genome was purified as described by Kovács et al. . A random library of DNA fragments was prepared and some clones were sequenced. To address the identification of the new isolate, a genome fragment encompassing traditionally studied viral genes like the hexon and the protease gene was targeted for amplification and sequencing. Two fragments of approximately 3,500 bp covering the whole coding sequence for the hexon, protease, DBP genes, and part of the 100 kDa gene were then amplified with FideliTaq (Amersham Biosciences) following the manufacturer's instructions. The fragments were purified with the QIAquick PCR purification kit (QIAgen) and sequenced with the ABI PRISM Dye Terminator Cycle Sequencing Ready Reaction kit (Applied Biosystems) using primer walking. The sequence was deposited in the GenBank with accession number EU293065.
The CynAdV genome contained a 2,775 bp-long hexon gene, a 615 bp-long protease gene, and a 1,320 bp-long DBP gene. The overall base composition of the 7,481 bp sequence was 21.43% A, 21.61% T, 26.55% G, and 30.41% C. The GC content was 56.96%, similar to those reported for the complete genome of HAdV-52 (55.1%) and SAdV-1 (55%). Interestingly, no intergenic region was observed between the hexon and protease genes, suggesting that the CynAdV genome size might be smaller than other SAdVs. When the predicted amino acid sequences were used as queries in BLASTp searches, CynAdV showed an overall amino acid identity of 75-88% for the hexon protein, a structural protein with long conserved domains (Table 1). Amino acid identity was reduced in the non-structural proteins, protease and DBP. The CynAdV protease consisted of a protein of 205 residues and showed 75-80% amino acid identity to most human and simian adenovirus protease sequences, whereas the DBP consisted of 493 residues and showed 46-59% amino acid identity to most adenoviral DBP sequences (Table 1).
Recombination among adenovirus serotypes generates new strains (intermediate types) that are difficult to classify. Due to selective pressure of neutralizing antibodies and viral escape, recombination events involving the hexon are common, as has recently been described [14, 15]. The hexon protein in HAdVs can be divided into 4 conserved (C1-4) and 3 variable regions (V1-3) . It contains the type-specific neutralization epitope (ε determinant) which consists of 6 hypervariable regions located in L1 and an extra one located in L2 . High identity of these variable regions with other HAdVs lineages could indicate recombination in this particular gene. To assess whether recombination occurred in the newly isolated monkey AdV hexon gene, we performed bootscan analysis of the prototypic hexon sequences available using SimPlot software . Hexon gene sequences were aligned using ClustalW and the CynAdV sequence was analyzed against those of HAdV and SAdV types grouped according to species. Grouping reduces effects of variation in individual sequences. However, considering the similarity between CynAdV and SAdV-7 and members of the species HAdV-G, these sequences were left ungrouped. The results of the bootscan analysis failed to identify recombination events (Figure 2).
In order to classify the new isolate, the evolutionary relationships of CynAdV were assessed in a phylogenetic analysis including all publicly available genomes of HAdVs and SAdVs. Sequences were aligned with MAFFT  using standard gap parameters and 100 refinement iterations. The alignment was then trimmed to match the coding regions of the genes sequenced for CynAdV: hexon, protease, DBP, and 100 kDa genes. Tree searches were performed on each gene individually and on the combined dataset (coding sequences of all 4 genes concatenated) using maximum likelihood methods as implemented in the program RAxML 7.0.4 . For all matrices, the evolutionary model used was the GTRGAMMA and 20 search replicates were performed. In the analysis of the combined dataset, model parameters were obtained for each gene separately (partition wise). Statistical significance of clades was obtained with 100 bootstrap searches. A second analysis of the combined dataset was done after removing highly variable regions including many gaps that could cause low accuracy in the alignment.
The resulting trees (combined dataset in Figure 3, individual gene trees in Figure 4) recovered all adenovirus species groups as in Roy et al. . The combined dataset tree had high support for all previously described HAdVs species and for most relationships among groups. In this tree, CynAdV clustered the closest with a clade including species SAdV-3 and SAdV-6 with moderately high bootstrap support (88%). The support for this relationship increased (94%) when highly variable regions were removed (Figure 3). Individual gene trees (Figure 4) disagreed on the relationships among species HAdV-A, F, G, the clade formed by types SAdV-3 and SAdV-6, and CynAdV. The conflicting relationships, however, did not receive high bootstrap support in any of the gene trees. The only exception was the DBP gene tree that agreed and showed high statistical support for the relationships of CynAdV depicted in the combined dataset tree. To check whether high levels of variation in DBP (as illustrated by low sequence identity across species, Table 1) could be biasing the combined analysis , we reanalyzed this gene using this time its amino acid sequence. The resulting tree is very similar to the one based on nucleotides, suggesting that this gene is not highly affected by sequence bias and that it is probably a good marker to use in the classification of adenoviruses (Figure 5).
Phylogenetic signal conflict among genes was assessed with the incongruence-length difference (ILD) test  as implemented in the software PAUP*. The IDL test showed significant conflict among the phylogenetic signal of individual genes (p = 0.001). One explanation for this conflict is nucleotide substitution saturation, which may blur phylogenetic signal in individual genes . This explanation is supported by the increase in bootstrap percentages when highly variable regions were removed from the combined dataset. An alternative explanation is recombination between an ancestor of the CynAdV virus and other viruses of related species groups A, F, or G. The low bootstrap percentages obtained for CynAdV in individual gene analyses (except for DBP), however, precludes the formulation of a specific hypothesis for the recombination event. As shown before, the test for recombination involving the hexon gene was negative.
In a previous study , three different lineages of adenovirus isolated from cynomolgus macaques were identified. Interestingly, these three lineages were also associated to SAdV-3 and SAdV-6 (isolated from Rhesus monkeys) in trees based on the DBP gene. Nevertheless, these relationships were not observed for two of the newly described varieties (SAdV-49, and SAdV-50) in trees obtained with other genes. In fact, their results were very similar to ours in that the relationships of the lineages classified as species A, F, and G were largely discordant. Unfortunately, since the sequences of the other cynomolgus macaque's adenoviruses - SAdV-48, SAdV-49, and SAdV-50 - were not publicly available, they could not be included in the phylogenetic analyses presented here.
In the past years, sequencing of the AdV hexon gene has proved to be sufficient for typing new isolates [11, 16]. However, phylogenetic data inferred from hexon gene analysis needs to be contrasted against other viral genes since selective pressure on the hexon (and fiber genes) may lead to recombination events and convergent evolution. In this study, analysis of a genomic region encompassing 4 different genes including the hexon gene distinguished the newly isolated monkey AdV as an independent clade close to the types SAdV-6 and SAdV-3. The trees shown here suggest that SAdV-6 may belong to species SAdV-A together with SAdV-3. However, CynAdV seems to adequately differ from these to belong to a separate (not yet established) species. This study represents an important contribution to the understanding of the diversity and evolutionary history of primate adenoviruses.
Swenson PDWG, Allard A, Hierholzer JC: Adenoviruses. Manual of Clinical Microbiology. Edited by: Murray PR BE, Jorgensen JH, Pfaller MA. 2003, Yolken RH: ASM Press, 2: 1404-1417.
Jones MS, Harrach B, Ganac RD, Gozum MM, Dela Cruz WP, Riedel B, Pan C, Delwart EL, Schnurr DP: New adenovirus species found in a patient presenting with gastroenteritis. J Virol. 2007, 81: 5978-5984. 10.1128/JVI.02650-06.
Bailey A, Mautner V: Phylogenetic relationships among adenovirus serotypes. Virology. 1994, 205: 438-452. 10.1006/viro.1994.1664.
Kovacs GM, Harrach B, Zakhartchouk AN, Davison AJ: Complete genome sequence of simian adenovirus 1: an Old World monkey adenovirus with two fiber genes. J Gen Virol. 2005, 86: 1681-1686. 10.1099/vir.0.80757-0.
Benkο MHB, Russell WC: Family Adenoviridae. Virus Taxonomy VIIth Report of the International Committee on Taxonomy of Viruses. Edited by: van Regenmortel MHV FC, Bishop DHL, Carstens EB, Estes MK, Lemon SM, Maniloff J, Mayo MA, McGeoch DJ, Pringle CR. 2000, Wickner RB: Academic Press, 227-238.
Roy S, Gao G, Clawson DS, Vandenberghe LH, Farina SF, Wilson JM: Complete nucleotide sequences and genome organization of four chimpanzee adenoviruses. Virology. 2004, 324: 361-372. 10.1016/j.virol.2004.03.047.
Wigand R, Mauss M, Adrian T: Chimpanzee adenoviruses are related to four subgenera of human adenoviruses. Intervirology. 1989, 30: 1-9. 10.1159/000150117.
Roy S, Vandenberghe LH, Kryazhimskiy S, Grant R, Calcedo R, Yuan X, Keough M, Sandhu A, Wang Q, Medina-Jaszek CA, et al: Isolation and characterization of adenoviruses persistently shed from the gastrointestinal tract of non-human primates. PLoS Pathog. 2009, 5: e1000503-10.1371/journal.ppat.1000503.
Wevers D, Leendertz FH, Scuda N, Boesch C, Robbins MM, Head J, Ludwig C, Kuhn J, Ehlers B: A novel adenovirus of Western lowland gorillas (Gorilla gorilla gorilla). Virol J. 7: 303-10.1186/1743-422X-7-303.
Heim A, Ebnet C, Harste G, Pring-Akerblom P: Rapid and quantitative detection of human adenovirus DNA by real-time PCR. J Med Virol. 2003, 70: 228-239. 10.1002/jmv.10382.
Madisch I, Harste G, Pommer H, Heim A: Phylogenetic analysis of the main neutralization and hemagglutination determinants of all human adenovirus prototypes as a basis for molecular classification and taxonomy. J Virol. 2005, 79: 15265-15276. 10.1128/JVI.79.24.15265-15276.2005.
Puig M, Jofre J, Lucena F, Allard A, Wadell G, Girones R: Detection of adenoviruses and enteroviruses in polluted waters by nested PCR amplification. Appl Environ Microbiol. 1994, 60: 2963-2970.
Kovacs GM, Davison AJ, Zakhartchouk AN, Harrach B: Analysis of the first complete genome sequence of an Old World monkey adenovirus reveals a lineage distinct from the six human adenovirus species. J Gen Virol. 2004, 85: 2799-2807. 10.1099/vir.0.80225-0.
Lukashev AN, Ivanova OE, Eremeeva TP, Iggo RD: Evidence of frequent recombination among human adenoviruses. J Gen Virol. 2008, 89: 380-388. 10.1099/vir.0.83057-0.
Walsh MP, Chintakuntlawar A, Robinson CM, Madisch I, Harrach B, Hudson NR, Schnurr D, Heim A, Chodosh J, Seto D, Jones MS: Evidence of molecular evolution driven by recombination events influencing tropism in a novel human adenovirus that causes epidemic keratoconjunctivitis. PLoS One. 2009, 4: e5635-10.1371/journal.pone.0005635.
Ebner K, Pinsker W, Lion T: Comparative sequence analysis of the hexon gene in the entire spectrum of human adenovirus serotypes: phylogenetic, taxonomic, and clinical implications. J Virol. 2005, 79: 12635-12642. 10.1128/JVI.79.20.12635-12642.2005.
Gall JG, Crystal RG, Falck-Pedersen E: Construction and characterization of hexon-chimeric adenoviruses: specification of adenovirus serotype. J Virol. 1998, 72: 10260-10264.
Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, Ingersoll R, Sheppard HW, Ray SC: Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999, 73: 152-160.
Katoh K, Kuma K, Toh H, Miyata T: MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005, 33: 511-518. 10.1093/nar/gki198.
Stamatakis A: RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006, 22: 2688-2690. 10.1093/bioinformatics/btl446.
Jeffroy O, Brinkmann H, Delsuc F, Philippe H: Phylogenomics: the beginning of incongruence?. Trends Genet. 2006, 22: 225-231. 10.1016/j.tig.2006.02.003.
Farris JS, Kallersjo M, Kluge AG, Bult C: Constructing a significance test for incongruence. Syst Biol. 1995, 44: 2-
Delsuc F, Brinkmann H, Philippe H: Phylogenomics and the reconstruction of the tree of life. Nat Rev Genet. 2005, 6: 361-375. 10.1038/nrg1603.
The authors thank the Serveis Científico-Tècnics of the University of Barcelona for the sequencing of genomic products, in particular Dr. Ramon Seminago.
The authors declare that they have no competing interests.
CMM and SB collected the sample and isolated the virus. AH cultured the virus and generated purified stocks. CMM and AH purified the viral genome, conceived the amplification and sequencing strategy, and analyzed the sequences obtained. FCA performed the phylogenetic analysis. CMM, FCA and RG analyzed the data and drafted the manuscript. All authors read and approved the final manuscript.