Contrasting selective patterns across the segmented genome of bluetongue virus in a global reassortment hotspot

Analyses and Insights into Genetic Reassortment and Natural Selection as Key Drivers of Piscine orthoreovirus Evolution

Piscine orthoreovirus (PRV) is a pathogen that causes heart and skeletal muscle inflammation in Salmo salar and has also been linked to circulatory disorders in other farmed salmonids, such as Oncorhynchus kisutch and Oncorhynchus mykiss. The virus has a segmented, double-stranded RNA genome, which makes it possible to undergo genetic reassortment and increase its genomic diversity through point mutations. In this study, genetic reassortment in PRV was assessed using the full genome sequences available in public databases. This study used full genome sequences that were concatenated and genome-wide reassortment events, and phylogenetic analyses were performed using the recombination/reassortment detection program version 5 (RDP5 V 5.5) software. Additionally, each segment was aligned codon by codon, and overall mean distance and selection was tested using the Molecular Evolutionary Genetics Analysis X software, version 10.2 (MEGA X version 10.2). The results showed that there were 17 significant reassortment events in 12 reassortant sequences, involving genome exchange between low and highly virulent genotypes. PRV sequences from different salmonid host species did not appear to limit the reassortment. This study found that PRV frequently undergoes reassortment events to increase the diversity of its segmented genome, leading to antigenic variation and increased virulence. This study also noted that to date, no reassortment events have been described between PRV-1 and PRV-3 genotypes. However, the number of complete genomic sequences within each genotype is uneven. This is important because PRV-3 induces cross-protection against PRV-1, making it a potential vaccine candidate.

In Situ Hybridization (RNAscope) Detection of Bluetongue Virus Serotypes 10 and 17 in Experimentally Co-Infected Culicoides sonorensis

Bluetongue virus (BTV) is a segmented, double-stranded RNA virus transmitted by Culicoides biting midges. Infection of domestic and wild ruminants with BTV can result in a devastating disease and significant economic losses. As a virus with a segmented genome, reassortment among the BTV serotypes that have co-infected a host may increase genetic diversity, which can alter BTV transmission dynamics and generate epizootic events. The objective of this study was to determine the extent of dissemination and characterize the tropism of BTV serotypes 10 and 17 in co-infected Culicoides sonorensis. Midges were exposed to both BTV serotypes via blood meal and processed for histologic slides 10 days after infection. An in situ hybridization approach was employed using the RNAscope platform to detect the nucleic acid segment 2 of both serotypes. Observations of the mosaic patterns in which serotypes did not often overlap suggest that co-infection at the cellular level may not be abundant with these two serotypes in C. sonorensis. This could be a consequence of superinfection exclusion. Understanding BTV co-infection and its biological consequences will add an important dimension to the modeling of viral evolution and emergence.

Widespread Reassortment Contributes to Antigenic Shift in Bluetongue Viruses from South Africa

Bluetongue (BT), a viral disease of ruminants, is endemic throughout South Africa, where outbreaks of different serotypes occur. The predominant serotypes can differ annually due to herd immunity provided by annual vaccinations using a live attenuated vaccine (LAV). This has led to both wild-type and vaccine strains co-circulating in the field, potentially leading to novel viral strains due to reassortment and recombination. Little is known about the molecular evolution of the virus in the field in South Africa. The purpose of this study was to investigate the genetic diversity of field strains of BTV in South Africa and to provide an initial assessment of the evolutionary processes shaping BTV genetic diversity in the field. Complete genomes of 35 field viruses belonging to 11 serotypes, collected from different regions of the country between 2011 and 2017, were sequenced. The sequences were phylogenetically analysed in relation to all the BTV sequences available from GenBank, including the LAVs and reference strains, resulting in the analyses and reassortment detection of 305 BTVs. Phylogenomic analysis indicated a geographical selection of the genome segments, irrespective of the serotype. Based on the initial assessment of the current genomic clades that circulate in South Africa, the selection for specific clades is prevalent in directing genome segment reassortment, which seems to exclude the vaccine strains and in multiple cases involves Segment-2 resulting in antigenic shift.

