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.

Apis mellifera filamentous virus from a honey bee gut microbiome survey in Hungary.

In Hungary, as part of a nationwide, climatically balanced survey for a next-generation sequencing-based study of the honey bee (Apis mellifera) gut microbiome, repeated sampling was carried out during the honey production season (March and May 2019). Among other findings, the presence of Apis mellifera filamentous virus (AmFV) was detected in all samples, some at very high levels. AmFV-derived reads were more abundant in the March samples than in the May samples. In March, a higher abundance of AmFV-originated reads was identified in samples collected from warmer areas compared to those collected from cooler areas. A lower proportion of AmFV-derived reads were identified in samples collected in March from the wetter areas than those collected from the drier areas. AmFV-read abundance in samples collected in May showed no significant differences between groups based on either environmental temperature or precipitation. The AmFV abundance correlated negatively with Bartonella apihabitans, Bartonella choladocola, and positively with Frischella perrara, Gilliamella apicola, Gilliamella sp. ESL0443, Lactobacillus apis, Lactobacillus kullabergensis, Lactobacillus sp. IBH004. De novo metagenome assembly of four samples resulted in almost the complete AmFV genome. According to phylogenetic analysis based on DNA polymerase, the Hungarian strains are closest to the strain CH-05 isolated in Switzerland.

Genetic Factors Associated with the Development of Neuropathy in Type 2 Diabetes.

Neuropathy is a serious and frequent complication of type 2 diabetes (T2DM). This study was carried out to search for genetic factors associated with the development of diabetic neuropathy by whole exome sequencing. For this study, 24 patients with long-term type 2 diabetes with neuropathy and 24 without underwent detailed neurological assessment and whole exome sequencing. Cardiovascular autonomic function was evaluated by cardiovascular reflex tests. Heart rate variability was measured by the triangle index. Sensory nerve function was estimated by Neurometer and Medoc devices. Neuropathic symptoms were characterized by the neuropathy total symptom score (NTSS). Whole exome sequencing (WES) was performed on a Thermo Ion GeneStudio S5 system determining the coding sequences of approximately 32,000 genes comprising 50 million base pairs. Variants were detected by Ion Reporter software and annotated using ANNOVAR, integrating database information from dbSNP, ClinVar, gnomAD, and OMIM. Integrative genomics viewer (IGV) was used for visualization of the mapped reads. We have identified genetic variants that were significantly associated with increased (22-49-fold) risk of neuropathy (rs2032930 and rs2032931 of recQ-mediated genome instability protein 2 (RMI2) gene), rs604349 of myosin binding protein H like (MYBPHL) gene and with reduced (0.07-0.08-fold) risk (rs917778 of multivesicular body subunit 12B (MVB12B) and rs2234753 of retinoic acid X receptor alpha (RXRA) genes). The rs2032930 showed a significant correlation with current perception thresholds measured at 5 Hz and 250 Hz for n. medianus (p = 0.042 and p = 0.003, respectively) and at 5 Hz for n. peroneus (p = 0.037), as well as the deep breath test (p = 0.022) and the NTSS (p = 0.023). The rs2032931 was associated with current perception thresholds (p = 0.003 and p = 0.037, respectively), deep breath test (p = 0.022), and NTSS (p = 0.023). The rs604349 correlated with values measured at 2000 (p = 0.049), 250 (p = 0.018), and 5 Hz (p = 0.005) for n. medianus, as well as warm perception threshold measured by Medoc device (p = 0.042). The rs2234753 showed correlations with a current perception threshold measured at 2000 Hz for n. medianus (p = 0.020), deep breath test (p = 0.040), and NTSS (p = 0.003). There was a significant relationship between rs91778 and cold perception threshold (p = 0.013). In our study, genetic variants have been identified that may have an impact on the risk of neuropathy developing in type 2 diabetic patients. These results could open up new opportunities for early preventive measures and might provide targets for new drug developments in the future.

Ixodes ricinus tick bacteriome alterations based on a climatically representative survey in Hungary.

