The Evolutionary Dynamics of Bluetongue Virus

Bluetongue virus (BTV) is a midge-borne member of the genus Orbivirus that causes an eponymous debilitating livestock disease of great agricultural impact and which has expanded into Europe in recent decades. Reassortment among the ten segments comprising the double-stranded (ds) RNA genome of BTV has played an important role in generating the epidemic strains of this virus in Europe. In this study, we investigated the dynamics of BTV genome segment evolution utilizing time-structured data sets of complete sequences from four segments, totalling 290 sequences largely sampled from ruminant hosts. Our analysis revealed that BTV genome segments generally evolve under strong purifying selection and at substitution rates that are generally lower (mean rates of ~0.5–7 × 10⁻⁴ nucleotide substitutions per site, per year) than vector-borne positive-sense viruses with single-strand (ss) RNA genomes. These also represent the most robust estimates of the nucleotide substitution rate in a dsRNA virus generated to date. Additionally, we determined that patterns of geographic structure and times to most recent common ancestor differ substantially between each segment, including a relatively recent origin for the diversity of segment 10 within the past millennium. Together, these findings demonstrate the effect of reassortment to decouple the evolutionary dynamics of BTV genome segments.     

Heavily glycosylated, highly fit SIVMne variants continue to diversify and undergo selection after transmission to a new host and they elicit early antibody dependent cellular responses but delayed neutralizing antibody responses

The evolutionary success of La Crosse virus (LACV, family Bunyaviridae) is due to its ability to adapt to changing conditions through intramolecular genetic changes and segment reassortment. Vertical transmission of LACV in mosquitoes increases the potential for segment reassortment. Studies were conducted to determine if segment reassortment was occurring in naturally infected Aedes triseriatus from Wisconsin and Minnesota in 2000, 2004, 2006 and 2007. Mosquito eggs were collected from various sites in Wisconsin and Minnesota. They were reared in the laboratory and adults were tested for LACV antigen by immunofluorescence assay. RNA was isolated from the abdomen of infected mosquitoes and portions of the small (S), medium (M) and large (L) viral genome segments were amplified by RT-PCR and sequenced. Overall, the viral sequences from 40 infected mosquitoes and 5 virus isolates were analyzed. Phylogenetic and linkage disequilibrium analyses revealed that approximately 25% of infected mosquitoes and viruses contained reassorted genome segments, suggesting that LACV segment reassortment is frequent in nature.     

Analysis of mixed infection of sheep with bluetongue virus serotypes 10 and 17: evidence for genetic reassortment in the vertebrate host.

Two seronegative sheep were infected intravenously with 10(9) PFU each of bluetongue virus (BTV) serotype 10 and BTV serotype 17. One animal experienced a mild bluetongue-like disease, and both experienced a short-duration viremia and developed neutralizing immune responses to both virus serotypes. Progeny virus was isolated from venous blood from each animal by using conditions in which reassortment could not have occurred during isolation. Electropherotypes were determined for the progeny viruses from the infected sheep, yielding strikingly similar results for the two animals. In both sheep, serotype 10 dominated among the progeny, accounting for 92% of the progeny. Serotype 17 was rarely isolated and accounted for 3% of the progeny analyzed. The remaining 5% of the progeny clones were reassortant and derived genome segments from both serotypes 10 and 17. Analysis of the parental origin of genome segments in the small number of reassortant progeny analyzed suggested that selection of specific genome segments may have occurred in the infected sheep. These data indicate that reassortment of genome segments occurs, at low frequency, in sheep mixedly infected with BTV.     

Genetic reassortants for identification of the genome segment coding for the bluetongue virus hemagglutinin.

Two bluetongue virus (BTV) serotypes isolated in Australia and two selected reassortants derived from cells coinfected with these viruses have been used to identify the gene coding for the virus hemagglutinin. The parent viruses had characteristic hemagglutination patterns: BTV type 20 agglutinated sheep erythrocytes only; and BTV type 21 agglutinated sheep, bovine, human, and goose erythrocytes. Analysis of the two virus clones that had reassorted in genes coding for the outer capsid polypeptides demonstrated that hemagglutination and hemagglutination inhibition are functions associated with the outer capsid protein (VP2), which is encoded by genome segment 2.     