Genome-Wide Reassortment Analysis of Influenza A H7N9 Viruses Circulating in China during 2013-2019

Reassortment with the H9N2 virus gave rise to the zoonotic H7N9 avian influenza virus (AIV), which caused more than five outbreak waves in humans, with high mortality. The frequent exchange of genomic segments between H7N9 and H9N2 has been well-documented. However, the reassortment patterns have not been described and are not yet fully understood. Here, we used phylogenetic analyses to investigate the patterns of intersubtype and intrasubtype/intralineage reassortment across the eight viral segments. The H7N9 virus and its progeny frequently exchanged internal genes with the H9N2 virus but rarely with the other AIV subtypes. Before beginning the intrasubtype/intralineage reassortment analyses, five Yangtze River Delta (YRD A-E) and two Pearl River Delta (PRD A-B) clusters were divided according to the HA gene phylogeny. The seven reset segment genes were also nomenclatured consistently. As revealed by the tanglegram results, high intralineage reassortment rates were determined in waves 2–3 and 5. Additionally, the clusters of PB2 c05 and M c02 were the most dominant in wave 5, which could have contributed to the onset of the largest H7N9 outbreak in 2016–2017. Meanwhile, a portion of the YRD-C cluster (HP H7N9) inherited their PB2, PA, and M segments from the co-circulating YRD-E (LP H7N9) cluster during wave 5. Untanglegram results revealed that the reassortment rate between HA and NA was lower than HA with any of the other six segments. A multidimensional scaling plot revealed a robust genetic linkage between the PB2 and PA genes, indicating that they may share a co-evolutionary history. Furthermore, we observed relatively more robust positive selection pressure on HA, NA, M2, and NS1 proteins. Our findings demonstrate that frequent reassortment, particular reassorted patterns, and adaptive mutations shaped the H7N9 viral genetic diversity and evolution. Increased surveillance is required immediately to better understand the current state of the HP H7N9 AIV.

Bluetongue Research at a Crossroads: Modern Genomics Tools Can Pave the Way to New Insights

Bluetongue virus (BTV) is an arthropod-borne, segmented double-stranded RNA virus that can cause severe disease in both wild and domestic ruminants. BTV evolves via several key mechanisms, including the accumulation of mutations over time and the reassortment of genome segments.Additionally, BTV must maintain fitness in two disparate hosts, the insect vector and the ruminant. The specific features of viral adaptation in each host that permit host-switching are poorly characterized. Limited field studies and experimental work have alluded to the presence of these phenomena at work, but our understanding of the factors that drive or constrain BTV's genetic diversification remains incomplete. Current research leveraging novel approaches and whole genome sequencing applications promises to improve our understanding of BTV's evolution, ultimately contributing to the development of better predictive models and management strategies to reduce future impacts of bluetongue epizootics.

Rotavirus A Genome Segments Show Distinct Segregation and Codon Usage Patterns

Reassortment of the Rotavirus A (RVA) 11-segment dsRNA genome may generate new genome constellations that allow RVA to expand its host range or evade immune responses. Reassortment may also produce phylogenetic incongruities and weakly linked evolutionary histories across the 11 segments, obscuring reassortment-specific epistasis and changes in substitution rates. To determine the co-segregation patterns of RVA segments, we generated time-scaled phylogenetic trees for each of the 11 segments of 789 complete RVA genomes isolated from mammalian hosts and compared the segments’ geodesic distances. We found that segments 4 (VP4) and 9 (VP7) occupied significantly different tree spaces from each other and from the rest of the genome. By contrast, segments 10 and 11 (NSP4 and NSP5/6) occupied nearly indistinguishable tree spaces, suggesting strong co-segregation. Host-species barriers appeared to vary by segment, with segment 9 (VP7) presenting the weakest association with host species. Bayesian Skyride plots were generated for each segment to compare relative genetic diversity among segments over time. All segments showed a dramatic decrease in diversity around 2007 coinciding with the introduction of RVA vaccines. To assess selection pressures, codon adaptation indices and relative codon deoptimization indices were calculated with respect to different host genomes. Codon usage varied by segment with segment 11 (NSP5) exhibiting significantly higher adaptation to host genomes. Furthermore, RVA codon usage patterns appeared optimized for expression in humans and birds relative to the other hosts examined, suggesting that translational efficiency is not a barrier in RVA zoonosis.

