Inferring the rate and time-scale of dengue virus evolution

Future issues in RNA virus evolution

The study of RNA virus evolution has developed rapidly during the past 30 years. This review outlines some important recent findings, as well as a number of the remaining major challenges, particularly those that might explain why RNA viruses are the most important class of emerging diseases, yet often have difficulties adapting to sustained transmission cycles in new hosts. The author emphasizes the relevance of research on the underlying dynamics of mutation, fitness landscapes and the constraints to viral adaptation, as well as the evolution of recombination and reassortment. It is also suggested that a combination of theoretical, experimental and comparative approaches is essential for future studies of viral evolution, coupled with new genome sequence data on intrahost genetic variation.

Invasion and maintenance of dengue virus type 2 and type 4 in the Americas

Dengue virus type 4 (DENV-4) was first reported in the Americas in 1981, where it caused epidemics of dengue fever throughout the region. In the same year, the region's first epidemic of dengue hemorrhagic fever was reported, caused by an Asian strain of dengue virus type 2 (DENV-2) that was distinct from the American subtype circulating previously. Despite the importance of these epidemics, little is known about the rates or determinants of viral spread among island and mainland populations or their directions of movement. We employed a Bayesian coalescent approach to investigate the transmission histories of DENV-2 and DENV-4 since their introduction in 1981 and a parsimony method to assess patterns of strain migration. For both viruses there was an initial invasion phase characterized by an exponential increase in the number of DENV lineages, after which levels of genetic diversity remained constant despite reported fluctuations in DENV-2 and DENV-4 activity. Strikingly, viral lineage numbers increased far more rapidly for DENV-4 than DENV-2, indicative of a more rapid rate of exponential population growth in DENV-4 or a higher rate of geographic dispersal, allowing this virus to move more effectively among localities. We propose that these contrasting dynamics may reflect underlying differences in patterns of host immunity. Despite continued gene flow along particular transmission routes, the overall extent of viral traffic was less than expected under panmixis. Hence, DENV in the Americas has a clear geographic structure that maintains viral diversity between outbreaks.

DNA microarray technique for detection and identification of seven flaviviruses pathogenic for man

A flavivirus microarray was developed for detection and identification of yellow fever (YF), West Nile, Japanese encephalitis (JE), and the dengue 1-4 viruses, which are causing severe human disease all over the world. The microarray was based on 500-nucleotide probe fragments from five different parts of the seven viral genomes. A low-stringent amplification method targeting the corresponding regions of the viral genomic RNA was developed and combined with hybridization to the microarray for detection and identification. For distinction of the generated virus-specific fluorescence-patterns a fitting analysis procedure was adapted. The method was verified as functional for all seven flaviviruses and the strategy for the amplification, combined with the long probes, provided a high tolerance for smaller genetic variability, most suitable for these rapidly changing RNA viruses. A potentially high detection and identification capacity was proven on diverged strains of West Nile and dengue viruses. The lower limit for detection was equivalent, or better, when compared to routinely used RT-PCR methods. The performance of the method was verified on human patient samples containing dengue viruses, or normal human serum spiked with YF or JE viruses. The results demonstrated the ability of the flavivirus microarray to screen simultaneously a sample for several viruses in parallel, in combination with a good lower limit of detection.

Biological transmission of arboviruses: reexamination of and new insights into components, mechanisms, and unique traits as well as their evolutionary trends

Among animal viruses, arboviruses are unique in that they depend on arthropod vectors for transmission. Field research and laboratory investigations related to the three components of this unique mode of transmission, virus, vector, and vertebrate host, have produced an enormous amount of valuable information that may be found in numerous publications. However, despite many reviews on specific viruses, diseases, or interests, a systematic approach to organizing the available information on all facets of biological transmission and then to interpret it in the context of the evolutionary process has not been attempted before. Such an attempt in this review clearly demonstrates tremendous progress made worldwide to characterize the viruses, to comprehend disease transmission and pathogenesis, and to understand the biology of vectors and their role in transmission. The rapid progress in molecular biologic techniques also helped resolve many virologic puzzles and yielded highly valuable data hitherto unavailable, such as characterization of virus receptors, the genetic basis of vertebrate resistance to viral infection, and phylogenetic evidence of the history of host range shifts in arboviruses. However, glaring gaps in knowledge of many critical subjects, such as the mechanism of viral persistence and the existence of vertebrate reservoirs, are still evident. Furthermore, with the accumulated data, new questions were raised, such as evolutionary directions of virus virulence and of host range. Although many fundamental questions on the evolution of this unique mode of transmission remained unresolved in the absence of a fossil record, available observations for arboviruses and the information derived from studies in other fields of the biological sciences suggested convergent evolution as a plausible process. Overall, discussion of the diverse range of theories proposed and observations made by many investigators was found to be highly valuable for sorting out the possible mechanism(s) of the emergence of arboviral diseases.

