Global Genomic Variations of Pathogenic Bacteria Salmonella Reading

US Centers for Disease Control and Prevention (CDC) estimates food-borne pathogenic Salmonella bacteria cause about 1.35 million infections, 26,500 hospitalizations, and 420 deaths in the U.S. every year. Outbreaks of multidrug-resistant Salmonella Reading linked mainly to raw turkey products and alfalfa sprouts have accounted for human illness including mortality. Therefore, we downloaded and compared genome sequences of 897 Salmonella Reading isolated from 27 resources in 11 countries at 21 year-points from the National Center for Biotechnology Information database. Isolates were clustered into 10 clades which consisted of 3 major and 7 minor or single node clades. Although 3 major clades included both Europe and American isolates, one clade consisted of 605 North American isolates out of 614 isolates. Evolutionary distance is more related to the continent than the source of isolation. The host source, continent (North America and Europe) and phylogenetic clade were related to the prevalence of isolates encoding Antimicrobial Resistance Genes (ARGs). Prevalence of prophages was greater in bovine and swine isolates than poultry and human isolates and the least prevalence was found in human isolates. Between continents, the prevalence of phage was greater in North American isolates than European. The diversity of virulence factors in swine isolates differed from poultry isolate while no difference was found among continents. In conclusion, evolutionary distance is related to isolation host source rather than the continent, and genome features were distinguished by host and cluster. Our genomic analysis implies that Salmonella Reading evolved indeHow to cite this paper: Lee, W.H., Jeon, Y., Chung, J.S., Shim, W.J., Shin, M.J., Han, S.G., Baek, S., Kim, J., Kim, J. and Chung, P.S. (2021) Global Genomic Variations of Pathogenic Bacteria Salmonella Reading. Computational Molecular Bioscience, 11, 57-68. https://doi.org/10.4236/cmb.2021.113004 Received: September 1, 2021 Accepted: September 27, 2021 Published: September 30, 2021 Copyright © 2021 by author(s) and Scientific Research Publishing Inc. This work is licensed under the Creative Commons Attribution International License (CC BY 4.0). http://creativecommons.org/licenses/by/4.0/ Open Access


Introduction
Salmonella, an important causative microorganism of food-borne diseases, has extensive attention for worldwide high incidence [1]. US Centers for Disease Control and Prevention (CDC) estimates food-borne pathogenic Salmonella bacteria cause about 1.35 million infections, 26,500 hospitalizations, and 420 deaths in the U.S. every year. Non-typhoidal Salmonella spp. is a leading cause of human food poisoning worldwide and is responsible for 93 million infections annually [2]. Outbreaks of Multidrug-Resistant (MDR) Salmonella Reading linked mainly to raw turkey products and Alfalfa sprouts have accounted for human illness including mortality. Since S. Reading was identified from a water supply in Reading, England in 1916, it has been reported from various sources including domestic animals, poultry, produce and environmental samples. Recently, two large outbreaks of S. Reading were reported. CDC [3] reported an outbreak from November 2017 through March 2019 and it was linked to living turkeys and raw turkey products, but no single source product or company was attri- MDR S. Newport, S. Typhimurium, and other Salmonella serovars and phage types have been associated with increased morbidity and mortality in humans and animals [4]. The emergence of bacteria ARGs occurs during the evolutionary process and vertical gene transfer [5]. The acquired resistance includes genetic exchanges through horizontal gene transfer with the acquisition of ARGs on mobile genetic elements [6]. The therapeutic and non-therapeutic uses of antibiotics in humans and animals have increased selective pressure [7]. Furthermore, pathogenic bacteria including Salmonella have been selected for ARGs due to the constant exposure to antibiotics [8].
Considering the recent multiple outbreaks in North America and its multidrug resistance, S. Reading might need to be better understood its ecology and evolution. This study aims to investigate the phylogenomic distance of and ARGs and genomic features in publicly available widespread S. Reading genomes. To achieve the purpose of this study, a comprehensive genomic investigation was performed to distinguish genomic changes related to host source and isolation locations.  (3), bovine (44), canine (1), companion animal (3), environment (12), feed (11), food (8), human (91), laboratory (1), ovine (2), poultry (510), swine (87), wild animal (1) and no data (123).

