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.

Genomic Epidemiology and Evolution of Fowl Adenovirus 1.

Fowl adenovirus 1 (FAdV-1) is the main cause of gizzard erosion in chickens. Whole genome sequencing and sequence analyses of 32 FAdV-1 strains from a global collection provided evidence that multiple recombination events have occurred along the entire genome. In gene-wise phylogenies, only the adenoviral pol gene formed a tree topology that corresponded to whole genome-based phylogeny. Virus genetic features that were clearly connected to gizzard erosion were not identified in our analyses. However, some genome variants tended to be more frequently identified from birds with gizzard erosion and strains isolated from healthy birds or birds with non-specific pathologies tended to form common clusters in multiple gene phylogenies. Our data show that the genetic diversity is greater, and the evolutionary mechanisms are more complex within FAdV-1 than previously thought. The implications of these findings for viral pathogenesis and epidemiology await further investigation.

Marked Genotype Diversity among Reoviruses Isolated from Chicken in Selected East-Central European Countries.

The concern that the vaccines currently used against Avian orthoreovirus (ARV) infections are less efficient in the field justifies the need for the close monitoring of circulating ARV strains. In this study, we collected necropsy samples from various chicken breeds and tested for ARV by virus isolation, RT-PCR assay and sequence analysis. ARVs were isolated from birds showing runting-stunting syndrome, uneven growth, lameness or increased mortality, with relative detection rates of 38%, 35%, 6% and 25%, respectively. Partial σC gene sequences were determined for nearly 90% of ARV isolates. The isolates could be classified into one of the major genetic clusters. Interestingly, cluster 2 and cluster 5 were isolated from vaccinated broiler breeders, while clusters 1 to 4 were isolated from unvaccinated broilers. The isolates shared less than 75% amino acid identities with the vaccine strains (range, 44.3-74.6%). This study reaffirms the global distribution of the major genetic clusters of ARVs in chicken. The diversity of ARV strains isolated from unvaccinated broilers was greater than those detected from vaccinated animals, however, the relative importance of passive and active immunity on the selection of novel strains in different chicken breeds needs to be better understood.

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.

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.

A novel gyrovirus in a common pheasant (Phasianus colchicus) with poult enteritis and mortality syndrome.

A novel gyrovirus was detected in an intestinal specimen of a common pheasant that died due to poult enteritis and mortality syndrome. The genome of the pheasant-associated gyrovirus (PAGyV) is 2353 nucleotides (nt) long and contains putative genes for the VP1, VP2, and VP3 proteins in an arrangement that is typical for gyroviruses. Gyrovirus-specific motifs were identified in both the coding region and the intergenic region of the PAGyV genome. The VP1 of PAGyV shares up to 67.6% pairwise nt sequence identity with reference sequences and forms a distinct branch in the phylogenetic tree. Thus, according to the recently described species demarcation criteria, PAGyV belongs to a novel species in the genus Gyrovirus, family Anelloviridae, for which we propose the name "Gyrovirus phaco 1".

Mitochondrial sequence diversity reveals the hybrid origin of invasive gibel carp (Carassius gibelio) populations in Hungary.