Strategically examining the full-genome of dengue virus type 3 in clinical isolates reveals its mutation spectra

BACKGROUND

Previous studies presented the quasispecies spectrum of the envelope region of dengue virus type 3 (DENV-3) from either clinical specimens or field-caught mosquitoes. However, the extent of sequence variation among full genomic sequences of DENV within infected individuals remains largely unknown.

RESULTS

Instead of arbitrarily choosing one genomic region in this study, the full genomic consensus sequences of six DENV-3 isolates were used to locate four genomic regions that had a higher potential of sequence heterogeneity at capsid-premembrane (C-prM), envelope (E), nonstructural protein 3 (NS3), and NS5. The extent of sequence heterogeneity revealed by clonal sequencing was genomic region-dependent, whereas the NS3 and NS5 had lower sequence heterogeneity than C-prM and E. Interestingly, the Phylogenetic Analysis by Maximum Likelihood program (PAML) analysis supported that the domain III of E region, the most heterogeneous region analyzed, was under the influence of positive selection.

CONCLUSION

This study confirmed previous reports that the most heterogeneous region of the dengue viral genome resided at the envelope region, of which the domain III was under positive selection pressure. Further studies will need to address the influence of these mutations on the overall fitness in different hosts (i.e., mosquito and human) during dengue viral transmission.

Evolutionary timescale of rabies virus adaptation to North American bats inferred from the substitution rate of the nucleoprotein gene

Throughout North America, rabies virus (RV) is endemic in bats. Distinct RV variants exist that are closely associated with infection of individual host species, such that there is little or no sustained spillover infection away from the primary host. Using Bayesian methodology, nucleotide substitution rates were estimated from alignments of partial nucleoprotein (N) gene sequences of nine distinct bat RV variants from North America. Substitution rates ranged from 2.32 x 10(-4) to 1.38 x 10(-3) substitutions per site per year. A maximum-likelihood (ML) molecular clock model was rejected for only two of the nine datasets. In addition, using sequences from bat RV variants across the Americas, the evolutionary rate for the complete N gene was estimated to be 2.32 x 10(-4). This rate was used to scale trees using Bayesian and ML methods, and the time of the most recent common ancestor for current bat RV variant diversity in the Americas was estimated to be 1660 (range 1267-1782) and 1651 (range 1254-1773), respectively. Our reconstructions suggest that RV variants currently associated with infection of bats from Latin America (Desmodus and Tadarida) share the earliest common ancestor with the progenitor RV. In addition, from the ML tree, times were estimated for the emergence of the three major lineages responsible for bat rabies cases in North America. Adaptation to infection of the colonial bat species analysed (Eptesicus fuscus, Myotis spp.) appears to have occurred much quicker than for the solitary species analysed (Lasionycteris noctivagans, Pipistrellus subflavus, Lasiurus borealis, Lasiurus cinereus), suggesting that the process of virus adaptation may be dependent on host biology.

Overcoming antigenic diversity and improving vaccines using DNA shuffling and screening technologies

