Molecular evolution of poxviruses

Previous restriction fragment length polymorphism analysis divided variola virus (VARV) strains into two subtypes, one of which included West African and South American isolates. This allowed a dating to be introduced for the first time in estimation of the VARV evolution rate. The results were used to analyze the molecular evolution of the total family Poxviridae. Comparisons of the known nucleotide sequences were performed for the extended conserved central genome region in 42 orthopoxvirus strains and for the eight genes of multisubunit RNA polymerase in 65 viruses belonging to various genera of the family Poxviridae. Using the Bayesian dating method, the mutation accumulation rate of poxviruses was estimated at (1.7–8.8) × 10^?6 nucleotide substitutions per site per year. Computations showed that the modern poxvirus genera started diverging from an ancestral virus more than 200 thousand years ago and that an ancestor of the genus Orthopoxvirus emerged 131 ± 45 thousand years ago. The other genera of mammalian poxviruses with a low GC content diverged approximately 110–90 thousand years ago. The independent evolution of VARV started 3.4 ± 0.8 thousand years ago. It was shown with the example of VARV and the monkeypox virus (MPXV) that divergent evolution of these orthopoxviruses started and the West African subtypes of VARV and MPXV were formed as geographical conditions changed to allow isolation of West African animals from other African regions.     

The time scale in poxvirus evolution

Unlike vertebrates and RNA-containing viruses, the objective estimate of molecular clock for DNA-containing viruses was so far absent. An extended central conservative genomic region of orthopoxviruses (about 102 kbp) and the sequence of DNA polymerase gene (about 3 kbp) of the viruses belonging to various genera from the family Poxviridae were analyzed. During this analysis, the known dating of variola virus (VARV) transfer from West Africa to South America (XVI century) and our own data on close phylogenetic relations between the modem West African and South American VARV isolates were used. As a result of this work, it was calculated for the first time that the rate of mutation accumulation in these DNA-containing viruses amounted to 0.9-1.2 x 10(-6) substitutions per site per year. The poxviruses started separating from the ancestor virus to form the modem genera approximately 500 thousand years ago; the ancestor of the genus Orthopoxvirus separated about 300 thousand years ago; and its division into the modem studied species took place approximately 14 thousand years ago.     

Time scale of Poxvirus evolution

Unlike in vertebrates and RNA viruses, the molecular clock has not been estimated so far for DNA viruses. The extended conserved central region (102 kb) of the orthopoxvirus genome and the DNA polymerase gene (3 kb) were analyzed in viruses representing several genera of the family Poxviridae. Analysis was based on the known dating of the variola virus (VARV) transfer from Western Africa to South America and previous data on the phylogenetic relatedness of modern West African and South American isolates of VARV. The mutation accumulation rate was for the first time estimated for these DNA viruses at (0.9–1.2) × 10^?6 substitutions per site per year. It was assumed that poxviruses diverged from an ancestor approximately 500,000 years ago to form the recent species and that the ancestor of the genus Orthopoxvirus emerged approximately 300,000 years ago and gave origin to the modern species approximately 14,000 years ago.     

Pandemic potential of poxviruses: From an ancient killer causing smallpox to the surge of monkeypox

Smallpox caused by the variola virus (VARV) was one of the greatest infectious killers of mankind. Historical records trace back smallpox for at least a millennium while phylogenetic analysis dated the ancestor of VARV circulating in the 20th century into the 19th century. The discrepancy was solved by the detection of distinct VARV sequences first in 17th-century mummies and then in human skeletons dated to the 7th century. The historical records noted marked variability in VARV virulence which scientists tentatively associated with gene losses occurring when broad-host poxviruses narrow their host range to a single host. VARV split from camel and gerbil poxviruses and had no animal reservoir, a prerequisite for its eradication led by WHO. The search for residual pockets of VARV led to the discovery of the monkeypox virus (MPXV); followed by the detection of endemic smallpox-like monkeypox (mpox) disease in Africa. Mpox is caused by less virulent clade 2 MPXV in West Africa and more virulent clade 1 MPXV in Central Africa. Exported clade 2 mpox cases associated with the pet animal trade were observed in 2003 in the USA. In 2022 a world-wide mpox epidemic infecting more than 80,000 people was noted, peaking in August 2022 although waning rapidly. The cases displayed particular epidemiological characteristics affecting nearly exclusively young men having sex with men (MSM). In contrast, mpox in Africa mostly affects children by non-sexual transmission routes possibly from uncharacterized animal reservoirs. While African children show a classical smallpox picture, MSM mpox cases show few mostly anogenital lesions, low-hospitalization rates and 140 fatal cases worldwide. MPXV strains from North America and Europe are closely related, derived from clade 2 African MPXV. Distinct transmission mechanisms are more likely causes for the epidemiological and clinical differences between endemic African cases and the 2022 epidemic cases than viral traits.     

