Next Article in Journal
Synthetic Hexanucleotides as a Tool to Overcome Excessive Neutrophil Activation Caused by CpG-Containing Oligonucleotides
Next Article in Special Issue
Prevalence of Salmonella spp. and Staphylococcus aureus in Chicken Meat and Pork from Cambodian Markets
Previous Article in Journal
Unique Attributes of the Laurel Wilt Fungal Pathogen, Raffaelea lauricola, as Revealed by Metabolic Profiling
Previous Article in Special Issue
Virulence Characterization of Listeria monocytogenes, Listeria innocua, and Listeria welshimeri Isolated from Fish and Shrimp Using In Vivo Early Zebrafish Larvae Models and Molecular Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genomic Characterization of Salmonella typhimurium DT104 Strains Associated with Cattle and Beef Products

1
Produce Safety and Microbiology Research Unit, Western Regional Research Center, Agricultural Research Service, U.S. Department of Agriculture, Albany, CA 94710, USA
2
Department of Biology, California State University-Northridge, Northridge, CA 91330, USA
3
School of Animal and Comparative Biomedical Sciences, University of Arizona, Tucson, AZ 85721, USA
*
Author to whom correspondence should be addressed.
Pathogens 2021, 10(5), 529; https://doi.org/10.3390/pathogens10050529
Submission received: 17 March 2021 / Revised: 4 April 2021 / Accepted: 22 April 2021 / Published: 27 April 2021

Abstract

:
Salmonella enterica subsp. enterica serovar Typhimurium DT104, a multidrug-resistant phage type, has emerged globally as a major cause of foodborne outbreaks particularly associated with contaminated beef products. In this study, we sequenced three S. Typhimurium DT104 strains associated with a 2009 outbreak caused by ground beef, including the outbreak source strain and two clinical strains. The goal of the study was to gain a stronger understanding of the genomics and genomic epidemiology of highly clonal S. typhimurium DT104 strains associated with bovine sources. Our study found no single nucleotide polymorphisms (SNPs) between the ground beef source strain and the clinical isolates from the 2009 outbreak. SNP analysis including twelve other S. typhimurium strains from bovine and clinical sources, including both DT104 and non-DT104, determined DT104 strains averaged 55.0 SNPs between strains compared to 474.5 SNPs among non-DT104 strains. Phylogenetic analysis separated the DT104 strains from the non-DT104 strains, but strains did not cluster together based on source of isolation even within the DT104 phage type. Pangenome analysis of the strains confirmed previous studies showing that DT104 strains are missing the genes for the allantoin utilization pathway, but this study confirmed that the genes were part of a deletion event and not substituted or disrupted by the insertion of another genomic element. Additionally, cgMLST analysis revealed that DT104 strains with cattle as the source of isolation were quite diverse as a group and did not cluster together, even among strains from the same country. Expansion of the analysis to 775 S. typhimurium ST19 strains associated with cattle from North America revealed diversity between strains, not limited to just among DT104 strains, which suggests that the cattle environment is favorable for a diverse group of S. typhimurium strains and not just DT104 strains.

1. Introduction

Salmonella is a major foodborne pathogen, annually resulting in 1.4 million cases in the United States [1] and over 90 million on a global scale [2]. Salmonella enterica subsp. enterica serovar Typhimurium (Salmonella typhimurium) was the most common serovar fifteen years ago, and it is currently the third most common serovar in the United States [3] with many different types of food products, including beef products, associated with outbreaks. In fact, beef products are a critically important vehicle for S. typhimurium transmission to humans. Ground beef in particular has been implicated in an estimated 45% of all U.S. Salmonella outbreaks during 2002–2011 [4]. Since 1973, there have been 16 S. typhimurium outbreaks associated with beef products in the U.S., with ground beef accounting for 37.5% of those outbreaks [4]. While beef consumption in the U.S. is four times higher than the world average [5], this does suggest surveillance of beef products for S. typhimurium is critically important to the safety of consumers on a global scale.
S. typhimurium is highly clonal, and several types of subtyping methods, including pulsed field gel electrophoresis (PFGE), are not always effective at distinguishing strains [6], which prior to whole genome sequencing (WGS) had made epidemiological investigations more complicated. Several retrospective studies conducted using the higher resolution of WGS have begun to improve the genomic epidemiology of S. typhimurium, although most of these studies have concentrated on using SNP analysis to compare the strains [7,8,9,10]. In the examination of genomes from five different outbreaks, one study found that for four outbreaks, epidemiologically linked strain only varied by one or two SNPs. In a fifth outbreak, the epidemiologically linked strains varied by up to 12 SNPs [9]. Whereas, another retrospective study of S. typhimurium DT8 strains from an outbreak found only three SNP differences between outbreak-linked strains, but those unrelated to the outbreak had up to 342 SNPs [11]. As the genome epidemiology of Salmonella continues to be investigated and expanded, it becomes clear that detailed studies of different species, phage types, sources, and many other variables need to be examined to truly understand the genomic epidemiology of Salmonella on a global scale.
S. typhimurium phage type DT104 has emerged as a common phage type accounting for up to 50% of all the S. typhimurium cases in many countries [12]. DT104 strains have also been linked with numerous outbreaks around the world associated with beef products [13,14,15], including ground beef [16,17]. S. typhimurium DT104 is a multidrug-resistant definitive phage type that first emerged in the 1980s [18], and commonly has resistance to ampicillin, chloramphenicol, spectinomycin, streptomycin, sulphonamides, and tetracyclines [19]. Studies have shown that S. typhimurium DT104 strains cannot be distinguished by pulsed field gel electrophoresis (PFGE) [20], but multiple locus variable number of tandem repeat analysis (MLVA) has shown the ability to discriminate some strains of DT104 [21]. The advent of whole genome sequencing (WGS) allows for resolution down to the individual nucleotide, which has allowed for expansion of the global genomic epidemiology of DT104 strains. Two major studies have explored the genomics of DT104 strains from different hosts [22] and temporally and spatially different strains [23], but currently no studies have specifically examined DT104 associated with beef products. As cattle and beef products are a critical source of global DT104 infections, as well as S. typhimurium, it is vital to gain a better understanding of the global genomic epidemiology.
In 2009, S. typhimurium DT104-contaminated ground beef was determined to be the cause of an outbreak in Colorado, USA, leading to 14 clinical cases and the recall of over 466,000 pounds of ground beef [16]. The Colorado Department of Public Health and Environment (CDPHE) isolated several clinical isolates and an additional isolate from a ground beef patty obtained from the home of a patient. The goals of this study were to (1) sequence and compare two clinical isolates and ground beef source strains; (2) conduct a detailed genomic comparison of S. typhimurium clinical and bovine-associated DT104 and non-DT104 strains; (3) compare the ground beef outbreak strains against a global DT104 population; (4) compare these bovine-associated DT104 strains to other bovine-associated S. typhimurium from North America, which will provide a detailed genomic characterization of DT104 strains associated with cattle and beef products.

2. Results

2.1. 2009 Colorado Ground Beef Outbreak DT104 Strains

The complete genome of S. typhimurium DT104 strain RM9437 (ground beef) from the 2009 Colorado outbreak is 5,030,429 bp (composed of 4,936,499 bp chromosome and 93,930 bp plasmid), and encodes for 4569 CDSs, 7 rRNA operons, and 83 tRNAs. The genome contains eight prophages, including six complete and two incomplete prophages (Table 1). The plasmid pRM9437 is >99% identical to pLT2, the described virulence plasmid found in many S. typhimurium strains. It possesses spvRABCD and pefBACD (fimbriae) and srgABC (involved in fimbriae biogenesis). This version of the plasmid does not possess antibiotic resistance genes. The strain also possesses the Salmonella Genomic Island 1 (SGI1), which is the multidrug resistance element possessing resistance genes to ampicillin, chloramphenicol, streptomycin, sulfonamides, and tetracycline (ACSSuT resistance) on the chromosome and described previously [24,25]. Based on the sequences of seven housekeeping genes used for multilocus sequence typing (MLST) of S. enterica, the sequence type for these outbreak strains was determined to be ST19. Two clinical strains from the outbreak, RM9435 (clinical) and RM9436 (clinical), were used for manual high-quality SNP curation between clinical strains and ground beef source strain from the same outbreak, which revealed that there are no SNPs between these three strains. Therefore, all additional SNP analysis only involved RM9437.

2.2. SNP Analysis between S. typhimurium DT104 and Non-DT104 Strains

Thirteen S. typhimurium ST19 strains with complete genomes were selected for further SNP analysis (Table 1). These strains included seven from bovine sources (three DT104 and four non-DT104) and six from clinical sources (four DT104 and two non-DT104). Using CSI Phylogeny SNP analysis software and RM9437 as the reference genome, we determined that the seven DT104 strains had significantly less SNPs compared to the six non-DT104 strains (Figure 1). This finding was strengthened by performing SNP analysis of all the strains against each other, which found a total of 1781 high-quality SNPs between all strains with an average of 440.9 SNPs between any two strains. However, the DT104 strains possessed an average of 55.0 SNPs compared to 474.5 SNPs for non-DT104 strains (Table 2). Phylogenetic analysis of the thirteen ST19 strains reinforced this separation of the DT104 and non-DT104 strains, as all of the DT104 strains clustered together, and demonstrated a common ancestor deviation from non-DT104 strains (Figure S1). Yet the source of isolation did not have a similar impact on the relationship of the S. typhimurium strains regardless of the phage type, as neither bovine nor clinical strains clustered together on the maximum-likelihood (ML) tree even within the separate phage type clusters. While it would be expected that clinical strains would not cluster together unless they had a common infectious host source, it was interesting that the cattle strains did not cluster together, which suggests there is diversity among S. typhimurium strains associated with cattle. Furthermore, the bovine isolates did have significantly more SNPs between the strains than the clinical strains, as the bovine strains averaged 511.8 SNPs compared to an average of 344.9 SNPs for the clinical strains (Table 2).

