SEARCH · Search PubMed

Results for “RNA virus evolution”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

PubMed2026

Transcriptome mining and comparative genomics reveal 36 putative novel marafivirus species and conserved evolution of the marafibox regulatory element.

BACKGROUND: Marafiviruses are plant-infecting RNA viruses associated with several economically important crops, but their genomic diversity remains incompletely characterized. OBJECTIVE: This study aimed to identify previously unrecognized marafivirus genomes and investigate their genomic features and evolutionary relationships. METHODS: Publicly available plant transcriptome datasets were systematically mined to detect marafivirus-like sequences. Recovered genomes were analyzed using comparative sequence analysis, phylogenetic reconstruction, and genome organization characterization. RESULTS: A total of 62 marafivirus-like genomes were recovered from 33 independent sources representing diverse plant hosts. Polyprotein-based comparative and phylogenetic analyses grouped these genomes into 36 lineages likely representing novel species. All newly identified viruses clustered within the Marafivirus clade. Genome organization analysis revealed conserved polyprotein architecture and widespread presence of the marafibox promoter element. Conservation of additional open reading frames among closely related isolates aided identification of potentially functional genes. CONCLUSION: These findings substantially expand the known diversity of marafiviruses and demonstrate the effectiveness of transcriptome mining for discovering previously unrecognized plant viruses.

Phylogeny↗

PubMed2025

Analysis of Leishbuviridae from Trypanosomatids.

Over the last decade, considerable progress has been made in unraveling RNA virus diversity. This has contributed to our understanding of the evolution of these viruses, which include emerging zoonotic human pathogens. Current success has been greatly facilitated by the development of next-generation sequencing platforms instrumental for meta-transcriptomic studies. However, due to the rapid evolution of RNA viruses, there are numerous "blind spots" waiting to be explored; one of those is the RNA virome of unicellular eukaryotes. Here, we present the pipeline, which has been successfully used to characterize various types of RNA viruses, including Leishbuviridae (Bunyaviricetes, Hareavirales) in the parasitic flagellates of the family Trypanosomatidae. The pipeline relies on axenic in vitro cell culture and double-stranded RNA enrichment, followed by direct RNA-sequencing. A detailed procedure description starting from the initial total RNA preparation to the final assembly of the viral segments is provided.

High-Throughput Nucleotide Sequencing↗

PubMed2025

A positive-sense single-stranded RNA virus acquired a negative-sense open reading frame through recombination.

Although positive- and negative-sense single-stranded RNA viruses are ubiquitous in nature, there is currently no evidence of recombination or reassortment between viruses with these two major forms of genome organization. Here, we describe the discovery of brine shrimp virga-like virus 1 (BSVV1), a novel positive-sense single-stranded RNA virus with a recombinant genome structure derived from two viral phyla with differing genome organizations. The genome of BSVV1 comprises three open reading frames (ORFs). ORF1 resembles the RNA-dependent RNA polymerase of Ips virga-like virus 1 (a positive-sense RNA virus), while ORF2, transcribed in the positive orientation, is related to the glycoprotein of Hubei bunya-like virus 10 and other negative-sense RNA viruses. The predicted ORF3 was unique to BSVV1 without known homologs identified. The presence of the three protein products was verified by mass spectrometry. Notably, our analysis also revealed that BSVV1 is geographically widespread and found in brine shrimp from at least eight countries on four continents. In addition, BSVV1 was successfully cultured and proliferated to high viral loads during brine shrimp development. In sum, we provide compelling evidence of an ancient recombination event between negative- and positive-sense single-stranded RNA viruses, enriching our understanding of the evolution of genome structures in RNA viruses.

Open Reading Frames↗

PubMed2026

Selection profiles in RNA viruses reflect the characteristics of viruses more than individual proteins.

Proteins that are exposed on the surface of a virus are frequently subject to strong selection to escape from neutralizing antibodies. To investigate whether surface-exposed (SE) and non-exposed (NE) proteins encoded by RNA viruses exhibit different patterns of evolution under selection, we analyzed 244 protein-coding genes from 28 species of RNA viruses representing 15 taxonomic families. First, we show that gene-wide rates of non-synonymous (dN) and synonymous (dS) substitutions do not differentiate between SE and NE proteins. To incorporate variation in substitution rates among codon sites, we inferred the posterior distribution over a fixed grid of dN and dS rates for each alignment. This 'evolutionary fingerprint' provides a common framework for comparing the selection profiles of non-homologous genes. Next, we computed the Wasserstein distance for every pair of fingerprints, which is analogous to amount of work required to reshape one distribution to another. After compensating for differences in genetic variation among alignments, we found a small but significant difference between the fingerprints of SE and NE proteins (PERMANOVA, P = 0.03). However, we observed larger and more significant effects of whether the virus is enveloped (P < 10-5) and the interaction between these factors (P=6.9×10-4). The latter effects were driven by high levels of purifying selection in capsid proteins of Picornaviruses. Furthermore, greater amounts of variation in fingerprints were explained by significant differences among virus families and modes of transmission (P < 10-5). These results imply the pattern of selection on a virus protein is shaped more by characteristics of the virus than the protein itself.

RNA Viruses↗

PubMed2026

Evolution of virulence of a plant RNA virus in developmental stage-structured host populations.

Natural host populations are age-structured, and developmental stages differ in susceptibility and within-host pathogen dynamics, potentially imposing distinct selective pressures on viruses. However, the evolutionary consequences of host age structure remain poorly understood. We experimentally evolved turnip mosaic potyvirus for 5 passages in Arabidopsis thaliana populations spanning 7 demographic regimes, from juvenile- to mature-dominated cohorts. We quantified disease progression, symptom severity, and viral load, cross-inoculated evolved lineages across host stages to construct infection matrices, and performed whole-population sequencing at passages 1 and 5. Disease traits changed markedly with passage, demography, and their interaction. Disease progression evolved faster in older populations, whereas symptom severity was independent of median age, indicating demographic reweighting of virulence components. Viral load increased across passages and positively correlated with severity, linking within-host fitness to symptoms. Cross-inoculation assays revealed a modular infection network: juvenile-evolved lineages specialized on juvenile hosts, whereas lineages from intermediate and older populations were more generalist. Genomically, we detected both parallel and demography-specific adaptations, including recurrent changes in the viral protein VPg (involved in translation, replication, and host interactions) as well as synonymous variants showing consistent or opposing selection across host population stage structures. Overall, host age structure emerges as a major ecological driver of virulence evolution, shaping tradeoffs between disease progression and severity and determining specialization versus generalism. These results integrate phenotypic and genomic responses and suggest that manipulating crop age structure could steer virus evolution toward less damaging outcomes.

Virulence↗

PubMed2026

Discovery of the order 'Quisvirales' redefines the evolution of RNA replication and transcription in the phylum Pisuviricota.

Genome replication in positive-stranded RNA (ssRNA+) viruses is mediated by cognate enzymes, including ubiquitous RNA-dependent RNA polymerase (RdRp). In ssRNA+ viruses with multiple open reading frames (ORFs) in their genomes, replication often is accompanied by synthesis of subgenomic RNAs (transcription) for expression of 3'-proximal ORFs. In addition, all ssRNA+ viruses with genomes larger than ~7 kb encode helicases, linking helicases to RNA genome expansion. Helicases are essential ATPases that unwind nucleic acids and are classified into six recognized superfamilies (SF1-SF6). In the phylum Pisuviricota that includes important pathogens, helicases of SF1-SF3 are integrated into multi-enzyme replicase polyprotein(s) including 3C(-like) protease (3CLpro) and RdRp. Here, large-scale mining of invertebrate metatranscriptomes and targeted genome sequence assembly uncovered six spider-associated ssRNA+ viruses that, based on their conserved 3CLpro-RdRp module in replicase polyproteins, genome size (20-22 kb), and phylogeny, form a family-like cluster in a putative order, named 'Quisvirales'. Quisviruses have similar genome and replicase architectures to enveloped coronaviruses and other nidoviruses. Notably, quisviruses encode ORFs 1a and 1b with predicted -1 programmed ribosomal frameshifting elements in the ORF1a/b overlap region. Using an original mapping approach for detecting chimeric sequencing reads, we obtained evidence that 3'-proximal ORFs are expressed via 5'-coterminal, leader-containing subgenomic RNAs. This suggests that the quisvirus subgenomic RNAs are generated through discontinuous transcription-a mechanism otherwise exclusively found in nidoviruses among the many ssRNA+ virus orders that synthesize subgenomic RNAs. Striking differences between nido- and quisviruses are, however, the RdRp being the only common core ORF1b-encoded enzyme and the replacement of the nidovirus SF1 helicase by a novel superfamily helicase. This quisvirus SF7 helicase, like the Picornavirales SF3 helicase, comprises an AAA+ (ATPase-like) domain typical for ring-forming helicases and thus must play an essential role in replication. The discovery of the order 'Quisvirales' demonstrates that viruses employing large replicase polyproteins of nidovirus-like complexity and discontinuous transcription may have evolved repeatedly from an 3CLpro-RdRp-encoding ancestor.

