Revisiting the origins of the Sobemovirus genus: A case for ancient origins of plant viruses.

The discrepancy between short- and long-term rate estimates, known as the time-dependent rate phenomenon (TDRP), poses a challenge to extrapolating evolutionary rates over time and reconstructing evolutionary history of viruses. The TDRP reveals a decline in evolutionary rate estimates with the measurement timescale, explained empirically by a power-law rate decay, notably observed in animal and human viruses. A mechanistic evolutionary model, the Prisoner of War (PoW) model, has been proposed to address TDRP in viruses. Although TDRP has been studied in animal viruses, its impact on plant virus evolutionary history remains largely unexplored. Here, we investigated the consequences of TDRP in plant viruses by applying the PoW model to reconstruct the evolutionary history of sobemoviruses, plant pathogens with significant importance due to their impact on agriculture and plant health. Our analysis showed that the Sobemovirus genus dates back over four million years, indicating an ancient origin. We found evidence that supports deep host jumps to Poaceae, Fabaceae, and Solanaceae occurring between tens to hundreds of thousand years ago, followed by specialization. Remarkably, the TDRP-corrected evolutionary history of sobemoviruses was extended far beyond previous estimates that had suggested their emergence nearly 9,000 years ago, a time coinciding with the Neolithic period in the Near East. By incorporating sequences collected through metagenomic analyses, the resulting phylogenetic tree showcases increased genetic diversity, reflecting a deep history of sobemovirus species. We identified major radiation events beginning between 4,600 to 2,000 years ago, which aligns with the Neolithic period in various regions, suggesting a period of rapid diversification from then to the present. Our findings make a case for the possibility of deep evolutionary origins of plant viruses.

EvoLaps: a web interface to visualize continuous phylogeographic reconstructions.

BACKGROUND: Phylogeographic reconstructions serve as a basis to understand the spread and evolution of pathogens. Visualization of these reconstructions often lead to complex graphical representations which are difficult to interpret. RESULT: We present EvoLaps, a user-friendly web interface to visualize phylogeographic reconstructions based on the analysis of latitude/longitude coordinates with various clustering levels. EvoLaps also produces transition diagrams that provide concise and easy to interpret summaries of phylogeographic reconstructions. CONCLUSION: The main contribution of EvoLaps is to assemble known numerical and graphical methods/tools into a user-friendly interface dedicated to the visualization and edition of evolutionary scenarios based on continuous phylogeographic reconstructions. EvoLaps is freely usable at www.evolaps.org .

ICTV Virus Taxonomy Profile: Solemoviridae 2021.

The family Solemoviridae includes viruses with icosahedral particles (26-34 nm in diameter) assembled on T=3 symmetry with a 4-6 kb positive-sense, monopartite, polycistronic RNA genome. Transmission of members of the genera Sobemovirus and Polemovirus occurs via mechanical wounding, vegetative propagation, insect vectors or abiotically through soil; members of the genera Polerovirus and Enamovirus are transmitted by specific aphids. Most solemoviruses have a narrow host range. This is a summary of the International Committee on Taxonomy of Viruses (ICTV) Report on the family Solemoviridae, which is available at ictv.global/report/solemoviridae.

Identification of a Hypervirulent Pathotype of Rice yellow mottle virus: A Threat to Genetic Resistance Deployment in West-Central Africa.

