Modeling multipartite virus evolution: the genome formula facilitates rapid adaptation to heterogeneous environments†

Replication-independent change in the frequencies of distinct genome segments of a multipartite virus during its transit within aphid vectors

Multipartite viruses exhibit a fragmented genome composed of several nucleic acid segments individually packaged in distinct viral particles. The genome of all species of the genus Nanovirus holds eight segments, which accumulate at a very specific and reproducible relative frequency in the host plant tissues. In a given host species, the steady state pattern of the segments' relative frequencies is designated the genome formula and is thought to have an adaptive function through the modulation of gene expression. Nanoviruses are aphid-transmitted circulative non-propagative viruses, meaning that the virus particles are internalized into the midgut cells, transferred to the hemolymph, and then to the saliva, with no replication during this transit. Unexpectedly, a previous study on the faba bean necrotic stunt virus revealed that the genome formula changes after ingestion by aphids. We investigate here the possible mechanism inducing this change by first comparing the relative segment frequencies in different compartments of the aphid. We show that changes occur both in the midgut lumen and in the secreted saliva but not in the gut, salivary gland, or hemolymph. We further establish that the viral particles differentially resist physicochemical variations, in particular pH, ionic strength, and/or type of salt, depending on the encapsidated segment. We thus propose that the replication-independent genome formula changes within aphids are not adaptive, contrary to changes occurring in plants, and most likely reflect a fortuitous differential degradation of virus particles containing distinct segments when passing into extra-cellular media such as gastric fluid or saliva.

IMPORTANCE

The genome of multipartite viruses is composed of several segments individually packaged into distinct viral particles. Each segment accumulates at a specific frequency that depends on the host plant species and regulates gene expression. Intriguingly, the relative frequencies of the genome segments also change when the octopartite faba bean necrotic stunt virus (FBNSV) is ingested by aphid vectors, despite the present view that this virus travels through the aphid gut and salivary glands without replicating. By monitoring the genomic composition of FBNSV populations during the transit in aphids, we demonstrate here that the changes take place extracellularly in the gut lumen and in the saliva. We further show that physicochemical factors induce differential degradation of viral particles depending on the encapsidated segment. We propose that the replication-independent changes within the insect vector are not adaptive and result from the differential stability of virus particles containing distinct segments according to environmental parameters.

Robust Approaches to the Quantitative Analysis of Genome Formula Variation in Multipartite and Segmented Viruses

When viruses have segmented genomes, the set of frequencies describing the abundance of segments is called the genome formula. The genome formula is often unbalanced and highly variable for both segmented and multipartite viruses. A growing number of studies are quantifying the genome formula to measure its effects on infection and to consider its ecological and evolutionary implications. Different approaches have been reported for analyzing genome formula data, including qualitative description, applying standard statistical tests such as ANOVA, and customized analyses. However, these approaches have different shortcomings, and test assumptions are often unmet, potentially leading to erroneous conclusions. Here, we address these challenges, leading to a threefold contribution. First, we propose a simple metric for analyzing genome formula variation: the genome formula distance. We describe the properties of this metric and provide a framework for understanding metric values. Second, we explain how this metric can be applied for different purposes, including testing for genome-formula differences and comparing observations to a reference genome formula value. Third, we re-analyze published data to illustrate the applications and weigh the evidence for previous conclusions. Our re-analysis of published datasets confirms many previous results but also provides evidence that the genome formula can be carried over from the inoculum to the virus population in a host. The simple procedures we propose contribute to the robust and accessible analysis of genome-formula data.

Conserved untranslated regions of multipartite viruses: Natural markers of novel viral genomic components and tags of viral evolution

