Analysis of the northern pitcher plant (Sarracenia purpurea L.) phytotelm bacteriome throughout a temperate region growing season.

The insectivorous Northern Pitcher Plant, Sarracenia purpurea, recruits a dynamic biotic community in the rainwater collected by its pitcher-shaped leaves. Insect capture and degradation within the pitcher fluid (phytotelma) has been well documented as a mechanism for supplementing the plant's nitrogen, phosphorous, and micronutrient requirements. Metagenomic studies have shown a diverse microbiome in this phytotelm environment, including taxa that contribute metabolically to prey digestion. In this investigation, we used high-throughput 16S rDNA sequencing and bioinformatics to analyze the S. purpurea phytotelm bacteriome as it changes through the growing season (May-September) in plants from the north-central region of the species' native range. Additionally, we used molecular techniques to detect and quantify bacterial nitrogenase genes (nifH) in all phytotelm samples to explore the hypothesis that diazotrophy is an additional mechanism of supplying biologically available nitrogen to S. purpurea. The results of this study indicate that while prokaryote diversity remains relatively stable in plants at different locations within our region, diversity changes significantly as the growing season progresses. Furthermore, nifH genes were detected at biologically significant concentrations in one hundred percent of samples, suggesting that nitrogen fixation may be an important contributor to the S. purpurea nutrient budget.

Slaughtering processes impact microbial communities and antimicrobial resistance genes of pig carcasses.

Several steps in the abattoir can influence the presence of microbes and associated resistance genes (ARGs) on the animal carcasses used for further meat processing. We investigated how these processes influence the resistome-microbiome of groups of pigs with different on-farm antimicrobial exposure status, from the moment they entered the abattoir until the end of carcass processing. Using a targeted enrichment metagenomic approach, we identified 672 unique ARGs conferring resistance to 43 distinct AMR classes from pooled skin (N = 42) and carcass swabs (N = 63) collected sequentially before, during, and after the slaughter process and food safety interventions. We observed significant variations in the resistome and microbial profiles of pigs before and after slaughter, as well as a significant decline in ARG counts, diversity, and microbial DNA load during slaughter and carcass processing, irrespective of prior antimicrobial treatments on the farm. These results suggest that existing interventions in the abattoir are effective in reducing not only the pathogen load but also the overall bacterial burden, including ARGs on pork carcasses. Concomitant with reductions in microbial and ARG counts, we observed an increase in the relative abundance of non-drug-specific ARGs, such as those conferring resistance to metals and biocides, and in particular mercury. Using a strict colocalization procedure, we found that most mercury ARGs were associated with genomes from the Pseudomonadaceae and Enterobacteriaceae families. Collectively, these findings demonstrate that slaughter and processing practices within the abattoir can shape the microbial and ARG profiles of pork carcasses during the transition from living muscle to meat.

Scalable de novo classification of antibiotic resistance of Mycobacterium tuberculosis.

MOTIVATION: World Health Organization estimates that there were over 10 million cases of tuberculosis (TB) worldwide in 2019, resulting in over 1.4 million deaths, with a worrisome increasing trend yearly. The disease is caused by Mycobacterium tuberculosis (MTB) through airborne transmission. Treatment of TB is estimated to be 85% successful, however, this drops to 57% if MTB exhibits multiple antimicrobial resistance (AMR), for which fewer treatment options are available. RESULTS: We develop a robust machine-learning classifier using both linear and nonlinear models (i.e. LASSO logistic regression (LR) and random forests (RF)) to predict the phenotypic resistance of Mycobacterium tuberculosis (MTB) for a broad range of antibiotic drugs. We use data from the CRyPTIC consortium to train our classifier, which consists of whole genome sequencing and antibiotic susceptibility testing (AST) phenotypic data for 13 different antibiotics. To train our model, we assemble the sequence data into genomic contigs, identify all unique 31-mers in the set of contigs, and build a feature matrix M, where M[i, j] is equal to the number of times the ith 31-mer occurs in the jth genome. Due to the size of this feature matrix (over 350 million unique 31-mers), we build and use a sparse matrix representation. Our method, which we refer to as MTB++, leverages compact data structures and iterative methods to allow for the screening of all the 31-mers in the development of both LASSO LR and RF. MTB++ is able to achieve high discrimination (F-1 >80%) for the first-line antibiotics. Moreover, MTB++ had the highest F-1 score in all but three classes and was the most comprehensive since it had an F-1 score >75% in all but four (rare) antibiotic drugs. We use our feature selection to contextualize the 31-mers that are used for the prediction of phenotypic resistance, leading to some insights about sequence similarity to genes in MEGARes. Lastly, we give an estimate of the amount of data that is needed in order to provide accurate predictions. AVAILABILITY: The models and source code are publicly available on Github at https://github.com/M-Serajian/MTB-Pipeline.

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.