Rice yellow mottle virus (RYMV) causes high losses to rice production in Africa. Several sources of varietal high resistance are available but the emergence of virulent pathotypes that are able to overcome one or two resistance alleles can sometimes occur. Both resistance spectra and viral adaptability have to be taken into account to develop sustainable rice breeding strategies against RYMV. In this study, we extended previous resistance spectrum analyses by testing the rymv1-4 and rymv1-5 alleles that are carried by the rice accessions Tog5438 and Tog5674, respectively, against isolates that are representative of RYMV genetic and pathogenic diversity. Our study revealed a hypervirulent pathotype, named thereafter pathotype T', that is able to overcome all known sources of high resistance. This pathotype, which is spatially localized in West-Central Africa, appears to be more abundant than previously suspected. To better understand the adaptive processes of pathotype T', molecular determinants of resistance breakdown were identified via Sanger sequencing and validated through directed mutagenesis of an infectious clone. These analyses confirmed the key role of convergent nonsynonymous substitutions in the central part of the viral genome-linked protein to overcome RYMV1-mediated resistance. In addition, deep-sequencing analyses revealed that resistance breakdown does not always coincide with fixed mutations. Actually, virulence mutations that are present in a small proportion of the virus population can be sufficient for resistance breakdown. Considering the spatial distribution of RYMV strains in Africa and their ability to overcome the RYMV resistance genes and alleles, we established a resistance-breaking risk map to optimize strategies for the deployment of sustainable and resistant rice lines in Africa.

Alternative mutational pathways, outside the VPg, of rice yellow mottle virus to overcome eIF(iso)4G-mediated rice resistance under strong genetic constraints.

The adaptation of rice yellow mottle virus (RYMV) to rymv1-mediated resistance has been reported to involve mutations in the viral genome-linked protein (VPg). In this study, we analysed several cases of rymv1-2 resistance breakdown by an isolate with low adaptability. Surprisingly, in these rarely occurring resistance-breaking (RB) genotypes, mutations were detected outside the VPg, in the ORF2a/ORF2b overlapping region. The causal role of three mutations associated with rymv1-2 resistance breakdown was validated via directed mutagenesis of an infectious clone. In resistant plants, these mutations increased viral accumulation as efficiently as suboptimal RB mutations in the VPg. Interestingly, these mutations are located in a highly conserved, but unfolded, domain. Altogether, our results indicate that under strong genetic constraints, a priori unfit genotypes can follow alternative mutational pathways, i.e. outside the VPg, to overcome rymv1-2 resistance.

Sites under positive selection modulate the RNA silencing suppressor activity of rice yellow mottle virus movement protein P1.

RNA silencing is a eukaryotic mechanism for RNA-based gene regulation that plays an essential role in diverse biological processes, such as defence against viral infections. The P1 of rice yellow mottle virus (RYMV) is a movement protein and displays RNA silencing suppression activity with variable efficiency, depending on the origin of the isolates. In this study, the positive selection pressure acting on the P1 protein gene was assessed. A site-by-site analysis of the dN/dS ratio was performed and 18 positively selected sites were identified. Four of these were mutated, and the ability to suppress RNA silencing was evaluated for the resulting mutants in a transient expression assay. All mutations affected quantitatively RNA silencing suppression, one caused a significant decrease in the activity and three significantly increased it. This work demonstrates, for what is to the best of our knowledge the first time, that the RYMV gene encoding the P1 RNA silencing suppressor is under adaptive evolution.

Historical contingencies modulate the adaptability of Rice yellow mottle virus.

The rymv1-2 and rymv1-3 alleles of the RYMV1 resistance to Rice yellow mottle virus (RYMV), coded by an eIF(iso)4G1 gene, occur in a few cultivars of the Asiatic (Oryza sativa) and African (O. glaberrima) rice species, respectively. The most salient feature of the resistance breaking (RB) process is the converse genetic barrier to rymv1-2 and rymv1-3 resistance breakdown. This specificity is modulated by the amino acid (glutamic acid vs. threonine) at codon 49 of the Viral Protein genome-linked (VPg), a position which is adjacent to the virulence codons 48 and 52. Isolates with a glutamic acid (E) do not overcome rymv1-3 whereas those with a threonine (T) rarely overcome rymv1-2. We found that isolates with T49 had a strong selective advantage over isolates with E49 in O. glaberrima susceptible cultivars. This explains the fixation of the mutation T49 during RYMV evolution and accounts for the diversifying selection estimated at codon 49. Better adapted to O. glaberrima, isolates with T49 are also more prone than isolates with E49 to fix rymv1-3 RB mutations at codon 52 in resistant O. glaberrima cultivars. However, subsequent genetic constraints impaired the ability of isolates with T49 to fix rymv1-2 RB mutations at codons 48 and 52 in resistant O. sativa cultivars. The origin and role of the amino acid at codon 49 of the VPg exemplifies the importance of historical contingencies in the ability of RYMV to overcome RYMV1 resistance.

