Validation of high rates of nucleotide substitution in geminiviruses: phylogenetic evidence from East African cassava mosaic viruses

A novel recombinant variant of latent membrane protein 1 from Epstein Barr virus in Argentina denotes phylogeographical association

Epstein Barr virus (EBV) infection in Argentina occurs at an early age and occasionally develops infectious mononucleosis (IM). EBV is also related with lymphomas. LMP1, the viral oncoprotein is polymorphic and is used to define viral variants.

The diversification of begomovirus populations is predominantly driven by mutational dynamics

Begomoviruses (single-stranded DNA plant viruses) are responsible for serious agricultural threats. Begomovirus populations exhibit a high degree of within-host genetic variation and evolve as quickly as RNA viruses. Although the recombination-prone nature of begomoviruses has been extensively demonstrated, the relative contribution of recombination and mutation to the genetic variation of begomovirus populations has not been assessed. We estimated the genetic variability of begomovirus datasets from around the world. An uneven distribution of genetic variation across the length of the cp and rep genes due to recombination was evident from our analyses. To estimate the relative contributions of recombination and mutation to the genetic variability of begomoviruses, we mapped all substitutions over maximum likelihood trees and counted the number of substitutions on branches which were associated with recombination (η_r) and mutation (η_μ). In addition, we also estimated the per generation relative rates of both evolutionary mechanisms (r/μ) to express how frequently begomovirus genomes are affected by recombination relative to mutation. We observed that the composition of genetic variation in all begomovirus datasets was dominated by mutation. Additionally, the low correlation between the estimates indicated that the relative contributions of recombination and mutation are not necessarily a function of their relative rates. Our results show that, although a considerable fraction of the genetic variation levels could be assigned to recombination, it was always lower than that due to mutation, indicating that the diversification of begomovirus populations is predominantly driven by mutational dynamics.

Sequence-based taxonomic framework for the classification of uncultured single-stranded DNA viruses of the family Genomoviridae

With the advent of metagenomics approaches, a large diversity of known and unknown viruses has been identified in various types of environmental, plant, and animal samples. One such widespread virus group is the recently established family Genomoviridae which includes viruses with small (∼2-2.4 kb), circular ssDNA genomes encoding rolling-circle replication initiation proteins (Rep) and unique capsid proteins. Here, we propose a sequence-based taxonomic framework for classification of 121 new virus genomes within this family. Genomoviruses display ∼47% sequence diversity, which is very similar to that within the well-established and extensively studied family Geminiviridae (46% diversity). Based on our analysis, we establish a 78% genome-wide pairwise identity as a species demarcation threshold. Furthermore, using a Rep sequence phylogeny-based analysis coupled with the current knowledge on the classification of geminiviruses, we establish nine genera within the Genomoviridae family. These are Gemycircularvirus (n = 73), Gemyduguivirus (n = 1), Gemygorvirus (n = 9), Gemykibivirus (n = 29), Gemykolovirus (n = 3), Gemykrogvirus (n = 3), Gemykroznavirus (n = 1), Gemytondvirus (n = 1), Gemyvongvirus (n = 1). The presented taxonomic framework offers rational classification of genomoviruses based on the sequence information alone and sets an example for future classification of other groups of uncultured viruses discovered using metagenomics approaches.

Genetic variation and population structure of Cucumber green mottle mosaic virus

Cucumber green mottle mosaic virus (CGMMV) is a single-stranded, positive sense RNA virus infecting cucurbitaceous plants. In recent years, CGMMV has become an important pathogen of cucurbitaceous crops including watermelon, pumpkin, cucumber and bottle gourd in China, causing serious losses to their production. In this study, we surveyed CGMMV infection in various cucurbitaceous crops grown in Zhejiang Province and in several seed lots purchased from local stores with the dot enzyme-linked immunosorbent assay (dot-ELISA), using a CGMMV specific monoclonal antibody. Seven CGMMV isolates obtained from watermelon, grafted watermelon or oriental melon samples were cloned and sequenced. Identity analysis showed that the nucleotide identities of the seven complete genome sequences ranged from 99.2 to 100%. Phylogenetic analysis of seven CGMMV isolates as well as 24 other CGMMV isolates from the GenBank database showed that all CGMMV isolates could be grouped into two distinct monophyletic clades according to geographic distribution, i.e. Asian isolates for subtype I and European isolates for subtype II, indicating that population diversification of CGMMV isolates may be affected by geographical distribution. Site variation rate analysis of CGMMV found that the overall variation rate was below 8% and mainly ranged from 2 to 5%, indicating that the CGMMV genomic sequence was conservative. Base substitution type analysis of CGMMV showed a mutational bias, with more transitions (A↔G and C↔T) than transversions (A↔C, A↔T, G↔C and G↔T). Most of the variation occurring in the CGMMV genome resulted in non-synonymous substitutions, and the variation rate of some sites was higher than 30% because of this mutational bias. Selection constraint analysis of CGMMV ORFs showed strong negative selection acting on the replication-associated protein, similar to what occurs for other plant RNA viruses. Finally, potential recombination analysis identified isolate Ec as a recombinant with a low degree of confidence.

Divergent evolutionary and epidemiological dynamics of cassava mosaic geminiviruses in Madagascar

BACKGROUND

Cassava mosaic disease (CMD) in Madagascar is caused by a complex of at least six African cassava mosaic geminivirus (CMG) species. This provides a rare opportunity for a comparative study of the evolutionary and epidemiological dynamics of distinct pathogenic crop-infecting viral species that coexist within the same environment. The genetic and spatial structure of CMG populations in Madagascar was studied and Bayesian phylogeographic modelling was applied to infer the origins of Madagascan CMG populations within the epidemiological context of related populations situated on mainland Africa and other south western Indian Ocean (SWIO) islands.

RESULTS

The isolation and analysis of 279 DNA-A and 117 DNA-B sequences revealed the presence in Madagascar of four prevalent CMG species (South African cassava mosaic virus, SACMV; African cassava mosaic virus, ACMV; East African cassava mosaic Kenya virus, EACMKV; and East African cassava mosaic Cameroon virus, EACMCV), and of numerous CMG recombinants that have, to date, only ever been detected on this island. SACMV and ACMV, the two most prevalent viruses, displayed low degrees of genetic diversity and have most likely been introduced to the island only once. By contrast, EACMV-like CMG populations (consisting of East African cassava mosaic virus, EAMCKV, EACMCV and complex recombinants of these) were more diverse, more spatially structured, and displayed evidence of at least three independent introductions from mainland Africa. Although there were no statistically supported virus movement events between Madagascar and the other SWIO islands, at least one mainland African ACMV variant likely originated in Madagascar.

CONCLUSIONS

Our study highlights both the complexity of CMD in Madagascar, and the distinct evolutionary and spatial dynamics of the different viral species that collectively are associated with this disease. Given that more distinct CMG species and recombinants have been found in Madagascar than any other similarly sized region of the world, the risks of recombinant CMG variants emerging on this island are likely to be higher than elsewhere. Evidence of an epidemiological link between Madagascan and mainland African CMGs suggests that the consequences of such emergence events could reach far beyond the shores of this island.

Grapevine red blotch-associated virus is Present in Free-Living Vitis spp. Proximal to Cultivated Grapevines

Red blotch is an emerging disease of grapevine associated with grapevine red blotch-associated virus (GRBaV). The virus spreads with infected planting stocks but no vector of epidemiological significance has been conclusively identified. A vineyard block of red-blotch-affected Vitis vinifera 'Cabernet franc' clone 214 was observed in California, with a clustering of infected, symptomatic vines focused along one edge of the field proximal to a riparian habitat with free-living Vitis spp. No genetic heterogeneity was observed in a 587-nucleotide region of the GRBaV genome in a population of 44 Cabernet franc clone 214 isolates. By contrast, genetic differences were observed in isolates from other cultivars and clones growing in adjacent blocks. GRBaV was confirmed infecting four free-living vines, two of which were shown to be V. californica × V. vinifera hybrids. The genomes of three free-living GRBaV vine isolates and seven from V. vinifera cultivars were compared; free-living vine isolates were shown to be more similar to each other and a 'Merlot' isolate than to the other cultivated vine isolates. The finding that GRBaV is present in free-living Vitis spp. indicates the virus can be spread by natural (nonhuman-mediated) means, and we hypothesize that in-field spread of GRBaV is occurring.

