Temporal dynamics of the fecal microbiome in female pigs from early life through estrus, parturition, and weaning of the first litter of piglets

BACKGROUND

Age-associated changes in the gastrointestinal microbiome of young pigs have been robustly described; however, the temporal dynamics of the fecal microbiome of the female pig from early life to first parity are not well understood. Our objective was to describe microbiome and antimicrobial resistance dynamics of the fecal microbiome of breeding sows from early life through estrus, parturition and weaning of the first litter of piglets (i.e., from 3 to 53 weeks of age).

RESULTS

Our analysis revealed that fecal bacterial populations in developing gilts undergo changes consistent with major maturation milestones. As the pigs progressed towards first estrus, the fecal bacteriome shifted from Rikenellaceae RC9 gut group- and UCG-002-dominated enterotypes to Treponema- and Clostridium sensu stricto 1-dominated enterotypes. After first estrus, the fecal bacteriome stabilized, with minimal changes in enterotype transition and associated microbial diversity from estrus to parturition and subsequent weaning of first litter piglets. Unlike bacterial communities, fecal fungal communities exhibited low diversity with high inter- and intra-pig variability and an increased relative abundance of certain taxa at parturition, including Candida spp. Counts of resistant fecal bacteria also fluctuated over time, and were highest in early life and subsequently abated as the pigs progressed to adulthood.

CONCLUSIONS

This study provides insights into how the fecal microbial community and antimicrobial resistance in female pigs change from three weeks of age throughout their first breeding lifetime. The fecal bacteriome enterotypes and diversity are found to be age-driven and established by the time of first estrus, with minimal changes observed during subsequent physiological stages, such as parturition and lactation, when compared to the earlier age-related shifts. The use of pigs as a model for humans is well-established, however, further studies are needed to understand how our results compare to the human microbiome dynamics. Our findings suggest that the fecal microbiome exhibited consistent changes across individual pigs and became more diverse with age, which is a beneficial characteristic for an animal model system.

Effect of castration timing and weaning strategy on the taxonomic and functional profile of ruminal bacteria and archaea of beef calves

BACKGROUND

Beef cattle experience several management challenges across their lifecycle. Castration and weaning, two major interventions in the early life of beef cattle, can have a substantial impact on animal performance. Despite the key role of the rumen microbiome on productive traits of beef cattle, the effect of castration timing and weaning strategy on this microbial community has not been formally described. We assessed the effect of four castration time windows (at birth, turnout, pre-weaning and weaning) and two weaning strategies (fence-line and truck transportation) on the rumen microbiome in a randomized controlled study with 32 male calves across 3 collection days (i.e., time points). Ruminal fluid samples were submitted to shotgun metagenomic sequencing and changes in the taxonomic (microbiota) and functional profile (metagenome) of the rumen microbiome were described.

RESULTS

Using a comprehensive yet stringent taxonomic classification approach, we identified 10,238 unique taxa classified under 40 bacterial and 7 archaeal phyla across all samples. Castration timing had a limited long-term impact on the rumen microbiota and was not associated with changes in alpha and beta diversity. The interaction of collection day and weaning strategy was associated with changes in the rumen microbiota, which experienced a significant decrease in alpha diversity and shifts in beta diversity within 48 h post-weaning, especially in calves abruptly weaned by truck transportation. Calves weaned using a fence-line weaning strategy had lower relative abundance of Bacteroides, Lachnospira, Fibrobacter and Ruminococcus genera compared to calves weaned by truck transportation. Some genes involved in the hydrogenotrophic methanogenesis pathway (fwdB and fwdF) had higher relative abundance in fence-line-weaned calves post-weaning. The antimicrobial resistance gene tetW consistently represented more than 50% of the resistome across time, weaning and castration groups, without significant changes in relative abundance.

CONCLUSIONS

Within the context of this study, castration timing had limited long-term effects on the rumen microbiota, while weaning strategy had short-term effects on the rumen microbiota and methane-associated metagenome, but not on the rumen resistome.

Age influences the temporal dynamics of microbiome and antimicrobial resistance genes among fecal bacteria in a cohort of production pigs

BACKGROUND

The pig gastrointestinal tract hosts a diverse microbiome, which can serve to select and maintain a reservoir of antimicrobial resistance genes (ARG). Studies suggest that the types and quantities of antimicrobial resistance (AMR) in fecal bacteria change as the animal host ages, yet the temporal dynamics of AMR within communities of bacteria in pigs during a full production cycle remains largely unstudied.

RESULTS

A longitudinal study was performed to evaluate the dynamics of fecal microbiome and AMR in a cohort of pigs during a production cycle; from birth to market age. Our data showed that piglet fecal microbial communities assemble rapidly after birth and become more diverse with age. Individual piglet fecal microbiomes progressed along similar trajectories with age-specific community types/enterotypes and showed a clear shift from E. coli/Shigella-, Fusobacteria-, Bacteroides-dominant enterotypes to Prevotella-, Megaspheara-, and Lactobacillus-dominated enterotypes with aging. Even when the fecal microbiome was the least diverse, the richness of ARGs, quantities of AMR gene copies, and counts of AMR fecal bacteria were highest in piglets at 2 days of age; subsequently, these declined over time, likely due to age-related competitive changes in the underlying microbiome. ARGs conferring resistance to metals and multi-compound/biocides were detected predominately at the earliest sampled ages.

CONCLUSIONS

