Experimental observations of rapid Maize streak virus evolution reveal a strand-specific nucleotide substitution bias

Twenty years of evolution and diversification of digitaria streak virus in Digitaria setigera

Within the family Geminiviridae, the emergence of new species results from their high mutation and recombination rates. In this study, we report the variability and evolution of digitaria streak virus (DSV), a mastrevirus isolated in 1986 from the grass Digitaria setigera in an island of the Vanuatu archipelago. Viral DNA of DSV samples was amplified from D. setigera specimens, derived from the naturally infected original plant, which were propagated in different laboratories in France and Italy for more than 20 years. From the consensus sequences, the nucleotide substitution rate was estimated for the period between a sample and the original sequence published in 1987, as well as for the period between samples. In addition, the intra-host genetic complexity and diversity of 8 DSV populations with a total of 165 sequenced haplotypes was characterized. The evolutionary rate of DSV was estimated to be between 1.13 × 10⁻⁴ and 9.87 × 10⁻⁴ substitutions/site/year, within the ranges observed in other single-stranded DNA viruses and RNA viruses. Bioinformatic analyses revealed high variability and heterogeneity in DSV populations, which confirmed that mutant spectra are continuously generated and are organized as quasispecies. The analysis of polymorphisms revealed nucleotide substitution biases in viral genomes towards deamination and oxidation of single-stranded DNA. The differences in variability in each of the genomic regions reflected a dynamic and modular evolution in the mutant spectra that was not reflected in the consensus sequences. Strikingly, the most variable region of the DSV genome, encoding the movement protein, showed rapid fixation of the mutations in the consensus sequence and a concomitant dN/dS ratio of 6.130, which suggests strong positive selection in this region. Phylogenetic analyses revealed a possible divergence in three genetic lineages from the original Vanuatu DSV isolate.

Identification of the Viral Determinant of Hypovirulence and Host Range in Sclerotiniaceae of a Genomovirus Reconstructed from the Plant Metagenome

Uncharacterized viral genomes that encode circular replication-associated proteins of single-stranded DNA viruses have been discovered by metagenomics/metatranscriptomics approaches. Some of these novel viruses are classified in the newly formed family Genomoviridae. Here, we determined the host range of a novel genomovirus, SlaGemV-1, through the transfection of Sclerotinia sclerotiorum with infectious clones. Inoculating with the rescued virions, we further transfected Botrytis cinerea and Monilinia fructicola, two economically important members of the family Sclerotiniaceae, and Fusarium oxysporum. SlaGemV-1 causes hypovirulence in S. sclerotiorum, B. cinerea, and M. fructicola. SlaGemV-1 also replicates in Spodoptera frugiperda insect cells but not in Caenorhabditis elegans or plants. By expressing viral genes separately through site-specific integration, the replication protein alone was sufficient to cause debilitation. Our study is the first to demonstrate the reconstruction of a metagenomically discovered genomovirus without known hosts with the potential of inducing hypovirulence, and the infectious clone allows for studying mechanisms of genomovirus-host interactions that are conserved across genera. IMPORTANCE Little is known about the exact host range of widespread genomoviruses. The genome of soybean leaf-associated gemygorvirus-1 (SlaGemV-1) was originally assembled from a metagenomic/metatranscriptomic study without known hosts. Here, we rescued SlaGemV-1 and found that it could infect three important plant-pathogenic fungi and fall armyworm (S. frugiperda Sf9) insect cells but not a model nematode, C. elegans, or model plant species. Most importantly, SlaGemV-1 shows promise for inducing hypovirulence of the tested fungal species in the family Sclerotiniaceae, including Sclerotinia sclerotiorum, Botrytis cinerea, and Monilinia fructicola. The viral determinant of hypovirulence was further identified as replication initiation protein. As a proof of concept, we demonstrate that viromes discovered in plant metagenomes can be a valuable genetic resource when novel viruses are rescued and characterized for their host range.

Population diversity of cassava mosaic begomoviruses increases over the course of serial vegetative propagation