Cassava Ivorian bacilliform virus is a member of the genus Anulavirus.

The complete genomic sequence of Cassava Ivorian bacilliform virus (CIBV) is described. The virus has a genomic organization similar to that of pelargonium zonate spot virus (PZSV), the type member of the genus Anulavirus, but it is most closely related to a second, recently described, anulavirus, Amazon lily mild mottle virus (ALiMMV).

Modeling the velocity of evolving lineages and predicting dispersal patterns.

Accurate estimation of the dispersal velocity or speed of evolving organisms is no mean feat. In fact, existing probabilistic models in phylogeography or spatial population genetics generally do not provide an adequate framework to define velocity in a relevant manner. For instance, the very concept of instantaneous speed simply does not exist under one of the most popular approaches that models the evolution of spatial coordinates as Brownian trajectories running along a phylogeny [30]. Here, we introduce a new family of models - the so-called "Phylogenetic Integrated Velocity" (PIV) models - that use Gaussian processes to explicitly model the velocity of evolving lineages instead of focusing on the fluctuation of spatial coordinates over time. We describe the properties of these models and show an increased accuracy of velocity estimates compared to previous approaches. Analyses of West Nile virus data in the U.S.A. indicate that PIV models provide sensible predictions of the dispersal of evolving pathogens at a one-year time horizon. These results demonstrate the feasibility and relevance of predictive phylogeography in monitoring epidemics in time and space.

How fast are viruses spreading in the wild?

Genomic data collected from viral outbreaks can be exploited to reconstruct the dispersal history of viral lineages in a two-dimensional space using continuous phylogeographic inference. These spatially explicit reconstructions can subsequently be used to estimate dispersal metrics allowing to unveil the dispersal dynamics and evaluate the capacity to spread among hosts. Heterogeneous sampling intensity of genomic sequences can however impact the accuracy of dispersal insights gained through phylogeographic inference. In our study, we implement a simulation framework to evaluate the robustness of three dispersal metrics - a lineage dispersal velocity, a diffusion coefficient, and an isolation-by-distance signal metric - to the sampling effort. Our results reveal that both the diffusion coefficient and isolation-by-distance signal metrics appear to be robust to the number of samples considered for the phylogeographic reconstruction. We then use these two dispersal metrics to compare the dispersal pattern and capacity of various viruses spreading in animal populations. Our comparative analysis reveals a broad range of isolation-by-distance patterns and diffusion coefficients mostly reflecting the dispersal capacity of the main infected host species but also, in some cases, the likely signature of rapid and/or long-distance dispersal events driven by human-mediated movements through animal trade. Overall, our study provides key recommendations for the lineage dispersal metrics to consider in future studies and illustrates their application to compare the spread of viruses in various settings.

Genetic diversity and epidemic histories of rice yellow mottle virus in Ghana.

Rice yellow mottle virus (RYMV) has persisted as a major biotic constraint to rice production in Africa. However, no data on RYMV epidemics were available in Ghana, although it is an intensive rice-producing country. Surveys were performed from 2010 to 2020 in eleven rice-growing regions of Ghana. Symptom observations and serological detections confirmed that RYMV is circulating in most of these regions. Coat protein gene and complete genome sequencings revealed that RYMV in Ghana almost exclusively belongs to the strain S2, one of the strains covering the largest area in West Africa. We also detected the presence of the S1ca strain which is being reported for the first time outside its area of origin. These results suggested a complex epidemiological history of RYMV in Ghana and a recent expansion of S1ca to West Africa. Phylogeographic analyses reconstructed at least five independent RYMV introductions in Ghana for the last 40 years, probably due to rice cultivation intensification in West Africa leading to a better circulation of RYMV. In addition to identifying some routes of RYMV dispersion in Ghana, this study contributes to the epidemiological surveillance of RYMV and helps to design disease management strategies, especially through breeding for rice disease resistance.