The fecal microbiome and resistome-along with evaluated descriptors of phenotypic antimicrobial susceptibility of fecal bacteria-among a cohort of pigs, demonstrated opposing trajectories in diversity primarily driven by the aging of pigs.

Evaluating the effects of antimicrobial drug use on the ecology of antimicrobial resistance and microbial community structure in beef feedlot cattle

Introduction

Use of antimicrobial drugs (AMDs) in food producing animals has received increasing scrutiny because of concerns about antimicrobial resistance (AMR) that might affect consumers. Previously, investigations regarding AMR have focused largely on phenotypes of selected pathogens and indicator bacteria, such as Salmonella enterica or Escherichia coli. However, genes conferring AMR are known to be distributed and shared throughout microbial communities. The objectives of this study were to employ target-enriched metagenomic sequencing and 16S rRNA gene amplicon sequencing to investigate the effects of AMD use, in the context of other management and environmental factors, on the resistome and microbiome in beef feedlot cattle.

Methods

This study leveraged samples collected during a previous longitudinal study of cattle at beef feedlots in Canada. This included fecal samples collected from randomly selected individual cattle, as well as composite-fecal samples from randomly selected pens of cattle. All AMD use was recorded and characterized across different drug classes using animal defined daily dose (ADD) metrics.

Results

Overall, fecal resistome composition was dominated by genes conferring resistance to tetracycline and macrolide-lincosamide-streptogramin (MLS) drug classes. The diversity of bacterial phyla was greater early in the feeding period and decreased over time in the feedlot. This decrease in diversity occurred concurrently as the microbiome represented in different individuals and different pens shifted toward a similar composition dominated by Proteobacteria and Firmicutes. Some antimicrobial drug exposures in individuals and groups were associated with explaining a statistically significant proportion of the variance in the resistome, but the amount of variance explained by these important factors was very small (<0.6% variance each), and smaller than associations with other factors measured in this study such as time and feedlot ID. Time in the feedlot was associated with greater changes in the resistome for both individual animals and composite pen-floor samples, although the proportion of the variance associated with this factor was small (2.4% and 1.2%, respectively).

Discussion

Results of this study are consistent with other investigations showing that, compared to other factors, AMD exposures did not have strong effects on antimicrobial resistance or the fecal microbial ecology of beef cattle.

ggMOB: Elucidation of genomic conjugative features and associated cargo genes across bacterial genera using genus-genus mobilization networks

Horizontal gene transfer mediated by conjugation is considered an important evolutionary mechanism of bacteria. It allows organisms to quickly evolve new phenotypic properties including antimicrobial resistance (AMR) and virulence. The frequency of conjugation-mediated cargo gene exchange has not yet been comprehensively studied within and between bacterial taxa. We developed a frequency-based network of genus-genus conjugation features and candidate cargo genes from whole-genome sequence data of over 180,000 bacterial genomes, representing 1,345 genera. Using our method, which we refer to as ggMOB, we revealed that over half of the bacterial genomes contained one or more known conjugation features that matched exactly to at least one other genome. Moreover, the proportion of genomes containing these conjugation features varied substantially by genus and conjugation feature. These results and the genus-level network structure can be viewed interactively in the ggMOB interface, which allows for user-defined filtering of conjugation features and candidate cargo genes. Using the network data, we observed that the ratio of AMR gene representation in conjugative versus non-conjugative genomes exceeded 5:1, confirming that conjugation is a critical force for AMR spread across genera. Finally, we demonstrated that clustering genomes by conjugation profile sometimes correlated well with classical phylogenetic structuring; but that in some cases the clustering was highly discordant, suggesting that the importance of the accessory genome in driving bacterial evolution may be highly variable across both time and taxonomy. These results can advance scientific understanding of bacterial evolution, and can be used as a starting point for probing genus-genus gene exchange within complex microbial communities that include unculturable bacteria. ggMOB is publicly available under the GNU licence at https://ruiz-hci-lab.github.io/ggMOB/.

Target-enriched long-read sequencing (TELSeq) contextualizes antimicrobial resistance genes in metagenomes

BACKGROUND

Metagenomic data can be used to profile high-importance genes within microbiomes. However, current metagenomic workflows produce data that suffer from low sensitivity and an inability to accurately reconstruct partial or full genomes, particularly those in low abundance. These limitations preclude colocalization analysis, i.e., characterizing the genomic context of genes and functions within a metagenomic sample. Genomic context is especially crucial for functions associated with horizontal gene transfer (HGT) via mobile genetic elements (MGEs), for example antimicrobial resistance (AMR). To overcome this current limitation of metagenomics, we present a method for comprehensive and accurate reconstruction of antimicrobial resistance genes (ARGs) and MGEs from metagenomic DNA, termed target-enriched long-read sequencing (TELSeq).

RESULTS

Using technical replicates of diverse sample types, we compared TELSeq performance to that of non-enriched PacBio and short-read Illumina sequencing. TELSeq achieved much higher ARG recovery (>1,000-fold) and sensitivity than the other methods across diverse metagenomes, revealing an extensive resistome profile comprising many low-abundance ARGs, including some with public health importance. Using the long reads generated by TELSeq, we identified numerous MGEs and cargo genes flanking the low-abundance ARGs, indicating that these ARGs could be transferred across bacterial taxa via HGT.

CONCLUSIONS

TELSeq can provide a nuanced view of the genomic context of microbial resistomes and thus has wide-ranging applications in public, animal, and human health, as well as environmental surveillance and monitoring of AMR. Thus, this technique represents a fundamental advancement for microbiome research and application. Video abstract.