Inferences from tip-calibrated phylogenies: a review and a practical guide

Molecular dating of phylogenetic trees is a growing discipline using sequence data to co‐estimate the timing of evolutionary events and rates of molecular evolution. All molecular‐dating methods require converting genetic divergence between sequences into absolute time. Historically, this could only be achieved by associating externally derived dates obtained from fossil or biogeographical evidence to internal nodes of the tree. In some cases, notably for fast‐evolving genomes such as viruses and some bacteria, the time span over which samples were collected may cover a significant proportion of the time since they last shared a common ancestor. This situation allows phylogenetic trees to be calibrated by associating sampling dates directly to the sequences representing the tips (terminal nodes) of the tree. The increasing availability of genomic data from ancient DNA extends the applicability of such tip‐based calibration to a variety of taxa including humans, extinct megafauna and various microorganisms which typically have a scarce fossil record. The development of statistical models accounting for heterogeneity in different aspects of the evolutionary process while accommodating very large data sets (e.g. whole genomes) has allowed using tip‐dating methods to reach inferences on divergence times, substitution rates, past demography or the age of specific mutations on a variety of spatiotemporal scales. In this review, we summarize the current state of the art of tip dating, discuss some recent applications, highlight common pitfalls and provide a ‘how to’ guide to thoroughly perform such analyses.

Genome Hotspots for Nucleotide Substitutions and the Evolution of Influenza A (H1N1) Human Strains

In recent years a number of studies have brought attention to the role of positive selection during the evolution of antigenic escape by influenza strains. Particularly, the identification of positively selected sites within antigenic domains of viral surface proteins has been used to suggest that the evolution of viral–host receptor binding specificity is driven by selection. Here we show that, following the 1918 outbreak, the antigenic sites of the hemagglutinin (HA) viral surface protein and the stalk region of neuraminidase became substitution hotspots. The hotspots show similar patterns of nucleotide substitution bias at synonymous and nonsynonymous sites. Such bias imposes directionality in amino acid replacements that can influence signals of selection at antigenic sites. Our results suggest that the high accumulation of substitutions within the antigenic sites of HA can explain not only cases of antigenic escape by antigenic drift but also lead to occasional episodes of viral extinction.

Two Novel DNAs That Enhance Symptoms and Overcome CMD2 Resistance to Cassava Mosaic Disease

Cassava mosaic begomoviruses (CMBs) cause cassava mosaic disease (CMD) across Africa and the Indian subcontinent. Like all members of the geminivirus family, CMBs have small, circular single-stranded DNA genomes. We report here the discovery of two novel DNA sequences, designated SEGS-1 and SEGS-2 (forsequencesenhancinggeminivirussymptoms), that enhance symptoms and break resistance to CMD. The SEGS are characterized by GC-rich regions and the absence of long open reading frames. Both SEGS enhanced CMD symptoms in cassava (Manihot esculentaCrantz) when coinoculated withAfrican cassava mosaic virus(ACMV),East African cassava mosaic Cameroon virus(EACMCV), orEast African cassava mosaic virus-Uganda(EACMV-UG). SEGS-1 also overcame resistance of a cassava landrace carrying the CMD2 resistance locus when coinoculated with EACMV-UG. Episomal forms of both SEGS were detected in CMB-infected cassava but not in healthy cassava. SEGS-2 episomes were also found in virions and whiteflies. SEGS-1 has no homology to geminiviruses or their associated satellites, but the cassava genome contains a sequence that is 99% identical to full-length SEGS-1. The cassava genome also includes three sequences with 84 to 89% identity to SEGS-2 that together encompass all of SEGS-2 except for a 52-bp region, which includes the episomal junction and a 26-bp sequence related to alphasatellite replication origins. These results suggest that SEGS-1 is derived from the cassava genome and facilitates CMB infection as an integrated copy and/or an episome, while SEGS-2 was originally from the cassava genome but now is encapsidated into virions and transmitted as an episome by whiteflies.

IMPORTANCE

Cassava is a major crop in the developing world, with its production in Africa being second only to maize. CMD is one of the most important diseases of cassava and a serious constraint to production across Africa. CMD2 is a major CMD resistance locus that has been deployed in many cassava cultivars through large-scale breeding programs. In recent years, severe, atypical CMD symptoms have been observed occasionally on resistant cultivars, some of which carry the CMD2 locus, in African fields. In this report, we identified and characterized two DNA sequences, SEGS-1 and SEGS-2, which produce similar symptoms when coinoculated with cassava mosaic begomoviruses onto a susceptible cultivar or a CMD2-resistant landrace. The ability of SEGS-1 to overcome CMD2 resistance and the transmission of SEGS-2 by whiteflies has major implications for the long-term durability of CMD2 resistance and underscore the need for alternative sources of resistance in cassava.

The recombination mediator RAD51D promotes geminiviral infection

To study a possible role for homologous recombination in geminivirus replication, we challenged Arabidopsis recombination gene knockouts by Euphorbia yellow mosaic virus infection. Our results show that the RAD51 paralog RAD51D, rather than RAD51 itself, promotes viral replication at early stages of infection. Blot hybridization analyses of replicative intermediates using one- and two-dimensional gels and deep sequencing point to an unexpected facet of recombination-dependent replication, the repair by single-strand annealing (SSA) during complementary strand replication. A significant decrease of both intramolecular, yielding defective DNAs and intermolecular recombinant molecules between the two geminiviral DNA components (A, B) were observed in the absence of RAD51D. By contrast, DNA A and B reacted differentially with the generation of inversions. A model to implicate single-strand annealing recombination in geminiviral recombination-dependent replication is proposed.

Diversity, Mutation and Recombination Analysis of Cotton Leaf Curl Geminiviruses

The spread of cotton leaf curl disease in China, India and Pakistan is a recent phenomenon. Analysis of available sequence data determined that there is a substantial diversity of cotton-infecting geminiviruses in Pakistan. Phylogenetic analyses indicated that recombination between two major groups of viruses, cotton leaf curl Multan virus (CLCuMuV) and cotton leaf curl Kokhran virus (CLCuKoV), led to the emergence of several new viruses. Recombination detection programs and phylogenetic analyses showed that CLCuMuV and CLCuKoV are highly recombinant viruses. Indeed, CLCuKoV appeared to be a major donor virus for the coat protein (CP) gene, while CLCuMuV donated the Rep gene in the majority of recombination events. Using recombination free nucleotide datasets the substitution rates for CP and Rep genes were determined. We inferred similar nucleotide substitution rates for the CLCuMuV-Rep gene (4.96X10^-4) and CLCuKoV-CP gene (2.706X10^-4), whereas relatively higher substitution rates were observed for CLCuMuV-CP and CLCuKoV-Rep genes. The combination of sequences with equal and relatively low substitution rates, seemed to result in the emergence of viral isolates that caused epidemics in Pakistan and India. Our findings also suggest that CLCuMuV is spreading at an alarming rate, which can potentially be a threat to cotton production in the Indian subcontinent.

Insight into virus encapsulation mechanism through in silico interaction study between coat protein and RNA operator loops of Sesbania mosaic virus

Viruses are parasite by nature and they are responsible for many diseases. Inhibitor development is very difficult for viruses due to their rapid mutative nature.

Host ecology determines the dispersal patterns of a plant virus