Serological Cross-Reactions between Expressed VP2 Proteins from Different Bluetongue Virus Serotypes

Bluetongue (BT) is a severe and economically important disease of ruminants that is widely distributed around the world, caused by the bluetongue virus (BTV). More than 28 different BTV serotypes have been identified in serum neutralisation tests (SNT), which, along with geographic variants (topotypes) within each serotype, reflect differences in BTV outer-capsid protein VP2. VP2 is the primary target for neutralising antibodies, although the basis for cross-reactions and serological variations between and within BTV serotypes is poorly understood. Recombinant BTV VP2 proteins (rVP2) were expressed in Nicotiana benthamiana, based on sequence data for isolates of thirteen BTV serotypes (primarily from Europe), including three 'novel' serotypes (BTV-25, -26 and -27) and alternative topotypes of four serotypes. Cross-reactions within and between these viruses were explored using rabbit anti-rVP2 sera and post BTV-infection sheep reference-antisera, in I-ELISA (with rVP2 target antigens) and SNT (with reference strains of BTV-1 to -24, -26 and -27). Strong reactions were generally detected with homologous rVP2 proteins or virus strains/serotypes. The sheep antisera were largely serotype-specific in SNT, but more cross-reactive by ELISA. Rabbit antisera were more cross-reactive in SNT, and showed widespread, high titre cross-reactions against homologous and heterologous rVP2 proteins in ELISA. Results were analysed and visualised by antigenic cartography, showing closer relationships in some, but not all cases, between VP2 topotypes within the same serotype, and between serotypes belonging to the same 'VP2 nucleotype'.

Genetic and phylogenetic characterization of polycistronic dsRNA segment-10 of bluetongue virus isolates from India between 1985 and 2011

The smallest polycistronic dsRNA segment-10 (S10) of bluetongue virus (BTV) encodes NS3/3A and putative NS5. The S10 sequence data of 46 Indian BTV field isolates obtained between 1985 and 2011 were determined and compared with the cognate sequences of global BTV strains. The largest ORF on S10 encodes NS3 (229 aa) and an amino-terminal truncated form of the protein (NS3A) and a putative NS5 (50-59 aa) due to alternate translation initiation site. The overall mean distance of the global NS3 was 0.1106 and 0.0269 at nt and deduced aa sequence, respectively. The global BTV strains formed four major clusters. The major cluster of Indian BTV strains was closely related to the viruses reported from Australia and China. A minor sub-cluster of Indian BTV strains were closely related to the USA strains and a few of the Indian strains were similar to the South African reference and vaccine strains. The global trait association of phylogenetic structure indicates the evolution of the global BTV S10 was not homogenous but rather represents a moderate level of geographical divergence. There was no evidence of an association between the virus and the host species, suggesting a random spread of the viruses. Conflicting selection pressure on the alternate coding sequences of the S10 was evident where NS3/3A might have evolved through strong purifying (negative) selection and NS5 through a positive selection. The presence of multiple positively selected codons on the putative NS5 may be advantageous for adaptation of the virus though their precise role is unknown.

Phylodynamics and Codon Usage Pattern Analysis of Broad Bean Wilt Virus 2

Broad bean wilt virus 2 (BBWV-2), which belongs to the genus Fabavirus of the family Secoviridae, is an important pathogen that causes damage to broad bean, pepper, yam, spinach and other economically important ornamental and horticultural crops worldwide. Previously, only limited reports have shown the genetic variation of BBWV2. Meanwhile, the detailed evolutionary changes, synonymous codon usage bias and host adaptation of this virus are largely unclear. Here, we performed comprehensive analyses of the phylodynamics, reassortment, composition bias and codon usage pattern of BBWV2 using forty-two complete genome sequences of BBWV-2 isolates together with two other full-length RNA1 sequences and six full-length RNA2 sequences. Both recombination and reassortment had a significant influence on the genomic evolution of BBWV2. Through phylogenetic analysis we detected three and four lineages based on the ORF1 and ORF2 nonrecombinant sequences, respectively. The evolutionary rates of the two BBWV2 ORF coding sequences were 8.895 × 10⁻⁴ and 4.560 × 10⁻⁴ subs/site/year, respectively. We found a relatively conserved and stable genomic composition with a lower codon usage choice in the two BBWV2 protein coding sequences. ENC-plot and neutrality plot analyses showed that natural selection is the key factor shaping the codon usage pattern of BBWV2. Strong correlations between BBWV2 and broad bean and pepper were observed from similarity index (SiD), codon adaptation index (CAI) and relative codon deoptimization index (RCDI) analyses. Our study is the first to evaluate the phylodynamics, codon usage patterns and adaptive evolution of a fabavirus, and our results may be useful for the understanding of the origin of this virus.

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.