IMPORTANCE: Climate-sensitive disease vectors, such as ticks, respond to the environment with changes in their microbiome. These changes can affect the emergence or re-emergence of various vector-borne pathogens, such as the causative agent of Lyme borreliosis (LB) or tick-borne encephalitis. This aspect is particularly emphasized in light of climate change. The climatically representative assessment of microbiome differences in various developmental stages of the most common Central European tick species, Ixodes ricinus, deepens our understanding of the potential climatic factors behind microbial relative abundance and interaction changes. This knowledge can support the development of novel disease vector control strategies.

Procalcitonin-guided antibiotic therapy may shorten length of treatment and may improve survival-a systematic review and meta-analysis.

BACKGROUND: Appropriate antibiotic (AB) therapy remains a challenge in the intensive care unit (ICU). Procalcitonin (PCT)-guided AB stewardship could help optimize AB treatment and decrease AB-related adverse effects, but firm evidence is still lacking. Our aim was to compare the effects of PCT-guided AB therapy with standard of care (SOC) in critically ill patients. METHODS: We searched databases CENTRAL, Embase and Medline. We included randomized controlled trials (RCTs) comparing PCT-guided AB therapy (PCT group) with SOC reporting on length of AB therapy, mortality, recurrent and secondary infection, ICU length of stay (LOS), hospital LOS or healthcare costs. Due to recent changes in sepsis definitions, subgroup analyses were performed in studies applying the Sepsis-3 definition. In the statistical analysis, a random-effects model was used to pool effect sizes. RESULTS: We included 26 RCTs (n = 9048 patients) in the quantitative analysis. In comparison with SOC, length of AB therapy was significantly shorter in the PCT group (MD - 1.79 days, 95% CI: -2.65, - 0.92) and was associated with a significantly lower 28-day mortality (OR 0.84, 95% CI: 0.74, 0.95). In Sepsis-3 patients, mortality benefit was more pronounced (OR 0.46 95% CI: 0.27, 0.79). Odds of recurrent infection were significantly higher in the PCT group (OR 1.36, 95% CI: 1.10, 1.68), but there was no significant difference in the odds of secondary infection (OR 0.81, 95% CI: 0.54, 1.21), ICU and hospital length of stay (MD - 0.67 days 95% CI: - 1.76, 0.41 and MD - 1.23 days, 95% CI: - 3.13, 0.67, respectively). CONCLUSIONS: PCT-guided AB therapy may be associated with reduced AB use, lower 28-day mortality but higher infection recurrence, with similar ICU and hospital length of stay. Our results render the need for better designed studies investigating the role of PCT-guided AB stewardship in critically ill patients.

A biallelic multiple nucleotide length polymorphism explains functional causality at 5p15.33 prostate cancer risk locus.

To date, single-nucleotide polymorphisms (SNPs) have been the most intensively investigated class of polymorphisms in genome wide associations studies (GWAS), however, other classes such as insertion-deletion or multiple nucleotide length polymorphism (MNLPs) may also confer disease risk. Multiple reports have shown that the 5p15.33 prostate cancer risk region is a particularly strong expression quantitative trait locus (eQTL) for Iroquois Homeobox 4 (IRX4) transcripts. Here, we demonstrate using epigenome and genome editing that a biallelic (21 and 47 base pairs (bp)) MNLP is the causal variant regulating IRX4 transcript levels. In LNCaP prostate cancer cells (homozygous for the 21 bp short allele), a single copy knock-in of the 47 bp long allele potently alters the chromatin state, enabling de novo functional binding of the androgen receptor (AR) associated with increased chromatin accessibility, Histone 3 lysine 27 acetylation (H3K27ac), and ~3-fold upregulation of IRX4 expression. We further show that an MNLP is amongst the strongest candidate susceptibility variants at two additional prostate cancer risk loci. We estimated that at least 5% of prostate cancer risk loci could be explained by functional non-SNP causal variants, which may have broader implications for other cancers GWAS. More generally, our results underscore the importance of investigating other classes of inherited variation as causal mediators of human traits.

Development of a versatile LCM-Seq method for spatial transcriptomics of fluorescently tagged cholinergic neuron populations.