2.3. Gene Content Analysis of S. typhimurium DT104 and Non-DT104 Strains

Core genome analysis of the thirteen strains found that there were 4300 core genes, 877 shell genes (present in 1 ≤ strains > 12), and 587 cloud genes (present in only one strain) among all strains, which was fairly consistent for non-DT104 strains (4331 core genes), bovine strains (4308 core genes), and clinical strains (4350 core genes). However, DT104 strains as a group had appreciably more core genes at 4512 compared to these other groups, which also resulted in fewer shell genes (300) and cloud genes (49) compared to the rest of the overall group (Supplementary Table S1). Fifty-three (53) genes unique to all the DT104 strains were identified in the pan-GWAS analysis, which were composed entirely of genes from prophages, although based on the PHASTER analysis there was some variation in the specific prophage that contained the genes between the different DT104 strains. For example, the 53 DT104 unique genes in S. typhimurium strain DT104 included 38 genes that were present in an intact Aeromonas phage vB AsaM-56 (as identified by PHASTER analysis), and the other 15 genes were present in an intact Salmonella phage ST64B (as identified by PHASTER analysis; Supplementary Table S2), whereas the 53 unique genes in S. typhimurium str. RM9437 were broken down as 38 genes in an intact Edwardsiella phage GF-2 (as identified by PHASTER) and the other 15 genes were present in an intact Salmonella phage 118970 sal3 (as identified by PHASTER). The difference in identified intact prophages containing the DT104 unique genes is likely due to PHASTER analysis, but could be worth investigating in the future. On the other hand, none of the 12 genes found to be unique to the non-DT104 strains or missing from DT104 strains were associated with prophages, and the genes included 11 genes in the allantoin–glyoxylate metabolic pathway and one hypothetical protein that is 87 amino acids long (Supplementary Table S3). These 11 allantoin–glyoxylate metabolic pathway genes seem to have been deleted in the DT104 strains as there is nothing inserted or replacing the genes in the genomes (Supplementary Figure S2A). No unique genes were identified for either the bovine strains or the clinical strains during the pan-GWAS analysis.

2.4. Prophage and Antibiotic Resistance Analysis

Phage type DT104 strains are considered a significant public health concern as they are frequently multidrug resistant, therefore the intact prophage profile and antibiotic resistance gene profile for all the strains were examined. Similar to the SNP analysis, DT104 and non-DT104 strains were found to have very different profiles for both the prophages and antibiotic resistance genes (Figure 2). Yet, the phage Gifsy-2 was common to all of the strains in the study and 11/13 (84.6%) possessed the phage Gifsy-1. Additionally, 5/7 (71.4%) of the DT104 strains contained the Edwardsiella phage GF-2 (defined by PHASTER), but none of the non-DT104 strains had this prophage. The Enterobacteria phage ST104, which is commonly present in DT104 strains [28], was also found in 6/7 (85.7%) of the DT104 strains. Interestingly, prophage ST104 was completely missing (not incomplete or questionable either) from the DT104 strain CDC 2011K-1702, but the prophage was present in the non-DT104 strain CDC 2010K-1587. Further analysis shows that CDC 2011K-1702 has the prophage Enterobacteria phage P4 inserted at the site of Enterobacteria phage ST104 insertion in the other DT104 strains. Nucleotide alignment of the prophage ST104 sequence from CDC 2010K-1587 (non-DT104) against the prophage sequences from DT104 strains (Figure S2B) showed only 67.44% identity. Alignment of the prophage ST104 sequences from DT104 strains against each other showed 99.99% identity. Bovine strains also had quite a bit of prophage variation, as none of the bovine strains shared an identical prophage pattern. There were 13 different prophages identified between the seven bovine strains, but only two prophages were shared by all the strains. The other 11 prophages were present in at less than 50% of the strains. None of the clinical strains had identical prophage patterns. The six clinical strains had a total of 10 different prophages identified, but only the phage Gifsy-2 was present in all the clinical strains.
While there was some diversity in the antibiotic resistance genes present in the different S. typhimurium genomes, it was not as distinctive as the prophage diversity among the strains. In fact, all 13 of the strains contained the aac(6′)-laa gene that encodes for aminoglycoside resistance (Figure 2). According to the Comprehensive Antibiotic Resistance Database (CARD) [29], this gene is present in 91.3% of Salmonella enterica genomes. However, there were some differences in antibiotic resistance genes between DT104 strains and non-DT104 strains. For example, none of the non-DT104 strains contained the previously described SGI1 encoding ACSSuT resistance. For the DT104 strains, 5/7 (71.4%) had the entire genomic island. As a group, the bovine strains did not share an antibiotic profile. This was also the case for the clinical strains.

2.5. Global Genomic Epidemiology of DT104 Strains

We performed hierarchical clustering (HC) of the core genome multilocus sequence types (cgMLST), which provides a relationship between the three Colorado outbreak strains (RM9435, RM9436, and RM9437) and a global set of DT104 strains at different HC levels. At the hierarchal clustering 100 (HC100) level, the differences between strains’ cgMLST alleles are 100 or fewer, and this level has been associated with strains from long-term endemic persistence [30]. At this level, DT104 strains from bovine, swine, and food sources were distributed across the minimum spanning tree (MST) based on the cgMLST. These strains were also intermixed with human clinical strains across the entire MST, whereas DT104 strains from avian sources formed a specific cluster with only a few clinical strains and a single bovine strain. The U.S. outbreak strains clustered together as an individual node and away from most other bovine DT104 strains, but the clinical and ground beef forming a single node confirms the outbreak relationship of these strains. The country of origin does not seem to impact the clustering of the DT104 strains even for similar sources of isolation (Figure 3). To examine the relationship of these global DT104 strains a little closer, strains RM9435, RM9436, and RM9437, 25 bovine strains, and a random selection of 72 strains representing other sources of isolation were selected for pangenome analysis using Anvi’o software. The analysis found 4362 core gene clusters, 621 accessory gene clusters, and 424 singleton gene clusters among all 100 strains. Similar to cgMLST analysis, when clustering the strains based on presence/absence of the gene clusters, the swine and bovine strains did not all cluster together and were intermixed with the clinical strains (Figure 4). Examination of genes present or absent from DT104 strains associated with cattle found that there were no genes present in only cattle-associated DT104 strains and absent in non-cattle DT104 strains or vice versa. However, there were 11 genes that were statistically significant in presence or absence from cattle-associated DT104 strains when compared to non-cattle-associated strains. Only two genes were more commonly present in cattle DT104 strains than in non-cattle-associated DT104 strains. The beta-lactamase PSE-1 (pse1) gene was present in all 26 bovine DT104 strains examined, but only present in 75.0% of the 72 non-cattle-associated DT104 strains (p-value = 0.0026). The tRNA-Ala gene was present in 19.2% of cattle-associated DT104 strains compared to only 1.4% of non-cattle-associated DT104 strains (p-value = 0.0047). Whereas, the other nine genes were statistically more likely to be missing from cattle-associated DT104 strains than from non-cattle DT104 strains, including five hypothetical genes encoding proteins, the regulatory protein Rop, mobilization protein MbeC, endoribonuclease ToxN, and dihydropteroate synthase (Table 3).

2.6. Genomic Epidemiology of Cattle-Associated S. typhimurium Strains

To further examine the genomic epidemiology of bovine-associated S. typhimurium strains, the fifteen ST19 S. typhimurium strains were compared to both bovine-associated DT104 and non-DT104 ST19 S. typhimurium strains from North America. The MST based on the cgMLST analysis demonstrated quite a bit of diversity among the various bovine strains as there were not many large nodes. However, this analysis did separate the fifteen DT104 and non-DT104 strains from each other. All DT104 strains branched off the same major node and away from six non-DT104 strains (Figure 5). The four non-DT104 bovine strains were also spread nearly across the entire MST. Overall, the results suggest that S. typhimurium strains associated with different bovine sources are extremely diverse even within a single geographical region, particularly as S. typhimurium is a highly clonal group as a whole.

3. Discussion