Protein-RNA linkage and post-translational modifications of two sobemovirus VPgs.

Sobemoviruses possess a viral genome-linked protein (VPg) attached to the 5' end of viral RNA. VPg is processed from the viral polyprotein. In the current study, Cocksfoot mottle virus (CfMV) and Rice yellow mottle virus (RYMV) VPgs were purified from virions and analysed by mass spectrometry. The cleavage sites in the polyprotein and thereof the termini of VPg were experimentally proven. The lengths of the mature VPgs were determined to be 78 and 79 aa residues, respectively. The amino acid residues covalently linked to RNA in the two VPgs were, surprisingly, not conserved; it is a tyrosine at position 5 of CfMV VPg and serine at position 1 of RYMV VPg. Phosphorylations were identified in CfMV and RYMV VPgs with two positionally similar locations T20/S14 and S71/S72, respectively. RYMV VPg contains an additional phosphorylation site at S41.

Rivers and landscape ecology of a plant virus, Rice yellow mottle virus along the Niger Valley.

To investigate the spread of Rice yellow mottle virus (RYMV) along the Niger River, regular sampling of virus isolates was conducted along 500 km of the Niger Valley in the Republic of Niger and was complemented by additional sampling in neighbouring countries in West Africa and Central Africa. The spread of RYMV into and within the Republic of Niger was inferred as a continuous process using a Bayesian statistical framework applied previously to reconstruct its dispersal history in West Africa, East Africa, and Madagascar. The spatial resolution along this section of the Niger River was the highest implemented for RYMV and possibly for any plant virus. We benefited from the results of early field surveys of the disease for the validation of the phylogeographic reconstruction and from the well-documented history of rice cultivation changes along the Niger River for their interpretation. As a prerequisite, the temporal signal of the RYMV data sets was revisited in the light of recent methodological advances. The role of the hydrographic network of the Niger Basin in RYMV spread was examined, and the link between virus population dynamics and the extent of irrigated rice was assessed. RYMV was introduced along the Niger River in the Republic of Niger in the early 1980s from areas to the southwest of the country where rice was increasingly grown. Viral spread was triggered by a major irrigation scheme made of a set of rice perimeters along the river valley. The subsequent spatial and temporal host continuity and the inoculum build-up allowed for a rapid spread of RYMV along the Niger River, upstream and downstream, over hundreds of kilometres, and led to the development of severe epidemics. There was no evidence of long-distance dissemination of the virus through natural water. Floating rice in the main meanders of the Middle Niger did not contribute to virus dispersal from West Africa to Central Africa. RYMV along the Niger River is an insightful example of how agricultural intensification favours pathogen emergence and spread.

Direct interaction between the Rice yellow mottle virus (RYMV) VPg and the central domain of the rice eIF(iso)4G1 factor correlates with rice susceptibility and RYMV virulence.

The adaptation of Rice yellow mottle virus (RYMV) to recessive resistance mediated by the rymv1-2 allele has been reported as a model to study the emergence and evolution of virulent variants. The resistance and virulence factors have been identified as eukaryotic translation initiation factor eIF(iso)4G1 and viral genome-linked protein (VPg), respectively, but the molecular mechanisms involved in their interaction are still unknown. In this study, we demonstrated a direct interaction between RYMV VPg and the central domain of rice eIF(iso)4G1 both in vitro, using recombinant proteins, and in vivo, using a yeast two-hybrid assay. Insertion of the E309K mutation in eIF(iso)4G1, conferring resistance in planta, strongly diminished the interaction with avirulent VPg. The efficiency of the major virulence mutations at restoring the interaction with the resistance protein was assessed. Our results explain the prevalence of virulence mutations fixed during experimental evolution studies and are consistent with the respective viral RNA accumulation levels of avirulent and virulent isolates. Our results also explain the origin of the residual multiplication of wild-type isolates in rymv1-2-resistant plants and the role of genetic context in the poor adaptability of the S2/S3 strain. Finally, the strategies of RYMV and members of family Potyviridae to overcome recessive resistance were compared.

