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.

Untranslated Regions of a Segmented Kindia Tick Virus Genome Are Highly Conserved and Contain Multiple Regulatory Elements for Viral Replication

Novel segmented tick-borne RNA viruses belonging to the group of Jingmenviruses (JMVs) are widespread across Africa, Asia, Europe, and America. In this work, we obtained whole-genome sequences of two Kindia tick virus (KITV) isolates and performed modeling and the functional annotation of the secondary structure of 5' and 3' UTRs from JMV and KITV viruses. UTRs of various KITV segments are characterized by the following points: (1) the polyadenylated 3' UTR; (2) 5' DAR and 3' DAR motifs; (3) a highly conserved 5'-CACAG-3' pentanucleotide; (4) a binding site of the La protein; (5) multiple UAG sites providing interactions with the MSI1 protein; (6) three homologous sequences in the 5' UTR and 3' UTR of segment 2; (7) the segment 2 3' UTR of a KITV/2017/1 isolate, which comprises two consecutive 40 nucleotide repeats forming a Y-3 structure; (8) a 35-nucleotide deletion in the second repeat of the segment 2 3' UTR of KITV/2018/1 and KITV/2018/2 isolates, leading to a modification of the Y-3 structure; (9) two pseudoknots in the segment 2 3' UTR; (10) the 5' UTR and 3' UTR being represented by patterns of conserved motifs; (11) the 5'-CAAGUG-3' sequence occurring in early UTR hairpins. Thus, we identified regulatory elements in the UTRs of KITV, which are characteristic of orthoflaviviruses. This suggests that they hold functional significance for the replication of JMVs and the evolutionary similarity between orthoflaviviruses and segmented flavi-like viruses.

Conserved Sequences in the 5' and 3' Untranslated Regions of Jingmenvirus Group Representatives

The Jingmenvirus group (JVG), with members such as Jingmen tick virus (JMTV), Alongshan virus (ALSV), Yanggou tick virus (YGTV), and Takachi virus (TAKV), is drawing attention due to evidence of it causing disease in humans and its unique genome architecture. In the current work, complete untranslated regions (UTRs) of four strains of ALSV and eight strains of YGTV were obtained. An analysis of these sequences, as well as JVG sequences from GenBank, uncovered several regions within viral UTRs that were highly conserved for all the segments and viruses. Bioinformatics predictions suggested that the UTRs of all the segments of YGTV, ALSV, and JMTV could form similar RNA structures. The most notable feature of these structures was a stable stem-loop with one (5' UTR) or two (3' UTR) AAGU tetraloops on the end of a hairpin.

Distribution and Characterisation of Tick-Borne Flavi-, Flavi-like, and Phenuiviruses in the Chelyabinsk Region of Russia

In this work, we presented data from a two-year study of flavi-, flavi-like, and phenuiviruses circulation in the population of ixodid ticks in the Chelyabinsk region. We isolated three tick-borne encephalitis virus (TBEV) strains from I. persulcatus, which was not detected in the ticks of the genus Dermacentor. The virus prevalence ranged from 0.66% to 2.28%. The Yanggou tick virus (YGTV) is widespread in steppe and forest-steppe zones and is mainly associated with ticks of the genus Dermacentor. We isolated 26 strains from D. reticulatus, D. marginatus, and I. persulcatus ticks in the HAE/CTVM8 tick cell line. The virus prevalence ranged from 1.58% to 4.18% in D. reticulatus, ranged from 0.78% to 3.93% in D. marginatus, and was 0.66% in I. persulcatus. There was combined focus of TBEV and YGTV in the territory of the Chelyabinsk region. The Alongshan virus (ALSV) was found to be associated with I. persulcatus ticks and is spread in forest zone. We detected 12 amplicons and isolated 7 strains of ALSV in tick cells. The virus prevalence ranged from 1.13% to 6.00%. The phlebovirus Gomselga and unclassified phenuivirus Stavropol were associated with I. persulcatus and D. reticulatus ticks, respectively. Virus prevalence of the unclassified phenuivirus Stavropol in the Chelyabinsk region is lower than that in neighbouring regions.

Molecular characterisation of a new tenuivirus from Festuca sp

A novel virus with a quadruple genome of negative-sense, single-stranded RNA was identified by high-throughput sequencing (HTS) in a grass sample from Saxony-Anhalt, Germany, and tentatively called Festuca stripe-associated virus (FSaV). The genome of FSaV consists of four segments and a total of 16,535 nucleotides (nt) which encode seven open reading frames (ORF). FSaV shares highest nt identity (between 72.84% to 80.74%) to Iranian wheat stripe virus (IWSV) and rice hoja blanca virus (RHBV). Additionally, pairwise comparisons between the amino acid sequences of the ORFs on the genome of FSaV and the corresponding ones on the genomes of the members of the Tenuvirus genus showed that FSaV shared 83.17% and 90.85% (amino acid) aa identity to IWSV. Moreover, the non-coding intergenic regions (ncIR) shared only between 49.5% to 60.87% nt identity to the corresponding regions on the IWSV genome. Based on the ICTV species demarcation, the results suggest that FSaV may represent a new species of the genus Tenuivirus. Plastid sequence analysis of the HTS data showed that the original host is a member of the genus Festuca most likely the species Festuca pratensis.