Diversity of Transmission Outcomes Following Co-Infection of Sheep with Strains of Bluetongue Virus Serotype 1 and 8

Bluetongue virus (BTV) causes an economically important disease, bluetongue (BT), in susceptible ruminants and is transmitted primarily by species of Culicoides biting midges (Diptera: Ceratopogonidae). Since 2006, northern Europe has experienced multiple incursions of BTV through a variety of routes of entry, including major outbreaks of strains of BTV serotype 8 (BTV-8) and BTV serotype 1 (BTV-1), which overlapped in distribution within southern Europe. In this paper, we examined the variation in response to coinfection with strains of BTV-1 and BTV-8 using an in vivo transmission model involving Culicoides sonorensis, low passage virus strains, and sheep sourced in the United Kingdom. In the study, four sheep were simultaneously infected using BTV-8 and BTV-1 intrathoracically inoculated C. sonorensis and co-infections of all sheep with both strains were established. However, there were significant variations in both the initiation and peak levels of virus RNA detected throughout the experiment, as well as in the infection rates in the C. sonorensis that were blood-fed on experimentally infected sheep at peak viremia. This is discussed in relation to the potential for reassortment between these strains in the field and the policy implications for detection of BTV strains.

Ecological Dynamics Impacting Bluetongue Virus Transmission in North America

Bluetongue virus (BTV) is an arbovirus transmitted to domestic and wild ruminants by certain species of Culicoides midges. The disease resulting from infection with BTV is economically important and can influence international trade and movement of livestock, the economics of livestock production, and animal welfare. Recent changes in the epidemiology of Culicoides-transmitted viruses, notably the emergence of exotic BTV genotypes in Europe, have demonstrated the devastating economic consequences of BTV epizootics and the complex nature of transmission across host-vector landscapes. Incursions of novel BTV serotypes into historically enzootic countries or regions, including the southeastern United States (US), Israel, Australia, and South America, have also occurred, suggesting diverse pathways for the transmission of these viruses. The abundance of BTV strains and multiple reassortant viruses circulating in Europe and the US in recent years demonstrates considerable genetic diversity of BTV strains and implies a history of reassortment events within the respective regions. While a great deal of emphasis is rightly placed on understanding the epidemiology and emergence of BTV beyond its natural ecosystem, the ecological contexts in which BTV maintains an enzootic cycle may also be of great significance. This review focuses on describing our current knowledge of ecological factors driving BTV transmission in North America. Information presented in this review can help inform future studies that may elucidate factors that are relevant to longstanding and emerging challenges associated with prevention of this disease.

Deciphering type-specific neutralizing antibodies to bluetongue virus in goats of Andaman and Nicobar Islands, India