BACKGROUND: Invasive gibel carp, Carassius gibelio (Bloch, 1782) has become well-established in the Hungarian waters and now are spreading in the European waters. On major concern now is the potential hybridization between gibel carp and the other invasive species in the Carassius auratus complex (CAC), which may further accelerate the spread of the whole invasive species complex. The identification of gibel carp and their hybrids is difficult because of its morphological similarity to the other species in CAC. Here we carry out a genomic assessment to understand the history of gibel carp invasion and its phylogenetic relationship with the other species in CAC. Three loci of the mitochondrial genome (D-loop, CoI, Cytb) were used to determine the phylogenetic origin of individuals and relarionship among six gibel carp populations and the other species in the CAC. METHODOLGY: A total of 132 gibel carp samples from six locations in Southern Transdanubia (Hungary) were collected after phenotypic identification to measure the genetic diversity within and among gibel carp populations of Southern Transdanubia (Hungary). The genetic background was examined by the sequences of the mitochondrial genome: D-loop, Cytochrome c oxidase I (CoI) and Cytochrome b (Cytb). Mitochondrial genetic markers are excellent tools for phylogenetic studies because they are maternally inherited. Successfully identified haplotypes were aligned and with reference sequences in nucleotide databases (i.e., NCBI-BLAST: National Centre for Biotechnology Information and BOLD: Barcode of Life Data System). The phylogenetic relationships among gibel carp populations were then analyzed together with the reference sequences to understand the relationship and the level of hybridization with the species in CAC. RESULTS: Among the 132 aligned D-loop sequences 22 haplotypes were identified. Further examination of representative individuals of the 22 haplotypes, six Cytb and four CoI sequences were detected. The largest number of haplotypes of all three loci were found in Lake Balaton, the largest shallow lake in Central Europe. Based on the NCBI-BLAST alignment of the D-loop, haplotypes of Carassius auratus auratus and Carassius a. buergeri in CAC were identified in the C. gibelio samples. Further analysis of haplotypes with the other two mitochondrial markers confirmed the occurrence of intragenus hybridization of C. gibelio in the Hungarian waters. CONCLUSION: By using three mitochondrial markers (D-loop, Cytb, CoI), we genomically characterized a gibel carp-complex in Hungarian waters and assessed the C. gibelio phylogenetic status between them. Hybrid origin of locally invasive Carassius taxon was detected in Hungary. It points out that invasive species are not only present in Hungary but reproduce with each other in the waters, further accelerating their spread.

Genomic Epidemiology and Evolution of Duck Hepatitis A Virus.

Duck hepatitis A virus (DHAV), an avian picornavirus, causes high-mortality acute disease in ducklings. Among the three serotypes, DHAV-1 is globally distributed, whereas DHAV-2 and DHAV-3 serotypes are chiefly restricted to Southeast Asia. In this study, we analyzed the genomic evolution of DHAV-1 strains using extant GenBank records and genomic sequences of 10 DHAV-1 strains originating from a large disease outbreak in 2004-2005, in Hungary. Recombination analysis revealed intragenotype recombination within DHAV-1 as well as intergenotype recombination events involving DHAV-1 and DHAV-3 strains. The intergenotype recombination occurred in the VP0 region. Diversifying selection seems to act at sites of certain genomic regions. Calculations estimated slightly lower rates of evolution of DHAV-1 (mean rates for individual protein coding regions, 5.6286 × 10-4 to 1.1147 × 10-3 substitutions per site per year) compared to other picornaviruses. The observed evolutionary mechanisms indicate that whole-genome-based analysis of DHAV strains is needed to better understand the emergence of novel strains and their geographical dispersal.

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.

Increased spatial resolution of sampling in the Carpathian basin helps to understand the phylogeny of central European stream-dwelling gudgeons.

BACKGROUND: Phylogenetic studies of widespread European fish species often do not completely cover their entire distribution area, and some areas are often excluded from analyses than others. For example, Carpathian stocks are often omitted from these surveys or are under-represented in the samples. However, this area served as an extra-Mediterranean refugia for many species; therefore, it is assumed that fish stocks here may show special phylogenetic features. For this reason, increased spatial resolution of sampling, namely revealing genetic information from unexamined Carpathian areas within the range of doubtful taxa, may help us better understand their phylogenetic features. To test this hypothesis, a phylogenetic investigation using a partial mtCR sequence data was conducted on 56 stream-dwelling freshwater fish (Gobio spp.) individuals collected from 11 rivers of the data-deficient Southeastern Carpathian area. Moreover, we revieved the available phylogenetic data of Middle-Danubian stream-dwelling gudgeon lineages to delineate their distribution in the area. RESULTS: Seven out of the nine detected haplotypes were newly described, suggesting the studied area hosts distinct and diverse Gobio stocks. Two valid species (G. obtusirostris, G. gobio), and a haplogroup with doubtful phylogenetic position" G. sp. 1" were detected in the area, showing a specific spatial distribution pattern. The distribution of the detected lineages in the Middle-Danubian area correspond with recent and paleo hydrogeographic features, at the same time mainly on their bordering areas show considerable overlap. CONCLUSIONS: Despite the relatively limited geographic range of the study, our results provide important information which contributes to a better understanding of the phylogenetic, taxonomic and distribution features of Central European gudgeons. The genetically confirmed distribution data of the indicated lineages corresponds well with the recent and near-recent hydrogeographic features of the area, and assumes several hybrid zones in the Carpathian Basin. Additionally, the results show that the middle and lower Danubian watershed cannot be excluded from the range of G. gobio. Moreover, the" G. sp. 1", is slightly differentiated but phylogenetically distinct entity, and is the only Gobio taxa thus far detected in the middle and lower Tisza-basin. However, further investigations are necessary to clarify the taxonomic position of this group.