Cassava mosaic disease (CMD) represents a serious threat to cassava, a major root crop for more than 300 million Africans. CMD is caused by single-stranded DNA begomoviruses that evolve rapidly, making it challenging to develop durable disease resistance. In addition to the evolutionary forces of mutation, recombination and reassortment, factors such as climate, agriculture practices and the presence of DNA satellites may impact viral diversity. To gain insight into the factors that alter and shape viral diversity in planta, we used high-throughput sequencing to characterize the accumulation of nucleotide diversity after inoculation of infectious clones corresponding to African cassava mosaic virus (ACMV) and East African cassava mosaic Cameroon virus (EACMCV) in the susceptible cassava landrace Kibandameno. We found that vegetative propagation had a significant effect on viral nucleotide diversity, while temperature and a satellite DNA did not have measurable impacts in our study. EACMCV diversity increased linearly with the number of vegetative propagation passages, while ACMV diversity increased for a time and then decreased in later passages. We observed a substitution bias toward C→T and G→A for mutations in the viral genomes consistent with field isolates. Non-coding regions excluding the promoter regions of genes showed the highest levels of nucleotide diversity for each genome component. Changes in the 5' intergenic region of DNA-A resembled the sequence of the cognate DNA-B sequence. The majority of nucleotide changes in coding regions were non-synonymous, most with predicted deleterious effects on protein structure, indicative of relaxed selection pressure over six vegetative passages. Overall, these results underscore the importance of knowing how cropping practices affect viral evolution and disease progression.

A Novel SARS-CoV-2 Viral Sequence Bioinformatic Pipeline Has Found Genetic Evidence That the Viral 3' Untranslated Region (UTR) Is Evolving and Generating Increased Viral Diversity

An unprecedented amount of SARS-CoV-2 sequencing has been performed, however, novel bioinformatic tools to cope with and process these large datasets is needed. Here, we have devised a bioinformatic pipeline that inputs SARS-CoV-2 genome sequencing in FASTA/FASTQ format and outputs a single Variant Calling Format file that can be processed to obtain variant annotations and perform downstream population genetic testing. As proof of concept, we have analyzed over 229,000 SARS-CoV-2 viral sequences up until November 30, 2020. We have identified over 39,000 variants worldwide with increased polymorphisms, spanning the ORF3a gene as well as the 3' untranslated (UTR) regions, specifically in the conserved stem loop region of SARS-CoV-2 which is accumulating greater observed viral diversity relative to chance variation. Our analysis pipeline has also discovered the existence of SARS-CoV-2 hypermutation with low frequency (less than in 2% of genomes) likely arising through host immune responses and not due to sequencing errors. Among annotated non-sense variants with a population frequency over 1%, recurrent inactivation of the ORF8 gene was found. This was found to be present in the newly identified B.1.1.7 SARS-CoV-2 lineage that originated in the United Kingdom. Almost all VOC-containing genomes possess one stop codon in ORF8 gene (Q27∗), however, 13% of these genomes also contains another stop codon (K68∗), suggesting that ORF8 loss does not interfere with SARS-CoV-2 spread and may play a role in its increased virulence. We have developed this computational pipeline to assist researchers in the rapid analysis and characterization of SARS-CoV-2 variation.

Molecular phylogenetics and evolutionary analysis of a highly recombinant begomovirus, Cotton leaf curl Multan virus, and associated satellites

Cotton leaf curl Multan virus (CLCuMuV) and its associated satellites are a major part of the cotton leaf curl disease (CLCuD) caused by the begomovirus species complex. Despite the implementation of potential disease management strategies, the incessant resurgence of resistance-breaking variants of CLCuMuV imposes a continuous threat to cotton production. Here, we present a focused effort to map the geographical prevalence, genomic diversity, and molecular evolutionary endpoints that enhance disease complexity by facilitating the successful adaptation of CLCuMuV populations to the diversified ecosystems. Our results demonstrate that CLCuMuV populations are predominantly distributed in China, while the majority of alphasatellites and betasatellites exist in Pakistan. We demonstrate that together with frequent recombination, an uneven genetic variation mainly drives CLCuMuV and its satellite's virulence and evolvability. However, the pattern and distribution of recombination breakpoints greatly vary among viral and satellite sequences. The CLCuMuV, Cotton leaf curl Multan alphasatellite, and Cotton leaf curl Multan betasatellite populations arising from distinct regions exhibit high mutation rates. Although evolutionarily linked, these populations are independently evolving under strong purifying selection. These findings will facilitate to comprehensively understand the standing genetic variability and evolutionary patterns existing among CLCuMuV populations across major cotton-producing regions of the world.

Intra-host evolution of the ssDNA virus tomato severe rugose virus (ToSRV)