Geographical and Tick-Dependent Distribution of Flavi-Like Alongshan and Yanggou Tick Viruses in Russia

The genus Flavivirus includes related, unclassified segmented flavi-like viruses, two segments of which have homology with flavivirus RNA-dependent RNA polymerase NS5 and RNA helicase-protease NS3. This group includes such viruses as Jingmen tick virus, Alongshan virus, Yanggou tick virus and others. We detected the Yanggou tick virus in Dermacentor nuttalli and Dermacentor marginatus ticks in two neighbouring regions of Russia. The virus prevalence ranged from 0.5% to 8.0%. We detected RNA of the Alongshan virus in 44 individuals or pools of various tick species in eight regions of Russia. The virus prevalence ranged from 0.6% to 7.8%. We demonstrated the successful replication of the Yanggou tick virus and Alongshan virus in IRE/CTVM19 and HAE/CTVM8 tick cell lines without a cytopathic effect. According to the phylogenetic analysis, we divided the Alongshan virus into two groups: an Ixodes persulcatus group and an Ixodes ricinus group. In addition, the I. persulcatus group can be divided into European and Asian subgroups. We found amino acid signatures specific to the I. ricinus and I. persulcatus groups and also distinguished between the European and Asian subgroups of the I. persulcatus group.

In Vitro Reassortment between Endemic Bluetongue Viruses Features Global Shifts in Segment Frequencies and Preferred Segment Combinations

Bluetongue virus (BTV) is an arthropod-borne pathogen that is associated with sometimes severe disease in both domestic and wild ruminants. Predominantly transmitted by Culicoides spp. biting midges, BTV is composed of a segmented, double-stranded RNA genome. Vector expansion and viral genetic changes, such as reassortment between BTV strains, have been implicated as potential drivers of ongoing BTV expansion into previously BTV-free regions. We used an in vitro system to investigate the extent and flexibility of reassortment that can occur between two BTV strains that are considered enzootic to the USA, BTV-2 and BTV-10. Whole genome sequencing (WGS) was coupled with plaque isolation and a novel, amplicon-based sequencing approach to quantitate the viral genetic diversity generated across multiple generations of in vitro propagation. We found that BTV-2 and BTV-10 were able to reassort across multiple segments, but that a preferred BTV-2 viral backbone emerged in later passages and that certain segments were more likely to be found in reassortant progeny. Our findings indicate that there may be preferred segment combinations that emerge during BTV reassortment. Moreover, our work demonstrates the usefulness of WGS and amplicon-based sequencing approaches to improve understanding of the dynamics of reassortment among segmented viruses such as BTV.

Frustration and fidelity in influenza genome assembly

The genome of the influenza virus consists of eight distinct single-stranded RNA segments, each encoding proteins essential for the viral life cycle. When the virus infects a host cell, these segments must be replicated and packaged into new budding virions. The viral genome is assembled with remarkably high fidelity: experiments reveal that most virions contain precisely one copy of each of the eight RNA segments. Cell-biological studies suggest that genome assembly is mediated by specific reversible and irreversible interactions between the RNA segments and their associated proteins. However, the precise inter-segment interaction network remains unresolved. Here, we computationally predict that tree-like irreversible interaction networks guarantee high-fidelity genome assembly, while cyclic interaction networks lead to futile or frustrated off-pathway products. We test our prediction against multiple experimental datasets. We find that tree-like networks capture the nearest-neighbour statistics of RNA segments in packaged virions, as observed by electron tomography. Just eight tree-like networks (of a possible 262 144) optimally capture both the nearest-neighbour data and independently measured RNA–RNA binding and co-localization propensities. These eight do not include the previously proposed hub-and-spoke and linear networks. Rather, each predicted network combines hub-like and linear features, consistent with evolutionary models of interaction gain and loss.

Trans-Acting RNA-RNA Interactions in Segmented RNA Viruses

RNA viruses represent a large and important group of pathogens that infect a broad range of hosts. Segmented RNA viruses are a subclass of this group that encode their genomes in two or more molecules and package all of their RNA segments in a single virus particle. These divided genomes come in different forms, including double-stranded RNA, coding-sense single-stranded RNA, and noncoding single-stranded RNA. Genera that possess these genome types include, respectively, Orbivirus (e.g., Bluetongue virus), Dianthovirus (e.g., Red clover necrotic mosaic virus) and Alphainfluenzavirus (e.g., Influenza A virus). Despite their distinct genomic features and diverse host ranges (i.e., animals, plants, and humans, respectively) each of these viruses uses trans-acting RNA–RNA interactions (tRRIs) to facilitate co-packaging of their segmented genome. The tRRIs occur between different viral genome segments and direct the selective packaging of a complete genome complement. Here we explore the current state of understanding of tRRI-mediated co-packaging in the abovementioned viruses and examine other known and potential functions for this class of RNA–RNA interaction.