Metagenomic Identification of Novel Eukaryotic Viruses with Small DNA Genomes in Pheasants.

A panel of intestinal samples collected from common pheasants (Phasianus colchicus) between 2008 and 2017 was used for metagenomic investigation using an unbiased enrichment protocol and different bioinformatic pipelines. The number of sequence reads in the metagenomic analysis ranged from 1,419,265 to 17,507,704 with a viral sequence read rate ranging from 0.01% to 59%. When considering the sequence reads of eukaryotic viruses, RNA and DNA viruses were identified in the samples, including but not limited to coronaviruses, reoviruses, parvoviruses, and CRESS DNA viruses (i.e., circular Rep-encoding single-stranded DNA viruses). Partial or nearly complete genome sequences were reconstructed of at least three different parvoviruses (dependoparvovirus, aveparvovirus and chaphamaparvovirus), as well as gyroviruses and diverse CRESS DNA viruses. Generating information of virus diversity will serve as a basis for developing specific diagnostic tools and for structured epidemiological investigations, useful to assess the impact of these novel viruses on animal health.

Spatiotemporal Distribution of PRRSV-1 Clades in Hungary with a Focus on the Era of Disease Eradication.

Porcine reproductive and respiratory syndrome (PRRS) is the cause of the most severe economic losses in the pig industry worldwide. PRRSV is extremely diverse in Europe, which poses a significant challenge to disease control within a country or any region. With the combination of phylogenetic reconstruction and network analysis, we aimed to uncover the major routes of the dispersal of PRRSV clades within Hungary. In brief, by analyzing >2600 ORF5 sequences, we identified at least 12 clades (including 6 clades within lineage 1 and 3 clades within lineage 3) common in parts of Western Europe (including Denmark, Germany and the Netherlands) and identified 2 novel clades (designated X1 and X2). Of interest, some genetic clades unique to other central European countries, such as the Czech Republic and Poland, were not identified. The pattern of PRRSV clade distribution is consistent with the route of the pig trade among countries, showing that most of the identified clades were introduced from Western Europe when fatteners were transported to Hungary. As a result of rigorous implementation of the national eradication program, the swine population was declared officially free from PRRSV. This map of viral diversity and clade distribution will serve as valuable baseline information for the maintenance of PRRSV-free status in the post-eradication era.

A Snapshot on the Genomic Epidemiology of Turkey Reovirus Infections, Hungary.

Reovirus infections in turkeys are associated with arthritis and lameness. Viral genome sequence data are scarce, which makes an accurate description of the viral evolution and epidemiology difficult. In this study, we isolated and characterized turkey reoviruses from Hungary. The isolates were identified in 2016; these isolates were compared with earlier Hungarian turkey reovirus strains and turkey reoviruses isolated in the 2010s in the United States. Gene-wise sequence and phylogenetic analyses identified the cell-receptor binding protein and the main neutralization antigen, σC, to be the most conserved. The most genetically diverse gene was another surface antigen coding gene, µB. This gene was shown to undergo frequent reassortment among chicken and turkey origin reoviruses. Additional reassortment events were found primarily within members of the homologous turkey reovirus clade. Our data showed evidence for low variability among strains isolated from independent outbreaks, a finding that suggests a common source of turkey reoviruses in Hungarian turkey flocks. Given that commercial vaccines are not available, identification of the source of these founder virus strains would permit a more efficient prevention of disease outbreaks before young birds are settled to fattening facilities.

Genetic diversity among reptilian orthoreoviruses isolated from pet snakes and lizards.