Properties of the recombinant TNF-binding proteins from variola, monkeypox, and cowpox viruses are different

Tumor necrosis factor (TNF), a potent proinflammatory and antiviral cytokine, is a critical extracellular immune regulator targeted by poxviruses through the activity of virus-encoded family of TNF-binding proteins (CrmB, CrmC, CrmD, and CrmE). The only TNF-binding protein from variola virus (VARV), the causative agent of smallpox, infecting exclusively humans, is CrmB. Here we have aligned the amino acid sequences of CrmB proteins from 10 VARV, 14 cowpox virus (CPXV), and 22 monkeypox virus (MPXV) strains. Sequence analyses demonstrated a high homology of these proteins. The regions homologous to cd00185 domain of the TNF receptor family, determining the specificity of ligand-receptor binding, were found in the sequences of CrmB proteins. In addition, a comparative analysis of the C-terminal SECRET domain sequences of CrmB proteins was performed. The differences in the amino acid sequences of these domains characteristic of each particular orthopoxvirus species were detected. It was assumed that the species-specific distinctions between the CrmB proteins might underlie the differences in these physicochemical and biological properties. The individual recombinant proteins VARV-CrmB, MPXV-CrmB, and CPXV-CrmB were synthesized in a baculovirus expression system in insect cells and isolated. Purified VARV-CrmB was detectable as a dimer with a molecular weight of 90 kDa, while MPXV- and CPXV-CrmBs, as monomers when fractioned by non-reducing SDS-PAGE. The CrmB proteins of VARV, MPXV, and CPXV differed in the efficiencies of inhibition of the cytotoxic effects of human, mouse, or rabbit TNFs in L929 mouse fibroblast cell line. Testing of CrmBs in the experimental model of LPS-induced shock using SPF BALB/c mice detected a pronounced protective effect of VARV-CrmB. Thus, our data demonstrated the difference in anti-TNF activities of VARV-, MPXV-, and CPXV-CrmBs and efficiency of VARV-CrmB rather than CPXV- or MPXV-CrmBs against LPS-induced mortality in mice.     

Dating the origin of the African human T-cell lymphotropic virus type-i (HTLV-I) subtypes

To investigate the origin of the African PTLV-I virus, we phylogenetically analyzed the available HTLV-I and STLV-I strains. We also attempted to date the presumed interspecies transmissions that resulted in the African HTLV-I subtypes. Molecular-clock analysis was performed using the Tamura-Nei substitution model and gamma distributed rate heterogeneity based on the maximum-likelihood topology of the combined long-terminal-repeat and env third-codon-position sequences. Since the molecular clock was not rejected and no evidence for saturation was found, a constant rate of evolution at these positions for all 33 HTLV-I and STLV-I strains was reasonably assumed. The spread of PTLV-I in Africa is estimated to have occurred at least 27,300 +/- 8,200 years ago. Using the available strains, the HTLV-If subtype appears to have emerged within the last 3,000 years, and the HTLV-Ia, HTLV-Ib, HTLV-Id, and HTLV-Ie subtypes appear to have diverged between 21,100 and 5,300 years ago. Interspecies transmissions, most probably simian to human, must have occurred around that time and probably continued later. When the synonymous and nonsynonymous substitution ratios were compared, it was clear that purifying selection was the driving force for PTLV-I evolution in the env gene, irrespective of the host species. Due to the small number of strains in some of the investigated groups, these data on selective pressure should be taken with caution.     

Molecular evolution of the tick-borne encephalitis and Powassan viruses

