BackgroundShotgun 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.ResultsThe 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.ConclusionsThese 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.Electronic supplementary materialThe online version of this article (10.1186/s40168-017-0361-8) contains supplementary material, which is available to authorized users.

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.

Microbiome studies in animal science using 16S rRNA gene sequencing have become increasingly common in recent years as sequencing costs continue to fall and bioinformatic tools become more powerful and user-friendly. The combination of molecular biology, microbiology, microbial ecology, computer science, and bioinformatics—in addition to the traditional considerations when conducting an animal science study—makes microbiome studies sometimes intimidating due to the intersection of different fields. The objective of this review is to serve as a jumping-off point for those animal scientists less familiar with 16S rRNA gene sequencing and analyses and to bring up common issues and concerns that arise when planning an animal microbiome study from design through analysis. This review includes an overview of 16S rRNA gene sequencing, its advantages, and its limitations; experimental design considerations such as study design, sample size, sample pooling, and sample locations; wet lab considerations such as field handing, microbial cell lysis, low biomass samples, library preparation, and sequencing controls; and computational considerations such as identification of contamination, accounting for uneven sequencing depth, constructing diversity metrics, assigning taxonomy, differential abundance testing, and, finally, data availability. In addition to general considerations, we highlight some special considerations by species and sample type.

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.

No abstract