Resistome analysis of Escherichia coli isolates from layers in Hungary.

The authors aimed to investigate eight strains of Escherichia coli (E. coli) strains from Hungarian layer flocks for antimicrobial resistance genes (ARG), using metagenomic methods. The strains were isolated from cloacal swabs of healthy adult layers. This study employed shotgun sequencing-based genetic and bioinformatic analysis along with determining phenotypic minimum inhibitory concentrations. A total of 59 ARGs were identified in the eight E. coli isolates, carrying ARGs against 15 groups of antibiotics. Among these, 28 ARGs were identified as transferable. Specifically, four ARGs were plasmid-derived, 18 ARGs were phage-derived and an additional six ARGs were predicted to be mobile, contributing to their mobility and potential spread between bacteria.

Phylogeny of Transferable Oxazolidinone Resistance Genes and Homologs.

Oxazolidinone resistance, especially transmissible resistance, is a major public health concern, and the origin of this resistance mechanism is not yet resolved. This study aims to delve into the phylogenetic origin of the transmissible oxazolidinone resistance mechanisms conferring cross-resistance to other drugs of human and veterinary importance. The amino acid sequences of the five cfr ribosomal methylases and optrA and poxtA were used as queries in searches against 219,549 bacterial proteomes in the NCBI RefSeq database. Hits with >40% amino acid identity and >80% query coverage were aligned, and phylogenetic trees were reconstructed. All five cfr genes yielded highly similar trees, with rlmN housekeeping ribosomal methylases located basal to the sister groups of S-adenosyl-methionine-dependent methyltransferases from various Deltaproteobacteria and Actinomycetia, including antibiotic-producing Streptomyces species, and the monophyletic group of cfr genes. The basal branches of the latter contained paenibacilli and other soil bacteria; they then could be split into the clades [cfr(C):cfr(E)] and [[cfr:cfr(B)]:cfr(D)], always with different Bacillaceae in their stems. Lachnospiraceae were encountered in the basal branches of both optrA and poxtA trees. The ultimate origin of the cfr genes is the rlmN housekeeping ribosomal methylases, which evolved into a suicide-avoiding methylase in antibiotic producers; a soil organism (Lachnospiraceae, Paenibacilli) probably acted as a transfer organism into pathogenic bacteria. In the case of optrA, the porcine pathogenic Streptococcus suis was present in all branches, while the proteins closest to poxtA originated from Clostridia.

Genomic epidemiology of antifungal resistance in human and avian isolates of Candida albicans: a pilot study from the One Health perspective.

Stress-induced genomic changes in Candida albicans contribute to the adaptation of this species to various environmental conditions. Variations of the genome composition of animal-origin C. albicans strains are largely unexplored and drug resistance or other selective pressures driving the evolution of these yeasts remained an intriguing question. Comparative genome analysis was carried out to uncover chromosomal aneuploidies and regions with loss of heterozygosity (LOH), two mechanisms that manage genome plasticity. We detected aneuploidy only in human isolates. Bird-derived isolates showed LOH in genes commonly associated with antifungal drug resistance similar to human isolates. Our study suggests that environmental fungicide usage might exert selective pressure on C. albicans infecting animals, thus contributing to the spread of potentially resistant strains between different hosts.

In Vitro Microevolution and Co-Selection Assessment of Amoxicillin and Cefotaxime Impact on Escherichia coli Resistance Development.

The global spread of antimicrobial resistance has become a prominent issue in both veterinary and public health in the 21st century. The extensive use of amoxicillin, a beta-lactam antibiotic, and consequent resistance development are particularly alarming in food-producing animals, with a focus on the swine and poultry sectors. Another beta-lactam, cefotaxime, is widely utilized in human medicine, where the escalating resistance to third- and fourth-generation cephalosporins is a major concern. The aim of this study was to simulate the development of phenotypic and genotypic resistance to beta-lactam antibiotics, focusing on amoxicillin and cefotaxime. The investigation of the minimal inhibitory concentrations (MIC) of antibiotics was performed at 1×, 10×, 100×, and 1000× concentrations using the modified microbial evolution and growth arena (MEGA-plate) method. Our results indicate that amoxicillin significantly increased the MIC values of several tested antibiotics, except for oxytetracycline and florfenicol. In the case of cefotaxime, this increase was observed in all classes. A total of 44 antimicrobial resistance genes were identified in all samples. Chromosomal point mutations, particularly concerning cefotaxime, revealed numerous complex mutations, deletions, insertions, and single nucleotide polymorphisms (SNPs) that were not experienced in the case of amoxicillin. The findings suggest that, regarding amoxicillin, the point mutation of the acrB gene could explain the observed MIC value increases due to the heightened activity of the acrAB-tolC efflux pump system. However, under the influence of cefotaxime, more intricate processes occurred, including complex amino acid substitutions in the ampC gene promoter region, increased enzyme production induced by amino acid substitutions and SNPs, as well as mutations in the acrR and robA repressor genes that heightened the activity of the acrAB-tolC efflux pump system. These changes may contribute to the significant MIC increases observed for all tested antibiotics. The results underscore the importance of understanding cross-resistance development between individual drugs when choosing clinical alternative drugs. The point mutations in the mdtB and emrR genes may also contribute to the increased activity of the mdtABC-tolC and emrAB-tolC pump systems against all tested antibiotics. The exceptionally high mutation rate induced by cephalosporins justifies further investigations to clarify the exact mechanism behind.

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.