Associations between early lactation intramammary infections and udder health and performance during the first 180 days in milk in first-lactation organic dairy cows.

Prior data from our group showed that first-lactation cows under organic management in United States have a high prevalence of Staphylococcus aureus, Streptococcus spp., and Staphylococcus chromogenes intramammary infections (IMI) in early lactation. Nonetheless, the relationship between IMI, udder health, and milk production in organically reared primiparous cows remains elusive. The objectives of this observational study were to investigate the relationship between presence and persistence of IMI in the first 35 d in milk (DIM) and somatic cell count (SCC) and milk production during the first 6 mo of lactation on first-lactation organic dairy cows. The analysis included a total of 1,348 composite milk samples collected during the first 35 DIM that were submitted for milk culture and 1,674 Dairy Herd Improvement Association (DHIA) tests during the first 180 DIM from 333 heifers in 4 organic dairy farms, enrolled between February 2019 and January 2020. The association between IMI in the first 35 DIM and new high SCC (SCC > 200,000 cells/mL) and milk production during the first 6 mo of lactation was investigated using Cox proportional hazards regression and mixed linear regression, respectively. The association between IMI persistence (harboring the same microorganism as reported by the laboratory for 2 or more samples) in the first 35 DIM and number of DHIA tests with high SCC during the first 6 mo of lactation was modeled using negative binomial regression. The presence of IMI by Staph. aureus (hazard ratio [HR] [95% confidence interval {CI}]: 3.35 [2.64, 4.25]) or Streptococcus spp. (HR [95% CI]: 2.25 [2.12, 2.39]) during the first 35 DIM was associated with an increased risk of new high SCC during the first 6 mo of lactation. Milk production was reduced when Streptococcus spp. were identified in milk samples. However, there was no evidence of a difference in milk production in Staph. aureus IMI. Isolation of non-aureus staphylococci and mammaliicocci was related to a mild increase in the hazards of high SCC (HR [95% CI]: 1.34 [0.97, 1.85]) and a decrease in milk production during one or more postpartum tests. Presence of gram-negative or Streptococcus-like organisms IMI was not associated with either high SCC or milk production. Presence of Bacillus IMI was associated with a lower hazard of new high SCC (HR [95% CI]: 0.45 [0.30, 0.68]), and higher milk production during the first 180 d of lactation (overall estimate [95% CI]: 1.7 kg/d [0.3, 3.0]). The persistence of IMI in the first 35 DIM was associated with the number of tests with high SCC during the lactation for all microorganisms except for Staphylococcus chromogenes. Therefore, our results suggest that the persistence of IMI in the first 35 DIM could be an important factor to understand the association between IMI detected in early lactation and lactational SCC and milk production in organic dairy heifers. Our study described associations between IMI, udder health, and milk production in first-lactation organic dairy cows that are consistent with findings from conventional dairy farms.

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.

Evaluating variations in metabolic profiles during the dry period related to the time of hyperketonemia onset in dairy cows.

