Developments in high-throughput next generation sequencing (NGS) technology have rapidly advanced the understanding of overall microbial ecology as well as occurrence and diversity of specific genes within diverse environments. In the present study, we compared the ability of varying sequencing depths to generate meaningful information about the taxonomic structure and prevalence of antimicrobial resistance genes (ARGs) in the bovine fecal microbial community. Metagenomic sequencing was conducted on eight composite fecal samples originating from four beef cattle feedlots. Metagenomic DNA was sequenced to various depths, D1, D0.5 and D0.25, with average sample read counts of 117, 59 and 26 million, respectively. A comparative analysis of the relative abundance of reads aligning to different phyla and antimicrobial classes indicated that the relative proportions of read assignments remained fairly constant regardless of depth. However, the number of reads being assigned to ARGs as well as to microbial taxa increased significantly with increasing depth. We found a depth of D0.5 was suitable to describe the microbiome and resistome of cattle fecal samples. This study helps define a balance between cost and required sequencing depth to acquire meaningful results.

For a One-Health investigation of antimicrobial resistance (AMR) inEnterococcusspp., isolates from humans and beef cattle along with abattoirs, manured fields, natural streams, and wastewater from both urban and cattle feedlot sources were collected over two years. Species identification ofEnterococcusrevealed distinct associations across the continuum. Of the 8430 isolates collected,Enterococcus faeciumandEnterococcus faecaliswere the main species in urban wastewater (90%) and clinical human isolates (99%);Enterococcus hiraepredominated in cattle (92%) and feedlot catch-basins (60%), whereas natural streams harbored environmentalEnterococcusspp. Whole-genome sequencing ofE. faecalis(n = 366 isolates) andE. faecium(n = 342 isolates), revealed source clustering of isolates, indicative of distinct adaptation to their respective environments. Phenotypic resistance to tetracyclines and macrolides encoded bytet(M)anderm(B)respectively, was prevalent amongEnterococcusspp. regardless of source. ForE. faeciumfrom cattle, resistance to β-lactams and quinolones was observed among 3% and 8% of isolates respectively, compared to 76% and 70% of human clinical isolates. Clinical vancomycin-resistantE. faeciumexhibited high rates of multi-drug resistance, with resistance to all β-lactam, macrolides, and quinolones tested. Differences in the AMR profiles among isolates reflected antimicrobial use practices in each sector of the One-Health continuum.

Background: In fungal plant pathogens, genome rearrangements followed by selection pressure for adaptive traits have facilitated the co-evolutionary arms race between hosts and their pathogens.Pyrenophora tritici-repentis(Ptr) has emerged recently as a foliar pathogen of wheat worldwide and its populations consist of isolates that vary in their ability to produce combinations of different necrotrophic effectors. These effectors play vital roles in disease development. Here, we sequenced the genomes of a global collection (40 isolates) of Ptr to gain insights into its gene content and genome rearrangements.Results: A comparative genome analysis revealed an open pangenome, with an abundance of accessory genes (~ 57%) reflecting Ptr’s adaptability. A clear distinction between pathogenic and non-pathogenic genomes was observed in size, gene content, and phylogenetic relatedness. Chromosomal rearrangements and structural organization, specifically around effector coding genes, were detailed using long-read assemblies (PacBio RS II) generated in this work in addition to previously assembled genomes. We also discovered the involvement of large mobile elements associated with Ptr’s effectors:ToxA, the gene encoding for the necrosis effector, was found as a single copy within a 143-kb ‘Starship’ transposon (dubbed ‘Horizon’) with a clearly defined target site and target site duplications. ‘Horizon’ was located on different chromosomes in different isolates, indicating mobility, and the previously described ToxhAT transposon (responsible for horizontal transfer ofToxA) was nested within this newly identified Starship. Additionally,ToxB, the gene encoding the chlorosis effector, was clustered as three copies on a 294-kb element, which is likely a different putative ‘Starship’ (dubbed ‘Icarus’) in a ToxB-producing isolate.ToxBand its putative transposon were missing from theToxBnon-coding reference isolate, but the homologtoxband ‘Icarus’ were both present in a different non-coding isolate. This suggests thatToxBmay have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption ofToxBvertical inheritance. Finally, the genome architecture of Ptr was defined as ‘one-compartment’ based on calculated gene distances and evolutionary rates.Conclusions: These findings together reflect on the highly plastic nature of the Ptr genome which has likely helped to drive its worldwide adaptation and has illuminated the involvement of giant transposons in facilitating the evolution of virulence in Ptr.