Viral, bacterial and parasitic pathogens have evolved multiple strategies to evade the immune response, facilitate transmission and establish chronic infections. One of the underlying strategies that pathogens have evolved is antigenic variation of immune response targets that reduce the affinity of antigen binding to antibodies and major histocompatability complex class I and II receptors. Vaccine candidates generally target a limited number of these antigen variants or combine antigens from several variants to include in multivalent vaccine formulations. DNA shuffling and screening technologies, also known as MolecularBreeding (Maxygen, Inc.) directed molecular evolution, have been successfully used to identify and develop novel and chimaeric vaccine candidates capable of inducing immune responses that recognise and control multiple antigenic variants. DNA shuffling and screening strategies also select vaccine candidates with improved immunogenicity, increased expression as recombinant polypeptides and improved growth of whole viruses in cell culture. As DNA shuffling and screening strategies can be applied to many pathogens, there remain numerous applications of DNA shuffling to solve challenging problems in vaccine process development and manufacture.

Emerging flaviviruses: the spread and resurgence of Japanese encephalitis, West Nile and dengue viruses

Mosquito-borne flaviviruses provide some of the most important examples of emerging and resurging diseases of global significance. Here, we describe three of them: the resurgence of dengue in tropical and subtropical areas of the world, and the spread and establishment of Japanese encephalitis and West Nile viruses in new habitats and environments. These three examples also illustrate the complexity of the various factors that contribute to their emergence, resurgence and spread. Whereas some of these factors are natural, such as bird migration, most are due to human activities, such as changes in land use, water impoundments and transportation, which result in changed epidemiological patterns. The three examples also show the ease with which mosquito-borne viruses can spread to and colonize new areas, and the need for continued international surveillance and improved public health infrastructure to meet future emerging disease threats.

The molecular epidemiology of dengue virus serotype 4 in Bangkok, Thailand

Dengue represents a major public health problem in Thailand, with all four viral serotypes co-circulating. Dengue virus serotype 4 (DENV-4) is the least frequently sampled serotype, although one that is often associated with hemorrhagic fever during secondary infection. To determine the evolutionary forces shaping the genetic diversity of DENV-4, and particularly whether its changing prevalence could be attributed to instances of adaptive evolution in the viral genome, we undertook a large-scale molecular epidemiological analysis of DENV-4 in Bangkok, Thailand, using both E gene and complete coding region sequences. This analysis revealed extensive genetic diversity within a single locality at a single time, including the discovery of a new and divergent genotype of DENV-4, as well as a pattern of continual lineage turnover. We also recorded the highest average rate of evolutionary change for this serotype, at 1.072 x 10(-3) nucleotide substitutions per site, per year. However, despite this abundant genetic variation, there was no evidence for adaptive evolution in any gene, codon, or lineage of DENV-4, with the highest rate of nonsynonymous substitution observed in NS2A. Consequently, the rapid turnover of DENV-4 lineages through time is most likely the consequence of a high rate of deleterious mutation in the viral genome coupled to seasonal fluctuations in the size of the vector population.

Explosive Radiations and the Reliability of Molecular Clocks: Island Endemic Radiations as a Test Case

The reliability of molecular clocks has been questioned for several key evolutionary radiations on the basis that the clock might run fast in explosive radiations. Molecular date estimates for the radiations of metazoan phyla (the Cambrian explosion) and modern orders of mammals and birds are in many cases twice as old as the palaeontological evidence would suggest. Could some aspect of explosive radiations speed the molecular clock, making molecular date estimates too old? Here we use 19 independent instances of recent explosive radiations of island endemic taxa as a model system for testing the proposed influence of rapid adaptive radiation on the rate of molecular evolution. These radiations are often characterized by many of the potential mechanisms for fast rates in explosive radiations--such as small population size, elevated speciation rate, rapid rate of morphological change, release from previous ecological constraints, and adaptation to new niches--and represent a wide variety of species, islands, and genes. However, we find no evidence of a consistent increase in rates in island taxa compared to their mainland relatives, and therefore find no support for the hypothesis that the molecular clock runs fast in explosive radiations.

Sequences of flavivirus-related RNA viruses persist in DNA form integrated in the genome of Aedes spp. mosquitoes

