Genotype classification and pathogenicity of infectious bursal disease virus circulating in vaccinated broiler chicken farms.

This study investigated the genotype classification and pathogenicity of infectious bursal disease virus (IBDV) circulating in vaccinated broiler chicken farms in Egypt. A total of 150 samples were collected from 30 vaccinated commercial broiler chicken farms and pooled into 30 working samples. IBDV was tested using reverse transcriptase polymerase chain reaction (RT-PCR) amplification of the hypervariable region of the viral protein 2 (hvVP2) and the VP1 gene 5' extremity. Both RT-PCR fragments were sequenced from six samples, and then the obtained nucleotide sequences were analyzed. The IBDV genotypes were identified using nucleotide sequences. Five sequences of the six strains examined were classified as genotype A3B2 for the highly virulent segments A and B (vv-A/vv-B IBDV). Interestingly, this study identified and classified a novel segment-reassortant strain as the A1B2 genotype. Specifically, it involved the segment reassortment of classical virulent segment A (cv-A) with vv-B producing cv-A/vv-B reassortant IBDV. Subsequently, we compared the pathogenicity of reassortant (cv-A/vv-B) IBDV and vvIBDV strains identified in this study. Both strains developed typical IBD clinical signs, postmortem lesions, histopathology, immunohistochemistry, and lesion scores, which were more severe in vvIBDV than reassortant IBDV. In conclusion, this is the first report of the genotype classification based on both genome segments (hvVP2 and VP1) with pathogenicity of IBDV circulating in vaccinated broiler chicken farms and this pathogenicity is more severe in vvIBDV strain than a novel reassortant IBDV strain.

Molecular characterisation of IBDV isolates in Turkey revealed reassortant strains.

1. Infectious bursal disease (IBD) is an acute, highly contagious viral disease of chickens caused by a virus (IBDV) which has a bi-segmented, double-stranded RNA genome. It has five viral proteins in its structure; the VP1 gene is encoded in segment B and the other four are in segment A.2. In this study, bursae of Fabricius and spleen samples taken from chickens suspected of having clinical or subclinical IBD from a total of 50 chicken flocks located in different geographical regions of Turkey were examined.3. The RT-PCR analysis of the VP2 gene showed that 30 of the 50 samples (60%) tested positive. Eight positive isolates were chosen and RT-PCR was performed to amplify the VP1 gene.4. The study showed that reassortant field strains that cause clinical or subclinical disease are currently circulating in broiler flocks across Turkey.

Genomic diversity and evolutionary dynamics of Influenza A viruses in Colombian swine: implications for one health surveillance and control.

Influenza A viruses (IAV) impose significant respiratory disease burdens in both swine and humans worldwide, with frequent human-to-swine transmission driving viral evolution in pigs and highlighting the risk at the animal-human interface. Therefore, a comprehensive One Health approach (interconnection among human, animal, and environmental health) is needed for IAV prevention, control, and response. Animal influenza genomic surveillance remains limited in many Latin American countries, including Colombia. To address this gap, we genetically characterized 170 swine specimens from Colombia (2011-2017). Whole genome sequencing revealed a predominance of pandemic-like H1N1 lineage, with a minority belonging to H3N2 and H1N2 human seasonal-like lineage and H1N1 early classical swine lineages. Significantly, we have identified reassortant and recombinant viruses (H3N2, H1N1) not previously reported in Colombia. This suggests a broad genotypic viral diversity, likely resulting from reassortment between classical endemic viruses and new introductions established in Colombia's swine population (e.g. the 2009 H1N1 pandemic). Our study highlights the importance of a One Health approach in disease control, particularly in an ecosystem where humans are a main source of IAV to swine populations, and emphasizes the need for continued surveillance and enhanced biosecurity measures. The co-circulation of multiple subtypes in regions with high swine density facilitates viral exchange, underscoring the importance of monitoring viral evolution to inform vaccine selection and public health policies locally and globally.