Diversification of rice yellow mottle virus and related viruses spans the history of agriculture from the neolithic to the present.

The mechanisms of evolution of plant viruses are being unraveled, yet the timescale of their evolution remains an enigma. To address this critical issue, the divergence time of plant viruses at the intra- and inter-specific levels was assessed. The time of the most recent common ancestor (TMRCA) of Rice yellow mottle virus (RYMV; genus Sobemovirus) was calculated by a Bayesian coalescent analysis of the coat protein sequences of 253 isolates collected between 1966 and 2006 from all over Africa. It is inferred that RYMV diversified approximately 200 years ago in Africa, i.e., centuries after rice was domesticated or introduced, and decades before epidemics were reported. The divergence time of sobemoviruses and viruses of related genera was subsequently assessed using the age of RYMV under a relaxed molecular clock for calibration. The divergence time between sobemoviruses and related viruses was estimated to be approximately 9,000 years, that between sobemoviruses and poleroviruses approximately 5,000 years, and that among sobemoviruses approximately 3,000 years. The TMRCA of closely related pairs of sobemoviruses, poleroviruses, and luteoviruses was approximately 500 years, which is a measure of the time associated with plant virus speciation. It is concluded that the diversification of RYMV and related viruses has spanned the history of agriculture, from the Neolithic age to the present.

Theme and variations in the evolutionary pathways to virulence of an RNA plant virus species.

The diversity of a highly variable RNA plant virus was considered to determine the range of virulence substitutions, the evolutionary pathways to virulence, and whether intraspecific diversity modulates virulence pathways and propensity. In all, 114 isolates representative of the genetic and geographic diversity of Rice yellow mottle virus (RYMV) in Africa were inoculated to several cultivars with eIF(iso)4G-mediated Rymv1-2 resistance. Altogether, 41 virulent variants generated from ten wild isolates were analyzed. Nonconservative amino acid replacements at five positions located within a stretch of 15 codons in the central region of the 79-aa-long protein VPg were associated with virulence. Virulence substitutions were fixed predominantly at codon 48 in most strains, whatever the host genetic background or the experimental conditions. There were one major and two isolate-specific mutational pathways conferring virulence at codon 48. In the prevalent mutational pathway I, arginine (AGA) was successively displaced by glycine (GGA) and glutamic acid (GAA). Substitutions in the other virulence codons were displaced when E48 was fixed. In the isolate-specific mutational pathway II, isoleucine (ATA) emerged and often later coexisted with valine (GTA). In mutational pathway III, arginine, with the specific S2/S3 strain codon usage AGG, was displaced by tryptophane (TGG). Mutational pathway I never arose in the widely spread West African S2/S3 strain because G48 was not infectious in the S2/S3 genetic context. Strain S2/S3 least frequently overcame resistance, whereas two geographically localized variants of the strain S4 had a high propensity to virulence. Codons 49 and 26 of the VPg, under diversifying selection, are candidate positions in modulating the genetic barriers to virulence. The theme and variations in the evolutionary pathways to virulence of RYMV illustrates the extent of parallel evolution within a highly variable RNA plant virus species.

Intrinsic disorder in Viral Proteins Genome-Linked: experimental and predictive analyses.

BACKGROUND: VPgs are viral proteins linked to the 5' end of some viral genomes. Interactions between several VPgs and eukaryotic translation initiation factors eIF4Es are critical for plant infection. However, VPgs are not restricted to phytoviruses, being also involved in genome replication and protein translation of several animal viruses. To date, structural data are still limited to small picornaviral VPgs. Recently three phytoviral VPgs were shown to be natively unfolded proteins. RESULTS: In this paper, we report the bacterial expression, purification and biochemical characterization of two phytoviral VPgs, namely the VPgs of Rice yellow mottle virus (RYMV, genus Sobemovirus) and Lettuce mosaic virus (LMV, genus Potyvirus). Using far-UV circular dichroism and size exclusion chromatography, we show that RYMV and LMV VPgs are predominantly or partly unstructured in solution, respectively. Using several disorder predictors, we show that both proteins are predicted to possess disordered regions. We next extend theses results to 14 VPgs representative of the viral diversity. Disordered regions were predicted in all VPg sequences whatever the genus and the family. CONCLUSION: Based on these results, we propose that intrinsic disorder is a common feature of VPgs. The functional role of intrinsic disorder is discussed in light of the biological roles of VPgs.

