The pangenome of the wheat pathogen Pyrenophora tritici-repentis reveals novel transposons associated with necrotrophic effectors ToxA and ToxB.

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 of ToxA) 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. ToxB and its putative transposon were missing from the ToxB non-coding reference isolate, but the homolog toxb and 'Icarus' were both present in a different non-coding isolate. This suggests that ToxB may have been mobile at some point during the evolution of the Ptr genome which is contradictory to the current assumption of ToxB vertical 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.

Comparison of Strategies for Isolating Anaerobic Bacteria from the Porcine Intestine.

The isolation of bacteria that represent the diversity of autochthonous taxa in the gastrointestinal tract is necessary to fully ascertain their function, but the majority of bacterial species inhabiting the intestines of mammals are fastidious and thus challenging to isolate. The goal of the current study was to isolate a diverse assemblage of anaerobic bacteria from the intestine of pigs as a model animal and to comparatively examine various novel and traditional isolation strategies. Methods used included long-term enrichments, direct plating, a modified ichip method, as well as ethanol and tyndallization treatments of samples to select for endospore-forming taxa. A total of 234 taxa (91 previously uncultured) comprising 80 genera and 7 phyla were isolated from mucosal and luminal samples from the ileum, cecum, ascending colon, and spiral colon removed from animals under anesthesia. The diversity of bacteria isolated from the large intestine was less than that detected by next-generation sequence analysis. Long-term enrichments yielded the greatest diversity of recovered bacteria (Shannon's index [SI] = 4.7). Methods designed to isolate endospore-forming bacteria produced the lowest diversity (SI ≤ 2.7), with tyndallization yielding lower diversity than the ethanol method. However, the isolation frequency of previously uncultured bacteria was highest for ethanol-treated samples (41.9%) and the ichip method (32.5%). The goal of recovering a diverse collection of enteric bacteria was achieved. Importantly, the study findings demonstrate that it is necessary to use a combination of methods in concert to isolate bacteria that are representative of the diversity within the intestines of mammals.IMPORTANCE This work determined that using a combination of anaerobic isolation methods is necessary to increase the diversity of bacteria recovered from the intestines of monogastric mammals. Direct plating methods have traditionally been used to isolate enteric bacteria, and recent methods (e.g., diffusion methods [i.e., ichip] or differential isolation of endospore-forming bacteria) have been suggested to be superior at increasing diversity, including the recovery of previously uncultured taxa. We showed that long-term enrichment of samples using a variety of media isolated the most diverse and novel bacteria. Application of the ichip method delivered a diversity of bacteria similar to those of enrichment and direct plating methods. Methods that selected for endospore-forming bacteria generated collections that differed in composition from those of other methods with reduced diversity. However, the ethanol treatment frequently isolated novel bacteria. By using a combination of methods in concert, a diverse collection of enteric bacteria was generated for ancillary experimentation.

Salmonella enterica Serovar Typhimurium Temporally Modulates the Enteric Microbiota and Host Responses To Overcome Colonization Resistance in Swine.

Salmonella enterica serovar Typhimurium is a prevalent incitant of enteritis in human beings and nonhuman animals. It has been proposed that host defense responses incited by Salmonella allow the bacterium to overcome colonization resistance. Piglets (n = 24) were orally inoculated with S. enterica serovar Typhimurium DT104 or buffer alone, and the host and microbial responses were temporally examined at the acute (2 days postinoculation [dpi]), subacute (6 dpi), and recovery (10 dpi) stages of salmonellosis. At the acute stage of disease, body temperatures were elevated, and feed consumption and weight gain were reduced. The densities of Salmonella associated with the gut mucosa decreased over time, with higher densities of the bacterium in the ileum and the large intestine. Moreover, substantive histopathological changes were observed as a function of time, with prominent epithelial injury and neutrophil infiltration observed at 2 dpi. Correspondingly, a variety of host metrics were temporally affected in piglets with salmonellosis (e.g., TNFα, IFNγ, PR39, ßD2, iNOS, IL8, REGIIIγ). The enteric microbiota was characterized using culture-independent and -dependent methods in concert, and taxon- and location-specific changes to the microbiota were observed in infected piglets. Bacteroides spp. (e.g., Bacteroides uniformis, Bacteroides fragilis), Streptococcus spp. (e.g., Streptococcus gallolyticus), and various Gammaproteobacteria were highly associated with inflamed tissues, while bacteria within the Ruminococcaceae and Veillonellaceae families were mainly associated with healthy mucosae. In conclusion, the study findings showed that S Typhimurium incited temporal and spatial modifications to the swine autochthonous microbiota, and to host defense responses, that were consistent with overcoming colonization resistance to incite salmonellosis in swine.IMPORTANCE Limited information is available on host and enteric microbiota responses incited by Salmonella enterica serovar Typhimurium in swine and on possible mechanisms by which the bacterium overcomes colonization resistance to incite salmonellosis. Temporal characterization of a variety of host metrics in piglets (e.g., physiological, histopathological, and immunological) showed the importance of studying the progression of salmonellosis. A number of host responses integrally associated with disease development were identified. Utilization of next-generation sequence analysis to characterize the enteric microbiota was found to lack sufficient resolution; however, culture-dependent and -independent methods in combination identified taxon- and location-specific changes to bacterial communities in infected piglets. The study identified bacterial and host responses associated with salmonellosis, which will be beneficial in understanding colonization resistance and in the development of effective alternatives to antibiotics to mitigate salmonellosis.

Comparative genomics of multidrug-resistant Enterococcus spp. isolated from wastewater treatment plants.

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.

Application of culturomics to characterize diverse anaerobic bacteria from the gastrointestinal tract of broiler chickens in relation to environmental reservoirs.

Characterization of the microbiota of chickens is of current interest. The goals of the current study were to apply anaerobic isolation methods to comprehensively isolate and identify bacteria from the gastrointestinal tract of chickens and their environment. Bacterial communities within the drinking water were dominated by Escherichia, whereas communities in litter were more representative of the cecum. The crop and small intestine (jejunum and ileum) were dominated by Lactobacillus and Enterococcus spp., and the cecum was dominated by Proteus spp. The collection of bacteria isolated was dominated by Enterococcus spp., Escherichia/Shigella spp., Lactobacillus spp., and Proteus spp.; however, many rare taxa were observed. These included members of the Clostridiales and Clostridium spp., which were commonly isolated from the ileum and cecum. Bacteria isolated by enrichment and direct plating differed. The selective de Man-Rogosa-Sharpe agar was commonly associated with the isolation of Lactobacillus spp. and yielded the lowest diversity of all methods utilized. Increased diversity and frequency of Clostridium spp. was observed in enrichments of blood and mucus or by plating on Columbia agar supplemented with 10% blood and gentamicin. The bacteria isolated from this study provide source material for genomic and functional studies in chicken hosts.

Biogeography and genotypic diversity of Metarhizium brunneum and Metarhizium robertsii in northwestern North American soils.

The biogeography and genotype diversity of Metarhizium species in northwestern North American soils was examined; 20 ecoregions were sampled, including 58 agricultural and 80 natural habitat subsites, and areas that were glaciated during the Pleistocene epoch. One hundred and twenty-nine isolates of M. brunneum, 26 isolates of M. robertsii, four isolates of M. guizhouense, one isolate of M. flavoviride, and 55 isolates of Beauveria were recovered. Metarhizium and Beauveria species were isolated in diverse ecoregions within the study area, but a trend for increased isolation of Metarhizium species in western regions of the study area was observed. Consistent with this observation, the prevalence of M. brunneum and M. robertsii decreased at higher elevations, and the opposite was true for Beauveria. Both M. brunneum and M. robertsii were more commonly isolated from agricultural and natural habitat subsites, and considerable genotypic diversity was observed in both habitats and within the same subsite. Metarhizium robertsii, but not M. brunneum, was more commonly isolated from nonglaciated locations; however, less diversity and richness was observed for M. brunneum recovered from glaciated versus nonglaciated locations consistent with insular biogeography. The study has implications for microbial control strategies in the region.

Characterization of the Hoof Bacterial Communities of Active Digital Dermatitis Lesions in Feedlot Cattle.

Digital dermatitis (DD) is a costly hoof infection, causing lameness and pain in feedlot cattle. DD lesions can develop nonlinearly through a series of clinical stages, which can be classified by Dopfer's M-stage scoring system. This widely adopted lesion scoring system recognizes five DD stages, where M1 (early lesion), M2 (acute ulcerative lesion), and M4.1 (chronic proliferative lesion with new developing lesion) are considered active but separate stages of the disease. This study assessed the skin surface microbiota of the active DD lesions of feedlot cattle. The DD lesions from three commercial feedlots were swabbed and then scored according to Dopfer's M-stage scoring system. Swab samples were collected from 12 M2- and 15 M4.1-stage lesions. A total of 21 control swab samples from healthy contralateral feet (DD control) were classified as stage M0. An additional six skin swabs (M0) were collected from completely healthy (CH control) cattle with no lesions. The bacterial communities of active DD lesions (M2 and M4.1) and healthy skin (M0) were profiled using 16S amplicon sequencing. Diversity analyses showed that the hoof bacterial communities of M2 and M4.1 lesions were each distinct from those of M0 skin. However, the bacterial communities between the two active lesion stages were not different from each other. A significant increase in the relative abundance of Spirochaetota and Fusobacteriota and an overall decrease in bacterial diversity contributed to the altered bacterial communities in M2 and M4.1 lesions compared to those of healthy skin (M0). Although stages M2 and M4.1 are considered clinically different stages, the lesion-associated bacterial community is similar between the two active stages.

Characterization of the hoof bacterial communities in feedlot cattle affected with digital dermatitis, foot rot or both using a surface swab technique.

BACKGROUND: Lameness is defined as altered or abnormal gait due to dysfunction of the locomotor system, and is a health issue of feedlot cattle, having major economic, labour, and welfare implications. Digital dermatitis (DD-a lesion of the plantar surface of the foot) and foot rot (FR-affects the interdigital cleft) are common infectious causes of lameness in feedlots. These hoof lesions can occur alone or in combination (DD + FR) in the same hoof. A total of 208 hoof swabs were collected from three commercial feedlots located in southern Alberta. Every lesion sample was matched with a corresponding control skin sample taken from a healthy contralateral foot. Control skin samples were also collected from cattle with no lesion on any feet. Bacterial communities of three types of hoof lesions (DD, DD + FR, FR) and healthy skin were profiled using 16S amplicon sequencing. RESULTS: Alpha diversity analysis revealed a lower bacterial diversity on DD and FR lesions compared to control skin. Beta diversity analysis showed that bacterial communities of DD, FR, and DD + FR lesions were distinct from those of the control skin. While the impact of feedlot was minimal, lesion type contributed to 22% of the variation observed among bacterial communities (PERMANOVA-R = 0.22, P < 0.01). Compared to the corresponding control skin, there were 11, 12, and 3 differentially abundant (DA) bacterial genera in DD, DD + FR, and FR lesions, respectively. CONCLUSIONS: The bacterial community description of a DD + FR lesion is a novel finding. Not only did lesions lead to altered bacterial communities when compared to healthy skin, but the composition of those communities also differed depending on the hoof lesion. The 16S amplicon sequencing of surface swabs has significant value as a research tool in separating different hoof lesions and can provide additional insights to the polybacterial etiology of DD and FR in feedlot cattle.

Elucidation of Physiological, Transcriptomic and Metabolomic Salinity Response Mechanisms in Medicago sativa.

Alfalfa (Medicago sativa L.) is a widely grown perennial leguminous forage crop with a number of positive attributes. However, despite its moderate ability to tolerate saline soils, which are increasing in prevalence worldwide, it suffers considerable yield declines under these growth conditions. While a general framework of the cascade of events involved in plant salinity response has been unraveled in recent years, many gaps remain in our understanding of the precise molecular mechanisms involved in this process, particularly in non-model yet economically important species such as alfalfa. Therefore, as a means of further elucidating salinity response mechanisms in this species, we carried out in-depth physiological assessments of M. sativa cv. Beaver, as well as transcriptomic and untargeted metabolomic evaluations of leaf tissues, following extended exposure to salinity (grown for 3-4 weeks under saline treatment) and control conditions. In addition to the substantial growth and photosynthetic reductions observed under salinity treatment, we identified 1233 significant differentially expressed genes between growth conditions, as well as 60 annotated differentially accumulated metabolites. Taken together, our results suggest that changes to cell membranes and walls, cuticular and/or epicuticular waxes, osmoprotectant levels, antioxidant-related metabolic pathways, and the expression of genes encoding ion transporters, protective proteins, and transcription factors are likely involved in alfalfa's salinity response process. Although some of these alterations may contribute to alfalfa's modest salinity resilience, it is feasible that several may be disadvantageous in this context and could therefore provide valuable targets for the further improvement of tolerance to this stress in the future.

Identification of Differential Drought Response Mechanisms in Medicago sativa subsp. sativa and falcata through Comparative Assessments at the Physiological, Biochemical, and Transcriptional Levels.

Alfalfa (Medicago sativa L.) is an extensively grown perennial forage legume, and although it is relatively drought tolerant, it consumes high amounts of water and depends upon irrigation in many regions. Given the progressive decline in water available for irrigation, as well as an escalation in climate change-related droughts, there is a critical need to develop alfalfa cultivars with improved drought resilience. M. sativa subsp. falcata is a close relative of the predominantly cultivated M. sativa subsp. sativa, and certain accessions have been demonstrated to exhibit superior performance under drought. As such, we endeavoured to carry out comparative physiological, biochemical, and transcriptomic evaluations of an as of yet unstudied drought-tolerant M. sativa subsp. falcata accession (PI 641381) and a relatively drought-susceptible M. sativa subsp. sativa cultivar (Beaver) to increase our understanding of the molecular mechanisms behind the enhanced ability of falcata to withstand water deficiency. Our findings indicate that unlike the small number of falcata genotypes assessed previously, falcata PI 641381 may exploit smaller, thicker leaves, as well as an increase in the baseline transcriptional levels of genes encoding particular transcription factors, protective proteins, and enzymes involved in the biosynthesis of stress-related compounds. These findings imply that different falcata accessions/genotypes may employ distinct drought response mechanisms, and the study provides a suite of candidate genes to facilitate the breeding of alfalfa with enhanced drought resilience in the future.

Inconsistent PCR detection of Shiga toxin-producing Escherichia coli: Insights from whole genome sequence analyses.

Shiga toxin-producing Escherichia coli (STEC) have been linked to food-borne disease outbreaks. As PCR is routinely used to screen foods for STEC, it is important that factors leading to inconsistent detection of STEC by PCR are understood. This study used whole genome sequencing (WGS) to investigate causes of inconsistent PCR detection of stx1, stx2, and serogroup-specific genes. Fifty strains isolated from Alberta feedlot cattle from three different studies were selected with inconsistent or consistent detection of stx and serogroup by PCR. All isolates were initially classified as STEC by PCR. Sequencing was performed using Illumina MiSeq® with sample library by Nextera XT. Virtual PCRs were performed using Geneious and bacteriophage content was determined using PHASTER. Sequencing coverage ranged from 47 to 102x, averaging 74x, with sequences deposited in the NCBI database. Eleven strains were confirmed by WGS as STEC having complete stxA and stxB subunits. However, truncated stx fragments occurred in twenty-two other isolates, some having multiple stx fragments in the genome. Isolates with complete stx by WGS had consistent stx1 and stx2 detection by PCR, although one also having a stx2 fragment had inconsistent stx2 PCR. For all STEC and 18/39 non-STEC, serogroups determined by PCR agreed with those determined by WGS. An additional three WGS serotypes were inconclusive and two isolates were Citrobacter spp. Results demonstrate that stx fragments associated with stx-carrying bacteriophages in the E. coli genome may contribute to inconsistent detection of stx1 and stx2 by PCR. Fourteen isolates had integrated stx bacteriophage but lacked complete or fragmentary stx possibly due to partial bacteriophage excision after sub-cultivation or other unclear mechanisms. The majority of STEC isolates (7/11) did not have identifiable bacteriophage DNA in the contig(s) where stx was located, likely increasing the stability of stx in the bacterial genome and its detection by PCR.

Whole-Genome Draft Assemblies of Difficult-to-Classify Escherichia coli O157 and Non-O157 Isolates from Feces of Canadian Feedlot Cattle.

Forty-eight Escherichia coli strains were chosen due to variable detection of stx or serogroup by PCR. Although all strains were initially determined to be Shiga toxin-producing Escherichia coli (STEC), their genomes revealed 11 isolates carrying stx 1a, stx 1b, stx 2a, and/or stx 2b Assembled genome sizes varied between 4,667,418 and 5,556,121 bp, with N 50 values between 79,648 and 294,166 bp and G+C contents between 50.3% and 51.4%.

Comparison of biochemical and genotypic speciation methods for vancomycin-resistant enterococci isolated from urban wastewater treatment plants.

Enterococci species in wastewater including Enterococcus faecalis, Enterococcus faecium, Enterococcus casseliflavus and Enterococcus gallinarum isolates (n = 308) with low or high level vancomycin resistance were determined and compared using a phenotypic method (RapID™ STR system), 16S rRNA sequencing, and multi-locus (atpA, groESL, and pheS) sequence analysis (MLSA). Error rates for the RapID™ STR system were E. faecalis (15.9%), E. faecium (21.5%), and E. casseliflavus/E. gallinarum (56.9%) when referenced to the consensus of all methods tested. Comparison of single nucleotide polymorphism (SNP) distances and phylogenetic trees suggested that the groESL locus delineated species more effectively than other loci. The groESL locus was the most reliable loci for the correct identification of Enterococcus spp., including E. casseliflavus and E. gallinarum, with high congruence compared to the consensus (Adjusted Rand Index = 0.954; Adjusted Wallace Co-efficient = 0.941). All of the methods were compared to whole genome sequencing, which acted as a gold standard, for the isolates from this study and those downloaded from NCBI.

Serotyping and antimicrobial resistance of Mannheimia haemolytica strains from European cattle with bovine respiratory disease.

The objective of this study was to assess the prevalence of three serotypes, A1, A2, and A6 in 98 M. haemolytica isolates collected from clinical BRD cases in European cattle and assess their antimicrobial resistance profiles. Isolates were characterized by serotyping (plate agglutination and serotype specific PCR) and antimicrobial susceptibility testing. The study identified a predominance of serotypes A1 (59%) and A6 (22%) in European M. haemolytica isolates exhibiting a relatively low level of antimicrobial resistance. A comprehensive understanding of the relative prevalence of different M. haemolytica serotypes in Europe informs a targeted approach for vaccine design against BRD.

Characterization of Non-O157 Escherichia coli from Cattle Faecal Samples in the North-West Province of South Africa.

Escherichia coli are commensal bacteria in the gastrointestinal tract of mammals, but some strains have acquired Shiga-toxins and can cause enterohemorrhagic diarrhoea and kidney failure in humans. Shiga-toxigenic E. coli (STEC) strains such as E. coli O157:H7 and some non-O157 strains also contain other virulence traits, some of which contribute to their ability to form biofilms. This study characterized non-O157 E. coli from South African cattle faecal samples for their virulence potential, antimicrobial resistance (AMR), biofilm-forming ability, and genetic relatedness using culture-based methods, pulsed-field gel electrophoresis (PFGE), and whole genome sequencing (WGS). Of 80 isolates screened, 77.5% (62/80) possessed Shiga-toxins genes. Of 18 antimicrobials tested, phenotypic resistance was detected against seven antimicrobials. Resistance ranged from 1.3% (1/80) for ampicillin-sulbactam to 20% (16/80) for tetracycline. Antimicrobial resistance genes were infrequently detected except for tetA, which was found in 31.3% (25/80) and tetB detected in 11.3% (9/80) of isolates. Eight biofilm-forming associated genes were detected in STEC isolates (n = 62) and two non-STEC strains. Prevalence of biofilm genes ranged from 31.3% (20/64) for ehaAß passenger to 100% for curli structural subunit (csgA) and curli regulators (csgA and crl). Of the 64 STEC and multi-drug resistant isolates, 70.3% (45/64) and 37.5% (24/64) formed strong biofilms on polystyrene at 22 and 37 °C, respectively. Of 59 isolates screened by PFGE, 37 showed unique patterns and the remaining isolates were grouped into five clusters with a ≥90% relatedness. In silico serotyping following WGS on a subset of 24 non-O157 STEC isolates predicted 20 serotypes comprising three novel serotypes, indicating their diversity as potential pathogens. These findings show that North West South African cattle harbour genetically diverse, virulent, antimicrobial-resistant and biofilm-forming non-O157 E. coli. Biofilm-forming ability may increase the likelihood of persistence of these pathogens in the environment and facilitate their dissemination, increasing the risk of cross contamination or establishment of infections in hosts.

Quantification and Multidrug Resistance Profiles of Vancomycin-Resistant Enterococci Isolated from Two Wastewater Treatment Plants in the Same Municipality.

Wastewater treatment plants (WWTPs) are points of control for the environmental dissemination of antimicrobial resistant bacteria. Vancomycin-resistant enterococci (VRE) were used as indicators of antimicrobial resistance (AMR) in two WWTPs (biologically aerated filter (BAF) and conventional activated sludge (CAS)) in the same municipality. The removal and abundance of enterococci and VRE as well as the species and antimicrobial resistance profiles of VRE were assessed. Enterococci and VRE from the primary and final effluents were enumerated. Results were assessed from an ecological context. VRE was not selected for by either WWTP but the BAF system outperformed the CAS system for the removal of enterococci/VRE. Enterococcus faecalis (n = 151), E. faecium (n = 94) and E. casseliflavus/E. gallinarum (n = 59) were the dominant VRE species isolated. A decrease in levofloxacin resistance in enterococci was observed in the BAF WWTP. An increase in nitrofurantoin resistant (p < 0.001) and a decrease in quinupristin/dalfopristin (p = 0.003) and streptomycin (p = 0.022) resistant enterococci were observed in the CAS WWTP, corresponding to a shift of VRE from E. faecalis to E. faecium. Wastewater treatment processes can be managed to limit the dissemination of antimicrobial resistance determinants into the surrounding environment.

Impact of sequencing depth on the characterization of the microbiome and resistome.

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.