Flavivirus-related sequences have been discovered in the dsDNA genome of Aedes albopictus and Aedes aegypti mosquitoes, demonstrating for the first time an integration into a eukaryotic genome of a multigenic sequence from an RNA virus that replicates without a recognized DNA intermediate. In the Aedes albopictus C6/36 cell line, an open reading frame (ORF) of 1557 aa with protease/helicase and polyprotein processing domains characteristic of flaviviruses was identified. It is closely related to NS1-NS4A genes of the Cell Fusing Agent and Kamiti River virus and the corresponding mRNAs were detected. Integrated sequences homologous to the envelope, NS4B and polymerase genes of flaviviruses were identified. Overall, approximately two-thirds of a flavivirus-like genome were characterized. In the Aedes aegypti A20 cell line, a 492 aa ORF related to the polymerase of the Cell Fusing Agent and Kamiti River virus was identified. These flavivirus-related integrated DNA sequences were detected in laboratory-bred and wild Aedes albopictus and Aedes aegypti mosquitoes, demonstrating that their discovery is not an artefact resulting from the manipulation of mosquito cell lines, since they exist under natural conditions. This finding has major implications regarding evolution, as it represents an entirely different mechanism by which genetic diversity may be generated in eukaryotic cells distinct from accepted processes.

Phylogeography and molecular evolution of dengue 2 in the Caribbean basin, 1981-2000

We sequenced the envelope (E) genes of 59 DEN-2 isolates collected from ten Caribbean islands, six South American countries, and two Central American countries between 1981 and 2000, a period characterized by hyperendemicity and increased incidence of severe dengue. Fifty-two isolates belonged to "American/Asian" subtype IIIb, possessing a characteristic polar residue at envelope aa position 390 (N [n = 48] or S [n = 4]) common to that group. Six isolates from Trinidad (1981), Honduras (1991 [4]), and El Salvador (1987) fell into the "Native American" subtype V (D at aa 390), and one from Honduras (1986) belonged to "Asian" subtype I. The data suggest that after its first isolation in the Caribbean in 1981, genotype IIIb spread throughout the Americas and effectively replaced subtype V throughout the Caribbean basin. The strain also evolved into several distinct lineages, based on substitutions in the E glycoprotein (amino acids 91 and 131), two of which were still in circulation in 2000. Interestingly, a molecular clock did not fit the data well, suggesting that other sources of rate variation, such as differential selection or differences in effective population sizes, may exist among lineages. Our results indicate the importance of large temporal- and geographical-scale phylogenetic studies in understanding disease dynamics, particularly where replacements between regions can occur.

Molecular analysis of the dengue virus type 1 and 2 in Brazil based on sequences of the genomic envelope-nonstructural protein 1 junction region

The genomic sequences of the Envelope-Non-Structural protein 1 junction region (E/NS1) of 84 DEN-1 and 22 DEN-2 isolates from Brazil were determined. Most of these strains were isolated in the period from 1995 to 2001 in endemic and regions of recent dengue transmission in São Paulo State. Sequence data for DEN-1 and DEN-2 utilized in phylogenetic and split decomposition analyses also include sequences deposited in GenBank from different regions of Brazil and of the world. Phylogenetic analyses were done using both maximum likelihood and Bayesian approaches. Results for both DEN-1 and DEN-2 data are ambiguous, and support for most tree bipartitions are generally poor, suggesting that E/NS1 region does not contain enough information for recovering phylogenetic relationships among DEN-1 and DEN-2 sequences used in this study. The network graph generated in the split decomposition analysis of DEN-1 does not show evidence of grouping sequences according to country, region and clades. While the network for DEN-2 also shows ambiguities among DEN-2 sequences, it suggests that Brazilian sequences may belong to distinct subtypes of genotype III.

A phylogenetic reconstruction of the epidemiological history of canine rabies virus variants in Colombia

