Identification of a potential interspecies reassortant rotavirus G and avastrovirus 2 co-infection from black-headed gull (Chroicocephalus ridibundus) in Hungary

Abstract

The black-headed gull is the most common nesting gull species in Hungary. Based on the lifestyle and feeding habits of the black-headed gull, which is highly adapted to the human environment, they can be carriers and spreaders of potential human and other animal pathogens. Between 2014 and 2018 within the framework of the “Life Bird Ringing program” a total of 7 faecal samples were collected from gulls and one sample (MR04) was randomly selected for viral metagenomics and mass sequencing. 95.4% and 4% of the reads were classified into family Seadornaviridae and Astroviridae, respectively, and then were verified by RT-PCR method. In this study, the complete genome of a potential interspecies reassortant rotavirus (RV) strain gull/MR04_RV/HUN/2014 (PP239049-PP239059) and the partial ORF1ab, complete ORF2 of a novel avian nephritis virus strain gull/MR04_AAstV/HUN/2014 (PP239060) was discussed. The strain gull/MR04_RV/HUN/2014 was closely related to rotavirus G (RVG) viruses based on the proteins VP1–VP3, VP6, NSP2, NSP3, and NSP5, but it was more related to the human rotavirus B (RVB) strain Bang373 based on the NSP1, NSP4 and VP7, VP4 proteins, which is assumed to be the result of reassortment between different RVG-RVB rotavirus species. The strain gull/MR04_AAstV/HUN/2014 belonged to the genus Avastrovirus species avastrovirus 2 (AAstV-2) and is related to members of group 6 of avian nephritis viruses (ANVs), but based on the genetic distances it may be the first representative of a separate group. Additional gull samples were found to be negative by RT-PCR. Gulls, which are well adapted to the human environment, could potentially spread enterically transmitted viral pathogens like interspecies reassortant rotaviruses (RVG/RVB), but further molecular surveillance is needed to explore more deeply the viral communities of gulls or other related species adapted to human environments.

Introduction