Reassortment between Two Serologically Unrelated Bluetongue Virus Strains Is Flexible and Can Involve any Genome Segment

Coinfection of a cell by two different strains of a segmented virus can give rise to a “reassortant” with phenotypic characteristics that might differ from those of the parental strains. Bluetongue virus (BTV) is a double-stranded RNA (dsRNA) segmented virus and the cause of bluetongue, a major infectious disease of livestock. BTV exists as at least 26 different serotypes (BTV-1 to BTV-26). Prompted by the isolation of a field reassortant between BTV-1 and BTV-8, we systematically characterized the process of BTV reassortment. Using a reverse genetics approach, our study clearly indicates that any BTV-1 or BTV-8 genome segment can be rescued in the heterologous “backbone.” To assess phenotypic variation as a result of reassortment, we examined viral growth kinetics and plaque sizes in in vitro experiments and virulence in an experimental mouse model of bluetongue disease. The monoreassortants generated had phenotypes that were very similar to those of the parental wild-type strains both in vitro and in vivo. Using a forward genetics approach in cells coinfected with BTV-1 and BTV-8, we have shown that reassortants between BTV-1 and BTV-8 are generated very readily. After only four passages in cell culture, we could not detect wild-type BTV-1 or BTV-8 in any of 140 isolated viral plaques. In addition, most of the isolated reassortants contained heterologous VP2 and VP5 structural proteins, while only 17% had homologous VP2 and VP5 proteins. Our study has shown that reassortment in BTV is very flexible, and there is no fundamental barrier to the reassortment of any genome segment. Given the propensity of BTV to reassort, it is increasingly important to have an alternative classification system for orbiviruses.     

Disruption of Specific RNA-RNA Interactions in a Double-Stranded RNA Virus Inhibits Genome Packaging and Virus Infectivity

Bluetongue virus (BTV) causes hemorrhagic disease in economically important livestock. The BTV genome is organized into ten discrete double-stranded RNA molecules (S1-S10) which have been suggested to follow a sequential packaging pathway from smallest to largest segment during virus capsid assembly. To substantiate and extend these studies, we have investigated the RNA sorting and packaging mechanisms with a new experimental approach using inhibitory oligonucleotides. Putative packaging signals present in the 3’untranslated regions of BTV segments were targeted by a number of nuclease resistant oligoribonucleotides (ORNs) and their effects on virus replication in cell culture were assessed. ORNs complementary to the 3’ UTR of BTV RNAs significantly inhibited virus replication without affecting protein synthesis. Same ORNs were found to inhibit complex formation when added to a novel RNA-RNA interaction assay which measured the formation of supramolecular complexes between and among different RNA segments. ORNs targeting the 3’UTR of BTV segment 10, the smallest RNA segment, were shown to be the most potent and deletions or substitution mutations of the targeted sequences diminished the RNA complexes and abolished the recovery of viable viruses using reverse genetics. Cell-free capsid assembly/RNA packaging assay also confirmed that the inhibitory ORNs could interfere with RNA packaging and further substitution mutations within the putative RNA packaging sequence have identified the recognition sequence concerned. Exchange of 3’UTR between segments have further demonstrated that RNA recognition was segment specific, most likely acting as part of the secondary structure of the entire genomic segment. Our data confirm that genome packaging in this segmented dsRNA virus occurs via the formation of supramolecular complexes formed by the interaction of specific sequences located in the 3’ UTRs. Additionally, the inhibition of packaging in-trans with inhibitory ORNs suggests this that interaction is a bona fide target for the design of compounds with antiviral activity.     

Phylogenetic comparison of the S3 gene of United States prototype strains of bluetongue virus with that of field isolates from California.