Structural similarity of human papillomavirus E4 and polyomaviral VP4 exhibited by genomic analysis of the common kestrel (Falco tinnunculus) polyomavirus.

Polyomaviruses are widely distributed viruses of birds that may induce developmental deformities and internal organ disorders primarily in nestlings. In this study, polyomavirus sequence was detected in kidney and liver samples of a common kestrel (Falco tinnunculus) that succumbed at a rescue station in Hungary. The amplified 5025 nucleotide (nt) long genome contained the early (large and small T antigen, LTA and STA) and late (viral proteins, VP1, VP2, VP3) open reading frames (ORFs) typical for polyomaviruses. One of the additional putative ORFs (named VP4) showed identical localization with the VP4 and ORF-X of gammapolyomaviruses, but putative splicing sites could not be found in its sequence. Interestingly, the predicted 123 amino acid (aa) long protein sequence showed the highest similarity with human papillomavirus E4 early proteins in respect of the aa distribution and motif arrangement implying similar functions. The LTA of the kestrel polyomavirus shared <59.2% nt and aa pairwise identity with the LTA sequence of other polyomaviruses and formed a separated branch in the phylogenetic tree among gammapolyomaviruses. Accordingly, the kestrel polyomavirus may be the first member of a novel species within the Gammapolyomavirus genus, tentatively named Gammapolyomavirus faltin.

Antimicrobial resistance genes and associated mobile genetic elements in Lactobacillales from various sources.

Lactobacillales are commonly used in food products and as probiotics in animal and human medicine. Despite being generally recognized as safe, lactic acid bacteria may harbor a variety of antimicrobial resistance genes (ARGs), which may be transferable to human or veterinary pathogens, thus, may pose veterinary and public health concerns. This study investigates the resistome of Lactobacillales. A total of 4,286 whole-genome sequences were retrieved from NCBI RefSeq database. We screened ARGs in whole genome sequences and assessed if they are transmissible by plasmid transfer or by linkage to integrative mobile genetic elements. In the database, 335 strains were found to carry at least one ARG, and 194 strains carried at least one potentially transferable ARG. The most prevalent transferable ARG were tetM and tetW conferring antibiotic resistance to tetracycline. This study highlights the importance of the One Health concept by demonstrating the potential for Lactobacillales, commonly used in food products, to serve as reservoirs and vectors for ARGs.

In Vitro Microevolution and Co-Selection Assessment of Florfenicol Impact on Escherichia coli Resistance Development.

The issue of antimicrobial resistance is becoming an increasingly serious challenge in both human and veterinary medicine. Prudent antimicrobial use in veterinary medicine is warranted and supported by international guidelines, with the Antimicrobial Advice Ad Hoc Expert Group (AMEG) placing particular emphasis on the critically important group B antimicrobials. These antimicrobials are commonly employed, especially in the poultry and swine industry. The impact of florfenicol, a veterinary antibiotic, was studied on the resistance development of Escherichia coli. The aim of the study was to investigate the effect of the use of florfenicol on the development of phenotypic and genomic resistances, not only to the drug itself but also to other drugs. The minimum inhibitory concentrations (MICs) of the antibiotics were investigated at 1×, 10×, 100× and 1000× concentrations using the adapted Microbial Evolution and Growth Arena (MEGA-plate) method. The results demonstrate that florfenicol can select for resistance to fluoroquinolone antibiotics (167× MIC value increase) and cephalosporins (67× MIC value increase). A total of 44 antimicrobial resistance genes were identified, the majority of which were consistent across the samples. Chromosomal point mutations, including alterations in resistance-associated and regulatory genes (acrB, acrR, emrR and robA), are thought to trigger multiple drug efflux pump activations, leading to phenotypically increased resistance. The study underscores the impact of florfenicol and its role in the development of antimicrobial resistance, particularly concerning fluoroquinolone antibiotics and cephalosporins. This study is the first to report florfenicol's dose-dependent enhancement of other antibiotics' MICs, linked to mutations in SOS-box genes (mdtABC-tolC, emrAB-tolC and acrAB-tolC) and increased multidrug efflux pump genes. Mutations in the regulatory genes acrR, emrR and rpbA support the possibility of increased gene expression. The results are crucial for understanding antimicrobial resistance and its development, highlighting the promising potential of in vitro evolutionary and coselection studies for future research.

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.