To evaluate and quantify the evolutionary dynamics of the bipartite begomovirus tomato severe rugose virus (ToSRV) in a cultivated and a non-cultivated host, plants of tomato and Nicandra physaloides were biolistically inoculated with an infectious clone and systemically infected leaves were sampled at 30, 75 and 120 days after inoculation. Total DNA was extracted and sequenced in the Illumina HiSeq 2000 platform. The datasets were trimmed with the quality score limit set to 0.01, and the assembly was performed using the infectious clone sequence as reference. SNPs were filtered using a minimum p-value of 0.001 and the sum frequencies were used to calculate the deviation from the original clone sequence. Nucleotide substitution rates were calculated for the two DNA components in both hosts: 1.73 × 10^-3 and 3.07 × 10^-4 sub/site/year for the DNA-A and DNA-B, respectively, in N. physaloides, and 8.05 × 10^-4 and 7.02 × 10^-5 sub/site/year the for DNA-A and DNA-B, respectively, in tomato. These values are in the same range of those estimated for viruses with single-stranded RNA genomes and for other begomoviruses. Strikingly, the number of substitutions decreased over time, suggesting the presence of bottlenecks during systemic infection. Determination of Shannon's entropy indicated different patterns of variation in the DNA-A and the DNA-B, suggesting distinct evolutionary forces acting upon each component.

Geographic Distribution and Genetic Diversity of Rice Stripe Mosaic Virus in Southern China

Rice stripe mosaic virus (RSMV) transmitted by the leafhopper Recilia dorsalis is a tentative new species in the genus Cytorhabdovirus identified recently in South China. To explore its geographic distribution and genetic diversity, field investigation and viral whole-genome sequencing were conducted in this study. The results indicated that RSMV was present in the rice samples collected across southern China. Twelve representative samples from different geographical regions were selected for viral whole-genome sequencing and the viral genome variation was analyzed in combination with a previously reported RSMV isolate. Identity analysis showed that the genome sequences of 13 RSMV isolates were highly conserved with nucleotide identities over 99.4%. There was a strong negative selection pressure during the evolution of RSMV with more transitions (72.08%) than transversions (27.92%) found between the RSMV isolates. Among the seven genes encoded by RSMV, the P gene was the most variable, followed by N, M, L, and G; the P3 and P6 amino acid sequences were not found to be mutated and no mutations were found in the non-coding region. A phylogenetic tree based on the RSMV whole-genome nucleotide sequences revealed that all RSMV isolates clustered in two groups based on geographical origin. Notably, the L proteins of the Guangxi and Hainan isolates had five and one specific amino acid sites, respectively, suggesting that the L gene has undergone environmental adaptive variation during the dispersal of RSMV.

Eukaryotic Circular Rep-Encoding Single-Stranded DNA (CRESS DNA) Viruses: Ubiquitous Viruses With Small Genomes and a Diverse Host Range

While single-stranded DNA (ssDNA) was once thought to be a relatively rare genomic architecture for viruses, modern metagenomics sequencing has revealed circular ssDNA viruses in most environments and in association with diverse hosts. In particular, circular ssDNA viruses encoding a homologous replication-associated protein (Rep) have been identified in the majority of eukaryotic supergroups, generating interest in the ecological effects and evolutionary history of circular Rep-encoding ssDNA viruses (CRESS DNA) viruses. This review surveys the explosion of sequence diversity and expansion of eukaryotic CRESS DNA taxonomic groups over the last decade, highlights similarities between the well-studied geminiviruses and circoviruses with newly identified groups known only through their genome sequences, discusses the ecology and evolution of eukaryotic CRESS DNA viruses, and speculates on future research horizons.

The diversity, evolution and epidemiology of plant viruses: A phylogenetic view

During the past four decades, the scientific community has seen an exponential advance in the number, sophistication, and quality of molecular techniques and bioinformatics tools for the genetic characterization of plant virus populations. Predating these advances, the field of Phylogenetics has significantly contributed to understand important aspects of plant virus evolution. This review aims at summarizing the impact of Phylogenetics in the current knowledge on three major aspects of plant virus evolution that have benefited from the development of phylogenetic inference: (1) The identification and classification of plant virus diversity. (2) The mechanisms and forces shaping the evolution of plant virus populations. (3) The understanding of the interaction between plant virus evolution, epidemiology and ecology. The work discussed here highlights the important role of phylogenetic approaches in the study of the dynamics of plant virus populations.

Differential Shape of Geminivirus Mutant Spectra Across Cultivated and Wild Hosts With Invariant Viral Consensus Sequences

