Evolution of Salmonella enterica serotype Typhimurium driven by anthropogenic selection and niche adaptation.

Salmonella enterica serotype Typhimurium (S. Typhimurium) is a leading cause of gastroenteritis and bacteraemia worldwide, and a model organism for the study of host-pathogen interactions. Two S. Typhimurium strains (SL1344 and ATCC14028) are widely used to study host-pathogen interactions, yet genotypic variation results in strains with diverse host range, pathogenicity and risk to food safety. The population structure of diverse strains of S. Typhimurium revealed a major phylogroup of predominantly sequence type 19 (ST19) and a minor phylogroup of ST36. The major phylogroup had a population structure with two high order clades (α and ß) and multiple subclades on extended internal branches, that exhibited distinct signatures of host adaptation and anthropogenic selection. Clade α contained a number of subclades composed of strains from well characterized epidemics in domesticated animals, while clade ß contained multiple subclades associated with wild avian species. The contrasting epidemiology of strains in clade α and ß was reflected by the distinct distribution of antimicrobial resistance (AMR) genes, accumulation of hypothetically disrupted coding sequences (HDCS), and signatures of functional diversification. These observations were consistent with elevated anthropogenic selection of clade α lineages from adaptation to circulation in populations of domesticated livestock, and the predisposition of clade ß lineages to undergo adaptation to an invasive lifestyle by a process of convergent evolution with of host adapted Salmonella serotypes. Gene flux was predominantly driven by acquisition and recombination of prophage and associated cargo genes, with only occasional loss of these elements. The acquisition of large chromosomally-encoded genetic islands was limited, but notably, a feature of two recent pandemic clones (DT104 and monophasic S. Typhimurium ST34) of clade α (SGI-1 and SGI-4).

Genome Variation and Molecular Epidemiology of Salmonella enterica Serovar Typhimurium Pathovariants.

Salmonella enterica serovar Typhimurium is one of approximately 2,500 distinct serovars of the genus Salmonella but is exceptional in its wide distribution in the environment, livestock, and wild animals. S Typhimurium causes a large proportion of nontyphoidal Salmonella (NTS) infections, accounting for a quarter of infections, second only to S. enterica serovar Enteritidis in incidence. S Typhimurium was once considered the archetypal broad-host-range Salmonella serovar due to its wide distribution in livestock and wild animals, and much of what we know of the interaction of Salmonella with the host comes from research using a small number of laboratory strains of the serovar (LT2, SL1344, and ATCC 14028). But it has become clear that these strains do not reflect the genotypic or phenotypic diversity of S Typhimurium. Here, we review the epidemiological record of S Typhimurium and studies of the host-pathogen interactions of diverse strains of S Typhimurium. We present the concept of distinct pathovariants of S Typhimurium that exhibit diversity of host range, distribution in the environment, pathogenicity, and risk to food safety. We review recent evidence from whole-genome sequencing that has revealed the extent of genomic diversity of S Typhimurium pathovariants, the genomic basis of differences in the level of risk to human and animal health, and the molecular epidemiology of prominent strains. An improved understanding of the impact of genome variation of bacterial pathogens on pathogen-host and pathogen-environment interactions has the potential to improve quantitative risk assessment and reveal how new pathogens evolve.

Genomic analysis of clinical Aeromonas isolates reveals genetic diversity but little evidence of genetic determinants for diarrhoeal disease.

Aeromonas spp. are associated with a number of infectious syndromes in humans including gastroenteritis and dysentery. Our understanding of the genetic diversity, population structure, virulence determinants and antimicrobial resistance of the genus has been limited by a lack of sequenced genomes linked to metadata. We performed a comprehensive analysis of the whole genome sequences of 447 Aeromonas isolates from children in Karachi, Pakistan, with moderate-to-severe diarrhoea (MSD) and from matched controls without diarrhoea that were collected as part of the Global Enteric Multicenter Study (GEMS). Human-associated Aeromonas isolates exhibited high species diversity and extensive antimicrobial and virulence gene content. Aeromonas caviae, A. dhankensis, A. veronii and A. enteropelogenes were all significantly associated with MSD in at least one cohort group. The maf2 and lafT genes that encode components of polar and lateral flagella, respectively, exhibited a weak association with isolates originating from cases of gastroenteritis.

Population structure of Salmonella enterica Typhi in Harare, Zimbabwe (2012-19) before typhoid conjugate vaccine roll-out: a genomic epidemiology study.