To better define the molecular epidemiology of bluetongue virus (BTV) infection, the genetic characteristics and phylogenetic relationships of the S3 genes of the five U.S. prototype strains of BTV, the commercially available serotype 10 modified live virus vaccine, and 18 field isolates of BTV serotypes 10, 11, 13, and 17 obtained in California during 1980, 1981, 1989, and 1990 were determined. With the exception of the S3 gene of the U.S. prototype strain of BTV serotype 2 (BTV 2), these viruses had an overall sequence homology of between 95 and 100%. Phylogenetic analyses segregated the prototype U.S. BTV 2 strain to a unique branch (100% bootstrap value), whereas the rest of the viruses clustered in two main monophyletic groups that were not correlated with their serotype, year of isolation, or geographical origin. The lack of consistent association between S3 gene sequence and virus serotype likely is a consequence of reassortment of BTV gene segments during natural mixed infections of vertebrate and invertebrate hosts. The prototype strain of BTV 13, which is considered an introduction to the U.S. like BTV 2, presents an S3 gene which is highly homologous to those of some isolates of BTV 10 and especially to that of the vaccine strain. This finding strongly suggests that the U.S. prototype strain of BTV 13 is a natural reassortant. The different topologies of the phylogenetic trees of the L2 and S3 genes of the various viruses indicate that these two genome segments evolve independently. We conclude that the S3 gene segment of populations of BTV in California is formed by different consensus sequences which cocirculate and which cannot be grouped by serotype.     

Widespread Recombination,Reassortment, and Transmission of Unbalanced Compound Viral Genotypes in Natural Arenavirus Infections

Arenaviruses are one of the largest families of human hemorrhagic fever viruses and are known to infect both mammals and snakes. Arenaviruses package a large (L) and small (S) genome segment in their virions. For segmented RNA viruses like these, novel genotypes can be generated through mutation, recombination, and reassortment. Although it is believed that an ancient recombination event led to the emergence of a new lineage of mammalian arenaviruses, neither recombination nor reassortment has been definitively documented in natural arenavirus infections. Here, we used metagenomic sequencing to survey the viral diversity present in captive arenavirus-infected snakes. From 48 infected animals, we determined the complete or near complete sequence of 210 genome segments that grouped into 23 L and 11 S genotypes. The majority of snakes were multiply infected, with up to 4 distinct S and 11 distinct L segment genotypes in individual animals. This S/L imbalance was typical: in all cases intrahost L segment genotypes outnumbered S genotypes, and a particular S segment genotype dominated in individual animals and at a population level. We corroborated sequencing results by qRT-PCR and virus isolation, and isolates replicated as ensembles in culture. Numerous instances of recombination and reassortment were detected, including recombinant segments with unusual organizations featuring 2 intergenic regions and superfluous content, which were capable of stable replication and transmission despite their atypical structures. Overall, this represents intrahost diversity of an extent and form that goes well beyond what has been observed for arenaviruses or for viruses in general. This diversity can be plausibly attributed to the captive intermingling of sub-clinically infected wild-caught snakes. Thus, beyond providing a unique opportunity to study arenavirus evolution and adaptation, these findings allow the investigation of unintended anthropogenic impacts on viral ecology, diversity, and disease potential.     

Evolution of novel reassortant A/H3N2 influenza viruses in North American swine and humans, 2009-2011