Reovirus infections in reptiles are frequently detected and associated with various clinical diseases; yet, our knowledge about their genetic diversity and evolutionary relationships remains limited. In this study, we characterize at the genomic level five reptile origin orthoreovirus strains isolated from exotic snakes and lizards in Hungary and Germany. The genomic organization of the study strains was similar to that of the representative strains of reptile origin reoviruses belonging to species Reptilian orthoreovirus and Testudine orthoreovirus. Additionally, all five study strains clustered with the bush viper origin reference Reptilian orthoreovirus strain, 47/02. The nucleotide sequence divergence among strains fell from 56.64 to 99.36%. Based on genome segment constellations two well separated groups were observed, which may represent two genetic lineages of reptilian orthoreoviruses we tentatively referred here as genogroups, classifying two squamata origin strains with available whole genome sequences into genogroup I (GGI) and four strains into genogroup II (GGII). The representative GGI and GGII Reptilian orthoreovirus strains are characterized by moderate-to-high nucleotide and amino acid similarities within genogroups (range, 69.45 to 99.36% and 74.64 to 100.00%), whereas lower nucleotide and amino acid similarities (range, 56.64 to 77.24% and 54.53 to 93.85%) and different structures of the bicistronic S1 segment were found between genogroups. Further studies are needed to explore the genomic diversity among reptilian reoviruses of squamata origin; this would be critical to establish a robust classification system for these viruses and to see if interaction among members of distinct lineages may result in viable progenies with novel genetic features.

Corrigendum: Genetic diversity of imported PRRSV-2 strains, 2005-2020, Hungary.

[This corrects the article DOI: 10.3389/fvets.2022.986850.].

Genetic Diversity of Actinobacillus pleuropneumoniae Serovars in Hungary.

A total of 114 Actinobacillus pleuropneumoniae isolates from porcine hemorrhagic necrotic pleuropneumonia were characterized by the examination of biotype, serovar, antibiotic resistance genes, and genes of toxin production. Pulsed-field gel electrophoresis was used to analyze their genetic relationship, which identified 16 clusters. Serovar 2 (50 isolates), serovar 13 (25 isolates), serovar 9 (11 isolates), and serovar 16 (7 isolates) were the most frequent serovars. Serovar 2 formed nine distinguishable clusters; serovar 13 and serovar 16 were less diverse, exhibiting two potentially related subclusters; serovar 9 was represented by a single cluster. Remarkably small differences were seen in the core genome when nine representative isolates of serovar 13 were subjected to whole-genome sequencing. Tetracycline resistance was relatively frequent in the two clusters of serovar 13; one of them was also frequently resistant against beta-lactams. Resistance in other serovars was sporadic. All isolates carried the apxIV gene. The toxin profiles of serovar 2 were characterized by the production of ApxII and ApxIII toxins, except for a small cluster of three isolates: serovar 9 and serovar 16 isolates produced ApxI and ApxII toxins. Serovar 13 carried apxII and apxIBD genes, indicating the production of the ApxII toxin, but not of ApxI or ApxIII. The unusually high frequency and low diversity of serovar 13 are not explained by its virulence properties, but the high frequency of resistance to beta-lactams and tetracyclines may have played a role in its spread. The emergence of serovar 16 may be facilitated by its high virulence, also explaining its high clonality.

Genetic diversity of imported PRRSV-2 strains, 2005-2020, Hungary.

Porcine reproductive and respiratory syndrome virus 2 (PRRSV-2) remains sporadic in Europe. In this study, we investigated the molecular epidemiology of PRRSV-2 infections encompassing 15 years in Hungary. Partial (423 bp long) ORF5 sequences (n = 44) from 20 Hungarian pig herds were analyzed. The study strains fell into two genetic lineages, L1 and L5, being L5 strains more prevalent (88.6 vs. 11.4%). Pairwise sequence identities within Hungarian representative PRRSV-2 strains ranged between 84.7 to 100% (nucleotide, nt) and 85 to 100% (amino acid, aa). When compared with reference strains, identity values fell between 87 and 100% (L1, nt 87-91%, aa 87-93%, reference strain IAF-exp91; L5, nt 87-100%, aa 88-100%, reference strain Ingelvac MLV). Epidemiologic examination implied that the majority of L5 strains were imported repeatedly from other European countries where Ingelvac MLV was approved for routine use. The emergence of L1 strains was thought to be associated with a single introduction and subsequent dissemination between pig farms of a large integrator. Results presented here contribute to a better understanding of the epizootiology of PRRSV-2 infections and shed light on the genetic diversity of viral strains in non-endemic countries.

Candidate 'Avian orthoreovirus B': An emerging waterfowl pathogen in Europe and Asia?

A fusogenic virus was isolated from a flock of breeder Pekin ducks in 2019, Hungary. The affected flock experienced a marked decrease in egg production. Histopathological lesions were seen in the oviduct and in the lungs of birds sent for diagnostic investigation. The fusogenic agent was characterized as an orthoreovirus by viral metagenomics. The assembled viral genome was composed of 10 genomic segments and was 23,433 nucleotides (nt) in length. The study strain, designated Reo/HUN/DuckDV/2019, shared low-to-medium gene-wise sequence identity with avian orthoreovirus strains from galliform and anseriform birds (nt, 38.90%-72.33%) as well as with representative strains of neoavian orthoreoviruses (nt, 40.07%-68.23%). On the contrary, the study strain shared 86.48%-95.01% pairwise nt sequence identities with recent German and Chinese reovirus isolates, D2533/6 and Ych, respectively. Phylogenetic analysis clustered all three unusual waterfowl pathogens on a monophyletic branch, indicating a common evolutionary origin of Reo/HUN/DuckDV/2019 with these enigmatic orthoreoviruses described over the past few years. The finding that a candidate new orthoreovirus species, tentatively called Avian orthoreovirus B, was isolated in recent years in Europe and Asia in moribund ducks seems an alarming sign that needs to be better evaluated by extending laboratory diagnosis of viral pathogens in countries where the waterfowl industry is important.

Novel Lineage of Infectious Bronchitis Virus from Sub-Saharan Africa Identified by Random Amplification and Next-Generation Sequencing of Viral Genome.

Avian infectious bronchitis (IB) is among the major viral respiratory and reproductive diseases of chickens caused by Avian coronavirus. In the African continent, IB was first described in countries located in the Mediterranean basin. In other parts of the continent, the epidemiological situation of IB remains unclear. In this study, the complete genome sequences of five IBV strains, originating from the sub-Saharan area were determined. Phylogenetic analysis based on the full-length S1 sequences identified three lineages (GI-14, GI-16, and GI-19) common in Africa and revealed that a strain, D2334/11/2/13/CI, isolated in Ivory Coast may represent a novel lineage within genotype GI. The maximum inter- and intragenotype sequence identities between this strain and other IBVs were 67.58% and 78.84% (nucleotide) and 64.44% and 78.6% (amino acid), respectively. The whole-genome nucleotide identity of the novel variant shared the highest values with a reference Belgian nephropathogenic strain (B1648, 92.4%) and with another study strain from Ivory Coast (D2334/12/2/13/CI, 94.6%). This study illustrates the importance of epidemiological monitoring of IBV in sub-Saharan Africa, as the area may serve as a focal point for newly emerging viral lineages.

Genetic characterisation of a novel reptarenavirus detected in a dead pet red-tailed boa (Boa constrictor).

Boid inclusion body disease (BIBD) is a severe and transmissible disease of snakes worldwide. Reptarenaviruses have been identified as the aetiological agents of BIBD. We determined the almost complete genome sequence of an arenavirus detected in a female red-tailed boa that had succumbed in a private collection in Hungary. We used a combination of next generation sequencing and Sanger sequencing methods. Based on the analysis of the obtained sequence data, the virus, tentatively named Coldvalley virus, seemed to belong to the Reptarenavirus genus of the Arenaviridae family. This classification was confirmed by the genome structure (bisegmented single-stranded RNA) characteristic of the genera Mammarenavirus and Reptarenavirus. The pairwise comparison of the nucleotide and amino acid sequences, as well as the topology of the maximum likelihood phylogenetic trees, suggested that the newly-characterised Coldvalley virus can be classified into the species Rotterdam reptarenavirus.