Since its isolation in 1966 in Kenya, rice yellow mottle virus (RYMV) has been reported throughout Africa resulting in one of the economically most important tropical plant emerging diseases. A thorough understanding of RYMV evolution and dispersal is critical to manage viral spread in tropical areas that heavily rely on agriculture for subsistence. Phylogenetic analyses have suggested a relatively recent expansion, perhaps driven by the intensification of agricultural practices, but this has not yet been examined in a coherent statistical framework. To gain insight into the historical spread of RYMV within Africa rice cultivations, we analyse a dataset of 300 coat protein gene sequences, sampled from East to West Africa over a 46-year period, using Bayesian evolutionary inference. Spatiotemporal reconstructions date the origin of RMYV back to 1852 (1791-1903) and confirm Tanzania as the most likely geographic origin. Following a single long-distance transmission event from East to West Africa, separate viral populations have been maintained for about a century. To identify the factors that shaped the RYMV distribution, we apply a generalised linear model (GLM) extension of discrete phylogenetic diffusion and provide strong support for distances measured on a rice connectivity landscape as the major determinant of RYMV spread. Phylogeographic estimates in continuous space further complement this by demonstrating more pronounced expansion dynamics in West Africa that are consistent with agricultural intensification and extensification. Taken together, our principled phylogeographic inference approach shows for the first time that host ecology dynamics have shaped the historical spread of a plant virus.

The contribution of translesion synthesis polymerases on geminiviral replication

Geminiviruses multiply primarily in the plant phloem, but never in meristems. Their Rep protein can activate DNA synthesis in differentiated cells. However, when their single-stranded DNA is injected into the phloem by insects, no Rep is present for inducing initial complementary strand replication. Considering a contribution of translesion synthesis (TLS) polymerases in plants, four of them (Polη, Polζ, Polκ, Rev1) are highly and constitutively expressed in differentiated tissues like the phloem. Two geminiviruses (Euphorbia yellow mosaic virus, Cleome leaf crumple virus), inoculated either biolistically or by whiteflies, replicated in Arabidopsis thaliana mutant lines of these genes to the same extent as in wild type plants. Comparative deep sequencing of geminiviral DNAs, however, showed a high exchange rate (10(-4)-10(-3)) similar to the phylogenetic variation described before and a significant difference in nucleotide substation rates if Polη and Polζ were absent, with a differential response to the viral DNA components.

Vector-Enabled Metagenomic (VEM) Surveys Using Whiteflies (Aleyrodidae) Reveal Novel Begomovirus Species in the New and Old Worlds

Whitefly-transmitted viruses belonging to the genus Begomovirus (family Geminiviridae) represent a substantial threat to agricultural food production. The rapid evolutionary potential of these single-stranded DNA viruses combined with the polyphagous feeding behavior of their whitefly vector (Bemisia tabaci) can lead to the emergence of damaging viral strains. Therefore, it is crucial to characterize begomoviruses circulating in different regions and crops globally. This study utilized vector-enabled metagenomics (VEM) coupled with high-throughput sequencing to survey begomoviruses directly from whiteflies collected in various locations (California (USA), Guatemala, Israel, Puerto Rico, and Spain). Begomoviruses were detected in all locations, with the highest diversity identified in Guatemala where up to seven different species were identified in a single field. Both bipartite and monopartite viruses were detected, including seven new begomovirus species from Guatemala, Puerto Rico, and Spain. This begomovirus survey extends the known diversity of these highly damaging plant viruses. However, the new genomes described here and in the recent literature appear to reflect the outcome of interactions between closely-related species, often resulting from recombination, instead of unique, highly divergent species.

Complexity of begomovirus and betasatellite populations associated with chilli leaf curl disease in India

Chilli, which encompasses several species in the genus Capsicum, is widely consumed throughout the world. In the Indian subcontinent, production of chilli is constrained due to chilli leaf curl disease (ChiLCD) caused by begomoviruses. Despite the considerable economic consequences of ChiLCD on chilli cultivation in India, there have been scant studies of the genetic diversity and structure of the begomoviruses that cause this disease. Here we report on a comprehensive survey across major chilli-growing regions in India. Analysis of samples collected in the survey indicates that ChiLCD-infected plants are associated with a complex of begomoviruses (including one previously unreported species) with a diverse group of betasatellites found in crops and weeds. The associated betasatellites neither enhanced the accumulation of the begomovirus components nor reduced the incubation period in Nicotiana benthamiana. The ChiLCD-associated begomoviruses induced mild symptoms on Capsicum spp., but both the level of helper virus that accumulated and the severity of symptoms were increased in the presence of cognate betasatellites. Interestingly, most of the begomoviruses were found to be intra-species recombinants. The betasatellites possess high nucleotide variability, and recombination among them was also evident. The nucleotide substitution rates were determined for the AV1 gene of begomoviruses (2.60 × 10- 3 substitutions site- 1 year- 1) and the βC1 gene of betasatellites [chilli leaf curl betasatellite (ChiLCB), 2.57 × 10- 4 substitution site- 1 year- 1; tomato leaf curl Bangladesh betasatellite (ToLCBDB), 5.22 × 10- 4 substitution site- 1 year- 1]. This study underscores the current understanding of Indian ChiLCD-associated begomoviruses and also demonstrates the crucial role of betasatellites in severe disease development in Capsicum spp.

The effect of genetic structure on molecular dating and tests for temporal signal

‘Dated‐tip’ methods of molecular dating use DNA sequences sampled at different times, to estimate the age of their most recent common ancestor. Several tests of ‘temporal signal’ are available to determine whether data sets are suitable for such analysis. However, it remains unclear whether these tests are reliable.

We investigate the performance of several tests of temporal signal, including some recently suggested modifications. We use simulated data (where the true evolutionary history is known), and whole genomes of methicillin‐resistant Staphylococcus aureus (to show how particular problems arise with real‐world data sets).

We show that all of the standard tests of temporal signal are seriously misleading for data where temporal and genetic structures are confounded (i.e. where closely related sequences are more likely to have been sampled at similar times). This is not an artefact of genetic structure or tree shape per se, and can arise even when sequences have measurably evolved during the sampling period. More positively, we show that a ‘clustered permutation’ approach introduced by Duchêne et al. (Molecular Biology and Evolution, 32, 2015, 1895) can successfully correct for this artefact in all cases and introduce techniques for implementing this method with real data sets.

The confounding of temporal and genetic structures may be difficult to avoid in practice, particularly for outbreaks of infectious disease, or when using ancient DNA. Therefore, we recommend the use of ‘clustered permutation’ for all analyses. The failure of the standard tests may explain why different methods of dating pathogen origins have reached such wildly different conclusions.

Novel circular single-stranded DNA viruses identified in marine invertebrates reveal high sequence diversity and consistent predicted intrinsic disorder patterns within putative structural proteins

Viral metagenomics has recently revealed the ubiquitous and diverse nature of single-stranded DNA (ssDNA) viruses that encode a conserved replication initiator protein (Rep) in the marine environment. Although eukaryotic circular Rep-encoding ssDNA (CRESS-DNA) viruses were originally thought to only infect plants and vertebrates, recent studies have identified these viruses in a number of invertebrates. To further explore CRESS-DNA viruses in the marine environment, this study surveyed CRESS-DNA viruses in various marine invertebrate species. A total of 27 novel CRESS-DNA genomes, with Reps that share less than 60.1% identity with previously reported viruses, were recovered from 21 invertebrate species, mainly crustaceans. Phylogenetic analysis based on the Rep revealed a novel clade of CRESS-DNA viruses that included approximately one third of the marine invertebrate associated viruses identified here and whose members may represent a novel family. Investigation of putative capsid proteins (Cap) encoded within the eukaryotic CRESS-DNA viral genomes from this study and those in GenBank demonstrated conserved patterns of predicted intrinsically disordered regions (IDRs), which can be used to complement similarity-based searches to identify divergent structural proteins within novel genomes. Overall, this study expands our knowledge of CRESS-DNA viruses associated with invertebrates and explores a new tool to evaluate divergent structural proteins encoded by these viruses.