BACKGROUND: The continued emergence of Salmonella enterica serovar Typhi, with ever increasing antimicrobial resistance, necessitates the use of vaccines in endemic countries. A typhoid fever outbreak in Harare, Zimbabwe, in 2018 from a multidrug resistant S Typhi with additional resistance to ciprofloxacin was the catalyst for the introduction of a typhoid conjugate vaccine programme. We aimed to investigate the emergence and evolution of antimicrobial resistance of endemic S Typhi in Zimbabwe and to determine the population structure, gene flux, and sequence polymorphisms of strains isolated before a typhoid conjugate vaccine programme to provide a baseline for future evaluation of the effect of the vaccination programme. METHODS: In this genomic epidemiology study, we used short-read whole-genome sequencing of S Typhi isolated from clinical cases of typhoid fever in Harare, Zimbabwe, between Jan 1, 2012, and Feb 9, 2019, to determine the S Typhi population structure, gene flux, and sequence polymorphisms and reconstructed the evolution of antimicrobial resistance. Maximum likelihood time-scaled phylogenetic trees of Zimbabwe isolates in the context of global isolates obtained from the National Center for Biotechnology Information were constructed to infer spread and emergence of antimicrobial resistance. FINDINGS: The population structure of S Typhi in Harare, Zimbabwe, from 2012 to 2019 was dominated by multidrug resistant genotype 4.3.1.1.EA1 (H58) that spread to Zimbabwe from neighbouring countries in around 2009 (95% credible interval 2008·5-2010·0). Acquisition of an IncN plasmid carrying antimicrobial resistance genes including a qnrS gene and a mutation in the quinolone resistance determining region of gyrA gene contributed to non-susceptibility and resistance to quinolone antibiotics. A minority population of antimicrobial susceptible S Typhi genotype 3.3.1 strains were present throughout. INTERPRETATION: The currently dominant S Typhi population is genotype 4.3.1.1 that spread to Zimbabwe and acquired additional antimicrobial resistance though acquisition of a plasmid and mutation in the gyrA gene. This study provides a baseline population structure for future evaluation of the effect of the typhoid conjugate vaccine programme in Harare. FUNDING: Bill & Melinda Gates Foundation and the Biotechnology and Biological Sciences Research Council Institute Strategic Programme.

A genomic appraisal of invasive Salmonella Typhimurium and associated antibiotic resistance in sub-Saharan Africa.

Invasive non-typhoidal Salmonella (iNTS) disease manifesting as bloodstream infection with high mortality is responsible for a huge public health burden in sub-Saharan Africa. Salmonella enterica serovar Typhimurium (S. Typhimurium) is the main cause of iNTS disease in Africa. By analysing whole genome sequence data from 1303 S. Typhimurium isolates originating from 19 African countries and isolated between 1979 and 2017, here we show a thorough scaled appraisal of the population structure of iNTS disease caused by S. Typhimurium across many of Africa's most impacted countries. At least six invasive S. Typhimurium clades have already emerged, with ST313 lineage 2 or ST313-L2 driving the current pandemic. ST313-L2 likely emerged in the Democratic Republic of Congo around 1980 and further spread in the mid 1990s. We observed plasmid-borne as well as chromosomally encoded fluoroquinolone resistance underlying emergences of extensive-drug and pan-drug resistance. Our work provides an overview of the evolution of invasive S. Typhimurium disease, and can be exploited to target control measures.

Single-cell genomics reveals population structures from in vitro evolutionary studies of Salmonella.

Single-cell DNA sequencing has the potential to reveal detailed hierarchical structures in evolving populations of cells. Single cell approaches are increasingly used to study clonal evolution in human ageing and cancer but have not yet been deployed to study evolving clonal microbial populations. Here, we present an approach for single bacterial genomic analysis for in vitro evolution experiments using FACS isolation of individual bacteria followed by whole-genome amplification and sequencing. We apply this to the experimental evolution of a hypermutator strain of Salmonella in response to antibiotic stress (ciprofloxacin). By analysing sequence polymorphisms in individual cells from populations we identified the presence and prevalence of sub-populations which have acquired polymorphisms in genes previously demonstrated to be associated with ciprofloxacin susceptibility. We were also able to identify that the population exposed to antibiotic stress was able to develop resistance whilst maintaining diversity. This population structure could not be resolved from bulk sequence data, and our results show how high-throughput single-cell sequencing can enhance experimental studies of bacterial evolution.

Increased phage resistance through lysogenic conversion accompanying emergence of monophasic Salmonella Typhimurium ST34 pandemic strain.