Single-cell transcriptomics are powerful tools to define neuronal cell types based on co-expressed gene clusters. Limited RNA input in these technologies necessarily compromises transcriptome coverage and accuracy of differential expression analysis. We propose that bulk RNA-Seq of neuronal pools defined by spatial position offers an alternative strategy to overcome these technical limitations. We report a laser-capture microdissection (LCM)-Seq method which allows deep transcriptome profiling of fluorescently tagged neuron populations isolated with LCM from histological sections of transgenic mice. Mild formaldehyde fixation of ZsGreen marker protein, LCM sampling of ∼300 pooled neurons, followed by RNA isolation, library preparation and RNA-Seq with methods optimized for nanogram amounts of moderately degraded RNA enabled us to detect ∼15,000 different transcripts in fluorescently labeled cholinergic neuron populations. The LCM-Seq approach showed excellent accuracy in quantitative studies, allowing us to detect 2891 transcripts expressed differentially between the spatially defined and clinically relevant cholinergic neuron populations of the dorsal caudate-putamen and medial septum. In summary, the LCM-Seq method we report in this study is a versatile, sensitive, and accurate bulk sequencing approach to study the transcriptome profile and differential gene expression of fluorescently tagged neuronal populations isolated from transgenic mice with high spatial precision.

Antimicrobial resistance gene lack in tick-borne pathogenic bacteria.

Tick-borne infections, including those of bacterial origin, are significant public health issues. Antimicrobial resistance (AMR), which is one of the most pressing health challenges of our time, is driven by specific genetic determinants, primarily by the antimicrobial resistance genes (ARGs) of bacteria. In our work, we investigated the occurrence of ARGs in the genomes of tick-borne bacterial species that can cause human infections. For this purpose, we processed short/long reads of 1550 bacterial isolates of the genera Anaplasma (n = 20), Bartonella (n = 131), Borrelia (n = 311), Coxiella (n = 73), Ehrlichia (n = 13), Francisella (n = 959) and Rickettsia (n = 43) generated by second/third generation sequencing that have been freely accessible at the NCBI SRA repository. From Francisella tularensis, 98.9% of the samples contained the FTU-1 beta-lactamase gene. However, it is part of the F. tularensis representative genome as well. Furthermore, 16.3% of them contained additional ARGs. Only 2.2% of isolates from other genera (Bartonella: 2, Coxiella: 8, Ehrlichia: 1, Rickettsia: 2) contained any ARG. We found that the odds of ARG occurrence in Coxiella samples were significantly higher in isolates related to farm animals than from other sources. Our results describe a surprising lack of ARGs in these bacteria and suggest that Coxiella species in farm animal settings could play a role in the spread of AMR.

A survey on antimicrobial resistance genes of frequently used probiotic bacteria, 1901 to 2022.

BackgroundAntimicrobial resistance (AMR) is caused by AMR determinants, mainly genes (ARGs) in the bacterial genome. Bacteriophages, integrative mobile genetic elements (iMGEs) or plasmids can allow ARGs to be exchanged among bacteria by horizontal gene transfer (HGT). Bacteria, including bacteria with ARGs, can be found in food. Thus, it is conceivable that in the gastrointestinal tract, bacteria from the gut flora could take up ARGs from food.AimThe study objective was to gain insight into the ARG set carried by commonly used probiotic bacteria that may enter the human body with non-fermented foods, fermented foods, or probiotic dietary supplements (FFPs) and to assess ARG mobility.MethodsNext generation sequencing whole genome data from 579 isolates of 12 commonly employed probiotic bacterial species were collected from a public repository. Using bioinformatical tools, ARGs were analysed and linkage with mobile genetic elements assessed.ResultsResistance genes were found in eight bacterial species. The ratios of ARG positive/negative samples per species were: Bifidobacterium animalis (65/0), Lactiplantibacillus plantarum (18/194), Lactobacillus delbrueckii (1/40), Lactobacillus helveticus (2/64), Lactococcus lactis (74/5), Leucoconstoc mesenteroides (4/8), Levilactobacillus brevis (1/46), Streptococcus thermophilus (4/19). In 66% (112/169) of the ARG-positive samples, at least one ARG could be linked to plasmids or iMGEs. No bacteriophage-linked ARGs were found.ConclusionThe finding of potentially mobile ARGs in probiotic strains for human consumption raises awareness of a possibility of ARG HGT in the gastrointestinal tract. In addition to existing recommendations, screening FFP bacterial strains for ARG content and mobility characteristics might be considered.

Evaluation of the effects of Lake Hévíz sulfur thermal water on skin microbiome in plaque psoriasis: An open label, pilot study.