The presence of antibodies to bluetongue virus (BTV) and the viral antigen is reported recently from the Andaman and Nicobar Islands, a group of islands at the juncture of the Bay of Bengal and the Andaman Sea. A retrospective study was conducted to investigate the presence of neutralizing antibodies to different BTV serotypes in the seroconverted goats of the Islands. Thirty six samples out of 186 serum samples tested were selected on the basis of high antibody titre as predicted in an indirect ELISA. Each of the selected serum samples was used for neutralization of six BTV serotypes (BTV-1, BTV-2, BTV-9, BTV-10, BTV-16 and BTV-23), the most commonly reported serotypes in India. Out of 36 serum samples used in the neutralization study, neutralizing antibodies could be determined in 15 samples. The neutralizing antibodies to BTV-10 were found in more number of the serum samples followed by BTV-1, BTV-2 and BTV-23 and BTV-9 and BTV-16. Many of the serum samples could neutralize more than one BTV serotypes indicating possible widespread superinfections by multiple BTV serotypes in goats in the Islands. Majority of the serum samples used in the neutralization study could not neutralize any of the six BTV serotypes commonly reported from India indicating possible circulation of other BTV serotypes yet to confirm. The present study reveals circulation of multiple BTV serotypes in Andaman and Nicobar Islands where there was no such report available earlier. The findings are laudable as the baseline information for further investigations to identify and characterize the virus and competent vectors and for implementing appropriate suitable control strategies for bluetongue in the Islands and the nearby territories.

Origin of Bluetongue Virus Serotype 8 Outbreak in Cyprus, September 2016

In September 2016, clinical signs, indicative of bluetongue, were observed in sheep in Cyprus. Bluetongue virus serotype 8 (BTV-8) was detected in sheep, indicating the first incursion of this serotype into Cyprus. Following virus propagation, Nextera XT DNA libraries were sequenced on the MiSeq instrument. Full-genome sequences were obtained for five isolates CYP2016/01-05 and the percent of nucleotide sequence (% nt) identity between them ranged from 99.92% to 99.95%, which corresponded to a few (2-5) amino acid changes. Based on the complete coding sequence, the Israeli ISR2008/13 (98.42-98.45%) was recognised as the closest relative to CYP2016/01-05. However, the phylogenetic reconstruction of CYP2016/01-05 revealed that the possibility of reassortment in several segments: 4, 7, 9 and 10. Based on the available sequencing data, the incursion BTV-8 into Cyprus most likely occurred from the neighbouring countries (e.g., Israel, Lebanon, Syria, or Jordan), where multiple BTV serotypes were co-circulating rather than from Europe (e.g., France) where a single BTV-8 serotype was dominant. Supporting this hypothesis, atmospheric dispersion modelling identified wind-transport events during July-September that could have allowed the introduction of BTV-8 infected midges from Lebanon, Syria or Israel coastlines into the Larnaca region of Cyprus.

Prospects of Next-Generation Vaccines for Bluetongue

Bluetongue (BT) is a haemorrhagic disease of wild and domestic ruminants with a huge economic worldwide impact on livestock. The disease is caused by BT-virus transmitted by Culicoides biting midges and disease control without vaccination is hardly possible. Vaccination is the most feasible and cost-effective way to minimize economic losses. Marketed BT vaccines are successfully used in different parts of the world. Inactivated BT vaccines are efficacious and safe but relatively expensive, whereas live-attenuated vaccines are efficacious and cheap but are unsafe because of under-attenuation, onward spread, reversion to virulence, and reassortment events. Both manufactured BT vaccines do not enable differentiating infected from vaccinated animals (DIVA) and protection is limited to the respective serotype. The ideal BT vaccine is a licensed, affordable, completely safe DIVA vaccine, that induces quick, lifelong, broad protection in all susceptible ruminant species. Promising vaccine candidates show improvement for one or more of these main vaccine standards. BTV protein vaccines and viral vector vaccines have DIVA potential depending on the selected BTV antigens, but are less effective and likely more costly per protected animal than current vaccines. Several vaccine platforms based on replicating BTV are applied for many serotypes by exchange of serotype dominant outer shell proteins. These platforms based on one BTV backbone result in attenuation or abortive virus replication and prevent disease by and spread of vaccine virus as well as reversion to virulence. These replicating BT vaccines induce humoral and T-cell mediated immune responses to all viral proteins except to one, which could enable DIVA tests. Most of these replicating vaccines can be produced similarly as currently marketed BT vaccines. All replicating vaccine platforms developed by reverse genetics are classified as genetic modified organisms. This implies extensive and expensive safety trails in target ruminant species, and acceptance by the community could be hindered. Nonetheless, several experimental BT vaccines show very promising improvements and could compete with marketed vaccines regarding their vaccine profile, but none of these next generation BT vaccines have been licensed yet.