Novel H3N2 influenza viruses (H3N2v) containing seven genome segments from swine lineage triple-reassortant H3N2 viruses and a 2009 pandemic H1N1 (H1N1pdm09) matrix protein segment (pM) were isolated from 12 humans in the United States between August and December 2011. To understand the evolution of these novel H3N2 viruses in swine and humans, we undertook a phylogenetic analysis of 674 M sequences and 388 HA and NA sequences from influenza viruses isolated from North American swine during 2009-2011, as well as HA, NA, and M sequences from eight H3N2v viruses isolated from humans. We identified 34 swine influenza viruses (termed rH3N2p) with the same combination of H3, N2, and pM segments as the H3N2v viruses isolated from humans. Notably, these rH3N2p viruses were generated in swine via reassortment events between H3N2 viruses and the pM segment approximately 4 to 10 times since 2009. The pM segment has also reassorted with multiple distinct lineages of H1 virus, especially H1δ viruses. Importantly, the N2 segment of all H3N2v viruses isolated from humans is derived from a genetically distinct N2 lineage that has circulated in swine since being acquired by reassortment with seasonal human H3N2 viruses in 2001-2002, rather than from the N2 that is associated with the 1998 H3N2 swine lineage. The identification of this N2 variant may have implications for influenza vaccine design and the potential pandemic threat of H3N2v to human age groups with differing levels of prior exposure and immunity.     

Genome segment reassortment between two serotypes of bluetongue virus in a natural host.

The occurrence of genome segment reassortment between two antigenically related orbiviruses was demonstrated in cattle. Individual virus clones were isolated by cell culture plaque assays directly from the blood of a calf infected with two serotypes of bluetongue virus. The majority (89%) of progeny viruses isolated from the calf represented reassortant viruses. A minimum of six genome segments participated in reassortment, with 16 unique reassortant constellations being identified. Such genome segment reassortment between unique, though antigenically related, orbiviruses has undoubtedly played a major role in generating the extensive phenotypic and genotypic diversity that is characteristic of this serogroup.     

GEOGRAPHIC DIFFERENCES IN SEXUAL REASSORTMENT IN RNA PHAGE

The genetic structure of natural bacteriophage populations is poorly understood. Recent metagenomic studies suggest that phage biogeography is characterized by frequent migration. Using virus samples mostly isolated in Southern California, we recently showed that very little population structure exists in segmented RNA phage of the Cystoviridae family due to frequent segment reassortment (sexual genetic mixis) between unrelated virus individuals. Here we use a larger genetic dataset to examine the structure of Cystoviridae phage isolated from three geographic locations in Southern New England. We document extensive natural variation in the physical sizes of RNA genome segments for these viruses. In addition, consistent with earlier findings, our phylogenetic analyses and calculations of linkage disequilibrium (LD) show no evidence of within‐segment recombination in wild populations. However, in contrast to the prior study, our analysis finds that reassortment of segments between individual phage plays a lesser role among cystoviruses sampled in New England, suggesting that the evolutionary importance of genetic mixis in Cystoviridae phage may vary according to geography. We discuss possible explanations for these conflicting results across the studies, such as differing local ecology and its impact on phage growth, and geographic differences in selection against hybrid phage genotypes.     

Multiple genotypes of influenza B virus circulated between 1979 and 2003

The segmented genome of influenza B virus allows exchange of gene segments between cocirculating strains. Through this process of reassortment, diversity is generated by the mixing of genes between viruses that differ in one or more gene segments. Phylogenetic and evolutionary analyses of all 11 genes of 31 influenza B viruses isolated from 1979 to 2003 were used to study the evolution of whole genomes. All 11 genes diverged into two new lineages prior to 1987. All genes except the NS1 gene were undergoing linear evolution, although the rate of evolution and the degree to which nucleotide changes translated into amino acid changes varied between lineages and by gene. Frequent reassortment generated 14 different genotypes distinct from the gene constellation of viruses circulating prior to 1979. Multiple genotypes cocirculated in some locations, and a sequence of reassortment events over time could not be established. The surprising diversity of the viruses, unrestricted mixing of lineages, and lack of evidence for coevolution of gene segments do not support the hypothesis that the reassortment process is driven by selection for functional differences.     

Evolution and Phylogenetic Analysis of Full-Length VP3 Genes of Eastern Mediterranean Bluetongue Virus Isolates