Whole-genome analysis of a rare G15P[21] group A rotavirus detected at a dairy cattle farm.

Group A rotaviruses (RVAs) are a major cause of severe enteritis in humans and animals. RVAs have been identified in several animal species and their genetic diversity, the segmented nature of their RNA genome and the ability to spill over from one species to another can generate new RVA strains. In this study, we investigated the genome constellations of an unusual, rare, bovine RVA strain, G15P[21], identified from a farm with neonatal diarrhoea of calves in 2006. In parallel, the genome constellations of other RVA strains with different G/P types identified from the same farm in the same time span (2006-2008) were analysed. The genome constellation of strain K53 was G15-P[21]-I2-R2-C2-M2-A13-N2-T9-E2-H3 and was similar, overall, to that of the other bovine RVA strains (G6/10-P[11]-I2-R2-C2-M2-A13-N2-T6-E2-H3) with the exception of the NSP3 segment (T9 vs T6). This study describes RVA genomes with different genotype combinations isolated at a farm and also contributes to the understanding of the diversity and evaluation of rotavirus in a global context.

Genomic Diversity of CRESS DNA Viruses in the Eukaryotic Virome of Swine Feces.

Replication-associated protein (Rep)-encoding single-stranded DNA (CRESS DNA) viruses are a diverse group of viruses, and their persistence in the environment has been studied for over a decade. However, the persistence of CRESS DNA viruses in herds of domestic animals has, in some cases, serious economic consequence. In this study, we describe the diversity of CRESS DNA viruses identified during the metagenomics analysis of fecal samples collected from a single swine herd with apparently healthy animals. A total of nine genome sequences were assembled and classified into two different groups (CRESSV1 and CRESSV2) of the Cirlivirales order (Cressdnaviricota phylum). The novel CRESS DNA viral sequences shared 85.8-96.8% and 38.1-94.3% amino acid sequence identities for the Rep and putative capsid protein sequences compared to their respective counterparts with extant GenBank record. Data presented here show evidence for simultaneous infection of swine herds with multiple novel CRESS DNA viruses, including po-circo-like viruses and fur seal feces-associated circular DNA viruses. Given that viral genomes with similar sequence and structure have been detected in swine fecal viromes from independent studies, investigation of the association between presence of CRESS DNA viruses and swine health conditions seems to be justified.

Recombination Events Shape the Genomic Evolution of Infectious Bronchitis Virus in Europe.

Infectious bronchitis of chicken is a high morbidity and mortality viral disease affecting the poultry industry worldwide; therefore, a better understanding of this pathogen is of utmost importance. The primary aim of this study was to obtain a deeper insight into the genomic diversity of field infectious bronchitis virus (IBV) strains using phylogenetic and recombination analysis. We sequenced the genome of 20 randomly selected strains from seven European countries. After sequencing, we created a genome sequence data set that contained 36 European origin field isolates and 33 vaccine strains. When analyzing these 69 IBV genome sequences, we identified 215 recombination events highlighting that some strains had multiple recombination breaking points. Recombination hot spots were identified mostly in the regions coding for non-structural proteins, and multiple recombination hot spots were identified in the nsp2, nsp3, nsp8, and nsp12 coding regions. Recombination occurred among different IBV genotypes and involved both field and vaccine IBV strains. Ninety percent of field strains and nearly half of vaccine strains showed evidence of recombination. Despite the low number and the scattered geographical and temporal origin of whole-genome sequence data collected from European Gammacoronaviruses, this study underlines the importance of recombination as a major evolutionary mechanism of IBVs.

Genomic characterization of avian and neoavian orthoreoviruses detected in pheasants.