S. typhimurium strains are highly clonal, and many traditional subtyping methods such as pulsed field gel electrophoresis (PFGE) have trouble resolving strains during an outbreak [6]. This was particularly true with DT104 strains, as prior to the advent of whole genome sequencing (WGS) only MLVA was shown to have some resolution to differentiate DT104 strains [21,31]. Using WGS, this study investigated the level of clonality of DT104 strains within an outbreak, compared bovine and clinical DT104 strains to non-DT104 strains from within the same sequence type (ST19), additional DT104 strains from around the world, and cattle-associated S. typhimurium strains from North America. S. typhimurium DT104 strains are a common cause of disease in many parts of the world, and improving the genomic epidemiology of this phage type of S. typhimurium provides input information for public health laboratories and clarifies the role cattle and beef products may play in human infection. The SNP analysis between source and clinical strains from an outbreak revealed no SNP differences, which was confirmed at the allele level via cgMLST analysis. This suggests that due to the high level of clonality, there will be a limited number of SNPs between clinical DT104 strains during an outbreak, but further outbreak investigations are needed to confirm these results. However, this information can be used for improving source attribution of DT104 strains, and ultimately to improve source-tracking to assist in rapidly conducting future outbreak investigations to minimize the public’s exposure to contaminated food products.
The advent of WGS has dramatically improved the ability of public health agencies around the world to source track foodborne outbreaks [32]. Nevertheless, the continued development and improvement of WGS as a critical public health tool for outbreak source tracking for pathogens is vital to improving global food safety, particularly understanding SNP variation between different serovars and phage types of Salmonella. SNP analysis of the bovine and clinical strains within the ST19 group averaged 440.9 SNPs. That average was slightly higher than was found in another study (380.6 SNPs) [33]; however, that study examined fewer strains than our study. Furthermore, our study found that DT104 strains had significantly less SNPs than non-DT104 ST19 strains as a whole, with a range of 30 to 85 and an average of 55.0 SNPs between strains. Non-DT104 strains had a range of 10 to 646 and an average of 474.5 SNPs. Mather et al. found that out of 359 DT104 strains, there was a range of 0 to 167 SNPs between the strains [22], which even on the high end of the range is significantly lower than that of most non-DT104 strains from the ST19 group. Taken together, these data provide genomic evidence that DT104 strains have a higher level of clonality when compared with other ST19 S. typhimurium strains. Thus, they may require special consideration during analysis for outbreak source tracking.
The detailed analysis of DT104 strains from bovine sources and clinical sources against non-DT104 strains found that not all DT104 strains contained the SGI1 genomic island that encodes multidrug resistance, which has been described before for certain DT104 strains. For example, DT104 strains have been described before with a partial deletion of SGI1 that lacked flo and tet genes [34], which is similar to the DT104 strain CDC 2011K-1702 that was used in this study. Although, we also found that DT104 strain SA972816 from China was missing the entire 13 kb multidrug resistance section of the SGI1 genomic island, which is not as common in DT104 strains, but has been described previously [23,35]. The pangenome analysis using Roary/Scoary of the 98 DT104 strains to determine DT104 cattle unique genes found that all the cattle-associated DT104 strains (DT104 cattle-associated strain SA972816 was not included in this particular analysis) had the beta-lactamase PSE-1, but it was missing in 25.0% of the non-cattle-associated DT104 strains. The pse1 or blaP1 gene is in the SGI1 genomic island and encodes for ampicillin resistance [36], which suggests that 25.0% of the non-cattle-associated DT104 strains could be missing part if not all of the SGI1 genomic island. However, it seems much more common for cattle-associated DT104 strains to have the SGI1 genomic island that encodes for multidrug resistance, which might be due to the heavy use of antimicrobials in the beef industry, but would need further confirmation. We investigated the Enterobacteria phage ST104 from several DT104 strains as this is the phage that is believed to provide the DT104 phage type [28], and it remains unclear as to how DT104 strains CDC 2011K-1702 could be missing the phage and still be a DT104 strain. Several possibilities exist: (1) it is not a DT104 strain and was mis-identified by the CDC, but it did group with the other DT104 strains in the phylogenetic analysis although it was also missing the SGI1 genomic island; (2) multiple bacteriophage that are closely related could be responsible for the DT104 phage type. Currently, additional analysis is needed to determine the answer to this issue. However, the reason for the non-DT104 strain carrying the ST104 phage was revealed to be that there was only ~67% sequence homology, which suggests the phage may not be functional in this strain, thus impacting the phage type.
It has previously been described that the allantoin utilization loci is deleted in DT104 strains [37,38]. Yet, the previous types of analysis were not able to determine if the genes were disrupted due to insertion, deletion, or substitution in the DT104 genomes. Therefore, this study is the first to show that the allantoin utilization loci was deleted from the DT104 genome, which suggests it might not be needed for DT104 strains during colonization of common isolation sources like cattle and swine. In fact, Matiasovicova et al. investigated if deletion of allantoin utilization loci in S. typhimurium affected virulence in mice or colonization in chickens, since poultry have high levels of allantoin. The study found that the loss of the loci S. typhimurium does not appear to affect virulence in mice or colonization of chickens [38], so this may not be the reason for the rarity of DT104 strains in poultry. Thus, the reason for the deletion of the allantoin utilization loci in DT104 strains, whether in colonization, survival, virulence, or for another reason, still needs further investigation.
DT104 strains are commonly isolated from swine and cattle, yet not as frequently from poultry [39], and the results of this study demonstrate that it could be due to more diversity in the DT104 strains associated with cattle and swine compared to those associated with poultry. The cgMLST analysis of DT104 strains from EnteroBase found that those strains isolated from poultry were clustered together compared to cattle and swine strains. Cattle and swine strains were distributed across the MST, suggesting a significant amount of diversity among these strains. It should be noted that all of the poultry strains were isolated in Ireland, which could explain the lack of diversity in the poultry-associated strains. However, cattle strains as strains restricted to specific European regions (Germany, Ireland, and Northern Ireland) were spread across the tree and not clustered together by country. The diversity observed in this study was at the HC100 level that is associated with long-term endemic persistence in a host [30], therefore it suggests that cattle and swine are endemically colonized with DT104 strains for long time periods, which results in more frequent isolation and higher risk of human infection. Overall, we suggest that DT104 strains have established a niche in cattle and swine populations around the world that has resulted in numerous outbreaks associated particularly with contaminated beef products. Interestingly, analysis of DT104 strains to determine unique genes associated with cattle colonization found that there were no unique genes associated only with DT104 cattle strains, which suggests that DT104 strains colonizing cattle are more generalist and not specialist, thus probably resulting in increased diversity of the population associated with cattle across the world. Furthermore, when ST19 lineage S. typhimurium strains associated with cattle, including multiple DT104 strains, from a defined region (North America) were examined using cgMLST, there was also quite a bit of diversity in general among S. typhimurium strains. Additionally, the diversity of the cattle-associated DT104 strains was further supported as the multiple DT104 cattle strains were spread across the tree and not together in a tight cluster, even compared to non-DT104 strains.
In conclusion, this study was the first to explore the genomics of DT104 strains strictly associated with cattle and/or beef products from the perspective of understanding the genomics and genomic epidemiology for improving global food safety. The results of the study found that DT104 strains are more clonal than non-DT104 S. typhimurium strains even from the same ST, but this changes when the DT104 strains are associated with cattle. In fact, an association with cattle and/or beef products seems to increase the diversity of S. typhimurium strains in general even within limited geographical areas. Therefore, from a global food safety perspective, it should be understood that DT104 and non-DT104 S. typhimurium strains associated with cattle and/or beef products can be quite diverse. However, DT104 strains from outbreaks should have clinical strains and source strains with practically no SNPs, like that demonstrated for the 2009 Colorado outbreak strains in this study.

4. Materials and Methods

4.1. Strains