Issues associated with newly emerging viruses, their genetic diversity, and viral evolution in modern environments are currently attracting growing attention. In this study, a phylogenetic analysis was performed and the evolution rate was evaluated for such pathogenic flaviviruses endemic to Russia as tick-borne encephalitis virus (TBEV) and Powassan virus (PV). The analysis involved 47 nucleotide sequences of the TBEV genome region encoding protein E and 17 sequences of the PV NS5-encoding region. The nucleotide substitution rate was estimated as 1.4 × 10⁻⁴ and 5.4 × 10⁻⁵ substitutions per site per year for the E protein-encoding region of the TBEV genome and for the NS5 genome region of PV, respectively. The ratio of non-synonymous to synonymous nucleotide substitutions (dN/dS) in viral sequences was calculated as 0.049 for TBEV and 0.098 for PV. The highest dN/dS values of 0.201–0.220 were found in the subcluster of Russian and Canadian PV strains, and the lowest value of 0.024 was observed in the cluster of Russian and Chinese strains of the Far Eastern TBEV genotype. Evaluation of time intervals between the events of viral evolution showed that the European subtype of TBEV diverged from the common TBEV ancestor approximately 2750 years ago, while the Siberian and Far Eastern subtypes emerged approximately 2250 years ago. The PV was introduced into its natural foci of the Russian Primorskii krai only approximately 70 years ago; these strains were very close to Canadian PV strains. The pattern of PV evolution in North America was similar to the evolution of the Siberian and Far Eastern TBEV subtypes in Asia. The moments of divergence between major genetic groups of TBEV and PV coincide with historical periods of climate warming and cooling, suggesting that climate change was a key factor in the evolution of flaviviruses in past millennia.     

Molecular evolution of the tick-borne encephalitis and Powassan viruses

The problem of emerging viruses, their genetic diversity and viral evolution in nature are attracting more attention. The phylogenetic analysis and evaluationary rate estimation were made for pathogenic flaviviruses such as tick-borne encephalitis virus (TBEV) and Powassan (PV) circulated in natural foci in Russia. 47 nucleotide sequences of encoded protein E of the TBEV and 17 sequences of NS5 genome region of the PV have been used. It was found that the rate of accumulation of nucleotide substitutions for E genome region of TBEV was approximately 1.4 x 10(-4) and 5.4 x 10(-5) substitutions per site per year for NS5 genome region of PV. The ratio of non-synonymous nucleotide substitutions to synonymous substitution (dN/dS) for viral sequences were estimated of 0.049 for TBEV and 0.098 for PV. Maximum value dN/dS was 0.201-0.220 for sub-cluster of Russian and Canadian strains of PV and the minimum - 0.024 for cluster of Russian and Chinese strains of Far Eastern genotype TBEV. Evaluation of time intervals of evolutionary events associated with these viruses showed that European subtype TBEV are diverged from all-TBEV ancestor within approximately 2750 years and the Siberian and Far Eastern subtypes are emerged about 2250 years ago. The PV was introduced into natural foci of the Primorsky Krai of Russia only about 70 years ago and PV is a very close to Canadian strains of PV. Evolutionary picture for PV in North America is similar to evolution of Siberian and Far Eastern subtypes TBEV in Asia. The divergence time for main genetic groups of TBEV and PV are correlated with historical periods of warming and cooling. These allow to propose a hypothesis that climate changes were essential to the evolution of the flaviviruses in the past millenniums.     

Correlates of protective immunity in poxvirus infection: where does antibody stand?

Even though smallpox has been eradicated, the threat of accidental or intentional release has highlighted the fact there is little consensus about correlates of protective immunity or immunity against re-infection with the causative poxvirus, variola virus (VARV). As the existing vaccine for smallpox has unacceptable rates of side effects and complications, new vaccines are urgently needed. Surrogate animal models of VARV infection in humans, including vaccinia virus (VACV) and ectromelia virus (ECTV) infection in mice, monkeypox virus (MPXV) infection in macaques have been used as tools to dissect the immune response to poxviruses. Mousepox, caused by ECTV, a natural mouse pathogen, is arguably the best surrogate small-animal model, as it shares many aspects of virus biology, pathology and clinical features with smallpox in humans. The requirements for recovery from a primary ECTV infection have been well characterized and include type I and II interferons, natural killer cells, CD4T cells, CD8T cell effector function and antibody. From a vaccine standpoint, it is imperative that the requirements for recovery from secondary infection are also identified. We have investigated host immune parameters in response to a secondary ECTV infection, and have identified that interferon and CD8T cell effector functions are not essential; however, T- and B-cell interaction and antibody are absolutely critical for recovery from a secondary challenge. The central role of antibody has been also been identified in the secondary response to other poxviruses. These findings have important clinical implications and would greatly assist the design of therapeutic interventions and new vaccines for smallpox.     

Origin and evolution of influenza virus hemagglutinin genes