Salmonella enterica serovar Typhimurium (S. Typhimurium) comprises a group of closely related human and animal pathogens that account for a large proportion of all Salmonella infections globally. The epidemiological record of S. Typhimurium in Europe is characterized by successive waves of dominant clones, each prevailing for approximately 10-15 years before replacement. Succession of epidemic clones may represent a moving target for interventions aimed at controlling the spread and impact of this pathogen on human and animal health. Here, we investigate the relationship of phage sensitivity and population structure of S. Typhimurium using data from the Anderson phage typing scheme. We observed greater resistance to phage predation of epidemic clones circulating in livestock over the past decades compared to variants with a restricted host range implicating increased resistance to phage in the emergence of epidemic clones of particular importance to human health. Emergence of monophasic S. Typhimurium ST34, the most recent dominant multidrug-resistant clone, was accompanied by increased resistance to phage predation during clonal expansion, in part by the acquisition of the mTmII prophage that may have contributed to the fitness of the strains that replaced ancestors lacking this prophage.

Ecological niche adaptation of Salmonella Typhimurium U288 is associated with altered pathogenicity and reduced zoonotic potential.

The emergence of new bacterial pathogens is a continuing challenge for agriculture and food safety. Salmonella Typhimurium is a major cause of foodborne illness worldwide, with pigs a major zoonotic reservoir. Two phylogenetically distinct variants, U288 and ST34, emerged in UK pigs around the same time but present different risk to food safety. Here we show using genomic epidemiology that ST34 accounts for over half of all S. Typhimurium infections in people while U288 less than 2%. That the U288 clade evolved in the recent past by acquiring AMR genes, indels in the virulence plasmid pU288-1, and accumulation of loss-of-function polymorphisms in coding sequences. U288 replicates more slowly and is more sensitive to desiccation than ST34 isolates and exhibited distinct pathogenicity in the murine model of colitis and in pigs. U288 infection was more disseminated in the lymph nodes while ST34 were recovered in greater numbers in the intestinal contents. These data are consistent with the evolution of S. Typhimurium U288 adaptation to pigs that may determine their reduced zoonotic potential.

Whole-genome epidemiology links phage-mediated acquisition of a virulence gene to the clonal expansion of a pandemic Salmonella enterica serovar Typhimurium clone.

Epidemic and pandemic clones of bacterial pathogens with distinct characteristics continually emerge, replacing those previously dominant through mechanisms that remain poorly characterized. Here, whole-genome-sequencing-powered epidemiology linked horizontal transfer of a virulence gene, sopE, to the emergence and clonal expansion of a new epidemic Salmonella enterica serovar Typhimurium (S. Typhimurium) clone. The sopE gene is sporadically distributed within the genus Salmonella and rare in S. enterica Typhimurium lineages, but was acquired multiple times during clonal expansion of the currently dominant pandemic monophasic S. Typhimurium sequence type (ST) 34 clone. Ancestral state reconstruction and time-scaled phylogenetic analysis indicated that sopE was not present in the common ancestor of the epidemic clade, but later acquisition resulted in increased clonal expansion of sopE-containing clones that was temporally associated with emergence of the epidemic, consistent with increased fitness. The sopE gene was mainly associated with a temperate bacteriophage mTmV, but recombination with other bacteriophage and apparent horizontal gene transfer of the sopE gene cassette resulted in distribution among at least four mobile genetic elements within the monophasic S. enterica Typhimurium ST34 epidemic clade. The mTmV prophage lysogenic transfer to other S. enterica serovars in vitro was limited, but included the common pig-associated S. enterica Derby (S. Derby). This may explain mTmV in S. Derby co-circulating on farms with monophasic S. Typhimurium ST34, highlighting the potential for further transfer of the sopE virulence gene in nature. We conclude that whole-genome epidemiology pinpoints potential drivers of evolutionary and epidemiological dynamics during pathogen emergence, and identifies targets for subsequent research in epidemiology and bacterial pathogenesis.

SGI-4 in Monophasic Salmonella Typhimurium ST34 Is a Novel ICE That Enhances Resistance to Copper.

A multi drug resistant Salmonella enterica 4,[5],12:i- of sequence type 34 (monophasic S. Typhimurium ST34) is a current pandemic clone associated with livestock, particularly pigs, and numerous outbreaks in the human population. A large genomic island, termed SGI-4, is present in the monophasic Typhimurium ST34 clade and absent from other S. Typhimurium strains. SGI-4 consists of 87 open reading frames including sil and pco genes previously implicated in resistance to copper (Cu) and silver, and multiple genes predicted to be involved in mobilization and transfer by conjugation. SGI-4 was excised from the chromosome, circularized, and transferred to recipient strains of S. Typhimurium at a frequency influenced by stress induced by mitomycin C, and oxygen tension. The presence of SGI-4 was associated with increased resistance to Cu, particularly but not exclusively under anaerobic conditions. The presence of silCBA genes, predicted to encode an RND family efflux pump that transports Cu from the periplasm to the external milieu, was sufficient to impart the observed enhanced resistance to Cu, above that commonly associated with S. Typhimurium isolates. The presence of these genes resulted in the absence of Cu-dependent induction of pco genes encoding multiple proteins linked to Cu resistance, also present on SGI-4, suggesting that the system effectively limits the Cu availability in the periplasm, but did not affect SodCI-dependent macrophage survival.

