Whole genome sequence analysis of Salmonella Typhi provides evidence of phylogenetic linkage between cases of typhoid fever in Santiago, Chile in the 1980s and 2010-2016 ========================================================================================================================================================================= * Mailis Maes * Michael J. Sikorski * Megan E. Carey * Ellen E. Higginson * Zoe A. Dyson * Alda Fernandez * Pamela Araya * Sharon M. Tennant * Stephen Baker * Rosanna Lagos * Juan Carlos Hormazábal * Myron M. Levine * Gordon Dougan ## Abstract Typhoid fever epidemiology was investigated rigorously in Santiago, Chile during the 1980s, when *Salmonella enterica* serovar Typhi (*S*. Typhi) caused seasonal, hyperendemic disease. Targeted interventions reduced the annual typhoid incidence rates from 128-220 cases/105 population occurring between 1977-1984 to <8 cases/105 from 1992 onwards. As such, Santiago represents a contemporary example of the epidemiologic transition of an industrialized city from amplified hyperendemic typhoid fever to a period when typhoid is no longer endemic. We used whole genome sequencing (WGS) and phylogenetic analysis to compare the genotypes of *S*. Typhi cultured from acute cases of typhoid fever occurring in Santiago during the hyperendemic period of the 1980s (n=74) versus the nonendemic 2010s (n=80) when typhoid fever was rare. The genotype distribution between “historical” (1980s) isolates and “modern” (2011-2016) isolates was similar, with genotypes 3.5 and 2 comprising the majority of isolations, and 73/80 (91.3%) of modern isolates matching a genotype detected in the 1980s. Additionally, phylogenomically ‘ancient’ genotypes 1.1 and 1.2.1, uncommon in the global collections, were also detected in both eras, with a notable rise amongst the modern isolates. Thus, genotypes of *S*. Typhi causing acute illness in the modern nonendemic era match the genotypes circulating during the hyperendemic 1980s. The persistence of historical genotypes may be explained by chronic typhoid carriers originally infected during or before the 1980s. **Author summary** Studies of *Salmonella* Typhi (the cause of typhoid fever) rarely include isolates collected both before and after the interruption of hyperendemic transmission because this typically occurred decades before modern bacteria preservation methods. After interruption, it was assumed that sporadic cases and infrequent outbreaks were due to either chronic biliary carriers or importations, but this was difficult to characterize with low resolution bacterial typing methods. In Santiago, Chile, typhoid fever persisted at hyperendemic levels through the 1980s until organized control efforts in the 1980s and changes to wastewater policy in 1991 caused annual typhoid incidence to plummet. In this study, we used whole genome sequencing (WGS) to investigate whether recent sporadic cases occurring in Santiago in the 2010s were genomically similar to *S*. Typhi circulating in the 1980s, or dissimilar, possibly representing importations of *S*. Typhi from outside of Chile. We found concordance amongst *S*. Typhi genotypes between the 1980s and 2010s, and differences from genotypes circulating in Southeast Asia and Africa where typhoid remains hyperendemic. Our findings suggest that a proportion of modern, rare typhoid cases in Santiago are autochthonous, and that chronic carriers or another unknown reservoir likely contribute. Broadly, our findings corroborate the epidemiologic importance of long-term reservoirs of typhoid fever decades after typhoid elimination. Keywords * *S*. Typhi * Chile * phylogenetics * historical lineages * South America ## Introduction Typhoid fever is a systemic infection, caused by the bacterium *Salmonella enterica* subspecies *enterica* serovar Typhi (*S*. Typhi), which is spread through ingestion of food and water contaminated with human faeces containing *S*. Typhi. Typhoid fever in Santiago, Chile rose to hyper-endemic levels in the mid-1970s with a striking annual incidence of 177-220 cases per 105 population per year from 1977-1983, partially declined during the years 1984-1991 (51-128 cases/105 per year), and finally plummeted below 8 cases/105 per year after 1991, reaching nonendemic levels of ∼0.5/105 in 2012 [1]. During the hyperendemic years, the highest incidence rate was observed in school age children, and case numbers peaked in the rainless summer months [2]. To address this public health issue, in 1979 the Ministry of Health of Chile implemented a Typhoid Fever Control Program involving strategic preventive interventions including large-scale field trials of live oral typhoid vaccine Ty21a in >500,000 schoolchildren [3-6], and identification and treatment of chronic carriers among food handlers [7]. Collectively, these interventions led to a steady decline in typhoid cases between 1983 and 1990 [3-6, 8]. Epidemiologic and environmental bacteriology studies in the 1980s undertaken by the Chilean Typhoid Fever Control Program also incriminated the use of raw untreated sewage water for irrigation of vegetable crops during the rainless summer months as a key risk factor responsible for amplified transmission of *S*. Typhi in Santiago [8, 9]. In April 1991, the government of Chile abruptly prohibited the practice of irrigating crops with untreated sewage wastewater in Santiago [1, 8, 10, 11]. This action was taken following the enormous explosive epidemics of El Tor cholera in 1991 that occurred in neighbouring countries to the North, including in Peru, beginning in late January [12, 13], followed by Ecuador in February [13], and Colombia in March [13]. The pivotal stimulus for the intervention was a small outbreak of cholera in Santiago, Chile, in April, 1991 [1, 8, 10, 11]. Following imposition and strict enforcement of the prohibition of use of untreated sewage water for irrigation, the Santiago cholera epidemic quickly came to an end [10, 11], and the incidence of typhoid dropped precipitously in subsequent years [1, 8, 10]. From 2000-2014, the annual incidence of confirmed typhoid fever in Santiago has ranged from 0.20-1.18 cases/105 population. The current rare sporadic cases of typhoid in Santiago include illnesses among: travellers who visit countries with high incidence rates of typhoid; immigrants from typhoid-endemic countries; and transmission from Chilean chronic carriers who were originally infected decades earlier (pre-1991) during the era of hyperendemic transmission. Phylogenetic analysis of whole genome sequencing (WGS) data has become the gold standard for inferring relatedness between *S*. Typhi isolates [14]. WGS can also be used to assign genotypes, identify and characterize plasmids, compare gene content and homology, and identify molecular determinants of antibiotic resistance and virulence [15, 16]. Whereas WGS data from many global isolates are available, few bacterial collections include isolates from both before and after the interruption of amplified hyperendemic typhoid fever. Also, WGS-derived information concerning the genotypes isolated from typhoid fever cases in South American countries is scarce. The majority of publicly available sequenced isolates come from hyperendemic regions in South Asia and East Africa where H58 (4.3.1) genotypes dominate [17]. Herein we describe a genomic analysis of *S*. Typhi isolates from Santiago, Chile from two distinct epidemiologic periods including the early 1980s (N=74), when typhoid was hyperendemic, and from 2010-2016 (N=80), a quarter century after the interruption of amplified transmission of *S*. Typhi in Santiago. We show that the genotypic composition of the population structure is similar between these two periods, even though the local epidemiology has changed drastically. In addition, we describe a potential carrier-related cluster of the “ancient” 1.1 clade, and the identification of a new variant of plasmid pHCM2 carried by 2 isolates of another ancient genotype (1.2.1). ## Material and methods ### Study design In total, 74 “historical” isolates of *S*. Typhi from cases of acute typhoid fever in Santiago, Chile in the period 1981 to 1986, during the era of hyperendemic typhoid, were obtained from the culture collection of the Center for Vaccine Development and Global Health of the University of Maryland, Baltimore. “Modern” *S*. Typhi isolates (n=80) were strains collected from acute cases of typhoid fever and referred to the Instituto de Salud Pública de Chile (ISP, Institute of Public Health) by clinical microbiology laboratories from 2011-2016; these referred strains were generally accompanied by only limited demographic and clinical information. ### DNA extraction and sequencing DNA was isolated from the 154 *S*. Typhi isolates using the Norgen DNA extraction kit as per the manufacturer’s instructions. Genomic DNA was sequenced by Illumina HiSeq at the Wellcome Sanger Institute, generating 125 bp paired reads. Raw read data were deposited in the European Nucleotide Archive and individual accession numbers are listed in Table S1. ### Assembly, mapping, SNP-calling and genotype assignment Raw Illumina reads were assembled using Velvet v1.2 via the Wellcome Sanger Institute automated analysis pipeline [18]. Sequenced reads were mapped and annotated, single nucleotide polymorphisms (SNPs) were called against the *S*. Typhi CT18 reference genome [19] using SMALT v0.7.4 and subsequent variant calling was carried out as previously described [20]. Genotypes were assigned using GenoTyphi (v1.9.0) [14, 16]. ### Phylogenetic analysis Previously defined [21-23] recombinant regions such as prophages and plasmids were excluded manually, and any remaining recombinant regions were filtered using Gubbins (v2.5.0) [24]. The resultant core genome alignment was used to infer Maximum Likelihood (ML) phylogenies using RAxML (v8.2.8) [25], specifying a generalized time-reversible model and a gamma distribution to model site-specific rate variation (GTR+ G substitution model; GTRGAMMA in RAxML) with 100 bootstrap pseudoreplicates used to assess branch support. SNP distances for the core genome alignment of all the strains have also been calculated from this alignment using snp-dists package (available at: [https://github.com/tseemann/snp-dists](https://github.com/tseemann/snp-dists)). Resulting phylogenies were visualized with iTOL [26] and are available for interactive inspection on Microreact [27] ([https://microreact.org/project/j25cGgoVX8ixQCRstPVhgo-chile-1980-2010](https://microreact.org/project/j25cGgoVX8ixQCRstPVhgo-chile-1980-2010)). ### Molecular determination of AMR, plasmids and virulence genes Raw read data were screened using SRST2 v0.2.0 [15] with the ARG-ANNOT [28], PlasmidFinder [29], and *Salmonella* Virulence Factor DB databases [30] to detect known and putative antimicrobial resistance (AMR) genes, known plasmid replicons, and known virulence-associated genes, respectively. Homology cutoffs of 90% nucleotide similarity were used for all SRST2 screens. Plasmid pHCM2 sequences were assembled with Unicycler [31], and the circularised sequences extracted from *de bruijn* assembly graphs with Bandage (v0.8.1) [32]. These were then compared with the pHCM2 reference sequence (accession no: [AL513384](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=AL513384&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom)) using Artemis-ACT [33] and visualized with EasyFig [34]. ## Results ### *S*. Typhi genotypes in Santiago, Chile Genotypic analysis indicated that most of the isolates from Santiago in this study belonged to genotypes 2 (n=43, 27.9%), 3.5 (n=51, 33.1%) and 1.1 (n=22, 14.3%) (Table 1). While these three genotypes were the most prevalent, we also detected 11 additional distinct genotypes across the two time periods. These include three previously described importations of isolates belonging to H58 genotypes 4.3.1.1 and 4.3.1.2 [35]. View this table: [Table 1.](http://medrxiv.org/content/early/2022/01/24/2022.01.19.22269577/T1) Table 1. S. Typhi genotypes detected in isolates from acute clinical cases in Santiago, Chile during the hyperendemic 1980s versus the nonendemic 2010s The genomic population structure of *S*. Typhi in Santiago appeared highly consistent between the hyperendemic and post-hyperendemic eras (Fig 1). Both rare “ancient” genotypes 1.1 and 1.2.1, as well as more commonly found genotypes 2 and 3.5, were frequently present across both study periods. Despite several decades having passed since the cessation of the era of hyperendemic transmission, most of the post-2010 sporadic *S*. Typhi isolates shared bootstrap-supported phylogenetic clusters with the 1980 hyperendemic isolates. Only within genotype 1.1 did recent and historical sequences form independent clusters. The intermingling of historical and recent sequences suggests that the 2010s isolates originated from progenitors within the same pool of *S*. Typhi circulating in the 1980s. ![Fig. 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/24/2022.01.19.22269577/F1.medium.gif) [Fig. 1.](http://medrxiv.org/content/early/2022/01/24/2022.01.19.22269577/F1) Fig. 1. Genotypes and molecular determinates of AMR and plasmids in the Santiago collection. SNP based maximum likelihood core genome phylogeny indicating the genotypes of all the Santiago *S*. Typhi isolates is shown by the branch colour as per the inset legend. A diverse range of 14 genotypes were detected including 1.1, 1.2.1, 2, 2.0.2, 2.2, 2.3.2, 2.3.3, 2.3.4, 3.1, 3.3, 3.5, 4.1, 4.3.1.1, 4.3.1.2. Most of the isolates belonged to genotype 1.1, 2 and 3.5. Year of isolation, SNPs in the QRDR region associated with reduced susceptibility to fluoroquinolones as well as plasmids and detected resistance genes are shown in the columns adjacent to the tree. Bootstrap confidence of 80 or above is indicated with a black circle on the branch. Amongst the 151 non-H58 isolates, few antimicrobial resistance mechanisms were detected. One genotype 2.0.0 isolate carried a *tetB* gene, known to encode a mechanism of resistance to tetracyclines (albeit not an antibiotic of clinical relevance for treatment of typhoid fever). Two genotype 2.2 isolates and one genotype 4.1 isolate had single non-synonymous point mutations in the Quinolone Resistance Determining Region (QRDR) of genes *gyrA* and *gyrB*, associated with reduced susceptibility to fluoroquinolones [36]; in the late 1980s, ciprofloxacin was the antibiotic of choice for treating patients with typhoid fever. Only 17 candidate plasmids were detected across all the isolates, including a Colicin plasmid, three related to plasmid pHCM2 and 13 of incompatibility type IncX1 (Fig 1). ### Global transmission dynamics To better understand how this local population structure fits into the global context of *S*. Typhi, a maximum likelihood phylogenetic tree was generated using a total of 2,013 non-H58 *S*. Typhi genomes from approximately 70 countries (Fig 2). The *S*. Typhi isolates from Santiago of genotypes 3.5, 2, 1.1 and 1.2.1 clustered together and were genetically distinct from those sequenced previously from other countries. The other genotypes present in Santiago were mostly found mostly on long branches, indicating significant differentiation from their closest relatives in the global context phylogeny. ![Fig. 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/24/2022.01.19.22269577/F2.medium.gif) [Fig. 2.](http://medrxiv.org/content/early/2022/01/24/2022.01.19.22269577/F2) Fig. 2. A global context maximum likelihood core genome phylogeny incorporating publicly available non-H58 and Santiago *S*. Typhi isolates. SNP based core genome maximum likelihood tree of 2013 global isolates, branches are coloured by genotype as per the inset legend; the ring indicates country of isolation as per the inset legend; Santiago isolates are highlighted in black and the Colombian isolates in red. ### “Ancient” *S*. Typhi clades endemic to Chile Two of the clades common in our dataset, 1.1 and 1.2, are found closer to the root of the global *S*. Typhi tree and thus can be considered to comprise “ancient” clades, which are detected with very low global frequency. Of the 1.1 isolates (n=22), two were from the 1980s collection (Fig 3). Only six other isolates, 0.3% of the global collection described by Wong et al [14], belonged to clade 1.1. Interestingly, a cluster of closely related genotype 1.1 *S*. Typhi isolates was detected in the years 2011 (1 isolate) and 2012 (19 isolates). These isolates also exhibit lower genetic variation (median SNP distance of 8) than those observed among the two 1980s 1.1 isolates (119 SNPs). The smallest SNP distance between a 1980s 1.1 isolate and a 2010s 1.1 isolate was 79 SNPs. This suggests the 2011-2012 isolates have a recent common ancestor. ![Fig. 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/24/2022.01.19.22269577/F3.medium.gif) [Fig. 3.](http://medrxiv.org/content/early/2022/01/24/2022.01.19.22269577/F3) Fig. 3. *S*. Typhi clade 1.1 phylogeny showing sequences carrying *ssrA* mutation (*ssrA-* G639E). Phylogenetic core genome analysis of all *S*. Typhi isolates of genotype 1.1 identified from the global collection. The majority of these isolates were collected in Santiago (black bar), two isolates belonged to the 1980s collection appearing ancestral to the 2010s collection. The twenty isolates from 2011-2012 showed an interesting cluster with a median SNP distance of 8 SNPs which is unlikely to have evolved among acute cases, where the mutation rate is calculated to be < 1 SNP/genome/year, in one year. Also shown is an alignment of a section of the *ssrA* gene, showing the conserved SNP at position 1648986 of the *S*. Typhi CT18 reference genome. * the precise isolation year of one of the 1980s isolates was unclear. We sought to determine if there was a specific virulence gene associated with the 2011-2012 genotype 1.1 cluster, using SRST2 with the *Salmonella* VFDB. Although there were no additional or missing virulence factor-associated genes detected, this screening did reveal that all of the genotype 1.1 isolates from 2011-2012 cluster harboured a SNP in the *ssrA* gene. The alignment of this gene revealed the same missense mutation (G to A at position 1648968 of CT18, causing a glycine to glutamic acid amino acid change) in the *ssrA* gene (*ssrA-* G639E), (Fig 3). The other ancient clade present amongst the Santiago isolates was genotype 1.2.1. This was the only clade in our datasets (except from one recently imported H58 isolate) in which isolates harboured pHCM2 plasmids. pHCM2 plasmids are found in ∼15% of *S*. Typhi and are generally highly conserved (Fig 4A). To our knowledge based on the data available in the *S*. Typhi public collection on Pathogenwatch: ([https://pathogen.watch/collection/1y3cdsjq55hf-public-genomes-accessedDec2021](https://pathogen.watch/collection/1y3cdsjq55hf-public-genomes-accessedDec2021)) pHCM2 plasmids have not been detected in ancient clades 0.0 to 1.2 previously. Two pHCM2 plasmids were found among our isolates from genotype 1.2.1; one in a 1984 isolate (ERR2710544) and the other from 2013 (ERR4289185). These plasmids were highly similar to each other, but significantly different to the pHCM2 found on *S*. Typhi reference CT18 (accession no: [AL513384](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=AL513384&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom)), and to another novel pHCM2 found in a 1973 isolate from Thailand [37]. The pHCM2 sequences described here differed by >2000 SNPs (98.4% nucleotide identity, 89% coverage) from the reference (Fig 4). ![Fig. 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/01/24/2022.01.19.22269577/F4.medium.gif) [Fig. 4.](http://medrxiv.org/content/early/2022/01/24/2022.01.19.22269577/F4) Fig. 4. Genetic comparison of cryptic plasmid pHCM2-like sequences. Global representative collection of pHCM2 aligned to the reference pHCM2 from CT18. (A) Presence and absence of genes located on pHCM2. Red indicates presence of the genes located on the reference pHCM2 sequence from CT18, while blue indicates absence of the genes compared to the reference. Nucleotide BLAST comparison of the 2 Chilean assembled pHCM2-like plasmids (ERR271054, ERR4289185) against the reference pHCM2 plasmid from *S*. Typhi CT18 (accession no:[AL513384](http://medrxiv.org/lookup/external-ref?link_type=GEN&access_num=AL513384&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom)) (B). Gene annotations are shown for the reference pHCM2 from CT18, with the shading intensity indicating percent nucleotide homology between the plasmid sequences. ## Discussion This study utilizes WGS analyses to compare *S*. Typhi isolates from patients who developed acute typhoid fever in Santiago during the modern era (i.e., 2011-2016, decades post-interruption of amplified typhoid transmission in 1991) with isolates from acute cases during the hyper-endemic era of the 1980s. Using phylogenetics and genome-wide screening, we compared the genotypes, AMR determinants, plasmid content of *S*. Typhi isolates between these two time periods. Whilst comparisons of isolates from the 1980’s with isolates from 1990 (post vaccinating interventions-but pre-irrigation ban, when incidences were still high 47 cases /105 [1]) have been undertaken previously with pulse field gel electrophoresis (PFGE) [38], to our knowledge, this is the first WGS-based analysis comparing *S*. Typhi collected during a period of hyper-endemicity and decades after typhoid control. The observation that the 1980s and 2010s collections of *S*. Typhi isolates do not cluster independently within our commonly found genotypes (2, 3.5 and 1.2.1) suggests that many of the sporadic cases detected in the 2010s share common ancestry with *S*. Typhi cases from the 1980s. These may, in turn, be linked to shedding of *S*. Typhi from chronic typhoid carriers. This is expected decades after hyper-endemic typhoid is controlled, as chronic carriers will continue to shed sporadically until they are treated or die out of the population [39]. Although there have been several detected introductions of the globally dominant H58 *S*. Typhi clade (4.3.1) into Santiago, as occurs with most industrialized countries in North America, Europe, and Australia, this was not a frequent genotype found in the sporadic modern Chilean cases. This is distinct from countries such as the United Kingdom where imported cases of H58 *S*. Typhi have become common among travel-associated typhoid cases in recent years [40, 41]. *S*. Typhi WGS data from a recent study in Colombia found genotypes 2, 2.5 and 3.5 to be the local endemic genotypes [42]. Genotype 2.5 was not detected among the Santiago isolates; however, genotype 3.5 and 2 were detected in Santiago, and were the nearest genomic neighbours of the Colombian isolates on the global tree. Even though the sequenced isolates from Colombia and Santiago, Chile both contained isolates form the 2010s and historical isolates, they clustered separately on the global tree, suggesting that there has not been obvious transmission between the two countries. More representative WGS data from other Central and South American countries are needed to determine if this is the case with other unsampled endemic countries in the area. One of the more interesting findings from this dataset was the cluster of 20 highly related strains of genotype 1.1 in 2011-2012. These isolates differed by a median of 8 SNPs in the core genome. For context, the average core genome mutation rate in *S*. Typhi is approximately 1 substitution/genome/year every year [21, 43-45]. As such, it is unlikely that these isolates originated from a clonal source and sequentially acquired this genetic variation in <2 years. Conversely, these strains are clearly highly related and were all collected during a short time span. A potential explanation to account for this cluster is that it was an outbreak caused by bacteria shed from a long-term carrier. *S*. Typhi colonizing the gallbladder or other long-term physiologic niche of a chronic typhoid carrier may evolve and mutate over time into a heterogeneous, yet closely related, community of bacteria [46]. Thus, the closely related nature of these isolates occurring mainly over a single year and during a period when typhoid was no longer hyperendemic in Santiago could be explained by an outbreak instigated by a chronic carrier who was shedding multiple highly related strains. If this hypothetical carrier was working as a food handler, intermittently contaminating food vehicles consumed by subjects during 2011-12, this could result in the cluster of closely related isolates we observe. The cessation of cases could then be explained by a change of employment or residence, a cholecystectomy procedure, a course of antibiotics that concentrate in the gallbladder (e.g., ciprofloxacin) prescribed for another infection or surgery, or death of the hypothetical carrier at the end of 2012. Notably, the only detectable difference in virulence factor-associated genes from this potential carrier-related 2011-2012 cluster is the *ssrA* missense mutation. The *ssrA/B* two component response regulator is a principle regulator of *Salmonella* pathogenicity island-2 (SPI-2) [47, 48]. SPI-2 has been shown to be insignificant in *S*. Typhi macrophage survival capability and *S*. Typhi’s infection efficiency of humanized mice [49, 50]. The *ssrA* gene has been linked to biofilm formation in *Salmonella* and persistence on gallstones [51, 52]. This observation may be worth additional investigation to investigate the role of SPI-2 and *ssrA* in Typhoid carriage. Our analysis identified a novel variant of cryptic plasmid pHCM2 in two sequences from genotype 1.2.1. A pHCM2 type plasmid was present in the *S*. Typhi reference isolate CT18 from Vietnam [53] and the plasmid family was described in some detail [54]. Interestingly, any role for pHCM2 in *S*. Typhi biology remains elusive. We know that almost all of the pHCM2 plasmids in the global *S*. Typhi collection are highly conserved in terms of overall nucleotide diversity. However, the Chilean pHCM2 from clade 1.2.1 is substantially different from most other family members. The only other reported pHCM2 plasmid with a similar level of variation is harboured in a *S*. Typhi Thailand isolate from 1973 [55]. In conclusion, this genomics analysis suggests that asymptomatic typhoid carriers (or other unknown reservoirs) are likely responsible for continuing sporadic typhoid cases in Chile. The ancient clades present in the country have persisted. We recommend that all sporadic typhoid cases and clusters should be evaluated using epidemiological outbreak investigations and WGS as routine surveillance. This will help ensure the detection of any potential introduction or development of drug-resistant *S*. Typhi and may lead to the identification and treatment of carriers, which would further limit ongoing disease transmission. ## Data Availability All Raw Sequencing files are available from the ENA database (accession numbers are listed in Table S1). ## Authorship declaration All the Co-authors contributed to manuscript editing ## Authors contributions **Conception and design of experiments** Myron M. Levine, Rosanna Lagos, **Sample collection and DNA extractions** Alda Fernandez, Pamela Araya, Rosanna Lagos, Juan Carlos Hormazabal **Analyzed the data** Mailis Maes, Megan E. Carey, Ellen E. Higginson, **Contributed reagents/materials/analysis tools** Zoe A. Dyson, Stephen Baker, Michael J. Sikorski, Sharon M. Tennant, **Funding acquisition**, Myron M. Levine, Gordon Dougan **Writing original draft** Mailis Maes, Michael J. Sikorski, Megan E. Carey,, Ellen E. Higginson **Writing review and editing** Mailis Maes, Michael J. Sikorski, Megan E. Carey, Ellen E. Higginson,, Myron M. Levine, Zoe A. Dyson, Sharon M. Tennant ## Funding This work was supported by a grant [OPP1161058] from the Bill & Melinda Gates Foundation, this work was also supported by the Gates Foundation (TyVAC) [OPP1151153] and the Wellcome Trust (STRATAA) [106158 and 098051]. MM was funded by the NIHR Cambridge Biomedical Research Centre and NIHR AMR Research Capital Funding Scheme [NIHR200640]. ZAD was supported by a grant funded by the Wellcome Trust (STRATAA; 106158/Z/14/Z), and received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement TyphiNET (#845681). MJS received research support from the NIH NIAID grants F30AI156973 and U19AI110820 as well as NIDDK T32DK067872. The views expressed are those of the authors and not necessarily those of the NIHR or the Department of Health and Social Care. The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. ## Data availability Raw sequence data have been submitted to the European Nucleotide Archive (ENA) and accession numbers of individual isolates are listed in the appendix. ## Competing interest The authors have declared that no competing interests exist ## Acknowledgements We express our thanks to the typhoid fever patients whose samples were analysed and included in this paper, as well as the local public health personnel in Chile. Furthermore, we wish to thank Nicholas Thomson for access to the Sanger analysis pipelines. * Received January 19, 2022. * Revision received January 19, 2022. * Accepted January 21, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/) ## References 1. 1.Marco C, Delgado I, Vargas C, Munoz X, Bhutta ZA, Ferreccio C. Typhoid Fever in Chile 1969-2012: Analysis of an Epidemic and Its Control. Am J Trop Med Hyg. 2018;99(3_Suppl):26–33. doi:10.4269/ajtmh.18-0125 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4269/ajtmh.18-0125&link_type=DOI) 2. 2.Levine MM BR, Ferreccio C, Clements ML, Lanata C, Sear S, Morris JG, Cisneros L, Germanier R. Interventions to control endemic typhoid fever: field studies in Santiago, Chile. Control and Eradication of Infectious Diseases, An International Symposium. 1985;Series 1:37-53. 3. 3.Black RE, Levine MM, Ferreccio C, Clements ML, Lanata C, Rooney J, et al. Efficacy of one or two doses of Ty21a Salmonella typhi vaccine in enteric-coated capsules in a controlled field trial. Vaccine. 1990;8(1):81–4. doi:10.1016/0264-410X(90)90183-M [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0264-410X(90)90183-M&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2180234&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990CM94100016&link_type=ISI) 4. 4.Levine MM, Ferreccio C, Ortiz E, Cryz S. Comparison of enteric-coated capsules and liquid formulation of Ty21a typhoid vaccine in randomised controlled field trial. The Lancet. 1990;336(8720):891–4. doi;10.1016/0140-6736(90)92266-K [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0140-6736(90)92266-K&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1976928&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1990EC45600001&link_type=ISI) 5. 5.Levine MM, Ferreccio C, Black RE, Germanier R. Large-scale field trial of Ty21a live oral typhoid vaccine in enteric-coated capsule formulation. Lancet. 1987;1(8541):1049–52. doi:10.1016/s0140-6736(87)90480-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(87)90480-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2883393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1987H180100001&link_type=ISI) 6. 6.Ferreccio C, Levine MM, Rodriguez H, Contreras R. Comparative efficacy of two, three, or four doses of TY21a live oral typhoid vaccine in enteric-coated capsules: a field trial in an endemic area. J Infect Dis. 1989;159(4):766–9. doi:10.1093/infdis/159.4.766 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/159.4.766&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=2647863&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1989U006300026&link_type=ISI) 7. 7.Ferreccio C, Morris JG, Jr.., Valdivieso C, Prenzel I, Sotomayor V, Drusano GL, et al. Efficacy of ciprofloxacin in the treatment of chronic typhoid carriers. J Infect Dis. 1988;157(6):1235–9. doi:10.1093/infdis/157.6.123 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/157.6.1235&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3373023&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1988N738500017&link_type=ISI) 8. 8.Gauld JS, Hu H, Klein DJ, Levine MM. Typhoid fever in Santiago, Chile: Insights from a mathematical model utilizing venerable archived data from a successful disease control program. PLoS Negl Trop Dis. 2018;12(9):e0006759. doi:10.1371/journal.pntd.0006759 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0006759&link_type=DOI) 9. 9.Sears SD, Ferreccio C, Levine MM, Cordano AM, Monreal J, Black RE, et al. The use of Moore swabs for isolation of Salmonella typhi from irrigation water in Santiago, Chile. J Infect Dis. 1984;149(4):640–2. doi:10.1093/infdis/149.4.640 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/149.4.640&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=6373964&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 10. 10.Levine M. South America: The return of cholera. The Lancet. 1991;338(8758):45–6. doi:10.1016/0140-6736(91)90024-J [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0140-6736(91)90024-J&link_type=DOI) 11. 11.Medina E. [Epidemic of cholera in Chile; 1991]. Rev Med Chil. 1991;119(8):943–56. [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1844781&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 12. 12.Swerdlow DL, Mintz ED, Rodriguez M, Tejada E, Ocampo C, Espejo L, et al. Waterborne transmission of epidemic cholera in Trujillo, Peru: lessons for a continent at risk. Lancet. 1992;340(8810):28–33. doi:10.1016/0140-6736(92)92432-f [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/0140-6736(92)92432-F&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1351608&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992JC35400012&link_type=ISI) 13. 13.Tauxe RV, Blake PA. Epidemic cholera in Latin America. Jama. 1992;267(10):1388–90. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jama.1992.03480100098039&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1740864&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992HG67800036&link_type=ISI) 14. 14.Wong VK, Baker S, Connor TR, Pickard D, Page AJ, Dave J, et al. An extended genotyping framework for Salmonella enterica serovar Typhi, the cause of human typhoid. Nature communications. 2016;7(1):1–11. doi: 10.1038/ncomms12827 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms12827&link_type=DOI) 15. 15.Inouye M, Dashnow H, Raven LA, Schultz MB, Pope BJ, Tomita T, et al. SRST2: Rapid genomic surveillance for public health and hospital microbiology labs. Genome Med. 2014;6(11):90. doi:10.1186/s13073-014-0090-6 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1186/s13073-014-0090-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25422674&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 16. 16.Dyson ZA, Holt KE. Five years of GenoTyphi: updates to the global Salmonella Typhi genotyping framework. The Journal of Infectious Diseases. 2021. doi:10.1093/infdis/jiab414 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/infdis/jiab414&link_type=DOI) 17. 17.Argimon S, Yeats CA, Goater RJ, Abudahab K, Taylor B, Underwood A, et al. A global resource for genomic predictions of antimicrobial resistance and surveillance of Salmonella Typhi at pathogenwatch. Nat Commun. 2021;12(1):2879. doi:10.1038/s41467-021-23091-2 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-23091-2&link_type=DOI) 18. 18.Page AJ, De Silva N, Hunt M, Quail MA, Parkhill J, Harris SR, et al. Robust high-throughput prokaryote de novo assembly and improvement pipeline for Illumina data. Microb Genom. 2016;2(8):e000083. doi:10.1099/mgen.0.000083 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/mgen.0.000083&link_type=DOI) 19. 19.Parkhill J, Dougan G, James K, Thomson N, Pickard D, Wain J, et al. Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature. 2001;413(6858):848–52. doi:10.1038/35101607 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/35101607&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11677608&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171750200047&link_type=ISI) 20. 20.Kariuki S, Mbae C, Van Puyvelde S, Onsare R, Kavai S, Wairimu C, et al. High relatedness of invasive multi-drug resistant non-typhoidal Salmonella genotypes among patients and asymptomatic carriers in endemic informal settlements in Kenya. PLoS Negl Trop Dis. 2020;14(8):e0008440. doi:10.1371/journal.pntd.0008440 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0008440&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 21. 21.Britto CD, Dyson ZA, Duchene S, Carter MJ, Gurung M, Kelly DF, et al. Laboratory and molecular surveillance of paediatric typhoidal Salmonella in Nepal: Antimicrobial resistance and implications for vaccine policy. PLoS neglected tropical diseases. 2018;12(4):e0006408. doi: 10.1371/journal.pntd.0006408 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0006408&link_type=DOI) 22. 22.Rahman SIA, Dyson ZA, Klemm EJ, Khanam F, Holt KE, Chowdhury EK, et al. Population structure and antimicrobial resistance patterns of Salmonella Typhi isolates in urban Dhaka, Bangladesh from 2004 to 2016. PLoS neglected tropical diseases. 2020;14(2):e0008036. doi: 10.1371/journal.pntd.0008036 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0008036&link_type=DOI) 23. 23.Wong VK, Holt KE, Okoro C, Baker S, Pickard DJ, Marks F, et al. Molecular Surveillance Identifies Multiple Transmissions of Typhoid in West Africa. PLoS Negl Trop Dis. 2016;10(9):e0004781. doi:10.1371/journal.pntd.0004781 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0004781&link_type=DOI) 24. 24.Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, Bentley SD, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43(3):e15. doi:10.1093/nar/gku1196 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gku1196&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25414349&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 25. 25.Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90. doi:10.1093/bioinformatics/btl446 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btl446&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16928733&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000241629600016&link_type=ISI) 26. 26.Letunic I, Bork P. Interactive Tree Of Life (iTOL): an online tool for phylogenetic tree display and annotation. Bioinformatics. 2007;23(1):127–8. doi:10.1093/bioinformatics/btl529 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btl529&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17050570&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000243060300021&link_type=ISI) 27. 27.Argimón S, Abudahab K, Goater RJE, Fedosejev A, Bhai J, Glasner C, et al. Microreact: visualizing and sharing data for genomic epidemiology and phylogeography. Microb Genom. 2016;2(11):e000093. doi:10.1099/mgen.0.000093 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/mgen.0.000093&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28348833&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 28. 28.Gupta SK, Padmanabhan BR, Diene SM, Lopez-Rojas R, Kempf M, Landraud L, et al. ARG-ANNOT, a new bioinformatic tool to discover antibiotic resistance genes in bacterial genomes. Antimicrob Agents Chemother. 2014;58(1):212–20. doi:10.1128/aac.01310-13 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjg6IjU4LzEvMjEyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDEvMjQvMjAyMi4wMS4xOS4yMjI2OTU3Ny5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 29. 29.Carattoli A, Zankari E, Garcia-Fernandez A, Voldby Larsen M, Lund O, Villa L, et al. In silico detection and typing of plasmids using PlasmidFinder and plasmid multilocus sequence typing. Antimicrob Agents Chemother. 2014;58(7):3895–903. doi:10.1128/AAC.02412-14 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYWFjIjtzOjU6InJlc2lkIjtzOjk6IjU4LzcvMzg5NSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAxLzI0LzIwMjIuMDEuMTkuMjIyNjk1NzcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 30. 30.Chen L, Yang J, Yu J, Yao Z, Sun L, Shen Y, et al. VFDB: a reference database for bacterial virulence factors. Nucleic Acids Res. 2005;33(Database issue):D325–8. doi:10.1093/nar/gki008 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gki177&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15608208&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000226524300067&link_type=ISI) 31. 31.Wick RR, Judd LM, Gorrie CL, Holt KE. Unicycler: Resolving bacterial genome assemblies from short and long sequencing reads. PLoS Comput Biol. 2017;13(6):e1005595. doi:10.1371/journal.pcbi.100559 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pcbi.1005595&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28594827&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 32. 32.Wick RR, Schultz MB, Zobel J, Holt KE. Bandage: interactive visualization of de novo genome assemblies. Bioinformatics. 2015;31(20):3350–2. doi:10.1093/bioinformatics/btv383 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv383&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26099265&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 33. 33.Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, Parkhill J. ACT: the Artemis Comparison Tool. Bioinformatics. 2005;21(16):3422–3. doi:10.1093/bioinformatics/bti553 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/bti553&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15976072&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000231360600013&link_type=ISI) 34. 34.Sullivan MJ, Petty NK, Beatson SA. Easyfig: a genome comparison visualizer. Bioinformatics. 2011;27(7):1009–10. doi:10.1093/bioinformatics/btr0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btr039&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21278367&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000289162000018&link_type=ISI) 35. 35.Maes M, Dyson ZA, Higginson EE, Fernandez A, Araya P, Tennant SM, et al. Multiple Introductions of Salmonella enterica Serovar Typhi H58 with Reduced Fluoroquinolone Susceptibility into Chile. Emerg Infect Dis. 2020;26(11):2736–40. doi:10.3201/eid2611.201676 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3201/eid2611.201676&link_type=DOI) 36. 36.Dyson ZA, Klemm EJ, Palmer S, Dougan G. Antibiotic resistance and typhoid. Clinical Infectious Diseases. 2019;68(Supplement_2):S165–S70. doi: 10.1093/cid/ciy111 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciy1111&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 37. 37.Dyson ZA, Thanh DP, Bodhidatta L, Mason CJ, Srijan A, Rabaa MA, et al. Whole genome sequence analysis of Salmonella Typhi isolated in Thailand before and after the introduction of a national immunization program. PLoS neglected tropical diseases. 2017;11(1):e0005274. doi: 10.1371/journal.pntd.0005274 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0005274&link_type=DOI) 38. 38.Fica AE, Prat-Miranda S, Fernandez-Ricci A, D’Ottone K, Cabello FC. Epidemic typhoid in Chile: analysis by molecular and conventional methods of Salmonella typhi strain diversity in epidemic (1977 and 1981) and nonepidemic (1990) years. J Clin Microbiol. 1996;34(7):1701–7. doi:10.1128/jcm.34.7.1701-1707.1996 [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNtIjtzOjU6InJlc2lkIjtzOjk6IjM0LzcvMTcwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzAxLzI0LzIwMjIuMDEuMTkuMjIyNjk1NzcuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 39. 39.Ames WR, Robins M. Age and Sex as Factors in the Development of the Typhoid Carrier State, and a Method for Estimating Carrier Prevalence. Am J Public Health Nations Health. 1943;33(3):221–30. doi:10.2105/ajph.33.3.221 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2105/AJPH.33.3.221&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18015749&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 40. 40.Chattaway MA, Gentle A, Nair S, Tingley L, Day M, Mohamed I, et al. Phylogenomics and antimicrobial resistance of Salmonella Typhi and Paratyphi A, B and C in England, 2016-2019. Microb Genom. 2021;7(8). doi:10.1099/mgen.0.000633 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/mgen.0.000633&link_type=DOI) 41. 41.Ingle DJ, Nair S, Hartman H, Ashton PM, Dyson ZA, Day M, et al. Informal genomic surveillance of regional distribution of Salmonella Typhi genotypes and antimicrobial resistance via returning travellers. PLoS Negl Trop Dis. 2019;13(9):e0007620. doi:10.1371/journal.pntd.0007620 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0007620&link_type=DOI) 42. 42.Guevara PD, Maes M, Thanh DP, Duarte C, Rodriguez EC, Montaño LA, et al. A genomic snapshot of Salmonella enterica serovar Typhi in Colombia. PLoS Negl Trop Dis. 2021;15(9):e0009755. doi:10.1371/journal.pntd.0009755 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0009755&link_type=DOI) 43. 43.Wong VK, Baker S, Pickard DJ, Parkhill J, Page AJ, Feasey NA, et al. Phylogeographical analysis of the dominant multidrug-resistant H58 clade of Salmonella Typhi identifies inter-and intracontinental transmission events. Nature genetics. 2015;47(6):632–9. doi: 10.1038/ng.3281 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3281&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25961941&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 44. 44.Dyson ZA, Malau E, Horwood PF, Ford R, Siba V, Yoannes M, et al. Whole genome sequence analysis of Salmonella Typhi in Papua New Guinea reveals an established population of genotype 2.1.7 sensitive to antimicrobials. medRxiv. 2021:2021.09.27.21264209. 45. 45.Kariuki S, Dyson ZA, Mbae C, Ngetich R, Kavai SM, Wairimu C, et al. Multiple introductions of multidrug-resistant typhoid associated with acute infection and asymptomatic carriage, Kenya. Elife. 2021;10. doi:10.7554/eLife.67852 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.67852&link_type=DOI) 46. 46.Thanh Duy P, Thieu NTV, Nguyen Thi Nguyen T, Ngoc Dan Thanh H, Dongol S, Karkey A, et al. Gallbladder carriage generates genetic variation and genome degradation in Salmonella Typhi. PLoS Pathog. 2020;16(10):e1008998. doi:10.1371/journal.ppat.1008998 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1008998&link_type=DOI) 47. 47.Johnson R, Mylona E, Frankel G. Typhoidal Salmonella: Distinctive virulence factors and pathogenesis. Cell Microbiol. 2018;20(9):e12939. doi:10.1111/cmi.12939 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/cmi.12939&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30030897&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 48. 48.Murret-Labarthe C, Kerhoas M, Dufresne K, Daigle F. New Roles for Two-Component System Response Regulators of Salmonella enterica Serovar Typhi during Host Cell Interactions. Microorganisms. 2020;8(5). doi:10.3390/microorganisms8050722 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/microorganisms8050722&link_type=DOI) 49. 49.Forest CG, Ferraro E, Sabbagh SC, Daigle F. Intracellular survival of Salmonella enterica serovar Typhi in human macrophages is independent of Salmonella pathogenicity island (SPI)-2. Microbiology (Reading). 2010;156(Pt 12):3689–98. doi:10.1099/mic.0.041624-0 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1099/mic.0.041624-0&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20817644&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000285806200018&link_type=ISI) 50. 50.Karlinsey JE, Stepien TA, Mayho M, Singletary LA, Bingham-Ramos LK, Brehm MA, et al. Genome-wide Analysis of Salmonella enterica serovar Typhi in Humanized Mice Reveals Key Virulence Features. Cell Host Microbe. 2019;26(3):426-34.e6. doi:10.1016/j.chom.2019.08.001 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.chom.2019.08.001&link_type=DOI) 51. 51.Desai SK, Winardhi RS, Periasamy S, Dykas MM, Jie Y, Kenney LJ. The horizontally-acquired response regulator SsrB drives a Salmonella lifestyle switch by relieving biofilm silencing. Elife. 2016;5. doi:10.7554/eLife.10747 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.10747&link_type=DOI) 52. 52.Marshall JM, Flechtner AD, La Perle KM, Gunn JS. Visualization of extracellular matrix components within sectioned Salmonella biofilms on the surface of human gallstones. PLoS One. 2014;9(2):e89243. doi:10.1371/journal.pone.0089 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0089243&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24551241&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) 53. 53.Parkhill J, Dougan G, James KD, Thomson NR, Pickard D, Wain J, et al. Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature. 2001;413(6858):848–52. doi:10.1038/35101607 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/35101607&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11677608&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171750200047&link_type=ISI) 54. 54.Kidgell C, Pickard D, Wain J, James K, Diem Nga LT, Diep TS, et al. Characterisation and distribution of a cryptic Salmonella typhi plasmid pHCM2. Plasmid. 2002;47(3):159–71. doi:10.1016/s0147-619x(02)00013 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0147-619X(02)00013-6&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12151231&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F01%2F24%2F2022.01.19.22269577.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000177455300001&link_type=ISI) 55. 55.Dyson ZA, Thanh DP, Bodhidatta L, Mason CJ, Srijan A, Rabaa MA, et al. Whole Genome Sequence Analysis of Salmonella Typhi Isolated in Thailand before and after the Introduction of a National Immunization Program. PLoS Negl Trop Dis. 2017;11(1):e0005274. doi:10.1371/journal.pntd.0005 [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pntd.0005274&link_type=DOI)