Multiple layers of chimerism in a single-stranded DNA virus discovered by deep sequencing

Viruses with single-stranded (ss) DNA genomes infect hosts in all three domains of life and include many medically, ecologically, and economically important pathogens. Recently, a new group of ssDNA viruses with chimeric genomes has been discovered through viral metagenomics. These chimeric viruses combine capsid protein genes and replicative protein genes that, respectively, appear to have been inherited from viruses with positive-strand RNA genomes, such as tombusviruses, and ssDNA genomes, such as circoviruses, nanoviruses or geminiviruses. Here, we describe the genome sequence of a new representative of this virus group and reveal an additional layer of chimerism among ssDNA viruses. We show that not only do these viruses encompass genes for capsid proteins and replicative proteins that have distinct evolutionary histories, but also the replicative genes themselves are chimeras of functional domains inherited from viruses of different families. Our results underscore the importance of horizontal gene transfer in the evolution of ssDNA viruses and the role of genetic recombination in the emergence of novel virus groups.

The Performance of the Date-Randomization Test in Phylogenetic Analyses of Time-Structured Virus Data

Rates and timescales of viral evolution can be estimated using phylogenetic analyses of time-structured molecular sequences. This involves the use of molecular-clock methods, calibrated by the sampling times of the viral sequences. However, the spread of these sampling times is not always sufficient to allow the substitution rate to be estimated accurately. We conducted Bayesian phylogenetic analyses of simulated virus data to evaluate the performance of the date-randomization test, which is sometimes used to investigate whether time-structured data sets have temporal signal. An estimate of the substitution rate passes this test if its mean does not fall within the 95% credible intervals of rate estimates obtained using replicate data sets in which the sampling times have been randomized. We find that the test sometimes fails to detect rate estimates from data with no temporal signal. This error can be minimized by using a more conservative criterion, whereby the 95% credible interval of the estimate with correct sampling times should not overlap with those obtained with randomized sampling times. We also investigated the behavior of the test when the sampling times are not uniformly distributed throughout the tree, which sometimes occurs in empirical data sets. The test performs poorly in these circumstances, such that a modification to the randomization scheme is needed. Finally, we illustrate the behavior of the test in analyses of nucleotide sequences of cereal yellow dwarf virus. Our results validate the use of the date-randomization test and allow us to propose guidelines for interpretation of its results.

Recombination in viruses: mechanisms, methods of study, and evolutionary consequences

Recombination is a pervasive process generating diversity in most viruses. It joins variants that arise independently within the same molecule, creating new opportunities for viruses to overcome selective pressures and to adapt to new environments and hosts. Consequently, the analysis of viral recombination attracts the interest of clinicians, epidemiologists, molecular biologists and evolutionary biologists. In this review we present an overview of three major areas related to viral recombination: (i) the molecular mechanisms that underlie recombination in model viruses, including DNA-viruses (Herpesvirus) and RNA-viruses (Human Influenza Virus and Human Immunodeficiency Virus), (ii) the analytical procedures to detect recombination in viral sequences and to determine the recombination breakpoints, along with the conceptual and methodological tools currently used and a brief overview of the impact of new sequencing technologies on the detection of recombination, and (iii) the major areas in the evolutionary analysis of viral populations on which recombination has an impact. These include the evaluation of selective pressures acting on viral populations, the application of evolutionary reconstructions in the characterization of centralized genes for vaccine design, and the evaluation of linkage disequilibrium and population structure.

Bioinformatic analysis reveals genome size reduction and the emergence of tyrosine phosphorylation site in the movement protein of New World bipartite begomoviruses

Begomovirus (genus Begomovirus, family Geminiviridae) infection is devastating to a wide variety of agricultural crops including tomato, squash, and cassava. Thus, understanding the replication and adaptation of begomoviruses has important translational value in alleviating substantial economic loss, particularly in developing countries. The bipartite genome of begomoviruses prevalent in the New World and their counterparts in the Old World share a high degree of genome homology except for a partially overlapping reading frame encoding the pre-coat protein (PCP, or AV2). PCP contributes to the essential functions of intercellular movement and suppression of host RNA silencing, but it is only present in the Old World viruses. In this study, we analyzed a set of non-redundant bipartite begomovirus genomes originating from the Old World (N = 28) and the New World (N = 65). Our bioinformatic analysis suggests ∼ 120 nucleotides were deleted from PCP's proximal promoter region that may have contributed to its loss in the New World viruses. Consequently, genomes of the New World viruses are smaller than the Old World counterparts, possibly compensating for the loss of the intercellular movement functions of PCP. Additionally, we detected substantial purifying selection on a portion of the New World DNA-B movement protein (MP, or BC1). Further analysis of the New World MP gene revealed the emergence of a putative tyrosine phosphorylation site, which likely explains the increased purifying selection in that region. These findings provide important information about the strategies adopted by bipartite begomoviruses in adapting to new environment and suggest future in planta experiments.

Evolution of cocirculating varicella-zoster virus genotypes during a chickenpox outbreak in Guinea-Bissau

Varicella-zoster virus (VZV), a double-stranded DNA alphaherpesvirus, is associated with seasonal outbreaks of varicella in nonimmunized populations. Little is known about whether these outbreaks are associated with a single or multiple viral genotypes and whether new mutations rapidly accumulate during transmission. Here, we take advantage of a well-characterized population cohort in Guinea-Bissau and produce a unique set of 23 full-length genome sequences, collected over 7 months from eight households. Comparative sequence analysis reveals that four distinct genotypes cocirculated among the population, three of which were present during the first week of the outbreak, although no patients were coinfected, which indicates that exposure to infectious virus from multiple sources is common during VZV outbreaks. Transmission of VZV was associated with length polymorphisms in the R1 repeat region and the origin of DNA replication. In two cases, these were associated with the formation of distinct lineages and point to the possible coevolution of these loci, despite the lack of any known functional link in VZV or related herpesviruses. We show that these and all other sequenced clade 5 viruses possess a distinct R1 repeat motif that increases the acidity of an ORF11p protein domain and postulate that this has either arisen or been lost following divergence of the major clades. Thus, sequencing of whole VZV genomes collected during an outbreak has provided novel insights into VZV biology, transmission patterns, and (recent) natural history.

IMPORTANCE VZV is a highly infectious virus and the causative agent of chickenpox and shingles, the latter being particularly associated with the risk of painful complications. Seasonal outbreaks of chickenpox are very common among young children, yet little is known about the dynamics of the virus during person-to-person to transmission or whether multiple distinct viruses seed and/or cocirculate during an outbreak. In this study, we have sequenced chickenpox viruses from an outbreak in Guinea-Bissau that are supported by detailed epidemiological data. Our data show that multiple different virus strains seeded and were maintained throughout the 6-month outbreak period and that viruses transmitted between individuals accumulated new mutations in specific genomic regions. Of particular interest is the potential coevolution of two distinct parts of the genomes and our calculations of the rate of viral mutation, both of which increase our understanding of how VZV evolves over short periods of time in human populations.

A novel papillomavirus in Adélie penguin (Pygoscelis adeliae) faeces sampled at the Cape Crozier colony, Antarctica