Influenza A, B, and C viruses are the etiological agents of influenza. Hemagglutinin (HA) is the major envelope glycoprotein of influenza A and B viruses, and hemagglutinin-esterase (HE) in influenza C viruses is a protein homologous to HA. Because influenza A virus pandemics in humans appear to occur when new subtypes of HA genes are introduced from aquatic birds that are known to be the natural reservoir of the viruses, an understanding of the origin and evolution of HA genes is of particular importance. We therefore conducted a phylogenetic analysis of HA and HE genes and showed that the influenza A and B virus HA genes diverged much earlier than the divergence between different subtypes of influenza A virus HA genes. The rate of amino acid substitution for A virus HAs from duck, a natural reservoir, was estimated to be 3.19 x 10(-4) per site per year, which was slower than that for human and swine A virus HAs but similar to that for influenza B and C virus HAs (HEs). Using this substitution rate from the duck, we estimated that the divergences between different subtypes of A virus HA genes occurred from several thousand to several hundred years ago. In particular, the earliest divergence time was estimated to be about 2,000 years ago. Also, the A virus HA gene diverged from the B virus HA gene about 4,000 years ago and from the C virus HE gene about 8,000 years ago. These time estimates are much earlier than the previous ones.     

Identification of Wild-Derived Inbred Mouse Strains Highly Susceptible to Monkeypox Virus Infection for Use as Small Animal Models