Psoriasis is a chronic inflammatory skin disease. It is associated with changes in skin microbiome. The aim of this study was to evaluate how Lake Hévíz sulfur thermal water influences the composition of microbial communities that colonizes skin in patients with psoriasis. Our secondary objective was to investigate the effects of balneotherapy on disease activity. In this open label study, participants with plaque psoriasis underwent 30-min therapy sessions in Lake Hévíz, at a temperature of 36 °C, five times a week for 3 weeks. The skin microbiome samples were collected by swabbing method from two different areas (lesional skin-psoriatic plaque and non-lesional skin). From 16 patients, 64 samples were processed for a 16S rRNA sequence-based microbiome analysis. Outcome measures were alpha-diversity (Shannon, Simpson, and Chao1 indexes), beta-diversity (Bray-Curtis metric), differences in genus level abundances, and Psoriasis Area and Severity Index (PASI). Skin microbiome samples were collected at baseline, and immediately after treatment. Based on the visual examination of the employed alpha- and beta-diversity measures, no systematic difference based on sampling timepoint or sample location could be revealed in these regards. Balneotherapy in the unaffected area significantly increased the level of Leptolyngbya genus, and significantly decreased the level of Flavobacterium genus. A similar trend was revealed by the results of the psoriasis samples, but the differences were not statistically significant. In patients with mild psoriasis, a significant improvement was observed in PASI scores.

Quantitative Analysis of Inflammatory Uterine Lesions of Pregnant Gilts with Digital Image Analysis Following Experimental PRRSV-1 Infection.

Reproductive disorders caused by porcine reproductive and respiratory syndrome virus-1 are not yet fully characterized. We report QuPath-based digital image analysis to count inflammatory cells in 141 routinely, and 35 CD163 immunohistochemically stained endometrial slides of vaccinated or unvaccinated pregnant gilts inoculated with a high or low virulent PRRSV-1 strain. To illustrate the superior statistical feasibility of the numerical data determined by digital cell counting, we defined the association between the number of these cells and endometrial, placental, and fetal features. There was strong concordance between the two manual scorers. Distributions of total cell counts and endometrial and placental qPCR results differed significantly between examiner1's endometritis grades. Total counts' distribution differed significantly between groups, except for the two unvaccinated. Higher vasculitis scores were associated with higher endometritis scores, and higher total cell counts were expected with high vasculitis/endometritis scores. Cell number thresholds of endometritis grades were determined. A significant correlation between fetal weights and total counts was shown in unvaccinated groups, and a significant positive correlation was found between these counts and endometrial qPCR results. We revealed significant negative correlations between CD163+ counts and qPCR results of the unvaccinated group infected with the highly virulent strain. Digital image analysis was efficiently applied to assess endometrial inflammation objectively.

Genomic Evidence for Direct Transmission of mecC-MRSA between a Horse and Its Veterinarian.

Methicillin-resistant Staphylococcus aureus bearing the mecC gene (mecC-MRSA) has been reported from animals and humans in recent years. This study describes the first mecC-MRSA isolates of human and equine origin in Hungary (two isolates from horses and one from a veterinarian, who treated one of the infected horses, but was asymptomatic). MRSA isolates were identified by cultivation and PCR detection of the species-specific spa gene and mecA/mecC methicillin resistance genes. The isolates were characterized by antibiotic susceptibility testing, MLST, spa, SCCmec typing, PFGE and whole genome sequencing (WGS). All three isolates belonged to the ST130-t843-SCCmec XI genotype, and carried the mecC and blaZ genes. Apart from beta-lactam drugs, they were sensitive to all tested antibiotics. The isolates of the infected horse and its veterinarian had the same PFGE pulsotype and showed only slight differences with WGS. Hence, this is the first description of direct transmission of a mecC-carrying MRSA between a horse and its veterinarian. The emergence of mecC in the country highlights the importance of the appropriate diagnostics in MRSA identification.

Comparison of the short-term complications of TTA-rapid and modified cTTA procedures.