Bluetongue virus (BTV) is the ‘type’ species of the genus Orbivirus within the family Reoviridae. The BTV genome is composed of ten linear segments of double-stranded RNA (dsRNA), each of which codes for one of ten distinct viral proteins. Previous phylogenetic comparisons have evaluated variations in genome segment 3 (Seg-3) nucleotide sequence as way to identify the geographical origin (different topotypes) of BTV isolates. The full-length nucleotide sequence of genome Seg-3 was determined for thirty BTV isolates recovered in the eastern Mediterranean region, the Balkans and other geographic areas (Spain, India, Malaysia and Africa). These data were compared, based on molecular variability, positive-selection-analysis and maximum-likelihood phylogenetic reconstructions (using appropriate substitution models) to 24 previously published sequences, revealing their evolutionary relationships. These analyses indicate that negative selection is a major force in the evolution of BTV, restricting nucleotide variability, reducing the evolutionary rate of Seg-3 and potentially of other regions of the BTV genome. Phylogenetic analysis of the BTV-4 strains isolated over a relatively long time interval (1979–2000), in a single geographic area (Greece), showed a low level of nucleotide diversity, indicating that the virus can circulate almost unchanged for many years. These analyses also show that the recent incursions into south-eastern Europe were caused by BTV strains belonging to two different major-lineages: representing an ‘eastern’ (BTV-9, -16 and -1) and a ‘western’ (BTV-4) group/topotype. Epidemiological and phylogenetic analyses indicate that these viruses originated from a geographic area to the east and southeast of Greece (including Cyprus and the Middle East), which appears to represent an important ecological niche for the virus that is likely to represent a continuing source of future BTV incursions into Europe.     

Reassortment Patterns in Swine Influenza Viruses

Three human influenza pandemics occurred in the twentieth century, in 1918, 1957, and 1968. Influenza pandemic strains are the results of emerging viruses from non-human reservoirs to which humans have little or no immunity. At least two of these pandemic strains, in 1957 and in 1968, were the results of reassortments between human and avian viruses. Also, many cases of swine influenza viruses have reportedly infected humans, in particular, the recent H1N1 influenza virus of swine origin, isolated in Mexico and the United States. Pigs are documented to allow productive replication of human, avian, and swine influenza viruses. Thus it has been conjectured that pigs are the “mixing vessel” that create the avian-human reassortant strains, causing the human pandemics. Hence, studying the process and patterns of viral reassortment, especially in pigs, is a key to better understanding of human influenza pandemics. In the last few years, databases containing sequences of influenza A viruses, including swine viruses, collected since 1918 from diverse geographical locations, have been developed and made publicly available. In this paper, we study an ensemble of swine influenza viruses to analyze the reassortment phenomena through several statistical techniques. The reassortment patterns in swine viruses prove to be similar to the previous results found in human viruses, both in vitro and in vivo, that the surface glycoprotein coding segments reassort most often. Moreover, we find that one of the polymerase segments (PB1), reassorted in the strains responsible for the last two human pandemics, also reassorts frequently.     

Intragenic Recombination as a Mechanism of Genetic Diversity in Bluetongue Virus