Hyperketonemia (HYK) in early lactation can have a different impact on health and productivity depending on the timing of HYK onset. While specific metabolites measured during the dry period may serve as biomarkers of HYK, the correlations between metabolites represent a challenge for the use of metabolic profiles dataset, and little has been explored on HYK. This exploratory cohort study aimed a) to characterize the correlations among metabolites measured during the late dry period in dairy cows, and b) to identify biomarkers in the late dry period associated with the onset of HYK at the first (wk1) and second (wk2) week of lactation. Individual blood samples from 440 Holstein dairy cows were collected at 21 ± 3 days before expected parturition. From each sample, 36 different metabolites were measured in serum and plasma. Hyperketonemia was diagnosed in wk1 and wk2 of lactation based on the blood concentration of beta-hydroxybutyrate (BHB > 1.2 mmol/L). Principal component analysis (PCA) was performed to reduce metabolites to a smaller number of uncorrelated components. Multivariable logistic regression models were applied to assess the associations between principal components (PC) and HYK at wk1 only (HYK+ wk1), wk2 only (HYK+ wk2), or both weeks (HYK+ wk1-2). The incidence of HYK was 16.2% in the first week, 13.0% in the second week, and 21.2% within the first two weeks of lactation. The results of PCA highlighted 10 PCs from which two were associated with HYK+ wk1 as compared with cows without HYK during the first two weeks of lactation (non-HYK); the PC a2 led by bilirubin and non-esterified fatty acids (OR = 1.29; 95%CI: 1.02-1.68), and the PC a5 led by alkaline phosphatase (ALP) and gamma-glutamyl transferase (GGT) (OR = 2.77; 95%CI: 1.61-4.97). There was no evidence of an association between any PC and HYK+ wk2 (vs. non-HYK cows). Cows with elevated PC a5 (led by ALP and GGT) in the dry period were 3.18 times more likely to be HYK+ wk1 than HYK+ wk2 (OR: 3.18, 95%CI: 1.34-8.73; P = 0.013). Overall, the main hypothesis generated by our exploratory study suggests that cows with biomarkers of liver dysfunction (ALP, GGT, bilirubin) assessed by PCA at 3 weeks before calving are more likely to develop HYK during the first week of lactation compared to the second week. In addition, results suggest that cows with HYK in both of the first two weeks of lactation had an overall metabolic disbalance during the onset of the late dry period, which based on PCs, encompass biomarkers related to glucogenic and ketogenic metabolic pathways as well as liver dysfunction and fatty liver. Further research is needed to determine the underlying mechanisms associated with the different adaptations between cows that develop HYK during the first and second week of lactation.

Investigation of intramammary infections in primiparous cows during early lactation on organic dairy farms.

Previous studies have shown that organically raised dairy cows have an increased prevalence of Staphylococcus aureus compared with conventionally raised dairy cows. However, little information exists about the dynamics of intramammary infection (IMI) in primiparous cows during early lactation on organic dairy farms. The objective of this study was to describe the IMI dynamics of primiparous cows on certified organic farms during early lactation. This longitudinal study enrolled 503 primiparous cows from 5 organic dairy farms from February 2019 to January 2020. Quarter-level milk samples were collected aseptically on a weekly basis during the first 5 wk of lactation. Samples were pooled by cow and time point into composite samples inside a sterilized laminar hood and submitted for microbiological culture. For each of the different microorganisms identified, we estimated the prevalence in each postpartum sample, period prevalence (PP), cumulative incidence, and persistence of IMI. Logistic regression models were used to investigate whether the prevalence of IMI differed by farm or sampling time points and whether IMI persistence differed between detected microorganisms. Our findings revealed a high prevalence of Staphylococcus aureus (PP = 18.9%), non-aureus staphylococci and closely related mammaliicoccal species (PP = 52.1%), and Streptococcus spp. and Streptococcus-like organisms (PP = 32.1%) within the study population. The prevalence of these microorganisms varied significantly between farms. Staphylococcus aureus and Staphylococcus chromogenes exhibited significantly higher IMI persistence compared with other detected bacterial taxa, confirming the divergent epidemiological behavior in terms of IMI chronicity across different microorganisms. This study improves our understanding of the epidemiology of mastitis-causing pathogens in organically raised primiparous cows, which can be used to tailor mastitis control plans for this unique yet growing subpopulation of dairy cows.

