Comparative genomics of Salmonella enterica serovar Montevideo reveals lineage-specific gene differences that may influence ecological niche association

Salmonella enterica serovar Montevideo has been linked to recent foodborne illness outbreaks resulting from contamination of products such as fruits, vegetables, seeds and spices. Studies have shown that Montevideo also is frequently associated with healthy cattle and can be isolated from ground beef, yet human salmonellosis outbreaks of Montevideo associated with ground beef contamination are rare. This disparity fuelled our interest in characterizing the genomic differences between Montevideo strains isolated from healthy cattle and beef products, and those isolated from human patients and outbreak sources. To that end, we sequenced 13 Montevideo strains to completion, producing high-quality genome assemblies of isolates from human patients (n=8) or from healthy cattle at slaughter (n=5). Comparative analysis of sequence data from this study and publicly available sequences (n=72) shows that Montevideo falls into four previously established clades, differentially occupied by cattle and human strains. The results of these analyses reveal differences in metabolic islands, environmental adhesion determinants and virulence factors within each clade, and suggest explanations for the infrequent association between bovine isolates and human illnesses.

Comparative Genomics of All Three Campylobacter sputorum Biovars and a Novel Cattle-Associated C. sputorum Clade

Campylobacter sputorum is a nonthermotolerant campylobacter that is primarily isolated from food animals such as cattle and sheep. C. sputorum is also infrequently associated with human illness. Based on catalase and urease activity, three biovars are currently recognized within C. sputorum: bv. sputorum (catalase negative, urease negative), bv. fecalis (catalase positive, urease negative), and bv. paraureolyticus (catalase negative, urease positive). A multi-locus sequence typing (MLST) method was recently constructed for C. sputorum. MLST typing of several cattle-associated C. sputorum isolates suggested that they are members of a divergent C. sputorum clade. Although catalase positive, and thus technically bv. fecalis, the taxonomic position of these strains could not be determined solely by MLST. To further characterize C. sputorum, the genomes of four strains, representing all three biovars and the divergent clade, were sequenced to completion. Here we present a comparative genomic analysis of the four C. sputorum genomes. This analysis indicates that the three biovars and the cattle-associated strains are highly related at the genome level with similarities in gene content. Furthermore, the four genomes are strongly syntenic with one or two minor inversions. However, substantial differences in gene content were observed among the three biovars. Finally, although the strain representing the cattle-associated isolates was shown to be C. sputorum, it is possible that this strain is a member of a novel C. sputorum subspecies; thus, these cattle-associated strains may form a second taxon within C. sputorum.

Comparative Genomic Analysis Identifies a Campylobacter Clade Deficient in Selenium Metabolism

The nonthermotolerant Campylobacter species C. fetus, C. hyointestinalis, C. iguaniorum, and C. lanienae form a distinct phylogenetic cluster within the genus. These species are primarily isolated from foraging (swine) or grazing (e.g., cattle, sheep) animals and cause sporadic and infrequent human illness. Previous typing studies identified three putative novel C. lanienae-related taxa, based on either MLST or atpA sequence data. To further characterize these putative novel taxa and the C. fetus group as a whole, 76 genomes were sequenced, either to completion or to draft level. These genomes represent 26 C. lanienae strains and 50 strains of the three novel taxa. C. fetus, C. hyointestinalis and C. iguaniorum genomes were previously sequenced to completion; therefore, a comparative genomic analysis across the entire C. fetus group was conducted (including average nucleotide identity analysis) that supports the initial identification of these three novel Campylobacter species. Furthermore, C. lanienae and the three putative novel species form a discrete clade within the C. fetus group, which we have termed the C. lanienae clade. This clade is distinguished from other members of the C. fetus group by a reduced genome size and distinct CRISPR/Cas systems. Moreover, there are two signature characteristics of the C. lanienae clade. C. lanienae clade genomes carry four to ten unlinked and similar, but nonidentical, flagellin genes. Additionally, all 76 C. lanienae clade genomes sequenced demonstrate a complete absence of genes related to selenium metabolism, including genes encoding the selenocysteine insertion machinery, selenoproteins, and the selenocysteinyl tRNA.

The impact of the bovine faecal microbiome on Escherichia coli O157:H7 prevalence and enumeration in naturally infected cattle

AIMS

The objective of this study was to determine if the faecal microbiome has an association with Escherichia coli O157:H7 prevalence and enumeration.

METHODS AND RESULTS