Comparing patterns and scales of plant virus phylogeography: Rice yellow mottle virus in Madagascar and in continental Africa.

Rice yellow mottle virus (RYMV) in Madagascar Island provides an opportunity to study the spread of a plant virus disease after a relatively recent introduction in a large and isolated country with a heterogeneous host landscape ecology. Here, we take advantage of field survey data on the occurrence of RYMV disease throughout Madagascar dating back to the 1970s, and of virus genetic data from ninety-four isolates collected since 1989 in most regions of the country to reconstruct the epidemic history. We find that the Malagasy isolates belong to a unique recombinant strain that most likely entered Madagascar through a long-distance introduction from the most eastern part of mainland Africa. We infer the spread of RYMV as a continuous process using a Bayesian statistical framework. In order to calibrate the time scale in calendar time units in this analysis, we pool the information about the RYMV evolutionary rate from several geographical partitions. Whereas the field surveys and the phylogeographic reconstructions both point to a rapid southward invasion across hundreds of kilometers throughout Madagascar within three to four decades, they differ on the inferred origin location and time of the epidemic. The phylogeographic reconstructions suggest a lineage displacement and unveil a re-invasion of the northern regions that may have remained unnoticed otherwise. Despite ecological differences that could affect the transmission potential of RYMV in Madagascar and in mainland Africa, we estimate similar invasion and dispersal rates. We could not identify environmental factors that have a relevant impact on the lineage dispersal velocity of RYMV in Madagascar. This study highlights the value and complementarity of (historical) nongenetic and (more contemporaneous) genetic surveillance data for reconstructing the history of spread of plant viruses.

Complete Genome Sequences of Rice Yellow Mottle Virus Isolates from the Federal Democratic Republic of Ethiopia.

The full-length genomes of two isolates of Rice yellow mottle virus from Ethiopia were sequenced. A comparison with 28 sequences from East Africa showed that they clustered within a new strain named S4et, related to the S4mg and S4ug strains found in the Lake Victoria Basin and Madagascar, respectively.

Genetic diversity and silencing suppression effects of Rice yellow mottle virus and the P1 protein.

BACKGROUND: PTGS (post-transcriptional gene silencing) is used to counter pathogenic invasions, particularly viruses. In return, many plant viruses produce proteins which suppress silencing directed against their RNA. The diversity of silencing suppression at the species level in natural hosts is unknown. RESULTS: We investigated the functional diversity of silencing suppression among isolates of the African RYMV (Rice yellow mottle virus) in rice. The RYMV-P1 protein is responsible for cell-to-cell movement and is a silencing suppressor. Transgenic gus-silencing rice lines were used to investigate intra-specific and serogroup silencing suppression diversity at two different levels: that of the virion and the P1 silencing suppressor protein. Our data provide evidence that silencing suppression is a universal phenomenon for RYMV species. However, we found considerable diversity in their ability to suppress silencing which was not linked to RYMV phylogeny, or pathogenicity. At the level of the silencing suppressor P1 alone, we found similar results to those previously found at the virion level. In addition, we showed that cell-to-cell movement of P1 was crucial for the efficiency of silencing suppression. Mutagenesis of P1 demonstrated a strong link between some amino acids and silencing suppression features with, one on the hand, the conserved amino acids C95 and C64 involved in cell-to-cell movement and the strength of suppression, respectively, and on the other hand, the non conserved F88 was involved in the strength of silencing suppression. CONCLUSION: We demonstrated that intra-species diversity of silencing suppression is highly variable and by mutagenesis of P1 we established the first link between silencing suppression and genetic diversity. These results are potentially important for understanding virus-host interactions.