Historically, canine rabies in Colombia has been caused by two geographically distinct canine variants of rabies virus (RV) which between 1992 and 2002 accounted for approximately 95% of Colombian rabies cases. Genetic variant 1 (GV1) has been isolated up until 1997 in the Central Region and the Department of Arauca, and is now considered extinct through a successful vaccination program. Genetic variant 2 (GV2) has been isolated from the northern Caribbean Region and continues to circulate at present. Here we have analyzed two sets of sequence data based upon either a 147 nucleotide region of the glycoprotein (G) gene or a 258 nucleotide region that combines a fragment of the non-coding intergenic region and a fragment of the polymerase gene. Using both maximum likelihood (ML) and Markov chain Monte Carlo (MCMC) methods we have estimated the time of the most recent common ancestor (MRCA) of the two variants to be between 1983 and 1988. Reconstructions of the population history suggest that GV2 has been circulating in Colombia since the 1960s and that GV1 evolved as a separate lineage from GV2. Estimations of the effective population size at present show the GV2 outbreak to be approximately 20 times greater than that of GV1. Demographic reconstructions were unable to detect a decrease in population size concurrent with the elimination of GV1. We find a raised rate of nucleotide substitution for GV1 gene sequences when compared to that of GV2, although all estimates have wide confidence limits. We demonstrate that phylogenetic reconstructions and sequence analysis can be used to support incidence data from the field in the assessment of RV epidemiology.

Origins, evolution, and vector/host coadaptations within the genus Flavivirus

Although viruses in the genus Flavivirus share complex antigenic interrelationships, they can be divided into four phylogenetic/ecological groups: two mosquito-borne groups, a tick-borne group, and nonvectored viruses. These divisions largely reflect the selective constraints imposed on the viruses by the vertebrate hosts, the invertebrate vectors, and the associated ecologies. Phylogenetic trees based on the flavivirus genetic sequence show characteristic branching patterns that reflect these groupings. This review describes the evolution and possible origins of individual flaviviruses, correlating ecological and epidemiological characteristics with their phylogenies and geographic dispersal. It will also become apparent that many of the phylogenetic lineages that define species diverged relatively recently, and the subsequent dispersal and epidemiology of these viruses have therefore been significantly influenced by increasing human population densities and activities such as recreation, urbanization, land reclamation, transportation, and deforestation. This review also considers some of the likely implications of persistent/chronic infections in relation to virus dispersal and recombination between related flaviviruses on phylogenetic analysis and vaccine development strategies.

E/NS1 modifications of dengue 2 virus after serial passages in mammalian and/or mosquito cells

OBJECTIVE

Dengue viruses are routinely maintained in nature by transmission cycles involving the passage of virus between humans and AEDES mosquitoes. The number of dengue virus lineages has been increasing over time. The aim of this study was to identify the genetic diversity of dengue 2 virus serially transferred in mammalian and/or mosquito cells.

METHODS

The E/NS1 gene of dengue 2 virus variants derived from serial passages in Vero or C6/36 cells, or alternately in both cell systems, was amplified and sequenced in order to observe gene modification after serial passages.

RESULTS

Three nucleotides (two in E and one in NS1) or two amino acids (one each in E and NS1) changed in the virus that was continuously cultured in Vero cells for 20 passages, whereas four nucleotides (two each in E and NS1) or three amino acids (one in E and two in NS1) changed in the virus cultured for 30 passages. The genome of dengue 2 virus remained stable even when the virus was serially transferred in C6/36 cells for 30 generations. However, there was one amino acid substitution (E46 I-->V) resulting from a single nucleotide change in the E region of dengue 2 virus alternately transferred in C6/36 and Vero cells for either 20 or 30 passages. In addition, dengue 2 virus obtained from serially cultured Vero cells usually replicated better when it reinfected Vero cells, reflecting its high adaptation fitness to the host cell.

CONCLUSIONS

It is concluded that genetic changes of dengue 2 virus are constrained in Vero (mammalian) cells, resulting in a variety of genome-related quasispecies populations. Some populations of the virus are subsequently selected by and genetically (at least in the E/NS1 portion of the viral genome) maintained in C6/36 (mosquito) cells during replicative competition.

The modern molecular clock

The discovery of the molecular clock--a relatively constant rate of molecular evolution--provided an insight into the mechanisms of molecular evolution, and created one of the most useful new tools in biology. The unexpected constancy of rate was explained by assuming that most changes to genes are effectively neutral. Theory predicts several sources of variation in the rate of molecular evolution. However, even an approximate clock allows time estimates of events in evolutionary history, which provides a method for testing a wide range of biological hypotheses ranging from the origins of the animal kingdom to the emergence of new viral epidemics.