Pyrosequencing analysis of faecal microbiome was performed from feedlot cattle fed one of three diets: (i) 94 heifers fed low concentrate (LC) diet, (ii) 142 steers fed moderate concentrate (MC) diet, and (iii) 132 steers fed high concentrate (HC) diet. A total of 322 585 OTUs were calculated from 2,411,122 high-quality sequences obtained from 368 faecal samples. In the LC diet group, OTUs assigned to the orders Clostridiales and RF39 (placed within the class Mollicutes) were positively correlated with both E. coli O157:H7 prevalence and enumeration. In the MC diet group, OTUs assigned to Prevotella copri were positively correlated with both E. coli O157:H7 prevalence and enumeration, whereas OTUs assigned to Prevotella stercorea were negatively correlated with both E. coli O157:H7 prevalence and enumeration. In both the MC diet group and the HC diet group, OTUs assigned to taxa placed within Clostridiales were both positively and negatively correlated with both E. coli O157:H7 prevalence and enumeration. However, all correlations were weak. In both the MC diet group and the HC diet group, stepwise linear regression through backward elimination analyses indicated that these OTUs were significantly correlated (P < 0·001) with prevalence or enumeration, explaining as much as 50% of variability in E. coli O157:H7 prevalence or enumeration.

CONCLUSIONS

Individual colonic bacterial species have little impact on E. coli O157:H7 shedding but collectively groups of bacteria were strongly associated with pathogen shedding.

SIGNIFICANCE AND IMPACT OF THE STUDY

Bacterial groups in the bovine colon may impact faecal shedding of the zoonotic pathogen E. coli O157:H7, and manipulation of the intestinal microbiota to alter these bacteria may reduce shedding of this pathogen and foodborne illnesses.

The effects of monensin in diets fed to finishing beef steers and heifers on growth performance and fecal shedding of Escherichia coli O157:H7

Two experiments were conducted to study the effects of monensin dose on growth performance and O157:H7 shedding in finishing beef cattle. In Exp. 1, 198 heifers (298 ± 1.1 kg BW) were allocated to 1 of 2 treatments consisting of 1) 200 mg/heifer daily of monensin and 2) 400 mg/heifer daily of monensin and fed for 151 d. In Exp. 2, 199 steers (430 ± 1.9 kg BW) were stratified by BW and allocated to 1 of 2 treatments consisting of 1) 0 mg/steer daily of monensin and 2) 400 mg/steer daily of monensin and fed for 128 d. For both experiments, there were 4 pen replicates per treatment. For Exp. 1 and Exp. 2, the model included the fixed effect of treatment for growth performance measures and the fixed effects of treatment, time, and treatment × time interaction, respectively, for O157:H7 shedding. In Exp. 1, final BW was 1.9% greater for heifers fed 400 mg/d monensin than for heifers fed 200 mg/d monensin ( = 0.05). Furthermore, ADG was 4.9% greater ( = 0.05) and G:F was 3.1% greater ( = 0.04) when the heifers were fed 400 mg/d monensin vs. 200 mg/d monensin. Pen prevalence for O157:H7 ( = 0.96) and the percentage of animals in the pen shedding O157:H7 at enumerable levels ( = 0.82) did not differ between heifers fed 200 mg/d monensin and heifers fed 400 mg/d monensin over the 4 sampling periods. For Exp. 2, steers fed the supplement containing monensin had a 1.9% greater final BW ( = 0.04) and a 5.2% greater ADG ( = 0.02) than steers fed a control supplement without monensin. No differences in DMI or G:F were noted across the treatments ( ≥ 0.14). O157:H7 percentage of enumerable cattle within the pen was greater for the steers fed monensin than the control steers not fed monensin than the control steers not fed monensin ( = 0.02) over the 4 sampling periods. However, the percentage of animals in the pen shedding O157:H7 (prevalence positive) did not differ between treatments ( = 0.18), nor did the average fecal counts ( = 0.45). In conclusion, feeding a higher dose (400 mg/d) of monensin improved final BW and ADG compared with a low dose of monensin or a no-monensin control in steers and heifers across multiple years. The percentage of animals shedding O157:H7 at enumerable levels was greater for steers fed the monensin supplement than for steers fed the control supplement, yet the presence of monensin, irrespective of the dose, did not affect the percentage of animals in the pen shedding O157:H7.

Evolution of a zoonotic pathogen: investigating prophage diversity in enterohaemorrhagic Escherichia coli O157 by long-read sequencing

Enterohaemorrhagic Escherichia coli (EHEC) O157 is a zoonotic pathogen for which colonization of cattle and virulence in humans is associated with multiple horizontally acquired genes, the majority present in active or cryptic prophages. Our understanding of the evolution and phylogeny of EHEC O157 continues to develop primarily based on core genome analyses; however, such short-read sequences have limited value for the analysis of prophage content and its chromosomal location. In this study, we applied Single Molecule Real Time (SMRT) sequencing, using the Pacific Biosciences long-read sequencing platform, to isolates selected from the main sub-clusters of this clonal group. Prophage regions were extracted from these sequences and from published reference strains. Genome position and prophage diversity were analysed along with genetic content. Prophages could be assigned to clusters, with smaller prophages generally exhibiting less diversity and preferential loss of structural genes. Prophages encoding Shiga toxin (Stx) 2a and Stx1a were the most diverse, and more variable compared to prophages encoding Stx2c, further supporting the hypothesis that Stx2c-prophage integration was ancestral to acquisition of other Stx types. The concept that phage type (PT) 21/28 (Stx2a+, Stx2c+) strains evolved from PT32 (Stx2c+) was supported by analysis of strains with excised Stx-encoding prophages. Insertion sequence elements were over-represented in prophage sequences compared to the rest of the genome, showing integration in key genes such as stx and an excisionase, the latter potentially acting to capture the bacteriophage into the genome. Prophage profiling should allow more accurate prediction of the pathogenic potential of isolates.