Whole genome sequencing of a rare rotavirus from archived stool sample demonstrates independent zoonotic origin of human G8P[14] strains in Hungary.

Genotype P[14] rotaviruses in humans are thought to be zoonotic strains originating from bovine or ovine host species. Over the past 30 years only few genotype P[14] strains were identified in Hungary totaling<0.1% of all human rotaviruses whose genotype had been determined. In this study we report the genome sequence and phylogenetic analysis of a human genotype G8P[14] strain, RVA/Human-wt/HUN/182-02/2001/G8P[14]. The whole genome constellation (G8-P[14]-I2-R2-C2-M2-A11-N2-T6-E2-H3) of this strain was shared with another Hungarian zoonotic G8P[14] strain, RVA/Human-wt/HUN/BP1062/2004/G8P[14], although phylogenetic analyses revealed the two rotaviruses likely had different progenitors. Overall, our findings indicate that human G8P[14] rotavirus detected in Hungary in the past originated from independent zoonotic events. Further studies are needed to assess the public health risk associated with infections by various animal rotavirus strains.

Development of a consensus reverse transcription PCR assay for the specific detection of tortoise picornaviruses.

Picornaviruses (PVs) of different terrestrial tortoise species, previously designated as Virus "X," have been frequently detected from various tissues by virus isolation in Terrapene heart cell culture as the preferred laboratory method for diagnosis. Here, we describe the development of 2 diagnostic reverse transcription (RT)-PCR-based assays for the identification and characterization of tortoise PVs belonging to the tentative genus Topivirus To test the novel diagnostic systems, PVs were isolated from swab and tissue samples collected in Germany, Italy, and Hungary between 2000 and 2013. All 25 tested isolates gave positive results with both novel consensus primer sets. Sequencing of the amplified products confirmed that all studied viruses were members of the new proposed genus Topivirus Phylogenetic analyses clearly distinguished 2 lineages within the genus. Based on sequence analysis, no association was observed between the geographic distribution and genetic relatedness. Furthermore, no strict host specificity was indicated. The PCR-based diagnosis may provide a time-saving and sensitive method to detect tortoise PVs, and evaluation of PV presence in these animals may help control virus spread.

Complete Genome Sequence of a Genotype G23P[37] Pheasant Rotavirus Strain Identified in Hungary.

We investigated the genomic properties of a rotavirus A strain isolated from diarrheic pheasant poults in Hungary in 2015. Sequence analyses revealed a shared genomic constellation (G23-P[37]-I4-R4-C4-M4-A16-N10-T4-E4-H4) and close relationship (range of nucleotide sequence similarity: VP2, 88%; VP1 and NSP4, 98%) with another pheasant rotavirus strain isolated previously in Germany.

Whole genome characterization of a chelonian orthoreovirus strain identifies significant genetic diversity and may classify reptile orthoreoviruses into distinct species.