Avian reoviruses are well-known pathogens seriously affecting the productivity of poultry industry. Game birds represent a small segment of the agricultural sector and much remained to be learnt about factors affecting productivity. Here we show that reovirus infections might occur in pheasants and demonstrate that reoviruses of pheasants are of diverse origin. The complete or coding-complete genomic sequences of two Hungarian reovirus strains, D1996/2/1 and Reo/HUN/Pheasant/216/2015, have been determined in this study. The strain D1996/2/1 was isolated in 2012 from birds with gizzard erosion, whereas the other strain was isolated in 2015 from diarrheic pheasant poults. Phylogenetic analyses showed that none of the Hungarian isolates shared common origin with a pheasant reovirus detected recently in the United States. Additionally, we found that Reo/HUN/Pheasant/216/2015 is a multi-reassortant reovirus within the species Avian orthoreovirus that shared genetic relationship with turkey reoviruses (σC), partridge reoviruses (λA, σB), and chicken reoviruses (λB, λC, µA, σA, and σNS), in the respective gene phylogenies, whereas two genes (µB and µNS) did not reveal any possible common ancestors. The other isolate, D1996/2/1, was found to be distantly related to previously described reoviruses raising the possibility that it might represent a novel orthoreovirus species or a new genogroup within the newly accepted species, Neoavian orthoreovirus. The genetic diversity among pheasant reoviruses could raise challenges for virus classification as well as for development of molecular diagnostic tools and vaccine based prevention and control measures.

Genomic evolution of avian polyomaviruses with a focus on budgerigar fledgling disease virus.

Gammapolyomaviruses may cause serious inflammatory diseases in a broad range of avian hosts. In this study we investigated genomic evolution of and selection constraint acting on avian polyomaviruses (APyVs). Our analyses suggested that goose haemorrhagic polyomavirus (GHPV) evolves more slowly (3.03 × 10-5 s/s/y mean evolutionary rate) than budgerigar fledgling disease virus (BFDV), finch polyomavirus (FPyV) and canary polyomavirus (CaPyV) (1.39 × 10-4 s/s/y, 2.63 × 10-4 s/s/y and 1.41 × 10-4 s/s/y mean evolutionary rate, respectively). In general, purifying selection seems to act on the protein coding regions of APyV genomes, although positive Darwinian selection was also predicted in a few positions (e.g., in the large tumor antigen coding region of BFDV and GHPV and in the capsid protein sequences of BFDV). The importance of these aa changes remains elusive. Overall, a better understanding of adaptive changes in the genome of APyVs requires additional data from various incidental hosts and reservoir species.

Phylogeny of Hungarian EBLV-1 strains using whole-genome sequence data.

European bat lyssavirus 1 (EBLV-1) is a widespread lyssavirus across Europe, whose epizootic cycle is linked to a few bat species. Occasionally, EBLV-1 infection may occur in domestic animals and humans. EBLV-1 can be classified into two subtypes, where subtype EBLV-1a shows a wide geographic distribution between France and Russia whereas subtype EBLV-1b is distributed between Spain and Poland. In this study, we determined the genome sequence of two recent EBLV-1a strains detected in Hungary and analysed their adaptive evolution and phylodynamics. The data set that included 100 EBLV-1 genome sequences identified positive selection at selected sites in genes coding for viral proteins (N, codon 18; P, 141 and 155; G, 244 and 488; L, 168, 980, 1597 and 1754). A major genetic clade containing EBLV-1a isolates from Hungary, Slovakia, Denmark and Poland was estimated to have diverged during the 19th century whereas the divergence of the most recent ancestor of Hungarian and Slovakian isolates dates back to 1950 (time span, 1930 to 1970). Phylogeographic analysis of the EBLV-1a genomic sequences demonstrated strong evidence of viral dispersal from Poland to Hungary. This new information indicates that additional migratory flyways may help the virus spread, a finding that supplements the general theory on a west-to-east dispersal of EBLV-1a strains. Long-distance migrant bats may mediate the dispersal of EBLV-1 strains across Europe; however, structured surveillance and extended genome sequencing would be needed to better understand the epizootiology of EBLV-1 infections in Europe.