Gulls belong to the family Laridae of birds (Aves), which are a group well adapted to the human environment found worldwide. The species of common black-headed gull (Chroicocephalus ridibundus) is the most common nesting gull species in Hungary (https://mme.hu/magyarorszagmadarai/madaradatbazis-larrid) and widespread species in Europe, too (http://datazone.birdlife.org/species/factsheet/black-headed-gull-larus-ridibundus) [1]. In terms of its lifestyle, it is a short-term migratory bird species, its domestic breeding population is 4,000–6,400 pairs (2015–2017) in Hungary, and its wintering area is the central and western part of the Mediterranean, and Western Europe. Despite its frequency, it is a protected species (least concern/LC) in Hungary. Regarding its diet, the species exhibits a highly adaptable feeding behaviour, provisioning its young with the most accessible food sources. In agricultural landscapes, it forages for insects, amphibians, invertebrates, and small mammals unearthed by fresh ploughing. Urban areas and their surroundings, however, offer a more abundant and easily accessible variety of resources, including waste from garbage dumps, carrion sites, and fish markets. This species may act as a carrier and vector for pathogens that pose risks to humans and other animals. To better understand these dynamics, next-generation sequencing (NGS) methods could be employed to identify and characterize potential pathogens, providing valuable insights into their transmission and impact [2].

Rotaviruses (RVs) (family Seadornaviridae) and astroviruses (family Astroviridae) are one of the most important and most frequently occurring viral pathogens in faeces which have the important pathogenic potential for humans and livestock (porcine/bovine RVs; porcine mamastro-, and chicken/turkey avastroviruses) [3,4]. The RV genome consists of 11 dsRNA segments that encode eight nonstructural proteins (RdRp, Cap and NSP1–NSP5/6), and six structural proteins (VP1–VP4, VP6, and VP7). Based on the analysis of conservative genome segments (e.g., Seg1) comparisons [5] and phylogenetic analysis of the intermediate capsid protein (VP6) (Seg6) [6], RVs can currently be classified into 9 groups A to J recognized by the International Committee on Taxonomy of Viruses (ICTV), where rotavirus A (RVA), C (RVC), D (RVD) and F (RVF) forming RVA-like clade (clade 1) and rotavirus B (RVB), G (RVG), H (RVH), I (RVI), J (RVJ) and forming RVB-like clade (clade 2) [3,5,7,8]. Additional strains are awaiting taxonomic classification representing novel RVL and RVK groups [9]. Species in RVA, RVD, RVF and RVG have been identified in bird species [8–10]. Several evolutionary mechanisms (mutation, recombination, and reassortment) have already been documented between strains of the same RV species [11–17], and recombination events between different species (interspecies) were also reported [18,19], but no information is available on reassortment between segments of different RV species.

Astroviruses have a positive sense, single-stranded RNA genome, having a genome-linked viral protein (VPg) on the 5’ terminal side and a polyadenylated tail (poly(A)) on the 3’ end. Astroviruses were initially considered as an etiological agent of enteric disease [20]; however, astroviruses have already been identified from infections causing central nervous system disease [21], hepatitis [22] and interstitial nephritis with growth retardation [23] in many mammalian and bird species. It seems that astroviruses have a broad host range and genetic diversity [24]. Currently, the members of the genus Mammastrovirus are infecting mammals and the members of the genus Avastrovirus could infect birds were distinguished [4,24]. Currently, species classification is based on comparative genetic analyses of the ORF2 capsid protein transcribed from sgRNA [25].

In this study, a faecal sample of a common black-headed gull was investigated by the NGS method and data was analyzed by an “in-house” data processing pipeline. We report the detection and complete genome characterization of a novel RV and an avastrovirus (AAstV) from a common black-headed gull in Hungary.

Materials and methods

Sample collection

A total of seven faecal samples were collected as part of the “Life Bird Ringing” program between 2014 and 2018 in Rétszilas, Hungary [26]. These included five samples from common black-headed gulls (Chroicocephalus ridibundus) (HA15097, MR03, MR04, MR05, and SH05705) and two samples from yellow-legged gulls (Larus michahellis) (RE00320 and RE00999). In Hungary, the annual number of ringed black-headed gulls typically ranges from a few hundred individuals, with the highest recorded number being 2,246. The breeding population sampled in Rétszilas consisted of 80–100 breeding pairs. In contrast, the yellow-legged gull has a significantly smaller population size, with fewer than 10 breeding pairs, a trend that is observed nationwide and is not limited to Rétszilas. The licensed ringer and their assistants entered the nesting colony for 3–4 hours on a suitable day in late May or early June, depending on the age of the chicks, to apply leg bands and collect faecal samples during defecation events. Approximately 2–5 grams of faecal samples were collected as a one-off from gull chicks on the nesting island during the ringing process. The samples were transferred into Eppendorf tubes and stored at −20°C for subsequent analysis.

Ethics statement

The black-headed gull (Chroicocephalus ridibundus) is listed as “Least Concern” by both eBird (https://ebird.org/), the International Union for Conservation of Nature’s Society of the Red List of Threatened Species (https://www.iucnredlist.org/) and the BirdLife International (https://www.birdlife.org/). The Hungarian Ornithological and Nature Association (MME) (https://mme.hu/) works internationally as a member of the BirdLife International Partnership and undertakes practical work (survey, observation and ringing) with valid permission to conserve Hungary’s biodiversity. During the “Life bird ringing” program, the faecal samples were collected non-invasively by qualified ornithologists with valid permission (14/3859-9/2012, PE-KTF/97-13/2017), preserving the physical fitness of the birds. The ornithologist’s currently valid certificate (TMF-1180/22/2003) was issued by the National Inspectorate for Environment, Nature and Water.

Sample preparation for viral metagenomics and next-generation sequencing

The faecal samples were homogenized mechanically with 0.1 M phosphate-buffered saline (35–40V/v%) and the total RNA was extracted with TRI Reagent (MRC, USA) according to the manufacturer’s protocol. 200 µL supernatant of faecal suspension of sample MR04 selected for viral metagenomics and next-generation sequencing (VM-NGS) analyses as described previously [27,28]. Briefly, the sample was first filtered through a 0.45 µm filter (Millipore), treated with Nuclease-mix (Turbo DNase, Ambion; Baseline-Zero. Epicentre; Benzonase, Novagen and RNase Promega). After total nucleic acid isolation (Quick-RNA Viral Kit, Zymo Research, Irvine, USA) and reverse transcription, DNA was amplified by random PCR. cDNA library was constructed by NEBNext Ultra II FS DNA Library Prep Kit for Illumina (NEB, Ipswitch, MA, USA). Illumina sequencing was performed on NovaSeq 6000 instrument (Illumina, San Diego, CA, USA) with a 2x151 run configuration [27,28]. The generated NGS sequence data was filtered, analyzed and classified using an in-house assembled bioinformatics pipeline which includes the Kaiju and Diamond algorithms [28]. Virus-related reads were verified using Megan and NCBI BlastX methods. Detailed bash and python scripts and settings used in this study for data filtering, creation and management of the Kaiju database, and Diamond-based validation are available upon request.

Sample preparation for RT-PCR “screening” and Sanger sequencing

RV and astrovirus-related reads were aligned (MUSCLE aligner, https://www.ebi.ac.uk/jdispatcher/msa/muscle) to the most closely related sequences found in the GenBank (by BLAST searches). The alignments were visualized using GeneDoc and Geneious Prime software. Based on the aligned reads, several read-specific primers were designed to characterize the full-length viral sequences using conventional RT-PCR reactions, primer-walking [29], terminal deoxyribonucleotide transferase-based [30] and adapter-ligation-based 5′/3′ RACE [31,32] techniques. Based on multiple sequence alignments of reference and study strains, generic primer-pairs for VP6 of RV (MR04-GVP6-SCR-R: 5′-TGACACCAWACYTTRTC-3′, MR04-GVP6-SCR-F: 5′-ATCAATGATTACAATGC-3′, W: A/T, Y: C/T, R: A/G) and sequence-specific primer-pair for AAstV (MR04-AvAst2-SCR-R: 5′-GTTCATTACCTCAGACACATG-3′, MR04-AvAst2-F: 5′-AAGAGCTGGATGGTGGCTT-3′) were designed for epidemiological investigations. PCR products were directly sequenced with BigDye Terminator v1.1 (Thermo-Fisher) and run on an automated Sanger sequencer (Genetic Analyzer 3500, Applied Biosystems, Hitachi, Japan). Oligonucleotides were ordered from Integrated DNA Technologies (IDT DNA, https://eu.idtdna.com).

Comparative and phylogenetic analysis of sequence data

Differences between nucleotide sequences were determined using Needleman-Wunsch global alignment (https://blast.ncbi.nlm.nih.gov/Blast.cgi), where identities (%), gaps (%) and NW Score were used for analysis and the calculated values were visualized by R/RStudio. Genetic clustering measurements were estimated using with R/RStudio with the help of “seqinr”, “msa”, “cluster”, “factoextra”, “randomcoloR”, “ggrepel” and “ape” libraries. Scripts for graphing NW data and cluster analysis of sequences are available on request. The nucleotide sequences of segments 6 and 11 were aligned by Multiple alignment program based on fast Fourier transform (MAFFT) v7.525 (2024/Mar/13; https://www.ebi.ac.uk/jdispatcher/msa/mafft) [33] and then analysed using IQ-TREE software using the following parameters: -m MFP, -alrt 1000, -T AUTO, -bb 100 [34–36]. According to the Bayesian Information Criterion (BIC) score the best-fit model for the VP6 dataset was TIM2 + F + R3 and for the NSP5 dataset was TIM2 + F + I + G4 models, respectively. Both the SH-aLRT test and the ultrafast bootstrap (UFBoot) values were set up to 1000 is each tree creation, and then the maximum likelihood consensus tree was used. The similarity analysis between the study strain gull/MR04-RV/HUN/2014 (RVG) and the ruddy turnstone rotavirus isolate MW05 (RVG) and human rotavirus B strain Bang373 (RVB) sequences included in the analysis was made using SimPlot 3.5.1 software [37]. For the analysis, the entire genomes of the study strain gull/MR04-RV/HUN/2014 genome segments (PP239049–PP239059) rotavirus, the closest relative ruddy turnstone rotavirus isolate MW05 (MH453863–MH453873) and human rotavirus B strain Bang373 (EU490415, EU490418, AY238384, AY238385, AY238388–AY238394) were concatenated and aligned using MAFFT method. The aligned nucleotide sequences were analysed with the SimPlot program based on the following setup: Window: 500 bp, Step: 20 bp, GapStrip: On, F84 (“Maximum Likelihood”), T/t: 2,0). Phylogenetic analysis of additional segments of RV, the deduced amino acid sequences were aligned using the MEGA11/Muscle method with default settings and the aligned amino acid sequences were tested with MEGA11/Find Best DNA/protein model search [38]. The statistical method determined by the lowest BIC scores was chosen as the basis for the phylogenetic analysis in each tree (S3 Fig). The phylogenetic analysis of astroviruses capsid sequences was aligned by the ClustalW method in MEGA11 and MAFFT alignment method, respectively (S4 Fig) and was tested with MEGA11/Find Best DNA/protein model search [38]. Maximum Likelihood statistical method, Lg with freqs (+F) model was chosen to estimate phylogenetic analysis. Rates among sites were calculated with Gamma Distributed with invariant sites (G + I) option and bootstrap (1000 replicates) was also used. After the pre-defined groups were formed, the mean group distances within groups were calculated with MEGA11 software [38]. The cluster analysis was performed using a script written in R (available for publication upon request) on avian nephritis viruses (ANV) sequences belonging to group 6.

Placement of sequences in a public database

The complete nucleotide sequences of all genomic segments of RV strain gull/MR04_RV/HUN/2014 (PP239049-PP239059) and the partial ORF1ab-complete ORF2 sequence of AAstV strain gull/MR04_AAstV/HUN/2014 (PP239060) were submitted to the GenBank database in National Library of Medicine (NIH) also shared with DNA DataBank of Japan (DDBJ) and the European Nucleotide Archive (ENA).

Results

After the bioinformatic analysis, a total of 1,305,096 virus-related reads were obtained which can be classified into virus families infecting protozoa (N = 504), bacteria (N = 24), plants (N = 18,108), fungi (N = 194), eucaryotic algae (N = 36), invertebrates (N = 220) and vertebrates (N = 1,286,010) (S1 Fig). The viral families infecting vertebrates (98.5% of the virus-related reads) were selected for further analysis where the reads can be potentially classified into family Sedoreoviridae (95.4%), Astroviridae (4%), less than 1% was the sum of families Coronaviridae, Picobirnaviridae, Caliciviridae, Retroviridae and Picornaviridae (S1 Fig.).

Characterization of RV strain gull/MR04_RV/HUN/2014

The determined complete genome of RV strain gull/MR04_RV/HUN/2014 (PP239049-PP239059) is 18,352 nucleotide (nt) long. The lengths and detailed descriptions of the segments were found in the S1 Table. The 5′ and 3′ non-coding regions (NCRs) of each segment have a length of 9 to 209 nt in the 5′ NCR and 30 to 112 nt in the 3′ NCR (S1 Table). At the extreme termini of each 5′/3′ NCRs of gull/MR04_RV/HUN/2014, conserved nucleotide sequence motifs were identified and analysed compared to the corresponding region of the reference strains. The common nucleotide sequences were seen in S2 Fig.

The amino acid (aa) sequences encoded by the gull/MR04_RV/HUN/2014 genome segments were compared with the corresponding homologues of the reference RV sequences belonging to groups A to D and F to J. The conservative motifs and “hotspots” analysed by Trojnar et al. [8] were looked for and aas that have changed compared to the group A RV isolate RVA/Simian-tc/ZAF/SA11-H96/1958/G3P5B[2] sequence in bold.

In the VP1, VP3 and NSP2, NSP3 proteins several conserved aa motifs could be detected with some aa changes. Of the six aa motifs described for RV VP1 protein, the motif A (YTD*VS*Q*WD*ASQHNT), the motif B (RYHGVAS*G*E*K*TTK*IGNSYANIALITTV), the motif C (THLRVDGD*D*NVVTAYTS), the motif D (MNARVKALASYTGLEMAKR), the motif E (KALASYTGLEMAK*RFIICGKIFERGA) and motif F (YPLDRLEAYSILPWP) could be found with minor/major changes and aa residues in motifs with a key role in catalysis are marked with an asterisk [8,39,40]. No aa residues interacting with the γ-phosphate of the incoming NTP and nucleosides in the P or N site could be found in motif F [39], but additional conservative motifs could be observed (AKxxSxPMx, xSxxFHVG, x represent any aa) at the N-terminal site of the protein.

After aligning the VP3 protein sequences, neither the KxTAMD nor the Kx[P/N]G motifs were found in the viral protein of the study strain [8], but aa motifs: GxxxE(S/T → GPEFES and LΩxL(S/T)NxxN (where Ω indicates an aromatic aa) → LYSISNTKN in the N7-MTase domain, aa motifs: SGHΦ (where Φ is a hydrophobic residue) → SGHL in the GTPase/RTPase domain and two HΦ(T/S)Φ motifs: HITL and HMTI in the VP3 C terminus (a predicted structural domain described in RVA) have been observed [41]. Further aa residues have also been detected in VP3 (one residue in N7-MTase domain: VRETVF → VRNTVF and four residues in the GTPase/RTPase domain: SRRSQFLR → RQVSIYNR; PMRWAK → PMRWSK and DRHSLH → EKHTRE) [41].

In the NSP2 protein, the conserved aa motif of the NTPase domain (KITxNxxxxxxDxxxxxxxAxxxxxxNxFAxIxHGxxHxR) was observed based on the structural similarity in RVA, RVC and RVD strains [8], of which the histidine triad (HIT motif: HΦHΦHΦΦ, Φ is a hydrophobic residue) were the conserved residues [40,42]. The modified HIT-like motif (HGHGHVRSV, catalytically active histidine bolded) was also detected in the study strain gull/MR04_RV/HUN/2014 NSP2.

Based on the NCBI Conserved Domain Search (CDS) the NSP3 protein of the strain gull/MR04_RV/HUN/2014 belongs to the NSP3_RV superfamily and has a similar structure to protein cd20714 (E-value: 6.20e-23). A total of 30 aa residues of the RNA-binding site were detected and within and after this site, three distinguished aa motifs (PKR, NLEKEVRMLR and NSTSDLKLS) could be identified compared to the reference protein.

Each nucleotide sequences of strain gull/MR04_RV/HUN/2014 segments were sent to the NCBI BlastX to search the translated nucleotide query in the clustered protein database (S2 Table). All segments have high (more than 85%) nucleotide identity to ruddy turnstone RV isolate MW05 (MH453863-MH453873) in RVG with significant query cover, but four segments have low sequence identity to other types of RVs (S2 Table). The segment 4 encoding the VP4 (VP5/VP8) protein(s) has 40.54% aa identity (query cover: 97%) to RVB species strain RVB/Goat-wt/USA/Minnesota-1/2016 (KY689690/ASV45170) as the best match, followed by the closest RVG strain chicken/03V0567/DEU/2003 (JQ920006/ AFL91895) with 39.45% identity (query cover: 96%). The genetic analysis of segments 5 (NSP1), 9 (VP7) and 10 (NSP4), although the segment identities were the highest with homologue segments of RVG strains RVG/turkey-wt/USA/Minnesota-1/2016 (KY689680/ASV45160), Shackleton virus (MT025062/QIS87933) and DRVG/G018-19/Australia/2019 (MW388702/UAJ21477), respectively, with a very low nt and aa identities (S2 Table).

The nt sequence analyses of the VP6 inner capsid protein of the study strain showed that it is genetically most closely related to members of RVG species, and the closest relative is a partial 570nt/190aa long VP6 sequence of the strain herring gull/H0110385/NL/2014 (KP057507, AKE33300) isolated from an FFPE cloacal bursa tissue of host Larus argentatus in the Netherlands in 2014 (Fig 1). The second closest related sequence is the Ruddy turnstone rotavirus isolate MW05 (MH453863- MH453873) isolated from a combined oropharyngeal and cloacal swab in the King Island, Australia in 2014 (Fig 1). The phylogenetic analysis of the further segments of the study strain, each segment showed a relative relationship with the corresponding segments of the isolate MW05, however, segments 4, 5, 9 and 10 are more homologous to members of the RVB group. The list of RV sequences included in this study is collected in S3 Table and the detailed phylogenetic analysis of each segment can be found in S3 Fig.

[Figure omitted. See PDF.]

Phylogenetic analysis of the internal capsid protein VP6 encoded segment 6 of strain gull/MR04_RV/HUN/2014 (PP239054, PP239059) rotavirus (RV) and all complete and partial amino acid coding nucleotide sequences of group G rotaviruses available in the NCBI GenBank database were used for the VP6 analysis. The nucleotide sequences were aligned using the Multiple Alignment using Fast Fourier Transform (MAFFT) v7.525 (2024/Mar/13) [33] and then analysed using IQ-TREE software with the following parameters: -m MFP, -alrt 1000, -T AUTO, -bb 100 [34,35,36]. According to the Bayesian Information Criterion (BIC) score the best-fit model for the VP6 dataset was TIM2 + F + R3 and for the NSP5 dataset was TIM2 + F + I + G4 models, respectively. Both the SH-aLRT test and the ultrafast bootstrap (UFBoot) values were set up to 1000 in each analysis, then the maximum likelihood consensus tree was used. The nucleotide, the amino acid accession identification numbers and the strain/isolate name of the sequences were indicated on the phylogenetic tree. The study strain gull/MR04_RV/HUN/2014 (PP239054, PP239059) sequence is marked in bold. [P]: partial sequence.

By aligning all the genome segments of the strain gull/MR04_RV/HUN/2014 to the corresponding segments of reference and the closest related RV strains using the Needleman-Wunsch alignment method, low similarity values were obtained, especially for those of the diverse four segments 4, 5, 9 and 10 mentioned above (Fig 2) [43].

[Figure omitted. See PDF.]

The main diagrams (segments 1–11) show the relationship between the reference or closest related sequences aligned to the segments of the study strain sequence using the Needleman-Wunsch (NW) Global Alignment method based on the NW scores (X-axis: NWS), the alignment GAP (Y-axis: GAP) and the percent of identity (points). The phylogenetic position of the study strain has also been indicated in the small figure. A detailed phylogenetic analysis of each segment can be found in S3 Fig.

Analysis of the merged segments of the strain gull/MR04_RV/HUN/2014, the ruddy turnstone RV isolate MW05 (RVG) and the human rotavirus B (RVB) strain Bang373 yielded surprising results. In the analysis, segments similar to the RVG strain homologated to the RVG reference sequence, while segments with low identity values (5: NSP1 and 10: NSP4, 9: VP7 and 4: VP4) showed surprising similarity to the corresponding segments of human rotavirus strain Bang373 (RVB), which may be an exciting example of reassortment between RV species (RVG/RVB) (Fig 3).

[Figure omitted. See PDF.]

37]. For the analysis, the entire genomes of the study strain gull/MR04-RV/HUN/2014 genome segments (PP239049-PP239059) rotavirus, the closest relative ruddy turnstone rotavirus isolate MW05 (MH453863-MH453873) and human rotavirus B strain Bang373 (EU490415, EU490418, AY238384, AY238385, AY238388- AY238394) were concatenated and aligned using Multiple Alignment using Fast Fourier Transform (MAFFT) method (https://www.ebi.ac.uk/jdispatcher/msa/mafft). The aligned nucleotide sequences were analysed with the SimPlot program based on the following setup: Window: 500 bp, Step: 20 bp, GapStrip: On, F84 (“Maximum Likelihood”), T/t: 2,0). The segments coding for non-structural and capsid proteins (with frames) and possible reassortant segments (with arrows) are marked in the figure. RVG: rotavirus group G, RVB: rotavirus group B, S1–S11: segment 1–11.

Characterization of Avastrovirus 2 (AAstV-2) strain gull/MR04_AAstV/HUN/2014

The determined partial genome of astrovirus strain gull/MR04_AAstV/HUN/2014 (PP239060) is 3,448nt long (partial ORF1b is 1,137nt/378aa and the full-length ORF2 is 2,139 nt/712aa) (Fig 4A). The partial ORF1b encoded polyprotein contains the putative aa motifs (DWTRFD, GNPSG, YGDD, and FGMWVK) of the catalytic domain of the astroviral RdRp enzyme [44,45]. The sub-genomic RNA nt sequence motif (UUUGGAGNGGNGGACCNAAN_4-11AUGNC) is identifiable in the study strain, the modified sequence of which is as follows UGAGNGGNGGACCNAAN₁₄AUGGC (N: any of four nucleotides, the ORF2 initiation codon is underlined). The ORF1b is non-overlapped with ORF2, a 24nt-long spacer is observed, but the “CCGA” AAstV pentamer [46] is not present. The conserved nt sequence of stem-loop II (s2m) mobile genetic elements could not be identified in either the partial ORF1b or the 3’UTR. The closest related sequence of the partial ORF1a in the GenBank database was an Astroviridae sp. isolate prf038ast1 (MT138011) sequence from a bird metagenome with 66.55% nt identity (query cover: 77%) and the ORF1b protein of ruddy turnstone astrovirus isolate Prototype-RtAstV (MK189093, QCP68856) sequence with 64.55% (query cover: 100%) aa identity.

[Figure omitted. See PDF.]

A) The genome organization of the partial sequence of gull/MR04_AAstV/HUN/2014 astrovirus was constructed based on the Duck astrovirus strain CPH (KJ020899). The conserved nucleotide sgRNA signal and the conserved amino acid motifs in the ORF1b and ORF2 proteins are highlighted. B: basic amino acids, A: acidic amino acids. B) The cluster analysis of the study strain gull/MR04_AAstV/HUN/2014 (PP239060) and the representative strains in AAstV-2 group 6. The cluster analysis was done with an R script (in this study) and “seqinr”, “msa”, “cluster”, “factoextra”, “randomcoloR”, “ggrepel” and “ape” packages using R/RStudio. The script of the analysis can be provided as requested. Arrow shows the most closely related virus sequence AAstV-2 isolate KH08-1279 (JX985649) with the BlastP identity and the Needleman-Wunsch (NW) Global alignment score results. The study strain is coloured red. The detailed phylogenetic analysis of avian astroviruses can be found in the appendix S4 Fig.

The ORF2 capsid protein of the strain gull/MR04_AAstV/HUN/2014 is likely to be an astrovirus capsid protein based on the NCBI CDS and BlastP identity analysis. The ORF2 protein is separated into a conserved particle assembly domain (AD) and a variable receptor-interaction domain (RID) (Fig 4A). Based on the NCBI CDS and a comparative sequence analysis with the reference/representative types of avian nephritis viruses (ANV) [46], the first 390 aa elements build the inner and outer ORF2 capsid core contain conservative elements/motifs while the spike of the ORF2 protein is variable and has no detected conserved motif and has low sequence identity. Examining the nt sequence encoding the ORF2 capsid protein, the closest related sequence was an Astroviridae sp. isolate tom152ast2 (MN920667) sequence from a bird metagenome with 66.51% nt identity (query cover: 52%) and the corresponding protein of Avastrovirus 2 isolate KH08-1279 (JX985649, AFW05403) from an Ardeola bird sp. with 52.76% aa identity (query cover: 83%) (Fig 4B). Based on the phylogenetic analysis performed on all capsid (complete and partial) proteins (N = 821) of genus Avastrovirus (taxid:249589) the study strain belongs to the avian nephritis virus (ANV) group in Avastrovirus 2 (AAstV-2) (S4 Fig). According to the genetic identity and the aa-based cluster analysis, the study strain is closely related to isolate KH08-1279 (JX985649, AFW05403) and strain red-necked avocet/MW08/Interior (MH453800, AXF38644) viruses which are classified in group 6 of ANVs (Fig 4B, S4 Fig) [46].

Using the screening primer pairs (MR04-GVP6-SCR-R/ MR04-GVP6-SCR-F; MR04-AvAst2-SCR-R/ MR04-AvAst2-F) neither RV nor astrovirus could be found in the further faecal samples of common black-headed (N=4) and yellow-legged (N = 2) gulls.

Discussion

The species-rich Laridae bird family includes the genera gulls, terns, noddies, skimmers and kittiwakes, which are primarily seabirds which are widespread in the world, but species that can also be found in Western Europe, e.g.,: European herring gull (Larus argentatus) and in Northern Europe, e.g.,: common gull (Larus canus) and in Central Europe, e.g.,: black-headed (Chroicocephalus ridibundus) and yellow-legged (Larus michahellis) gulls including Hungary [1]. Considering their lifestyle, these birds may come into close contact with humans, and the viral metagenomic analysis of their excrement can provide useful information about the prevalence of viruses excreted in faeces. Viral metagenomics using NGS methods can reveal many potential new pathogens, however, this discovery is only a meaningless and unmanageable data set in the absence of further detailed analysis [28,43].

In this study, after evaluating the mass sequencing data, we focused on viral sequences occurring in vertebrates, among which viral families Sedoreoviridae and Astroviridae were interesting for deeper investigation. Due to the low percentage of reads (less than 0.3% per virus family) and the uncertain taxonomic classification (less than 30% identity) of reads belonging to the families Coronaviridae, Picobirnaviridae, Retroiridae, Caliciviridae and Picornaviridae, were not analyzed in this study. The data obtained during mass sequencing was filtered using the Kaiju method based on the NCBI nr protein database (accessed 03.27.2023.) and then confirmed using the diamond method. The obtained reads were aligned to the closest related sequence, and then verified with sequence-specific primers, after we reported on a potential interspecies reassortant RVG strain gull/MR04_RV/HUN/2014 and an avian astrovirus (AAstV-II) strain gull/MR04_AAstV/HUN/2014.

The RV identified in the black-headed gull faecal sample is an exciting finding from several points of view. Firstly, according to the rule defined by the International Committee on Taxonomy of Viruses (ICTV) Sedoreoviridae study group, viruses whose VP6, encoding in segment 6, amino acid sequence has more than 53% identity belong to the same species within the RV genus (https://ictv.global/report/chapter/sedoreoviridae/sedoreoviridae/rotavirus) [3]. Based on this criterion, the black-headed gull RV strain gull/MR04_RV/HUN/2014 belongs to the genus Rotavirus and species Rotavirus G (RVG), but with a low (66%) nucleotide identity value. Secondly, however, most of the RV conserved motifs identified by Trojnar et al. [8] in the functionally important proteins (RdRp/VP1, Cap/VP3, NSP2) were present in the study RV strain, but most parts of the viral genome show significant differences from members of the RVG, especially the segments encoding the NSP1, NSP4, VP7 and VP4 proteins.

The first RVG species was identified from chickens (Gallus gallus) [47] as potentially causing runting and stunting syndrome [48] and then an RVG species was also described from turkeys (Meleagris gallopavo) [49]. To date, RVG species have been identified in many other species of wild birds, such as ducks (Anatidae), pigeons (Columbidae), partridges (Perdicinae), gulls (Laridae), and interestingly pigs (Sus scrofa domesticus) (Fig 1). Based on the analysis of segment 6 encoding the VP6 capsid protein, the closest relative was the RV identified from herring gull in the Netherlands in 2014. Unfortunately, only the partial VP6 sequence of isolate herring gull/H01-10385/NL/2014 is available in the GenBank database [50], so we could not use this isolate for a deeper analysis covering additional segments. Comparing the corresponding regions of viruses with a complete genome, the closest relative of the study strain was the RVG isolate MW05 (MH453868) described from ruddy turnstone waterfowl in Australia in 2014. Until August 2023, RVG types had been described only from avian species, but Dias and colleagues found RVG group rotavirus strains RVG/Pig-wt/BRA/BP52/2016 (OR067849) and RVG/Pig-wt/BRA/BP56/2016 (OR067850) in pigs in Brazil in 2016 (unpublished, only sequences are available). The RVG viruses identified in pigs are closely related to RVG types identified in chickens, including the reference sequence of the RVG group strain chicken/03V0567/DEU/2003 (NC_021587). Unfortunately, only the partial NSP5 portion of the two pig RVG strains mentioned above is available in the GenBank database, which we could not use for our VP6-based sequence analysis.

In a detailed analysis, the segments encoding additional structural and non-structural proteins, very low nt identities (with a high overlap value) were observed, in segments 4, 5, 9 and 10 compared to any of the known corresponding RV sequences. Among these four segments, segments 5 and 10 encode non-structural NSP1 and NSP4 proteins respectively, of which the nearly 500-aa-long NSP1 protein has a potential role in inhibiting the host’s innate immunity by breaking down the key factors required to activate host interferon production [51]. The aa similarity of the NSP1 protein to members of RVG proved to be at most 35%, however, the phylogenetic analysis of this segment with members of other RV groups showed a clear relationship with the human RVB strain Bang373 (RVB) (Fig 2). The NSP4 of segment 10 is a relatively small, 175-aa-long protein that functions as a viroporin [52]. The nt sequence of segment 10 had no significant match with any of the known RVG viruses and the aa level has the closest to RVG strain DRVG/G018-19/Australia/2019 (UAJ21477) but with an extremely low aa identity value. Interestingly, in the nt level the NSP4 protein of the study strain shows significant identity (with no significant query cover) with RVB strain RVB/Pig-wt/VNM/14176_8/NSP4 (KX362398) described in porcine [53] (S2 Table). The VP7 protein of segment 9 shows low nt and aa identity to RVG viruses, while VP4 encoded in segment 4 has no relationship to RVG viruses at either nt or aa level. In contrast, the VP4 is related to the corresponding protein of the RVB strain JN311 (KU562898) isolated from human stool in Bangladesh [54]. The low-level nt and aa identities observed in segments 4, 5, 9 and 10 of the study strain compared to RVG lead to the conclusion that these segments have originated from a currently unknown RV related distantly to RVB which could make the study RV strain a potential reassortant virus and may be able to infect and replicate successfully among in RVB viruses. Based on the analysis of the concatenated nucleotides of all segments of the study strain gull/MR04_RV/HUN/2014, and the closest relative partners ruddy turnstone RV isolate MW05 (RVG) and human rotavirus B strain Bang373 (RVB), it can be assumed that the RV identified in this study is, in fact, an interspecies (RVG/RVB) reassortant virus capable of infecting birds (RVG), rodents (RVB), pigs (RVB) and possibly the human (RVB) population.

Group B rotaviruses (RVBs) have been identified in various species, including humans [54; 55], swine [56–58], bovines [58], caprine [59], equines [58,60,61], lambs [62], and goats [62], and rodents [63]. Human RVB infections or outbreaks have been reported in Asia, specifically in Nepal [64], India [65,66], and Bangladesh [54,67], as well as in South Africa [68]. In other parts of the world, swine are considered the primary reservoir for RVBs. But in Europe, no human RVB infections or outbreaks have been documented. However, a few studies have identified RVBs in animal hosts, particularly swine [56–58], bovines [58], equines [58], lambs [62], and goats [62], albeit in limited geographic regions. Despite these findings, there is no formal epidemiological surveillance system in the EU/EEA for monitoring rotavirus infections or circulating strains, as noted by the European Centre for Disease Prevention and Control (ECDC) (https://www.ecdc.europa.eu/en/rotavirus-infection/facts). Most available data on RVBs in Europe come from isolated studies, highlighting the need for further research to better understand the significance of RVB in animal and human diarrheal diseases.

Point mutations, intragenic [13] and intergenic (inter-lineage/-sub-lineage) recombination [14,15,16] and reassortation (between single or multiple genome segments) [17] evolutionary mechanisms as a whole are the driving forces of strain variation that can facilitate the adaptation of the segmented RV with a double-stranded RNA genome to the host [11–13,69]. However, point mutations and intragenic and intergenic recombination are well-studied mechanisms, as a new observation, a recombination phenomenon was also described between segments of different RV species in two cases [18,19]. Both studies reported partial or complete replacement of the NSP3 protein (segment 7) of porcine RV belonging to the RVH group with the corresponding segment of porcine RV belonging to the RVC species, which was the first evidence of interspecies recombination. In 2017, two turkey RVs (strains Turkey-wt/USA/Minnesota-1/2016 and Turkey-wt/USA/Minnesota-2/2016) were described in the RVG group, a possible recombination event was detected in VP6 and NSP3 proteins compared to the member of the RVB human strains RVB/Human-wt/BAN/Bang117/2002 (GU391305) and RVB/Human-wt/JAP/MMR-B1/2007 (GU370059) corresponding partial regions [70].

The mechanism of reassortment has been observed between the same RV species, especially in group A RVs, is also a well-known phenomenon [71,72], where the two most important neutralization antigens VP4 (encoded on segment 4) and VP7 (encoded on segment 9) proteins, as well as NSP4 (encoded on segment 10) and NSP5/6 (encoded on segment 11) proteins, are most likely to be replaced [73]. The VP4 binds the virion to the host’s epithelial cell receptors and delivers it into the host cell [74]. Presumably, the binding of RVs and their entry into the host cell must be the result of multiple consecutive contacts between the outer capsid proteins VP4 and VP7 and cell receptors, in which a significant change can modify the virulence of the virus, and even widen the possible host spectrum or other host specificity can result [74]. The NSP4 protein is important for viral morphogenesis and has enterotoxin activity [73], and the NSP5/6 proteins are useful for tracing the species origin of RVs [73].

The molecular process of reassortment is not fully understood, but, likely, reassortment does not require physical proximity to the parental genomes and there is growing evidence to support the selective packaging model, that uses a signal-driven selective segment incorporation process to shape newly generated viral particles, as seen in the influenza virus [69,75–77]. The virus assortment is likely coordinated by RNA–RNA interactions and a selective packaging signal/ cis-acting element according to the strategy of the “Concerted model” [78]. The conserved sequences at the 5’ and 3’ ends of the rotavirus segments are so-called “cis-acting elements (CRE)” with different functional content essential for selective packaging [78], but the CRE at the 3’end serves as a critical polymerase recognition element and provides a template for genome replication [78].

In our study, the genetically and phylogenetically analysed strain gull/MR04_RV/HUN/2014 is closely related to RVG group viruses based on the proteins VP1-VP4, VP6, NSP2, NSP3, NSP5, while based on the non-structural proteins NSP1, NSP4 and the structural proteins VP7, VP4, it is more related to the human RVB strain Bang373. The RVG species described in this communication is presumably an interspecies reassortant virus, which has been supported by several indirect observations so far: a) the RVG and RVB virus species belong to the same RVB-like clade (clade 2) and are juxtaposed in the phylogenetic analysis of all 11 segments as seen in S3 Fig [6,9] and phylogenetically closely related viruses are genetically similar and, in case of a co-infection, are more likely to produce infectious virus particles by reassortation. b) Comparing the 5’ and 3’ consensus ends of the study strain, it is worth directing attention to segments 2 (VP2), 4 (VP4), 5 (NSP1) and 10 (NSP4) where a high degree match with the consensus elements of the RVB group was found (S2 Fig). c) Two cases of putative recombination events between turkey RVG and human RVB species have been reported, involving the segments encoding VP6 and NSP3 proteins [70], in the study strain there was not a pointwise or intermittent excision and replacement of certain parts of the segments, but a replacement of entire segments could have happened. d) Several cases of rotaviruses belonging to different RV groups have been described that were present in the same host, causing co-infection [79–81]. Interestingly, in a Chinese study, a survey of pig herds showed the co-occurrence of RV belonging to other species in many cases, where some type of co-infection was found in 40% of the examined samples [81]. Among the test results, the most surprising result was the 2% quadruple (RVA+RVB+RVC+RVH) RV infection [81]. e) Finally, the RVG type has been detected only in avian species but it has been described in 2 cases in pigs (unpublished, only nucleotide sequences are available).