The objective of this retrospective study was to determine the complications of the first 30 tibial tuberosity advancement rapid (TTA-rapid) and 30 modified circular tibial tuberosity advancement (mcTTA) procedures performed by our team, and to compare the results with the findings reported in the literature. Our research was based on 30 procedures in each group. All dogs were client-owned. Data were collected only for the study of cases that had a minimum follow-up period of 3 months. Intraoperative (IO) and postoperative (PO) complications were assessed, with the latter divided into two subgroups: major and minor. Results obtained for the TTA-rapid group: IO complications 23.3% (7/30), major PO complications 13.3% (4/30), minor PO complications 16.7% (5/30). Results of the mcTTA group: IO complications 0% (0/30), major PO complications 3.3% (1/30), minor PO complications 20% (6/30). Comparing the complication rates, we found that there was a significant difference between the two groups in the occurrence of IO complications (P = 0.01054); however, there was no significant difference in the incidence of major (P = 0.3533) and minor (P > 0.9999) PO complications between groups. Our results are consistent with the findings reported in the literature and suggest that both techniques are efficient and carry a relatively low complication rate.

Methodological and Biological Factors Influencing Global DNA Methylation Results Measured by LINE-1 Pyrosequencing Assay in Colorectal Tissue and Liquid Biopsy Samples.

Long interspersed nuclear element 1 (LINE-1) bisulfite pyrosequencing is a widely used technique for genome-wide methylation analyses. We aimed to investigate the effects of experimental and biological factors on its results to improve the comparability. LINE-1 bisulfite pyrosequencing was performed on colorectal tissue (n = 222), buffy coat (n = 39), and plasma samples (n = 9) of healthy individuals and patients with colorectal tumors. Significantly altered methylation was observed between investigated LINE-1 CpG positions of non-tumorous tissues (p ≤ 0.01). Formalin-fixed, paraffin-embedded biopsies (73.0 ± 5.3%) resulted in lower methylation than fresh frozen samples (76.1 ± 2.8%) (p ≤ 0.01). DNA specimens after long-term storage showed higher methylation levels (+3.2%, p ≤ 0.01). In blood collection tubes with preservatives, cfDNA and buffy coat methylation significantly changed compared to K3EDTA tubes (p ≤ 0.05). Lower methylation was detected in older (>40 years, 76.8 ± 1.7%) vs. younger (78.1 ± 1.0%) female patients (p ≤ 0.05), and also in adenomatous tissues with MTHFR 677CT, or 1298AC mutations vs. wild-type (p ≤ 0.05) comparisons. Based on our findings, it is highly recommended to consider the application of standard DNA samples in the case of a possible clinical screening approach, as well as in experimental research studies.

Correction: Papp, M.; Solymosi, N. Review and Comparison of Antimicrobial Resistance Gene Databases. Antibiotics 2022, 11, 339.

While preparing our next manuscript, we noticed a few errors in our published paper [...].

Natural diversity of the honey bee (Apis mellifera) gut bacteriome in various climatic and seasonal states.

As pollinators and producers of numerous human-consumed products, honey bees have great ecological, economic and health importance. The composition of their bacteriota, for which the available knowledge is limited, is essential for their body's functioning. Based on our survey, we performed a metagenomic analysis of samples collected by repeated sampling. We used geolocations that represent the climatic types of the study area over two nutritionally extreme periods (March and May) of the collection season. Regarding bacteriome composition, a significant difference was found between the samples from March and May. The samples' bacteriome from March showed a significant composition difference between cooler and warmer regions. However, there were no significant bacteriome composition differences among the climatic classes of samples taken in May. Based on our results, one may conclude that the composition of healthy core bacteriomes in honey bees varies depending on the climatic and seasonal conditions. This is likely due to climatic factors and vegetation states determining the availability and nutrient content of flowering plants. The results of our study prove that in order to gain a thorough understanding of a microbiome's natural diversity, we need to obtain the necessary information from extreme ranges within the host's healthy state.

Global DNA hypomethylation of colorectal tumours detected in tissue and liquid biopsies may be related to decreased methyl-donor content.