Genomic signatures of Mannheimia haemolytica that associate with the lungs of cattle with respiratory disease, an integrative conjugative element, and antibiotic resistance genes

BACKGROUND

Mannheimia haemolytica typically resides in cattle as a commensal member of the upper respiratory tract microbiome. However, some strains can invade their lungs and cause respiratory disease and death, including those with multi-drug resistance. A nucleotide polymorphism typing system was developed for M. haemolytica from the genome sequences of 1133 North American isolates, and used to identify genetic differences between isolates from the lungs and upper respiratory tract of cattle with and without clinical signs of respiratory disease.

RESULTS

A total of 26,081 nucleotide polymorphisms were characterized after quality control filtering of 48,403 putative polymorphisms. Phylogenetic analyses of nucleotide polymorphism genotypes split M. haemolytica into two major genotypes (1 and 2) that each were further divided into multiple subtypes. Multiple polymorphisms were identified with alleles that tagged genotypes 1 or 2, and their respective subtypes. Only genotype 2 M. haemolytica associated with the lungs of diseased cattle and the sequence of a particular integrative and conjugative element (ICE). Additionally, isolates belonging to one subtype of genotype 2 (2b), had the majority of antibiotic resistance genes detected in this study, which were assorted into seven combinations that ranged from 1 to 12 resistance genes.

CONCLUSIONS

Typing of diverse M. haemolytica by nucleotide polymorphism genotypes successfully identified associations with diseased cattle lungs, ICE sequence, and antibiotic resistance genes. Management of cattle by their carriage of M. haemolytica could be an effective intervention strategy to reduce the prevalence of respiratory disease and supplemental needs for antibiotic treatments in North American herds.

Short-term evolution of Shiga toxin-producing Escherichia coli O157:H7 between two food-borne outbreaks

Shiga toxin-producing Escherichia coli (STEC) O157:H7 is a public health threat and outbreaks occur worldwide. Here, we investigate genomic differences between related STEC O157:H7 that caused two outbreaks, eight weeks apart, at the same restaurant. Short-read genome sequencing divided the outbreak strains into two sub-clusters separated by only three single-nucleotide polymorphisms in the core genome while traditional typing identified them as separate phage types, PT8 and PT54. Isolates did not cluster with local strains but with those associated with foreign travel to the Middle East/North Africa. Combined long-read sequencing approaches and optical mapping revealed that the two outbreak strains had undergone significant microevolution in the accessory genome with prophage gain, loss and recombination. In addition, the PT54 sub-type had acquired a 240 kbp multi-drug resistance (MDR) IncHI2 plasmid responsible for the phage type switch. A PT54 isolate had a general fitness advantage over a PT8 isolate in rich medium, including an increased capacity to use specific amino acids and dipeptides as a nitrogen source. The second outbreak was considerably larger and there were multiple secondary cases indicative of effective human-to-human transmission. We speculate that MDR plasmid acquisition and prophage changes have adapted the PT54 strain for human infection and transmission. Our study shows the added insights provided by combining whole-genome sequencing approaches for outbreak investigations.

Disinfectant and Antimicrobial Susceptibility Profiles of the Big Six Non-O157 Shiga Toxin-Producing Escherichia coli Strains from Food Animals and Humans

The disinfectant and antimicrobial susceptibility profiles of 138 non-O157 Shiga toxin-producing Escherichia coli strains (STECs) from food animals and humans were determined. Antimicrobial resistance (AMR) was moderate (39.1% of strains) in response to 15 antimicrobial agents. Animal strains had a lower AMR prevalence (35.6%) than did human strains (43.9%) but a higher prevalence of the resistance profile GEN-KAN-TET. A decreasing prevalence of AMR was found among animal strains from serogroups O45 > O145 > O121 > O111 > O26 > O103 and among human strains from serogroups O145 > O103 > O26 > O111 > O121 > O45. One animal strain from serogroups O121 and O145 and one human strain from serogroup O26 had extensive drug resistance. A high prevalence of AMR in animal O45 and O121 strains and no resistance or a low prevalence of resistance in human strains from these serogroups suggests a source other than food animals for human exposure to these strains. Among the 24 disinfectants evaluated, all strains were susceptible to triclosan. Animal strains had a higher prevalence of resistance to chlorhexidine than did human strains. Both animal and human strains had a similar low prevalence of low-level benzalkonium chloride resistance, and animal and human strains had similar susceptibility profiles for most other disinfectants. Benzyldimethylammonium chlorides and C10AC were the primary active components in disinfectants DC&R and P-128, respectively, against non-O157 STECs. A disinfectant FS512 MIC ≥ 8 μg/ml was more prevalent among animal O121 strains (61.5%) than among human O121 strains (25%), which may also suggest a source of human exposure to STEC O121 other than food animals. Bacterial inhibition was not dependent solely on pH but was correlated with the presence of dissociated organic acid species and some undissociated acids.