Geminiviruses (family Geminiviridae) possess single-stranded circular DNA genomes that are replicated by cellular polymerases in plant host cell nuclei. In their hosts, geminivirus populations behave as ensembles of mutant and recombinant genomes, known as viral quasispecies. This favors the emergence of new geminiviruses with altered host range, facilitating new or more severe diseases or overcoming resistance traits. In warm and temperate areas several whitefly-transmitted geminiviruses of the genus Begomovirus cause the tomato yellow leaf curl disease (TYLCD) with significant economic consequences. TYLCD is frequently controlled in commercial tomatoes by using the dominant Ty-1 resistance gene. Over a 45 day period we have studied the diversification of three begomoviruses causing TYLCD: tomato yellow leaf curl virus (TYLCV), tomato yellow leaf curl Sardinia virus (TYLCSV) and tomato yellow leaf curl Malaga virus (TYLCMaV, a natural recombinant between TYLCV and TYLCSV). Viral quasispecies resulting from inoculation of geminivirus infectious clones were examined in plants of susceptible tomato (ty-1/ty-1), heterozygous resistant tomato (Ty-1/ty-1), common bean, and the wild reservoir Solanum nigrum. Differences in virus fitness across hosts were observed while viral consensus sequences remained invariant. However, the complexity and heterogeneity of the quasispecies were high, especially in common bean and the wild host. Interestingly, the presence or absence of the Ty-1 allele in tomato did not lead to differences in begomovirus mutant spectra. However, the fitness decrease of TYLCSV and TYLCV in tomato at 45 dpi might be related to an increase in CP (Coat protein) mutation frequency. In Solanum nigrum the recombinant TYLCMaV, which showed lower fitness than TYLCSV, at 45 dpi actively explored Rep (Replication associated protein) ORF but not the overlapping C4. Our results underline the importance of begomovirus mutant spectra during infections. This is especially relevant in the wild reservoir of the viruses, which has the potential to maintain highly diverse mutant spectra without modifying their consensus sequences.

Possible risks posed by single-stranded DNA viruses of pigs associated with xenotransplantation

Routine large-scale xenotransplantation from pigs to humans is getting closer to clinical reality owing to several state-of-the-art technologies, especially the ability to rapidly engineer genetically defined pigs. However, using pig organs in humans poses risks including unwanted cross-species transfer of viruses and adaption of these pig viruses to the human organ recipient. Recent developments in the field of virology, including the advent of metagenomic techniques to characterize entire viromes, have led to the identification of a plethora of viruses in many niches. Single-stranded DNA (ssDNA) viruses are the largest group prevalent in virome studies in mammals. Specifically, the ssDNA viral genomes are characterized by a high rate of nucleotide substitution, which confers a proclivity to adapt to new hosts and cross-species barriers. Pig-associated ssDNA viruses include torque teno sus viruses (TTSuV) in the Anelloviridae family, porcine parvoviruses (PPV), and porcine bocaviruses (PBoV) both in the family of Parvoviridae, and porcine circoviruses (PCV) in the Circoviridae family, some of which have been confirmed to be pathogenic to pigs. The risks of these viruses for the human recipient during xenotransplantation procedures are relatively unknown. Based on the scant knowledge available on the prevalence, predilection, and pathogenicity of pig-associated ssDNA viruses, careful screening and monitoring are required. In the case of positive identification, risk assessments and strategies to eliminate these viruses in xenotransplantation pig stock may be needed.

Evolution of a wild-plant tobamovirus passaged through an exotic host: Fixation of mutations and increased replication

Tobamovirus is a group of viruses that have become serious pathogens of crop plants. As part of a study informing risk of wild plant virus spill over to crops, we investigated the capacity of a solanaceous-infecting tobamovirus from an isolated indigenous flora to adapt to new exotic hosts. Yellow tailflower mild mottle virus (YTMMV) (genus Tobamovirus, family Virgaviridae) was isolated from a wild plant of yellow tailflower (Anthocercis littoria, family Solanaceae) and initially passaged through a plant of Nicotiana benthamiana, then one of Nicotiana glutinosa where a single local lesion was used to inoculate a N. benthamiana plant. Sap from this plant was used as starting material for nine serial passages through three plant species. The virus titre was recorded periodically, and 85% of the virus genome was sequenced at each passage for each host. Six polymorphic sites were found in the YTMMV genome across all hosts and passages. At five of these, the alternate alleles became fixed in the viral genome until the end of the experiment. Of these five alleles, one was a non-synonymous mutation (U1499C) that occurred only when the virus replicated in tomato. The mutant isolate harbouring U1499C, designated YTMMV-δ, increased its titre over passages in tomato and outcompeted the wild-type isolate when both were co-inoculated to tomato. That YTMMV-δ had greater reproductive fitness in an exotic host than did the wild type isolate suggests YTMMV evolution is influenced by host changes.

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.

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.

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.

The role of corchorus in spreading of tomato yellow leaf curl virus on tomato in Jeddah, Saudi Arabia