The K-mer antibiotic resistance gene variant analyzer (KARGVA).

Characterization of antibiotic resistance genes (ARGs) from high-throughput sequencing data of metagenomics and cultured bacterial samples is a challenging task, with the need to account for both computational (e.g., string algorithms) and biological (e.g., gene transfers, rearrangements) aspects. Curated ARG databases exist together with assorted ARG classification approaches (e.g., database alignment, machine learning). Besides ARGs that naturally occur in bacterial strains or are acquired through mobile elements, there are chromosomal genes that can render a bacterium resistant to antibiotics through point mutations, i.e., ARG variants (ARGVs). While ARG repositories also collect ARGVs, there are only a few tools that are able to identify ARGVs from metagenomics and high throughput sequencing data, with a number of limitations (e.g., pre-assembly, a posteriori verification of mutations, or specification of species). In this work we present the k-mer, i.e., strings of fixed length k, ARGV analyzer - KARGVA - an open-source, multi-platform tool that provides: (i) an ad hoc, large ARGV database derived from multiple sources; (ii) input capability for various types of high-throughput sequencing data; (iii) a three-way, hash-based, k-mer search setup to process data efficiently, linking k-mers to ARGVs, k-mers to point mutations, and ARGVs to k-mers, respectively; (iv) a statistical filter on sequence classification to reduce type I and II errors. On semi-synthetic data, KARGVA provides very high accuracy even in presence of high sequencing errors or mutations (99.2 and 86.6% accuracy within 1 and 5% base change rates, respectively), and genome rearrangements (98.2% accuracy), with robust performance on ad hoc false positive sets. On data from the worldwide MetaSUB consortium, comprising 3,700+ metagenomics experiments, KARGVA identifies more ARGVs than Resistance Gene Identifier (4.8x) and PointFinder (6.8x), yet all predictions are below the expected false positive estimates. The prevalence of ARGVs is correlated to ARGs but ecological characteristics do not explain well ARGV variance. KARGVA is publicly available at https://github.com/DataIntellSystLab/KARGVA under MIT license.

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.

MEGARes and AMR++, v3.0: an updated comprehensive database of antimicrobial resistance determinants and an improved software pipeline for classification using high-throughput sequencing.

Antimicrobial resistance (AMR) is considered a critical threat to public health, and genomic/metagenomic investigations featuring high-throughput analysis of sequence data are increasingly common and important. We previously introduced MEGARes, a comprehensive AMR database with an acyclic hierarchical annotation structure that facilitates high-throughput computational analysis, as well as AMR++, a customized bioinformatic pipeline specifically designed to use MEGARes in high-throughput analysis for characterizing AMR genes (ARGs) in metagenomic sequence data. Here, we present MEGARes v3.0, a comprehensive database of published ARG sequences for antimicrobial drugs, biocides, and metals, and AMR++ v3.0, an update to our customized bioinformatic pipeline for high-throughput analysis of metagenomic data (available at MEGLab.org). Database annotations have been expanded to include information regarding specific genomic locations for single-nucleotide polymorphisms (SNPs) and insertions and/or deletions (indels) when required by specific ARGs for resistance expression, and the updated AMR++ pipeline uses this information to check for presence of resistance-conferring genetic variants in metagenomic sequenced reads. This new information encompasses 337 ARGs, whose resistance-conferring variants could not previously be confirmed in such a manner. In MEGARes 3.0, the nodes of the acyclic hierarchical ontology include 4 antimicrobial compound types, 59 resistance classes, 233 mechanisms and 1448 gene groups that classify the 8733 accessions.

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/.

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.

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.