Syotti: scalable bait design for DNA enrichment

MOTIVATION

Bait enrichment is a protocol that is becoming increasingly ubiquitous as it has been shown to successfully amplify regions of interest in metagenomic samples. In this method, a set of synthetic probes ('baits') are designed, manufactured and applied to fragmented metagenomic DNA. The probes bind to the fragmented DNA and any unbound DNA is rinsed away, leaving the bound fragments to be amplified for sequencing. Metsky et al. demonstrated that bait-enrichment is capable of detecting a large number of human viral pathogens within metagenomic samples.

RESULTS

We formalize the problem of designing baits by defining the Minimum Bait Cover problem, show that the problem is NP-hard even under very restrictive assumptions, and design an efficient heuristic that takes advantage of succinct data structures. We refer to our method as Syotti. The running time of Syotti shows linear scaling in practice, running at least an order of magnitude faster than state-of-the-art methods, including the method of Metsky et al. At the same time, our method produces bait sets that are smaller than the ones produced by the competing methods, while also leaving fewer positions uncovered. Lastly, we show that Syotti requires only 25 min to design baits for a dataset comprised of 3 billion nucleotides from 1000 related bacterial substrains, whereas the method of Metsky et al. shows clearly super-linear running time and fails to process even a subset of 17% of the data in 72 h.

AVAILABILITY AND IMPLEMENTATION

https://github.com/jnalanko/syotti.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evaluation of Contamination in Milk Samples Pooled From Independently Collected Quarters Within a Laboratory Setting

The primary objective of this observational study was to evaluate the prevalence of contamination from independently collected quarter-level milk samples pooled in a laboratory and subjected to bacterial culture. To address this objective, weekly quarter-level milk samples were collected longitudinally from a cohort of 503 primiparous cows from five organic dairy farms during the first 5 weeks after calving. Individual quarter milk samples were pooled in a laboratory using aseptic technique (“lab-pooled”) and subjected to bacterial culture. In the sample set of 2,006 lab-pooled milk samples, 207 (10.3%) were classified as contaminated using a standard definition (i.e., growth of three or more distinct microorganisms). Subsequent culturing of corresponding quarter-level milk samples revealed that many of the contaminated lab-pooled sample results (i.e., 46.7%) were the result of intramammary infections with different pathogens across the quarters, rather than actual contamination within any single quarter (i.e., “true contamination”). The odds of true contamination were lower when the lab-pooled sample exhibited growth of three microorganisms compared to more than 3 microorganisms. Our findings suggest that pooling of quarter samples within a laboratory setting may yield lower rates of contamination compared to those previously reported from samples composited on-farm, but that current cut-offs to define contamination may need to be evaluated for use with lab-pooled samples. Further investigation of use of lab-pooled samples may be warranted to reduce costs while still providing useful scientific insight.

A Longitudinal Investigation of the Effects of Age, Dietary Fiber Type and Level and Injectable Antimicrobials on the Fecal Microbiome and Antimicrobial Resistance of Finisher Pigs

Age and diet are among the factors that influence the community composition of the fecal microbiome. Additionally, antimicrobial use can alter the composition of bacterial communities. An 86-d study with finisher pigs aimed to evaluate age-related dynamics (d 98-177 of age), effects of types and levels of dietary fiber, and injectable antimicrobials on the fecal microbiome and antimicrobial resistance (AMR) was conducted. A total of 287 pigs, housed in 36 pens, with 7 to 8 pigs per pen, fed a corn grain and soybean meal-based basal diet, formulated to contain 8.7% neutral detergent fiber (NDF), were randomly assigned to one of three treatments: 1. basal diet with no supplement, 2. basal diet supplemented with 20% distillers dried grains with solubles (DDGS) formulated to contain 13.6% NDF, or 3. basal diet supplemented with 14.5% sugar beet pulp (SBP) formulated to contain 13.6% NDF. Five finisher pigs from each treatment group were selected randomly, and fecal samples were collected on d 98, 110, 144, and 177 of age. In addition, fecal samples were collected from pigs that were injected intramuscularly ceftiofur hydrochloride or penicillin G on d 1 and 3 along with pen-mate untreated controls on d 1. Fecal samples were subjected to 16S rRNA amplicon-based microbiome analysis and culture methods to quantify the abundance of total and AMR coliforms and enterococci populations. The alpha diversity, such as species richness, increased with age, and the overall bacterial composition changed with age (P =0.001) and diet (P = 0.001). Diet-associated shifts in the specific bacterial taxa were observed. The richness, diversity, and evenness of bacterial taxa did not differ between pigs that were injected with ceftiofur versus their untreated pen mates or by dietary treatments, but differed in pigs that received penicillin G injection. Both antimicrobial treatments contributed to changes in the overall fecal bacterial composition at the genus level. Collectively, the data demonstrate that both age and the diet (control vs. DDGS-, control vs. SBP- or DDGS- vs. SBP-based diets) were associated with overall bacterial community composition and the impact of age on variations in fecal microbiome composition was greater than the diet. Antibiotic treatment had minimal effect on bacterial diversity and relative abundance of taxa. Further, diets and antimicrobial treatment had minimal impact on the overall counts of AMR coliforms and enterococci populations in feces.

The microbiome of common bedding materials before and after use on commercial dairy farms

Background