Driving forces of continuing evolution of rotaviruses.

Rotaviruses are non-enveloped double-stranded RNA virus that causes acute diarrheal diseases in children (< 5 years). More than 90% of the global rotavirus infection in humans was caused by Rotavirus group A. Rotavirus infection has caused more than 200000 deaths annually and predominantly occurs in the low-income countries. Rotavirus evolution is indicated by the strain dynamics or the emergence of the unprecedented strain. The major factors that drive the rotavirus evolution include the genetic shift that is caused by the reassortment mechanism, either in the intra- or the inter-genogroup. However, other factors are also known to have an impact on rotavirus evolution. This review discusses the structure and types, epidemiology, and evolution of rotaviruses. This article also reviews other supplemental factors of rotavirus evolution, such as genetic reassortment, mutation rate, glycan specificity, vaccine introduction, the host immune responses, and antiviral drugs.

Molecular evolutionary analysis of novel NSP4 mono-reassortant G1P[8]-E2 rotavirus strains that caused a discontinuous epidemic in Japan in 2015 and 2018.

In the 2010s, several unusual rotavirus strains emerged, causing epidemics worldwide. This study reports a comprehensive molecular epidemiological study of rotaviruses in Japan based on full-genome analysis. From 2014 to 2019, a total of 489 rotavirus-positive stool specimens were identified, and the associated viral genomes were analyzed by next-generation sequencing. The genotype constellations of those strains were classified into nine patterns (G1P[8] (Wa), G1P[8]-E2, G1P[8] (DS-1), G2P[4] (DS-1), G3P[8] (Wa), G3P[8] (DS-1), G8P[8] (DS-1), G9P[8] (Wa), and G9P[8]-E2). The major prevalent genotype differed by year, comprising G8P[8] (DS-1) (37% of that year's isolates) in 2014, G1P[8] (DS-1) (65%) in 2015, G9P[8] (Wa) (72%) in 2016, G3P[8] (DS-1) (66%) in 2017, G1P[8]-E2 (53%) in 2018, and G9P[8] (Wa) (26%) in 2019. The G1P[8]-E2 strains (G1-P[8]-I1-R1-C1-M1-A1-N1-T1-E2-H1) isolated from a total of 42 specimens in discontinuous years (2015 and 2018), which were the newly-emerged NSP4 mono-reassortant strains. Based on the results of the Bayesian evolutionary analyses, G1P[8]-E2 and G9P[8]-E2 were hypothesized to have been generated from distinct independent inter-genogroup reassortment events. The G1 strains detected in this study were classified into multiple clusters, depending on the year of detection. A comparison of the predicted amino acid sequences of the VP7 epitopes revealed that the G1 strains detected in different years encoded VP7 epitopes harboring distinct mutations. These mutations may be responsible for immune escape and annual changes in the prevalent strains.

Genomic characterization of highly pathogenic avian influenza A H5N1 virus newly emerged in dairy cattle.

In March 2024, the emergence of highly pathogenic avian influenza (HPAI) A (H5N1) infections in dairy cattle was detected in the United Sates for the first time. We genetically characterize HPAI viruses from dairy cattle showing an abrupt drop in milk production, as well as from two cats, six wild birds, and one skunk. They share nearly identical genome sequences, forming a new genotype B3.13 within the 2.3.4.4b clade. B3.13 viruses underwent two reassortment events since 2023 and exhibit critical mutations in HA, M1, and NS genes but lack critical mutations in PB2 and PB1 genes, which enhance virulence or adaptation to mammals. The PB2 E627 K mutation in a human case associated with cattle underscores the potential for rapid evolution post infection, highlighting the need for continued surveillance to monitor public health threats.

Evolution of DS-1-like G8P[8] rotavirus A strains from Vietnamese children with acute gastroenteritis (2014-21): Adaptation and loss of animal rotavirus-derived genes during human-to-human spread.

Animal rotaviruses A (RVAs) are considered the source of emerging, novel RVA strains that have the potential to cause global spread in humans. A case in point was the emergence of G8 bovine RVA consisting of the P[8] VP4 gene and the DS-1-like backbone genes that appeared to have jumped into humans recently. However, it was not well documented what evolutionary changes occurred on the animal RVA-derived genes during circulation in humans. Rotavirus surveillance in Vietnam found that DS-1-like G8P[8] strains emerged in 2014, circulated in two prevalent waves, and disappeared in 2021. This surveillance provided us with a unique opportunity to investigate the whole process of evolutionary changes, which occurred in an animal RVA that had jumped the host species barrier. Of the 843 G8P[8] samples collected from children with acute diarrhoea in Vietnam between 2014 and 2021, fifty-eight strains were selected based on their distinctive electropherotypes of the genomic RNA identified using polyacrylamide gel electrophoresis. Whole-genome sequence analysis of those fifty-eight strains showed that the strains dominant during the first wave of prevalence (2014-17) carried animal RVA-derived VP1, NSP2, and NSP4 genes. However, the strains from the second wave of prevalence (2018-21) lost these genes, which were replaced with cognate human RVA-derived genes, thus creating strain with G8P[8] on a fully DS-1-like human RVA gene backbone. The G8 VP7 and P[8] VP4 genes underwent some point mutations but the phylogenetic lineages to which they belonged remained unchanged. We, therefore, propose a hypothesis regarding the tendency for the animal RVA-derived genes to be expelled from the backbone genes of the progeny strains after crossing the host species barrier. This study underlines the importance of long-term surveillance of circulating wild-type strains in order to better understand the adaptation process and the fate of newly emerging, animal-derived RVA among the human population. Further studies are warranted to disclose the molecular mechanisms by which spillover animal RVAs become readily transmissible among humans, and the roles played by the expulsion of animal-derived genes and herd immunity formed in the local population.

Genetic characterization and evidence for multiple reassortments of rotavirus A G3P[3] in dogs and cats in Thailand.

Rotavirus A (RVA) causes gastroenteritis in humans and animals. The zoonotic potential of RVA has been reported and raises major concerns, especially in animal-human interface settings. The study aimed to characterize and investigate the genetic diversity among RVAs in dogs and cats in Thailand. We collected 572 rectal swab samples from dogs and cats in Bangkok animal hospitals from January 2020 to June 2021. The one-step RT-PCR assay detected RVAs in 1.92% (11/572) of the samples, with 2.75% (8/290) in dogs and 1.06% (3/282) in cats. Two canine RVA and one feline RVA were subjected to whole genome sequencing. Our results showed that all three viruses were identified as RVA genotype G3P[3]. The genetic constellation of RVAs is unique for different species. For canine RVAs is G3-P [3]-I3-R3-C3-M3-A9-N2-T3-E3-H6, while Feline RVA is G3-P [3]-I8-R3-C3-M3-A9-N3-T3-E3-H6. Notably, both canine and feline RVAs contained the AU-1 genetic constellation with multiple reassortments. The results of phylogenetic, genetic, and bootscan analyses showed that canine RVAs may have reassorted from dog, human, and cat RVAs. While feline RVA was closely related to RVAs in humans, bats, and simians. This study provided genetic characteristics and diversity of RVAs in dogs and cats and suggested possible multiple reassortments, suggesting the zoonotic potential of the viruses. Thus, public health awareness should be raised regarding the zoonotic potential of RVAs in dogs and cats. Further studies on RVAs on a larger scale in dogs and cats in Thailand are needed.

IAVCP (Influenza A Virus Consensus and Phylogeny): Automatic Identification of the Genomic Sequence of the Influenza A Virus from High-Throughput Sequencing Data.

Recently, high-throughput sequencing of influenza A viruses has become a routine test. It should be noted that the extremely high diversity of the influenza A virus complicates the task of determining the sequences of all eight genome segments. For a fast and accurate analysis, it is necessary to select the most suitable reference for each segment. At the same time, there is no standardized method in the field of decoding sequencing results that allows the user to update the sequence databases to which the reads obtained by virus sequencing are compared. The IAVCP (influenza A virus consensus and phylogeny) was developed with the goal of automatically analyzing high-throughput sequencing data of influenza A viruses. Its goals include the extraction of a consensus genome directly from paired raw reads. In addition, the pipeline enables the identification of potential reassortment events in the evolutionary history of the virus of interest by analyzing the topological structure of phylogenetic trees that are automatically reconstructed.

Rift Valley Fever Phlebovirus Reassortment Study in Sheep.

Rift Valley fever (RVF) in ungulates and humans is caused by a mosquito-borne RVF phlebovirus (RVFV). Live attenuated vaccines are used in livestock (sheep and cattle) to control RVF in endemic regions during outbreaks. The ability of two or more different RVFV strains to reassort when co-infecting a host cell is a significant veterinary and public health concern due to the potential emergence of newly reassorted viruses, since reassortment of RVFVs has been documented in nature and in experimental infection studies. Due to the very limited information regarding the frequency and dynamics of RVFV reassortment, we evaluated the efficiency of RVFV reassortment in sheep, a natural host for this zoonotic pathogen. Co-infection experiments were performed, first in vitro in sheep-derived cells, and subsequently in vivo in sheep. Two RVFV co-infection groups were evaluated: group I consisted of co-infection with two wild-type (WT) RVFV strains, Kenya 128B-15 (Ken06) and Saudi Arabia SA01-1322 (SA01), while group II consisted of co-infection with the live attenuated virus (LAV) vaccine strain MP-12 and a WT strain, Ken06. In the in vitro experiments, the virus supernatants were collected 24 h post-infection. In the in vivo experiments, clinical signs were monitored, and blood and tissues were collected at various time points up to nine days post-challenge for analyses. Cell culture supernatants and samples from sheep were processed, and plaque-isolated viruses were genotyped to determine reassortment frequency. Our results show that RVFV reassortment is more efficient in co-infected sheep-derived cells compared to co-infected sheep. In vitro, the reassortment frequencies reached 37.9% for the group I co-infected cells and 25.4% for the group II co-infected cells. In contrast, we detected just 1.7% reassortant viruses from group I sheep co-infected with the two WT strains, while no reassortants were detected from group II sheep co-infected with the WT and LAV strains. The results indicate that RVFV reassortment occurs at a lower frequency in vivo in sheep when compared to in vitro conditions in sheep-derived cells. Further studies are needed to better understand the implications of RVFV reassortment in relation to virulence and transmission dynamics in the host and the vector. The knowledge learned from these studies on reassortment is important for understanding the dynamics of RVFV evolution.

Molecular Markers and Mechanisms of Influenza A Virus Cross-Species Transmission and New Host Adaptation.

Influenza A viruses continue to be a serious health risk to people and result in a large-scale socio-economic loss. Avian influenza viruses typically do not replicate efficiently in mammals, but through the accumulation of mutations or genetic reassortment, they can overcome interspecies barriers, adapt to new hosts, and spread among them. Zoonotic influenza A viruses sporadically infect humans and exhibit limited human-to-human transmission. However, further adaptation of these viruses to humans may result in airborne transmissible viruses with pandemic potential. Therefore, we are beginning to understand genetic changes and mechanisms that may influence interspecific adaptation, cross-species transmission, and the pandemic potential of influenza A viruses. We also discuss the genetic and phenotypic traits associated with the airborne transmission of influenza A viruses in order to provide theoretical guidance for the surveillance of new strains with pandemic potential and the prevention of pandemics.

Phylogenetic Analyses of Rotavirus A, B and C Detected on a Porcine Farm in South Africa.

Rotaviruses (RVs) are known to infect various avian and mammalian hosts, including swine. The most common RVs associated with infection in pigs are A, B, C and H (RVA-C; RVH). In this study we analysed rotavirus strains circulating on a porcine farm in the Western Cape province of South Africa over a two-year period. Whole genomes were determined by sequencing using Illumina MiSeq without prior genome amplification. Fifteen RVA genomes, one RVB genome and a partial RVC genome were identified. Phylogenetic analyses of the RVA data suggested circulation of one dominant strain (G5-P[6]/P[13]/P[23]-I5-R1-C1-M1-A8-N1-T7-E1-H1), typical of South African porcine strains, although not closely related to previously detected South African porcine strains. Reassortment with three VP4-encoding P genotypes was detected. The study also reports the first complete RVB genome (G14-P[5]-I13-R4-C4-M4-A10-T4-E4-H7) from Africa. The partial RVC (G6-P[5]-IX-R1-C1-MX-A9-N6-T6-EX-H7) strain also grouped with porcine strains. The study shows the continued circulation of an RVA strain, with a high reassortment rate of the VP4-encoding segment, on the porcine farm. Furthermore, incidents of RVB and RVC on this farm emphasize the complex epidemiology of rotavirus in pigs.

Genetic Reassortment in a Child Coinfected with Two Influenza B Viruses, B/Yamagata Lineage and B/Victoria-Lineage Strains.

We identified a child coinfected with influenza B viruses of B/Yamagata and B/Victoria lineages, in whom we analyzed the occurrence of genetic reassortment. Plaque purification was performed using a throat swab specimen from a 9-year-old child, resulting in 34 well-isolated plaques. The genomic composition of eight gene segments (HA, NA, PB1, PB2, PA, NP, M, and NS genes) for each plaque was determined at the lineage level. Of the 34 plaques, 21 (61.8%) had B/Phuket/3073/2013 (B/Yamagata)-like sequences in all gene segments, while the other 13 (38.2%) were reassortants with B/Texas/02/2013 (B/Victoria)-like sequences in 1-5 of the 8 segments. The PB1 segment had the most B/Victoria lineage genes (23.5%; 8 of 34 plaques), while PB2 and PA had the least (2.9%; 1 of 34 plaques). Reassortants with B/Victoria lineage genes in 2-5 segments showed the same level of growth as viruses with B/Yamagata lineage genes in all segments. However, reassortants with B/Victoria lineage genes only in the NA, PB1, NP, or NS segments exhibited reduced or undetectable growth. We demonstrated that various gene reassortments occurred in a child. These results suggest that simultaneous outbreaks of two influenza B virus lineages increase genetic diversity and could promote the emergence of new epidemic strains.

Avian Influenza outbreaks: Human infection risks for beach users - One health concern and environmental surveillance implications.

Despite its popularity for water activities, such as swimming, surfing, fishing, and rafting, inland and coastal bathing areas occasionally experience outbreaks of highly pathogenic avian influenza virus (HPAI), including A(H5N1) clade 2.3.4.4b. Asymptomatic infections and symptomatic outbreaks often impact many aquatic birds, which increase chances of spill-over events to mammals and pose concerns for public health. This review examined the existing literature to assess avian influenza virus (AIV) transmission risks to beachgoers and the general population. A comprehensive understanding of factors governing such crossing of the AIV host range is currently lacking. There is limited knowledge on key factors affecting risk, such as species-specific interactions with host cells (including binding, entry, and replication via viral proteins hemagglutinin, neuraminidase, nucleoprotein, and polymerase basic protein 2), overcoming host restrictions, and innate immune response. AIV efficiently transmits between birds and to some extent between marine scavenger mammals in aquatic environments via consumption of infected birds. However, the current literature lacks evidence of zoonotic AIV transmission via contact with the aquatic environment or consumption of contaminated water. The zoonotic transmission risk of the circulating A(H5N1) clade 2.3.4.4b virus to the general population and beachgoers is currently low. Nevertheless, it is recommended to avoid direct contact with sick or dead birds and to refrain from bathing in locations where mass bird mortalities are reported. Increasing reports of AIVs spilling over to non-human mammals have raised valid concerns about possible virus mutations that lead to crossing the species barrier and subsequent risk of human infections and outbreaks.

Development and application of a whole genome amplicon sequencing method for infectious salmon anemia virus (ISAV).

Infectious salmon anemia (ISA) is an infectious disease primarily affecting farmed Atlantic salmon, Salmo salar, which is caused by the ISA virus (ISAV). ISAV belongs to the Orthomyxoviridae family. The disease is a serious condition resulting in reduced fish welfare and high mortality. In this study, we designed an amplicon-based sequencing protocol for whole genome sequencing of ISAV. The method consists of 80 ISAV-specific primers that cover 92% of the virus genome and was designed to be used on an Illumina MiSeq platform. The sequencing accuracy was investigated by comparing sequences with previously published Sanger sequences. The sequences obtained were nearly identical to those obtained by Sanger sequencing, thus demonstrating that sequences produced by this amplicon sequencing protocol had an acceptable accuracy. The amplicon-based sequencing method was used to obtain the whole genome sequence of 12 different ISAV isolates from a small local epidemic in the northern part of Norway. Analysis of the whole genome sequences revealed that segment reassortment took place between some of the isolates and could identify which segments that had been reassorted.

Isolation, characterization and whole-genome analysis of G9 group a rotaviruses in China: Evidence for possible Porcine-Human interspecies transmission.

Group A rotaviruses (RVAs) are major causes of severe gastroenteritis in infants and young animals. To enhance our understanding of the relationship between human and animals RVAs, complete genome data are necessary. We screened 92 intestinal and stool samples from diarrheic piglets by RT‒PCR targeting the VP6 gene, revealing a prevalence of 10.9%. RVA was confirmed in two out of 5 calf samples. We successfully isolated two porcine samples using MA104 cell line. The full-length genetic constellation of the two isolates were determined to be G9-P[23]-I5-R1-C1-M1-A8-N1-T7-E1-H1, with close similarity to human Wa-like and porcine strains. Sequence analysis revealed the majority of genes were closely related to porcine and human RVAs. Phylogenetic analysis revealed that these isolates might have their ancestral origin from pigs, although some of their gene segments were related to human strains. This study reveals evidence of reassortment and possible interspecies transmission between pigs and humans in China.

Pathogenicity and transmission of novel highly pathogenic H7N2 variants originating from H7N9 avian influenza viruses in chickens.

The H7 subtype avian influenza viruses are circulating widely worldwide, causing significant economic losses to the poultry industry and posing a serious threat to human health. In 2019, H7N2 and H7N9 co-circulated in Chinese poultry, yet the risk of H7N2 remained unclear. We isolated and sequenced four H7N2 viruses from chickens, revealing them as novel reassortants with H7N9-derived HA, M, NS genes and H9N2-derived PB2, PB1, PA,NP, NA genes. To further explore the key segment of pathogenicity, H7N2-H7N9NA and H7N2-H9N2HA single-substitution were constructed. Pathogenicity study showed H7N2 isolates to be highly pathogenic in chickens, with H7N2-H7N9NA slightly weaker than H7N2-Wild type. Transcriptomic analysis suggested that H7N9-derived HA genes primarily drove the high pathogenicity of H7N2 isolates, eliciting a strong inflammatory response. These findings underscored the increased threat posed by reassorted H7N2 viruses to chickens, emphasizing the necessity of long-term monitoring of H7 subtype avian influenza viruses.

Novel Genotype of HA Clade 2.3.4.4b H5N8 Subtype High Pathogenicity Avian Influenza Virus Emerged at a Wintering Site of Migratory Birds in Japan, 2021/22 Winter.

Surveillance of avian influenza virus (AIV) was conducted in the 2021-2022 winter season at a wintering site of migratory Anatidae in Japan. An H5N8 subtype high pathogenicity AIV (HPAIV) with a unique gene constellation and four low pathogenicity AIVs (LPAIVs) were isolated from environmental samples. The genetic origin of the HPAIV (NK1201) was determined with whole-genome sequencing and phylogenetic analyses. Six of NK1201's eight genes were closely related to HA clade 2.3.4.4b H5N8 subtype HPAIVs, belonging to the G2a group, which was responsible for outbreaks in poultry farms in November 2021 in Japan. However, the remaining two genes, PB1 and NP, most closely matched those of the LPAIVs H7N7 and H1N8, which were isolated at the same place in the same 2021-2022 winter. No virus of the NK1201 genotype had been detected prior to the 2021-2022 winter, indicating that it emerged via genetic reassortment among HPAIV and LPAIVs, which were prevalent at the same wintering site. In addition, experimental infection in chickens indicated that NK1201 had slightly different infectivity compared to the reported infectivity of the representative G2a group H5N8 HPAIV, suggesting that the PB1 and NP genes derived from LPAIVs might have affected the pathogenicity of the virus in chickens. Our results directly demonstrate the emergence of a novel genotype of H5N8 HPAIV through gene reassortment at a wintering site. Analyses of AIVs at wintering sites can help to identify the emergence of novel HPAIVs, which pose risks to poultry, livestock, and humans.

The Number and Pattern of Viral Genomic Reassortments are not Necessarily Identifiable from Segment Trees.

Reassortment is an evolutionary process common in viruses with segmented genomes. These viruses can swap whole genomic segments during cellular co-infection, giving rise to novel progeny formed from the mixture of parental segments. Since large-scale genome rearrangements have the potential to generate new phenotypes, reassortment is important to both evolutionary biology and public health research. However, statistical inference of the pattern of reassortment events from phylogenetic data is exceptionally difficult, potentially involving inference of general graphs in which individual segment trees are embedded. In this paper, we argue that, in general, the number and pattern of reassortment events are not identifiable from segment trees alone, even with theoretically ideal data. We call this fact the fundamental problem of reassortment, which we illustrate using the concept of the "first-infection tree," a potentially counterfactual genealogy that would have been observed in the segment trees had no reassortment occurred. Further, we illustrate four additional problems that can arise logically in the inference of reassortment events and show, using simulated data, that these problems are not rare and can potentially distort our observation of reassortment even in small data sets. Finally, we discuss how existing methods can be augmented or adapted to account for not only the fundamental problem of reassortment, but also the four additional situations that can complicate the inference of reassortment.

Evolution of H6N6 viruses in China between 2014 and 2019 involves multiple reassortment events.

H6N6 avian influenza viruses (AIVs) have been widely detected in wild birds, poultry, and even mammals. Recently, H6N6 viruses were reported to be involved in the generation of H5 and H7 subtype viruses. To investigate the emergence, evolutionary pattern, and potential for an epidemic of H6N6 viruses, the complete genomes of 198 H6N6 viruses were analyzed, including 168 H6N6 viruses deposited in the NCBI and GISAID databases from inception to January 2019 and 30 isolates collected from China between November 2014 and January 2019. Using phylogenetic analysis, the 198 strains of H6N6 viruses were identified as 98 genotypes. Molecular clock analysis indicated that the evolution of H6N6 viruses in China was constant and not interrupted by selective pressure. Notably, the laboratory isolates reassorted with six subtype viruses: H6N2, H5N6, H7N9, H5N2, H4N2, and H6N8, resulting in nine novel H6N6 reassortment events. These results suggested that H6N6 viruses can act as an intermediary in the evolution of H5N6, H6N6, and H7N9 viruses. Animal experiments demonstrated that the 10 representative H6N6 viruses showed low pathogenicity in chickens and were capable of infecting mice without prior adaptation. Our findings suggest that H6N6 viruses play an important role in the evolution of AIVs, and it is necessary to continuously monitor and evaluate the potential epidemic of the H6N6 subtype viruses.