In this study we report the sequence and phylogenetic characterization of an orthoreovirus strain, CH1197/96, isolated from a spur-thighed tortoise (Testudo graeca) on chicken embryo fibroblast cells. The 23,957 bp long genome sequence was obtained by combined use of semiconductor and capillary sequencing. Although the genomic characterization showed that the virus was most similar to the bush viper reovirus strain, 47/02, and in phylogenies performed with all segments the two strains formed a monophyletic group, the nucleotide (48.4-70.3%) and amino acid (39.2-80.7%) sequence identity values were moderate between the two reptile origin reoviruses. Based on our results and existing classification criteria for the genus Orthoreovirus, the tortoise reovirus strain CH1197/96 might be the first representative of a novel reptilian origin Orthoreovirus species.

Coding-complete sequencing classifies parrot bornavirus 5 into a novel virus species.

In this study, we determined the sequence of the coding region of an avian bornavirus detected in a blue-and-yellow macaw (Ara ararauna) with pathological/histopathological changes characteristic of proventricular dilatation disease. The genomic organization of the macaw bornavirus is similar to that of other bornaviruses, and its nucleotide sequence is nearly identical to the available partial parrot bornavirus 5 (PaBV-5) sequences. Phylogenetic analysis showed that these strains formed a monophyletic group distinct from other mammalian and avian bornaviruses and in calculations performed with matrix protein coding sequences, the PaBV-5 and PaBV-6 genotypes formed a common cluster, suggesting that according to the recently accepted classification system for bornaviruses, these two genotypes may belong to a new species, provisionally named Psittaciform 2 bornavirus.

Ubiquiter circovirus sequences raise challenges in laboratory diagnosis: the case of honey bee and bee mite, reptiles, and free living amoebae.

Circoviruses of pigs and birds are established pathogens, however, the exact role of other, recently described circoviruses and circovirus-like viruses remains to be elucidated. The aim of this study was the detection of circoviruses in neglected host species, including honey bees, exotic reptiles and free-living amoebae by widely used broad-spectrum polymerase chain reaction (PCR) assays specific for the replication initiation protein coding gene of these viruses. The majority of sequences obtained from honey bees were highly similar to canine and porcine circoviruses, or, were distantly related to dragonfly cycloviruses. Other rep sequences detected in some honey bees, reptiles and amoebae showed similarities to various rep sequences deposited in the GenBank. Back-to-back PCR primers designed for the amplification of whole viral genomes failed to work that suggested the existence of integrated rep-like elements in many samples. Rolling circle amplification and exonuclease treatment confirmed the absence of small circular DNA genomes in the specimens analysed. In case of honey bees Varroa mite DNA contamination might be a source of the identified endogenous rep-like elements. The reptile and amoebae rep-like sequences were nearly identical with each other and with sequences detected in chimpanzee feces raising the possibility that detection of novel or unusual rep-like elements in some host species might originate from the microbial community of the host. Our results indicate that attention is needed when broad-spectrum rep gene specific polymerase chain reaction is chosen for laboratory diagnosis of circovirus infections.

Sequencing and phylogenetic analysis identifies candidate members of a new picornavirus genus in terrestrial tortoise species.

Near-complete genome sequences of seven picornavirus (PV) strains isolated from different terrestrial tortoise species were determined and characterized. The genome organization of the strains proved to be similar and displayed a typical PV layout, and the polyprotein-encoding regions showed low similarity to those of other PVs. The predicted regions of the tortoise PV genomes were related to the corresponding genome parts of viruses belonging to distinct genera, implying modular evolution of these novel viruses. Our results suggest that these tortoise PVs belong to a prototype species in a separate proposed genus in the family Picornaviridae, tentatively called Topivirus (Tortoise picornavirus).

Whole-genome sequencing of a green bush viper reovirus reveals a shared evolutionary history between reptilian and unusual mammalian orthoreoviruses.

In this study, we sequenced the whole genome of a reovirus isolated from a green bush viper (Atheris squamigera). The bush viper reovirus shared several features with other orthoreoviruses, including its genome organization. In phylogenetic analysis, this strain was monophyletic with Broome virus and baboon orthoreovirus, indicating that these viruses might have originated from a common ancestor. These new molecular data supplement previous information based mainly on biological properties of reptilian reoviruses, confirming their taxonomic position and broadening our knowledge of the evolution of members of the genus Orthoreovirus.