Mutations potentially associated with decreased susceptibility to fluoroquinolones, macrolides and lincomycin in Mycoplasma synoviae.

Mycoplasma synoviae is one of the economically most significant avian Mycoplasma species. It can cause great financial losses to the poultry industry by inducing respiratory diseases, infectious synovitis, or eggshell apex abnormalities. There are different approaches to control M. synoviae infection. Although antimicrobial therapy cannot replace long-term solutions, like eradication and vaccination, this strategy can be effective in the short term, as adequate antibiotic treatment can relieve economic losses through the attenuation of clinical signs and reduction of transmission. Using broth microdilution method, minimal inhibitory concentration (MIC) values to fourteen antibiotics related to eight antimicrobial groups were determined in 96 M. synoviae strains. Whole genome sequencing and sequence analysis revealed mutations potentially associated with decreased susceptibility to fluoroquinolones, macrolides and lincomycin. Molecular markers responsible for the high MICs to fluoroquinolones were found in the gyrA, gyrB, parC and parE genes. Besides, single nucleotide polymorphisms identified in genes encoding the 23S rRNA were found to be responsible for high MICs to the 50S inhibitor macrolides and lincomycin, while amino acid change in the 50S ribosomal protein L22 could be associated with decreased susceptibility to macrolides. The revealed mutations can contribute to the extension of knowledge about the genetic background of antibiotic resistance in M. synoviae. Moreover, the explored potentially resistance-related mutations may serve as targets for molecular biological assays providing data of antibiotic susceptibility prior to the laborious and time-consuming isolation of M. synoviae strains.

Publisher Correction: Occurrence and genetic diversity of CRESS DNA viruses in wild birds: a Hungarian study.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

The core genome multi-locus sequence typing of Mycoplasma anserisalpingitidis.

BACKGROUND: Mycoplasma anserisalpingitidis is a waterfowl pathogen that mainly infects geese, can cause significant economic losses and is present worldwide. With the advance of whole genome sequencing technologies, new methods are available for the researchers; one emerging methodology is the core genome Multi-Locus Sequence Typing (cgMLST). The core genome contains a high percentage of the coding DNA sequence (CDS) set of the studied strains. The cgMLST schemas are powerful genotyping tools allowing for the investigation of potential epidemics, and precise and reliable classification of the strains. Although whole genome sequences of M. anserisalpingitidis strains are available, to date, no cgMLST schema has been published for this species. RESULTS: In this study, Illumina short reads of 81 M. anserisalpingitidis strains were used, including samples from Hungary, Poland, Sweden, and China. Draft genomes were assembled with the SPAdes software and analysed with the online available chewBBACA program. User made modifications in the program enabled analysis of mycoplasmas and provided similar results as the conventional SeqSphere+ software. The threshold of the presence of CDS in the strains was set to 93% due to the quality of the draft genomes, resulting in the most accurate and robust schema. Three hundred thirty-one CDSs constituted our cgMLST schema (representing 42,77% of the whole CDS set of M. anserisalpingitidis ATCC BAA-2147), and a Neighbor joining tree was created using the allelic profiles. The correlation was observed between the strains' cgMLST profile and geographical origin; however, strains from the same integration but different locations also showed close relationship. Strains isolated from different tissue samples of the same animal revealed highly similar cgMLST profiles. CONCLUSIONS: The Neighbor joining tree from the cgMLST schema closely resembled the real-life spatial and temporal relationships of the strains. The incongruences between background data and the cgMLST profile in the strains from the same integration can be because of the higher probability of contacts between the flocks. This schema can help with the epidemiological investigation and can be used as a basis for further studies.