Papillomaviruses are epitheliotropic viruses that have circular dsDNA genomes encapsidated in non-enveloped virions. They have been found to infect a variety of mammals, reptiles and birds, but so far they have not been found in amphibians. Using a next-generation sequencing de novo assembly contig-informed recovery, we cloned and Sanger sequenced the complete genome of a novel papillomavirus from the faecal matter of Adélie penguins (Pygoscelis adeliae) nesting on Ross Island, Antarctica. The genome had all the usual features of a papillomavirus and an E9 ORF encoding a protein of unknown function that is found in all avian papillomaviruses to date. This novel papillomavirus genome shared ~60 % pairwise identity with the genomes of the other three known avian papillomaviruses: Fringilla coelebs papillomavirus 1 (FcPV1), Francolinus leucoscepus papillomavirus 1 (FlPV1) and Psittacus erithacus papillomavirus 1. Pairwise identity analysis and phylogenetic analysis of the major capsid protein gene clearly indicated that it represents a novel species, which we named Pygoscelis adeliae papillomavirus 1 (PaCV1). No evidence of recombination was detected in the genome of PaCV1, but we did detect a recombinant region (119 nt) in the E6 gene of FlPV1 with the recombinant region being derived from ancestral FcPV1-like sequences. Previously only paramyxoviruses, orthomyxoviruses and avian pox viruses have been genetically identified in penguins; however, the majority of penguin viral identifications have been based on serology or histology. This is the first report, to our knowledge, of a papillomavirus associated with a penguin species.

Pigeon circoviruses display patterns of recombination, genomic secondary structure and selection similar to those of beak and feather disease viruses

Pigeon circovirus (PiCV) has a ~2 kb genome circular ssDNA genome. All but one of the known PiCV isolates have been found infecting pigeons in various parts of the world. In this study, we screened 324 swab and tissue samples from Polish pigeons and recovered 30 complete genomes, 16 of which came from birds displaying no obvious pathology. Together with 17 other publicly available PiCV complete genomes sampled throughout the Northern Hemisphere and Australia, we find that PiCV displays a similar degree of genetic diversity to that of the related psittacine-infecting circovirus species, beak and feather disease virus (BFDV). We show that, as is the case with its pathology and epidemiology, PiCV also displays patterns of recombination, genomic secondary structure and natural selection that are generally very similar to those of BFDV. It is likely that breeding facilities play a significant role in the emergence of new recombinant PiCV variants and given that ~50 % of the domestic pigeon population is infected subclinically, all pigeon breeding stocks should be screened routinely for this virus.

The effect of mutation and selection on codon adaptation in Escherichia coli bacteriophage

Studying phage codon adaptation is important not only for understanding the process of translation elongation, but also for reengineering phages for medical and industrial purposes. To evaluate the effect of mutation and selection on phage codon usage, we developed an index to measure selection imposed by host translation machinery, based on the difference in codon usage between all host genes and highly expressed host genes. We developed linear and nonlinear models to estimate the C→T mutation bias in different phage lineages and to evaluate the relative effect of mutation and host selection on phage codon usage. C→T-biased mutations occur more frequently in single-stranded DNA (ssDNA) phages than in double-stranded DNA (dsDNA) phages and affect not only synonymous codon usage, but also nonsynonymous substitutions at second codon positions, especially in ssDNA phages. The host translation machinery affects codon adaptation in both dsDNA and ssDNA phages, with a stronger effect on dsDNA phages than on ssDNA phages. Strand asymmetry with the associated local variation in mutation bias can significantly interfere with codon adaptation in both dsDNA and ssDNA phages.

De novo reconstruction of consensus master genomes of plant RNA and DNA viruses from siRNAs

Virus-infected plants accumulate abundant, 21-24 nucleotide viral siRNAs which are generated by the evolutionary conserved RNA interference (RNAi) machinery that regulates gene expression and defends against invasive nucleic acids. Here we show that, similar to RNA viruses, the entire genome sequences of DNA viruses are densely covered with siRNAs in both sense and antisense orientations. This implies pervasive transcription of both coding and non-coding viral DNA in the nucleus, which generates double-stranded RNA precursors of viral siRNAs. Consistent with our finding and hypothesis, we demonstrate that the complete genomes of DNA viruses from Caulimoviridae and Geminiviridae families can be reconstructed by deep sequencing and de novo assembly of viral siRNAs using bioinformatics tools. Furthermore, we prove that this 'siRNA omics' approach can be used for reliable identification of the consensus master genome and its microvariants in viral quasispecies. Finally, we utilized this approach to reconstruct an emerging DNA virus and two viroids associated with economically-important red blotch disease of grapevine, and to rapidly generate a biologically-active clone representing the wild type master genome of Oilseed rape mosaic virus. Our findings show that deep siRNA sequencing allows for de novo reconstruction of any DNA or RNA virus genome and its microvariants, making it suitable for universal characterization of evolving viral quasispecies as well as for studying the mechanisms of siRNA biogenesis and RNAi-based antiviral defense.

A complete ancient RNA genome: identification, reconstruction and evolutionary history of archaeological Barley Stripe Mosaic Virus

The origins of many plant diseases appear to be recent and associated with the rise of domestication, the spread of agriculture or recent global movements of crops. Distinguishing between these possibilities is problematic because of the difficulty of determining rates of molecular evolution over short time frames. Heterochronous approaches using recent and historical samples show that plant viruses exhibit highly variable and often rapid rates of molecular evolution. The accuracy of estimated evolution rates and age of origin can be greatly improved with the inclusion of older molecular data from archaeological material. Here we present the first reconstruction of an archaeological RNA genome, which is of Barley Stripe Mosaic Virus (BSMV) isolated from barley grain ~750 years of age. Phylogenetic analysis of BSMV that includes this genome indicates the divergence of BSMV and its closest relative prior to this time, most likely around 2000 years ago. However, exclusion of the archaeological data results in an apparently much more recent origin of the virus that postdates even the archaeological sample. We conclude that this viral lineage originated in the Near East or North Africa, and spread to North America and East Asia with their hosts along historical trade routes.

Molecular variability and evolution of a natural population of tomato yellow leaf curl virus in Shanghai, China

Tomato yellow leaf curl virus (TYLCV), belonging to the genus Begomovirus of the family Geminiviridae, is emerging as the most destructive pathogen of tomato plants. Since the first report of TYLCV in Shanghai, China in 2006, TYLCV has spread rapidly to 13 provinces or autonomous regions of China. In this study, the molecular variability and evolution of TYLCV were monitored in Shanghai from its first upsurge in 2006 until 2010. Full-length genomic sequences of 26 isolates were obtained by rolling circle amplification. Sequence analysis showed that the intergenic region was the most variable, with a mean mutation rate of 4.81×10(-3) nucleotide substitutions per site per year. Genetic differentiation was found within isolates obtained from 2006, 2009, and 2010, though a linear increase in genetic diversity over time was not evident. Whilst significant parts of TYLCV genes were under negative selection, the C4 gene embedded entirely within the C1 gene had a tendency to undergo positive selection. Our results indicate that a mechanism of independent evolution of overlapping regions could apply to the natural population of TYLCV in Shanghai, China.

Cell tropism predicts long-term nucleotide substitution rates of mammalian RNA viruses

The high rates of RNA virus evolution are generally attributed to replication with error-prone RNA-dependent RNA polymerases. However, these long-term nucleotide substitution rates span three orders of magnitude and do not correlate well with mutation rates or selection pressures. This substitution rate variation may be explained by differences in virus ecology or intrinsic genomic properties. We generated nucleotide substitution rate estimates for mammalian RNA viruses and compiled comparable published rates, yielding a dataset of 118 substitution rates of structural genes from 51 different species, as well as 40 rates of non-structural genes from 28 species. Through ANCOVA analyses, we evaluated the relationships between these rates and four ecological factors: target cell, transmission route, host range, infection duration; and three genomic properties: genome length, genome sense, genome segmentation. Of these seven factors, we found target cells to be the only significant predictors of viral substitution rates, with tropisms for epithelial cells or neurons (P<0.0001) as the most significant predictors. Further, one-tailed t-tests showed that viruses primarily infecting epithelial cells evolve significantly faster than neurotropic viruses (P<0.0001 and P<0.001 for the structural genes and non-structural genes, respectively). These results provide strong evidence that the fastest evolving mammalian RNA viruses infect cells with the highest turnover rates: the highly proliferative epithelial cells. Estimated viral generation times suggest that epithelial-infecting viruses replicate more quickly than viruses with different cell tropisms. Our results indicate that cell tropism is a key factor in viral evolvability.