Infection with monkeypox virus (MPXV) causes disease manifestations in humans that are similar, although usually less severe, than those of smallpox. Since routine vaccination for smallpox ceased more than 30 years ago, there is concern that MPXV could be used for bioterrorism. Thus, there is a need to develop animal models to study MPXV infection. Accordingly, we screened 38 inbred mouse strains for susceptibility to MPXV. Three highly susceptible wild-derived inbred strains were identified, of which CAST/EiJ was further developed as a model. Using an intranasal route of infection with an isolate of the Congo Basin clade of MPXV, CAST/EiJ mice exhibited weight loss, morbidity, and death in a dose-dependent manner with a calculated 50% lethal dose (LD₅₀) of 680 PFU, whereas there were no deaths of BALB/c mice at a 10,000-fold higher dose. CAST/EiJ mice exhibited greater MPXV sensitivity when infected via the intraperitoneal route, with an LD₅₀ of 14 PFU. Both routes resulted in MPXV replication in the lung, spleen, and liver. Intranasal infection with an isolate of the less-pathogenic West African clade yielded an LD₅₀ of 7,600 PFU. The immune competence of CAST/EiJ mice was established by immunization with vaccinia virus, which induced antigen-specific T- and B-lymphocyte responses and fully protected mice from lethal doses of MPXV. The new mouse model has the following advantages for studying pathogenesis of MPXV, as well as for evaluation of potential vaccines and therapeutics: relative sensitivity to MPXV through multiple routes, genetic homogeneity, available immunological reagents, and commercial production.Monkeypox virus (MPXV), a member of the orthopoxvirus genus of the Chordopoxvirinae subfamily of the Poxviridae, was isolated in 1958 from lesions in a cynomolgous monkey that had been imported from Africa (27). The first human infections with MPXV were reported in 1972, and since then more than two thousand cases have been recorded, most in the Democratic Republic of the Congo and lesser numbers in West African countries (reviewed by Parker et al. [18]). The mortality from human monkeypox in the Congo is estimated to be 10% of infected individuals with clinical symptoms that mimic smallpox, which is caused by another member of the orthopoxvirus genus: variola virus. However, whereas the host range of variola virus is restricted to humans, serological studies indicate that MPXV naturally infects a large number of animal species, particularly squirrels and nonhuman primates. The sporadic occurrence of human monkeypox is thought to arise from close proximity and handling of infected animals. In this respect, a self-limited outbreak in the United States was traced to a shipment of West African rodents (19). Although monkeypox is a minor public health problem when compared historically to smallpox, the potential for expansion of the MPXV host range and adaptations to enhance human transmission make it prudent to continue careful surveillance. Moreover, the potential use of MPXV for bioterrorism has led to its inclusion as a select agent in the United States (http://www.selectagents.gov).Animal models are crucial for studying virus pathogenesis, and MXPV is no exception. Ground squirrels (22, 26), black-tailed prairie dogs (9, 11, 13, 30), and African dormice (23) are highly susceptible to MPXV. However, as experimental systems, each has limitations with regard to unavailability of commercial breeding, genetic heterogeneity and absence of immunological and other reagents. Laboratory mice, including BALB/c, C57BL/6, and several other mouse strains tested, were found to be resistant to MPXV disease unless impaired in innate or acquired immunity (10, 17, 24). In the present study, we tested a large group of distinct inbred strains of mice chosen for genetic diversity, inclusion of classical and wild-derived strains, and commercial availability. Of 38 inbred mouse strains tested, three wild-derived strains were highly susceptible to MPXV. One of these, CAST/EiJ, was further characterized with regard to MPXV strain sensitivity, route of inoculation, virus dissemination, immune response, and protection by vaccination and drug treatment.     

Archaic admixture in the human genome

One of the enduring questions in the evolution of our species surrounds the fate of 'archaic' forms of Homo. Did Neanderthals go extinct without interbreeding with modern humans 25-40 thousand years ago or are their genes present among modern-day Europeans? Recent work suggests that Neanderthals and an as yet unidentified archaic African population contributed to at least 5% of the modern European and West African gene pools, respectively. Extensive sequencing of Neanderthal and other archaic human nuclear DNA has the potential to answer this question definitively within the next few years.     

The genome of Yoka poxvirus

Yoka poxvirus was isolated almost four decades ago from a mosquito pool in the Central African Republic. Its classification as a poxvirus is based solely upon the morphology of virions visualized by electron microscopy. Here we describe sequencing of the Yoka poxvirus genome using a combination of Roche/454 and Illumina next-generation sequencing technologies. A single consensus contig of ～175 kb in length that encodes 186 predicted genes was generated. Multiple methods were used to show that Yoka poxvirus is most closely related to viruses in the Orthopoxvirus genus, but it is clearly distinct from previously described poxviruses. Collectively, the phylogenetic and genomic sequence analyses suggest that Yoka poxvirus is the prototype member of a new genus in the family Poxviridae.     

Comparative restriction enzyme analysis of the genome in variola virus strains from the Russian collection

Comparative RFLP analysis was for the first time performed for 21 variola virus (VARV) strains of the Russian collection with 20 amplicons covering the total VARV genome. The amplicons were synthesized in the long polymerase chain reaction. A database useful as a reference for identifying VARV strains was generated. VARV strains isolated in different geographical regions were compared and proved to vary mostly in variable genome regions. Each of the dendrograms constructed included three clusters of African, Asian, and VARV-alastrim isolates. The VARV-alastrim isolates differed to the greatest extent from the other strains. VARV strains isolated during an ecdemic variola burst in Moscow (1960) grouped with Asian isolates. Polymorphism of VARV strains was for the first time observed for a single variola burst with a few affected patients.     

Comparative Restriction Enzyme Analysis of the Genome in Variola Virus Strains from the Russian Collection

Comparative RFLP analysis was for the first time performed for 21 variola virus (VARV) strains of the Russian collection with 20 amplicons covering the total VARV genome. The amplicons were synthesized in the long polymerase chain reaction. A database useful as a reference for identifying VARV strains was generated. VARV strains isolated in different geographical regions were compared and proved to vary mostly in variable genome regions. Each of the dendrograms constructed included three clusters of African, Asian, and VARV-alastrim isolates. The VARV-alastrim isolates differed to the greatest extent from the other strains. VARV strains isolated during an ecdemic variola burst in Moscow (1960) grouped with Asian isolates. Polymorphism of VARV strains was for the first time observed for a single variola burst with a few affected patients.     

Poxvirus genome evolution by gene gain and loss

The poxviruses (Poxviridae) are a family of viruses with double-stranded DNA genomes and substantial numbers (often >200) of genes per genome. We studied the patterns of gene gain and loss over the evolutionary history of 17 poxvirus complete genomes. A phylogeny based on gene family presence/absence showed good agreement with families based on concatenated amino acid sequences of conserved single-copy genes. Gene duplications in poxviruses were often lineage specific, and the most extensively duplicated viral gene families were found in only a few of the genomes analyzed. A total of 34 gene families were found to include a member in at least one of the poxvirus genomes analyzed and at least one animal genome; in 16 (47%) of these families, there was evidence of recent horizontal gene transfer (HGT) from host to virus. Gene families with evidence of HGT included several involved in host immune defense mechanisms (the MHC class I, interleukin-10, interleukin-24, interleukin-18, the interferon gamma receptor, and tumor necrosis factor receptor II) and others (glutaredoxin and glutathione peroxidase) involved in resistance of cells to oxidative stress. Thus "capture" of host genes by HGT has been a recurrent feature of poxvirus evolution and has played an important role in adapting the virus to survive host antiviral defense mechanisms.     

The origin and evolution of Ebola and Marburg viruses

Molecular evolutionary analyses for Ebola and Marburg viruses were conducted with the aim of elucidating evolutionary features of these viruses. In particular, the rate of nonsynonymous substitutions for the glycoprotein gene of Ebola virus was estimated to be, on the average, 3.6 x 10(-5) per site per year. Marburg virus was also suggested to be evolving at a similar rate. Those rates were a hundred times slower than those of retroviruses and human influenza A virus, but were of the same order of magnitude as that of the hepatitis B virus. When these rates were applied to the degree of sequence divergence, the divergence time between Ebola and Marburg viruses was estimated to be more than several thousand years ago. Moreover, most of the nucleotide substitutions were transitions and synonymous for Marburg virus. This suggests that purifying selection has operated on Marburg virus during evolution.      

Transmissibility of the Monkeypox Virus Clades via Respiratory Transmission: Investigation Using the Prairie Dog-Monkeypox Virus Challenge System

Monkeypox virus (MPXV) is endemic within Africa where it sporadically is reported to cause outbreaks of human disease. In 2003, an outbreak of human MPXV occurred in the US after the importation of infected African rodents. Since the eradication of smallpox (caused by an orthopoxvirus (OPXV) related to MPXV) and cessation of routine smallpox vaccination (with the live OPXV vaccinia), there is an increasing population of people susceptible to OPXV diseases. Previous studies have shown that the prairie dog MPXV model is a functional animal model for the study of systemic human OPXV illness. Studies with this model have demonstrated that infected animals are able to transmit the virus to naive animals through multiple routes of exposure causing subsequent infection, but were not able to prove that infected animals could transmit the virus exclusively via the respiratory route. Herein we used the model system to evaluate the hypothesis that the Congo Basin clade of MPXV is more easily transmitted, via respiratory route, than the West African clade. Using a small number of test animals, we show that transmission of viruses from each of the MPXV clade was minimal via respiratory transmission. However, transmissibility of the Congo Basin clade was slightly greater than West African MXPV clade (16.7% and 0% respectively). Based on these findings, respiratory transmission appears to be less efficient than those of previous studies assessing contact as a mechanism of transmission within the prairie dog MPXV animal model.     

Recombinational history and molecular evolution of western equine encephalomyelitis complex alphaviruses.

Western equine encephalomyelitis (WEE) virus (Togaviridae: Alphavirus) was shown previously to have arisen by recombination between eastern equine encephalomyelitis (EEE)- and Sindbis-like viruses (C. S. Hahn, S. Lustig, E. G. Strauss, and J. H. Strauss, Proc. Natl. Acad. Sci. USA 85:5997-6001, 1988). We have now examined the recombinational history and evolution of all viruses belonging to the WEE antigenic complex, including the Buggy Creek, Fort Morgan, Highlands J, Sindbis, Babanki, Ockelbo, Kyzylagach, Whataroa, and Aura viruses, using nucleotide sequences derived from representative strains. Two regions of the genome were examined: sequences of 477 nucleotides from the C terminus of the E1 envelope glycoprotein gene which in WEE virus was derived from the Sindbis-like virus parent, and 517 nucleotide sequences at the C terminus of the nsP4 gene which in WEE virus was derived from the EEE-like virus parent. Trees based on the E1 region indicated that all members of the WEE virus complex comprise a monophyletic group. Most closely related to WEE viruses are other New World members of the complex: the Highlands J, Buggy Creek, and Fort Morgan viruses. More distantly related WEE complex viruses included the Old World Sindbis, Babanki, Ockelbo, Kyzylagach, and Whataroa viruses, as well as the New World Aura virus. Detailed analyses of 38 strains of WEE virus revealed at least 4 major lineages; two were represented by isolates from Argentina, one was from Brazil, and a fourth contained isolates from many locations in South and North America as well as Cuba. Trees based on the nsP4 gene indicated that all New World WEE complex viruses except Aura virus are recombinants derived from EEE- and Sindbis-like virus ancestors. In contrast, the Old World members of the WEE complex, as well as Aura virus, did not appear to have recombinant genomes. Using an evolutionary rate estimate (2.8 x 10(-4) substitutions per nucleotide per year) obtained from E1-3' sequences of WEE viruses, we estimated that the recombination event occurred in the New World 1,300 to 1,900 years ago. This suggests that the alphaviruses originated in the New World a few thousand years ago.