Temporal Genomic Phylogeny Reconstruction Indicates a Geospatial Transmission Path of Salmonella Cerro in the United States and a Clade-Specific Loss of Hydrogen Sulfide Production

Jasna Kovac, Kevin J. Cummings, Lorraine D. Rodriguez-Rivera, Laura M. Carroll, Anil Thachil, Martin Wiedmann

Salmonella Cerro has become one of the most prevalent Salmonella serotypes isolated from dairy cattle in several U.S. states, including New York where it represented 36% of all Salmonella isolates of bovine origin in 2015. This serotype is commonly isolated from dairy cattle with clinical signs of salmonellosis, including diarrhea and fever, although it has also been identified in herds without evidence of clinical disease or decreased production. To better understand the transmission patterns and drivers of its geographic spread, we have studied the genomic similarity and microevolution of S. Cerro isolates from the northeast U.S. and Texas. Eighty-three out of 86 isolates were confirmed as multilocus sequence type 367. We identified core genome SNPs in 57 upstate New York (NY), 2 Pennsylvania (PA), and 27 Texas S. Cerro isolates from dairy cattle, farm environments, raw milk, and one human clinical case and used them to construct a tip-dated phylogeny. S. Cerro isolates clustered in three distinct clades, including (i) clade I (n = 3; 2013) comprising isolates from northwest Texas (NTX), (ii) clade II (n = 14; 2009-2011, 2014) comprising isolates from NY, and (iii) clade III comprising isolates from NY, PA, and central Texas (CTX) in subclade IIIa (n = 45; 2008-2014), and only CTX isolates in subclade IIIb (n = 24; 2013). Temporal phylogenetic analysis estimated the divergence of these three clades from the most recent common ancestor in approximately 1980. The CTX clade IIIb was estimated to have evolved and diverged from the NY ancestor around 2004. Furthermore, gradual temporal loss of genes encoding a D-alanine transporter, involved in virulence, was observed. These genes were present in the isolates endemic to NTX clade I and were gradually lost in clades II and III. The virulence gene orgA, which is part of the Salmonella Pathogenicity Island 1, was lost in a subgroup of Texas isolates in clades I and IIIb. All S. Cerro isolates had an additional cytosine inserted in a cytosine-rich region of the virulence gene sopA, resulting in premature termination of translation likely responsible for loss of pathogenic capacity in humans. A group of closely related NY isolates was characterized by the loss of hydrogen sulfide production due to the truncation or complete loss of phsA. Our data suggest the ability of Salmonella to rapidly diverge and adapt to specific niches (e.g., bovine niche), and to modify virulence-related characteristics such as the ability to utilize tetrathionate as an alternative electron acceptor, which is commonly used to detect Salmonella. Overall, our results show that clinical outcome data and genetic data for S. Cerro isolates, such as truncations in virulence genes leading to novel pheno- and pathotypes, should be correlated to allow for accurate risk assessment.