Corchorus (Corchorus capsularis L. and Corchorus olitorius L.) is one of the most important fiber crops grown in tropical and subtropical regions throughout the world. Field survey was conducted and naturally infected leaf samples were collected from corchorus and tomato plants in Jeddah, Saudi Arabia. The causal virus was transmitted by whiteflies to tomato plants and begomovirus infection was confirmed by Polymerase chain reaction. The complete viral genome and associated betasatellites were amplified, cloned and sequenced from both corchorus and tomato samples. The genetic variability and phylogenetic relationships were determined for both isolates (corchorus and tomato). The complete genome sequences showed highest (99.5 % nt) similarity with tomato yellow leaf curl virus (TYLCV) and formed closest cluster with TYLCV-Tomato reported from Jizan and Al-Qasim, Saudi Arabia and betasatellites sequences showed highest similarity (99.8 % nt) with Tomato yellow leaf curl betasatellites-Jeddah followed by Tomato yellow leaf curl Oman betasatellites and formed closed cluster with TYLCV-Tomato. On the basis of results obtained from whiteflies transmission, sequence similarity and phylogenetic relationships; it is concluded that the identified virus could be a variant of TYLCV circulating in the Kingdom. The significance of this study demonstrated that the corchorus is serving as reservoir and alternative host and playing an important role in spreading the begomovirus associated disease in the Kingdom of Saudi Arabia.

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.

The global distribution of Banana bunchy top virus reveals little evidence for frequent recent, human-mediated long distance dispersal events

Banana bunchy top virus (BBTV; family Nanoviridae, genus Babuvirus) is a multi-component single-stranded DNA virus, which infects banana plants in many regions of the world, often resulting in large-scale crop losses. We analyzed 171 banana leaf samples from fourteen countries and recovered, cloned, and sequenced 855 complete BBTV components including ninety-four full genomes. Importantly, full genomes were determined from eight countries, where previously no full genomes were available (Samoa, Burundi, Republic of Congo, Democratic Republic of Congo, Egypt, Indonesia, the Philippines, and the USA [HI]). Accounting for recombination and genome component reassortment, we examined the geographic structuring of global BBTV populations to reveal that BBTV likely originated in Southeast Asia, that the current global hotspots of BBTV diversity are Southeast Asia/Far East and India, and that BBTV populations circulating elsewhere in the world have all potentially originated from infrequent introductions. Most importantly, we find that rather than the current global BBTV distribution being due to increases in human-mediated movements of bananas over the past few decades, it is more consistent with a pattern of infrequent introductions of the virus to different parts of the world over the past 1,000 years.

Single-stranded DNA viruses employ a variety of mechanisms for integration into host genomes

Single-stranded DNA (ssDNA) viruses are widespread in the environment and include economically, medically, and ecologically important pathogens. Recently, it has been discovered that ssDNA virus genomes are also prevalent in the chromosomes of their bacterial, archaeal, and eukaryotic hosts. Sequences originating from viruses of the families Parvoviridae, Circoviridae, and Geminiviridae are particularly widespread in the genomes of eukaryotes, where they are often fossilized as endogenous viral elements. ssDNA viruses have evolved diverse mechanisms to invade cellular genomes, and these principally vary between viruses infecting bacteria/archaea and eukaryotes. Filamentous bacteriophages (Inoviridae) use at least three major mechanisms of integration. Some of these phages encode integrases of serine or tyrosine recombinase superfamilies, while others utilize DDE transposases of the IS3, IS30, or IS110/IS492 families, whereas some inoviruses, and possibly certain members of the Microviridae, hijack the host XerCD recombination machinery. By contrast, eukaryotic viruses for integration rely on the endonuclease activity of their rolling-circle replication-initiation proteins, mimicking the mechanisms used by some bacterial transposons. Certain bacterial and eukaryotic ssDNA viruses have embraced a transposon-like means of propagation, with occasionally dramatic effects on host genome evolution. Here, we review the diversity of experimentally verified and hypothetical mechanisms of genome integration employed by ssDNA viruses, and consider the evolutionary implications of these processes, particularly in the emergence of novel virus groups.

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.

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.

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.

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.

Evolutionary time-scale of the begomoviruses: evidence from integrated sequences in the Nicotiana genome

Despite having single stranded DNA genomes that are replicated by host DNA polymerases, viruses in the family Geminiviridae are apparently evolving as rapidly as some RNA viruses. The observed substitution rates of geminiviruses in the genera Begomovirus and Mastrevirus are so high that the entire family could conceivably have originated less than a million years ago (MYA). However, the existence of geminivirus related DNA (GRD) integrated within the genomes of various Nicotiana species suggests that the geminiviruses probably originated >10 MYA. Some have even suggested that a distinct New-World (NW) lineage of begomoviruses may have arisen following the separation by continental drift of African and American proto-begomoviruses ∼110 MYA. We evaluate these various geminivirus origin hypotheses using Bayesian coalescent-based approaches to date firstly the Nicotiana GRD integration events, and then the divergence of the NW and Old-World (OW) begomoviruses. Besides rejecting the possibility of a<2 MYA OW-NW begomovirus split, we could also discount that it may have occurred concomitantly with the breakup of Gondwanaland 110 MYA. Although we could only confidently narrow the date of the split down to between 2 and 80 MYA, the most plausible (and best supported) date for the split is between 20 and 30 MYA--a time when global cooling ended the dispersal of temperate species between Asia and North America via the Beringian land bridge.