AAA+/RecA-like ATPase↗

PubMed2026

RNA Virus Diversity, Cross-Species Transmission, and Molecular Constraints in Two Closely Related Rat Species.

Viral infection involves co-evolution with hosts, yet the molecular determinants that constrain viral cross-species transmission remain poorly understood. Here, we established conspecific and heterospecific co-housing models for two closely related rat species, Rattus norvegicus (RN) and Rattus tanezumi (RT), both maintained in laboratory settings for over 10 generations, together with wild-caught RT individuals. Using meta-transcriptomic sequencing and population genomic analyses, we compared their RNA virus profiles and investigated the potential molecular constraints on cross-species viral transmission. From 63 rats, we characterized an extensive RNA virome comprising more than 600 viruses, including 7 zoonotic viruses, 29 viruses with cross-species transmission potential, and 335 novel viruses. Notably, the prevalence of Seoul orthohantavirus (SEOV) was significantly higher in RN than in RT. Population genomic analysis revealed that RN exhibited higher heterozygosity in Itgb3 (the gene encoding the SEOV receptor, β3-integrin) and Tlr7 (the gene encoding the receptor for viral ssRNA, Toll-like receptor 7) compared to RT. These genetic variations likely represent the molecular determinants responsible for the differential susceptibility to SEOV between the two species. Our findings clarify the diversity and prevalence of RNA viruses in closely related rodent species and highlight host genetic barriers that may influence zoonotic spillover risk.

Animals↗

PubMed2026

Nipah virus in the era of global connectivity: molecular evolution, transmission risk, and preparedness strategies.

Nipah virus (NiV) is a highly pathogenic zoonotic RNA virus belonging to the genus Henipavirus within the family Paramyxoviridae, representing a continuing global health concern due to its high case fatality rate and potential for epidemic expansion in the era of increasing international connectivity. The virus demonstrates strong evolutionary adaptability driven by the absence of proofreading mechanisms during RNA replication, enabling genetic diversification that may influence host range, virulence, and transmission dynamics. Molecular pathogenesis of NiV is primarily mediated through interaction of viral glycoproteins with ephrin-B2 and ephrin-B3 receptors, facilitating host cell entry, endothelial damage, and neuroinvasion. Immune evasion facilitated by the action of accessory proteins encoded by the P gene (P, V, W, and C) acts to suppress innate antiviral immunity through the inhibition of interferon induction and JAK/STAT signaling. Human-to-human transmission of Nipah virus remains limited, with epidemiological evidence indicating basic reproduction numbers generally below unity; however, respiratory involvement and healthcare-associated exposure may enhance cluster outbreaks. Global travel, ecological disruption, and fragmented surveillance systems contribute to spillover risk, particularly in South and Southeast Asia where fruit bats of the genus Pteropus serve as natural reservoirs. Despite advances in vaccine technology, including subunit, viral vector, mRNA-based platforms, and monoclonal antibody therapies, no licensed prophylactic or therapeutic agent is currently available for human use. Global preparedness remains challenged by the scarcity of high-containment biosafety facilities, limited research funding, and absence of integrated One Health surveillance networks. Ethical considerations surrounding wildlife population control further complicate disease mitigation strategies. Emerging genomic surveillance, artificial intelligence-assisted predictive modeling, and regional data-sharing frameworks are essential for early detection and response. Strengthening molecular research on viral-host interactions and transmission determinants will be critical for preventing future Nipah virus outbreaks in an increasingly interconnected world.

Genomic surveillance↗

PubMed2026

Creating bottom-up RNA transfer vehicles from synthetic protein assemblies.

Evolution guides biological systems to populate ecological niches, with viruses among the most successful examples of this principle. Viruses evolved over billions of years to efficiently transfer genetic information. Although viruses are highly diverse, most have converged towards remarkable similarity in the size and shape of their capsids1,2. By contrast, generative models for protein design enable the creation of protein architectures that are absent from nature3-5. Here we investigate whether protein assemblies designed by artificial intelligence can be functionalized to construct nucleic acid transport vehicles that are independent of evolutionary trajectories. By combining natural protein domains with synthetic protein assemblies, we create more than 100 bottom-up RNA transfer vehicles with unique sizes and shapes. These vehicles surpass the RNA transfer efficiency of widely used delivery vehicles by several orders of magnitude. In addition, we demonstrate that their tropism can be programmed by incorporation of computationally designed peptide binders and use them to deliver therapeutically relevant cargo RNAs into a wide range of cellular models. We show the in vivo biodistribution of one of these vehicles in a mouse at near-single-cell resolution, confirm its safety, and use it to perform a gene-editing treatment strategy for Duchenne muscular dystrophy in patient-derived cells and a pig. Our work demonstrates how proteins created by generative artificial intelligence can be harnessed for the rational engineering of RNA transport systems with the desired properties by overcoming the limitations of natural protein diversity.