All S. typhimurium strains sequenced (RM9435, RM9436, and RM9437) and/or used for the detailed comparative genomics or SNP analysis in this study are listed in Table 1 along with their respective accession number [22,26,27]. All strains sequenced in this study were grown on LB agar and/or LB broth at 37 °C unless otherwise stated and were sequenced as described below. Additionally, all global S. typhimurium strains utilized for cgMLST analysis in EnteroBase (http://enterobase.warwick.ac.uk/; accessed 11 February 2021) [30,40] including DT104 and North American bovine-associated strains were available on 11 February 2021 and 14 May 2020, respectively. Additionally, those strains listed above and the 100 DT104 strains selected for pangenome/DT104 cattle-specific analysis are also listed in Supplementary Table S4.

4.2. DNA Extraction

Genomic DNA for the three S. typhimurium strains (RM9435, RM9436, and RM9437) that were sequenced in this study were prepared from single colonies grown overnight in LB broth, and then extracted as previously described [41]. Briefly, bacteria were lysed with lysozyme (20 mg/mL in 50 mM Tris, pH 8.0; Fisher Scientific) and 10% SDS and then sequentially treated with RNase A (1 mg/mL at 37 °C for 24 h; Fisher Scientific, Waltham, MA, USA) and proteinase K (20 mg/mL at 37 °C for 4 h; Fisher BioReagents, Waltham, MA, USA). Extracted DNA was then precipitated with a sodium acetate/ethanol solution, purified by phenol/chloroform extraction, precipitated again with ethanol, and finally, re-suspended in Buffer EB (Qiagen; Hilden, Germany) for DNA sequencing library preparation. Prior to DNA library preparation, the DNA concentration was determined using a Qubit 3 fluorometer (Thermo-Scientific; Waltham, MA, USA) measuring at 485/530 nm and quality-assessed using a Nanodrop ND-1000 spectrophotometer (Thermo-Scientific; Waltham, MA, USA) measuring at 260/280 and 260/230 ratio wavelengths.

4.3. Genome Sequencing

DNA sequencing libraries were prepared using the KAPA Low-Throughput Library Preparation Kit with Standard PCR Amplification Module (Kapa Biosystems; Wilmington, MA, USA) following the manufacturer’s instructions except for the following modifications: 750 ng of DNA was sheared at 30 psi for 40 s and size selected to 700–770 bp following Illumina MiSeq protocols. Standard desalted TruSeq LT and PCR primers (Integrated DNA Technologies (IDT); Coralville, IA, USA) were used at 0.375 and 0.5 µM final concentrations, respectively. The PCR reaction was reduced to 3–5 cycles. Libraries were quantified using the KAPA Library Quantification Kit (Kapa Biosystems; Wilmington, MA, USA), except using a 10 µL volume and 90 s annealing/extension PCR, then pooled and normalized to 4 nM. Pooled libraries were re-quantified by digital droplet PCR (ddPCR) on a QX200 system (Bio-Rad; Hercules, CA, USA), using the Illumina TruSeq ddPCR Library Quantification Kit and following manufacturer’s protocols, except with an additional 2 min annealing/extension time. All libraries were sequenced on an Illumina MiSeq instrument (500-cycle v2 sequencing kit; Illumina; San Diego, CA, USA) following the manufacturer’s protocols. Sequencing resulted in a total of 2,477,178 paired-end (PE) 250 bp reads for RM9435 (~137× coverage); 1,265,954 PE 250 bp reads for RM9436 (~70× coverage), and 1,290,788 PE 250 bp reads for RM9437 (~71× coverage).

4.4. Assembly and Annotation

The PE sequence reads for S. typhimurium str. RM9437 were assembled de novo using the Roche Newbler assembler (v2.3) using default parameters, which resulted in a total of 70 contigs. Next, the 70 contigs for RM9437 were ordered by aligning against S. typhimurium str. DT104 (HF937208.1 and HF937209.1) using Mauve software [42]. The genome was closed in silico using the following steps: (1) identification of repeated contigs using the Perlscript Contig_extender3 [43] and determination of the number of reads per contigs length, and (2) contig extension using Illumina reads and Newbler contigs within Geneious software (v9.1) [44], which resulted in a single chromosome and a single plasmid sequence for S. typhimurium str. RM9437.
All coding DNA sequences (CDSs) were determined using the Glimmer 3.0 [45] plugin for Geneious software, and then predicted CDSs were compared against annotations for S. typhimurium LT2 genome (AE006468.1) and S. typhimurium DT104 (HF937208.1). Novel or unclear annotations were identified by using BLASTP against the National Center for Biotechnology Information (NCBI) nonredundant (nr) database with the parameters ≥80% identity across ≥75% of the query sequence for a confirmation hit. All tRNAs and rRNA loci were identified by using BLASTN again the S. typhimurium str. LT2 genome, and the final annotations were confirmed by comparing against the automated NCBI Prokaryotic Genome Annotation Pipeline (PGAP) annotation [46]. Reads for strains RM9435 and RM9436 were individually assembled in EnteroBase automatically upon submission, and both strains were only used for core genome multilocus sequence typing (cgMLST) and SNP analysis against RM9437 in the study.

4.5. Single Nucleotide Polymorphism (SNP) Analysis

SNP variations between S. typhimurium strains RM9435, RM9436, and RM9437 were identified by mapping the RM9435 or RM9436 reads against the reference genome of S. typhimurium str. RM9437 using Geneious Prime 2020.1.2 software. Variation/SNPs were found using the following parameters: minimum 30× coverage, >80% of reads contain variation, maximum variant p-value of 10−6, and minimum strand bias p-value of 10−5. Repeat regions were identified by BLASTN search of all SNP regions against the RM9437 chromosome. Those with multiple hits (≥50% identity across ≥50% of the query sequence) were examined to identify repeat regions, and any additional repeat regions were identified using the RepeatFinder plugin in Geneious Prime 2012.1.2 software using the following parameters: 100 bp minimum repeat length and 1% maximum mismatches. All repeat regions were then removed from the results of the SNP analysis.
The SNP matrix of the twelve complete genomes of the ST19 S. typhimurium strains lacking sequence reads was generated using CSI Phylogeny 1.4 [47] (https://cge.cbs.dtu.dk/services/CSIPhylogeny/; accessed on 14 May 2020) from the Center for Genomic Epidemiology. The complete genome sequence of RM9437 was included in the analysis as a reference genome to generate the SNP matrix, using the following parameters: 30× minimum depth at SNP positions; 25% minimum relative depth at SNP positions; 100 bp minimum distance between SNPs (pruning); 50 minimum SNP quality; 25 minimum read mapping quality; 1.96 minimum Z-score.
SNP analysis was also conducted among the 13 different bovine and clinical S. typhimurium strains using Parsnp software (v1.2) [48], and the SNP density of strains visualized using Gingr software (v1.2) as part of the Harvest suite [48]. A maximum-likelihood tree was generated by using the 1781 filtered high-quality SNPs identified by the program, which were concatenated and aligned using command line MUSCLE software (v3.8.31) with the default parameters [49,50], and a maximum-likelihood (ML) tree generated with MEGAX software (v10.1.8) [51]. First, a best-fit model was determined for the alignment data, then a maximum-likelihood (ML) tree was generated using the Tamura-Nei model [52] with Gamma distribution with 1000 pseudoreplicates.
Identified SNPs of the thirteen different ST19 S. typhimurium strains by either Parsnp or CSI Phylogeny are listed in Supplementary Table S5. A SNP map demonstrating the SNP differences among thirteen strains using the identified SNPs from CSI Phylogeny and S. typhimurium str. RM9437 as the reference genome was generated using BLAST Ring Image Generator (BRIG) software [53]. Locations of the forward and reverse coding DNA sequences (CDSs), prophage, and repeat regions in the RM9437 genome were also included for reference to location of SNPs.

4.6. Prophage and Antibiotic Resistance Genes

Prophage sequences in all 13 of the strains used in the detailed comparative genomic analysis (Table 1) were identified using the online tool PHASTER (www.phaster.ca; accessed 11 February, 2021) [54]. To compare the identified Enterobacteria phage ST104 prophage sequences of the non-DT104 strain CDC 2010K-1587 and the six DT104 strains that contain the prophage, the sequences were aligned using the MUSCLE plugin with default parameters and 10 iterations in Geneious Prime 2020.1.2 software. Whereas, antibiotic resistance genes present in the genome of the different S. typhimurium strains were identified using the online tool ResFinder (v4.0) [55] from the Center for Genomic Epidemiology.

4.7. Allele-Based Core Genome Multilocus Sequence Typing (cgMLST) Analysis

Allele-based genomic epidemiology analysis of the three Colorado DT104 outbreak strains (RM9437, RM9436, and RM9435) were compared to other DT104 strains from around the world, whereas the other twelve S. typhimurium DT104 and non-DT104 strains from bovine and clinical sources were compared to bovine-associated ST19 strains publicly available on EnteroBase on 11 February, 2021 or 14 May, 2020, respectively. An MST utilizing the improved minimal spanning tree algorithm MSTree V2 of GrapeTree [56] based on cgMLST allele differences of 2978 alleles was generated using the three DT104 Colorado outbreak strains sequenced for this study and 556 DT104 S. typhimurium strains comprising 330 different STs in EnteroBase at a hierarchal clustering level 100 (HC100). The cgMLST scheme utilized in Enterobase that was employed in this study has been previously described [30]. MST based on cgMLST allele differences was also generated using the fifteen ST19 strains and 775 ST19 S. typhimurium strains with 682 STs isolated from different bovine-associated sources in North America.

4.8. Genes Unique to DT104 and Cattle DT104 Strains Analysis

Genome-wide association studies (pan-GWAS) analysis was conducted on the thirteen ST19 S. typhimurium genomes (Table 1). First, each genome was annotated using Prokka (v1.13.3) [57] to the bacteria kingdom using default parameters. Next, the pan genome at 90% protein identity was determined for the thirteen strains using the program Roary (v3.12.0) [58] with the following options: create a multiFASTA alignment of core genes using PRANK and fast core gene alignment with MAFFT (Roary command: roary -f./DT104 -e -n -v -i 90./DT104/*.gff). The number of genes in the core genome, soft-core genome, cloud genes, and shell genes that were output from Roary was visualized using Roary.plots.py program. The entire Roary process was repeated at 90% protein identity for each of the following groups: DT104 strains, non-DT104 strains, bovine-associated strains, and clinical strains. Identification of those genes unique to DT104 strains compared to the other ST19 non-DT104 strains were identified using the output of Roary as input into the program Scoary (v1.6.16) [59], and the analysis was run with default parameters and traits file containing the phage type for each of the strains (Scoary command: scoary -g gene_presence_absence.csv -t DT104_traits.csv). The Scoary analysis was repeated to identify genes unique to non-DT104 strains, bovine strains, and clinical strains among the thirteen ST19 S. typhimurium strains used in this study. An expanded DT104 cattle-associated unique gene analysis was conducted by selecting the 2009 Colorado outbreak ground beef strain RM9437, 25 cattle-associated DT104 strains from EnteroBase, and an additional 72 random DT104 genomes including at least one representative from each source of isolation from each available country in EnteroBase (Supplementary Table S4; excluding RM9435 and RM9436). The remaining 50 strains were random selected and downloaded from the list of 556 DT104 strains in Enterobase. The 98 strains were annotated with Prokka and analyzed via Roary at 90% protein identity as described above. To determine the cattle-specific DT104 unique genes, the output of Roary was input into Scoary as described above with the traits file containing cattle/bovine/beef versus other for each DT104 strain.

4.9. Pangenome Analysis and Visualization with Anvi’o Software

The 100 DT104 strains described above were also used to visualize the pangenome by analysis with Anvi’o software (v6.2) [60] using the pangenome workflow [61]. Pangenome analysis was conducted using NCBI’s blastp program with a minimum bit score 0.5 for gene clustering and using the MCL algorithm [62] to identify gene clusters. The gene clusters determined in the analysis were binned based on the following parameters: (1) core genes–clusters present in 100% of the genomes; (2) accessory genes–gene clusters present in ≤99% of the genomes but more than one genome; (3) singleton genes–gene clusters present in only one genome maximum. The country and source of isolation were imported as additional layers in the database.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/pathogens10050529/s1. Figure S1: Phylogenetic analysis of thirteen ST19 S. typhimurium strains associated with either clinical or bovine sources. Figure S2: Genetic variation of DT104 strains compared to non-DT104 strains of the ST19 lineage. Supplementary Table S1: List of core genes for all 13 strains of ST19 lineage from bovine and clinical sources including DT104 and non-DT104 S. typhimurium strains and DT104 strains only, non-DT104 strains only, bovine strains and clinical strains. Supplementary Table S2: List of DT104 unique genes from Roary/Scoary analysis of 13 strains of ST19 lineage from bovine and clinical sources including DT104 and non-DT104 S. typhimurium. Supplementary Table S3: List of non-DT104 unique genes from Roary/Scoary analysis of strains of ST19 lineage from bovine and clinical sources including DT104 and non-DT104 S. typhimurium. Supplementary Table S4: Metadata of strains used for DT104 cgMLST, bovine cgMLST, and DT104 pangenome analyses. Supplementary Table S5: SNP analysis between 13 strains of ST19 lineage from bovine and clinical sources including DT104 and non-DT104 S. typhimurium strains.

Author Contributions

C.T.P. was responsible for conceptualization, methodology, validation, formal analysis, writing—original draft preparation, writing—review and editing; S.H. was responsible for investigation, validation, and writing—review and editing; A.A. and A.S.O. were responsible for formal analysis, investigation, writing—review and editing; K.K.C. was responsible for conceptualization, methodology, validation, formal analysis, investigation, writing—original draft preparation, visualization, and writing–review and editing, supervision, project administration and funding. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Technology and Research Initiative Fund (TRIF) provided to Kerry Cooper by the University of Arizona, and also the United States Department of Agriculture, Agricultural Research Service CRIS project 2030-42000-051-00D.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The sequence data presented or used in this study are openly available at either NCBI Genbank under the accession numbers listed in Table 1 or Enterobase as listed in Supplementary Table S4.

Acknowledgments

We thank Robert Mandrell and Anne Bates for curation of the Produce Safety and Microbiology Research Unit strain collection. We also thank Mary Cichon at the Colorado Department of Public Health and Environment for providing strains.

Conflicts of Interest

The authors declare that there are no conflicts of interest.

References

  1. Voetsch, A.C.; Van Gilder, T.J.; Angulo, F.J.; Farley, M.M.; Shallow, S.; Marcus, R.; Cieslak, P.R.; Deneen, V.C.; Tauxe, R.V.; Emerging Infections Program FoodNet Working Group. FoodNet estimate of the burden of illness caused by nontyphoidal Salmonella infections in the United States. Clin. Infect. Dis. 2004, 38 (Suppl. S3), S127–S134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Majowicz, S.E.; Musto, J.; Scallan, E.; Angulo, F.J.; Kirk, M.; O’Brien, S.J.; Jones, T.F.; Fazil, A.; Hoekstra, R.M.; International Collaboration on Enteric Disease ‘Burden of Illness Studies. The global burden of nontyphoidal Salmonella gastroenteritis. Clin. Infect. Dis. 2010, 50, 882–889. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. CDC. National Enteric Disease Surveillance: Salmonella Annual Report; Centers for Disease Control and Prevention: Atlanta, GA, USA, 2016. Available online: http://www.cdc.gov/nationalsurveillance/PDFs/NationalSalmSurveillOverview_508.pdf (accessed on 28 February 2021).
  4. Laufer, A.S.; Grass, J.; Holt, K.; Whichard, J.M.; Griffin, P.M.; Gould, L.H. Outbreaks of Salmonella infections attributed to beef—United States, 1973–2011. Epidemiol. Infect. 2015, 143, 2003–2013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. OECD.org. Meat Consumption (Indicator). Available online: https://data.oecd.org/agroutput/meat-consumption.htm (accessed on 12 July 2019).
  6. Trujillo, S.; Keys, C.E.; Brown, E.W. Evaluation of the taxonomic utility of six-enzyme pulsed-field gel electrophoresis in reconstructing Salmonella subspecies phylogeny. Infect. Genet. Evol. 2011, 11, 92–102. [Google Scholar] [CrossRef] [PubMed]
  7. Alexander, D.C.; Fitzgerald, S.F.; DePaulo, R.; Kitzul, R.; Daku, D.; Levett, P.N.; Cameron, A.D. Laboratory-Acquired Infection with Salmonella enterica Serovar Typhimurium Exposed by Whole-Genome Sequencing. J. Clin. Microbiol. 2016, 54, 190–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Gymoese, P.; Sorensen, G.; Litrup, E.; Olsen, J.E.; Nielsen, E.M.; Torpdahl, M. Investigation of Outbreaks of Salmonella enterica Serovar Typhimurium and Its Monophasic Variants Using Whole-Genome Sequencing, Denmark. Emerg. Infect. Dis. 2017, 23, 1631–1639. [Google Scholar] [CrossRef]
  9. Octavia, S.; Wang, Q.; Tanaka, M.M.; Kaur, S.; Sintchenko, V.; Lan, R. Delineating community outbreaks of Salmonella enterica serovar Typhimurium by use of whole-genome sequencing: Insights into genomic variability within an outbreak. J. Clin. Microbiol. 2015, 53, 1063–1071. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Phillips, A.; Sotomayor, C.; Wang, Q.; Holmes, N.; Furlong, C.; Ward, K.; Howard, P.; Octavia, S.; Lan, R.; Sintchenko, V. Whole genome sequencing of Salmonella typhimurium illuminates distinct outbreaks caused by an endemic multi-locus variable number tandem repeat analysis type in Australia, 2014. BMC Microbiol. 2016, 16, 211. [Google Scholar] [CrossRef] [Green Version]
  11. Ashton, P.M.; Peters, T.; Ameh, L.; McAleer, R.; Petrie, S.; Nair, S.; Muscat, I.; de Pinna, E.; Dallman, T. Whole Genome Sequencing for the Retrospective Investigation of an Outbreak of Salmonella typhimurium DT 8. PLoS Curr. 2015, 7. [Google Scholar] [CrossRef]
  12. Helms, M.; Ethelberg, S.; Molbak, K.; Group, D.T.S. International Salmonella typhimurium DT104 infections, 1992–2001. Emerg. Infect. Dis. 2005, 11, 859–867. [Google Scholar] [CrossRef] [PubMed]
  13. Ethelberg, S.; Sorensen, G.; Kristensen, B.; Christensen, K.; Krusell, L.; Hempel-Jorgensen, A.; Perge, A.; Nielsen, E.M. Outbreak with multi-resistant Salmonella typhimurium DT104 linked to carpaccio, Denmark, 2005. Epidemiol. Infect. 2007, 135, 900–907. [Google Scholar] [CrossRef] [PubMed]
  14. Kivi, M.; Hofhuis, A.; Notermans, D.W.; Wannet, W.J.; Heck, M.E.; Van De Giessen, A.W.; Van Duynhoven, Y.T.; Stenvers, O.F.; Bosman, A.; Van Pelt, W. A beef-associated outbreak of Salmonella typhimurium DT104 in The Netherlands with implications for national and international policy. Epidemiol. Infect. 2007, 135, 890–899. [Google Scholar] [CrossRef] [PubMed]
  15. Mindlin, M.J.; Lang, N.; Maguire, H.; Walsh, B.; Verlander, N.Q.; Lane, C.; Taylor, C.; Bishop, L.A.; Crook, P.D. Outbreak investigation and case-control study: Penta-resistant Salmonella typhimurium DT104 associated with biltong in London in 2008. Epidemiol. Infect. 2013, 141, 1920–1927. [Google Scholar] [CrossRef]
  16. United States Department of Agriculture. Colorado firm recalls ground beef products due to possible Salmonella contamination. In FSIS-RC-039-2009; United States Department of Agriculture: Washington, DC, USA, 2009. [Google Scholar]
  17. Dechet, A.M.; Scallan, E.; Gensheimer, K.; Hoekstra, R.; Gunderman-King, J.; Lockett, J.; Wrigley, D.; Chege, W.; Sobel, J.; Multistate Working, G. Outbreak of multidrug-resistant Salmonella enterica serotype Typhimurium Definitive Type 104 infection linked to commercial ground beef, northeastern United States, 2003–2004. Clin. Infect. Dis. 2006, 42, 747–752. [Google Scholar] [CrossRef]
  18. Threlfall, E.J. Epidemic Salmonella typhimurium DT 104—A truly international multiresistant clone. J. Antimicrob. Chemother. 2000, 46, 7–10. [Google Scholar] [CrossRef]
  19. Boyd, D.A.; Peters, G.A.; Ng, L.; Mulvey, M.R. Partial characterization of a genomic island associated with the multidrug resistance region of Salmonella enterica Typhymurium DT104. FEMS Microbiol. Lett. 2000, 189, 285–291. [Google Scholar] [CrossRef] [PubMed]
  20. Ribot, E.M.; Wierzba, R.K.; Angulo, F.J.; Barrett, T.J. Salmonella enterica serotype Typhimurium DT104 isolated from humans, United States, 1985, 1990, and 1995. Emerg. Infect. Dis. 2002, 8, 387–391. [Google Scholar] [CrossRef] [Green Version]
  21. Prendergast, D.M.; O’Grady, D.; Fanning, S.; Cormican, M.; Delappe, N.; Egan, J.; Mannion, C.; Fanning, J.; Gutierrez, M. Application of multiple locus variable number of tandem repeat analysis (MLVA), phage typing and antimicrobial susceptibility testing to subtype Salmonella enterica serovar Typhimurium isolated from pig farms, pork slaughterhouses and meat producing plants in Ireland. Food Microbiol. 2011, 28, 1087–1094. [Google Scholar] [CrossRef] [PubMed]
  22. Mather, A.E.; Reid, S.W.; Maskell, D.J.; Parkhill, J.; Fookes, M.C.; Harris, S.R.; Brown, D.J.; Coia, J.E.; Mulvey, M.R.; Gilmour, M.W.; et al. Distinguishable epidemics of multidrug-resistant Salmonella typhimurium DT104 in different hosts. Science 2013, 341, 1514–1517. [Google Scholar] [CrossRef] [Green Version]
  23. Leekitcharoenphon, P.; Hendriksen, R.S.; Le Hello, S.; Weill, F.X.; Baggesen, D.L.; Jun, S.R.; Ussery, D.W.; Lund, O.; Crook, D.W.; Wilson, D.J.; et al. Global Genomic Epidemiology of Salmonella enterica serovar Typhimurium DT104. Appl. Environ. Microbiol. 2016, 82, 2516–2526. [Google Scholar] [CrossRef] [Green Version]
  24. Briggs, C.E.; Fratamico, P.M. Molecular characterization of an antibiotic resistance gene cluster of Salmonella typhimurium DT104. Antimicrob. Agents Chemother. 1999, 43, 846–849. [Google Scholar] [CrossRef] [Green Version]
  25. Boyd, D.; Peters, G.A.; Cloeckaert, A.; Boumedine, K.S.; Chaslus-Dancla, E.; Imberechts, H.; Mulvey, M.R. Complete nucleotide sequence of a 43-kilobase genomic island associated with the multidrug resistance region of Salmonella enterica serovar Typhimurium DT104 and its identification in phage type DT120 and serovar Agona. J. Bacteriol. 2001, 183, 5725–5732. [Google Scholar] [CrossRef] [Green Version]
  26. Nguyen, S.V.; Harhay, D.M.; Bono, J.L.; Smith, T.P.; Fields, P.I.; Dinsmore, B.A.; Santovenia, M.; Kelley, C.M.; Wang, R.; Bosilevac, J.M.; et al. Complete, Closed Genome Sequences of 10 Salmonella enterica subsp. enterica Serovar Typhimurium Strains Isolated from Human and Bovine Sources. Genome Announc. 2016, 4. [Google Scholar] [CrossRef] [Green Version]
  27. Zhao, E.Y.; Bao, H.X.; Tang, L.; Zou, Q.H.; Liu, W.Q.; Zhu, D.L.; Chin, J.; Dong, Y.Y.; Li, Y.G.; Cao, F.L.; et al. Genomic comparison of Salmonella typhimurium DT104 with non-DT104 strains. Mol. Genet. Genom. 2013, 288, 549–557. [Google Scholar] [CrossRef]
  28. Tanaka, K.; Nishimori, K.; Makino, S.; Nishimori, T.; Kanno, T.; Ishihara, R.; Sameshima, T.; Akiba, M.; Nakazawa, M.; Yokomizo, Y.; et al. Molecular characterization of a prophage of Salmonella enterica serotype Typhimurium DT104. J. Clin. Microbiol. 2004, 42, 1807–1812. [Google Scholar] [CrossRef] [Green Version]
  29. Jia, B.; Raphenya, A.R.; Alcock, B.; Waglechner, N.; Guo, P.; Tsang, K.K.; Lago, B.A.; Dave, B.M.; Pereira, S.; Sharma, A.N.; et al. CARD 2017: Expansion and model-centric curation of the comprehensive antibiotic resistance database. Nucleic Acids Res. 2017, 45, D566–D573. [Google Scholar] [CrossRef]
  30. Zhou, Z.; Alikhan, N.F.; Mohamed, K.; Fan, Y.; Agama Study, G.; Achtman, M. The EnteroBase user’s guide, with case studies on Salmonella transmissions, Yersinia pestis phylogeny, and Escherichia core genomic diversity. Genome Res. 2020, 30, 138–152. [Google Scholar] [CrossRef] [Green Version]
  31. Lawson, A.J.; Desai, M.; O’Brien, S.J.; Davies, R.H.; Ward, L.R.; Threlfall, E.J. Molecular characterisation of an outbreak strain of multiresistant Salmonella enterica serovar Typhimurium DT104 in the UK. Clin. Microbiol. Infect. 2004, 10, 143–147. [Google Scholar] [CrossRef]
  32. Jackson, B.R.; Tarr, C.; Strain, E.; Jackson, K.A.; Conrad, A.; Carleton, H.; Katz, L.S.; Stroika, S.; Gould, L.H.; Mody, R.K.; et al. Implementation of Nationwide Real-time Whole-genome Sequencing to Enhance Listeriosis Outbreak Detection and Investigation. Clin. Infect. Dis. 2016, 63, 380–386. [Google Scholar] [CrossRef] [Green Version]
  33. Kingsley, R.A.; Kay, S.; Connor, T.; Barquist, L.; Sait, L.; Holt, K.E.; Sivaraman, K.; Wileman, T.; Goulding, D.; Clare, S.; et al. Genome and transcriptome adaptation accompanying emergence of the definitive type 2 host-restricted Salmonella enterica serovar Typhimurium pathovar. MBio 2013, 4, e00565-13. [Google Scholar] [CrossRef] [Green Version]
  34. Carattoli, A.; Filetici, E.; Villa, L.; Dionisi, A.M.; Ricci, A.; Luzzi, I. Antibiotic resistance genes and Salmonella genomic island 1 in Salmonella enterica serovar Typhimurium isolated in Italy. Antimicrob. Agents Chemother. 2002, 46, 2821–2828. [Google Scholar] [CrossRef] [Green Version]
  35. Baggesen, D.L.; Sandvang, D.; Aarestrup, F.M. Characterization of Salmonella enterica serovar typhimurium DT104 isolated from Denmark and comparison with isolates from Europe and the United States. J. Clin. Microbiol. 2000, 38, 1581–1586. [Google Scholar] [CrossRef] [Green Version]
  36. Ng, L.K.; Mulvey, M.R.; Martin, I.; Peters, G.A.; Johnson, W. Genetic characterization of antimicrobial resistance in Canadian isolates of Salmonella serovar Typhimurium DT104. Antimicrob. Agents Chemother. 1999, 43, 3018–3021. [Google Scholar] [CrossRef] [Green Version]
  37. Matiasovicova, J.; Adams, P.; Barrow, P.A.; Hradecka, H.; Malcova, M.; Karpiskova, R.; Budinska, E.; Pilousova, L.; Rychlik, I. Identification of putative ancestors of the multidrug-resistant Salmonella enterica serovar typhimurium DT104 clone harboring the Salmonella genomic island 1. Arch. Microbiol. 2007, 187, 415–424. [Google Scholar] [CrossRef]
  38. Matiasovicova, J.; Havlickova, H.; Sisak, F.; Pilousova, L.; Rychlik, I. allB, allantoin utilisation and Salmonella enterica serovar Enteritidis and Typhimurium colonisation of poultry and mice. Folia Microbiol. 2011, 56, 264–269. [Google Scholar] [CrossRef]
  39. Graziani, C.; Busani, L.; Dionisi, A.M.; Lucarelli, C.; Owczarek, S.; Ricci, A.; Mancin, M.; Caprioli, A.; Luzzi, I. Antimicrobial resistance in Salmonella enterica serovar Typhimurium from human and animal sources in Italy. Vet. Microbiol. 2008, 128, 414–418. [Google Scholar] [CrossRef] [PubMed]
  40. Alikhan, N.F.; Zhou, Z.; Sergeant, M.J.; Achtman, M. A genomic overview of the population structure of Salmonella. PLoS Genet. 2018, 14, e1007261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Parker, C.T.; Huynh, S.; Gorski, L.; Cooper, K.K.; Miller, W.G. Complete Genome Sequences of Two Outbreak Strains of Salmonella enterica subsp. enterica Serovar Thompson Associated with Cilantro. Genome Announc. 2015, 3. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Darling, A.E.; Mau, B.; Perna, N.T. progressiveMauve: Multiple genome alignment with gene gain, loss and rearrangement. PLoS ONE 2010, 5, e11147. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Merga, J.Y.; Winstanley, C.; Williams, N.J.; Yee, E.; Miller, W.G. Complete Genome Sequence of the Arcobacter butzleri Cattle Isolate 7h1h. Genome Announc. 2013, 1. [Google Scholar] [CrossRef] [Green Version]
  44. Kearse, M.; Moir, R.; Wilson, A.; Stones-Havas, S.; Cheung, M.; Sturrock, S.; Buxton, S.; Cooper, A.; Markowitz, S.; Duran, C.; et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics 2012, 28, 1647–1649. [Google Scholar] [CrossRef] [PubMed]
  45. Delcher, A.L.; Bratke, K.A.; Powers, E.C.; Salzberg, S.L. Identifying bacterial genes and endosymbiont DNA with Glimmer. Bioinformatics 2007, 23, 673–679. [Google Scholar] [CrossRef]
  46. Tatusova, T.; DiCuccio, M.; Badretdin, A.; Chetvernin, V.; Nawrocki, E.P.; Zaslavsky, L.; Lomsadze, A.; Pruitt, K.D.; Borodovsky, M.; Ostell, J. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res. 2016, 44, 6614–6624. [Google Scholar] [CrossRef]
  47. Kaas, R.S.; Leekitcharoenphon, P.; Aarestrup, F.M.; Lund, O. Solving the problem of comparing whole bacterial genomes across different sequencing platforms. PLoS ONE 2014, 9, e104984. [Google Scholar] [CrossRef] [Green Version]
  48. Treangen, T.J.; Ondov, B.D.; Koren, S.; Phillippy, A.M. The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecific microbial genomes. Genome Biol. 2014, 15, 524. [Google Scholar] [CrossRef] [Green Version]
  49. Edgar, R.C. MUSCLE: A multiple sequence alignment method with reduced time and space complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [Green Version]
  50. Edgar, R.C. MUSCLE: Multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [Green Version]
  51. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular Evolutionary Genetics Analysis across Computing Platforms. Mol. Biol. Evol. 2018, 35, 1547–1549. [Google Scholar] [CrossRef] [PubMed]
  52. Tamura, K.; Nei, M. Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol. Biol. Evol. 1993, 10, 512–526. [Google Scholar] [CrossRef]
  53. Alikhan, N.F.; Petty, N.K.; Ben Zakour, N.L.; Beatson, S.A. BLAST Ring Image Generator (BRIG): Simple prokaryote genome comparisons. BMC Genom. 2011, 12, 402. [Google Scholar] [CrossRef] [Green Version]
  54. Arndt, D.; Grant, J.R.; Marcu, A.; Sajed, T.; Pon, A.; Liang, Y.; Wishart, D.S. PHASTER: A better, faster version of the PHAST phage search tool. Nucleic Acids Res. 2016, 44, W16–W21. [Google Scholar] [CrossRef] [Green Version]
  55. Zankari, E.; Hasman, H.; Cosentino, S.; Vestergaard, M.; Rasmussen, S.; Lund, O.; Aarestrup, F.M.; Larsen, M.V. Identification of acquired antimicrobial resistance genes. J. Antimicrob. Chemother. 2012, 67, 2640–2644. [Google Scholar] [CrossRef]
  56. Zhou, Z.; Alikhan, N.F.; Sergeant, M.J.; Luhmann, N.; Vaz, C.; Francisco, A.P.; Carrico, J.A.; Achtman, M. GrapeTree: Visualization of core genomic relationships among 100,000 bacterial pathogens. Genome Res. 2018, 28, 1395–1404. [Google Scholar] [CrossRef] [Green Version]
  57. Seemann, T. Prokka: Rapid prokaryotic genome annotation. Bioinformatics 2014, 30, 2068–2069. [Google Scholar] [CrossRef]
  58. Page, A.J.; Cummins, C.A.; Hunt, M.; Wong, V.K.; Reuter, S.; Holden, M.T.; Fookes, M.; Falush, D.; Keane, J.A.; Parkhill, J. Roary: Rapid large-scale prokaryote pan genome analysis. Bioinformatics 2015, 31, 3691–3693. [Google Scholar] [CrossRef]
  59. Brynildsrud, O.; Bohlin, J.; Scheffer, L.; Eldholm, V. Rapid scoring of genes in microbial pan-genome-wide association studies with Scoary. Genome Biol. 2016, 17, 238. [Google Scholar] [CrossRef] [Green Version]
  60. Eren, A.M.; Esen, O.C.; Quince, C.; Vineis, J.H.; Morrison, H.G.; Sogin, M.L.; Delmont, T.O. Anvi’o: An advanced analysis and visualization platform for ‘omics data. PeerJ 2015, 3, e1319. [Google Scholar] [CrossRef] [PubMed]
  61. Delmont, T.O.; Eren, A.M. Linking pangenomes and metagenomes: The Prochlorococcus metapangenome. PeerJ 2018, 6, e4320. [Google Scholar] [CrossRef] [Green Version]
  62. Van Dongen, S.; Abreu-Goodger, C. Using MCL to extract clusters from networks. Methods Mol. Biol. 2012, 804, 281–295. [Google Scholar] [CrossRef]
Figure 1. SNP map between ST19 S. typhimurium strains. Location of SNPs from different bovine and clinical S. typhimurium strains compared to RM9437. Beginning in the inner circle and moving outwards: (1) RM9437 repeat regions; (2) RM9437 prophage regions; (3) RM9437 forward CDSs; (4) RM9437 reverse CDSs; bovine strains: (5) SA972816 (DT104); (6) USMARC-1808 (DT104); (7) USMARC-1810; (8) USMARC-1880; (9) USMARC-1896; (10) USMARC-1898; human clinical strains: (11) DT104 (DT104); (12) CDC H2662 (DT104); (13) CDC 2009K-1640 (DT104); (14) CDC 2011K-1702; (15) CDC 2009K-2059; (16) CDC 2010K-1587.
Figure 1. SNP map between ST19 S. typhimurium strains. Location of SNPs from different bovine and clinical S. typhimurium strains compared to RM9437. Beginning in the inner circle and moving outwards: (1) RM9437 repeat regions; (2) RM9437 prophage regions; (3) RM9437 forward CDSs; (4) RM9437 reverse CDSs; bovine strains: (5) SA972816 (DT104); (6) USMARC-1808 (DT104); (7) USMARC-1810; (8) USMARC-1880; (9) USMARC-1896; (10) USMARC-1898; human clinical strains: (11) DT104 (DT104); (12) CDC H2662 (DT104); (13) CDC 2009K-1640 (DT104); (14) CDC 2011K-1702; (15) CDC 2009K-2059; (16) CDC 2010K-1587.
Pathogens 10 00529 g001
Figure 2. Phylogenetic relationship of S. typhimurium strains. Maximum-likelihood (ML) tree of all the bovine and clinical strains used in the study based on 1798 SNPs identified between the 13 strains by Parsnp software using S. typhimurium str. RM9437 as the reference strain. ML tree generated with MEGAX software using Tamura–Nei model with a discrete Gamma distribution with 1000 pseudoreplicates. Antibiotic resistance gene and prophage profiles for each of the strains were identified with ResFinder and PHASTER, respectively (blue—present; yellow—absent). Note: * Includes only bacteriophage identified as intact by PHASTER.
Figure 2. Phylogenetic relationship of S. typhimurium strains. Maximum-likelihood (ML) tree of all the bovine and clinical strains used in the study based on 1798 SNPs identified between the 13 strains by Parsnp software using S. typhimurium str. RM9437 as the reference strain. ML tree generated with MEGAX software using Tamura–Nei model with a discrete Gamma distribution with 1000 pseudoreplicates. Antibiotic resistance gene and prophage profiles for each of the strains were identified with ResFinder and PHASTER, respectively (blue—present; yellow—absent). Note: * Includes only bacteriophage identified as intact by PHASTER.
Pathogens 10 00529 g002
Figure 3. cgMLST analysis of S. typhimurium DT104 strains. Core genome multilocus sequencing type (cgMLST) minimum spanning tree (MST) based on 2978 alleles of 556 DT104 strains representing 330 sequence types (STs) present in Enterobase on 11 February 2021. The numbers represent number of alleles between strains, and the black arrows mark the three United States strains sequenced in this study. Nodes are colored according to source type, while grey nodes (ND/others) and white nodes (missing) had no source type metadata in Enterobase. Logarithmic scale represents number of alleles.
Figure 3. cgMLST analysis of S. typhimurium DT104 strains. Core genome multilocus sequencing type (cgMLST) minimum spanning tree (MST) based on 2978 alleles of 556 DT104 strains representing 330 sequence types (STs) present in Enterobase on 11 February 2021. The numbers represent number of alleles between strains, and the black arrows mark the three United States strains sequenced in this study. Nodes are colored according to source type, while grey nodes (ND/others) and white nodes (missing) had no source type metadata in Enterobase. Logarithmic scale represents number of alleles.
Pathogens 10 00529 g003
Figure 4. Pangenome of 100 DT104 strains. Visualization of pangenome using Anvi’o software for 100 DT104 strains from Enterobase including the 3 Colorado outbreak strains (RM9435, RM9436, and RM9437), 25 bovine-associated strains, and 72 random selected strains from other sources. Partial rings represent comparison of DT104 genomes, including core genes present in all 100 strains, accessory genes present in ≥2 strains but ≤99 strains, and singleton genes that are only present in a single DT104 genome. Color of the ring is based on the source of isolation/source type. Order of the rings from outer most to inner most: (1) number of DT104 genomes contributing to gene cluster marked in strain ring; (2) number of genes in the gene cluster; (3) maximum number of paralogs in gene cluster; (3) geometric homogeneity identity of gene cluster; (4) functional homogeneity identity of gene cluster; (5) combined homogeneity identity of gene cluster; S. typhimurium strains: (6) R13; (7) BD-DU-2006-02-006862; (8) AH-00-5945; (9) BD-DU-2006-01-008520; (10) JM-08-31; (11) P5368550; (12) H041040309; (13) JM-05-18; (14) JM-04-98; (15) JM-04-96; (16) JM-04-88; (17) JM-04-66; (18) JM-04-33; (19) JM-04-32; (20) JM-04-26; (21) JM-04-16; (22) SE04-41; (23) DP-J28; (24) DP-J21; (25) DP-J20; (26) JE-S09-001058; (27) S-20031116-1; (28) AH-02-7174; (29) P5372472; (30) CK-95-00014475; (31) R11; (32) DP-A73; (33) AH-03-7513; (34) NL-A3; (35) JE-S09-002301; (36) DP-A58; (37) BD-OM-2007-04-003477; (38) MC-05-0834; (39) JM-01-36; (40) R9; (41) JE-4963; (42) JM-06-45; (43) JE-S09-002300; (44) R4; (45) A1519-04; (46) JM-02-29; (47) DP-A68; (48) JE-4657; (49) R5; (50) BD-OM-2007-07-005180; (51) MC-04-0302; (52) JM-04-300; (53) BD-DU-2006-07-000043; (54) R12; (55) MC-04-0395; (56) BD-DU-2006-06-004373; (57) R7; (58) JE-S09-002567; (59) BD-OM-2007-07-005216; (60) R10; (61) JE-S09-001062; (62) JM-08-01; (63) MC-05-0863; (64) JE-S09-002695; (65) BD-OM-2006-10-009515; (66) DP-A60; (67) DP-A43; (68) DP-A59; (69) JE-4960; (70) JM-01-80; (71) RM9435; (72) RM09437; (73) RM9436; (74) JM-03-73; (75) JE-3481; (76) SE04-155; (77) MC-04-0146; (78) SE04-154; (79) JE-S09-000879; (80) R8; (81) MC-05-0747; (82) BD-DU-2006-10-005166; (83) JM-07-08; (84) SM-DU-2007-09-007135; (85) SE06-42; (86) SM-DU-2007-03-009861; (87) AH-06-8974; (88) NL-I12; (89) AH-03-7508; (90) SE06-34; (91) SE06-33; (92) BD-PM-2007-09-006137; (93) NL-I11; (94) AH-03-7454; (95) R3; (96) NL-R64; (97) JE-S09-001059; (98) DP-D21; (99) JE-843; (100) JE-1363; (101) JE-S09-002574; (102) JE-S09-002079; (103) JE-S09-002077; (104) JE-S09-002078; (105) S-20030672-1.
Figure 4. Pangenome of 100 DT104 strains. Visualization of pangenome using Anvi’o software for 100 DT104 strains from Enterobase including the 3 Colorado outbreak strains (RM9435, RM9436, and RM9437), 25 bovine-associated strains, and 72 random selected strains from other sources. Partial rings represent comparison of DT104 genomes, including core genes present in all 100 strains, accessory genes present in ≥2 strains but ≤99 strains, and singleton genes that are only present in a single DT104 genome. Color of the ring is based on the source of isolation/source type. Order of the rings from outer most to inner most: (1) number of DT104 genomes contributing to gene cluster marked in strain ring; (2) number of genes in the gene cluster; (3) maximum number of paralogs in gene cluster; (3) geometric homogeneity identity of gene cluster; (4) functional homogeneity identity of gene cluster; (5) combined homogeneity identity of gene cluster; S. typhimurium strains: (6) R13; (7) BD-DU-2006-02-006862; (8) AH-00-5945; (9) BD-DU-2006-01-008520; (10) JM-08-31; (11) P5368550; (12) H041040309; (13) JM-05-18; (14) JM-04-98; (15) JM-04-96; (16) JM-04-88; (17) JM-04-66; (18) JM-04-33; (19) JM-04-32; (20) JM-04-26; (21) JM-04-16; (22) SE04-41; (23) DP-J28; (24) DP-J21; (25) DP-J20; (26) JE-S09-001058; (27) S-20031116-1; (28) AH-02-7174; (29) P5372472; (30) CK-95-00014475; (31) R11; (32) DP-A73; (33) AH-03-7513; (34) NL-A3; (35) JE-S09-002301; (36) DP-A58; (37) BD-OM-2007-04-003477; (38) MC-05-0834; (39) JM-01-36; (40) R9; (41) JE-4963; (42) JM-06-45; (43) JE-S09-002300; (44) R4; (45) A1519-04; (46) JM-02-29; (47) DP-A68; (48) JE-4657; (49) R5; (50) BD-OM-2007-07-005180; (51) MC-04-0302; (52) JM-04-300; (53) BD-DU-2006-07-000043; (54) R12; (55) MC-04-0395; (56) BD-DU-2006-06-004373; (57) R7; (58) JE-S09-002567; (59) BD-OM-2007-07-005216; (60) R10; (61) JE-S09-001062; (62) JM-08-01; (63) MC-05-0863; (64) JE-S09-002695; (65) BD-OM-2006-10-009515; (66) DP-A60; (67) DP-A43; (68) DP-A59; (69) JE-4960; (70) JM-01-80; (71) RM9435; (72) RM09437; (73) RM9436; (74) JM-03-73; (75) JE-3481; (76) SE04-155; (77) MC-04-0146; (78) SE04-154; (79) JE-S09-000879; (80) R8; (81) MC-05-0747; (82) BD-DU-2006-10-005166; (83) JM-07-08; (84) SM-DU-2007-09-007135; (85) SE06-42; (86) SM-DU-2007-03-009861; (87) AH-06-8974; (88) NL-I12; (89) AH-03-7508; (90) SE06-34; (91) SE06-33; (92) BD-PM-2007-09-006137; (93) NL-I11; (94) AH-03-7454; (95) R3; (96) NL-R64; (97) JE-S09-001059; (98) DP-D21; (99) JE-843; (100) JE-1363; (101) JE-S09-002574; (102) JE-S09-002079; (103) JE-S09-002077; (104) JE-S09-002078; (105) S-20030672-1.
Pathogens 10 00529 g004
Figure 5. cgMLST analysis of North American S. typhimurium bovine strains from ST19 lineage. Core genome multilocus sequence type (cgMLST) minimum spanning tree (MST) based on 2978 alleles of 775 S. typhimurium strains isolated from different bovine sources in North America that belong to the ST19 lineage and represent 682 core genome sequence types (ST) present in Enterobase on 14 May 2020. The numbers represent number of alleles between strains, and the fifteen strains used in this study are highlighted in yellow. Nodes are colored according to source details, and white nodes had the metadata missing. The arrows represent the two non-DT104 clinical strains, whereas the star marks the node with all seven DT104 strains. Logarithmic scale represents number of alleles.
Figure 5. cgMLST analysis of North American S. typhimurium bovine strains from ST19 lineage. Core genome multilocus sequence type (cgMLST) minimum spanning tree (MST) based on 2978 alleles of 775 S. typhimurium strains isolated from different bovine sources in North America that belong to the ST19 lineage and represent 682 core genome sequence types (ST) present in Enterobase on 14 May 2020. The numbers represent number of alleles between strains, and the fifteen strains used in this study are highlighted in yellow. Nodes are colored according to source details, and white nodes had the metadata missing. The arrows represent the two non-DT104 clinical strains, whereas the star marks the node with all seven DT104 strains. Logarithmic scale represents number of alleles.
Pathogens 10 00529 g005
Table 1. Salmonella typhimurium strains used in this study.
Table 1. Salmonella typhimurium strains used in this study.
StrainChromosome Size (bp)Number of PlasmidsProphage *DT104Collection YearCountry of IsolationIsolation SourceAccession NumberReference
RM94374,936,49918 (6)Yes2009USA: ColoradoGround BeefCP012985This study
RM94364,936,49918 (6)Yes2009USA: ColoradoHuman StoolSRR8660931This study
RM94354,936,49918 (6)Yes2009USA: ColoradoHuman StoolSRR8660930This study
CDC H26624,891,16516 (4)Yes1997USAHuman StoolCP014979[26]
CDC 2009K-16404,933,70718 (5)Yes2009USA: BIFSCo Region 2Human StoolCP014975[26]
SA9728164,891,92328 (5)Yes2002ChinaBovine StoolCP007484[27]
CDC 2011K-17024,906,321110 (5)Yes2011USA: BIFSCo Region 7Human UrineCP014967[26]
USMARC-18084,936,89417 (5)Yes2005USA: BIFSCo Region 8Bovine Post-eviscerationCP014969[26]
DT1044,933,63119 (5)Yes1988EnglandHuman stoolHF937208.1[22]
CDC 2009K-20594,823,79306 (3)No2009USA: BIFSCo Region 2Human StoolCP014983[26]
CDC 2010K-15874,799,39846 (2)No2010USA: BIFSCo Region 8Human StoolCP014965[26]
USMARC-18104,927,14509 (6)No2005USA: BIFSCo Region 5Bovine Pre-eviscerationCP014982[26]
USMARC-18804,815,20806 (3)No2003USA: BIFSCo Region 5Bovine Pre-eviscerationCP014981[26]
USMARC-18964,856,44017 (4)No2011USA: BIFSCo Region 2Bovine fat trimCP014977[26]
USMARC-18984,784,38535 (2)No2007USA: BIFSCo Region 3Ground BeefCP014971[26]
* Number in parentheses represent number of intact prophage sequences.
Table 2. Total number of SNPs between different strains of S. typhimurium from human and bovine sources.
Table 2. Total number of SNPs between different strains of S. typhimurium from human and bovine sources.
RM9437CDC 2009K-1640CDC 2011K-1702CDC H2662SA972816USMARC-1808DT104USMARC-1810USMARC-1880USMARC-1896USMARC-1898CDC 2009K-2059CDC 2010K-1587
RM9437 548560306761694593627584617650
CDC 2009K-164054 7148365347682579615570605638
CDC 2011K-17028571 51676872709606642599632665
CDC H2662604851 444549686583619576609642
SA97281630366744 4943678575611566601634
USMARC-18086753684549 54691588624581614647
DT104614772494354 681578614571604637
USMARC-1810694682709686678691681 591631584621646
USMARC-1880593579606583575588578591 512367502529
USMARC-1896627615642619611624614631512 50510305
USMARC-1898584570599576566581571584367505 495524
CDC 2009K-205961760563260960161460462150210495 295
CDC 2010K-1587650638665642634647637646529305524295
Table 3. Cattle-associated S. typhimurium DT104 genes *.
Table 3. Cattle-associated S. typhimurium DT104 genes *.
Prokka Annotation% of Bovine-Associated DT104 Genomes with Gene% of Non-Bovine-Associated DT104 Genomes # with Genep-Value
Hypothetical protein0%33.3%0.0003
Beta-lactamase PSE-1 (pse1)100%75.0%0.0026
Regulatory protein Rop3.8%33.3%0.0031
Hypothetical protein3.8%33.3%0.0031
Hypothetical protein3.8%31.9%0.0032
Hypothetical protein3.8%31.9%0.0032
tRNA-Ala(tgc)19.2%1.4%0.0047
Hypothetical protein3.8%30.5%0.0059
Mobilization protein MbeC3.8%30.5%0.0059
Endoribonuclease ToxN3.8%30.5%0.0059
Dihydropteroate synthase3.8%29.2%0.0062
* Genes listed are those statistically significant (p-value < 0.01). # Non-bovine sources include avian, clinical, ovine, equine, reptile, swine, water, environment, food, and primate. Note: 2009 Colorado outbreak strains RM9435 and RM9436 were not used in the analysis as these clinical strains are identical to the ground beef strain RM9437 and from the same outbreak.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Parker, C.T.; Huynh, S.; Alexander, A.; Oliver, A.S.; Cooper, K.K. Genomic Characterization of Salmonella typhimurium DT104 Strains Associated with Cattle and Beef Products. Pathogens 2021, 10, 529. https://doi.org/10.3390/pathogens10050529

AMA Style

Parker CT, Huynh S, Alexander A, Oliver AS, Cooper KK. Genomic Characterization of Salmonella typhimurium DT104 Strains Associated with Cattle and Beef Products. Pathogens. 2021; 10(5):529. https://doi.org/10.3390/pathogens10050529

Chicago/Turabian Style

Parker, Craig T., Steven Huynh, Aaron Alexander, Andrew S. Oliver, and Kerry K. Cooper. 2021. "Genomic Characterization of Salmonella typhimurium DT104 Strains Associated with Cattle and Beef Products" Pathogens 10, no. 5: 529. https://doi.org/10.3390/pathogens10050529

APA Style

Parker, C. T., Huynh, S., Alexander, A., Oliver, A. S., & Cooper, K. K. (2021). Genomic Characterization of Salmonella typhimurium DT104 Strains Associated with Cattle and Beef Products. Pathogens, 10(5), 529. https://doi.org/10.3390/pathogens10050529

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop