Antimicrobial resistance has become an imminent concern for public health. As methods for detection and characterization of antimicrobial resistance move from targeted culture and polymerase chain reaction to high throughput metagenomics, appropriate resources for the analysis of large-scale data are required. Currently, antimicrobial resistance databases are tailored to smaller-scale, functional profiling of genes using highly descriptive annotations. Such characteristics do not facilitate the analysis of large-scale, ecological sequence datasets such as those produced with the use of metagenomics for surveillance. In order to overcome these limitations, we present MEGARes (https://megares.meglab.org), a hand-curated antimicrobial resistance database and annotation structure that provides a foundation for the development of high throughput acyclical classifiers and hierarchical statistical analysis of big data. MEGARes can be browsed as a stand-alone resource through the website or can be easily integrated into sequence analysis pipelines through download. Also via the website, we provide documentation for AmrPlusPlus, a user-friendly Galaxy pipeline for the analysis of high throughput sequencing data that is pre-packaged for use with the MEGARes database.

The specific antimicrobial resistance (AMR) decreases that can be expected from reducing antimicrobial (AM) use in U.S. beef production have not been defined. To address this data gap, feces were recovered from 36 lots of “raised without antibiotics” (RWA) and 36 lots of “conventional” (CONV) beef cattle. Samples ( n = 719) were collected during harvest and distributed over a year. AMR was assessed by (i) the culture of six AM-resistant bacteria (ARB), (ii) quantitative PCR (qPCR) for 10 AMR genes (ARGs), (iii) a qPCR array of 84 ARGs, and (iv) metagenomic sequencing. Generally, AMR levels were similar, but some were higher in CONV beef cattle. The prevalence of third-generation cephalosporin-resistant (3GC ^r ) Escherichia coli was marginally different between production systems (CONV, 47.5%; RWA, 34.8%; P = 0.04), but the seasonal effect (summer, 92.8%; winter, 48.3%; P < 0.01) was greater. Erythromycin-resistant (ERY ^r ) Enterococcus sp. concentrations significantly differed between production systems (CONV, 1.91 log ₁₀ CFU/g; RWA, 0.73 log ₁₀ CFU/g; P < 0.01). Levels of aadA1 , ant ( 6 )- I , bla _ACI , erm (A), erm (B), erm (C), erm (F), erm (Q), tet (A), tet (B), tet (M), and tet (X) ARGs were higher ( P < 0.05) in the CONV system. Aggregate abundances of all 43 ARGs detected by metagenomic sequencing and the aggregate abundances of ARGs in the aminoglycoside, β-lactam, macrolide-lincosamide-streptogramin B (MLS), and tetracycline AM classes did not differ (log ₂ fold change < 1.0) between CONV and RWA systems. These results suggest that further reductions of AM use in U.S. beef cattle production may not yield significant AMR reductions beyond MLS and tetracycline resistance. IMPORTANCE The majority of antimicrobial (AM) use in the United States is for food-animal production, leading to concerns that typical AM use patterns during “conventional” (CONV) beef cattle production in the United States contribute broadly to antimicrobial resistance (AMR) occurrence. In the present study, levels of AMR were generally similar between CONV and “raised without antibiotics” (RWA) cattle. Only a limited number of modest AMR increases was observed in CONV cattle, primarily involving macrolide-lincosamide-streptogramin B (MLS) and tetracycline resistance. Macrolides (tylosin) and tetracyclines (chlortetracycline) are administered in-feed for relatively long durations to reduce liver abscesses. To ensure judicious AM use, the animal health, economic, and AMR impacts of shorter duration in-feed administration of these AMs should be examined. However, given the modest AMR reductions observed, further reductions of AM use in U.S. beef cattle production may not yield significant AMR reductions beyond MLS and tetracycline resistance.

Shotgun metagenomic sequencing is increasingly utilized as a tool to evaluate ecological-level dynamics of antimicrobial resistance and virulence, in conjunction with microbiome analysis. Interest in use of this method for environmental surveillance of antimicrobial resistance and pathogenic microorganisms is also increasing. In published metagenomic datasets, the total of all resistance- and virulence-related sequences accounts for < 1% of all sequenced DNA, leading to limitations in detection of low-abundance resistome-virulome elements. This study describes the extent and composition of the low-abundance portion of the resistome-virulome, using a bait-capture and enrichment system that incorporates unique molecular indices to count DNA molecules and correct for enrichment bias. The use of the bait-capture and enrichment system significantly increased on-target sequencing of the resistome-virulome, enabling detection of an additional 1441 gene accessions and revealing a low-abundance portion of the resistome-virulome that was more diverse and compositionally different than that detected by more traditional metagenomic assays. The low-abundance portion of the resistome-virulome also contained resistance genes with public health importance, such as extended-spectrum betalactamases, that were not detected using traditional shotgun metagenomic sequencing. In addition, the use of the bait-capture and enrichment system enabled identification of rare resistance gene haplotypes that were used to discriminate between sample origins. These results demonstrate that the rare resistome-virulome contains valuable and unique information that can be utilized for both surveillance and population genetic investigations of resistance. Access to the rare resistome-virulome using the bait-capture and enrichment system validated in this study can greatly advance our understanding of microbiome-resistome dynamics.

Metagenomic investigations have the potential to provide unprecedented insights into microbial ecologies, such as those relating to antimicrobial resistance (AMR). We characterized the microbial resistome in livestock operations raising cattle conventionally (CONV) or without antibiotic exposures (RWA) using shotgun metagenomics. Samples of feces, wastewater from catchment basins, and soil where wastewater was applied were collected from CONV and RWA feedlot and dairy farms. After DNA extraction and sequencing, shotgun metagenomic reads were aligned to reference databases for identification of bacteria (Kraken) and antibiotic resistance genes (ARGs) accessions (MEGARes). Differences in microbial resistomes were found across farms with different production practices (CONV vs. RWA), types of cattle (beef vs. dairy), and types of sample (feces vs. wastewater vs. soil). Feces had the greatest number of ARGs per sample (mean = 118 and 79 in CONV and RWA, respectively), with tetracycline efflux pumps, macrolide phosphotransferases, and aminoglycoside nucleotidyltransferases mechanisms of resistance more abundant in CONV than in RWA feces. Tetracycline and macrolide–lincosamide–streptogramin classes of resistance were more abundant in feedlot cattle than in dairy cow feces, whereas the β-lactam class was more abundant in dairy cow feces. Lack of congruence between ARGs and microbial communities (procrustes analysis) suggested that other factors (e.g., location of farms, cattle source, management practices, diet, horizontal ARGs transfer, and co-selection of resistance), in addition to antimicrobial use, could have impacted resistome profiles. For that reason, we could not establish a cause–effect relationship between antimicrobial use and AMR, although ARGs in feces and effluents were associated with drug classes used to treat animals according to farms’ records (tetracyclines and macrolides in feedlots, β-lactams in dairies), whereas ARGs in soil were dominated by multidrug resistance. Characterization of the “resistance potential” of animal-derived and environmental samples is the first step toward incorporating metagenomic approaches into AMR surveillance in agricultural systems. Further research is needed to assess the public-health risk associated with different microbial resistomes.

The objective was to examine effects of treating commercial beef feedlot cattle with therapeutic doses of tulathromycin, a macrolide antimicrobial drug, on changes in the fecal resistome and microbiome using shotgun metagenomic sequencing. Two pens of cattle were used, with all cattle in one pen receiving metaphylaxis treatment (800 mg subcutaneous tulathromycin) at arrival to the feedlot, and all cattle in the other pen remaining unexposed to parenteral antibiotics throughout the study period. Fecal samples were collected from 15 selected cattle in each group just prior to treatment (Day 1), and again 11 days later (Day 11). Shotgun sequencing was performed on isolated metagenomic DNA, and reads were aligned to a resistance and a taxonomic database to identify alignments to antimicrobial resistance (AMR) gene accessions and microbiome content. Overall, we identified AMR genes accessions encompassing 9 classes of AMR drugs and encoding 24 unique AMR mechanisms. Statistical analysis was used to identify differences in the resistome and microbiome between the untreated and treated groups at both timepoints, as well as over time. Based on composition and ordination analyses, the resistome and microbiome were not significantly different between the two groups on Day 1 or on Day 11. However, both the resistome and microbiome changed significantly between these two sampling dates. These results indicate that the transition into the feedlot—and associated changes in diet, geography, conspecific exposure, and environment—may exert a greater influence over the fecal resistome and microbiome of feedlot cattle than common metaphylactic antimicrobial drug treatment.

Treatment of food-producing animals with antimicrobial drugs (AMD) is controversial because of concerns regarding promotion of antimicrobial resistance (AMR). To investigate this concern, resistance genes in metagenomic bovine fecal samples during a clinical trial were analyzed to assess the impacts of treatment on beef feedlot cattle resistomes. Four groups of cattle were exposed, using a 2-by-2 factorial design, to different regimens of antimicrobial treatment. Injections of ceftiofur crystalline-free acid (a third-generation cephalosporin) were used to treat all cattle in treatment pens or only a single animal, and either chlortetracycline was included in the feed of all cattle in a pen or the feed was untreated. On days 0 and 26, respectively, pre- and posttrial fecal samples were collected, and resistance genes were characterized using shotgun metagenomics. Treatment with ceftiofur was not associated with changes to β-lactam resistance genes. However, cattle fed chlortetracycline had a significant increase in relative abundance of tetracycline resistance genes. There was also an increase of an AMR class not administered during the study, which is a possible indicator of coselection of resistance genes. Samples analyzed in this study had previously been evaluated by culture characterization ( Escherichia coli and Salmonella ) and quantitative PCR (qPCR) of metagenomic fecal DNA, which allowed comparison of results with this study. In the majority of samples, genes that were selectively enriched through culture and qPCR were not identified through shotgun metagenomic sequencing in this study, suggesting that changes previously documented did not reflect changes affecting the majority of bacterial genetic elements found in the predominant fecal resistome. IMPORTANCE Despite significant concerns about public health implications of AMR in relation to use of AMD in food animals, there are many unknowns about the long- and short-term impact of common uses of AMD for treatment, control, and prevention of disease. Additionally, questions commonly arise regarding how to best measure and quantify AMR genes in relation to public health risks and how to determine which genes are most important. These data provide an introductory view of the utility of using shotgun metagenomic sequencing data as an outcome for clinical trials evaluating the impact of using AMD in food animals.