Journal Article↗

PubMed2025

Molecular Evolution and Zoonotic Potential of Muju Virus (Orthohantavirus puumalaense) in Craseomys regulus, Republic of Korea.

Orthohantavirus puumalaense causes hemorrhagic fever with renal syndrome in Europe, with Puumala virus (PUUV) as its primary representative. Muju virus (MUJV), harbored by Craseomys regulus, an Arvicolinae rodent species endemic to the Republic of Korea (ROK), is also a genotype of O. puumalaense. However, their genomic diversity and zoonotic potential remain largely unknown. To investigate their prevalence, 185 voles were collected from 23 regions of the ROK between 2012 and 2023. Serological assays detected anti-PUUV immunoglobulin G antibodies in five samples (3.1%), whereas reverse-transcription polymerase chain reaction confirmed MUJV RNA in identical specimens (2.7%). Amplicon-based nanopore sequencing facilitates near-complete genome recovery, enabling high-resolution comparative analysis. Phylogenetic analysis revealed distinct genetic lineages in Gangwon and Jeollabuk Provinces. Evolutionary rate estimates indicated greater sequence divergence in the S and L segments than in the M segment. A zoonotic risk assessment revealed that most MUJV variants exhibited moderate-to-high spillover potential. The molecular detection of MUJV in Cheorwon, Gangwon Province, expands its known geographic range and provides the first molecular evidence of MUJV circulation in this region. These findings highlight the need for continued surveillance and seroprevalence studies of MUJV to assess its potential for human exposure and public health relevance in the ROK.

Animals↗

PubMed2027

Optimized Amplicon Strategy for Long-Read Sequencing of the Chikungunya Virus Genome.

Chikungunya virus (CHIKV) is a positive-sense RNA alphavirus transmitted to humans primarily by Aedes aegypti and Aedes albopictus mosquitoes. Its global circulation and significant public health impact underscore the need to better understand the molecular mechanisms driving CHIKV pathogenesis and transmission. Although robust molecular biology methods exist for CHIKV genome sequencing, a major limitation for surveillance and research is the inability to determine whether two nucleotide variations co-occur within the same viral genome when they are separated beyond the span of typical short-read designs. Here, we describe an optimized approach for processing CHIKV RNA samples that generates large amplicons suitable for long-read nanopore sequencing. This protocol enables amplification of the complete CHIKV genome in only two or three amplicons and facilitates detection of co-occurring nucleotide variations across 4-7.5 kb within the same molecule, thereby simplifying sequencing workflows and improving resolution in studies of viral evolution.

Chikungunya virus↗

PubMed2025

Deep Sequencing Reveals Dual Evolution of SARS-CoV-2: Insights Into Defective Genomes From Wuhan-Hu-1 Variants to Omicron Subvariants.

SARS-CoV-2 has evolved from early variants dominating the first (B.1.5, B.1.1) and second (B.1.177) pandemic waves, which exhibited a higher frequency of minority mutants with deletions leading to Defective Viral Genomes (DVGs) in the spike region near the S1/S2 cleavage site than the Alpha, Beta, and Delta variants. The emergence of Omicron has significantly altered the dominant variant profile, with Omicron subvariants now representing 100% of circulating viruses. To monitor the evolution and adaptation of Omicron in the human population, a deep-sequencing study was performed in RNA samples of BA.1, BA.1.1, BA.2, BA.5, BQ.1.1, XBB.1.5 and BA.2.86 Omicron subvariants. The findings reveal two occurrences of similar evolutionary patterns within SARS-CoV-2 characterized by a shift from a significant to a very low production of DVGs. This event suggests that DVGs might play a role in the virus's spread and adaptation for persistence in infected humans.

SARS-CoV-2↗

PubMed2026

Ramu stunt virus genome reveals previously unreported segments and nucleocapsid domain duplication in Mechlorovirus.