Large genomic differences between Moraxella bovoculi isolates acquired from the eyes of cattle with infectious bovine keratoconjunctivitis versus the deep nasopharynx of asymptomatic cattle

Moraxella bovoculi is a recently described bacterium that is associated with infectious bovine keratoconjunctivitis (IBK) or “pinkeye” in cattle. In this study, closed circularized genomes were generated for seven M. bovoculi isolates: three that originated from the eyes of clinical IBK bovine cases and four from the deep nasopharynx of asymptomatic cattle. Isolates that originated from the eyes of IBK cases profoundly differed from those that originated from the nasopharynx of asymptomatic cattle in genome structure, gene content and polymorphism diversity and consequently placed into two distinct phylogenetic groups. These results suggest that there are genetically distinct strains of M. bovoculi that may not associate with IBK.

Sequence of Colonization Determines the Composition of Mixed Biofilms by Escherichia coli O157:H7 and O111:H8 Strains

Bacterial biofilms are one of the potential sources of cross-contamination in food processing environments. Shiga toxin-producing Escherichia coli (STEC) O157:H7 and O111:H8 are important foodborne pathogens capable of forming biofilms, and the coexistence of these two STEC serotypes has been detected in various food samples and in multiple commercial meat plants throughout the United States. Here, we investigated how the coexistence of these two STEC serotypes and their sequence of colonization could affect bacterial growth competition and mixed biofilm development. Our data showed that E. coli O157:H7 strains were able to maintain a higher cell percentage in mixed biofilms with the co-inoculated O111:H8 companion strains, even though the results of planktonic growth competition were strain dependent. On solid surfaces with preexisting biofilms, the sequence of colonization played a critical role in determining the composition of the mixed biofilms because early stage precolonization significantly affected the competition results between the E. coli O157:H7 and O111:H8 strains. The precolonizer of either serotype was able to outgrow the other serotype in both planktonic and biofilm phases. The competitive interactions among the various STEC serotypes would determine the composition and structure of the mixed biofilms as well as their potential risks to food safety and public health, which is largely influenced by the dominant strains in the mixtures. Thus, the analysis of mixed biofilms under various conditions would be of importance to determine the nature of mixed biofilms composed of multiple microorganisms and to help implement the most effective disinfection operations accordingly.

Canonical Single Nucleotide Polymorphisms (SNPs) for High-Resolution Subtyping of Shiga-Toxin Producing Escherichia coli (STEC) O157:H7

The objective of this study was to develop a canonical, parsimoniously-informative SNP panel for subtyping Shiga-toxin producing Escherichia coli (STEC) O157:H7 that would be consistent with epidemiological, PFGE, and MLVA clustering of human specimens. Our group had previously identified 906 putative discriminatory SNPs, which were pared down to 391 SNPs based on their prevalence in a test set. The 391 SNPs were screened using a high-throughput form of TaqMan PCR against a set of clinical isolates that represent the most diverse collection of O157:H7 isolates from outbreaks and sporadic cases examined to date. Another 30 SNPs identified by others were also screened using the same method. Two additional targets were tested using standard TaqMan PCR endpoint analysis. These 423 SNPs were reduced to a 32 SNP panel with the almost the same discriminatory value. While the panel partitioned our diverse set of isolates in a manner that was consistent with epidemiological data and PFGE and MLVA phylogenies, it resulted in fewer subtypes than either existing method and insufficient epidemiological resolution in 10 of 47 clusters. Therefore, another round of SNP discovery was undertaken using comparative genomic resequencing of pooled DNA from the 10 clusters with insufficient resolution. This process identified 4,040 potential SNPs and suggested one of the ten clusters was incorrectly grouped. After its removal, there were 2,878 SNPs, of which only 63 were previously identified and 438 occurred across multiple clusters. Among highly clonal bacteria like STEC O157:H7, linkage disequilibrium greatly limits the number of parsimoniously informative SNPs. Therefore, it is perhaps unsurprising that our panel accounted for the potential discriminatory value of numerous other SNPs reported in the literature. We concluded published O157:H7 SNPs are insufficient for effective epidemiological subtyping. However, the 438 multi-cluster SNPs we identified may provide the additional information required.