Bovine mastitis is one of the most economically important diseases affecting dairy cows. The choice of bedding material has been identified as an important risk factor contributing to the development of mastitis. However, few reports examine both the culturable and nonculturable microbial composition of commonly used bedding materials, i.e., the microbiome. Given the prevalence of nonculturable microbes in most environments, this information could be an important step to understanding whether and how the bedding microbiome acts as a risk factor for mastitis. Therefore, our objective was to characterize the microbiome composition and diversity of bedding material microbiomes, before and after use.

Methods

We collected 88 bedding samples from 44 dairy farms in the U.S. Unused (from storage pile) and used (out of stalls) bedding materials were collected from four bedding types: new sand (NSA), recycled manure solids (RMS), organic non-manure (ON) and recycled sand (RSA). Samples were analyzed using 16S rRNA sequencing of the V3–V4 region.

Results

The overall composition as well as the counts of several microbial taxa differed between bedding types, with Proteobacteria, Actinobacteria, Bacteroidetes and Firmicutes dominating across all types. Used bedding contained a significantly different microbial composition than unused bedding, but the magnitude of this difference varied by bedding type, with RMS bedding exhibiting the smallest difference. In addition, positive correlations were observed between 16S rRNA sequence counts of potential mastitis pathogens (bacterial genera) and corresponding bedding bacterial culture data.

Conclusion

Our results strengthen the role of bedding as a potential source of mastitis pathogens. The consistent shift in the microbiome of all bedding types that occurred during use by dairy cows deserves further investigation to understand whether this shift promotes pathogen colonization and/or persistence, or whether it can differentially impact udder health outcomes. Future studies of bedding and udder health may be strengthened by including a microbiome component to the study design.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-022-00171-2.

Evaluation of the Impact of Antimicrobial Use Protocols in Porcine Reproductive and Respiratory Syndrome Virus-Infected Swine on Phenotypic Antimicrobial Resistance Patterns

A longitudinal study was conducted to assess the impact of different antimicrobial exposures of nursery-phase pigs on patterns of phenotypic antimicrobial resistance (AMR) in fecal indicator organisms throughout the growing phase. Based on practical approaches used to treat moderate to severe porcine reproductive and respiratory syndrome virus (PRRSV)-associated secondary bacterial infections, two antimicrobial protocols of differing intensities of exposure [44.1 and 181.5 animal-treatment days per 1000 animal days at risk (ATD)] were compared with a control group with minimal antimicrobial exposure (2.1 ATD). Litter-matched pigs (n = 108) with no prior antimicrobial exposure were assigned randomly to the treatment groups. Pen fecal samples were collected nine times during the wean-to-finish period and cultured for Escherichia coli and Enterococcus spp. Antimicrobial-susceptibility testing was conducted using NARMS Gram-negative and Gram-positive antibiotic panels. Despite up to 65-fold difference in ATD, few and modest differences were observed between groups and over time. Resistance patterns at marketing overall remained similar to those observed at weaning, prior to any antimicrobial exposures. Those differences observed could not readily be reconciled with the patterns of antimicrobial exposure. Resistance of E. coli to streptomycin was higher in the group exposed to 44.1 ATD, but no aminoglycosides were used. In all instances where resistances differed between time points, the higher resistance occurred early in the trial prior to any antimicrobial exposures. These minimal impacts on AMR despite substantially different antimicrobial exposures point to the lack of understanding of the drivers of AMR at the population level and the likely importance of factors other than antimicrobial exposure. IMPORTANCE Despite a recognized need for more longitudinal studies to assess the effects of antimicrobial use on resistance in food animals, they remain sparse in the literature, and most longitudinal studies of pigs have been observational. The current experimental study had the advantages of greater control of potential confounding, precise measurement of antimicrobial exposures which differed markedly between groups and tracking of pigs until market age. Overall, resistance patterns were remarkably stable between the treatment groups over time, and the differences observed could not be readily reconciled with the antimicrobial exposures, indicating the likely importance of other determinants of AMR at the population level.

Beyond Antimicrobial Use: A Framework for Prioritizing Antimicrobial Resistance Interventions

Antimicrobial resistance (AMR) is a threat to animal and human health. Antimicrobial use has been identified as a major driver of AMR, and reductions in use are a focal point of interventions to reduce resistance. Accordingly, stakeholders in human health and livestock production have implemented antimicrobial stewardship programs aimed at reducing use. Thus far, these efforts have yielded variable impacts on AMR. Furthermore, scientific advances are prompting an expansion and more nuanced appreciation of the many nonantibiotic factors that drive AMR, as well as how these factors vary across systems, geographies, and contexts. Given these trends, we propose a framework to prioritize AMR interventions. We use this framework to evaluate the impact of interventions that focus on antimicrobial use. We conclude by suggesting that priorities be expanded to include greater consideration of host-microbial interactions that dictate AMR, as well as anthropogenic and environmental systems that promote dissemination of AMR.

Metagenome SNP calling via read-colored de Bruijn graphs

MOTIVATION

Metagenomics refers to the study of complex samples containing of genetic contents of multiple individual organisms and, thus, has been used to elucidate the microbiome and resistome of a complex sample. The microbiome refers to all microbial organisms in a sample, and the resistome refers to all of the antimicrobial resistance (AMR) genes in pathogenic and non-pathogenic bacteria. Single-nucleotide polymorphisms (SNPs) can be effectively used to 'fingerprint' specific organisms and genes within the microbiome and resistome and trace their movement across various samples. However, to effectively use these SNPs for this traceability, a scalable and accurate metagenomics SNP caller is needed. Moreover, such an SNP caller should not be reliant on reference genomes since 95% of microbial species is unculturable, making the determination of a reference genome extremely challenging. In this article, we address this need.

