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

A Distinct Arabidopsis Latent Virus 1 Isolate Was Found in Wild Brassica hirta Plants and Bees, Suggesting the Potential Involvement of Pollinators in Virus Spread

During our search for aphid-pathogenic viruses, a comovirus was isolated from wild asymptomatic Brassica hirta (white mustard) plants harboring a dense population of Brevicoryne brassicae aphids. The transmission-electron-microscopy visualization of purified virions revealed icosahedral particles. The virus was mechanically transmitted to plants belonging to Brassicaceae, Solanaceae, Amaranthaceae, and Fabaceae families, showing unique ringspot symptoms only on B. rapa var. perviridis plants. The complete viral genome, comprised of two RNA segments, was sequenced. RNA1 and RNA2 contained 5921 and 3457 nucleotides, respectively, excluding the 3' terminal poly-adenylated tails. RNA1 and RNA2 each had one open-reading frame encoding a polyprotein of 1850 and 1050 amino acids, respectively. The deduced amino acids at the Pro-Pol region, delineated between a conserved CG motif of 3C-like proteinase and a GDD motif of RNA-dependent RNA polymerase, shared a 96.5% and 90% identity with the newly identified Apis mellifera-associated comovirus and Arabidopsis latent virus 1 (ArLV1), respectively. Because ArLV1 was identified early in 2018, the B. hirta comovirus was designated as ArLV1-IL-Bh. A high-throughput-sequencing-analyses of the extracted RNA from managed honeybees and three abundant wild bee genera, mining bees, long-horned bees, and masked bees, sampled while co-foraging in a Mediterranean ecosystem, allowed the assembly of ArLV1-IL-Bh, suggesting pollinators' involvement in comovirus spread in weeds.

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.

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.

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.

Temporal Dynamic of the Ratio between Monopartite Begomoviruses and Their Associated Betasatellites in Plants, and Its Modulation by the Viral Gene βC1

The begomovirus-betasatellite complex constantly threatens crops in Asia. However, the quantitative relationship between begomoviruses and betasatellites remains largely unknown. The quantities of tobacco curly shoot virus (TbCSV) and its betasatellite (TbCSB) and their ratio varied significantly in initial infection, and thereafter, the ratio tended to become constant. The TbCSB/TbCSV ratio in agrobacteria inoculum significantly affected that in plants in the initial infection but not thereafter. Null-mutation of βC1 that encodes a multifunctional protein important for pathogenesis in TbCSB significantly reduced the TbCSB/TbCSV ratio in plants. Viral inoculum plants with higher TbCSB/TbCSV ratios promoted whitefly transmission of the virus. The expression of AV1 encoded by TbCSV, βC1 encoded by TbCSB and the βC1/AV1 ratio varied significantly in the initial infection and thereafter the ratio tended to become constant. Additionally, the temporal dynamics of the ratio between another begomovirus and its betasatellite was similar to that of TbCSV and was positively regulated by βC1. These results indicate that the ratio between monopartite begomoviruses and betasatellites tend to become constant as infection progresses, and is modulated by βC1, but a higher betasatellite/begomovirus ratio in virally inoculated plants promotes virus transmission by whiteflies. Our findings provide novel insights into the association between begomoviruses and betasatellites.

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.