The co-circulation of different RV species in the same host species can facilitate the exchange of genes or even segments between different RV strains. Considering the lifestyle of the black-headed gulls described in this study, it is likely that they could have been infected with the RVB virus in human (contaminated with faeces) waste (food?) or in an intermediate host (swine?), which could have recombined during pathogenesis with the RVG virus typically found in birds.

The astrovirus identified in the analysed black-headed gull faecal sample phylogenetically belongs to the genus Avastrovirus species Avastrovirus 2 (AAstV-2) and avian nephritis viruses (ANVs). Unfortunately, due to the small amount of faecal sample collected from black-headed gull and the importance of the RV presented in this study, it was not possible to determine the complete genome of the strain gull/MR04_AAstV/HUN/2014. Unfortunately, only the partial ORF1b and the complete ORF2 capsid sequences, which exhibited properties characteristic of the AAstV (the nt sequence of the conserved sgRNA promoter between ORF1b/ORF2, the key conserved aa motifs characteristic of the polymerase enzyme) was analyzed. The classification of AAstVs is currently (since 2009) based on the hosts from which they have been isolated and the re-grouping of the known sequences is currently underway, with analysis of the entire capsid protein of the virus likely to be one of the pillars (https://ictv.global/report_9th/RNApos/Astroviridae). Currently, there are three official AAstV species (AAstV 1-3) and within these, many viral sequences await classification [4]. Based on the survey by Kariithi et al. [46], the AAstV-2 species, avian nephritis viruses (ANVs) group is subdivided into 11 additional sub-genotypes (groups 1–11). The nt genetic distance within each virus sub-genotype group varies between 0.204 and 0.284. The mean value calculated in group 6 is 0.418 [24,25,46], which raises the possibility of forming further groups, and we calculate that four more groups (group 6/A: 0.36, group 6/B: 0.43, group 6/C: 0.19, and group 6/D: 0.22) could be formed and visualized in S4 Fig, among which the gull/MR04_AAstV/HUN/2014 study strain presented here (group 6/B) would also occupy an important place.

The astrovirus included in the study is phylogenetically related to strain KH08-1279 virus isolated from pond herons (Ardeola sp.) in China and strain isolate AstV/Red-necked Avocet/MW08/Interior virus isolated from red-necked avocet (Recurvirostra novaehollandiae) in Australia. Both astroviruses were described in the faecal samples of waterfowl whose related species, e.g.,: little egret (Egretta garzetta), and pied avocet (Recurvirostra avosetta) also have known populations in Hungary, and this observation assumes that Hungarian herons, gulipans, and mostly gulls can be potential carriers and spreaders of AAstVs due to their migratory and aquatic lifestyle.

In conclusion, the wild bird populations are greatly affected by urbanization, which narrows the original habitats of waterfowl, forcing the majority of species to adapt to urban lifestyles. Paradoxically, gulls have been able to adapt to urban life, where landfills and polluted waterways allow them an expanded food source and at the same time transmission of human pathogens. Due to the very small number of study samples, likely, one piece of the whole puzzle regarding the diversity of RVs and astroviruses is found in nature. Further molecular surveillance for gulls or other related species adapted to human environments (crows, pigeons, etc.) and even more waterfowl can provide material for the discovery of already known or novel RVs, astroviruses and other novel viruses with RNA genomes.

Supporting information

S1 Fig. Classification of reads derived from viral metagenomic analysis.

https://doi.org/10.1371/journal.pone.0317400.s001

(DOCX)

S2 Fig. Comparison of 5’- and 3’-terminal non-coding regions (5’/3’ NCRs) sequences.

https://doi.org/10.1371/journal.pone.0317400.s002

(DOCX)

S3 Fig. Phylogenetic analysis of the corresponding segments of the study strain, closest related representative sequences and reference rotavirus strains.

https://doi.org/10.1371/journal.pone.0317400.s003

(DOCX)

S4 Fig. Phylogenetic and Cluster analysis of avastroviruses.

https://doi.org/10.1371/journal.pone.0317400.s004

(DOCX)

S1 Table. The genome organization of group G rotavirus strain gull/MR04-RV/HUN/2014.

https://doi.org/10.1371/journal.pone.0317400.s005

(DOCX)

S2 Table. Nucleotide and amino acid identity comparison of strain gull/MR04-RV/HUN/2014.

https://doi.org/10.1371/journal.pone.0317400.s006

(DOCX)

S3 Table. List of rotavirus sequences included in the analysis.

https://doi.org/10.1371/journal.pone.0317400.s007

(DOCX)

Acknowledgments

The research was performed in collaboration with the Genomics and Bioinformatics Core Facility at the Szentágothai Research Centre of the University of Pécs. Bioinformatics infrastructure was supported by ELIXIR Hungary (https://elixir-hungary.org).

References

1. 1. Martín-Vélez V, Montalvo T, Afán I, Sánchez-Márquez A, Aymí R, Figuerola J, et al. Gulls living in cities as overlooked seed dispersers within and outside urban environments. Sci Total Environ. 2022;823:153535. pmid:35104514

* View Article

* PubMed/NCBI

* Google Scholar

2. 2. Wille M, Eden J-S, Shi M, Klaassen M, Hurt AC, Holmes EC. Virus-virus interactions and host ecology are associated with RNA virome structure in wild birds. Mol Ecol. 2018;27(24):5263–78. pmid:30375075

* View Article

* PubMed/NCBI

* Google Scholar

3. 3. Matthijnssens J, Attoui H, Bányai K, Brussaard CPD, Danthi P, Del Vas M, et al. ICTV virus taxonomy profile: sedoreoviridae 2022. J Gen Virol. 2022;103(10):001782. pmid:36215107

* View Article

* PubMed/NCBI

* Google Scholar

4. 4. Cortez V, Meliopoulos VA, Karlsson EA, Hargest V, Johnson C, Schultz-Cherry S. Astrovirus biology and pathogenesis. Annu Rev Virol. 2017;4(1):327–48. pmid:28715976

* View Article

* PubMed/NCBI

* Google Scholar

5. 5. Ogden KM, Johne R, Patton JT. Rotavirus RNA polymerases resolve into two phylogenetically distinct classes that differ in their mechanism of template recognition. Virology. 2012;431(1–2):50–7. pmid:22687427

* View Article

* PubMed/NCBI

* Google Scholar

6. 6. Kindler E, Trojnar E, Heckel G, Otto PH, Johne R. Analysis of rotavirus species diversity and evolution including the newly determined full-length genome sequences of rotavirus F and G. Infect Genet Evol. 2013;14:58–67. pmid:23237956

* View Article

* PubMed/NCBI

* Google Scholar

7. 7. Matthijnssens J, Otto PH, Ciarlet M, Desselberger U, Van Ranst M, Johne R. VP6-sequence-based cutoff values as a criterion for rotavirus species demarcation. Arch Virol. 2012;157(6):1177–82. pmid:22430951

* View Article

* PubMed/NCBI

* Google Scholar

8. 8. Trojnar E, Otto P, Roth B, Reetz J, Johne R. The genome segments of a group D rotavirus possess group A-like conserved termini but encode group-specific proteins. J Virol. 2010;84(19):10254–65. pmid:20631147

* View Article

* PubMed/NCBI

* Google Scholar

9. 9. Johne R, Schilling-Loeffler K, Ulrich RG, Tausch SH. Whole genome sequence analysis of a prototype strain of the novel putative rotavirus species L. Viruses. 2022;14(3):462. pmid:35336869

* View Article

* PubMed/NCBI

* Google Scholar

10. 10. Dhama K, Saminathan M, Karthik K, Tiwari R, Shabbir MZ, Kumar N, et al. Avian rotavirus enteritis – an updated review. Vet Q. 2015;35(3):142–58. pmid:25917772

* View Article

* PubMed/NCBI

* Google Scholar

11. 11. Jere KC, Mlera L, Page NA, van Dijk AA, O’Neill HG. Whole genome analysis of multiple rotavirus strains from a single stool specimen using sequence-independent amplification and 454® pyrosequencing reveals evidence of intergenotype genome segment recombination. Infect Genet Evol. 2011;11(8):2072–82. pmid:22019521

* View Article

* PubMed/NCBI

* Google Scholar

12. 12. Jere KC, Chaguza C, Bar-Zeev N, Lowe J, Peno C, Kumwenda B, et al. Emergence of double- and triple-gene reassortant G1P[8] rotaviruses possessing a DS-1-like backbone after rotavirus vaccine introduction in Malawi. J Virol. 2018;92(3):e01246-17. pmid:29142125

* View Article

* PubMed/NCBI

* Google Scholar

13. 13. Hoxie I, Dennehy JJ. Intragenic recombination influences rotavirus diversity and evolution. Virus Evol. 2020;6(1):vez059. pmid:31949920

* View Article

* PubMed/NCBI

* Google Scholar

14. 14. Suzuki Y, Gojobori T, Nakagomi O. Intragenic recombinations in rotaviruses. FEBS Lett. 1998;427(2):183–7. pmid:9607308

* View Article

* PubMed/NCBI

* Google Scholar

15. 15. Parra GI, Espínola EE, Amarilla AA, Stupka J, Martinez M, Zunini M, et al. Diversity of group A rotavirus strains circulating in Paraguay from 2002 to 2005: detection of an atypical G1 in South America. J Clin Virol. 2007;40(2):135–41. pmid:17720620

* View Article

* PubMed/NCBI

* Google Scholar

16. 16. Phan TG, Okitsu S, Maneekarn N, Ushijima H. Evidence of intragenic recombination in G1 rotavirus VP7 genes. J Virol. 2007;81(18):10188–94. pmid:17609273

* View Article

* PubMed/NCBI

* Google Scholar

17. 17. Matthijnssens J, Rahman M, Van Ranst M. Two out of the 11 genes of an unusual human G6P[6] rotavirus isolate are of bovine origin. J Gen Virol. 2008;89(Pt 10):2630–5. pmid:18796733

* View Article

* PubMed/NCBI

* Google Scholar

18. 18. Suzuki T, Inoue D. Full genome-based genotyping system for rotavirus H and detection of potential gene recombination in nonstructural protein 3 between porcine rotavirus H and rotavirus C. J Gen Virol. 2018;99(12):1582–9. pmid:30328806

* View Article

* PubMed/NCBI

* Google Scholar

19. 19. Krasnikov N, Yuzhakov A. Interspecies recombination in NSP3 gene in the first porcine rotavirus H in Russia identified using nanopore-based metagenomic sequencing. Front Vet Sci. 2023;10:1302531. pmid:38116510

* View Article

* PubMed/NCBI

* Google Scholar

20. 20. Gabbay YB, Leite JPG, Oliveira DS, Nakamura LS, Nunes MRT, Mascarenhas JDP, et al. Molecular epidemiology of astrovirus type 1 in Belém, Brazil, as an agent of infantile gastroenteritis, over a period of 18 years (1982–2000): identification of two possible new lineages. Virus Res. 2007;129(2):166–74. pmid:17714822

* View Article

* PubMed/NCBI

* Google Scholar

21. 21. Finkbeiner SR, Holtz LR, Jiang Y, Rajendran P, Franz CJ, Zhao G, et al. Human stool contains a previously unrecognized diversity of novel astroviruses. Virol J. 2009;6:161. pmid:19814825

* View Article

* PubMed/NCBI

* Google Scholar

22. 22. Tan Z, Zhai H, Sun R, Sun Z, Xie R, Zhang L, et al. Complete genome sequence and phylogenetic analysis of a goose astrovirus isolate in China. Braz J Microbiol. 2023;54(1):427–34. pmid:36327040

* View Article

* PubMed/NCBI

* Google Scholar

23. 23. Yin L, Zhou Q, Mai K, Huang J, Yan Z, Wei X, et al. Isolation and characterization of a novel chicken astrovirus in China. Poult Sci. 2021;100(9):101363. pmid:34352410

* View Article

* PubMed/NCBI

* Google Scholar

24. 24. Donato C, Vijaykrishna D. The broad host range and genetic diversity of mammalian and avian astroviruses. Viruses. 2017;9(5):102. pmid:28489047

* View Article

* PubMed/NCBI

* Google Scholar

25. 25. De Benedictis P, Schultz-Cherry S, Burnham A, Cattoli G. Astrovirus infections in humans and animals - molecular biology, genetic diversity, and interspecies transmissions. Infect Genet Evol. 2011;11(7):1529–44. pmid:21843659

* View Article

* PubMed/NCBI

* Google Scholar

26. 26. Reuter G, Boros Á, Mátics R, Kapusinszky B, Delwart E, Pankovics P. A novel avian-like hepatitis E virus in wild aquatic bird, little egret (Egretta garzetta), in Hungary. Infect Genet Evol. 2016;46:74–7. pmid:27876615

* View Article

* PubMed/NCBI

* Google Scholar

27. 27. Boros Á, Albert M, Urbán P, Herczeg R, Gáspár G, Balázs B, et al. Unusual “Asian-origin” 2c to 2b point mutant canine parvovirus (Parvoviridae) and canine astrovirus (Astroviridae) co-infection detected in vaccinated dogs with an outbreak of severe haemorrhagic gastroenteritis with high mortality rate in Hungary. Vet Res Commun. 2022;46(4):1355–61. pmid:36129562

* View Article

* PubMed/NCBI

* Google Scholar

28. 28. Boros Á, Pankovics P, László Z, Urbán P, Herczeg R, Gáspár G, et al. The genomic and epidemiological investigations of enteric viruses of domestic caprine (Capra hircus) revealed the presence of multiple novel viruses related to known strains of humans and ruminant livestock species. Microbiol Spectr. 2023;11(6):e0253323. pmid:37823638

* View Article

* PubMed/NCBI

* Google Scholar

29. 29. Sverdlov E, Azhikina T. Primer walking. In: eLS, (Ed.). 2005. Available from:

* View Article

* Google Scholar

30. 30. Boros Á, Pankovics P, Simmonds P, Reuter G. Novel positive-sense, single-stranded RNA (+ssRNA) virus with di-cistronic genome from intestinal content of freshwater carp (Cyprinus carpio). PLoS One. 2011;6(12):e29145. pmid:22195010

* View Article

* PubMed/NCBI

* Google Scholar

31. 31. Wakuda M, Pongsuwanna Y, Taniguchi K. Complete nucleotide sequences of two RNA segments of human picobirnavirus. J Virol Methods. 2005;126(1–2):165–9. pmid:15847933

* View Article

* PubMed/NCBI

* Google Scholar

32. 32. Reuter G, Boros A, Delwart E, Pankovics P. Novel seadornavirus (family Reoviridae) related to Banna virus in Europe. Arch Virol. 2013;158(10):2163–7. pmid:23624680

* View Article

* PubMed/NCBI

* Google Scholar

33. 33. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6. pmid:28968734

* View Article

* PubMed/NCBI

* Google Scholar

34. 34. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530–4. pmid:32011700

* View Article

* PubMed/NCBI

* Google Scholar

35. 35. Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9. pmid:28481363

* View Article

* PubMed/NCBI

* Google Scholar

36. 36. Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35(2):518–22. pmid:29077904

* View Article

* PubMed/NCBI

* Google Scholar

37. 37. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73(1):152–60. pmid:9847317

* View Article

* PubMed/NCBI

* Google Scholar

38. 38. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38(7):3022–7. pmid:33892491

* View Article

* PubMed/NCBI

* Google Scholar

39. 39. McDonald SM, Aguayo D, Gonzalez-Nilo FD, Patton JT. Shared and group-specific features of the rotavirus RNA polymerase reveal potential determinants of gene reassortment restriction. J Virol. 2009;83(12):6135–48. pmid:19357162

* View Article

* PubMed/NCBI

* Google Scholar

40. 40. Vásquez-del Carpió R, Morales JL, Barro M, Ricardo A, Spencer E. Bioinformatic prediction of polymerase elements in the rotavirus VP1 protein. Biol Res. 2006;39(4):649–59. pmid:17657346

* View Article

* PubMed/NCBI

* Google Scholar

41. 41. Ogden KM, Snyder MJ, Dennis AF, Patton JT. Predicted structure and domain organization of rotavirus capping enzyme and innate immune antagonist VP3. J Virol. 2014;88(16):9072–85. pmid:24899176

* View Article

* PubMed/NCBI

* Google Scholar

42. 42. Vasquez-Del Carpio R, Gonzalez-Nilo FD, Riadi G, Taraporewala ZF, Patton JT. Histidine triad-like motif of the rotavirus NSP2 octamer mediates both RTPase and NTPase activities. J Mol Biol. 2006;362(3):539–54. pmid:16934294

* View Article

* PubMed/NCBI

* Google Scholar

43. 43. Pankovics P, Boros Á, László Z, Szekeres S, Földvári G, Altan E, et al. Genome characterization, prevalence and tissue distribution of astrovirus, hepevirus and norovirus among wild and laboratory rats (Rattus norvegicus) and mice (Mus musculus) in Hungary. Infect Genet Evol. 2021;93:104942. pmid:34044191

* View Article

* PubMed/NCBI

* Google Scholar

44. 44. Koonin EV. The phylogeny of RNA-dependent RNA polymerases of positive-strand RNA viruses. J Gen Virol. 1991;72(Pt 9):2197–206. pmid:1895057

* View Article

* PubMed/NCBI

* Google Scholar

45. 45. Méndez E, Murillo A, Velázquez R, Burnham A, Arias CF. Replication cycle of astroviruses. Astrovirus Res Ideas. 2012;7:19–45.

* View Article

* Google Scholar

46. 46. Kariithi HM, Volkening JD, Chiwanga GH, Pantin-Jackwood MJ, Msoffe PLM, Suarez DL. Genome sequences and characterization of chicken astrovirus and avian nephritis virus from tanzanian live bird markets. Viruses. 2023;15(6):1247. pmid:37376547

* View Article

* PubMed/NCBI

* Google Scholar

47. 47. McNulty MS, Todd D, Allan GM, McFerran JB, Greene JA. Epidemiology of rotavirus infection in broiler chickens: recognition of four serogroups. Arch Virol. 1984;81(1–2):113–21. pmid:6331344

* View Article

* PubMed/NCBI

* Google Scholar

48. 48. Otto P, Liebler-Tenorio EM, Elschner M, Reetz J, Löhren U, Diller R. Detection of rotaviruses and intestinal lesions in broiler chicks from flocks with runting and stunting syndrome (RSS). Avian Dis. 2006;50(3):411–8. pmid:17039842

* View Article

* PubMed/NCBI

* Google Scholar

49. 49. Kang SY, Nagaraja KV, Newman JA. Primary isolation and identification of avian rotaviruses from turkeys exhibiting signs of clinical enteritis in a continuous MA 104 cell line. Avian Dis. 1986;30(3):494–9. pmid:3021099

* View Article

* PubMed/NCBI

* Google Scholar

50. 50. Bodewes R, van Run PRWA, Schürch AC, Koopmans MPG, Osterhaus ADME, Baumgärtner W, et al. Virus characterization and discovery in formalin-fixed paraffin-embedded tissues. J Virol Methods. 2015;214:54–9. pmid:25681526

* View Article

* PubMed/NCBI

* Google Scholar

51. 51. Hou G, Zeng Q, Matthijnssens J, Greenberg HB, Ding S. Rotavirus NSP1 contributes to intestinal viral replication, pathogenesis, and transmission. mBio. 2021;12(6):e0320821. pmid:34903043

* View Article

* PubMed/NCBI

* Google Scholar

52. 52. Pham T, Perry JL, Dosey TL, Delcour AH, Hyser JM. The rotavirus NSP4 viroporin domain is a calcium-conducting ion channel. Sci Rep. 2017;7:43487. pmid:28256607

* View Article

* PubMed/NCBI

* Google Scholar

53. 53. Phan MVT, Anh PH, Cuong NV, Munnink BBO, van der Hoek L, My PT, et al. Unbiased whole-genome deep sequencing of human and porcine stool samples reveals circulation of multiple groups of rotaviruses and a putative zoonotic infection. Virus Evol. 2016;2(2):vew027. pmid:28748110

* View Article

* PubMed/NCBI

* Google Scholar

54. 54. Ahmed MU, Kobayashi N, Wakuda M, Sanekata T, Taniguchi K, Kader A, et al. Genetic analysis of group B human rotaviruses detected in Bangladesh in 2000 and 2001. J Med Virol. 2004;72(1):149–55. pmid:14635024

* View Article

* PubMed/NCBI

* Google Scholar

55. 55. Hung T, Chen GM, Wang CG, Chou ZY, Chao TX, Ye WW, et al. Rotavirus-like agent in adult non-bacterial diarrhoea in China. Lancet. 1983;2(8358):1078–9. pmid:6138617

* View Article

* PubMed/NCBI

* Google Scholar

56. 56. Smitalova R, Rodak L, Smid B, Psikal I. Detection of nongroup A rotaviruses in faecal samples of pigs in the Czech Republic. Vet Med (Praha). 2009;54(1):12–8.

* View Article

* Google Scholar

57. 57. Vidal A, Martín-Valls GE, Tello M, Mateu E, Martín M, Darwich L. Prevalence of enteric pathogens in diarrheic and non-diarrheic samples from pig farms with neonatal diarrhea in the North East of Spain. Vet Microbiol. 2019;237:108419. pmid:31585655

* View Article

* PubMed/NCBI

* Google Scholar

58. 58. Otto PH, Rosenhain S, Elschner MC, Hotzel H, Machnowska P, Trojnar E, et al. Detection of rotavirus species A, B and C in domestic mammalian animals with diarrhoea and genotyping of bovine species A rotavirus strains. Vet Microbiol. 2015;179(3–4):168–76. pmid:26223422

* View Article

* PubMed/NCBI

* Google Scholar

59. 59. Chen F, Knutson TP, Ciarlet M, Sturos M, Marthaler DG. Complete genome characterization of a rotavirus B (RVB) strain identified in Alpine goat kids with enteritis reveals inter-species transmission with RVB bovine strains. J Gen Virol. 2018;99(4):457–63. pmid:29517476

* View Article

* PubMed/NCBI

* Google Scholar

60. 60. Uprety T, Sreenivasan CC, Hause BM, Li G, Odemuyiwa SO, Locke S, et al. Identification of a ruminant origin group B rotavirus associated with diarrhea outbreaks in foals. Viruses. 2021;13(7):1330. pmid:34372536

* View Article

* PubMed/NCBI

* Google Scholar

61. 61. Sreenivasan CC, Naveed A, Uprety T, Soni S, Jacob O, Adam E, et al. Epidemiological investigation of equine rotavirus B outbreaks in horses in central Kentucky. Vet Microbiol. 2024;298:110278. pmid:39437661

* View Article

* PubMed/NCBI

* Google Scholar

62. 62. Muñoz M, Alvarez M, Lanza I, Cármenes P. Role of enteric pathogens in the aetiology of neonatal diarrhoea in lambs and goat kids in Spain. Epidemiol Infect. 1996;117(1):203–11. pmid:8760970

* View Article

* PubMed/NCBI

* Google Scholar

63. 63. Boey K, Shiokawa K, Avsaroglu H, Rajeev S. Seroprevalence of rodent pathogens in wild rats from the island of St. Kitts, West Indies. Animals (Basel). 2019;9(5):228. pmid:31083284

* View Article

* PubMed/NCBI

* Google Scholar

64. 64. Alam MM, Pun SB, Gauchan P, Yokoo M, Doan YH, Tran TNH, et al. The first identification of rotavirus B from children and adults with acute diarrhoea in Kathmandu, Nepal. Trop Med Health. 2013;41(3):129–34. pmid:24155654

* View Article

* PubMed/NCBI

* Google Scholar

65. 65. Lahon A, Maniya NH, Tambe GU, Chinchole PR, Purwar S, Jacob G, et al. Group B rotavirus infection in patients with acute gastroenteritis from India: 1994–1995 and 2004–2010. Epidemiol Infect. 2013;141(5):969–75. pmid:22813354

* View Article

* PubMed/NCBI

* Google Scholar

66. 66. Joshi MS, Lole KS, Barve US, Salve DS, Ganorkar NN, Chavan NA, et al. Investigation of a large waterborne acute gastroenteritis outbreak caused by group B rotavirus in Maharashtra state, India. J Med Virol. 2019;91(10):1877–81. pmid:31276221

* View Article

* PubMed/NCBI

* Google Scholar

67. 67. Sanekata T, Ahmed MU, Kader A, Taniguchi K, Kobayashi N. Human group B rotavirus infections cause severe diarrhea in children and adults in Bangladesh. J Clin Microbiol. 2003;41(5):2187–90. pmid:12734276

* View Article

* PubMed/NCBI

* Google Scholar

68. 68. Page NA, Netshikweta R, Tate JE, Madhi SA, Parashar UD, Groome MJ, et al. Microorganisms detected in intussusception cases and controls in children <3 years in South Africa from 2013 to 2017. Open Forum Infect Dis. 2023;10(9):ofad458. pmid:37720699

* View Article

* PubMed/NCBI

* Google Scholar

69. 69. Pérez-Losada M, Arenas M, Galán JC, Palero F, González-Candelas F. Recombination in viruses: mechanisms, methods of study, and evolutionary consequences. Infect Genet Evol. 2015;30:296–307. pmid:25541518

* View Article

* PubMed/NCBI

* Google Scholar

70. 70. Chen F, Knutson TP, Porter RE, Ciarlet M, Mor SK, Marthaler DG. Genome characterization of Turkey Rotavirus G strains from the United States identifies potential recombination events with human Rotavirus B strains. J Gen Virol. 2017;98(12):2931–6. pmid:29168675

* View Article

* PubMed/NCBI

* Google Scholar

71. 71. Martella V, Ciarlet M, Pratelli A, Arista S, Terio V, Elia G, et al. Molecular analysis of the VP7, VP4, VP6, NSP4, and NSP5/6 genes of a buffalo rotavirus strain: identification of the rare P[3] rhesus rotavirus-like VP4 gene allele. J Clin Microbiol. 2003;41(12):5665–75. pmid:14662959

* View Article

* PubMed/NCBI

* Google Scholar

72. 72. Matthijnssens J, Rahman M, Martella V, Xuelei Y, De Vos S, De Leener K, et al. Full genomic analysis of human rotavirus strain B4106 and lapine rotavirus strain 30/96 provides evidence for interspecies transmission. J Virol. 2006;80(8):3801–10. pmid:16571797

* View Article

* PubMed/NCBI

* Google Scholar

73. 73. Estes MK. Rotaviruses and their replication. In: Knipe DM, Howley PM, Griffin DE, Lamb RA, Martin MA, Roizman B, Straus SE, editors. Fields virology. 4th ed. Philadelphia, PA: Lippincott/The Williams & Wilkins Co. 2004. p. 1747–86.

74. 74. Hoshino Y, Kapikian AZ. Classification of rotavirus VP4 and VP7 serotypes. Arch Virol Suppl. 1996;12:99–111. pmid:9015107

* View Article

* PubMed/NCBI

* Google Scholar

75. 75. Simon-Loriere E, Holmes EC. Why do RNA viruses recombine?. Nat Rev Microbiol. 2011;9(8):617–26. pmid:21725337

* View Article

* PubMed/NCBI

* Google Scholar

76. 76. Marsh GA, Hatami R, Palese P. Specific residues of the influenza A virus hemagglutinin viral RNA are important for efficient packaging into budding virions. J Virol. 2007;81(18):9727–36. pmid:17634232

* View Article

* PubMed/NCBI

* Google Scholar

77. 77. Essere B, Yver M, Gavazzi C, Terrier O, Isel C, Fournier E, et al. Critical role of segment-specific packaging signals in genetic reassortment of influenza A viruses. Proc Natl Acad Sci U S A. 2013;110(40):E3840-8. pmid:24043788

* View Article

* PubMed/NCBI

* Google Scholar

78. 78. McDonald SM, Patton JT. Assortment and packaging of the segmented rotavirus genome. Trends Microbiol. 2011;19(3):136–44. pmid:21195621

* View Article

* PubMed/NCBI

* Google Scholar

79. 79. Miranda ARM, da Silva Mendes G, Santos N. Rotaviruses A and C in dairy cattle in the state of Rio de Janeiro, Brazil. Braz J Microbiol. 2022;53(3):1657–63. pmid:35478312

* View Article

* PubMed/NCBI

* Google Scholar

80. 80. Dias JBL, Pinheiro MS, Petrucci MP, Travassos CEPF, Mendes GS, Santos N. Rotavirus A and D circulating in commercial chicken flocks in southeastern Brazil. Vet Res Commun. 2024;48(2):743–8. pmid:37878188

* View Article

* PubMed/NCBI

* Google Scholar

81. 81. Shi K, Zhou H, Feng S, He J, Li B, Long F, et al. Development of a quadruplex RT-qPCR for the detection of porcine rotaviruses and the phylogenetic analysis of porcine RVH in China. Pathogens. 2023;12(9):1091. pmid:37764899

* View Article

* PubMed/NCBI

* Google Scholar

Web References

1. . Magyar Madártani és Természetvédelmi Egyesület (MME). [cited on 2024 Mar 25. ]. Available from: https://mme.hu/magyarorszagmadarai/madaradatbazis-larrid

* View Article

* Google Scholar

2. . BirdLife International. [cited on 2024 Mar 25. ]. Available from: http://datazone.birdlife.org/species/factsheet/black-headed-gull-larus-ridibundus

* View Article

* Google Scholar

3. . International Committee on Taxonomy of Viruses: ICTV. [cited on 2024 Mar 25. ]. Available from: https://ictv.global/report/chapter/sedoreoviridae/sedoreoviridae/rotavirus

* View Article

* Google Scholar

4. . International Committee on Taxonomy of Viruses: ICTV. [cited on 2024 Mar 25. ]. Available from: https://ictv.global/report_9th/RNApos/Astroviridae

* View Article

* Google Scholar

5. . Integrated DNA Technologies (IDT DNA). Available from: https://eu.idtdna.com

* View Article

* Google Scholar

6. . National Center for Biotechnology Information, Needleman-Wunsch Global Align. Available from: https://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE_TYPE=BlastSearch&PROG_DEF=blastn&BLAST_PROG_DEF=blastn&BLAST_SPEC=GlobalAln&LINK_LOC=BlastHomeLink

* View Article

* Google Scholar

7. . Multiple Sequence Alignment (MSA) – MUSCLE. Available from: https://www.ebi.ac.uk/jdispatcher/msa/muscle

* View Article

* Google Scholar

8. . European Centre for Disease Prevention and Control (ECDC), Disease factsheet about rotavirus. [cited on 2024 Dec 10. ]. Available from: https://www.ecdc.europa.eu/en/rotavirus-infection/facts

* View Article

* Google Scholar

Citation: Pankovics P, Takáts K, Urbán P, Mátics R, Reuter G, Boros Á (2025) Identification of a potential interspecies reassortant rotavirus G and avastrovirus 2 co-infection from black-headed gull (Chroicocephalus ridibundus) in Hungary. PLoS ONE 20(3): e0317400. https://doi.org/10.1371/journal.pone.0317400

About the Authors:

Péter Pankovics

Roles: Formal analysis, Investigation, Methodology, Visualization, Writing – original draft

E-mail: [email protected]

Affiliation: Department of Medical Microbiology and Immunology, Medical School, University of Pécs, Pécs, Hungary

ORICD: https://orcid.org/0000-0001-9419-2522

Károly Takáts

Roles: Formal analysis

Affiliation: Department of Medical Microbiology and Immunology, Medical School, University of Pécs, Pécs, Hungary

Péter Urbán

Roles: Data curation, Formal analysis

Affiliation: János Szentágothai Research Centre of the University of Pécs, Bioinformatics Research Group, Genomics and Bioinformatics Core Facility, Pécs, Hungary

Róbert Mátics

Roles: Resources

Affiliations: Hungarian Nature Research Society, Ajka, Hungary, Department of Behavioural Science, Medical School, University of Pécs, Pécs, Hungary

Gábor Reuter

Roles: Funding acquisition, Writing – review & editing

Affiliation: Department of Medical Microbiology and Immunology, Medical School, University of Pécs, Pécs, Hungary

Ákos Boros

Roles: Funding acquisition, Methodology, Supervision, Writing – review & editing

Affiliation: Department of Medical Microbiology and Immunology, Medical School, University of Pécs, Pécs, Hungary

[/RAW_REF_TEXT]

References

1. Martín-Vélez V, Montalvo T, Afán I, Sánchez-Márquez A, Aymí R, Figuerola J, et al. Gulls living in cities as overlooked seed dispersers within and outside urban environments. Sci Total Environ. 2022;823:153535. pmid:35104514

2. Wille M, Eden J-S, Shi M, Klaassen M, Hurt AC, Holmes EC. Virus-virus interactions and host ecology are associated with RNA virome structure in wild birds. Mol Ecol. 2018;27(24):5263–78. pmid:30375075

3. Matthijnssens J, Attoui H, Bányai K, Brussaard CPD, Danthi P, Del Vas M, et al. ICTV virus taxonomy profile: sedoreoviridae 2022. J Gen Virol. 2022;103(10):001782. pmid:36215107

4. Cortez V, Meliopoulos VA, Karlsson EA, Hargest V, Johnson C, Schultz-Cherry S. Astrovirus biology and pathogenesis. Annu Rev Virol. 2017;4(1):327–48. pmid:28715976

5. Ogden KM, Johne R, Patton JT. Rotavirus RNA polymerases resolve into two phylogenetically distinct classes that differ in their mechanism of template recognition. Virology. 2012;431(1–2):50–7. pmid:22687427

6. Kindler E, Trojnar E, Heckel G, Otto PH, Johne R. Analysis of rotavirus species diversity and evolution including the newly determined full-length genome sequences of rotavirus F and G. Infect Genet Evol. 2013;14:58–67. pmid:23237956

7. Matthijnssens J, Otto PH, Ciarlet M, Desselberger U, Van Ranst M, Johne R. VP6-sequence-based cutoff values as a criterion for rotavirus species demarcation. Arch Virol. 2012;157(6):1177–82. pmid:22430951

8. Trojnar E, Otto P, Roth B, Reetz J, Johne R. The genome segments of a group D rotavirus possess group A-like conserved termini but encode group-specific proteins. J Virol. 2010;84(19):10254–65. pmid:20631147

9. Johne R, Schilling-Loeffler K, Ulrich RG, Tausch SH. Whole genome sequence analysis of a prototype strain of the novel putative rotavirus species L. Viruses. 2022;14(3):462. pmid:35336869

10. Dhama K, Saminathan M, Karthik K, Tiwari R, Shabbir MZ, Kumar N, et al. Avian rotavirus enteritis – an updated review. Vet Q. 2015;35(3):142–58. pmid:25917772

11. Jere KC, Mlera L, Page NA, van Dijk AA, O’Neill HG. Whole genome analysis of multiple rotavirus strains from a single stool specimen using sequence-independent amplification and 454® pyrosequencing reveals evidence of intergenotype genome segment recombination. Infect Genet Evol. 2011;11(8):2072–82. pmid:22019521

12. Jere KC, Chaguza C, Bar-Zeev N, Lowe J, Peno C, Kumwenda B, et al. Emergence of double- and triple-gene reassortant G1P[8] rotaviruses possessing a DS-1-like backbone after rotavirus vaccine introduction in Malawi. J Virol. 2018;92(3):e01246-17. pmid:29142125

13. Hoxie I, Dennehy JJ. Intragenic recombination influences rotavirus diversity and evolution. Virus Evol. 2020;6(1):vez059. pmid:31949920

14. Suzuki Y, Gojobori T, Nakagomi O. Intragenic recombinations in rotaviruses. FEBS Lett. 1998;427(2):183–7. pmid:9607308

15. Parra GI, Espínola EE, Amarilla AA, Stupka J, Martinez M, Zunini M, et al. Diversity of group A rotavirus strains circulating in Paraguay from 2002 to 2005: detection of an atypical G1 in South America. J Clin Virol. 2007;40(2):135–41. pmid:17720620

16. Phan TG, Okitsu S, Maneekarn N, Ushijima H. Evidence of intragenic recombination in G1 rotavirus VP7 genes. J Virol. 2007;81(18):10188–94. pmid:17609273

17. Matthijnssens J, Rahman M, Van Ranst M. Two out of the 11 genes of an unusual human G6P[6] rotavirus isolate are of bovine origin. J Gen Virol. 2008;89(Pt 10):2630–5. pmid:18796733

18. Suzuki T, Inoue D. Full genome-based genotyping system for rotavirus H and detection of potential gene recombination in nonstructural protein 3 between porcine rotavirus H and rotavirus C. J Gen Virol. 2018;99(12):1582–9. pmid:30328806

19. Krasnikov N, Yuzhakov A. Interspecies recombination in NSP3 gene in the first porcine rotavirus H in Russia identified using nanopore-based metagenomic sequencing. Front Vet Sci. 2023;10:1302531. pmid:38116510

20. Gabbay YB, Leite JPG, Oliveira DS, Nakamura LS, Nunes MRT, Mascarenhas JDP, et al. Molecular epidemiology of astrovirus type 1 in Belém, Brazil, as an agent of infantile gastroenteritis, over a period of 18 years (1982–2000): identification of two possible new lineages. Virus Res. 2007;129(2):166–74. pmid:17714822

21. Finkbeiner SR, Holtz LR, Jiang Y, Rajendran P, Franz CJ, Zhao G, et al. Human stool contains a previously unrecognized diversity of novel astroviruses. Virol J. 2009;6:161. pmid:19814825

22. Tan Z, Zhai H, Sun R, Sun Z, Xie R, Zhang L, et al. Complete genome sequence and phylogenetic analysis of a goose astrovirus isolate in China. Braz J Microbiol. 2023;54(1):427–34. pmid:36327040

23. Yin L, Zhou Q, Mai K, Huang J, Yan Z, Wei X, et al. Isolation and characterization of a novel chicken astrovirus in China. Poult Sci. 2021;100(9):101363. pmid:34352410

24. Donato C, Vijaykrishna D. The broad host range and genetic diversity of mammalian and avian astroviruses. Viruses. 2017;9(5):102. pmid:28489047

25. De Benedictis P, Schultz-Cherry S, Burnham A, Cattoli G. Astrovirus infections in humans and animals - molecular biology, genetic diversity, and interspecies transmissions. Infect Genet Evol. 2011;11(7):1529–44. pmid:21843659

26. Reuter G, Boros Á, Mátics R, Kapusinszky B, Delwart E, Pankovics P. A novel avian-like hepatitis E virus in wild aquatic bird, little egret (Egretta garzetta), in Hungary. Infect Genet Evol. 2016;46:74–7. pmid:27876615

27. Boros Á, Albert M, Urbán P, Herczeg R, Gáspár G, Balázs B, et al. Unusual “Asian-origin” 2c to 2b point mutant canine parvovirus (Parvoviridae) and canine astrovirus (Astroviridae) co-infection detected in vaccinated dogs with an outbreak of severe haemorrhagic gastroenteritis with high mortality rate in Hungary. Vet Res Commun. 2022;46(4):1355–61. pmid:36129562

28. Boros Á, Pankovics P, László Z, Urbán P, Herczeg R, Gáspár G, et al. The genomic and epidemiological investigations of enteric viruses of domestic caprine (Capra hircus) revealed the presence of multiple novel viruses related to known strains of humans and ruminant livestock species. Microbiol Spectr. 2023;11(6):e0253323. pmid:37823638

29. Sverdlov E, Azhikina T. Primer walking. In: eLS, (Ed.). 2005. Available from:

30. Boros Á, Pankovics P, Simmonds P, Reuter G. Novel positive-sense, single-stranded RNA (+ssRNA) virus with di-cistronic genome from intestinal content of freshwater carp (Cyprinus carpio). PLoS One. 2011;6(12):e29145. pmid:22195010

31. Wakuda M, Pongsuwanna Y, Taniguchi K. Complete nucleotide sequences of two RNA segments of human picobirnavirus. J Virol Methods. 2005;126(1–2):165–9. pmid:15847933

32. Reuter G, Boros A, Delwart E, Pankovics P. Novel seadornavirus (family Reoviridae) related to Banna virus in Europe. Arch Virol. 2013;158(10):2163–7. pmid:23624680

33. Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2019;20(4):1160–6. pmid:28968734

34. Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, et al. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol. 2020;37(5):1530–4. pmid:32011700

35. Kalyaanamoorthy S, Minh BQ, Wong TKF, von Haeseler A, Jermiin LS. ModelFinder: fast model selection for accurate phylogenetic estimates. Nat Methods. 2017;14(6):587–9. pmid:28481363

36. Hoang DT, Chernomor O, von Haeseler A, Minh BQ, Vinh LS. UFBoot2: improving the ultrafast bootstrap approximation. Mol Biol Evol. 2018;35(2):518–22. pmid:29077904

37. Lole KS, Bollinger RC, Paranjape RS, Gadkari D, Kulkarni SS, Novak NG, et al. Full-length human immunodeficiency virus type 1 genomes from subtype C-infected seroconverters in India, with evidence of intersubtype recombination. J Virol. 1999;73(1):152–60. pmid:9847317

38. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38(7):3022–7. pmid:33892491

39. McDonald SM, Aguayo D, Gonzalez-Nilo FD, Patton JT. Shared and group-specific features of the rotavirus RNA polymerase reveal potential determinants of gene reassortment restriction. J Virol. 2009;83(12):6135–48. pmid:19357162

40. Vásquez-del Carpió R, Morales JL, Barro M, Ricardo A, Spencer E. Bioinformatic prediction of polymerase elements in the rotavirus VP1 protein. Biol Res. 2006;39(4):649–59. pmid:17657346

41. Ogden KM, Snyder MJ, Dennis AF, Patton JT. Predicted structure and domain organization of rotavirus capping enzyme and innate immune antagonist VP3. J Virol. 2014;88(16):9072–85. pmid:24899176

42. Vasquez-Del Carpio R, Gonzalez-Nilo FD, Riadi G, Taraporewala ZF, Patton JT. Histidine triad-like motif of the rotavirus NSP2 octamer mediates both RTPase and NTPase activities. J Mol Biol. 2006;362(3):539–54. pmid:16934294

43. Pankovics P, Boros Á, László Z, Szekeres S, Földvári G, Altan E, et al. Genome characterization, prevalence and tissue distribution of astrovirus, hepevirus and norovirus among wild and laboratory rats (Rattus norvegicus) and mice (Mus musculus) in Hungary. Infect Genet Evol. 2021;93:104942. pmid:34044191

44. Koonin EV. The phylogeny of RNA-dependent RNA polymerases of positive-strand RNA viruses. J Gen Virol. 1991;72(Pt 9):2197–206. pmid:1895057

45. Méndez E, Murillo A, Velázquez R, Burnham A, Arias CF. Replication cycle of astroviruses. Astrovirus Res Ideas. 2012;7:19–45.

46. Kariithi HM, Volkening JD, Chiwanga GH, Pantin-Jackwood MJ, Msoffe PLM, Suarez DL. Genome sequences and characterization of chicken astrovirus and avian nephritis virus from tanzanian live bird markets. Viruses. 2023;15(6):1247. pmid:37376547

47. McNulty MS, Todd D, Allan GM, McFerran JB, Greene JA. Epidemiology of rotavirus infection in broiler chickens: recognition of four serogroups. Arch Virol. 1984;81(1–2):113–21. pmid:6331344

48. Otto P, Liebler-Tenorio EM, Elschner M, Reetz J, Löhren U, Diller R. Detection of rotaviruses and intestinal lesions in broiler chicks from flocks with runting and stunting syndrome (RSS). Avian Dis. 2006;50(3):411–8. pmid:17039842

49. Kang SY, Nagaraja KV, Newman JA. Primary isolation and identification of avian rotaviruses from turkeys exhibiting signs of clinical enteritis in a continuous MA 104 cell line. Avian Dis. 1986;30(3):494–9. pmid:3021099

50. Bodewes R, van Run PRWA, Schürch AC, Koopmans MPG, Osterhaus ADME, Baumgärtner W, et al. Virus characterization and discovery in formalin-fixed paraffin-embedded tissues. J Virol Methods. 2015;214:54–9. pmid:25681526

51. Hou G, Zeng Q, Matthijnssens J, Greenberg HB, Ding S. Rotavirus NSP1 contributes to intestinal viral replication, pathogenesis, and transmission. mBio. 2021;12(6):e0320821. pmid:34903043

52. Pham T, Perry JL, Dosey TL, Delcour AH, Hyser JM. The rotavirus NSP4 viroporin domain is a calcium-conducting ion channel. Sci Rep. 2017;7:43487. pmid:28256607

53. Phan MVT, Anh PH, Cuong NV, Munnink BBO, van der Hoek L, My PT, et al. Unbiased whole-genome deep sequencing of human and porcine stool samples reveals circulation of multiple groups of rotaviruses and a putative zoonotic infection. Virus Evol. 2016;2(2):vew027. pmid:28748110

54. Ahmed MU, Kobayashi N, Wakuda M, Sanekata T, Taniguchi K, Kader A, et al. Genetic analysis of group B human rotaviruses detected in Bangladesh in 2000 and 2001. J Med Virol. 2004;72(1):149–55. pmid:14635024

55. Hung T, Chen GM, Wang CG, Chou ZY, Chao TX, Ye WW, et al. Rotavirus-like agent in adult non-bacterial diarrhoea in China. Lancet. 1983;2(8358):1078–9. pmid:6138617

56. Smitalova R, Rodak L, Smid B, Psikal I. Detection of nongroup A rotaviruses in faecal samples of pigs in the Czech Republic. Vet Med (Praha). 2009;54(1):12–8.

57. Vidal A, Martín-Valls GE, Tello M, Mateu E, Martín M, Darwich L. Prevalence of enteric pathogens in diarrheic and non-diarrheic samples from pig farms with neonatal diarrhea in the North East of Spain. Vet Microbiol. 2019;237:108419. pmid:31585655

58. Otto PH, Rosenhain S, Elschner MC, Hotzel H, Machnowska P, Trojnar E, et al. Detection of rotavirus species A, B and C in domestic mammalian animals with diarrhoea and genotyping of bovine species A rotavirus strains. Vet Microbiol. 2015;179(3–4):168–76. pmid:26223422

59. Chen F, Knutson TP, Ciarlet M, Sturos M, Marthaler DG. Complete genome characterization of a rotavirus B (RVB) strain identified in Alpine goat kids with enteritis reveals inter-species transmission with RVB bovine strains. J Gen Virol. 2018;99(4):457–63. pmid:29517476

60. Uprety T, Sreenivasan CC, Hause BM, Li G, Odemuyiwa SO, Locke S, et al. Identification of a ruminant origin group B rotavirus associated with diarrhea outbreaks in foals. Viruses. 2021;13(7):1330. pmid:34372536

61. Sreenivasan CC, Naveed A, Uprety T, Soni S, Jacob O, Adam E, et al. Epidemiological investigation of equine rotavirus B outbreaks in horses in central Kentucky. Vet Microbiol. 2024;298:110278. pmid:39437661

62. Muñoz M, Alvarez M, Lanza I, Cármenes P. Role of enteric pathogens in the aetiology of neonatal diarrhoea in lambs and goat kids in Spain. Epidemiol Infect. 1996;117(1):203–11. pmid:8760970

63. Boey K, Shiokawa K, Avsaroglu H, Rajeev S. Seroprevalence of rodent pathogens in wild rats from the island of St. Kitts, West Indies. Animals (Basel). 2019;9(5):228. pmid:31083284

64. Alam MM, Pun SB, Gauchan P, Yokoo M, Doan YH, Tran TNH, et al. The first identification of rotavirus B from children and adults with acute diarrhoea in Kathmandu, Nepal. Trop Med Health. 2013;41(3):129–34. pmid:24155654

65. Lahon A, Maniya NH, Tambe GU, Chinchole PR, Purwar S, Jacob G, et al. Group B rotavirus infection in patients with acute gastroenteritis from India: 1994–1995 and 2004–2010. Epidemiol Infect. 2013;141(5):969–75. pmid:22813354

66. Joshi MS, Lole KS, Barve US, Salve DS, Ganorkar NN, Chavan NA, et al. Investigation of a large waterborne acute gastroenteritis outbreak caused by group B rotavirus in Maharashtra state, India. J Med Virol. 2019;91(10):1877–81. pmid:31276221

67. Sanekata T, Ahmed MU, Kader A, Taniguchi K, Kobayashi N. Human group B rotavirus infections cause severe diarrhea in children and adults in Bangladesh. J Clin Microbiol. 2003;41(5):2187–90. pmid:12734276

68. Page NA, Netshikweta R, Tate JE, Madhi SA, Parashar UD, Groome MJ, et al. Microorganisms detected in intussusception cases and controls in children <3 years in South Africa from 2013 to 2017. Open Forum Infect Dis. 2023;10(9):ofad458. pmid:37720699

69. Pérez-Losada M, Arenas M, Galán JC, Palero F, González-Candelas F. Recombination in viruses: mechanisms, methods of study, and evolutionary consequences. Infect Genet Evol. 2015;30:296–307. pmid:25541518

70. Chen F, Knutson TP, Porter RE, Ciarlet M, Mor SK, Marthaler DG. Genome characterization of Turkey Rotavirus G strains from the United States identifies potential recombination events with human Rotavirus B strains. J Gen Virol. 2017;98(12):2931–6. pmid:29168675

71. Martella V, Ciarlet M, Pratelli A, Arista S, Terio V, Elia G, et al. Molecular analysis of the VP7, VP4, VP6, NSP4, and NSP5/6 genes of a buffalo rotavirus strain: identification of the rare P[3] rhesus rotavirus-like VP4 gene allele. J Clin Microbiol. 2003;41(12):5665–75. pmid:14662959

72. Matthijnssens J, Rahman M, Martella V, Xuelei Y, De Vos S, De Leener K, et al. Full genomic analysis of human rotavirus strain B4106 and lapine rotavirus strain 30/96 provides evidence for interspecies transmission. J Virol. 2006;80(8):3801–10. pmid:16571797

73. Estes MK. Rotaviruses and their replication. In: Knipe DM, Howley PM, Griffin DE, Lamb RA, Martin MA, Roizman B, Straus SE, editors. Fields virology. 4th ed. Philadelphia, PA: Lippincott/The Williams & Wilkins Co. 2004. p. 1747–86.

74. Hoshino Y, Kapikian AZ. Classification of rotavirus VP4 and VP7 serotypes. Arch Virol Suppl. 1996;12:99–111. pmid:9015107

75. Simon-Loriere E, Holmes EC. Why do RNA viruses recombine?. Nat Rev Microbiol. 2011;9(8):617–26. pmid:21725337

76. Marsh GA, Hatami R, Palese P. Specific residues of the influenza A virus hemagglutinin viral RNA are important for efficient packaging into budding virions. J Virol. 2007;81(18):9727–36. pmid:17634232

77. Essere B, Yver M, Gavazzi C, Terrier O, Isel C, Fournier E, et al. Critical role of segment-specific packaging signals in genetic reassortment of influenza A viruses. Proc Natl Acad Sci U S A. 2013;110(40):E3840-8. pmid:24043788

78. McDonald SM, Patton JT. Assortment and packaging of the segmented rotavirus genome. Trends Microbiol. 2011;19(3):136–44. pmid:21195621

79. Miranda ARM, da Silva Mendes G, Santos N. Rotaviruses A and C in dairy cattle in the state of Rio de Janeiro, Brazil. Braz J Microbiol. 2022;53(3):1657–63. pmid:35478312

80. Dias JBL, Pinheiro MS, Petrucci MP, Travassos CEPF, Mendes GS, Santos N. Rotavirus A and D circulating in commercial chicken flocks in southeastern Brazil. Vet Res Commun. 2024;48(2):743–8. pmid:37878188

81. Shi K, Zhou H, Feng S, He J, Li B, Long F, et al. Development of a quadruplex RT-qPCR for the detection of porcine rotaviruses and the phylogenetic analysis of porcine RVH in China. Pathogens. 2023;12(9):1091. pmid:37764899