Genetic variability of East African cassava mosaic Cameroon virus under field and controlled environment conditions

Cassava geminiviruses occur in all cassava growing areas of Africa and are considered to be the most damaging vector-borne plant pathogens. At least seven species of these viruses have been identified. We investigated genetic variation in East African cassava mosaic cassava Cameroon virus (EACMCV) from naturally infected cassava and from experimentally infected Nicotiana benthamiana. Results showed that the populations of EACMCV in cassava and in N. benthamiana were genetically heterogeneous. Mutation frequencies in the order of 10(-4), comparable to that reported for plant RNA viruses, were observed in both hosts. We also produced an EACMCV mutant that induces reversion and second site mutations, thus suggesting that a high mutation frequency facilitates the maintenance of genome structure and function. This is direct experimental evidence showing that cassava geminiviruses exhibit a high mutation frequency and that a single clone quickly transforms into a collection of mutant sequences upon introduction into the host.

Global genetic diversity and geographical and host-species distribution of beak and feather disease virus isolates

Psittacine beak and feather disease (PBFD) has a broad host range and is widespread in wild and captive psittacine populations in Asia, Africa, the Americas, Europe and Australasia. Beak and feather disease circovirus (BFDV) is the causative agent. BFDV has an ∼2 kb single stranded circular DNA genome encoding just two proteins (Rep and CP). In this study we provide support for demarcation of BFDV strains by phylogenetic analysis of 65 complete genomes from databases and 22 new BFDV sequences isolated from infected psittacines in South Africa. We propose 94% genome-wide sequence identity as a strain demarcation threshold, with isolates sharing >94% identity belonging to the same strain, and strain subtypes sharing >98% identity. Currently, BFDV diversity falls within 14 strains, with five highly divergent isolates from budgerigars probably representing a new species of circovirus with three strains (budgerigar circovirus; BCV-A, -B and -C). The geographical distribution of BFDV and BCV strains is strongly linked to the international trade in exotic birds; strains with more than one host are generally located in the same geographical area. Lastly, we examined BFDV and BCV sequences for evidence of recombination, and determined that recombination had occurred in most BFDV and BCV strains. We established that there were two globally significant recombination hotspots in the viral genome: the first is along the entire intergenic region and the second is in the C-terminal portion of the CP ORF. The implications of our results for the taxonomy and classification of circoviruses are discussed.

High variability and rapid evolution of a nanovirus

Nanoviruses are multipartite single-stranded DNA (ssDNA) plant viruses that cause important diseases of leguminous crops and banana. Little has been known about the variability and molecular evolution of these viruses. Here we report on the variability of faba bean necrotic stunt virus (FBNSV), a nanovirus from Ethiopia. We found mutation frequencies of 7.52 x 10(-4) substitutions per nucleotide in a field population of the virus and 5.07 x 10(-4) substitutions per nucleotide in a laboratory-maintained population derived thereof. Based on virus propagation for a period of more than 2 years, we determined a nucleotide substitution rate of 1.78 x 10(-3) substitutions per nucleotide per year. This high molecular evolution rate places FBNSV, as a representative of the family Nanoviridae, among the fastest-evolving ssDNA viruses infecting plants or vertebrates.

Maize streak virus: an old and complex 'emerging' pathogen

Maize streak virus (MSV; Genus Mastrevirus, Family Geminiviridae) occurs throughout Africa, where it causes what is probably the most serious viral crop disease on the continent. It is obligately transmitted by as many as six leafhopper species in the Genus Cicadulina, but mainly by C. mbila Naudé and C. storeyi. In addition to maize, it can infect over 80 other species in the Family Poaceae. Whereas 11 strains of MSV are currently known, only the MSV-A strain is known to cause economically significant streak disease in maize. Severe maize streak disease (MSD) manifests as pronounced, continuous parallel chlorotic streaks on leaves, with severe stunting of the affected plant and, usuallly, a failure to produce complete cobs or seed. Natural resistance to MSV in maize, and/or maize infections caused by non-maize-adapted MSV strains, can result in narrow, interrupted streaks and no obvious yield losses. MSV epidemiology is primarily governed by environmental influences on its vector species, resulting in erratic epidemics every 3-10 years. Even in epidemic years, disease incidences can vary from a few infected plants per field, with little associated yield loss, to 100% infection rates and complete yield loss.

TAXONOMY