Factors influencing usage of antimicrobial drugs among pastoralists in Kenya.

Agricultural use of antimicrobials in food animal production may contribute to the global emergence of antimicrobial resistance (AMR). However, considerable gaps exist in research on the use of antimicrobial drugs (AMDs) in food animals in small-scale production systems in low- and middle-income countries, despite the minimal regulation of antimicrobials in such regions. The aim of this study was to identify factors that may influence AMD use in livestock among pastoral communities in Kenya. We collected data related to household and herd demographics, herd health, and herd management from 55 households in the Maasai Mara ecosystem, Kenya, between 2018 and 2019. We used multi-model logistic regression inference (supervised machine learning) to ascertain trends in AMD use within these households. AMD use in cattle was significantly associated with AMD use in sheep and goats (p = 0.05), implying that decisions regarding AMD use in cattle or sheep and goats were interdependent. AMD use in sheep and goats was negatively associated with vaccination against the foot and mouth disease (FMD) virus in cattle (OR = 0.06, 95% CI 0.01-0.67, p = 0.02). Less AMD use was observed for vaccine-preventable diseases like contagious ecthyma when households had access to state veterinarians (OR = 0.06, p = 0.05, 95% CI 0.004-0.96). Overall, decisions to use AMDs were associated with vaccine usage, occurrence of respiratory diseases, and access to animal health advice. This hypothesis-generating study suggests that applying community-centric methods may be necessary to understand the use of AMDs in pastoral communities.

The impacts of viral infection and subsequent antimicrobials on the microbiome-resistome of growing pigs.