Genome sequencing and comparative genomics provides insights on the evolutionary dynamics and pathogenic potential of different H-serotypes of Shiga toxin-producing Escherichia coli O104

BACKGROUND

Various H-serotypes of the Shiga toxin-producing Escherichia coli (STEC) O104, including H4, H7, H21, and H¯, have been associated with sporadic cases of illness and have caused food-borne outbreaks globally. In the U.S., STEC O104:H21 caused an outbreak associated with milk in 1994. However, there is little known on the evolutionary origins of STEC O104 strains, and how genotypic diversity contributes to pathogenic potential of various O104 H-antigen serotypes isolated from different ecological niches and/or geographical regions.

RESULTS

Two STEC O104:H21 (milk outbreak strain) and O104:H7 (cattle isolate) strains were shot-gun sequenced, and the genomes were closed. The intimin (eae) gene, involved in the attaching-effacing phenotype of diarrheagenic E. coli, was not found in either strain. Examining various O104 genome sequences, we found that two "complete" left and right end portions of the locus of enterocyte effacement (LEE) pathogenicity island were present in 13 O104 strains; however, the central portion of LEE was missing, where the eae gene is located. In O104:H4 strains, the missing central portion of the LEE locus was replaced by a pathogenicity island carrying the aidA (adhesin involved in diffuse adherence) gene and antibiotic resistance genes commonly carried on plasmids. Enteroaggregative E. coli-specific virulence genes and European outbreak O104:H4-specific stx2-encoding Escherichia P13374 or Escherichia TL-2011c bacteriophages were missing in some of the O104:H4 genome sequences available from public databases. Most of the genomic variations in the strains examined were due to the presence of different mobile genetic elements, including prophages and genomic island regions. The presence of plasmids carrying virulence-associated genes may play a role in the pathogenic potential of O104 strains.

CONCLUSIONS

The two strains sequenced in this study (O104:H21 and O104:H7) are genetically more similar to each other than to the O104:H4 strains that caused an outbreak in Germany in 2011 and strains found in Central Africa. A hypothesis on strain evolution and pathogenic potential of various H-serotypes of E. coli O104 strains is proposed.

Comparison of whole genome sequences from human and non-human Escherichia coli O26 strains

Shiga toxin-producing Escherichia coli (STEC) O26 is the second leading E. coli serogroup responsible for human illness outbreaks behind E. coli O157:H7. Recent outbreaks have been linked to emerging pathogenic O26:H11 strains harboring stx 2 only. Cattle have been recognized as an important reservoir of O26 strains harboring stx 1; however the reservoir of these emerging stx 2 strains is unknown. The objective of this study was to identify nucleotide polymorphisms in human and cattle-derived strains in order to compare differences in polymorphism derived genotypes and virulence gene profiles between the two host species. Whole genome sequencing was performed on 182 epidemiologically unrelated O26 strains, including 109 human-derived strains and 73 non-human-derived strains. A panel of 289 O26 strains (241 STEC and 48 non-STEC) was subsequently genotyped using a set of 283 polymorphisms identified by whole genome sequencing, resulting in 64 unique genotypes. Phylogenetic analyses identified seven clusters within the O26 strains. The seven clusters did not distinguish between isolates originating from humans or cattle; however, clusters did correspond with particular virulence gene profiles. Human and non-human-derived strains harboring stx 1 clustered separately from strains harboring stx 2, strains harboring eae, and non-STEC strains. Strains harboring stx 2 were more closely related to non-STEC strains and strains harboring eae than to strains harboring stx 1. The finding of human and cattle-derived strains with the same polymorphism derived genotypes and similar virulence gene profiles, provides evidence that similar strains are found in cattle and humans and transmission between the two species may occur.

Effect of proximity to a cattle feedlot on Escherichia coli O157:H7 contamination of leafy greens and evaluation of the potential for airborne transmission

The impact of proximity to a beef cattle feedlot on Escherichia coli O157:H7 contamination of leafy greens was examined. In each of 2 years, leafy greens were planted in nine plots located 60, 120, and 180 m from a cattle feedlot (3 plots at each distance). Leafy greens (270) and feedlot manure samples (100) were collected six different times from June to September in each year. Both E. coli O157:H7 and total E. coli bacteria were recovered from leafy greens at all plot distances. E. coli O157:H7 was recovered from 3.5% of leafy green samples per plot at 60 m, which was higher (P < 0.05) than the 1.8% of positive samples per plot at 180 m, indicating a decrease in contamination as distance from the feedlot was increased. Although E. coli O157:H7 was not recovered from air samples at any distance, total E. coli was recovered from air samples at the feedlot edge and all plot distances, indicating that airborne transport of the pathogen can occur. Results suggest that risk for airborne transport of E. coli O157:H7 from cattle production is increased when cattle pen surfaces are very dry and when this situation is combined with cattle management or cattle behaviors that generate airborne dust. Current leafy green field distance guidelines of 120 m (400 feet) may not be adequate to limit the transmission of E. coli O157:H7 to produce crops planted near concentrated animal feeding operations. Additional research is needed to determine safe set-back distances between cattle feedlots and crop production that will reduce fresh produce contamination.

