The reticuloendotheliosis viruses (REVs) comprise several closely related amphotropic retroviruses isolated from birds. These viruses exhibit several highly unusual characteristics that have not so far been adequately explained, including their extremely close relationship to mammalian retroviruses, and their presence as endogenous sequences within the genomes of certain large DNA viruses. We present evidence for an iatrogenic origin of REVs that accounts for these phenomena. Firstly, we identify endogenous retroviral fossils in mammalian genomes that share a unique recombinant structure with REVs—unequivocally demonstrating that REVs derive directly from mammalian retroviruses. Secondly, through sequencing of archived REV isolates, we confirm that contaminated Plasmodium lophurae stocks have been the source of multiple REV outbreaks in experimentally infected birds. Finally, we show that both phylogenetic and historical evidence support a scenario wherein REVs originated as mammalian retroviruses that were accidentally introduced into avian hosts in the late 1930s, during experimental studies of P. lophurae, and subsequently integrated into the fowlpox virus (FWPV) and gallid herpesvirus type 2 (GHV-2) genomes, generating recombinant DNA viruses that now circulate in wild birds and poultry. Our findings provide a novel perspective on the origin and evolution of REV, and indicate that horizontal gene transfer between virus families can expand the impact of iatrogenic transmission events.
Retroviruses are characterized by their ability to insert a DNA copy of their genome into the chromosomes of infected cells. Occasionally, retroviruses insert into “germline” cells and are subsequently inherited as host alleles called endogenous retroviruses (ERVs). Vertebrate genomes contain thousands of ERV sequences derived from ancient retroviruses, and these viral sequences serve as molecular “fossils” that can be used to explore how retroviruses have evolved over millions of years. Here we combine an analysis of the retroviral “fossil record” with a phylogenetic and historical investigation to determine the origin of a group of avian retroviruses called reticuloendotheliosis viruses (REVs). We present evidence to demonstrate that rather than arising from natural infections of birds, REVs are in fact derived from mammalian retroviruses that were accidentally introduced into avian hosts during experimental studies of a malaria parasite in the late 1930s. Remarkably, REVs have subsequently inserted into the genomes of two large DNA viruses that infect birds, generating chimeric viruses that now circulate naturally in poultry and wild birds.
Citation: Niewiadomska AM, Gifford RJ (2013) The Extraordinary Evolutionary History of the Reticuloendotheliosis Viruses. PLoS Biol 11(8): e1001642. doi:10.1371/journal.pbio.1001642
Academic Editor: Bill Sugden, University of Wisconsin-Madison, United States of America
Received: February 25, 2013; Accepted: July 19, 2013; Published: August 27, 2013
Copyright: © 2013 Niewiadomska, Gifford. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This study was funded by the Aaron Diamond Research Center. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Abbreviations: DIAV, duck infectious anemia virus; ERV, endogenous retrovirus; FWPV, fowlpox virus; GHV-2, gallid herpesvirus type 2; ICT, interclass transmission; LTR, long terminal repeat; Ma, million years ago; PCR, polymerase chain reaction; REV, reticuloendotheliosis virus; RT, reverse transcriptase; SNV, spleen necrosis virus; TM, transmembrane protein
The reticuloendotheliosis viruses (REVs) comprise several closely related amphotropic retroviruses (family Retroviridae) isolated from birds . The prototypic REV isolate was isolated from a turkey in 1957 . Subsequently, REV infections have been reported in a diverse range of gamebirds (order Galliformes) and waterfowl (order Anseriformes). Infection is associated with a range of disease syndromes, including anemia, immunosuppression, neoplasia, runting, and feathering abnormalities called “nakanuke.” The etiology of REV infection remains enigmatic—although antibodies to REV are widespread in poultry, REV outbreaks occur only sporadically and are relatively rare .
All retroviruses replicate their genomes via a DNA intermediate that is integrated into the nuclear DNA of the host cell and is referred to as a “provirus.” Occasionally, infection of germ cells allows retroviral proviruses to enter the host germline, so that they can be vertically inherited as host alleles, called endogenous retroviruses (ERVs) , a proportion of which end up becoming fixed in the germline. These ancestral retrovirus sequences represent retroviral “fossils” ,, and as such they support “paleovirological” investigations that seek to address the long-term, macroevolutionary history of interaction between hosts and retroviruses ,. In a previous study, phylogenetic analysis of retroviral polymerase (pol) gene sequences revealed that REV groups robustly within the Gammaretrovirus genus, and is closely related to an ERV in the genome of the short-beaked echidna (Tachyglossus aculeatus)—an-egg laying mammal found only in Australia and New Guinea . This discovery reinforced the conclusions of earlier, serological studies, which proposed REVs to have originated in mammals .
Curiously, sequences derived from REV have also been identified in the genomes of two large DNA viruses that naturally infect birds: fowlpox virus (FWPV), a poxvirus , (family Poxviridae), and gallid herpesvirus 2 (GHV-2), a herpesvirus (family Herpesviridae). FWPV infects poultry and wild birds throughout the world, and causes a mild-to-severe, slow developing disease (avian pox) characterized by the formation of proliferative external lesions (dry pox), and diphtheritic lesions in the digestive and respiratory tracts (wet pox) . GHV-2 is the causative agent of Marek's disease, a highly contagious disease of chickens and other galliform birds that is associated with a wide range of clinical syndromes, including neoplasia and paralysis . Clinical disease is not always apparent in infected birds, but mortality rates in susceptible flocks can be very high .
Contamination of both FWPV and Marek's disease vaccines with replication competent REV, leading to outbreaks of REV infection, has been reported on numerous distinct occasions ,. However, only remnant REV sequences, incapable of expressing retrovirus, have been identified in GHV-2 and FWPV vaccine strains (typically a “solo LTR” derived from the long terminal repeat (LTR) regions that flank the provirus) ,. By contrast, FWPV field strains containing near full-length REV proviruses appear to circulate naturally in unvaccinated birds ,–. Recently, a field strain of GHV-2 containing a novel REV LTR insertion was reported .
In this study, we used a combination of PCR-based and in silico screening to explore the origin and evolutionary history of the REV lineage, and to investigate the processes linking exogenous REV isolates with endogenous REV-related sequences in virus and animal genomes.
Paleovirological History of the REV Lineage
To investigate the deeper origins of the REVs, we screened avian and mammalian genome sequence databases (Table S1) for ERV sequences closely related to REV (Table 1). Screening of 42 mammalian genomes identified numerous ERV loci that disclosed highly significant similarity to one or more REV coding domains, but none that matched closely to REV across the entire coding region of the genome. We found that all mammalian ERVs exhibiting a high degree of sequence similarity to REV in the gag-pol domain exhibited no such similarity in env, and vice versa. This can be assumed to reflect the recombinant genome structure of REV ,, comprising a Gammaretrovirus gag-pol domain fused to an env domain that is more commonly associated with the Betaretrovirus genus (although it also occurs in some other Gammaretroviruses, also considered to be recombinants ). No ERV loci closely related to REV were detected in avian genomes. We did identify numerous avian ERVs that disclosed weak similarity to REV in pol (30–40% amino acid identity). However, phylogenetic analysis revealed these ERVs to be derived from ancient, highly degenerated ERV lineages that were clearly distinct from modern Gammaretroviruses (Figure 1).
Shaded boxes indicate taxa that are known to occur as exogenous retroviruses. Brackets to the right indicate major lineages (note: an integrated taxonomy of exogenous and ERVs has yet to be established by the International Committee on Taxonomy of Viruses, and the groupings shown here are propositional). Associations of retrovirus groups and individual retroviral taxa with avian and mammalian hosts are indicated, as shown in the key. The phylogeny shown was constructed using NJ and a multiple sequence alignment spanning 140 amino acid residues in the reverse transcriptase protein (RT), and is midpoint rooted for display purposes. To obtain putative protein sequences for ERVs, frameshifting indels were inferred and removed, and the resulting nucleotide sequence was conceptually translated. Asterisks indicate clades with bootstrap support >90% in both NJ and maximum likelihood (ML) trees, based on 1,000 bootstrap replicates. The scale bar indicates evolutionary distance in substitutions per site. Table S2 provides details of all the ERVs and exogenous retrovirus taxa shown in the phylogeny.
In phylogenies based on reverse transcriptase (RT), avian REV isolates cluster tightly with a previously described ERV sequence derived from the short-beaked echidna genome . During a polymerase chain reaction (PCR)–based investigation of ERV diversity in Malagasy mammals, we serendipitously identified additional ERV RT sequences that grouped within this clade, in the genomes of two Malagasy carnivore species: the ring-tailed mongoose (Galidia elegans) and the narrow-striped mongoose (Mungotictis decemlineata). We recovered near complete proviral genome sequences for all three REV-related ERVs (hereafter referred to as echidna-ERV, Galidia-ERV, and Mungotictis-ERV) (Figure 2a), revealing that they exhibit similarity to REV throughout the entire internal coding region of the genome. Crucially, echidna-ERV, Galidia-ERV, and Mungotictis-ERV grouped robustly with REV isolates in phylogenies constructed using both the pol and env coding domains (Figure 3a and 3b), establishing that they share a common, recombinant ancestor with these viruses. Thus, ERVs belonging to the REV-lineage do occur in the genomic fossil record of mammals, but as with certain other retrovirus groups, such as foamy viruses and lentiviruses, they are relatively rare ,.
The schematic in panel (a) shows the genome structure of REV and SNV, a near full-length REV insertion in the FWPV genome, and the mammalian ERVs Echidna-ERV and Galidia-ERV. Percentage sequence identity to SNV, at the amino acid level, is shown for the putative Gag, Pol, and Env polyproteins of Echidna-ERV and Galidia-ERV. Proviral coding regions that disclose homology to Gammaretroviruses are shown in green, whereas those that disclose homology to Betaretroviruses are shown in blue. ORFs flanking the REV insertion in FWPV are in yellow. Panel (b) summarizes the genomic data used to estimate the minimum age of REV-related ERV insertions in Malagasy carnivore genomes. A time-scaled Carnivora phylogeny (based on Nyakatura et al. ) is shown on the left, with Malagasy carnivores shaded. A corresponding schematic on the right shows the genomic locus at which an orthologous ERV insertion was identified in a subset of Malagasy carnivores. Boxes represent the env gene (blue) and 3′ LTR sequences (green = U3; dark grey = R; light grey = U5). The adjacent black line represents flanking genomic DNA, spanning 238 nucleotides, obtained from the striped mongoose (Mungotictis decemlineata) and ring-tailed mongoose (Galidia elegans) genomes in our study, and aligned to a homologous genomic region (lacking a proviral insertion) in the cat (Felis catus), dog (Canis familiaris), and ferret (Mustela furo) genomes. An orthologous ERV insertion was detected in M. decemlineata and G. elegans genomes, but not in the more distantly related Fossa (Cryptoprocta ferox), indicating that germline invasion occurred between 18 and 8 Ma. Genetic data indicate that all Malagasy carnivores are derived from a single founder population that colonized Madagascar ~19 Ma ; thus, invasion of the Malagasy carnivore germline occurred in Madagascar. The nucleotide sequence alignment on which the schematic in panel (b) is based on is shown in Figure S1. Abbreviations: RV, retrovirus; Kb, Kilobases; ORF, open reading frame; PBS, primer binding site; Pro, proline; Thr, threonine; LTR, long terminal repeat; U3, unique three prime region; R, repeat region; U5, Unique five prime region; RT, reverse transcriptase; SU, surface protein; TM, transmembrane protein; M.dec, Mungotictis decemlineata; G.ele, Galidia elegans.
Panels (a) and (b) show ML phylogenies constructed from alignments of Gamma- and Betaretrovirus protein sequences. The phylogeny in panel (a) was constructed from an alignment spanning 157 residues of the RT protein encoded by pol, whereas the phylogeny in panel (b) was constructed from an alignment spanning 153 residues of the TM domain in the polypeptide encoded by env. Asterisks on internal nodes indicate ML bootstrap support >95%, (based on 1,000 bootstrap replicates). Asterisks beside taxa names indicate ERV families identified in this study. Open triangles indicate ERV lineages for which env genes were not identified. Scale bars indicate evolutionary distance in substitutions per site. Brackets to the right indicate genus designations, and viruses previously identified as Gamma- and Betaretrovirus (γ-β) recombinants. Table S2 provides details of all the ERVs and exogenous retrovirus taxa shown in the phylogeny. Abbreviations: RV, retrovirus; MoMLV, Moloney murine leukemia virus; FeLV, feline leukemia virus; GaLV, gibbon ape leukemia virus; KoRV, koala retrovirus; BAEV, baboon endogenous virus; SMRV, squirrel monkey retrovirus; TvERV, Trichosurus vulpecula endogenous retrovirus; JSRV, Jaagsiekte sheep retrovirus; SRV, simian retrovirus; MMTV, mouse mammary tumor virus.
PCR results suggested that all three ERVs are low copy number (1–2 proviruses) in their host species. Along with other factors, such as the relatively short length of LTRs, this precluded the confident use of molecular clock-based approaches to date the echidna-ERV, Galidia-ERV, and Mungotictis-ERV insertions. Notably, however, internal coding regions in all three ERVs were relatively intact (although echidna-ERV has a large deletion in region of the env gene encoding the surface (SU) glycoprotein (Figure 2a)).
Using a ligation-mediated PCR method, we recovered matching flanking insertion sites for Galidia-ERV and Mungotictis-ERV, confirming that REV-like viruses occur as orthologous insertions in distinct Malagasy mongoose species. This finding indicates that REV-like viruses entered the germline Malagasy mammals prior to the divergence of Galidia and Mungotictis ~8 million years ago (Ma) , (Figure 2b). REV-related ERVs were not detected in the more distantly related fossa (Cryptoprocta ferox). Together these findings establish that the entire REV lineage—including both mammalian and avian isolates—derives from a common founder that was generated by recombination, and circulated among mammals during the Miocene Epoch (~23–5 Ma).
Origin of the Avian REVs
Given that the REV lineage clearly originates in mammals, we decided to investigate the origins of the avian REVs in greater detail. We reviewed all published reports of REV outbreaks, and sought to obtain any available archived samples (Table 2). Among REV isolates that had not previously been sequenced, we were only able to obtain one (duck infectious anemia virus (DIAV)), for which all samples had not been exhausted or destroyed. We sequenced the complete DIAV genome and constructed ML phylogenies using all available REV sequence data. Phylogenies were constructed using alignment partitions representing (i) a conserved region of the pol coding domain (Figure 4a), (ii) the complete internal coding region of the viral genome (Figure 4b), and (iii) LTR sequences (Figure 5).
ML phylogenies constructed using (a) an alignment spanning residues 183–481 of the Pol polyprotein (DIAV coordinates) and containing REV and mammalian gammaretroviruses sequences and (b) a nucleotide alignment of the entire internal coding region of full-length avian isolates. The tree in panel (b) indicates the number of strain-specific, nonsynonymous mutations estimated to have occurred in the nucleocapsid (NC), capsid (CA), matrix (MA), protease (PR), RT, RNase-H (RH), integrase (IN), surface (SU), and TM genes of the exogenous isolate HA9901. Asterisks on internal nodes indicate ML bootstrap support >95%. All trees are midpoint rooted for display purposes. Scale bars indicate evolutionary distance in substitutions per site. Taxa labels include sequence accession numbers, and in panel (b) two-letter ISO country codes enclosed by brackets indicating the country of sampling. Further details of REV sequences included in these trees can be found in Table S4.
ML phylogenies constructed using an alignment of REV LTR sequences. Asterisks on internal nodes indicate ML bootstrap support >95%. The phylogeny is midpoint rooted for display purposes. Scale bars indicate evolutionary distance in substitutions per site. Taxa labels include two-letter ISO country codes indicating the country of sampling.. Taxa labels include accession numbers and two-letter ISO country codes enclosed by brackets indicating country of sampling. Where appropriate, FWPV and GHV-2 strain designations are shown in bold. Further details of REV sequences used in the tree, and an alignment figure highlighting lineage-specific LTR indels, can be found in Table S4 and Figure S2, respectively.
All three phylogenies consistently disclosed three major lineages. The first was comprised of spleen necrosis virus (SNV) and DIAV. Both these viruses were isolated from ducks that were experimentally infected with Plasmodium lophurae (SNV in 1959  and DIAV in 1972 ). The report describing the isolation of DIAV concluded that P. lophurae stocks were the source of infection, and demonstrated contamination of stocks in five different laboratories. Sequencing revealed that SNV and DIAV are highly related (~98% nucleotide identity), despite being isolated 13 years apart, establishing that contaminated stocks have been the source of multiple outbreaks of retroviral infection in P. lophurae–infected ducks, dating back as far as 1959, and likely earlier –.
A second clade comprised the REV insertion in FWPV and exogenous REV isolates obtained independently in different countries, including the prototypic REV isolates isolated in the United States . In addition, LTR phylogenies revealed this clade to include insertions present in two distinct GHV-2 strains: an attenuated lab strain (RM1 ), and a field strain (GX0101 ). Virus in this clade exhibit remarkably little genetic variation overall, despite having apparently been maintained in the avian population for at least 50 years . It thus appears unlikely that the exogenous REV isolates in this clade are spreading primarily through horizontal transmission of infectious retrovirus (since this would be expected to generate greater nucleotide sequence diversity among isolates). Instead, phylogenies suggest that exogenous retroviruses are being expressed from a stable FWPV-REV vector that circulates among domestic and wild birds. This would explain why antibodies to REV are widespread in poultry (see Table S5), and why REV infections occur not only in association with contaminated vaccines, but also in wild birds and unvaccinated commercial flocks . Revealingly, several reports describe FWPV or undiagnosed pox-like infections occurring in bird populations shortly prior to the occurrence of REV outbreaks ,–.
The third clade comprised the exogenous REV isolate HA9901, from China , as well as LTR sequences obtained from the JM-Hi3 strain of GHV-2, and a REV plasmid (pREVA6 ). This clade is robustly supported in bootstrapped phylogenies, and the presence of unique, shared indels in LTRs provides further evidence of common ancestry (Figure S2). These observations establish that HA9901 shares a common history with pREVA6, which can ultimately be traced back to the prototypic REV specimen ,. Interestingly, HA9901 has acquired numerous nonsynonymous mutations, consistent with ongoing replication as an exogenous retrovirus (Figure 4b).
The data presented in this study unequivocally demonstrate that REVs derive from a retrovirus that circulated in ancestral mammals, and originated through recombination more than 8 Ma. Furthermore, the extremely low genetic diversity observed among all avian REV isolates and sequences indicates a very recent origin for REV in birds (Figure 4 and 5). In previous studies it has generally been assumed that the REVs are a group of bona fide avian retroviruses that circulate in wild bird populations. However, phylogenetic evidence indicates that successful transmission of retroviruses, poxviruses, and herpesviruses across host classes is extremely rare, if indeed it occurs at all ,,. While such “long-distance” transmission events, leading to productive virus replication in the new host, likely do occur at an appreciable frequency for these viruses (particularly, for example, when the recipient host is immunocompromised), unless the transmitted virus is able to spread efficiently from host-to-host in the new species, these instances will typically represent evolutionary dead-ends .
Since REVs clearly originate in mammals, and all avian REVs are highly related, the entire avian REV lineage almost certainly derives from a single founder. Phylogenies rooted on mammalian REVs unambiguously place the SNV/DIAV lineage in a basal position relative to the FWPV-REV and HA9901 clades (Figure 4a). This is most readily reconciled with a scenario wherein REVs originated in P. lophurae experiments, and subsequently inserted into the FWPV and GHV-2 genomes (Figure 6). Importantly, this hypothesis of REV origin and evolution is not only consistent with the REV phylogeny, but also with the entire recorded history of REV-associated disease (Figure 7, Table 2), accounting for the disappearance of the SNV/DIAV lineage since the 1980s (when P. lophurae stocks were exhausted—see below), and the limited genetic diversity observed among all avian REVs (since relatively few virus replication cycles would be expected to separate all isolates). Moreover, this scenario accounts for the anomaly of retroviral interclass transmission, because it occurs in an experimental context wherein a pathogen (P. lophurae) is being deliberately adapted to a foreign host species.
A schematic showing the three possible scenarios via which the ancestor of REV could have crossed from birds (class Aves) into mammals (class Mammalia), assuming a maximum of one inter-class transmission (ICT) event in total. For each of the three scenarios shown, the phylogenetic relationships between REV isolates that would be expected to arise as a result are indicated (all phylogenies are rooted on the mammalian ancestor of avian REVs). A REV founder strain could conceivably have been transmitted from mammals to birds after first inserting into the genome of FWPV (panel a) or GHV-2 (panel b). However, only a scenario in which the SNV and DIAV lineage were established first (panel c)—as would be expected to occur if P. lophurae contamination enabled the iatrogenic emergence of virus—is compatible with the relationships observed in rooted phylogenetic trees (see Figure 4a). Abbreviations: FWPV, fowlpox virus; GHV-2, gallid herpesvirus 2.
A schematic representation of REV evolutionary history is shown, summarizing our hypothesis regarding the origin and evolution of the three major avian REV lineages (SNV/DIAV, REV/FWPV-REV, and HA9901) from a mammalian retrovirus ancestor that originated in the Cenozoic Era. REV-associated events (i.e., outbreaks of REV-associated disease, isolation of new REV strains, or identification of REV-containing DNA virus strains) reported in the literature have been mapped onto this schematic, as indicated in the key. Numbers shown above key symbols refer to Table 2, where details of the associated publication or report can be found. The broken scale bar shows time in years A.D. to the right of the break and Ma to the left of the break. A shaded background region indicates the time window for invasion of FWPV genome following iatrogenic introduction into poultry (assuming that reports of REV sequences in FWPV vaccine strains lyophilized in 1949  are accurate). Abbreviations: REV, reticuloendotheliosis virus; SNV, spleen necrosis virus; DIAV, duck infectious anemia virus; FWPV, fowlpox virus; GHV-2, gallid herpesvirus 2; FWPV-REV, Fowlpox virus with REV insertion.
P. lophurae has only been isolated once, in June 1937, in the New York Zoological Park (now Bronx Zoo), by Lowell T. Coggeshall. Coggeshall, who was then working for the Rockefeller Foundation, was searching for a parasite that could serve as an experimental model system for malaria research. In 1935, Émile Brumpt of the Pasteur Institute had identified Plasmodium gallinaceum, a parasite causing malarial disease in poultry, during an excursion to Ceylon (now Sri Lanka) . However, P. gallinaceum could not be introduced to the United States due to strict quarantine regulations against importation of poultry pathogens . Reasoning that other avian species from the same geographic region might harbor a similar parasite, Coggeshall screened some South East Asian bird species that had been introduced to the New York Zoological Park in the 1920s by ornithologist Lee Saunders Crandall . This led to the identification of a plasmodium in the blood of a Borneo firebacked pheasant (Lophura igniti igniti), which proved transmissible to very young chickens .
Stocks of this parasite, designated Plasmodium lophurae, were maintained by serial passage in chicken, duck, and turkey chicks, with 25 passages reported as of 1938 . Published reports suggest that contaminating virus was present from an early stage; a 1941 study of P. lophurae noted that anemia in infected animals appeared to be decoupled from parasite replication, indicating the presence of a second infectious agent . A study a few years later confirmed the presence of an additional “filterable agent”—the cause of a lethal anemia—in P. lophurae–infected poultry , and in 1959 William Trager identified this agent as SNV . Subsequently, SNV-like viruses were isolated from P. lophurae–infected ducks on multiple distinct occasions (Table 2) ,[31. The role of P. lophurae stocks as a source of infection appears to have gone unappreciated prior to the isolation of DIAV in 1972 . But while the associated study concluded that “DIAV has been an unrecognized companion of P. lophurae for many years,” the assumption remained that the contaminating virus was a natural pathogen of ducks.
Research on P. lophurae effectively ceased in the 1980s, when stocks could no longer be replenished. The organism has never subsequently been identified, and thus remains an enigma in many respects. Expeditions to Borneo have been mounted with the express purpose of obtaining further isolates, but these failed to identify the parasite in populations of wild birds . Since P. lophurae stocks ran out, no further viruses belonging to the SNV/DIAV lineage have been isolated, consistent with the hypothesis that contaminated stocks were the principle reservoir of infection for these viruses.
It remains unclear whether the progenitors of avian REVs were present in the animal from which P. lophurae was originally obtained or were introduced from an external source during serial passage. However, since none of the mammalian species that might be considered likely sources of contamination in a lab environment (i.e., mouse, rat, rabbit, guinea pig) appear to harbor truly REV-like viruses in exogenous or endogenous forms, whereas more exotic mammalian species do, cross-species transmission or contamination within the setting of the zoological park is an attractive hypothesis. Notably, we have identified REV-related ERVs in mammalian groups (Malagasy carnivores and Australian monotremes) that inhabit highly distinct and relatively isolated biogeographic regions, separated from one another by large expanses of ocean. This suggests that infection has been widespread in the past and that chiropteran (bat) vectors were likely at least partly involved in the spread of virus.
It also remains unclear precisely when and how the REV insertions in FWPV and GHV-2 genome were generated. REV could presumably have spread from birds experimentally infected with P. lophurae and into the wider environment either before or after inserting into a DNA virus vector. Notably, research on malaria was prioritized in the United States during World War II, and P. lophurae stocks were distributed to laboratories throughout the country (see Table S6) for experimental vaccine and drug research. During this period the poultry industry was scaling rapidly, and the first avian virus vaccines were being commercially developed (including live FWPV vaccines, based on attenuated virus strains grown in embryonated eggs ,). REV sequences have been reported in FWPV vaccines lyophilized in 1949 , suggesting that insertion had already occurred by this time. Unfortunately, however, this inference is subject to some incertitude, since it is based solely on PCR from a single archived sample, and no lyophilized material remains for study (Table 2).
The creation of Mareks disease vaccines became a priority in the United States during the 1950s, in response to devastating outbreaks of an apparently new, acute form of the disease . However, effective vaccines were not produced until after the first avian cell culture systems were established in the 1960s. These in vitro systems were key to the eventual development of vaccines based on (i) attenuated GHV-2 strains and (ii) the closely related herpesvirus of turkeys (HVT). Both of these vaccines were later discovered to be contaminated with REV. In previous studies it has generally been assumed that REV insertions into the GHV-2 genome originated in the distant evolutionary past ,, (although it is recognized that at least some were generated recently during in vitro attenuation ,). By contrast, our data suggest that all REV insertions into GHV-2 have been generated recently.
In 1960s and 1970s REV provided an experimental model for retrovirologists , and was sometimes used to transform avian cells . Thus it is likely that the emergence of avian cell culture systems was accompanied by the spread of REV as a contaminant. Interestingly, dissemination of REV genetic material appears to be ongoing; REV is apparently being maintained as an insertion in naturally circulating FWPV-REV, and field strains of GHV-2 containing novel REV LTR insertions have recently been reported ,. Furthermore, we show that the recently described exogenous REV isolate HA9901  shares a common history with REV plasmid pREVA6, which was in turn derived from the original tissue sample from which prototypic REV strains were isolated  (Figure 5, Figure S2). Thus it appears that in China, REV-contaminated materials may have given rise to independently circulating infectious retrovirus. The processes driving REV dissemination warrant further exploration, as does the potential role of co-opted REV sequences in altering the in vivo properties of FWPV and/or GHV-2.
In conclusion, historical, phylogenetic, and paleovirological evidence supports a scenario wherein REVs originated as mammalian retroviruses that were iatrogenically introduced into avian hosts, and subsequently integrated into the FWPV and GHV-2 genomes, generating recombinant DNA viruses that now circulate in wild birds and poultry. These data provide the first evidence that horizontal gene transfer between virus families can expand the impact of iatrogenic transmission events, raising questions about the potential, unintended impacts of live, recombinant vector vaccines. Broader surveillance of viral genetic diversity should be prioritized, so that the unintended consequences of experimental procedures on viral ecology and evolution can be better assessed and limited.
Materials and Methods
Screening in Silico
PERL scripts and the BLAST+ program suite were used to perform in silico screening of sequence databanks for sequences homologous to REV. We screened complete and low coverage whole genome sequence data representing 10 avian and 42 mammalian species (Table S1) and all poxvirus and herpesvirus-derived sequence data available in GenBank as of July 1, 2012. The noncoding nucleotide sequences (LTR and leader) and translated open reading frames (ORFs) (Gag, Pol, Env) of REV (FJ439119.1) were used as “probes” for in silico screening. Sequences that matched probes with high statistical significance (i.e., expect (e) values 0.001) were extracted and compared to a library of reference retroviral genomes (see Table S2), again using BLAST. The results of this “reciprocal” BLAST were examined, and the phylogenetic relationships of ERV loci that disclosed higher similarity to REV than to any other retroviral reference were investigated using the neighbor joining (NJ) algorithm implemented in PAUP . NJ trees revealed that among all the ERV loci identified by screening, ERVs in the European hedgehog (Erinaceus europaeus) and cape hyrax (Procavia capensis) genomes were most closely related to REV in the gag-pol and env genes, respectively (unpublished data). The median reciprocal BLAST bit score for these two subsets of ERVs was used to establish a threshold bit score for discriminating REV-related coding sequences from those of other ERVs.
Tissue Samples, Virus, and Cell Culture
Frozen tissue samples from Malagasy carnivores (Cryptoprocta ferox, Galidia elegans, Mungotictis decemlineata) were obtained from the American Museum of Natural Historys cryogenic collection. Frozen spleen tissue samples were obtained from a deceased echidna (Tachyglossus aculeatus) at the Philadelphia Zoo. Chicken embryonic fibroblasts, SL-29 cells (ATCC#: CRL1590), were maintained in DMEM medium (Life technologies) supplemented with 5% fetal bovine serum, 5% tryptose phosphate broth, penicillin (100 U/ml), and streptomycin (100 mg/ml). An aliquot of 400 ul of DIAV (ATCC #: VR775) was inoculated onto 30% confluent SL-29 cells in six-well plates. Media was changed after 2 days, and cells were allowed to grow for a total of 5 days. After 5 days, cells were harvested and genomic DNA was extracted.
PCR and Sequencing
Genomic DNA was extracted from tissue samples and SL-29 cells using the AllPrep dual DNA/RNA extraction kit (QIAGEN). Initial PCR amplification of endogenous retroviral fragments was performed using PCR primers (Integrated DNA Technologies) directed against two highly conserved motifs in retroviral protease (PR) and RT proteins. After initial sequencing of this genomic region, a combination of gene-specific primers and degenerate primers were used to amplify the remaining regions of the REV genomes found in Galidia elegans, Mungotictis decemlineata, and Tachyglossus aculeatus. LTR regions and genomic insertion sites were amplified and cloned by ligation-mediated PCR, using the GenomeWalker Universal kit (Clontech). For complete genome sequencing of DIAV, primers were based on equivalent targets in REV and SNV and were used to amplify multiple overlapping regions of the DIAV genome. A list of primer sequences, the genomes on which they were used, and their coordinates (based on alignment to the DIAV genome) can be found in Table S3. Basic PCR conditions were used for almost all reactions (denaturation at 95°C for 2 min, followed by 30 cycles of 94°C for 15 s, 55°C for 30 s, and 68°C for 1 min, final elongation for 10 min), although annealing temperatures and elongation times varied depending on the primers used (details available on request). For all reactions, gel-resolved amplicons were excised from 1% agarose gels and purified using the Qiaquick kit (QIAGEN) before TA cloning into pCR2.1 (Life Technologies, La Jolla, CA) and sequencing. All sequence analysis was performed by the GeneWiz commercial sequencing facility (GeneWiz, South Plainfield, NJ). Sequences obtained in this study have been submitted to Genbank under the following accession numbers: DIAV (KF313137); Echidna-ERV (KF313136), and Galidia-ERV (KF313135).
Sequence Data and Phylogenetic Analyses
Retroviral “pan-genus” phylogenies were constructed from an alignment of the highly conserved RT and transmembrane (TM) peptides. Sequences derived from the retroviral reference library (Table S2) were included, as well as a selection of the best matching, uncharacterized ERVs from in silico and PCR screening. For both genes, ProtTest was used to select the best fitting amino acid substitution matrix from a range of 96 different combinations of models and rate heterogeneity parameters, based on the Akaike information criterion (AIC) . The best fitting model for RT was rtREV , with gamma distributed rate heterogeneity (rtREV+Γ); for TM it was HIVw .
Phylogenetic investigation of within-REV variation was conducted using both peptide and nucleotide sequence data. We obtained all published REV sequence data from Genbank. Sequences shorter than 100 bp were excluded. The location and year of sampling, and host species associations, were extracted from the Genbank file or from an associated publication (Table S4). All sequences were profile aligned to a full genome reference (SNV; DQ003951.1). So that the phylogenetic relationships of all available sequences could be investigated, phylogenies were constructed for a range of alignment partitions: (i) complete genome, (ii) LTRs, (iii) gag, (iv) pol, and (v) env. Each partition was examined for evidence of recombination using GARD  and SplitsTree . One full genome sequence (GQ375848.1) appeared to be recombinant and was subsequently removed from our dataset. We used ProtTest and ModelTest to select the best fitting amino acid and/or nucleotide substitution matrices for each alignment partition. The best fitting model for all nucleotide alignments was the general time reversible model  with a proportion of invariable sites and a gamma-shaped rate variation across sites (GTR+I+G). The best fitting models for amino acid alignments were; Gag, JTT+I; Pol, FLU+G; Env, JTT. The ML phylogeny was constructed using RAXML , with 1,000 nonparametric bootstrap replicates. A REV ancestral sequence was reconstructed using PAML .
To systematically review REV-related literature, electronic searches of PubMed/Medline, JSTOR, Mendley, Scopus, Web of Science, and WorldCat were conducted in July 2012. Keywords used to search databases were “Reticuloendotheliosis Virus,” “Duck Infectious Anemia Virus,” “Spleen Necrosis Virus,” and “Chick Sync[i/y]tial Virus.” We restricted our search to papers with titles and abstracts available in English. The following data were searched for in texts: year of virus isolation, virus association, origin of isolation, animal status, secondary disease association, place of isolation, and the methods of isolation or detection. A completed PRISMA checklist and flow diagram for this systematic literature review can be found in Text S1.
A nucleotide alignment of orthologous ERV insertion sites in the Galidia elegans (Galidia-ERV-1,2) and Mungotictis decemlineata (Mungotictis-ERV-1,2) genomes. The alignment illustrated spans the 3′LTR and 3′ end of env of the orthologous REV-related ERV insertion in these two species and 238 bp of flanking genomic DNA (shown in grey). Flanking DNA is shown aligned to a homologous genomic sequence identified in the cat (Felis catus), dog (Canis familiaris), panda (Ailuropoda melanoleuca), and ferret (Mustela furo) genomes.
An alignment of REV LTR sequences, showing the presence of unique shared indels (insertions and deletions) that support the monophyletic relationship of the three sequences highlighted in gray, which include (i) the HA9901 strain of REV, (ii) REV plasmid (pREVA6), and (iii) a REV LTR insertion present in the JM-Hi3 strain of GHV-2. Shared indels are indicated by boxes.
Avian and mammalian whole genome sequences screened for REV-related ERVs.
Retroviral reference sequences used in the study. Annotated reference sequences representing newly described ERVs have been made available online (http://saturn.adarc.org/paleo/).
Primer Coordinates. aThe coordinates of primers are shown, based on alignment to DIAV reference sequence (Accession Number KF313137). Abbreviations: F, forward; R, reverse; REV, reticuloendotheliosis virus; LTR, long terminal repeat.
Abbreviations: REV, reticuloendotheliosis virus; FWPV-REV, fowlpox virus with REV insertion; GHV-2-REV, gallid herpesvirus type 2 with REV insertion; LTR, long terminal repeat.
aState or prefecture/two-letter ISO country code.
Asterisks indicate studies sponsored by the US Office of Scientific Research and Development (OSRD).
PRISMA checklist and flow diagram.
We thank the American Museum of Natural History and Philadelphia Zoo for providing tissue samples, and Eric Delwart, Karel A. Schat, Irwin W. Sherman, Michael Skinner, Greg Towers, Michael Tristem, and Robin Weiss for helpful discussions.
The author(s) have made the following declarations about their contributions: Conceived and designed the experiments: AMN RJG. Performed the experiments: AMN RJG. Analyzed the data: AMN RJG. Contributed reagents/materials/analysis tools: RJG. Wrote the paper: RJG.
- 1. Purchase HG, Ludford C, Nazerian K, Cox HW (1973) A new group of oncogenic viruses: reticuloendotheliosis, chick syncytial, duck infectious anemia, and spleen necrosis viruses. J Natl Cancer Inst 51: 489–499.
- 2. Robinson FR, Twiehaus MJ (1974) Isolation of the avian reticuloendothelial virus (strain T). Avian Dis 18: 278–288. doi: 10.2307/1589142
- 3. Fadly AM, Zavala G, Witter RL (2011) Reticuloendotheliosis. In: Saif YM, editor. Diseases of poultry. New York: John Wiley Sons.
- 4. Weiss RA (2006) The discovery of endogenous retroviruses. Retrovirology 3: 67.
- 5. Stoye JP (2012) Studies of endogenous retroviruses reveal a continuing evolutionary saga. Nat Rev Microbiol 10: 395–406. doi: 10.1038/nrmicro2783
- 6. Gifford R, Tristem M (2003) The evolution, distribution and diversity of endogenous retroviruses. Virus Genes 26: 291–315. doi: 10.1023/a:1024455415443
- 7. Feschotte C, Gilbert C (2012) Endogenous viruses: insights into viral evolution and impact on host biology. Nat Rev Genet 13: 283–296. doi: 10.1038/nrg3199
- 8. Patel MR, Emerman M, Malik HS (2011) Paleovirology - ghosts and gifts of viruses past. Curr Opin Virol 1: 304–309. doi: 10.1016/j.coviro.2011.06.007
- 9. Martin J, Herniou E, Cook J, ONeill RW, Tristem M (1999) Interclass transmission and phyletic host tracking in murine leukemia virus-related retroviruses. J Virol 73: 2442–2449.
- 10. Barbacid M, Hunter E, Aaronson SA (1979) Avian reticuloendotheliosis viruses: evolutionary linkage with mammalian type C retroviruses. J Virol 30: 508–514.
- 11. Isfort RJ, Witter R, Kung HJ (1994) Retrovirus insertion into herpesviruses. Trends Microbiol 2: 174–177. doi: 10.1016/0966-842x(94)90668-8
- 12. Hertig C, Coupar BE, Gould AR, Boyle DB (1997) Field and vaccine strains of fowlpox virus carry integrated sequences from the avian retrovirus, reticuloendotheliosis virus. Virology 235: 367–376. doi: 10.1006/viro.1997.8691
- 13. Hansen W (2001) Avian pox. Field manual of wildlife diseases general field procedures and diseases of birds biological resources division information and technology report 1999–2001. U.S.A. pp. 163–170.
- 14. Morrow C, Fehler F (2004) Mareks disease a worldwide problem. Mareks disease: an evolving problem. London: Elsevier Academic Press. pp. 49–61.
- 15. Yuasa N, Yoshida I, Taniguchi T (1976) Isolation of a reticuloendotheliosis virus from chickens inoculated with Mareks disease vaccine. National Institute of Animal Health Quarterly 16: 141–151.
- 16. Singh P, Schnitzlein WM, Tripathy DN (2003) Reticuloendotheliosis virus sequences within the genomes of field strains of fowlpox virus display variability. Journal of Virology 77: 5855–5862. doi: 10.1128/jvi.77.10.5855-5862.2003
- 17. Diallo IS, Mackenzie MA, Spradbrow PB, Robinson WF (1998) Field isolates of fowlpox virus contaminated with reticuloendotheliosis virus. Avian Pathology: Journal of the WVPA 27: 60–66. doi: 10.1080/03079459808419275
- 18. Moore KM, Davis JR, Sato T, Yasuda A (2000) Reticuloendotheliosis virus (REV) long terminal repeats incorporated in the genomes of commercial fowl poxvirus vaccines and pigeon poxviruses without indication of the presence of infectious REV. Avian Diseases 44: 827–841. doi: 10.2307/1593055
- 19. Biswas SK, Jana C, Chand K, Rehman W, Mondal B (2011) Detection of fowl poxvirus integrated with reticuloendotheliosis virus sequences from an outbreak in backyard chickens in India. Vet Ital 47: 147–153.
- 20. Cui Z, Zhuang G, Xu X, Sun A, Su S (2010) Molecular and biological characterization of a Mareks disease virus field strain with reticuloendotheliosis virus LTR insert. Virus Genes 40: 236–243. doi: 10.1007/s11262-009-0437-z
- 21. Kewalramani VN, Panganiban AT, Emerman M (1992) Spleen necrosis virus, an avian immunosuppressive retrovirus, shares a receptor with the type D simian retroviruses. J Virol 66: 3026–3031.
- 22. Koo HM, Gu J, Varela-Echavarria A, Ron Y, Dougherty JP (1992) Reticuloendotheliosis type C and primate type D oncoretroviruses are members of the same receptor interference group. J Virol 66: 3448–3454.
- 23. van der Kuyl AC, Dekker JT, Goudsmit J (1999) Discovery of a new endogenous type C retrovirus (FcEV) in cats: evidence for RD-114 being an FcEV(Gag-Pol)/baboon endogenous virus BaEV(Env) recombinant. J Virol 73: 7994–8002.
- 24. Katzourakis A, Gifford RJ, Tristem M, Gilbert MT, Pybus OG (2009) Macroevolution of complex retroviruses. Science 325: 1512. doi: 10.1126/science.1174149
- 25. Katzourakis A, Tristem M, Pybus OG, Gifford RJ (2007) Discovery and analysis of the first endogenous lentivirus. Proc Natl Acad Sci U S A 104: 6261–6265. doi: 10.1073/pnas.0700471104
- 26. Yoder AD, Burns MM, Zehr S, Delefosse T, Veron G, et al. (2003) Single origin of Malagasy Carnivora from an African ancestor. Nature 421: 734–737. doi: 10.1038/nature01303
- 27. Nyakatura K, Bininda-Emonds OR (2012) Updating the evolutionary history of Carnivora (Mammalia): a new species-level supertree complete with divergence time estimates. BMC Biol 10: 12. doi: 10.1186/preaccept-5398900576110216
- 28. Trager W (1959) A new virus of ducks interfering with development of malaria parasite (Plasmodium lophurae). Proc Soc Exp Biol Med 101: 578–582. doi: 10.3181/00379727-101-25023
- 29. Ludford CG, Purchase HG, Cox HW (1972) Duck infectious anemia virus associated with Plasmodium lophurae. Exp Parasitol 31: 29–38. doi: 10.1016/0014-4894(72)90044-6
- 30. Dearborn EH (1946) Filtrable agents lethal for ducks. Proc Soc Exp Biol Med 63: 48. doi: 10.3181/00379727-63-15488
- 31. Li J, Calnek BW, Schat KA, Graham DL (1983) Pathogenesis of reticuloendotheliosis virus infection in ducks. Avian Diseases 27: 1090–1105. doi: 10.2307/1590209
- 32. Ludford CG, Corwin RM, Cox HW, Sheldon TA (1969) Resistance of ducks to a Plasmodium sp. induced by a filterable agent. Military Medicine 134: 1276–1283.
- 33. Witter RL, Li D, Jones D, Lee LF, Kung HJ (1997) Retroviral insertional mutagenesis of a herpesvirus: a Mareks disease virus mutant attenuated for oncogenicity but not for immunosuppression or in vivo replication. Avian Dis 41: 407–421. doi: 10.2307/1592197
- 34. Hanson J, Howell J (1979) Suspected reticuloendotheliosis in turkeys with cutaneous lesions reminiscent of fowl pox. Can Vet J 20: 161–164.
- 35. Dren CN, Saghy E, Glavits R, Ratz F, Ping J, et al. (1983) Lymphoreticular tumour in pen-raised pheasants associated with a reticuloendotheliosis like virus infection. Avian Pathol 12: 55–71. doi: 10.1080/03079458308436149
- 36. Motha MX, Egerton JR (1987) Outbreak of atypical fowlpox in chickens with persistent reticuloendotheliosis viraemia. Avian Pathology: Journal of the WVPA 16: 177–182. doi: 10.1080/03079458708436362
- 37. Drew ML, Wigle WL, Graham DL, Griffin CP, Silvy NJ, et al. (1998) Reticuloendotheliosis in captive greater and Attwaters prairie chickens. J Wildl Dis 34: 783–791. doi: 10.7589/0090-3558-34.4.783
- 38. Trampel DW, Pepper TM, Witter RL (2002) Reticuloendotheliosis in Hungarian partridge. Journal of Wildlife Diseases 38: 438–442. doi: 10.7589/0090-3558-38.2.438
- 39. Wang Y, Cui Z, Jiang S (2006) Sequence analysis for the complete proviral genome of reticuloendotheliosis virus Chinese strain HA9901. Sci China C Life Sci 49: 149–157. doi: 10.1007/s11427-006-0149-8
- 40. Franklin RB, Maldonado RL, Bose HR (1974) Isolation and characterization of reticuloendotheliosis virus transformed bone marrow cells. Intervirology 3: 342–352. doi: 10.1159/000149771
- 41. Burmester BR, Purchase HG (1979) The history of avian medicine in the United States. V. Insights into avian tumor virus research. Avian Dis 23: 1–29. doi: 10.2307/1589667
- 42. Hughes AL, Irausquin S, Friedman R (2010) The evolutionary biology of poxviruses. Infect Genet Evol 10: 50–59. doi: 10.1016/j.meegid.2009.10.001
- 43. McGeoch DJ, Gatherer D (2005) Integrating reptilian herpesviruses into the family herpesviridae. J Virol 79: 725–731. doi: 10.1128/jvi.79.2.725-731.2005
- 44. Parrish CR, Holmes EC, Morens DM, Park EC, Burke DS, et al. (2008) Cross-species virus transmission and the emergence of new epidemic diseases. Microbiol Mol Biol Rev 72: 457–470. doi: 10.1128/mmbr.00004-08
- 45. Brumpt É (1935) Paludisme aviaire: Plasmodium gallinaceum n. sp. de la poule domestique. C R Hebd Seances Acad Sci 200: 783–785.
- 46. Vanderberg JP (2009) Reflections on early malaria vaccine studies, the first successful human malaria vaccination, and beyond. Vaccine 27: 2–9. doi: 10.1016/j.vaccine.2008.10.028
- 47. Conway WG (1972) In Memoriam: Lee Saunders Crandall. Auk 89: 420–248. doi: 10.2307/4084215
- 48. Coggeshall LT (1938) Plasmodium lophurae, a new species of malaria pathogenic for the domestic fowl. Am J Hyg 27: 615–618.
- 49. Terzian LA (1941) Studies on Plasmodium lophurae, a malarial parasite in fowls. American Journal of Epidemiology 33: 1–22.
- 50. Hitchner SB (1996) Poultry vaccine laboratories in the United States: a historical perspective. Avian Dis 40: 255–265. doi: 10.2307/1592218
- 51. Goldhaft AD (1957) The golden egg. New York: Horizon Press.
- 52. Kim TJ, Tripathy DN (2001) Reticuloendotheliosis virus integration in the fowl poxvirus genome: not a recent event. Avian diseases 45: 663–669. doi: 10.2307/1592909
- 53. Witter RL, Sharma JM, Fadly AM (1986) Nonbursal lymphomas induced by nondefective reticuloendotheliosis virus. Avian Pathology: Journal of the WVPA 15: 467–486. doi: 10.1080/03079458608436308
- 54. Isfort R, Jones D, Kost R, Witter R, Kung HJ (1992) Retrovirus insertion into herpesvirus in vitro and in vivo. Proc Natl Acad Sci U S A 89: 991–995. doi: 10.1073/pnas.89.3.991
- 55. Kang CY, Temin HM (1974) Reticuloendotheliosis virus nucleic acid sequences in cellular DNA. J Virol 14: 1179–1188.
- 56. Franklin RB, Kang CY, Min-Min Wan K, Bose HR Jr (1977) Transformation of chick embryo fibroblasts by reticuloendotheliosis virus. Virology 83: 313–321. doi: 10.1016/0042-6822(77)90176-3
- 57. Swofford DL (1998) PAUP*. Phylogenetic analysis using parsimony (* and other methods). 4 ed. Sunderland, MA: Sinauer Associates.
- 58. Abascal F, Zardoya R, Posada D (2005) ProtTest: selection of best-fit models of protein evolution. Bioinformatics 21: 2104–2105. doi: 10.1093/bioinformatics/bti263
- 59. Dimmic MW, Rest JS, Mindell DP, Goldstein RA (2002) rtREV: an amino acid substitution matrix for inference of retrovirus and reverse transcriptase phylogeny. J Mol Evol 55: 65–73. doi: 10.1007/s00239-001-2304-y
- 60. Nickle DC, Heath L, Jensen MA, Gilbert PB, Mullins JI, et al. (2007) HIV-specific probabilistic models of protein evolution. PLoS One 2: e503 doi:10.1371/journal.pone.0000503.
- 61. Kosakovsky Pond SL, Posada D, Gravenor MB, Woelk CH, Frost SD (2006) GARD: a genetic algorithm for recombination detection. Bioinformatics 22: 3096–3098. doi: 10.1093/bioinformatics/btl474
- 62. Huson DH (1998) SplitsTree: analyzing and visualizing evolutionary data. Bioinformatics 14: 68–73. doi: 10.1093/bioinformatics/14.1.68
- 63. Rodriguez F, Oliver JL, Marin A, Medina JR (1990) The general stochastic model of nucleotide substitution. J Theor Biol 142: 485–501. doi: 10.1016/s0022-5193(05)80104-3
- 64. Stamatakis A (2006) RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics 22: 2688–2690. doi: 10.1093/bioinformatics/btl446
- 65. Yang Z (1997) PAML: a program package for phylogenetic analysis by maximum likelihood. Comput Appl Biosci 13: 555–556. doi: 10.1093/bioinformatics/13.5.555
- 66. Corwin RM, McGhee RB (1966) Anemia in ducklings treated with malarious plasma. Experimental Parasitology 18: 281–289. doi: 10.1016/0014-4894(66)90027-0
- 67. Cook MK (1969) Cultivation of a filterable agent associated with Mareks disease. J Natl Cancer Inst 43: 203–212.
- 68. Carlson HC, Seawright GL, Pettit JR (1974) Reticuloendotheliosis in Japanese quail. Avian Pathol 3: 169–175. doi: 10.1080/03079457408418111
- 69. Grimes TM, Purchase HG (1973) Reticuloendotheliosis in a duck. Aust Vet J 49: 466–471. doi: 10.1111/j.1751-0813.1973.tb09295.x
- 70. Paul PS, Werdin RE, Pomeroy BS (1978) Spontaneously occurring lymphoproliferative disease in ducks. Avian Diseases 22: 191–195. doi: 10.2307/1589526
- 71. Schat KA, Gonzalez J, Solorzano A, Avila E, Witter RL (1976) A lymphoproliferative disease in Japanese Quail. Avian Diseases 20: 153–161. doi: 10.2307/1589484
- 72. Sarma PS, Jain DK, Mishra NK, Vernon ML, Paul PS, et al. (1975) Isolation and characterization of viruses from natural outbreaks of reticuloendotheliosis in turkeys. J Natl Cancer Inst 54: 1355–1359.
- 73. Koyama H, Inoue K, Nagashima T, Ohwada Y, Saito Y (1975) Cause of “nakanuke” in chickens. I. Occurrence of “nakanuke” in chicken inoculated with the cells showed coexistence or C-type virus and turkey herpesvirus. The Kitasato Archives of Experimental Medicine 48: 83–91.
- 74. McDougall JS, Biggs PM, Shilleto RW (1978) A leukosis in turkeys associated with infection with reticuloendotheliosis virus. Avian Pathology: Journal of the WVPA 7: 557–568. doi: 10.1080/03079457808418310
- 75. Solomon JJ, Witter RL, Nazerian K (1976) Studies on the etiology of lymphomas in turkeys: isolation of reticuloendotheliosis in virus. Avian Diseases 20: 735–747. doi: 10.2307/1589454
- 76. Ratnamohan N, Grimes TM, Bagust, Spradbrow PB (1980) A transmissible chicken tumour associated with reticuloendotheliosis virus infection. Aust Vet J 56: 34–38. doi: 10.1111/j.1751-0813.1980.tb02538.x
- 77. Witter RL, Johnson DC (1985) Epidemiology of reticuloendotheliosis virus in broiler breeder flocks. Avian Diseases 29: 1140–1154. doi: 10.2307/1590468
- 78. Dren CN, Nemeth I, Gyorvari I, Horvath P (1988) Isolation and characterization of reticuloendotheliosis virus from turkeys. Acta Vet Hung 36: 265–282.
- 79. Dren CN, Nemeth I, Sari I, Ratz F, Glavits R, et al. (1988) Isolation of a reticuloendotheliosis-like virus from naturally occurring lymphoreticular tumours of domestic goose. Avian Pathol 17: 259–277. doi: 10.1080/03079458808436446
- 80. Ley DH, Ficken MD, Cobb DT, Witter RL (1989) Histomoniasis and reticuloendotheliosis in a wild turkey (Meleagris gallopavo) in North Carolina. J Wildl Dis 25: 262–265. doi: 10.7589/0090-3558-25.2.262
- 81. Hayes LE, Langheinrich KA, Witter RL (1992) Reticuloendotheliosis in a wild turkey (Meleagris gallopavo) from coastal Georgia. J Wildl Dis 28: 154–158. doi: 10.7589/0090-3558-28.1.154
- 82. Fadly AM, Witter RL, Smith EJ, Silva RF, Reed WM, et al. (1996) An outbreak of lymphomas in commercial broiler breeder chickens vaccinated with a fowlpox vaccine contaminated with reticuloendotheliosis virus. Avian Pathology: Journal of the WVPA 25: 35–47. doi: 10.1080/03079459608419118
- 83. Miller PE, Paul-Murphy J, Sullivan R, Cooley AJ, Dubielzig RR, et al. (1998) Orbital lymphosarcoma associated with reticuloendotheliosis virus in a peafowl. J Am Vet Med Assoc 213: 377–380.
- 84. Davidson I, Borenstein R (1999) Multiple infection of chickens and turkeys with avian oncogenic viruses: prevalence and molecular analysis. Acta Virologica 43: 136–142.
- 85. Peterson MJ, Aguirre R, Ferro PJ, Jones DA, Lawyer TA, et al. (2002) Infectious disease survey of Rio Grande wild turkeys in the Edwards Plateau of Texas. J Wildl Dis 38: 826–833. doi: 10.7589/0090-3558-38.4.826
- 86. Crespo R, Woolcock PR, Fadly AM, Hall C, Shivaprasad HL (2002) Characterization of T-cell lymphomas associated with an outbreak of reticuloendotheliosis in turkeys. Avian Pathol 31: 355–361. doi: 10.1080/03079450220141624
- 87. Garcia M, Narang N, Reed WM, Fadly AM (2003) Molecular characterization of reticuloendotheliosis virus insertions in the genome of field and vaccine strains of fowl poxvirus. Avian Diseases 47: 343–354. doi: 10.1637/0005-2086(2003)047[0343:mcorvi]2.0.co;2
- 88. Cheng WH, Huang YP, Wang CH (2006) Serological and virological surveys of reticuloendotheliosis in chickens in Taiwan. The Journal of Veterinary Medical Science/The Japanese Society of Veterinary Science 68: 1315–1320. doi: 10.1292/jvms.68.1315
- 89. Barbosa T, Zavala G, Cheng S, Villegas P (2007) Full genome sequence and some biological properties of reticuloendotheliosis virus strain APC-566 isolated from endangered Attwaters prairie chickens. Virus Research 124: 68–77. doi: 10.1016/j.virusres.2006.10.002
- 90. Prukner-Radovcic E, Luschow D, Ciglar Grozdanic I, Tisljar M, Mazija H, et al. (2006) Isolation and molecular biological investigations of avian poxviruses from chickens, a turkey, and a pigeon in Croatia. Avian Diseases 50: 440–444. doi: 10.1637/7506-012006r.1
- 91. Yu LJ, Cui ZZ (2006) [Sequence analysis of integration sites of reticuloendotheliosis virus LTR in fowlpox vaccine virus genomes]. Wei sheng wu xue bao = Acta microbiologica Sinica 46: 660–662.
- 92. Tadese T, Fitzgerald S, Reed WM (2008) Detection and differentiation of re-emerging fowlpox virus (FWPV) strains carrying integrated reticuloendotheliosis virus (FWPV-REV) by real-time PCR. Veterinary Microbiology 127: 39–49. doi: 10.1016/j.vetmic.2007.08.012
- 93. Cheng Z, Shi Y, Zhang L, Zhu G, Diao X, et al. (2007) Occurrence of reticuloendotheliosis in Chinese partridge. The Journal of Veterinary Medical Science/The Japanese Society of Veterinary Science 69: 1295–1298. doi: 10.1292/jvms.69.1295
- 94. Liu Q, Zhao J, Su J, Pu J, Zhang G, et al. (2009) Full genome sequences of two reticuloendotheliosis viruses contaminating commercial vaccines. Avian Dis 53: 341–346. doi: 10.1637/8579-010609-reg.1
- 95. Lin CY, Chen CL, Wang CC, Wang CH (2009) Isolation, identification, and complete genome sequence of an avian reticuloendotheliosis virus isolated from geese. Veterinary Microbiology 136: 246–249. doi: 10.1016/j.vetmic.2008.11.017
- 96. Awad AM, Abd El-Hamid HS, Abou Rawash AA, Ibrahim HH (2010) Detection of reticuloendotheliosis virus as a contaminant of fowl pox vaccines. Poult Sci 89: 2389–2395. doi: 10.3382/ps.2010-00899
- 97. Arathy DS, Tripathy DN, Sabarinath GP, Bhaiyat MI, Chikweto A, et al. (2010) Preliminary molecular characterization of a fowl poxvirus isolate in Grenada. Avian Diseases 54: 1081–1085. doi: 10.1637/9015-081709-resnote.1
- 98. Diallo IS, Taylor J, Gibson J, Hoad J, De Jong A, et al. (2010) Diagnosis of a naturally occurring dual infection of layer chickens with fowlpox virus and gallid herpesvirus 1 (infectious laryngotracheitis virus). Avian Pathology: Journal of the WVPA 39: 25–30. doi: 10.1080/03079450903447412