BACKGROUND: Antimicrobials are used in food-producing animals for purposes of preventing, controlling, and/or treating infections. In swine, a major driver of antimicrobial use is porcine reproductive and respiratory syndrome (PRRS), which is caused by a virus that predisposes infected animals to secondary bacterial infections. Numerous antimicrobial protocols are used to treat PRRS, but we have little insight into how these treatment schemes impact antimicrobial resistance (AMR) dynamics within the fecal microbiome of commercial swine. The aim of this study was to determine whether different PRRS-relevant antimicrobial treatment protocols were associated with differences in the fecal microbiome and resistome of growing pigs. To accomplish this, we used a metagenomics approach to characterize and compare the longitudinal wean-to-market resistome and microbiome of pigs challenged with PRRS virus and then exposed to different antimicrobial treatments, and a group of control pigs not challenged with PRRS virus and having minimal antimicrobial exposure. Genomic DNA was extracted from pen-level composite fecal samples from each treatment group and subjected to metagenomic sequencing and microbiome-resistome bioinformatic and statistical analysis. Microbiome-resistome profiles were compared over time and between treatment groups. RESULTS: Fecal microbiome and resistome compositions both changed significantly over time, with a dramatic and stereotypic shift between weaning and 9 days post-weaning (dpw). Antimicrobial resistance gene (ARG) richness and diversity were significantly higher at earlier time points, while microbiome richness and diversity were significantly lower. The post-weaning shift was characterized by transition from a Bacteroides-dominated enterotype to Lactobacillus- and Streptococcus-dominated enterotypes. Both the microbiome and resistome stabilized by 44 dpw, at which point the trajectory of microbiome-resistome maturation began to diverge slightly between the treatment groups, potentially due to physical clustering of the pigs. Challenge with PRRS virus seemed to correspond to the re-appearance of many very rare and low-abundance ARGs within the feces of challenged pigs. Despite very different antimicrobial exposures after challenge with PRRS virus, resistome composition remained largely similar between the treatment groups. Differences in ARG abundance between the groups were mostly driven by temporal changes in abundance that occurred prior to antimicrobial exposures, with the exception of ermG, which increased in the feces of treated pigs, and was significantly more abundant in the feces of these pigs compared to the pigs that did not receive post-PRRS antimicrobials. CONCLUSIONS: The fecal microbiome-resistome of growing pigs exhibited a stereotypic trajectory driven largely by weaning and physiologic aging of the pigs. Events such as viral illness, antimicrobial exposures, and physical grouping of the pigs exerted significant yet relatively minor influence over this trajectory. Therefore, the AMR profile of market-age pigs is the culmination of the life history of the individual pigs and the populations to which they belong. Disease status alone may be a significant driver of AMR in market-age pigs, and understanding the interaction between disease processes and antimicrobial exposures on the swine microbiome-resistome is crucial to developing effective, robust, and reproducible interventions to control AMR. Video Abstract.

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 (day 98 to 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 days 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 days 1 and 3 along with pen-mate-untreated controls on day 1. Fecal samples were subjected to 16S rRNA amplicon-based microbiome analysis and culture methods to quantify the abundance of total 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 vs. 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 the 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. Furthermore, diets and antimicrobial treatment had minimal impact on the overall counts of AMR coliforms and enterococci populations in feces.

Bacterial communities in the gut and the feces are strongly influenced by a number of factors, particularly the age of the animal and the diet. In addition, antibiotic administration routinely used to treat bacterial diseases can also affect the community composition. A study with finisher pigs was conducted to evaluate age-related changes, effects of types­distiller's dried grains with solubles (DGGS) or sugar beet pulp (SBP)­and levels of dietary fiber, and injectable antibiotics on the fecal bacterial composition and antibiotic resistance in fecal bacteria. Fecal samples were collected from five pigs in each of the three dietary treatment groups, control diet with no supplement or supplemented with DDGS or SBP, on days 98, 110, 144, and 177 of age and on days 1 and 3 after the first injection of antibiotics, ceftiofur or penicillin G. Samples were analyzed to identify the bacterial community composition and prevalence of antibiotic resistance in fecal bacteria. Data generated suggested that the overall bacterial composition changed with age and diet, and age appeared to have a greater impact than diet. Antibiotics had only a modest impact on the bacterial community and had minimum impact on antibiotic resistance of fecal bacteria.

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.

AMR-meta: a k-mer and metafeature approach to classify antimicrobial resistance from high-throughput short-read metagenomics data.

BACKGROUND: Antimicrobial resistance (AMR) is a global health concern. High-throughput metagenomic sequencing of microbial samples enables profiling of AMR genes through comparison with curated AMR databases. However, the performance of current methods is often hampered by database incompleteness and the presence of homology/homoplasy with other non-AMR genes in sequenced samples. RESULTS: We present AMR-meta, a database-free and alignment-free approach, based on k-mers, which combines algebraic matrix factorization into metafeatures with regularized regression. Metafeatures capture multi-level gene diversity across the main antibiotic classes. AMR-meta takes in reads from metagenomic shotgun sequencing and outputs predictions about whether those reads contribute to resistance against specific classes of antibiotics. In addition, AMR-meta uses an augmented training strategy that joins an AMR gene database with non-AMR genes (used as negative examples). We compare AMR-meta with AMRPlusPlus, DeepARG, and Meta-MARC, further testing their ensemble via a voting system. In cross-validation, AMR-meta has a median f-score of 0.7 (interquartile range, 0.2-0.9). On semi-synthetic metagenomic data-external test-on average AMR-meta yields a 1.3-fold hit rate increase over existing methods. In terms of run-time, AMR-meta is 3 times faster than DeepARG, 30 times faster than Meta-MARC, and as fast as AMRPlusPlus. Finally, we note that differences in AMR ontologies and observed variance of all tools in classification outputs call for further development on standardization of benchmarking data and protocols. CONCLUSIONS: AMR-meta is a fast, accurate classifier that exploits non-AMR negative sets to improve sensitivity and specificity. The differences in AMR ontologies and the high variance of all tools in classification outputs call for the deployment of standard benchmarking data and protocols, to fairly compare AMR prediction tools.