Ramu stunt virus (RmSV), a member of the genus Mechlorovirus within the family Phenuiviridae, was previously described as a six-segmented RNA virus infecting sugarcane. In this study, we re-examined type material and additional isolates using high-throughput sequencing and RT-PCR validation, revealing that RmSV possesses a nine-segmented genome, making it the largest reported in the Phenuiviridae. This expanded architecture includes duplicated RNA segments (RNA 2a and RNA 2b) encoding nucleocapsid-like proteins and two novel segments (RNA 7 and RNA 8). Comparative analysis showed that RNA 2a and 2b share about 84% amino acid identity, while RNA 5 encodes a third nucleocapsid homolog, indicating unprecedented domain redundancy. Structural modeling confirmed that all three nucleocapsid proteins maintain a conserved fold despite low sequence identity, with electrostatic mapping suggesting differential RNA-binding potential. Additionally, RNA 6 encodes a hypothetical protein structurally similar to the rice stripe virus disease-specific S-protein, implicating a role in symptom development. Transcript abundance analysis revealed RNA 6 as the most highly expressed segment across isolates. These findings revise the genomic composition of RmSV, highlight mechanisms of genome plasticity and adaptive evolution in plant-infecting bunyaviruses, and underscore practical implications for diagnostic assay design, resistance breeding, and biosecurity surveillance.

Genome, Viral↗

PubMed2024

Replication of single viruses across the kingdoms, Fungi, Plantae, and Animalia.

It is extremely rare that a single virus crosses host barriers across multiple kingdoms. Based on phylogenetic and paleovirological analyses, it has previously been hypothesized that single members of the family Partitiviridae could cross multiple kingdoms. Partitiviridae accommodates members characterized by their simple bisegmented double-stranded RNA genome; asymptomatic infections of host organisms; the absence of an extracellular route for entry in nature; and collectively broad host range. Herein, we show the replicability of single fungal partitiviruses in three kingdoms of host organisms: Fungi, Plantae, and Animalia. Betapartitiviruses of the phytopathogenic fungusRosellinia necatrix could replicate in protoplasts of the carrot (Daucus carota), Nicotiana benthamiana and Nicotiana tabacum, in some cases reaching a level detectable by agarose gel electrophoresis. Moreover, betapartitiviruses showed more robust replication than the tested alphapartitiviruses. One of the fungal betapartitiviruses, RnPV18, could persistently and stably infect carrot plants regenerated from virion-transfected protoplasts. Both alpha- and betapartitiviruses, although with different host preference, could replicate in two insect cell lines derived from the fall armyworm Spodoptera frugiperda and the fruit fly Drosophila melanogaster. Our results indicate the replicability of single partitiviruses in members of three kingdoms and provide insights into virus adaptation, host jumping, and evolution.

Animals↗

PubMed2025

Enrichment of G-to-U Substitution in SARS-CoV-2 Functional Regions and Its Compensation via Concurrent Mutations.

We surveyed single nucleotide variant (SNV) patterns from 5 903 647 complete SARS-CoV-2 genomes. Among 10 012 SNVs, APOBEC-mediated C-to-U (C > U) deamination was the most prevalent, followed by G > U and other RNA editing-related substitutions including (A > G, U > C, G > A). However, C > U mutations were less frequent in functional regions, for example, S protein, intrinsic disordered regions, and nonsynonymous mutations, where G > U were over-represented. Notably, G-loss substitutions rarely appeared together. Instead, G-gain mutations tended to more frequently co-occur with others, with a marked preference in the S protein, suggesting a compensatory mechanism for G loss in G > U mutations. The temporal patterns revealed C > U frequency declined until late 2021 then resurged in early 2022. Conversely, G > U steadily decreased, with a pronounced drop in January 2022, coinciding with reduced COVID-19 severity. Vaccinated individuals exhibited a slightly but significantly higher C > U frequency and a notably lower G > U frequency compared to the unvaccinated group. Additionally, cancer patients had higher G > U frequency than general patients during the same period. Interestingly, none of the C > U SNVs were uniquely identified in 2724 environmental samples. These findings suggest novel functional roles of G > U in COVID-19 symptoms, potentially linked to oxidative stress and reactive oxygen species, while C > U remains the dominant substitution, likely driven by host immune-mediated RNA editing.

SARS-CoV-2↗

PubMed2027

Nanopore Sequencing for Chikungunya Virus: Principles and Application.