The only virus species known to cause MSD is MSV, the type member of the Genus Mastrevirus in the Family Geminiviridae. In addition to the MSV-A strain, which causes the most severe form of streak disease in maize, 10 other MSV strains (MSV-B to MSV-K) are known to infect barley, wheat, oats, rye, sugarcane, millet and many wild, mostly annual, grass species. Seven other mastrevirus species, many with host and geographical ranges partially overlapping those of MSV, appear to infect primarily perennial grasses.

PHYSICAL PROPERTIES

MSV and all related grass mastreviruses have single-component, circular, single-stranded DNA genomes of approximately 2700 bases, encapsidated in 22 x 38-nm geminate particles comprising two incomplete T = 1 icosahedra, with 22 pentameric capsomers composed of a single 32-kDa capsid protein. Particles are generally stable in buffers of pH 4-8.

DISEASE SYMPTOMS

In infected maize plants, streak disease initially manifests as minute, pale, circular spots on the lowest exposed portion of the youngest leaves. The only leaves that develop symptoms are those formed after infection, with older leaves remaining healthy. As the disease progresses, newer leaves emerge containing streaks up to several millimetres in length along the leaf veins, with primary veins being less affected than secondary or tertiary veins. The streaks are often fused laterally, appearing as narrow, broken, chlorotic stripes, which may extend over the entire length of severely affected leaves. Lesion colour generally varies from white to yellow, with some virus strains causing red pigmentation on maize leaves and abnormal shoot and flower bunching in grasses. Reduced photosynthesis and increased respiration usually lead to a reduction in leaf length and plant height; thus, maize plants infected at an early stage become severely stunted, producing undersized, misshapen cobs or giving no yield at all. Yield loss in susceptible maize is directly related to the time of infection: infected seedlings produce no yield or are killed, whereas plants infected at later times are proportionately less affected.

DISEASE CONTROL

Disease avoidance can be practised by only planting maize during the early season when viral inoculum loads are lowest. Leafhopper vectors can also be controlled with insecticides such as carbofuran. However, the development and use of streak-resistant cultivars is probably the most effective and economically viable means of preventing streak epidemics. Naturally occurring tolerance to MSV (meaning that, although plants become systemically infected, they do not suffer serious yield losses) has been found, which has primarily been attributed to a single gene, msv-1. However, other MSV resistance genes also exist and improved resistance has been achieved by concentrating these within individual maize genotypes. Whereas true MSV immunity (meaning that plants cannot be symptomatically infected by the virus) has been achieved in lines that include multiple small-effect resistance genes together with msv-1, it has proven difficult to transfer this immunity into commercial maize genotypes. An alternative resistance strategy using genetic engineering is currently being investigated in South Africa.

USEFUL WEBSITES

http://www.mcb.uct.ac.za/MSV/mastrevirus.htm; http://www.danforthcenter.org/iltab/geminiviridae/geminiaccess/mastrevirus/Mastrevirus.htm.

Recombination, selection and clock-like evolution of Rice yellow mottle virus

The clock-like diversification of Rice yellow mottle virus (RYMV), a widespread RNA plant virus that infects rice in Africa, was tested following a three-step approach with (i) an exhaustive search of recombinants, (ii) a comprehensive assessment of the selective constraints over lineages, and (iii) a stepwise series of tests of the molecular clock hypothesis. The first evidence of recombination in RYMV was found in East Africa, in the region most favorable to co-infection. RYMV evolved under a pronounced purifying selection, but the selection pressure did vary among lineages. There was no phylogenetic evidence of transient deleterious mutations. ORF2b, which codes for the polymerase and is the most constrained ORF, tends to diversify clock-like. With the other ORFs and the full genome, the departure from the strict clock model was limited. This likely reflects the dominant conservative selection pressure and the clock-like fixation of synonymous mutations.

Dating the origins of the maize-adapted strain of maize streak virus, MSV-A

Maize streak virus (MSV), which causes maize streak disease (MSD), is one of the most serious biotic threats to African food security. Here, we use whole MSV genomes sampled over 30 years to estimate the dates of key evolutionary events in the 500 year association of MSV and maize. The substitution rates implied by our analyses agree closely with those estimated previously in controlled MSV evolution experiments, and we use them to infer the date when the maize-adapted strain, MSV-A, was generated by recombination between two grass-adapted MSV strains. Our results indicate that this recombination event occurred in the mid-1800s, ∼20 years before the first credible reports of MSD in South Africa and centuries after the introduction of maize to the continent in the early 1500s. This suggests a causal link between MSV recombination and the emergence of MSV-A as a serious pathogen of maize.

Experimental evidence indicating that mastreviruses probably did not co-diverge with their hosts

Background