Bluetongue (BT), caused by Bluetongue virus (BTV), is an economically important disease affecting sheep, deer, cattle, and goats. Since 1998, a series of BT outbreaks have spread across much of southern and central Europe. To study why the epidemiology of the virus happens to change, it is important to fully know the mechanisms resulting in its genetic diversity. Gene mutation and segment reassortment have been considered as the key forces driving the evolution of BTV. However, it is still unknown whether intragenic recombination can occur and contribute to the process in the virus. We present here several BTV groups containing mosaic genes to reveal that intragenic recombination can take place between the virus strains and play a potential role in bringing novel BTV lineages.Bluetongue (BT) is an economically significant disease that seriously threatens sheep, some species of deer, and to a lesser extent cattle and goats. As a vector-borne viral disease of ruminants, BT is endemic in tropical and subtropical countries (46). However, a series of BT outbreaks have spread across much of southern and central Europe since 1998 (29). Thus, it is of great importance to fully understand the molecular basis driving the change of its epidemiology so as to prevent or limit future BT pandemics.Bluetongue virus (BTV), the pathogen of BT, belongs to the Orbivirus genus of the Reoviridae family (46). The virus has a segmented double-stranded RNA (dsRNA) genome that is packaged in a nonenveloped, icosahedral particle (46). Its 10 dsRNA segments encode 11 proteins, VP1 to VP7 (encoded by segments 1, 2, 3, 4, 6, 9, and 7, respectively), NS1 to SN3 (encoded by segments 5, 8, and 10, respectively), and NS3A (encoded by segment 10) (46). Two structural proteins, VP2 and VP5, form the outer layer of the virion particle and are responsible for cell attachment and virus entry (18, 31, 32), neutralizing epitope (14, 21), and virus virulence (36). Both of them are highly variable and generate 24 serotypes of the virus (44). The inner layers contain VP1, VP3, VP4, VP6, and VP7, and form the “core” of the BTV capsid. VP1 and VP6 are involved in RNA replication as the RNA-dependent RNA polymerase (54) and helicase/NTPase, respectively (49). VP7 forms the surface of the core and functions during the entry of the core into insect cells (44) and also can react with “core neutralizing” antibodies as a major serogroup-specific antigen (32, 44). These core proteins and two nonstructural proteins, NS1 and NS2, are thought to be relatively conservative, so that antigenic cross-reaction can take place between different BTV strains and serotypes, whereas NS3/N3a is more variable than the other nonstructural or core proteins (46).The genetic diversity and variation in sequences of different BTV genome segments were initially identified by RNA oligonucleotide fingerprint analysis of BTV field samples (47). Until now, reassortment and dynamic gene mutation, regarded as the key factors responsible for the genetic diversity of BTV, have been studied in details (46). The two mechanisms can result in both genetic drift and genetic shift and contribute to BTV evolution (47). It has been revealed that high-frequency genome segment reassortment occurs readily between different BTV serotypes (16). Thus, segment reassortment is an important factor in generation of genetic diversity in orbivirus populations in nature (45). In addition, it has been shown that homologous recombination can also play a role in the genetic diversity and evolution of some RNA viruses (24, 33) and bring on virulent variants of these viruses at last (8, 56). Although homologous recombination has been observed in rotavirus, a member of the Reoviridae (39, 40), it is still unknown whether the intragenic recombination can occur and play a role in the generation of genetic diversity in orbivirus populations.To determine whether homologous recombination shaped the evolution of BTV and to provide some insights into the recombination itself in the virus, we analyzed roughly 690 complete segments of BTV deposited in GenBank to see whether some of them underwent intragenic recombination event. Several BTV groups isolated at different time points and in different countries were found containing the same (or similar) mosaic segments, demonstrating that intragenic recombination had occurred in the field and that these viruses with mosaic segments had become prevailing strains. That is, intragenic recombination can play a potential role in generating genetic diversity of BTV and exert its influence on the change of BTV epidemiology.     

Reassortment Networks and the evolution of pandemic H1N1 swine-origin influenza

Prior research developed Reassortment Networks to reconstruct the evolution of segmented viruses under both reassortment and mutation. We report their application to the swine-origin pandemic H1N1 virus (S-OIV). A database of all influenza A viruses, for which complete genome sequences were available in Genbank by October 2009, was created and dynamic programming was used to compute distances between all corresponding segments. A reassortment network was created to obtain the minimum cost evolutionary paths from all viruses to the exemplar S-OIV A/California/04/2009. This analysis took 35 hours on the Cray Extreme Multithreading (XMT) supercomputer, which has special hardware to permit efficient parallelization. Six specific H1N1/H1N2 bottleneck viruses were identified that almost always lie on minimum cost paths to S-OIV. We conjecture that these viruses are crucial to S-OIV evolution and worthy of careful study from a molecular biology viewpoint. In phylogenetics, ancestors are typically medians that have no functional constraints. In our method, ancestors are not inferred, but rather chosen from previously observed viruses along a path of mutation and reassortment leading to the target virus. This specificity and functional constraint render our results actionable for further experiments in vitro and in vivo.     