Viruses with split genomes are classified as being either segmented or multipartite based on whether their genomic segments occur within a single virion or across different virions. Despite variations in number and sequence during evolution, the genomic segments of many viruses are conserved within the untranslated regions (UTRs). In this study, we present a methodology that combines RNA sequencing with iterative BLASTn of UTRs (https://github.com/qq371260/Iterative-blast-v.1.0) to identify new viral genomic segments. Some novel multipartite-like viruses related to the phylum Kitrinoviricota were annotated using sequencing data from field plant samples and public databases. We identified potentially plant-infecting jingmen-related viruses (order Amarillovirales) and jivi-related viruses (order Martellivirales) with at least six genomic components. The number of RNA molecules associated with a genome of the novel viruses in the families Closteroviridae, Kitaviridae, and Virgaviridae within the order Martellivirales reached five. Several of these viruses seem to represent new taxa at the subgenus, genus, and family levels. The diversity of novel genomic components and the multiple duplication of proteins or protein domains within single or multiple genomic components emphasize the evolutionary roles of genetic recombination (horizontal gene transfer), reassortment, and deletion. The relatively conserved UTRs at the genome level might explain the relationships between monopartite and multipartite viruses, as well as how subviral agents such as defective RNAs and satellite viruses can coexist with their helper viruses.

New clades of viruses infecting the obligatory biotroph Bremia lactucae representing distinct evolutionary trajectory for viruses infecting oomycetes

Recent advances in high throughput sequencing (HTS) approaches allowed a broad exploration of viromes from different fungal hosts, unveiling a great diversity of mycoviruses with interesting evolutionary features. The word mycovirus historically applies also to viruses infecting oomycetes but most studies are on viruses infecting fungi, with less mycoviruses found and characterized in oomycetes, particularly in the obligatory biotrophs. We, here, describe the first virome associated to Bremia lactucae, the causal agent of lettuce downy mildew, which is an important biotrophic pathogen for lettuce production and a model system for the molecular aspects of the plant-oomycetes interactions. Among the identified viruses, we could detect (1) two new negative sense ssRNA viruses related to the yueviruses, (2) the first example of permuted RdRp in a virus infecting fungi/oomycetes, (3) a new group of bipartite dsRNA viruses showing evidence of recent bi-segmentation and concomitantly, a possible duplication event bringing a bipartite genome to tripartite, (4) a first representative of a clade of viruses with evidence of recombination between distantly related viruses, (5) a new open reading frame (ORF)an virus encoding for an RdRp with low homology to known RNA viruses, and (6) a new virus, belonging to riboviria but not conserved enough to provide a conclusive phylogenetic placement that shows evidence of a recombination event between a kitrinoviricota-like and a pisuviricota-like sequence. The results obtained show a great diversity of viruses and evolutionary mechanisms previously unreported for oomycetes-infecting viruses, supporting the existence of a large diversity of oomycetes-specific viral clades ancestral of many fungal and insect virus clades.

The genome formula of a multipartite virus is regulated both at the individual segment and the segment group levels

Differential accumulation of the distinct genome segments is a common feature of viruses with segmented genomes. The reproducible and specific pattern of genome segment accumulation within the host is referred to as the "genome formula". There is speculation and some experimental support for a functional role of the genome formula by modulating gene expression through copy number variations. However, the mechanisms of genome formula regulation have not yet been identified. In this study, we investigated whether the genome formula of the octopartite nanovirus faba bean necrotic stunt virus (FBNSV) is regulated by processes acting at the individual segment vs. viral population levels. We used a leaf infiltration system to show that the two most accumulated genome segments of the FBNSV possess a greater intrinsic accumulation capacity in Vicia faba tissues than the other segments. Nevertheless, processes acting at the individual segment level are insufficient to generate the genome formula, suggesting the involvement of additional mechanisms acting at the supra-segment level. Indeed, the absence of segments with important functions during systemic infection strongly modifies the relative frequency of the others, indicating that the genome formula is a property of the segment group. Together, these results demonstrate that the FBNSV genome formula is shaped by a complex process acting at both the individual segment and the segment group levels.

Open questions in the social lives of viruses

Social interactions among viruses occur whenever multiple viral genomes infect the same cells, hosts, or populations of hosts. Viral social interactions range from cooperation to conflict, occur throughout the viral world, and affect every stage of the viral lifecycle. The ubiquity of these social interactions means that they can determine the population dynamics, evolutionary trajectory, and clinical progression of viral infections. At the same time, social interactions in viruses raise new questions for evolutionary theory, providing opportunities to test and extend existing frameworks within social evolution. Many opportunities exist at this interface: Insights into the evolution of viral social interactions have immediate implications for our understanding of the fundamental biology and clinical manifestation of viral diseases. However, these opportunities are currently limited because evolutionary biologists only rarely study social evolution in viruses. Here, we bridge this gap by (1) summarizing the ways in which viruses can interact socially, including consequences for social evolution and evolvability; (2) outlining some open questions raised by viruses that could challenge concepts within social evolution theory; and (3) providing some illustrative examples, data sources, and conceptual questions, for studying the natural history of social viruses.

cmv1-Mediated Resistance to CMV in Melon Can Be Overcome by Mixed Infections with Potyviruses

Resistance to cucumber mosaic virus (CMV) strain LS in melon is controlled by the gene cmv1, which restricts phloem entry. In nature, CMV is commonly found in mixed infections, particularly with potyviruses, where a synergistic effect is frequently produced. We have explored the possibility that this synergism could help CMV-LS to overcome cmv1-mediated resistance. We demonstrate that during mixed infection with a potyvirus, CMV-LS is able to overcome cmv1-controlled resistance and develop a systemic infection and that this ability does not depend on an increased accumulation of CMV-LS in mechanically inoculated cotyledons. Likewise, during a mixed infection initiated by aphids, the natural vector of both cucumoviruses and potyviruses that can very efficiently inoculate plants with a low number of virions, CMV-LS also overcomes cmv1-controlled resistance. This indicates that in the presence of a potyvirus, even a very low amount of inoculum, can be sufficient to surpass the resistance and initiate the infection. These results indicate that there is an important risk for this resistance to be broken in nature as a consequence of mixed infections, and therefore, its deployment in elite cultivars would not be enough to ensure a long-lasting resistance.

Constant ratio between the genomic components of bipartite begomoviruses during infection and transmission

The genomic components of multipartite viruses are encapsidated in separate virus particles, and the frequencies of genomic components represent one of the key genetic features. Many begomoviruses of economic significance are bipartite, and the details of the association between their genomic components remain largely unexplored. We first analyzed the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of the squash leaf curl China virus (SLCCNV) in plants and found that while the quantities of DNA-A and DNA-B varied significantly during infection, the B/A ratio remained constant. We then found that changes in the B/A ratio in agrobacteria inoculum may significantly alter the B/A ratio in plants at 6 days post inoculation, but the differences disappeared shortly thereafter. We next showed that while the quantities of DNA-A and DNA-B among plants infected by agrobacteria, sap transmission and whitefly-mediated transmission differed significantly, the B/A ratios were similar. Further analysis of gene expression revealed that the ratio of the expression of genes encoded by DNA-A and DNA-B varied significantly during infection. Finally, we monitored the temporal dynamics of the quantities of DNA-A and DNA-B and the B/A ratio of another bipartite begomovirus, and a constant B/A ratio was similarly observed. Our findings highlight the maintenance of a constant ratio between the two genomic components of bipartite begomoviruses during infection and transmission, and provide new insights into the biology of begomoviruses.

A collection of Trichoderma isolates from natural environments in Sardinia reveals a complex virome that includes negative-sense fungal viruses with unprecedented genome organizations

Trichoderma genus includes soil-inhabiting fungi that provide important ecosystem services in their interaction with plants and other fungi, as well as biocontrol of fungal plant diseases. A collection of Trichoderma isolates from Sardinia has been previously characterized, but here we selected 113 isolates, representatives of the collection, and characterized their viral components. We carried out high-throughput sequencing of ribosome-depleted total RNA following a bioinformatics pipeline that detects virus-derived RNA-directed RNA polymerases (RdRps) and other conserved viral protein sequences. This pipeline detected seventeen viral RdRps with two of them corresponding to viruses already detected in other regions of the world and the remaining fifteen representing isolates of new putative virus species. Surprisingly, eight of them are from new negative-sense RNA viruses, a first in the genus Trichoderma. Among them is a cogu-like virus, closely related to plant-infecting viruses. Regarding the positive-sense viruses, we report the presence of an 'ormycovirus' belonging to a recently characterized group of bisegmented single-stranded RNA viruses with uncertain phylogenetic assignment. Finally, for the first time, we report a bisegmented member of Mononegavirales which infects fungi. The proteins encoded by the second genomic RNA of this virus were used to re-evaluate several viruses in the Penicillimonavirus and Plasmopamonavirus genera, here shown to be bisegmented and encoding a conserved polypeptide that has structural conservation with the nucleocapsid domain of rhabdoviruses.

Evaluation of sequencing and PCR-based methods for the quantification of the viral genome formula

Viruses show great diversity in their genome organization. Multipartite viruses package their genome segments into separate particles, most or all of which are required to initiate infection in the host cell. The benefits of such seemingly inefficient genome organization are not well understood. One hypothesised benefit of multipartition is that it allows for flexible changes in gene expression by altering the frequency of each genome segment in different environments, such as encountering different host species. The ratio of the frequency of segments is termed the genome formula (GF). Thus far, formal studies quantifying the GF have been performed for well-characterised virus-host systems in experimental settings using RT-qPCR. However, to understand GF variation in natural populations or novel virus-host systems, a comparison of several methods for GF estimation including high-throughput sequencing (HTS) based methods is needed. Currently, it is unclear how HTS-methods compare a golden standard, such as RT-qPCR. Here we show a comparison of multiple GF quantification methods (RT-qPCR, RT-digital PCR, Illumina RNAseq and Nanopore direct RNA sequencing) using three host plants (Nicotiana tabacum, Nicotiana benthamiana, and Chenopodium quinoa) infected with cucumber mosaic virus (CMV), a tripartite RNA virus. Our results show that all methods give roughly similar results, though there is a significant method effect on genome formula estimates. While the RT-qPCR and RT-dPCR GF estimates are congruent, the GF estimates from HTS methods deviate from those found with PCR. Our findings emphasize the need to tailor the GF quantification method to the experimental aim, and highlight that it may not be possible to compare HTS and PCR-based methods directly. The difference in results between PCR-based methods and HTS highlights that the choice of quantification technique is not trivial.

Gene Copy Number Variations at the Within-Host Population Level Modulate Gene Expression in a Multipartite Virus

Multipartite viruses have a segmented genome, with each segment encapsidated separately. In all multipartite virus species for which the question has been addressed, the distinct segments reproducibly accumulate at a specific and host-dependent relative frequency, defined as the ‘genome formula’. Here, we test the hypothesis that the multipartite genome organization facilitates the regulation of gene expression via changes of the genome formula and thus via gene copy number variations. In a first experiment, the faba bean necrotic stunt virus (FBNSV), whose genome is composed of eight DNA segments each encoding a single gene, was inoculated into faba bean or alfalfa host plants, and the relative concentrations of the DNA segments and their corresponding messenger RNAs (mRNAs) were monitored. In each of the two host species, our analysis consistently showed that the genome formula variations modulate gene expression, the concentration of each genome segment linearly and positively correlating to that of its cognate mRNA but not of the others. In a second experiment, twenty parallel FBNSV lines were transferred from faba bean to alfalfa plants. Upon host switching, the transcription rate of some genome segments changes, but the genome formula is modified in a way that compensates for these changes and maintains a similar ratio between the various viral mRNAs. Interestingly, a deep-sequencing analysis of these twenty FBNSV lineages demonstrated that the host-related genome formula shift operates independently of DNA-segment sequence mutation. Together, our results indicate that nanoviruses are plastic genetic systems, able to transiently adjust gene expression at the population level in changing environments, by modulating the copy number but not the sequence of each of their genes.

Effects of an alphasatellite on life cycle of the nanovirus Faba bean necrotic yellows virus

Nanoviruses are plant viruses with a multipartite single-stranded DNA (ssDNA) genome. Alphasatellites are commonly associated with nanovirus infections, but their putative impact on their helper viruses is unknown. In this study, we investigated the role of subterranean clover stunt alphasatellite 1 (hereafter named SCSA 1) on various important traits of faba bean necrotic yellows virus (FBNYV) in its host plant Vicia faba and aphid vector Acyrthosiphon pisum, including disease symptoms, viral accumulation and transmission. The results indicate that SCSA 1 does not affect the symptom severity nor the overall FBNYV accumulation in V. faba, but changes the relative amounts of its different genomic segments. Moreover, the association of SCSA 1 with FBNYV increases the rate of plant-to-plant transmission by a process seemingly unrelated to simple increase of the viral accumulation in the vector. These results represent the first study on the impact of an alphasatellite on the biology of its helper nanovirus. They suggest that SCSA 1 may benefit FBNYV, but the genericity of this conclusion is discussed and questioned. Importance Alphasatellites are circular single stranded DNA molecules frequently found in association with natural isolates of nanoviruses and some geminiviruse, the two ssDNA plant infecting virus families. While the implications of alphasatellite presence in geminivirus infections are relatively well documented, comparable studies on alphasatellites associated with nanoviruses are not available. Here we confirm that subterranean clover stunt alphasatellite 1 affects different traits of its helper nanovirus, faba bean necrotic yellows virus, both in the host plant and aphid vector. We show that the frequencies of the virus segments change in the presence of alphasatellite, in both plant and vector. We also confirm that while within-plant virus load and symptom are not affected by alphasatellite, the presence of alphasatellite decreases within-aphid virus load, but significantly increases virus transmission rate, so may confer a possible evolutionary advantage for the helper virus.

Nine viruses from eight lineages exhibiting new evolutionary modes that co-infect a hypovirulent phytopathogenic fungus

Mycoviruses are an important component of the virosphere, but our current knowledge of their genome organization diversity and evolution remains rudimentary. In this study, the mycovirus composition in a hypovirulent strain of Sclerotinia sclerotiorum was molecularly characterized. Nine mycoviruses were identified and assigned into eight potential families. Of them, six were close relatives of known mycoviruses, while the other three had unique genome organizations and evolutionary positions. A deltaflexivirus with a tripartite genome has evolved via arrangement and horizontal gene transfer events, which could be an evolutionary connection from unsegmented to segmented RNA viruses. Two mycoviruses had acquired a second helicase gene by two different evolutionary mechanisms. A rhabdovirus representing an independent viral evolutionary branch was the first to be confirmed to occur naturally in fungi. The major hypovirulence-associated factor, an endornavirus, was finally corroborated. Our study expands the diversity of mycoviruses and potential virocontrol agents, and also provides new insights into virus evolutionary modes including virus genome segmentation.

Identification of mycoviruses in phytopathogenic fungi is necessary for understanding the origin of viruses and developing virocontrol strategies to protect plants. Nine mycoviruses with RNA genomes were identified in a hypovirulent strain of Sclerotinia sclerotiorum and were classified into eight potential viral families, suggesting that the composition of mycoviral communities was complex in this single fungal strain. They included four previously characterized mycoviruses and three distant relatives of known mycoviruses, as well as the first reports of a deltaflexivirus with a tripartite genome, and a fungal rhabdovirus. In addition, we found an endornavirus associated with hypovirulence in a phytopathogenic fungus. Our study makes a significant contribution because it not only expands the diversity-related knowledge of mycoviruses and potential virocontrol agents, but also provides new insights into mycovirus evolution.

Ever-increasing viral diversity associated with the red imported fire ant Solenopsis invicta (Formicidae: Hymenoptera)

Background

Advances in sequencing and analysis tools have facilitated discovery of many new viruses from invertebrates, including ants. Solenopsis invicta is an invasive ant that has quickly spread worldwide causing significant ecological and economic impacts. Its virome has begun to be characterized pertaining to potential use of viruses as natural enemies. Although the S. invicta virome is the best characterized among ants, most studies have been performed in its native range, with less information from invaded areas.

Methods

Using a metatranscriptome approach, we further identified and molecularly characterized virus sequences associated with S. invicta, in two introduced areas, U.S and Taiwan. The data set used here was obtained from different stages (larvae, pupa, and adults) of S. invicta life cycle. Publicly available RNA sequences from GenBank’s Sequence Read Archive were downloaded and de novo assembled using CLC Genomics Workbench 20.0.1. Contigs were compared against the non-redundant protein sequences and those showing similarity to viral sequences were further analyzed.

Results

We characterized five putative new viruses associated with S. invicta transcriptomes. Sequence comparisons revealed extensive divergence across ORFs and genomic regions with most of them sharing less than 40% amino acid identity with those closest homologous sequences previously characterized. The first negative-sense single-stranded RNA virus genomic sequences included in the orders Bunyavirales and Mononegavirales are reported. In addition, two positive single-strand virus genome sequences and one single strand DNA virus genome sequence were also identified. While the presence of a putative tenuivirus associated with S. invicta was previously suggested to be a contamination, here we characterized and present strong evidence that Solenopsis invicta virus 14 (SINV-14) is a tenui-like virus that has a long-term association with the ant. Furthermore, based on virus sequence abundance compared to housekeeping genes, phylogenetic relationships, and completeness of viral coding sequences, our results suggest that four of five virus sequences reported, those being SINV-14, SINV-15, SINV-16 and SINV-17, may be associated to viruses actively replicating in the ant S. invicta.

Conclusions

The present study expands our knowledge about viral diversity associated with S. invicta in introduced areas with potential to be used as biological control agents, which will require further biological characterization.

The Genome Segments of Bluetongue Virus Differ in Copy Number in a Host-Specific Manner

The variation in viral gene frequencies remains a largely unexplored aspect of within-host genetics. This phenomenon is often considered to be specific to multipartite viruses. Multipartite viruses have segmented genomes, but in contrast to segmented viruses, their segments are each encapsidated alone in a virion. A main hypothesis explaining the evolution of multipartism is that, compared to segmented viruses, it facilitates the regulation of segment abundancy, and the genes the segments carry, within a host. These differences in gene frequencies could allow for expression regulation. Here, we show that wild populations of a segmented virus, bluetongue virus (BTV), also present unequal segment frequencies. BTV cycles between ruminants and Culicoides biting midges. As expected from a role in expression regulation, segment frequencies tended to show specific values that differed between ruminants and midges. Our results expand previous knowledge on gene frequency variation and call for studies on its role and conservation beyond multipartite viruses.

Genome segmentation is mainly thought to facilitate reassortment. Here, we show that segmentation can also allow differences in segment abundance in populations of bluetongue virus (BTV). BTV has a genome consisting in 10 segments, and its cycle primarily involves periodic alternation between ruminants and Culicoides biting midges. We have developed a reverse transcription-quantitative PCR (RT-qPCR) approach to quantify each segment in wild BTV populations sampled in both ruminants and midges during an epizootic. Segment frequencies deviated from equimolarity in all hosts. Interestingly, segment frequencies were reproducible and distinct between ruminants and biting midges. Beyond a putative regulatory role in virus expression, this phenomenon could lead to different evolution rates between segments.

IMPORTANCE The variation in viral gene frequencies remains a largely unexplored aspect of within-host genetics. This phenomenon is often considered to be specific to multipartite viruses. Multipartite viruses have segmented genomes, but in contrast to segmented viruses, their segments are each encapsidated alone in a virion. A main hypothesis explaining the evolution of multipartism is that, compared to segmented viruses, it facilitates the regulation of segment abundancy, and the genes the segments carry, within a host. These differences in gene frequencies could allow for expression regulation. Here, we show that wild populations of a segmented virus, bluetongue virus (BTV), also present unequal segment frequencies. BTV cycles between ruminants and Culicoides biting midges. As expected from a role in expression regulation, segment frequencies tended to show specific values that differed between ruminants and midges. Our results expand previous knowledge on gene frequency variation and call for studies on its role and conservation beyond multipartite viruses.