RESULTS

We present LueVari, a reference-free SNP caller based on the read-colored de Bruijn graph, an extension of the traditional de Bruijn graph that allows repeated regions longer than the k-mer length and shorter than the read length to be identified unambiguously. LueVari is able to identify SNPs in both AMR genes and chromosomal DNA from shotgun metagenomics data with reliable sensitivity (between 91% and 99%) and precision (between 71% and 99%) as the performance of competing methods varies widely. Furthermore, we show that LueVari constructs sequences containing the variation, which span up to 97.8% of genes in datasets, which can be helpful in detecting distinct AMR genes in large metagenomic datasets.

AVAILABILITY AND IMPLEMENTATION

Code and datasets are publicly available at https://github.com/baharpan/cosmo/tree/LueVari.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Effect of subclinical mastitis detected in the first month of lactation on somatic cell count linear scores, milk yield, fertility, and culling of dairy cows in certified organic herds

It is well established that subclinical mastitis (SCM), characterized by somatic cell count (SCC) >200,000 cells/mL, has a negative effect on the productivity, reproductive performance, and survivability of cows from conventional dairy herds. However, in organic herds, where the use of antimicrobial drugs is restricted for the treatment and control of intramammary infections (IMI) in dairy cows, little is known about the effect of SCM on performance and survivability. The objective of this study was to evaluate whether SCM diagnosed during the first month of lactation was associated with SCC linear score dynamics, milk production, fertility, and culling of dairy cows in USDA-certified organic herds. We collected data from 2 organic herds in New Mexico and Texas. A total of 1,511 cows that calved between June 2018 and May 2019 were included in the study and were followed until month 10 of the current lactation. Cows with SCC >200,000 cells/mL in the first month of lactation were considered to have SCM. We used mixed linear regression models accounting for repeated measures to assess the effect of SCM on monthly milk production and SCC linear scores. We used Cox proportional hazards models to evaluate the effect of SCM on the risk of pregnancy and culling. We considered parity, farm, previous gestation length, stillbirth, twinning, dystocia, and 2- and 3-way interactions as potential confounders. Cows diagnosed with SCM during the first month of lactation produced less milk than cows without SCM. Cows with SCM had elevated SCC linear scores during their previous lactation and throughout the subsequent months of lactation compared to cows without SCM. The effect of SCM on SCC linear scores was more pronounced in multiparous than primiparous cows. Subclinical mastitis during the first month of lactation did not affect the likelihood of pregnancy during the first 300 d in milk. Cows with SCM in the first month were more likely to die or be culled during the 300 d of lactation than cows without SCM. We observed that elevated SCC in the first month of lactation had detrimental effects on the milk yield and survivability of dairy cows in USDA organic herds, but it did not affect reproductive performance. We demonstrated that cows with SCM diagnosed in the first month of lactation continued to have elevated SCC linear scores throughout their entire lactation, and that elevated SCC was carried over from the previous lactation.

Investigating the cow skin and teat canal microbiomes of the bovine udder using different sampling and sequencing approaches

There is a need for standardized, efficient, and practical sampling methods to support large population-based studies of the internal and external epithelial microbiomes of the bovine udder. The primary objective of this study was to evaluate different sampling devices for the isolation of microbial DNA originating from the internal and external teat epithelium. Secondary objectives were to survey and compare the microbial diversity of external and teat canal epithelial microbiomes using amplicon and shotgun metagenomic sequencing approaches. To address these objectives, we enrolled a convenience sample of 24 Holstein dairy cows and collected samples from the external epithelium at the base of udder, the external teat barrel epithelium, the external teat apex epithelium, and the teat canal epithelium. Extracted DNA was quantified and subjected to PCR amplification of the V4 hypervariable region of the 16S rRNA gene and sequenced on the Illumina MiSeq platform (Illumina Inc., San Diego, CA). A subset of samples was subjected to a shallow shotgun metagenomic assay on the Illumina HiSeq platform. For samples collected from the external teat epithelium, we found that gauze squares consistently yielded more DNA than swabs, and Simpson's reciprocal index of diversity was higher for gauze than for swabs. The teat canal epithelial samples exhibited significantly lower diversity than the external sampling locations, but there were no significant differences in diversity between teat apex, teat barrel, and base of the udder samples. There were, however, differences in the microbial distribution and abundances of specific bacteria across external epithelial surfaces. The proportion of shotgun sequence reads classified as Bos taurus was highly variable between sampling locations, ranging from 0.33% in teat apex samples to 99.91% in teat canal samples. These results indicate that gauze squares should be considered for studying the microbiome of the external epithelium of the bovine udder, particularly if DNA yield must be maximized. Further, the relative proportion of host to non-host DNA present in samples collected from the internal and external teat epithelium should be considered when designing studies that utilize shotgun metagenomic sequencing.

Mobilization of Antibiotic Resistance: Are Current Approaches for Colocalizing Resistomes and Mobilomes Useful?