Reducing assembly complexity of microbial genomes with single-molecule sequencing

BACKGROUND

The short reads output by first- and second-generation DNA sequencing instruments cannot completely reconstruct microbial chromosomes. Therefore, most genomes have been left unfinished due to the significant resources required to manually close gaps in draft assemblies. Third-generation, single-molecule sequencing addresses this problem by greatly increasing sequencing read length, which simplifies the assembly problem.

RESULTS

To measure the benefit of single-molecule sequencing on microbial genome assembly, we sequenced and assembled the genomes of six bacteria and analyzed the repeat complexity of 2,267 complete bacteria and archaea. Our results indicate that the majority of known bacterial and archaeal genomes can be assembled without gaps, at finished-grade quality, using a single PacBio RS sequencing library. These single-library assemblies are also more accurate than typical short-read assemblies and hybrid assemblies of short and long reads.

CONCLUSIONS

Automated assembly of long, single-molecule sequencing data reduces the cost of microbial finishing to $1,000 for most genomes, and future advances in this technology are expected to drive the cost lower. This is expected to increase the number of completed genomes, improve the quality of microbial genome databases, and enable high-fidelity, population-scale studies of pan-genomes and chromosomal organization.

Comparative genomics of the Campylobacter lari group

The Campylobacter lari group is a phylogenetic clade within the epsilon subdivision of the Proteobacteria and is part of the thermotolerant Campylobacter spp., a division within the genus that includes the human pathogen Campylobacter jejuni. The C. lari group is currently composed of five species (C. lari, Campylobacter insulaenigrae, Campylobacter volucris, Campylobacter subantarcticus, and Campylobacter peloridis), as well as a group of strains termed the urease-positive thermophilic Campylobacter (UPTC) and other C. lari-like strains. Here we present the complete genome sequences of 11 C. lari group strains, including the five C. lari group species, four UPTC strains, and a lari-like strain isolated in this study. The genome of C. lari subsp. lari strain RM2100 was described previously. Analysis of the C. lari group genomes indicates that this group is highly related at the genome level. Furthermore, these genomes are strongly syntenic with minor rearrangements occurring only in 4 of the 12 genomes studied. The C. lari group can be bifurcated, based on the flagella and flagellar modification genes. Genomic analysis of the UPTC strains indicated that these organisms are variable but highly similar, closely related to but distinct from C. lari. Additionally, the C. lari group contains multiple genes encoding hemagglutination domain proteins, which are either contingency genes or linked to conserved contingency genes. Many of the features identified in strain RM2100, such as major deficiencies in amino acid biosynthesis and energy metabolism, are conserved across all 12 genomes, suggesting that these common features may play a role in the association of the C. lari group with coastal environments and watersheds.

Investigation of bacterial diversity in the feces of cattle fed different diets

The objective of this study is to investigate individual animal variation of bovine fecal microbiota including as affected by diets. Fecal samples were collected from 426 cattle fed 1 of 3 diets typically fed to feedlot cattle: 1) 143 steers fed finishing diet (83% dry-rolled corn, 13% corn silage, and 4% supplement), 2) 147 steers fed late growing diet (66% dry-rolled corn, 26% corn silage, and 8% supplement), and 3) 136 heifers fed early growing diet (70% corn silage and 30% alfalfa haylage). Bacterial 16S rRNA gene amplicons were determined from individual fecal samples using next-generation pyrosequencing technology. A total of 2,149,008 16S rRNA gene sequences from 333 cattle with at least 2,000 sequences were analyzed. Firmicutes and Bacteroidetes were dominant phyla in all fecal samples. At the genus level, Oscillibacter, Turicibacter, Roseburia, Fecalibacterium, Coprococcus, Clostridium, Prevotella, and Succinivibrio were represented by more than 1% of total sequences. However, numerous sequences could not be assigned to a known genus. Dominant unclassified groups were unclassified Ruminococcaceae and unclassified Lachnospiraceae that could be classified to a family but not to a genus. These dominant genera and unclassified groups differed (P < 0.001) with diets. A total of 176,692 operational taxonomic units (OTU) were identified in combination across all the 333 cattle. Only 2,359 OTU were shared across 3 diet groups. UniFrac analysis showed that bacterial communities in cattle feces were greatly affected by dietary differences. This study indicates that the community structure of fecal microbiota in cattle is greatly affected by diet, particularly between forage- and concentrate-based diets.