Nanopore sequencing is transforming viral genomics through real-time, portable, long-read analysis of RNA and DNA. Unlike traditional short-read platforms, it detects nucleotide sequences by measuring ionic current changes as nucleic acids pass through nanoscale pores, enabling direct single-molecule sequencing and base modification detection. Its simplicity, flexibility, and capacity for ultra-long reads make it ideal for resolving complex genomic regions, structural variants, and full viral genomes. These advantages have accelerated its use in pathogen surveillance and outbreak response, especially in resource-limited settings. For chikungunya virus (CHIKV), nanopore sequencing allows rapid, culture-independent recovery of complete genomes from clinical and vector samples, enabling real-time tracking of viral diversity, evolution, and spread. Experiences from Ebola, Zika, and COVID-19 have demonstrated the power of portable sequencing, now applied to CHIKV monitoring. Advances in tools such as Guppy, Dorado, Minimap2, and Medaka enhance read quality, consensus accuracy, and downstream analyses. Despite challenges in basecalling and error correction, robust quality control pipelines ensure reliable results. Ongoing improvements in chemistry, flow cell design, and machine learning will further enhance fidelity and throughput, establishing nanopore sequencing as a cornerstone of CHIKV genomic surveillance and epidemic preparedness.

Chikungunya virus↗

PubMed2026

Evolutionary history of Jamestown Canyon virus reveals complex multi-vector ecology.

Jamestown Canyon virus (JCV) is a historically understudied mosquito-borne virus of increasing concern in North America. We generated 658 whole-genome JCV sequences from northeast United States, including 84% (500/597) of all JCV-positive mosquitoes detected in Connecticut from 1997 to 2022. Then, we applied phylodynamic methods to demonstrate how mosquito phenology structures the maintenance and evolution of JCV. Our phylogenetic analyses estimate that JCV was introduced in the Northeast by at least the early 1700s, and the primary introductions of lineages A and B into Connecticut occurred during the mid-1800s to mid-1900s. Further, we estimate that JCV evolves at a rate of ∼3 × 10-5 substitutions per site per year (s/s/y), making it one of the slowest-evolving known RNA viruses, because the virus spends ∼10 months per year in evolutionary stasis while overwintering in mosquito eggs. To investigate ecological drivers of JCV spread in Connecticut, we paired discrete trait and continuous phylogeographic reconstructions with mosquito surveillance data. We estimate that JCV has a low diffusion rate of ∼30-60 km2/year, which is more similar to slow-moving tick-borne viruses than to other mosquito-borne viruses. We found that univoltine Aedes mosquitoes were likely to maintain the virus across years through overwintering in eggs, accounting for its slow evolution and dispersal, while multivoltine mosquitoes contributed to periodic bursts of spatial diffusion and amplification within seasons. We demonstrate the utility of dense sequencing and phylodynamics to disentangle complex transmission cycles, offering a framework for rapidly advancing our evolutionary and ecological knowledge of understudied viruses.

Animals↗

PubMed2026

Longitudinal whole-genome analysis of bluetongue virus identifies conserved serotype-specific genomes and distinct genomic constellations within a Colorado sheep flock (2021-2023).

Bluetongue virus (BTV) is a segmented double-stranded RNA virus of ruminants transmitted by Culicoides spp. biting midges. Although the genome consists of ten segments, classification into serotypes is primarily based on genome segment 2. However, reassortment among genomic segments is a major driver of BTV evolution and diversity. This study used longitudinal whole-genome sequencing to characterize BTV genomes collected from 2021 to 2023 within a single sheep flock in Colorado, where multiple serotypes co-circulate. Whole-genome sequences were generated from fourteen blood samples representing four serotypes: BTV-6, -11, -13, and -17. Longitudinal sampling identified multiple BTV serotypes within individual sheep across consecutive years. Tanglegram analysis comparing segment phylogenies to the segment 2 tree demonstrated incongruent topologies across all genomic segments, suggestive of reassortment or the circulation of distinct genomic constellations. Nucleotide-level comparisons revealed high sequence homology among same-serotype samples from the same year, while the greatest genetic divergence was observed among BTV-17 genomes collected in different years. Additionally, all BTV-13 genomes contained a previously undescribed nonsynonymous substitution in segment 10 predicted to extend the encoded protein by three amino acids. Together, these findings demonstrate that highly conserved BTV genomes and distinct genomic constellations can be detected at the flock level across multiple years. This longitudinal whole-genome approach reveals the genetic complexity of endemic BTV populations, including novel variants and genomic patterns consistent with reassortment that are lost with conventional serotyped-based approaches, highlighting the need to integrate whole-genome characterization into endemic BTV monitoring programs.

Animals↗