Serotype Prediction
In Silico serotype prediction was performed with the Salmonella In Silico Typing Resource [9]. Only isolates with a predicted serotype of Reading were included in downstream analyses.

Analysis of Genomes
1) Core genome (cg) SNPs analysis: Parsnp was used to construct a non-redundant core genome single nucleotide polymorphisms (SNPs) of all 897 isolates and S. enterica serovar Typhi strain CT18 was used as reference for genome alignment.
2) Phylogeny reconstruction: Maximum Likelihood phylogenic tree was generated using MEGA [10]. Number of bootstrap replications, substitution model, rates of sites, gaps/missing data treatment, ML Heuristic method were 100, Tamura-Nei model, Gamma distributed with Invariant sites (G + I) with 5 of Discrete Gamma Categories, partial deletion with 95% site coverage cutoff and Nearest-Neighbor-Interchange (NMI) with NJ/BioNJ initial tree, respectively.

Statistical Analysis
Prevalence of antimicrobial resistance genes (ARGs) and phages by host source and continent source were compared using Wilcoxon/Kruskal-Wallis Rank Sums Tests. Relative abundances of VFs by sources were compared using Non-metric multi-dimensional scaling (NMDS) on Bray-Curtis distances with vegan package in R (version 3.5.2).

Results and Discussion
A total of 121,449 sites of SNPs was found for core genomes of all 897 isolates.  Clade 10 includes 38 bovine, 11 poultry and 85 swine, and few of companion animal, environment, feed, and 8 ND. Poultry isolates dominate the Clade 1 and also identified in Clade 10 but not in Clade 3. Human isolates dominated in Clade 3 and also identified in Clades 1 and 10. Swine isolates dominated in Clade 10 and identified in Clade 1 but not in Clade 3. Two large outbreaks of S. Reading with raw turkey products in North America (23, 24) has been reported and in this study majority of poultry strains (497 isolates) are phylogenetically linked closely in Clade 1 but also two different clusters, Clades 6 and 10, were identified for North America poultry strains. Poultry stains clustered with swine stains in Clade 10 and are polyphyletic from stains in Clade 1 and at least two lineages existed for poultry stains. Among 3 major clusters, any distinct lineage was identified with continents, Europe and North America, and it might reflect the rapid dissemination of this pathogen worldwide through the increased human travels and the massive production and the broad distribution of food and feed resources.
Global emergence of pathogens in the same lineage may alert the spread of antimicrobial resistance genes and create new control challenges in food safety and public health on a global scale [13]. Although the ResFinder database included 22 groups of antimicrobial resistance, 7 resistance groups such as aminoglycoside, beta-lactamase, chloramphenicol, fosfomycin, quinolone, sulfonamide and tetracycline were mainly found in isolates. Prevalence of ARGs expressed as % of isolates containing ARGs were shown in Figures 2-4. When all isolates were included, host effects (p < 0.0001) on prevalence were observed for all ARGs (Figure 2(A)). Aminoglycoside, beta-lactamase and sulfonamide prevalence were greater in bovine and poultry than swine. Chloramphenicol prevalence was greater in bovine than in poultry and swine. Fosfomycin and tetracycline were greater in bovine and swine than poultry. Quinolone was greater in human than other hosts. Variation in AMR profile diversity both within and between food-production animal host species with higher AMR diversity of isolates from broiler compared to layer chickens, breeder compared to rearer and finisher pigs and beef compared to dairy cattle was reported and alternate sources of AMR bacteria and variation in selective evolutionary pressures within and between food-production animal host species populations was suggested [14].
Continent effects (p < 0.0001) were also observed on all ARGs (Figure 2(B)). Aminoglycoside, beta-lactamase and sulfonamide prevalence were greater in North America isolates than Europe and chloramphenicol, fosfomycin, quinolone and tetracycline were greater in Europe than North America isolates. Considering the evolution and dissemination of pathogens in the same lineages, horizontal gene transfer might occur locally and under different pressures with local  drug policies [15]. Clade 1 had greater (p < 0.0001) prevalence of aminoglycoside, beta-lactamase, and sulfonamide than Clade 3 (Figure 2(C)). Clade 3 had greater prevalence of chloramphenicol and quinolone than Clade 1. Clade 10 contained higher prevalence of fosfomycin and tetracycline than Clades 1 and 3.
Host sources were compared in each continent to exclude the continent effect on host effects (Figure 3). ARGs prevalence is similar in North America isolates to in all isolates except chloramphenicol and tetracycline (Figure 3(A)). Chloramphenicol prevalence in bovine was greater (p < 0.0001) than in human, and tetracycline prevalence were not significantly different between bovine and swine in North America isolates. In Europe isolates (Figure 3(B)), Human isolates had greater prevalence of aminoglycoside (p = 0.0111), beta-lactamase (p = 0.0214), chloramphenicol and quinolone (p < 0.0001), and sulfonamide (p = 0.0214) than swine, and swine had greater (p < 0.0001) prevalence of fosfomycin and tetracycline than human isolates. For human isolates (Figure 4(A)), Europe isolates had greater prevalence of chloramphenicol (p = 0.0038) and quinolone (p = 0.0019, and for swine isolates (Figure 4(B)) prevalence of aminoglycoside, beta-lactamase, sulfonamide were greater (p = 0.05) in North America isolates and the prevalence of tetracycline was greater (p < 0.0001) in Europe isolates.
Prevalence of bacteriophages (phages) was greater (p < 0.0001) in bovine and swine than poultry and human isolates ( Figure 5). Poultry also had the greater (p < 0.0001) prevalence of phage than human. Between continents, North America had greater (p = 0.0061) prevalence of phage than Europe ( Figure 5). Phages are viruses that target and destroy bacteria [16]. Research developing therapeutic phage use as antibiotics in animal industry is increasing due to their specificity, efficient bacterial lysis and their capability to self-replicate [17]. However, Salmonella phages appear to be abundant (up to 97% prevalence) and widespread in animal feces, wastewater, and the environment of food animals [18] and the  diversity of phages are thought to have interactions between Salmonella serovars and specific phages [19]. Salmonella phages had different host range that US phages showed a narrow host range, whereas Thai phages showed a wide phage host range [20].  [21]. The most of SPI is conserved in serovars and the different profile of virulence factors in S. Reading implies the loss or acquirement of genes to adapt to host sources and maintained in their lineage but not within the continent.
In conclusion, evolutionary distance is related to isolation host source rather than continent, and genome features were distinguished by host and cluster. Our genomic analysis implies that Salmonella Reading evolved independently to environments within its lineages.

Summary
Two large outbreaks of S. Reading with raw turkey products in North America have been reported, the study found that the majority of poultry strains are phylogenetically linked in multi-clades identified for North American poultry strains. Poultry stains clustered with swine stains and are polyphyletic and at least two lineages existed for poultry stains.
Aminoglycoside, beta-lactamase, chloramphenicol, fosfomycin, quinolone, sulfonamide and tetracycline ARGs were found in isolates. Bovine and poultry showed greater prevalence in aminoglycoside, beta-lactam and sulfonamide resistances than in swine. Chloramphenicol prevalence was greater in bovine than poultry and swine. Fosfomycin and tetracycline prevalences were greater in bovine and swine than poultry while quinolone was greater in humans than in any other hosts. Aminoglycoside, beta-lactam and sulfonamide prevalences were greater in North America than in Europe. Chloramphenicol, fosfomycin, quinolone and tetracycline were greater in Europe than in North America. Considering the evolution and dissemination of pathogens in the same lineages, horizontal gene transfer might occur locally under varying pressures associated with local drug policies.
Prevalence of bacteriophages in bovine and swine were greater than poultry and human. Poultry also had a greater prevalence of phage than humans. North America displayed a greater prevalence of phage than Europe. Europe strains showed higher diversity of virulence factors than North America stains in Clade From the study, the rapid worldwide dissemination of pathogens including S.
Reading can be attributed to increased human travel and large-scale and broad distribution of food and feed resources.