RNA viruses are the fastest evolving human pathogens, making their treatment and control difficult. Compared to DNA viruses, RNA viruses replicate with much lower fidelity, which can explain why RNA viruses evolve significantly faster than most DNA viruses. However, there is tremendous variation among the evolutionary rates of different RNA viruses, which is not explained by variation in mutation rates. Here we present a survey of mammalian RNA virus rates of evolution, and a comprehensive comparison of these rates to different properties of virus genomic architecture and ecology. We found that cell tropism is the most significant predictor of long-term rates of mammalian RNA virus evolution. For instance, viruses targeting epithelial cells evolve significantly faster than viruses that target neurons. Our results provide mechanistic insight into why viruses that infect respiratory and gastrointestinal epithelia have been difficult to control.

Mutability dynamics of an emergent single stranded DNA virus in a naïve host

Quasispecies variants and recombination were studied longitudinally in an emergent outbreak of beak and feather disease virus (BFDV) infection in the orange-bellied parrot (Neophema chrysogaster). Detailed health monitoring and the small population size (<300 individuals) of this critically endangered bird provided an opportunity to longitudinally track viral replication and mutation events occurring in a circular, single-stranded DNA virus over a period of four years within a novel bottleneck population. Optimized PCR was used with different combinations of primers, primer walking, direct amplicon sequencing and sequencing of cloned amplicons to analyze BFDV genome variants. Analysis of complete viral genomes (n = 16) and Rep gene sequences (n = 35) revealed that the outbreak was associated with mutations in functionally important regions of the normally conserved Rep gene and immunogenic capsid (Cap) gene with a high evolutionary rate (3.41×10⁻³ subs/site/year) approaching that for RNA viruses; simultaneously we observed significant evidence of recombination hotspots between two distinct progenitor genotypes within orange-bellied parrots indicating early cross-transmission of BFDV in the population. Multiple quasispecies variants were also demonstrated with at least 13 genotypic variants identified in four different individual birds, with one containing up to seven genetic variants. Preferential PCR amplification of variants was also detected. Our findings suggest that the high degree of genetic variation within the BFDV species as a whole is reflected in evolutionary dynamics within individually infected birds as quasispecies variation, particularly when BFDV jumps from one host species to another.

Brazilian begomovirus populations are highly recombinant, rapidly evolving, and segregated based on geographical location

The incidence of begomovirus infections in crop plants sharply increased in Brazil during the 1990s following the introduction of the invasive B biotype of the whitefly vector, Bemisia tabaci. It is believed that this biotype transmitted begomoviruses from noncultivated plants to crop species with greater efficiency than indigenous B. tabaci biotypes. Either through rapid host adaptation or selection pressure in genetically diverse populations of noncultivated hosts, over the past 20 years various previously unknown begomovirus species have became progressively more prevalent in cultivated species such as tomato. Here we assess the genetic structure of begomovirus populations infecting tomatoes and noncultivated hosts in southeastern Brazil. Between 2005 and 2010, we sampled and sequenced 126 DNA-A and 58 DNA-B full-length begomovirus components. We detected nine begomovirus species in tomatoes and eight in the noncultivated host samples, with four species common to both tomatoes and noncultivated hosts. Like many begomoviruses, most species are obvious interspecies recombinants. Furthermore, species identified in tomato have probable parental viruses from noncultivated hosts. While the population structures of five well-sampled viral species all displayed geographical subdivision, a noncultivated host-infecting virus was more genetically variable than the four predominantly tomato-infecting viruses.

Fulfilling Koch's postulates confirms the monopartite nature of tomato leaf deformation virus: a begomovirus native to the New World

The monopartite nature of the begomovirus tomato leaf deformation virus (ToLDeV) reported in Peru is demonstrated here. The DNA molecule cloned from an infected plant was shown to be fully infectious in tomatoes inducing leaf curling and stunted growth similar to that observed in field-infected plants. The viral DNA was reisolated from systemically infected tissues of inoculated plants, thus fulfilling Koch's postulates. ToLDeV was demonstrated, therefore, as the causal agent of the disease syndrome widespread in tomato crops in Peru. This virus was shown to be present throughout the major tomato-growing regions of this country, both in tomatoes and wild plants. Analyses of the sequences of 51 ToLDeV isolates revealed a significant genetic diversity with three major genetic types co-circulating in the population. A geographical segregation was observed which should be taken into account for virus control. Constraints to genetic divergence found for the C4 gene of ToLDeV isolates suggest a relevant function for this protein. The results obtained confirm ToLDeV as a monopartite begomovirus native to the New World, which is a significant finding for this region.

Synonymous site variation due to recombination explains higher genetic variability in begomovirus populations infecting non-cultivated hosts

Begomoviruses are ssDNA plant viruses that cause serious epidemics in economically important crops worldwide. Non-cultivated plants also harbour many begomoviruses, and it is believed that these hosts may act as reservoirs and as mixing vessels where recombination may occur. Begomoviruses are notoriously recombination-prone, and also display nucleotide substitution rates equivalent to those of RNA viruses. In Brazil, several indigenous begomoviruses have been described infecting tomatoes following the introduction of a novel biotype of the whitefly vector in the mid-1990s. More recently, a number of viruses from non-cultivated hosts have also been described. Previous work has suggested that viruses infecting non-cultivated hosts have a higher degree of genetic variability compared with crop-infecting viruses. We intensively sampled cultivated and non-cultivated plants in similarly sized geographical areas known to harbour either the weed-infecting Macroptilium yellow spot virus (MaYSV) or the crop-infecting Tomato severe rugose virus (ToSRV), and compared the molecular evolution and population genetics of these two distantly related begomoviruses. The results reinforce the assertion that infection of non-cultivated plant species leads to higher levels of standing genetic variability, and indicate that recombination, not adaptive selection, explains the higher begomovirus variability in non-cultivated hosts.

Base composition and translational selection are insufficient to explain codon usage bias in plant viruses

Viral codon usage bias may be the product of a number of synergistic or antagonistic factors, including genomic nucleotide composition, translational selection, genomic architecture, and mutational or repair biases. Most studies of viral codon bias evaluate only the relative importance of genomic base composition and translational selection, ignoring other possible factors. We analyzed the codon preferences of ssRNA (luteoviruses and potyviruses) and ssDNA (geminiviruses) plant viruses that infect translationally distinct monocot and dicot hosts. We found that neither genomic base composition nor translational selection satisfactorily explains their codon usage biases. Furthermore, we observed a strong relationship between the codon preferences of viruses in the same family or genus, regardless of host or genomic nucleotide content. Our results suggest that analyzing codon bias as either due to base composition or translational selection is a false dichotomy that obscures the role of other factors. Constraints such as genomic architecture and secondary structure can and do influence codon usage in plant viruses, and likely in viruses of other hosts.

Adaptive evolution by recombination is not associated with increased mutation rates in Maize streak virus

BACKGROUND