Antimicrobial resistance (AMR) poses a global human and animal health threat, and predicting AMR persistence and transmission remains an intractable challenge. Shotgun metagenomic sequencing can help overcome this by enabling characterization of AMR genes within all bacterial taxa, most of which are uncultivatable in laboratory settings. Shotgun sequencing, therefore, provides a more comprehensive glance at AMR "potential" within samples, i.e., the "resistome." However, the risk inherent within a given resistome is predicated on the genomic context of various AMR genes, including their presence within mobile genetic elements (MGEs). Therefore, resistome risk stratification can be advanced if AMR profiles are considered in light of the flanking mobilizable genomic milieu (e.g., plasmids, integrative conjugative elements (ICEs), phages, and other MGEs). Because such mediators of horizontal gene transfer (HGT) are involved in uptake by pathogens, investigators are increasingly interested in characterizing that resistome fraction in genomic proximity to HGT mediators, i.e., the "mobilome"; we term this "colocalization." We explored the utility of common colocalization approaches using alignment- and assembly-based techniques, on clinical (human) and agricultural (cattle) fecal metagenomes, obtained from antimicrobial use trials. Ordination revealed that tulathromycin-treated cattle experienced a shift in ICE and plasmid composition versus untreated animals, though the resistome was unaffected during the monitoring period. Contrarily, the human resistome and mobilome composition both shifted shortly after antimicrobial administration, though this rebounded to pre-treatment status. Bayesian networks revealed statistical AMR-MGE co-occurrence in 19 and 2% of edges from the cattle and human networks, respectively, suggesting a putatively greater mobility potential of AMR in cattle feces. Conversely, using Mobility Index (MI) and overlap analysis, abundance of de novo-assembled contigs supporting resistomes flanked by MGE increased shortly post-exposure within human metagenomes, though > 40 days after peak dose such contigs were rare (∼2%). MI was not substantially altered by antimicrobial exposure across all cattle metagenomes, ranging 0.5-4.0%. We highlight that current alignment- and assembly-based methods estimating resistome mobility yield contradictory and incomplete results, likely constrained by approach-specific data inputs, and bioinformatic limitations. We discuss recent laboratory and computational advancements that may enhance resistome risk analysis in clinical, regulatory, and commercial applications.

MEGARes 2.0: a database for classification of antimicrobial drug, biocide and metal resistance determinants in metagenomic sequence data