Disinfectant and antibiotic susceptibility profiles of Escherichia coli O157:H7 strains from cattle carcasses, feces, and hides and ground beef from the United States

The disinfectant and antibiotic susceptibility profiles of 344 Escherichia coli O157:H7 strains from cattle carcasses, feces, and hides and ground beef from the United States were determined. A low prevalence of antibiotic resistance was observed (14%). The highest prevalences of resistance were to sulfisoxazole (10.5%), tetracycline (9.9%), streptomycin (7%), and chloramphenicol (4.9%). Four strains were resistant to eight antibiotics (two strains from ground beef and one strain each from hide and preevisceration carcass swabs of cull cattle at harvest). Pulsed-field gel electrophoresis analysis of the E. coli O157:H7 strains revealed two major groups (designated 1 and 2) composed of 17 and 20 clusters, respectively. Clusters 1A, 1B, 1C, and 1G.1 were associated with multidrug-resistant strains. There was no observed correlation between disinfectant resistance and antibiotic resistance. Sixty-nine (20%) of the 344 strains were resistant to chlorhexidine or benzalkonium chloride or the MICs of benzyldimethyldodecylammonium chloride were elevated. Inducible resistance was observed at elevated concentrations of antibiotics (1.4%) and disinfectants (6.1%). The highest rate of disinfectant inducible resistance was to OdoBan, quaternary ammonium chlorides, and the surface disinfectants F25, FS512, and MG, which are used in dairies, restaurants, and food processing plants. High MICs (1,024 to 4,096 m g/ml) of acetic, lactic, and citric acids were found. The decreasing order of acid potency based on molar MICs (MICs(molar)) was acetic, citric, and lactic acid. The correlation of the concentration of dissociated organic acids and MICs(molar) strongly suggests that the observed inhibition of E. coli O157:H7 was primarily due to dissociated forms of the acids.

Biofilm formation by Shiga toxin-producing Escherichia coli O157:H7 and Non-O157 strains and their tolerance to sanitizers commonly used in the food processing environment

Shiga toxin-producing Escherichia coli (STEC) strains are important foodborne pathogens. Among these, E. coli O157:H7 is the most frequently isolated STEC serotype responsible for foodborne diseases. However, the non-O157 serotypes have been associated with serious outbreaks and sporadic diseases as well. It has been shown that various STEC serotypes are capable of forming biofilms on different food or food contact surfaces that, when detached, may lead to cross-contamination. Bacterial cells at biofilm stage also are more tolerant to sanitizers compared with their planktonic counterparts, which makes STEC biofilms a serious food safety concern. In the present study, we evaluated the potency of biofilm formation by a variety of STEC strains from serotypes O157:H7, O26:H11, and O111:H8; we also compared biofilm tolerance with two types of common sanitizers, a quaternary ammonium chloride-based sanitizer and chlorine. Our results demonstrated that biofilm formation by various STEC serotypes on a polystyrene surface was highly strain-dependent, whereas the two non-O157 serotypes showed a higher potency of pellicle formation at air-liquid interfaces on a glass surface compared with serotype O157:H7. Significant reductions of viable biofilm cells were achieved with sanitizer treatments. STEC biofilm tolerance to sanitization was strain-dependent regardless of the serotypes. Curli expression appeared to play a critical role in STEC biofilm formation and tolerance to sanitizers. Our data indicated that multiple factors, including bacterial serotype and strain, surface materials, and other environmental conditions, could significantly affect STEC biofilm formation. The high potential for biofilm formation by various STEC serotypes, especially the strong potency of pellicle formation by the curli-positive non-O157 strains with high sanitization tolerance, might contribute to bacterial colonization on food contact surfaces, which may result in downstream product contamination.

Comparison of bacterial communities in faeces of beef cattle fed diets containing corn and wet distillers' grain with solubles

AIM

The mammalian intestinal microflora has been shown to impact host physiology. In cattle, intestinal bacteria are also associated with faecal contamination of environmental sources and human illness via foodborne pathogens. Use of wet distillers' grains with solubles (WDGS) in cattle feed creates a gastrointestinal environment where some bacterial species are enriched. Here, we examine if a diet containing 40% WDGS results in fundamentally different microbial community structures.

METHODS AND RESULTS

The 20,002 16S r-RNA gene sequences from 20 beef cattle were analysed using Sanger sequencing methods. At the genus level, Prevotella (Gram negative) and Anaerobacter (Gram positive) were the most frequently occurring bacteria in our beef cattle faecal samples. Diet-associated differences in prevalence were noted for Prevotella but not Anaerobacter.

CONCLUSIONS

Diet affects community structure. Faecal communities of co-housed beef cattle are not identical.

SIGNIFICANCE AND IMPACT OF THE STUDY