Single-stranded (ss) DNA viruses in the family Geminiviridae are proving to be very useful in real-time evolution studies. The high mutation rate of geminiviruses and other ssDNA viruses is somewhat mysterious in that their DNA genomes are replicated in host nuclei by high fidelity host polymerases. Although strand specific mutation biases observed in virus species from the geminivirus genus Mastrevirus indicate that the high mutation rates in viruses in this genus may be due to mutational processes that operate specifically on ssDNA, it is currently unknown whether viruses from other genera display similar strand specific mutation biases. Also, geminivirus genomes frequently recombine with one another and an alternative cause of their high mutation rates could be that the recombination process is either directly mutagenic or produces a selective environment in which the survival of mutants is favoured. To investigate whether there is an association between recombination and increased basal mutation rates or increased degrees of selection favoring the survival of mutations, we compared the mutation dynamics of the MSV-MatA and MSV-VW field isolates of Maize streak virus (MSV; Mastrevirus), with both a laboratory constructed MSV recombinant, and MSV recombinants closely resembling MSV-MatA. To determine whether strand specific mutation biases are a general characteristic of geminivirus evolution we compared mutation spectra arising during these MSV experiments with those arising during similar experiments involving the geminivirus Tomato yellow leaf curl virus (Begomovirus genus).

RESULTS

Although both the genomic distribution of mutations and the occurrence of various convergent mutations at specific genomic sites indicated that either mutation hotspots or selection for adaptive mutations might elevate observed mutation rates in MSV, we found no association between recombination and mutation rates. Importantly, when comparing the mutation spectra of MSV and TYLCV we observed similar strand specific mutation biases arising predominantly from imbalances in the complementary mutations G → T: C → A.

CONCLUSIONS

While our results suggest that recombination does not strongly influence mutation rates in MSV, they indicate that high geminivirus mutation rates are at least partially attributable to increased susceptibility of all geminivirus genomes to oxidative damage while in a single stranded state.

East African cassava mosaic-like viruses from Africa to Indian ocean islands: molecular diversity, evolutionary history and geographical dissemination of a bipartite begomovirus

BACKGROUND

Cassava (Manihot esculenta) is a major food source for over 200 million sub-Saharan Africans. Unfortunately, its cultivation is severely hampered by cassava mosaic disease (CMD). Caused by a complex of bipartite cassava mosaic geminiviruses (CMG) species (Family: Geminivirideae; Genus: Begomovirus) CMD has been widely described throughout Africa and it is apparent that CMG's are expanding their geographical distribution. Determining where and when CMG movements have occurred could help curtail its spread and reveal the ecological and anthropic factors associated with similar viral invasions. We applied Bayesian phylogeographic inference and recombination analyses to available and newly described CMG sequences to reconstruct a plausible history of CMG diversification and migration between Africa and South West Indian Ocean (SWIO) islands.

RESULTS

The isolation and analysis of 114 DNA-A and 41 DNA-B sequences demonstrated the presence of three CMG species circulating in the Comoros and Seychelles archipelagos (East African cassava mosaic virus, EACMV; East African cassava mosaic Kenya virus, EACMKV; and East African cassava mosaic Cameroon virus, EACMCV). Phylogeographic analyses suggest that CMG's presence on these SWIO islands is probably the result of at least four independent introduction events from mainland Africa occurring between 1988 and 2009. Amongst the islands of the Comoros archipelago, two major migration pathways were inferred: One from Grande Comore to Mohéli and the second from Mayotte to Anjouan. While only two recombination events characteristic of SWIO islands isolates were identified, numerous re-assortments events were detected between EACMV and EACMKV, which seem to almost freely interchange their genome components.

CONCLUSIONS

Rapid and extensive virus spread within the SWIO islands was demonstrated for three CMG complex species. Strong evolutionary or ecological interaction between CMG species may explain both their propensity to exchange components and the absence of recombination with non-CMG begomoviruses. Our results suggest an important role of anthropic factors in CMGs spread as the principal axes of viral migration correspond with major routes of human movement and commercial trade. Finer-scale temporal analyses of CMGs to precisely scale the relative contributions of human and insect transmission to their movement dynamics will require further extensive sampling in the SWIO region.

A melting pot of Old World begomoviruses and their satellites infecting a collection of Gossypium species in Pakistan