Background: Wastewater treatment plants (WWTPs) are considered hotspots for the environmental dissemination of antimicrobial resistance (AMR) determinants. Vancomycin-Resistant Enterococcus (VRE) are candidates for gauging the degree of AMR bacteria in wastewater. Enterococcus faecalis and Enterococcus faecium are recognized indicators of fecal contamination in water. Comparative genomics of enterococci isolated from conventional activated sludge (CAS) and biological aerated filter (BAF) WWTPs was conducted. Results: VRE isolates, including E. faecalis (n = 24), E. faecium (n = 11), E. casseliflavus (n = 2) and E. gallinarum (n = 2) were selected for sequencing based on WWTP source, species and AMR phenotype. The pangenomes of E. faecium and E. faecalis were both open. The genomic fraction related to the mobilome was positively correlated with genome size in E. faecium (p < 0.001) and E. faecalis (p < 0.001) and with the number of AMR genes in E. faecium (p = 0.005). Genes conferring vancomycin resistance, including vanA and vanM (E. faecium), vanG (E. faecalis), and vanC (E. casseliflavus/E. gallinarum), were detected in 20 genomes. The most prominent functional AMR genes were efflux pumps and transporters. A minimum of 16, 6, 5 and 3 virulence genes were detected in E. faecium, E. faecalis, E. casseliflavus and E. gallinarum, respectively. Virulence genes were more common in E. faecalis and E. faecium, than E. casseliflavus and E. gallinarum. A number of mobile genetic elements were shared among species. Functional CRISPR/Cas arrays were detected in 13 E. faecalis genomes, with all but one also containing a prophage. The lack of a functional CRISPR/Cas arrays was associated with multi-drug resistance in E. faecium. Phylogenetic analysis demonstrated differential clustering of isolates based on original source but not WWTP. Genes related to phage and CRISPR/Cas arrays could potentially serve as environmental biomarkers. Conclusions: There was no discernible difference between enterococcal genomes from the CAS and BAF WWTPs. E. faecalis and E. faecium have smaller genomes and harbor more virulence, AMR, and mobile genetic elements than other Enterococcus spp.

Comparative knowledge of microbiomes and resistomes across environmental interfaces between animal production systems and urban settings is lacking. In this study, we executed a comparative analysis of the microbiota and resistomes of metagenomes from cattle feces, catch basin water, manured agricultural soil and urban sewage. Metagenomic DNA from composite fecal samples (FC; n = 12) collected from penned cattle at four feedlots in Alberta, Canada, along with water from adjacent catchment basins (CB; n = 13), soil (n = 4) from fields in the vicinity of one of the feedlots and urban sewage influent (SI; n = 6) from two municipalities were subjected to Illumina HiSeq2000 sequencing. Firmicutes exhibited the highest prevalence (40%) in FC, whereas Proteobacteria were most abundant in CB (64%), soil (60%) and SI (83%). Among sample types, SI had the highest diversity of antimicrobial resistance (AMR), and metal and biocide resistance (MBR) classes (13 & 15) followed by FC (10 & 8), CB (8 & 4), and soil (6 & 1). The highest antimicrobial resistant (AMR) gene (ARG) abundance was harboured by FC, whereas soil samples had a very small, but unique resistome which did not overlap with FC & CB resistomes. In the beef production system, tetracycline resistance predominated followed by macrolide resistance. The SI resistome harboured β-lactam, macrolide, tetracycline, aminoglycoside, fluoroquinolone and fosfomycin resistance determinants. Metal and biocide resistance accounted for 26% of the SI resistome with a predominance of mercury resistance. This study demonstrates an increasing divergence in the nature of the microbiome and resistome as the distance from the feedlot increases. Consistent with antimicrobial use, tetracycline and macrolide resistance genes were predominant in the beef production system. One of the feedlots contributed both conventional (raised with antibiotics) and natural (raised without antibiotics) pens samples. Although natural pen samples exhibited a microbiota composition that was similar to samples from conventional pens, their resistome was less complex. Similarly, the SI resistome was indicative of drug classes used in humans and the greater abundance of mercury resistance may be associated with contamination of municipal water with household and industrial products.