BACKGROUND: Hypomethylation of long interspersed nuclear element 1 (LINE-1) is characteristic of various cancer types, including colorectal cancer (CRC). Malfunction of several factors or alteration of methyl-donor molecules' (folic acid and S-adenosylmethionine) availability can contribute to DNA methylation changes. Detection of epigenetic alterations in liquid biopsies can assist in the early recognition of CRC. Following the investigations of a Hungarian colon tissue sample set, our goal was to examine the LINE-1 methylation of blood samples along the colorectal adenoma-carcinoma sequence and in inflammatory bowel disease. Moreover, we aimed to explore the possible underlying mechanisms of global DNA hypomethylation formation on a multi-level aspect. METHODS: LINE-1 methylation of colon tissue (n = 183) and plasma (n = 48) samples of healthy controls and patients with colorectal tumours were examined with bisulfite pyrosequencing. To investigate mRNA expression, microarray analysis results were reanalysed in silico (n = 60). Immunohistochemistry staining was used to validate DNA methyltransferases (DNMTs) and folate receptor beta (FOLR2) expression along with the determination of methyl-donor molecules' in situ level (n = 40). RESULTS: Significantly decreased LINE-1 methylation level was observed in line with cancer progression both in tissue (adenoma: 72.7 ± 4.8%, and CRC: 69.7 ± 7.6% vs. normal: 77.5 ± 1.7%, p ≤ 0.01) and liquid biopsies (adenoma: 80.0 ± 1.7%, and CRC: 79.8 ± 1.3% vs. normal: 82.0 ± 2.0%, p ≤ 0.01). However, no significant changes were recognized in inflammatory bowel disease cases. According to in silico analysis of microarray data, altered mRNA levels of several DNA methylation-related enzymes were detected in tumours vs. healthy biopsies, namely one-carbon metabolism-related genes-which met our analysing criteria-showed upregulation, while FOLR2 was downregulated. Using immunohistochemistry, DNMTs, and FOLR2 expression were confirmed. Moreover, significantly diminished folic acid and S-adenosylmethionine levels were observed in parallel with decreasing 5-methylcytosine staining in tumours compared to normal adjacent to tumour tissues (p ≤ 0.05). CONCLUSION: Our results suggest that LINE-1 hypomethylation may have a distinguishing value in precancerous stages compared to healthy samples in liquid biopsies. Furthermore, the reduction of global DNA methylation level could be linked to reduced methyl-donor availability with the contribution of decreased FOLR2 expression.

Modulation of Hepatic Insulin and Glucagon Signaling by Nutritional Factors in Broiler Chicken.

Influencing the endocrine metabolic regulation of chickens by nutritional factors might provide novel possibilities for improving animal health and productivity. This study was designed to evaluate the impact of dietary cereal type (wheat-based (WB) vs. maize-based (MB) diets), crude protein level (normal (NP) vs. lowered (LP)), and sodium (n-)butyrate (1.5 g/kg diet) supplementation (vs. no butyrate) on the responsiveness of hepatic glucagon receptor (GCGR), insulin receptor beta (IRß) and mammalian target of rapamycin (mTOR) in the phase of intensive growth of chickens. Liver samples of Ross 308 broiler chickens (Gallus gallus domesticus) were collected on day 21 for quantitative real-time polymerase chain reaction and Western blot analyses. Hepatic GCGR and mTOR gene expressions were up-regulated by WB and LP diet. GCGR and IRß protein level decreased in groups with butyrate supplementation; however, the quantity of IRß and mTOR protein increased in WB groups. Based on these data, the applied dietary strategies may be useful tools to modulate hepatic insulin and glucagon signaling of chickens in the period of intensive growth. The obtained results might contribute to the better understanding of glycemic control of birds and increase the opportunity of ameliorating insulin sensitivity, hence, improving the production parameters and the welfare of broilers.

Antimicrobial resistance determinants in silage.

Animal products may play a role in developing and spreading antimicrobial resistance in several ways. On the one hand, residues of antibiotics not adequately used in animal farming can enter the human body via food. However, resistant bacteria may also be present in animal products, which can transfer the antimicrobial resistance genes (ARG) to the bacteria in the consumer's body by horizontal gene transfer. As previous studies have shown that fermented foods have a meaningful ARG content, it is indicated that such genes may also be present in silage used as mass feed in the cattle sector. In our study, we aspired to answer what ARGs occur in silage and what mobility characteristics they have? For this purpose, we have analyzed bioinformatically 52 freely available deep sequenced silage samples from shotgun metagenome next-generation sequencing. A total of 16 perfect matched ARGs occurred 54 times in the samples. More than half of these ARGs are mobile because they can be linked to integrative mobile genetic elements, prophages or plasmids. Our results point to a neglected but substantial ARG source in the food chain.