Prospective of Genomics in Revealing Transmission, Reassortment and Evolution of Wildlife-Borne Avian Influenza A (H5N1) Viruses

The outbreak of highly pathogenic avian influenza (HPAI) H5N1 disease has led to significant loss of poultry and wild life and case fatality rates in humans of 60%. Wild birds are natural hosts for all avian influenza virus subtypes and over120 bird species have been reported with evidence of H5N1 infection. Influenza A viruses possess a segmented RNA genome and are characterized by frequently occurring genetic reassortment events, which play a very important role in virus evolution and the spread of novel gene constellations in immunologically naïve human and animal populations. Phylogenetic analysis of whole genome or sub-genomic sequences is a standard means for delineating genetic variation, novel reassortment events, and surveillance to trace the global transmission pathways. In this paper, special emphasis is given to the transmission and circulation of H5N1 among wild life populations, and to the reassortment events that are associated with inter-host transmission of the H5N1 viruses when they infect different hosts, such as birds, pigs and humans. In addition, we review the inter-subtype reassortment of the viral segments encoding inner proteins between the H5N1 viruses and viruses of other subtypes, such as H9N2 and H6N1. Finally, we highlight the usefulness of genomic sequences in molecular epidemiological analysis of HPAI H5N1 and the technical limitations in existing analytical methods that hinder them from playing a greater role in virological research.     

Sequential packaging of RNA genomic segments during the assembly of Bluetongue virus

Bluetongue virus (BTV), a member of the Orbivirus genus within the Reoviridae family, has a genome of 10 double-stranded RNA segments, with three distinct size classes. Although the packaging of the viral genome is evidently highly specific such that every virus particle contains a set of 10 RNA segments, the order and mechanism of packaging are not understood. In this study we have combined the use of a cell-free in vitro assembly system with a novel RNA–RNA interaction assay to investigate the mechanism of single-stranded (ss) RNAs packaging during nascent capsid assembly. Exclusion of single or multiple ssRNA segments in the packaging reaction or their addition in different order significantly altered the outcome and suggested a particular role for the smallest segment, S10. Our data suggests that genome packaging probably initiates with the smallest segment which triggers RNA–RNA interaction with other smaller segments forming a complex network. Subsequently, the medium to larger size ssRNAs are recruited until the complete genome is packaging into the capsid. The untranslated regions of the smallest RNA segment, S10, is critical for the instigation of this process. We suggest that the selective packaging observed in BTV may also apply to other members of the Reoviridae family.     

Joint Inference of Migration and Reassortment Patterns for Viruses with Segmented Genomes

The structured coalescent allows inferring migration patterns between viral subpopulations from genetic sequence data. However, these analyses typically assume that no genetic recombination process impacted the sequence evolution of pathogens. For segmented viruses, such as influenza, that can undergo reassortment this assumption is broken. Reassortment reshuffles the segments of different parent lineages upon a coinfection event, which means that the shared history of viruses has to be represented by a network instead of a tree. Therefore, full genome analyses of such viruses are complex or even impossible. Although this problem has been addressed for unstructured populations, it is still impossible to account for population structure, such as induced by different host populations, whereas also accounting for reassortment. We address this by extending the structured coalescent to account for reassortment and present a framework for investigating possible ties between reassortment and migration (host jump) events. This method can accurately estimate subpopulation dependent effective populations sizes, reassortment, and migration rates from simulated data. Additionally, we apply the new model to avian influenza A/H5N1 sequences, sampled from two avian host types, Anseriformes and Galliformes. We contrast our results with a structured coalescent without reassortment inference, which assumes independently evolving segments. This reveals that taking into account segment reassortment and using sequencing data from several viral segments for joint phylodynamic inference leads to different estimates for effective population sizes, migration, and clock rates. This new model is implemented as the Structured Coalescent with Reassortment package for BEAST 2.5 and is available at https://github.com/jugne/SCORE.