Microevolution of antimicrobial resistance and biofilm formation of Salmonella Typhimurium during persistence on pig farms.

Salmonella Typhimurium and its monophasic variant S. 4,[5],12:i:- are the dominant serotypes associated with pigs in many countries. We investigated their population structure on nine farms using whole genome sequencing, and their genotypic and phenotypic variation. The population structure revealed the presence of phylogenetically distinct clades consisting of closely related clones of S. Typhimurium or S. 4,[5],12:i:- on each pig farm, that persisted between production cycles. All the S. 4,[5],12:i:- strains carried the Salmonella genomic island-4 (SGI-4), which confers resistance to heavy metals, and half of the strains contained the mTmV prophage, harbouring the sopE virulence gene. Most clonal groups were highly drug resistant due to the presence of multiple antimicrobial resistance (AMR) genes, and two clades exhibited evidence of recent on-farm plasmid-mediated acquisition of additional AMR genes, including an IncHI2 plasmid. Biofilm formation was highly variable but had a strong phylogenetic signature. Strains capable of forming biofilm with the greatest biomass were from the S. 4,[5],12:i:- and S. Typhimurium DT104 clades, the two dominant pandemic clones found over the last 25 years. On-farm microevolution resulted in enhanced biofilm formation in subsequent production cycle.

Genome-wide Analysis of Salmonella enterica serovar Typhi in Humanized Mice Reveals Key Virulence Features.

Salmonella enterica serovar Typhi causes typhoid fever only in humans. Murine infection with S. Typhimurium is used as a typhoid model, but its relevance to human typhoid is limited. Non-obese diabetic-scid IL2rγnull mice engrafted with human hematopoietic stem cells (hu-SRC-SCID) are susceptible to lethal S. Typhi infection. In this study, we use a high-density S. Typhi transposon library in hu-SRC-SCID mice to identify virulence loci using transposon-directed insertion site sequencing (TraDIS). Vi capsule, lipopolysaccharide (LPS), and aromatic amino acid biosynthesis were essential for virulence, along with the siderophore salmochelin. However, in contrast to the murine S. Typhimurium model, neither the PhoPQ two-component system nor the SPI-2 pathogenicity island was required for lethal S. Typhi infection, nor was the CdtB typhoid toxin. These observations highlight major differences in the pathogenesis of typhoid and non-typhoidal Salmonella infections and demonstrate the utility of humanized mice for understanding the pathogenesis of a human-specific pathogen.

Y Chromosome Sequences Reveal a Short Beringian Standstill, Rapid Expansion, and early Population structure of Native American Founders.

The Americas were the last inhabitable continents to be occupied by humans, with a growing multidisciplinary consensus for entry 15-25 thousand years ago (kya) from northeast Asia via the former Beringia land bridge [1-4]. Autosomal DNA analyses have dated the separation of Native American ancestors from the Asian gene pool to 23 kya or later [5, 6] and mtDNA analyses to ∼25 kya [7], followed by isolation ("Beringian Standstill" [8, 9]) for 2.4-9 ky and then a rapid expansion throughout the Americas. Here, we present a calibrated sequence-based analysis of 222 Native American and relevant Eurasian Y chromosomes (24 new) from haplogroups Q and C [10], with four major conclusions. First, we identify three to four independent lineages as autochthonous and likely founders: the major Q-M3 and rarer Q-CTS1780 present throughout the Americas, the very rare C3-MPB373 in South America, and possibly the C3-P39/Z30536 in North America. Second, from the divergence times and Eurasian/American distribution of lineages, we estimate a Beringian Standstill duration of 2.7 ky or 4.6 ky, according to alternative models, and entry south of the ice sheet after 19.5 kya. Third, we describe the star-like expansion of Q-M848 (within Q-M3) starting at 15 kya [11] in the Americas, followed by establishment of substantial spatial structure in South America by 12 kya. Fourth, the deep branches of the Q-CTS1780 lineage present at low frequencies throughout the Americas today [12] may reflect a separate out-of-Beringia dispersal after the melting of the glaciers at the end of the Pleistocene.