Antimicrobial resistance (AMR) is a threat to global public health and the identification of genetic determinants of AMR is a critical component to epidemiological investigations. High-throughput sequencing (HTS) provides opportunities for investigation of AMR across all microbial genomes in a sample (i.e. the metagenome). Previously, we presented MEGARes, a hand-curated AMR database and annotation structure developed to facilitate the analysis of AMR within metagenomic samples (i.e. the resistome). Along with MEGARes, we released AmrPlusPlus, a bioinformatics pipeline that interfaces with MEGARes to identify and quantify AMR gene accessions contained within a metagenomic sequence dataset. Here, we present MEGARes 2.0 (https://megares.meglab.org), which incorporates previously published resistance sequences for antimicrobial drugs, while also expanding to include published sequences for metal and biocide resistance determinants. In MEGARes 2.0, the nodes of the acyclic hierarchical ontology include four antimicrobial compound types, 57 classes, 220 mechanisms of resistance, and 1,345 gene groups that classify the 7,868 accessions. In addition, we present an updated version of AmrPlusPlus (AMR ++ version 2.0), which improves accuracy of classifications, as well as expanding scalability and usability.

MEGARes 2.0: a database for classification of antimicrobial drug, biocide and metal resistance determinants in metagenomic sequence data

Antimicrobial resistance (AMR) is a threat to global public health and the identification of genetic determinants of AMR is a critical component to epidemiological investigations. High-throughput sequencing (HTS) provides opportunities for investigation of AMR across all microbial genomes in a sample (i.e. the metagenome). Previously, we presented MEGARes, a hand-curated AMR database and annotation structure developed to facilitate the analysis of AMR within metagenomic samples (i.e. the resistome). Along with MEGARes, we released AmrPlusPlus, a bioinformatics pipeline that interfaces with MEGARes to identify and quantify AMR gene accessions contained within a metagenomic sequence dataset. Here, we present MEGARes 2.0 (https://megares.meglab.org), which incorporates previously published resistance sequences for antimicrobial drugs, while also expanding to include published sequences for metal and biocide resistance determinants. In MEGARes 2.0, the nodes of the acyclic hierarchical ontology include four antimicrobial compound types, 57 classes, 220 mechanisms of resistance, and 1,345 gene groups that classify the 7,868 accessions. In addition, we present an updated version of AmrPlusPlus (AMR ++ version 2.0), which improves accuracy of classifications, as well as expanding scalability and usability.

A Cautionary Report for Pathogen Identification Using Shotgun Metagenomics; A Comparison to Aerobic Culture and Polymerase Chain Reaction for Salmonella enterica Identification

This study was conducted to compare aerobic culture, polymerase chain reaction (PCR), lateral flow immunoassay (LFI), and shotgun metagenomics for identification of Salmonella enterica in feces collected from feedlot cattle. Samples were analyzed in parallel using all four tests. Results from aerobic culture and PCR were 100% concordant and indicated low S. enterica prevalence (3/60 samples positive). Although low S. enterica prevalence restricted formal statistical comparisons, LFI and deep metagenomic sequencing results were discordant with these results. Specifically, metagenomic analysis using k-mer-based classification against the RefSeq database indicated that 11/60 of samples contained sequence reads that matched to the S. enterica genome and uniquely identified this species of bacteria within the sample. However, further examination revealed that plasmid sequences were often included with bacterial genomic sequence data submitted to NCBI, which can lead to incorrect taxonomic classification. To circumvent this classification problem, we separated all plasmid sequences included in bacterial RefSeq genomes and reassigned them to a unique taxon so that they would not be uniquely associated with specific bacterial species such as S. enterica. Using this revised database and taxonomic structure, we found that only 6/60 samples contained sequences specific for S. enterica, suggesting increased relative specificity. Reads identified as S. enterica in these six samples were further evaluated using BLAST and NCBI's nr/nt database, which identified that only 2/60 samples contained reads exclusive to S. enterica chromosomal genomes. These two samples were culture- and PCR-negative, suggesting that even deep metagenomic sequencing suffers from lower sensitivity and specificity in comparison to more traditional pathogen detection methods. Additionally, no sample reads were taxonomically classified as S. enterica with two other metagenomic tools, Metagenomic Intra-species Diversity Analysis System (MIDAS) and Metagenomic Phylogenetic Analysis 2 (MetaPhlAn2). This study re-affirmed that the traditional techniques of aerobic culture and PCR provide similar results for S. enterica identification in cattle feces. On the other hand, metagenomic results are highly influenced by the classification method and reference database employed. These results highlight the nuances of computational detection of species-level sequences within short-read metagenomic sequence data, and emphasize the need for cautious interpretation of such results.

Characterization of the Microbial Resistome in Conventional and "Raised Without Antibiotics" Beef and Dairy Production Systems

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.

Hierarchical Hidden Markov models enable accurate and diverse detection of antimicrobial resistance sequences

The characterization of antimicrobial resistance genes from high-throughput sequencing data has become foundational in public health research and regulation. This requires mapping sequence reads to databases of known antimicrobial resistance genes to determine the genes present in the sample. Mapping sequence reads to known genes is traditionally accomplished using alignment. Alignment methods have high specificity but are limited in their ability to detect sequences that are divergent from the reference database, which can result in a substantial false negative rate. We address this shortcoming through the creation of Meta-MARC, which enables detection of diverse resistance sequences using hierarchical, DNA-based Hidden Markov Models. We first describe Meta-MARC and then demonstrate its efficacy on simulated and functional metagenomic datasets. Meta-MARC has higher sensitivity relative to competing methods. This sensitivity allows for detection of sequences that are divergent from known antimicrobial resistance genes. This functionality is imperative to expanding existing antimicrobial gene databases.

Investigating Effects of Tulathromycin Metaphylaxis on the Fecal Resistome and Microbiome of Commercial Feedlot Cattle Early in the Feeding Period

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.

Enrichment allows identification of diverse, rare elements in metagenomic resistome-virulome sequencing

BACKGROUND

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.

RESULTS

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.

CONCLUSIONS

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.

Succinct colored de Bruijn graphs

MOTIVATION

In 2012, Iqbal et al. introduced the colored de Bruijn graph, a variant of the classic de Bruijn graph, which is aimed at 'detecting and genotyping simple and complex genetic variants in an individual or population'. Because they are intended to be applied to massive population level data, it is essential that the graphs be represented efficiently. Unfortunately, current succinct de Bruijn graph representations are not directly applicable to the colored de Bruijn graph, which requires additional information to be succinctly encoded as well as support for non-standard traversal operations.

RESULTS

Our data structure dramatically reduces the amount of memory required to store and use the colored de Bruijn graph, with some penalty to runtime, allowing it to be applied in much larger and more ambitious sequence projects than was previously possible.

AVAILABILITY AND IMPLEMENTATION

https://github.com/cosmo-team/cosmo/tree/VARI.

CONTACT

martin.muggli@colostate.edu.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

MEGARes: an antimicrobial resistance database for high throughput sequencing

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.

Characterization of the resistome in manure, soil and wastewater from dairy and beef production systems

It has been proposed that livestock production effluents such as wastewater, airborne dust and manure increase the density of antimicrobial resistant bacteria and genes in the environment. The public health risk posed by this proposed outcome has been difficult to quantify using traditional microbiological approaches. We utilized shotgun metagenomics to provide a first description of the resistome of North American dairy and beef production effluents, and identify factors that significantly impact this resistome. We identified 34 mechanisms of antimicrobial drug resistance within 34 soil, manure and wastewater samples from feedlot, ranch and dairy operations. The majority of resistance-associated sequences found in all samples belonged to tetracycline resistance mechanisms. We found that the ranch samples contained significantly fewer resistance mechanisms than dairy and feedlot samples, and that the resistome of dairy operations differed significantly from that of feedlots. The resistome in soil, manure and wastewater differed, suggesting that management of these effluents should be tailored appropriately. By providing a baseline of the cattle production waste resistome, this study represents a solid foundation for future efforts to characterize and quantify the public health risk posed by livestock effluents.

Use of Metagenomic Shotgun Sequencing Technology To Detect Foodborne Pathogens within the Microbiome of the Beef Production Chain

Foodborne illnesses associated with pathogenic bacteria are a global public health and economic challenge. The diversity of microorganisms (pathogenic and nonpathogenic) that exists within the food and meat industries complicates efforts to understand pathogen ecology. Further, little is known about the interaction of pathogens within the microbiome throughout the meat production chain. Here, a metagenomic approach and shotgun sequencing technology were used as tools to detect pathogenic bacteria in environmental samples collected from the same groups of cattle at different longitudinal processing steps of the beef production chain: cattle entry to feedlot, exit from feedlot, cattle transport trucks, abattoir holding pens, and the end of the fabrication system. The log read counts classified as pathogens per million reads for Salmonella enterica,Listeria monocytogenes,Escherichia coli,Staphylococcus aureus, Clostridium spp. (C. botulinum and C. perfringens), and Campylobacter spp. (C. jejuni,C. coli, and C. fetus) decreased over subsequential processing steps. Furthermore, the normalized read counts for S. enterica,E. coli, and C. botulinumwere greater in the final product than at the feedlots, indicating that the proportion of these bacteria increased (the effect on absolute numbers was unknown) within the remaining microbiome. From an ecological perspective, data indicated that shotgun metagenomics can be used to evaluate not only the microbiome but also shifts in pathogen populations during beef production. Nonetheless, there were several challenges in this analysis approach, one of the main ones being the identification of the specific pathogen from which the sequence reads originated, which makes this approach impractical for use in pathogen identification for regulatory and confirmation purposes.

Resistome diversity in cattle and the environment decreases during beef production

Antimicrobial resistant determinants (ARDs) can be transmitted from livestock systems through meat products or environmental effluents. The public health risk posed by these two routes is not well understood, particularly in non-pathogenic bacteria. We collected pooled samples from 8 groups of 1741 commercial cattle as they moved through the process of beef production from feedlot entry through slaughter. We recorded antimicrobial drug exposures and interrogated the resistome at points in production when management procedures could potentially influence ARD abundance and/or transmission. Over 300 unique ARDs were identified. Resistome diversity decreased while cattle were in the feedlot, indicating selective pressure. ARDs were not identified in beef products, suggesting that slaughter interventions may reduce the risk of transmission of ARDs to beef consumers. This report highlights the utility and limitations of metagenomics for assessing public health risks regarding antimicrobial resistance, and demonstrates that environmental pathways may represent a greater risk than the food supply.

Mannheimia haemolytica in feedlot cattle: prevalence of recovery and associations with antimicrobial use, resistance, and health outcomes

BACKGROUND

Mannheimia haemolytica is an important etiological agent in bovine respiratory disease.

OBJECTIVES

Explore risk factors for recovery of susceptible and resistant M. haemolytica in feedlot cattle and explore associations with health outcomes.

ANIMALS

Cattle (n = 5,498) from 4 feedlots sampled at arrival and later in feeding period.

METHODS

Susceptibility of M. haemolytica isolates tested for 21 antimicrobials. Records of antimicrobial use and health events analyzed using multivariable regression.

RESULTS

M. haemolytica recovered from 29% of cattle (1,596/5,498), 13.1% at arrival (95% CI, 12.3-14.1%), and 19.8% at second sampling (95% CI, 18.7-20.9%). Nearly half of study cattle received antimicrobial drugs (AMDs) parenterally, mostly as metaphylactic treatment at arrival. Individual parenteral AMD exposures were associated with decreased recovery of M. haemolytica (OR, 0.2; 95% CI, 0.02-1.2), whereas exposure in penmates was associated with increased recovery (OR, 1.5; 95% CI, 1.05-2.2). Most isolates were pan-susceptible (87.8%; 95% CI, 87.0-89.4%). AMD exposures were not associated with resistance to any single drug. Multiply-resistant isolates were rare (5.9%; 95% CI, 5.1-6.9%), but AMD exposures in pen mates were associated with increased odds of recovering multiply-resistant M. haemolytica (OR, 23.9; 95% CI, 8.4-68.3). Cattle positive for M. haemolytica on arrival were more likely to become ill within 10 days (OR, 1.7; 95% CI, 1.1-2.4).

CONCLUSIONS AND CLINICAL IMPORTANCE

Resistance generally was rare in M. haemolytica. Antimicrobial drug exposures in penmates increased the risk of isolating susceptible and multiply-resistant M. haemolytica, a finding that could be explained by contagious spread.

Antimicrobial Resistance in Escherichia coli Recovered from Feedlot Cattle and Associations with Antimicrobial Use

The objectives of this study were to estimate the prevalence of antimicrobial resistance (AMR) and to investigate the associations between exposures to antimicrobial drugs (AMDs) and AMR in fecal non-type specific Escherichia coli (NTSEC) recovered from a large population of feedlot cattle. Two-stage random sampling was used to select individually identified cattle for enrollment, which were sampled at arrival and then a second time later in the feeding period. Advanced regression techniques were used to estimate resistance prevalences, and to investigate associations between AMD exposures in enrolled cattle and penmates and AMR identified in NTSEC recovered from the second sample set. Resistance was most commonly detected to tetracycline, streptomycin, and sulfisoxazole, and was rarely identified for critically important AMDs. All cattle were exposed to AMDs in feed, and 45% were treated parenterally. While resistance prevalence generally increased during the feeding period, most AMD exposures were not significantly associated with AMR outcomes. Exposures of enrolled cattle to tetracycline were associated with increased resistance to tetracycline and trimethoprim sulfa, while beta-lactam exposures were associated with decreased likelihood of detecting streptomycin resistance. Pen-level AMD exposure measures were not associated with resistance outcomes. These findings suggest that tetracycline treatment of feedlot cattle can be associated with modest increases in risk for recovery of resistant NTSEC, but the numerous treatments with an advanced macrolide (tulathromycin) were not associated with detectable increases in resistance in NTSEC. All cattle were exposed to in-feed treatments of tetracycline and this could limit the ability to identify the full impact of these exposures, but these exposures varied for enrolled cattle varied, providing an opportunity to evaluate a dose response. While AMD exposures were not associated with detectably increased risks for resistance to critically important AMDs, rare resistance outcomes and infrequent exposure to other important AMDs (e.g., cephalosporins) limited our ability to rigorously investigate questions regarding factors that can influence resistance to these important AMDs.