CLCuD in southern Asia is caused by a complex of multiple begomoviruses (whitefly transmitted, single-stranded [ss]DNA viruses) in association with a specific ssDNA satellite; Cotton leaf curl Multan betasatellite (CLCuMuB). A further single ssDNA molecule, for which the collective name alphasatellites has been proposed, is also frequently associated with begomovirus-betasatellite complexes. Multan is in the center of the cotton growing area of Pakistan and has seen some of the worst problems caused by CLCuD. An exhaustive analysis of the diversity of begomoviruses and their satellites occurring in 15 Gossypium species (including G. hirsutum, the mainstay of Pakistan's cotton production) that are maintained in an orchard in the vicinity of Multan has been conducted using φ29 DNA polymerase-mediated rolling-circle amplification, cloning and sequence analysis. The non-cultivated Gossypium species, including non-symptomatic plants, were found to harbor a much greater diversity of begomoviruses and satellites than found in the cultivated G. hirsutum. Furthermore an African cassava mosaic virus (a virus previously only identified in Africa) DNA-A component and a Jatropha curcas mosaic virus (a virus occurring only in southern India) DNA-B component were identified. Consistent with earlier studies of cotton in southern Asia, only a single species of betasatellite, CLCuMuB, was identified. The diversity of alphasatellites was much greater, with many previously unknown species, in the non-cultivated cotton species than in G. hirsutum. Inoculation of newly identified components showed them to be competent for symptomatic infection of Nicotiana benthamiana plants. The significance of the findings with respect to our understanding of the role of host selection in virus diversity in crops and the geographical spread of viruses by human activity are discussed.

One misdated sequence of rabbit hemorrhagic disease virus prevents accurate estimation of its nucleotide substitution rate

BACKGROUND

The literature is ripe with phylogenetic estimates of nucleotide substitution rates, especially of measurably evolving species such as RNA viruses. However, it is not known how robust these rate estimates are to inaccuracies in the data, particularly in sampling dates that are used for molecular clock calibration. Here we report on the rate of evolution of the emerging pathogen Rabbit hemorrhagic disease virus (RHDV), which has significantly different rates of evolution for the same outer capsid (VP60) gene published in the literature. In an attempt to reconcile the conflicting data and further elucidate details of RHDV 's evolutionary history, we undertook fresh Bayesian analyses and employed jackknife control methods to produce robust substitution rate and time to most recent common ancestor (TMRCA) estimates for RHDV based on the VP60 and RNA-dependent RNA polymerase genes.

RESULTS

Through these control methods, we were able to identify a single misdated taxon, a passaged lab strain used for vaccine production, which was responsible for depressing the RHDV capsid gene's rate of evolution by 65%. Without this isolate, the polymerase and the capsid protein genes had nearly identical rates of evolution: 1.90x10-3 nucleotide substitutions/site/year, ns/s/y, (95% highest probability density (HPD) 1.25x10-3-2.55x10-3) and 1.91x10-3 ns/s/y (95% HPD 1.50x10-3-2.34x10-3), respectively.

CONCLUSIONS

After excluding the misdated taxon, both genes support a significantly higher substitution rate as well as a relatively recent emergence of RHDV, and obviate the need for previously hypothesized decades of unobserved diversification of the virus. The control methods show that using even one misdated taxon in a large dataset can significantly skew estimates of evolutionary parameters and suggest that it is better practice to use smaller datasets composed of taxa with unequivocal isolation dates. These jackknife controls would be useful for future tip-calibrated rate analyses that include taxa with ambiguous dates of isolation.

The VIRCA Project: virus resistant cassava for Africa

The VIRCA (Virus Resistant Cassava for Africa) project is a collaborative program between the Donald Danforth Plant Science Center, USA the National Crops Resources Research Institute, Uganda and the Kenya Agricultural Research Institute, Kenya. VIRCA is structured to include all aspects of the intellectual property, technology, regulatory, biosafety, quality control, communication and distribution components required for a GM crop development and delivery process. VIRCA's goal is to improve cassava for resistance to the viral diseases cassava brown streak disease (CBSD) and cassava mosaic disease (CMD) using pathogen-derived RNAi technology, and to field test, obtain regulatory approval for and deliver these products to small landholder farmers. During Phase I of the project, proof of concept was achieved by production and testing of virus resistant plants under greenhouse and confined field trials in East Africa. In VIRCA Phase II, two farmer-preferred varieties will be modified for resistance to CBSD and CMD, and lead events identified after molecular and field screening. In addition to delivery of royalty-free improved planting materials for farmers, VIRCA capacity building activities are enhancing indigenous capability for crop biotechnology in East Africa.

Tracking viral evolution during a disease outbreak: the rapid and complete selective sweep of a circovirus in the endangered Echo parakeet

Circoviruses are among the smallest and simplest of all viruses, but they are relatively poorly characterized. Here, we intensively sampled two sympatric parrot populations from Mauritius over a period of 11 years and screened for the circovirus Beak and feather disease virus (BFDV). During the sampling period, a severe outbreak of psittacine beak and feather disease, which is caused by BFDV, occurred in Echo parakeets. Consequently, this data set presents an ideal system for studying the evolution of a pathogen in a natural population and to understand the adaptive changes that cause outbreaks. Unexpectedly, we discovered that the outbreak was most likely caused by changes in functionally important regions of the normally conserved replication-associated protein gene and not the immunogenic capsid. Moreover, these mutations were completely fixed in the Echo parakeet host population very shortly after the outbreak. Several capsid alleles were linked to the replication-associated protein outbreak allele, suggesting that whereas the key changes occurred in the latter, the scope of the outbreak and the selective sweep may have been influenced by positive selection in the capsid. We found evidence for viral transmission between the two host populations though evidence for the invasive species as the source of the outbreak was equivocal. Finally, the high evolutionary rate that we estimated shows how rapidly new variation can arise in BFDV and is consistent with recent results from other small single-stranded DNA viruses.

Ultrafast evolution and loss of CRISPRs following a host shift in a novel wildlife pathogen, Mycoplasma gallisepticum

Measureable rates of genome evolution are well documented in human pathogens but are less well understood in bacterial pathogens in the wild, particularly during and after host switches. Mycoplasma gallisepticum (MG) is a pathogenic bacterium that has evolved predominantly in poultry and recently jumped to wild house finches (Carpodacus mexicanus), a common North American songbird. For the first time we characterize the genome and measure rates of genome evolution in House Finch isolates of MG, as well as in poultry outgroups. Using whole-genome sequences of 12 House Finch isolates across a 13-year serial sample and an additional four newly sequenced poultry strains, we estimate a nucleotide diversity in House Finch isolates of only ∼2% of ancestral poultry strains and a nucleotide substitution rate of 0.8-1.2×10(-5) per site per year both in poultry and in House Finches, an exceptionally fast rate rivaling some of the highest estimates reported thus far for bacteria. We also found high diversity and complete turnover of CRISPR arrays in poultry MG strains prior to the switch to the House Finch host, but after the invasion of House Finches there is progressive loss of CRISPR repeat diversity, and recruitment of novel CRISPR repeats ceases. Recent (2007) House Finch MG strains retain only ∼50% of the CRISPR repertoire founding (1994-95) strains and have lost the CRISPR-associated genes required for CRISPR function. Our results suggest that genome evolution in bacterial pathogens of wild birds can be extremely rapid and in this case is accompanied by apparent functional loss of CRISPRs.

Evidence of inter-component recombination, intra-component recombination and reassortment in banana bunchy top virus

Banana bunchy top virus (BBTV; family Nanoviridae, genus Babuvirus) is a multi-component, ssDNA virus, which causes widespread banana crop losses throughout tropical Africa and Australasia. We determined the full genome sequences of 12 BBTV isolates from the Kingdom of Tonga and analysed these together with previously determined BBTV sequences to show that reassortment and both inter- and intra-component recombination have all been relatively frequent occurrences during BBTV evolution. We found that whereas DNA-U3 components display evidence of complex inter- and intra-component recombination, all of the South Pacific DNA-R components have a common intra-component recombinant origin spanning the replication-associated protein gene. Altogether, the DNA-U3 and DNA-M components display a greater degree of inter-component recombination than the DNA-R, -S, -C and -M components. The breakpoint distribution of the inter-component recombination events reveals a primary recombination hotspot around the 5' side of the common region major and, in accordance with recombination hotspots detectable in related ssDNA viruses, a secondary recombination hotspot near the origin of virion-strand replication.

Recombination hotspots and host susceptibility modulate the adaptive value of recombination during maize streak virus evolution

BACKGROUND

Maize streak virus -strain A (MSV-A; Genus Mastrevirus, Family Geminiviridae), the maize-adapted strain of MSV that causes maize streak disease throughout sub-Saharan Africa, probably arose between 100 and 200 years ago via homologous recombination between two MSV strains adapted to wild grasses. MSV recombination experiments and analyses of natural MSV recombination patterns have revealed that this recombination event entailed the exchange of the movement protein - coat protein gene cassette, bounded by the two genomic regions most prone to recombination in mastrevirus genomes; the first surrounding the virion-strand origin of replication, and the second around the interface between the coat protein gene and the short intergenic region. Therefore, aside from the likely adaptive advantages presented by a modular exchange of this cassette, these specific breakpoints may have been largely predetermined by the underlying mechanisms of mastrevirus recombination. To investigate this hypothesis, we constructed artificial, low-fitness, reciprocal chimaeric MSV genomes using alternating genomic segments from two MSV strains; a grass-adapted MSV-B, and a maize-adapted MSV-A. Between them, each pair of reciprocal chimaeric genomes represented all of the genetic material required to reconstruct - via recombination - the highly maize-adapted MSV-A genotype, MSV-MatA. We then co-infected a selection of differentially MSV-resistant maize genotypes with pairs of reciprocal chimaeras to determine the efficiency with which recombination would give rise to high-fitness progeny genomes resembling MSV-MatA.

RESULTS

Recombinants resembling MSV-MatA invariably arose in all of our experiments. However, the accuracy and efficiency with which the MSV-MatA genotype was recovered across all replicates of each experiment depended on the MSV susceptibility of the maize genotypes used and the precise positions - in relation to known recombination hotspots - of the breakpoints required to re-create MSV-MatA. Although the MSV-sensitive maize genotype gave rise to the greatest variety of recombinants, the measured fitness of each of these recombinants correlated with their similarity to MSV-MatA.

CONCLUSIONS

The mechanistic predispositions of different MSV genomic regions to recombination can strongly influence the accessibility of high-fitness MSV recombinants. The frequency with which the fittest recombinant MSV genomes arise also correlates directly with the escalating selection pressures imposed by increasingly MSV-resistant maize hosts.

The evolutionary rate of citrus tristeza virus ranks among the rates of the slowest RNA viruses

Citrus tristeza virus (CTV) has been studied intensively at the molecular level. However, knowledge regarding the dynamics of its evolution is practically non-existent. In the past, diverse authors have referred to CTV as a highly variable virus, implying rapid evolution. Others have, in recent times, referred to CTV as an exceptionally slowly evolving virus. In this work, we used the capsid protein (CP) gene to estimate the rate of evolution. This was obtained from a large set of heterochronous CP gene sequences using a bayesian coalescent approach. The best-fitting evolutionary and population models pointed to an evolutionary rate of 1.58×10(-4) nt per site year(-1) (95 % highest posterior density, 1.73×10(-5)-3.16×10(-4) nt per site year(-1)). For an unbiased comparison with other plant and animal viruses, the evolutionary rate of synonymous substitutions was considered. In a series of 88 synonymous evolutionary rates, ranging from 5.2×10(-6) to 6.2×10(-2) nt per site year(-1), CTV ranks in the 10th percentile, embedded among the slowest animal RNA viruses. At the time of citrus dissemination to Europe and the New World, the major clades that led to the current phylogenetic groups were already defined, which may explain the absence nowadays of geographical speciation.