Despite the demonstration that geminiviruses, like many other single stranded DNA viruses, are evolving at rates similar to those of RNA viruses, a recent study has suggested that grass-infecting species in the genus Mastrevirus may have co-diverged with their hosts over millions of years. This "co-divergence hypothesis" requires that long-term mastrevirus substitution rates be at least 100,000-fold lower than their basal mutation rates and 10,000-fold lower than their observable short-term substitution rates. The credibility of this hypothesis, therefore, hinges on the testable claim that negative selection during mastrevirus evolution is so potent that it effectively purges 99.999% of all mutations that occur.

Results

We have conducted long-term evolution experiments lasting between 6 and 32 years, where we have determined substitution rates of between 2 and 3 × 10^-4 substitutions/site/year for the mastreviruses Maize streak virus (MSV) and Sugarcane streak Réunion virus (SSRV). We further show that mutation biases are similar for different geminivirus genera, suggesting that mutational processes that drive high basal mutation rates are conserved across the family. Rather than displaying signs of extremely severe negative selection as implied by the co-divergence hypothesis, our evolution experiments indicate that MSV and SSRV are predominantly evolving under neutral genetic drift.

Conclusion

The absence of strong negative selection signals within our evolution experiments and the uniformly high geminivirus substitution rates that we and others have reported suggest that mastreviruses cannot have co-diverged with their hosts.

Diverse circovirus-like genome architectures revealed by environmental metagenomics

Single-stranded DNA (ssDNA) viruses with circular genomes are the smallest viruses known to infect eukaryotes. The present study identified 10 novel genomes similar to ssDNA circoviruses through data-mining of public viral metagenomes. The metagenomic libraries included samples from reclaimed water and three different marine environments (Chesapeake Bay, British Columbia coastal waters and Sargasso Sea). All the genomes have similarities to the replication (Rep) protein of circoviruses; however, only half have genomic features consistent with known circoviruses. Some of the genomes exhibit a mixture of genomic features associated with different families of ssDNA viruses (i.e. circoviruses, geminiviruses and parvoviruses). Unique genome architectures and phylogenetic analysis of the Rep protein suggest that these viruses belong to novel genera and/or families. Investigating the complex community of ssDNA viruses in the environment can lead to the discovery of divergent species and help elucidate evolutionary links between ssDNA viruses.

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

Whitefly-transmitted geminiviruses are major pathogens of the important crop cassava in Africa. The intensive sampling and sequencing of cassava mosaic disease-causing viruses that occurred in the wake of a severe outbreak in Central Africa (1997-2002) allowed us to estimate the rate of evolution of this virus. East African cassava mosaic virus and related species are obligately bipartite (DNA-A and DNA-B segments), and these two genome segments have different evolutionary histories. Despite these phylogenetic differences, we inferred high rates of nucleotide substitution in both segments: mean rates of 1.60x10(-3) and 1.33x10(-4) substitutions site(-1) year(-1) for DNA-A and DNA-B, respectively. While similarly high substitution rates were found in datasets free of detectable recombination, only that estimated for the coat protein gene (AV1), for which an additional DNA-A sequence isolated in 1995 was available, was statistically robust. These high substitution rates also confirm that those previously estimated for the monopartite tomato yellow leaf curl virus (TYLCV) are representative of multiple begomoviruses. We also validated our rate estimates by comparing them with those depicting the emergence of TYLCV in North America. These results further support the notion that geminiviruses evolve as rapidly as many RNA viruses.

Rapid host adaptation by extensive recombination

Experimental investigations into virus recombination can provide valuable insights into the biochemical mechanisms and the evolutionary value of this fundamental biological process. Here, we describe an experimental scheme for studying recombination that should be applicable to any recombinogenic viruses amenable to the production of synthetic infectious genomes. Our approach is based on differences in fitness that generally exist between synthetic chimaeric genomes and the wild-type viruses from which they are constructed. In mixed infections of defective reciprocal chimaeras, selection strongly favours recombinant progeny genomes that recover a portion of wild-type fitness. Characterizing these evolved progeny viruses can highlight both important genetic fitness determinants and the contribution that recombination makes to the evolution of their natural relatives. Moreover, these experiments supply precise information about the frequency and distribution of recombination breakpoints, which can shed light on the mechanistic processes underlying recombination. We demonstrate the value of this approach using the small single-stranded DNA geminivirus, maize streak virus (MSV). Our results show that adaptive recombination in this virus is extremely efficient and can yield complex progeny genomes comprising up to 18 recombination breakpoints. The patterns of recombination that we observe strongly imply that the mechanistic processes underlying rolling circle replication are the prime determinants of recombination breakpoint distributions found in MSV genomes sampled from nature.