It is known that a diet of 40% corn-based WDGS increases the generic Escherichia coli in the faeces and enriches E. coli O157:H7. The results from the current study suggest that in addition to previously observed changes in E. coli, the entire bacterial community structure is different for animals fed 40% corn-based WDGS compared to a traditional corn-finishing diet.

Phylogeny of Shiga toxin-producing Escherichia coli O157 isolated from cattle and clinically ill humans

Cattle are a major reservoir for Shiga toxin-producing Escherichia coli O157 (STEC O157) and harbor multiple genetic subtypes that do not all associate with human disease. STEC O157 evolved from an E. coli O55:H7 progenitor; however, a lack of genome sequence has hindered investigations on the divergence of human- and/or cattle-associated subtypes. Our goals were to 1) identify nucleotide polymorphisms for STEC O157 genetic subtype detection, 2) determine the phylogeny of STEC O157 genetic subtypes using polymorphism-derived genotypes and a phage insertion typing system, and 3) compare polymorphism-derived genotypes identified in this study with pulsed field gel electrophoresis (PFGE), the current gold standard for evaluating STEC O157 diversity. Using 762 nucleotide polymorphisms that were originally identified through whole-genome sequencing of 189 STEC O157 human- and cattle-isolated strains, we genotyped a collection of 426 STEC O157 strains. Concatenated polymorphism alleles defined 175 genotypes that were tagged by a minimal set of 138 polymorphisms. Eight major lineages of STEC O157 were identified, of which cattle are a reservoir for seven. Two lineages regularly harbored by cattle accounted for the majority of human disease in this study, whereas another was rarely represented in humans and may have evolved toward reduced human virulence. Notably, cattle are not a known reservoir for E. coli O55:H7 or STEC O157:H⁻ (the first lineage to diverge within the STEC O157 serogroup), which both cause human disease. This result calls into question how cattle may have originally acquired STEC O157. The polymorphism-derived genotypes identified in this study did not surpass PFGE diversity assessed by BlnI and XbaI digestions in a subset of 93 strains. However, our results show that they are highly effective in assessing the evolutionary relatedness of epidemiologically unrelated STEC O157 genetic subtypes, including those associated with the cattle reservoir and human disease.

Bacterial community analysis of beef cattle feedlots reveals that pen surface is distinct from feces

The surface of beef cattle feedlot pens is commonly conceptualized as being packed uncomposted manure. Despite the important role that the feedlot pen may play in the transmission of veterinary and zoonotic pathogens, the bacterial ecology of feedlot surface material is not well understood. Our present study characterized the bacterial communities of the beef cattle feedlot pen surface material using 3647 full-length 16S rDNA sequences, and we compared the community composition of feedlot pens to the fecal source material. The feedlot surface composite was represented by members of the phylum Actinobacteria (42%), followed by Firmicutes (24%), Bacteroidetes (24%), and Proteobacteria (9%). The feedlot pen surface material bacterial communities were clearly distinct from those of the feces from animals in the same pen. Comparisons with previously published results of feces from the animals in the same pen reveal that, of 139 genera identified, only 25 were present in both habitats. These results indicate that, microbiologically, the feedlot pen surface material is separate and distinct from the fecal source material, suggesting that bacteria that originate in cattle feces face different selection pressures and survival challenges during their tenure in the feedlot pen, as compared to their residence in the gastrointestinal tract.

Virulence-associated and antibiotic resistance genes of microbial populations in cattle feces analyzed using a metagenomic approach

The bovine fecal microbiota impacts human food safety as well as animal health. Although the bacteria of cattle feces have been well characterized using culture-based and culture-independent methods, techniques have been lacking to correlate total community composition with community function. We used high throughput sequencing of total DNA extracted from fecal material to characterize general community composition and examine the repertoire of microbial genes present in beef cattle feces, including genes associated with antibiotic resistance and bacterial virulence. Results suggest that traditional 16S sequencing using "universal" primers to generate full-length sequence may under represent Acitinobacteria and Proteobacteria. Over eight percent (8.4%) of the sequences from our beef cattle fecal pool sample could be categorized as virulence genes, including a suite of genes associated with resistance to antibiotic and toxic compounds (RATC). This is a higher proportion of virulence genes found in Sargasso sea, chicken cecum, and cow rumen samples, but comparable to the proportion found in Antarctic marine derived lake, human fecal, and farm soil samples. The quantitative nature of metagenomic data, combined with the large number of RATC classes represented in samples from widely different habitats indicates that metagenomic data can be used to track relative amounts of antibiotic resistance genes in individual animals over time. Consequently, these data can be used to generate sample-specific and temporal antibiotic resistance gene profiles to facilitate an understanding of the ecology of the microbial communities in each habitat as well as the epidemiology of antibiotic resistant gene transport between and among habitats.