Carbapenem-resistant Escherichia coli in Black-headed gulls, the Danube, and human clinical samples: A One Health comparison of contemporary isolates.

OBJECTIVES: Our aim was to characterize and compare contemporary carbapenem-resistant Enterobacterales (CRE) isolates from gulls, the River Danube, and humans in Hungary, Budapest. METHODS: Multiresistant Enterobacterales were sought for in 227 gull faecal and 24 Danube water samples from 2019 to 2020. Eosin-methylene blue agar containing 2 mg/L cefotaxime and Colilert-test containing 10 mg/L cefotaxime were used for gull and water samples, respectively. Isolates were characterized by polymerase chain reactions (PCRs); acquired carbapenemase producers were further analysed by whole-genome sequencing, together with 21 Hungarian human CR Escherichia coli (CREc) isolates. RESULTS: Gull and water samples exhibited a CRE prevalence of 7.4% (9/122) and 6.7% (7/105), none and 5/12 water samples yielded CRE from 2019 and 2020, respectively; CRE were found only in samples taken downstream of Budapest. The dominant species was Escherichia coli and the most prevalent carbapenemase was blaNDM-1. High-risk CREc clones were found both in gulls (ST224, ST372, ST744) and the Danube (ST10, ST354, ST410); the closest associations were between ST410 from humans and the Danube, among ST1437 among gulls, and between ST1437 in gulls and the Danube (46, 0, and 22-24 allelic distances, respectively). Direct links between human and gull isolates were not demonstrated. CONCLUSION: The study demonstrates potential epidemiological links among humans, a river crossing a city, and urbanised birds, suggesting a local transmission network. Water bodies receiving influent wastewater, together with animals using such habitats, may serve as a local reservoir system for CRE, highlighting the importance of One Health in CRE transmission, even in a country with a low CRE prevalence in humans.

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.

Genome stability assessment of PRRS vaccine strain with new ARTIC-style sequencing protocol.

A tiling amplicon sequencing protocol was developed to analyse the genome sequence stability of the modified live PRRSV vaccine strain, Porcilis MLV. The backbone of the ARTIC-style protocol was formed by 34 individual primer pairs, which were divided into two primer pools. Primer pairs were designed to amplify 532 to 588 bp fragments of the corresponding genomic region. The amplicons are suitable for sequencing on Illumina DNA sequencers with available 600-cycle sequencing kits. The concentration of primer pairs in the pools was optimized to obtain a balanced sequencing depth along the genome. Deep sequencing data of three vaccine batches were also analysed. All three vaccine batches were very similar to each other, although they also showed single nucleotide variations (SNVs) affecting less than 1 % of the genome. In the three vaccine strains, 113 to 122 SNV sites were identified; at these sites, the minority variants represented a frequency range of 1 to 48.7 percent. Additionally, the strains within the batches contained well-known length polymorphisms; the genomes of these minority deletion mutants were 135 to 222 bp shorter than the variant with the complete genome. Our results show the usefulness of ARTIC-style protocols in the evaluation of the genomic stability of PRRS MLV strains.

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

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

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.

Rapid genotyping of targeted viral samples using Illumina short-read sequencing data.

The most important information about microorganisms might be their accurate genome sequence. Using current Next Generation Sequencing methods, sequencing data can be generated at an unprecedented pace. However, we still lack tools for the automated and accurate reference-based genotyping of viral sequencing reads. This paper presents our pipeline designed to reconstruct the dominant consensus genome of viral samples and analyze their within-host variability. We benchmarked our approach on numerous datasets and showed that the consensus genome of samples could be obtained reliably without further manual data curation. Our pipeline can be a valuable tool for fast identifying viral samples. The pipeline is publicly available on the project's GitHub page (https://github.com/laczkol/QVG).

A novel gammapolyomavirus in a great cormorant (Phalacrocorax carbo).

In this study, the complete genome of a novel polyomavirus detected in a great cormorant (Phalacrocorax carbo) was characterized. The 5133-bp-long genome of the cormorant polyomavirus has a genomic structure typical of members of the genus Gammapolyomavirus, family Polyomaviridae, containing open reading frames encoding the large and small tumor antigens, viral proteins 1, 2, and 3, and the X protein. The large tumor antigen of the cormorant polyomavirus shares 45.6-50.4% amino acid sequence identity with the homologous sequences of other gammapolyomaviruses. These data, together with results of phylogenetic analysis, suggest that this cormorant polyomavirus should be considered the first member of a new species within the genus Gammapolyomavirus, for which we propose the name "Phalacrocorax carbo polyomavirus 1".

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.