BACKGROUND
typhoid fever is a human systemic infection that is caused by salmonella enterica subsp.
enterica serovar typhi .
this human-restricted and highly adapted pathogen is transmitted via the oral-faecal route.
s. typhi is responsible for  <dig>  million infections and results in approximately  <dig>  deaths worldwide annually  <cit> .
the disease primarily causes acute systemic infection with life-threatening complications, and the recovering patient may develop into a chronic carrier state  <cit> .
typhoid is endemic, with periodic outbreaks and sporadic cases occurring in developing countries, particularly in southeast asia, south central asia, latin america and southern africa where sanitary conditions are poor  <cit> .
among the  <dig> states of malaysia, kelantan has a significantly higher incidence of typhoid fever  <cit> .
a large typhoid outbreak occurred in kelantan state, which resulted in  <dig> cases of infection and two deaths in a short period of 3 months   <cit> .
previously, pulsed-field gel electrophoresis  revealed close genetic relatedness between a s. typhi strain that was isolated from an asymptomatic carrier in  <dig> and a strain that originated from a patient during the  <dig> outbreak in kelantan .
human carriers are the main reservoirs of s. typhi transmission, but the genetic basis, and the underlying mechanisms in particular, are unclear  <cit> .
it has been suggested that a carrier strain will likely lack gene acquisition capabilities and have little fitness advantages compared with those strains causing symptomatic infections because the human reservoir is small and physiologically isolated  <cit> .
therefore, it is of great interest to know as to what extent these closely related strains differed or shared in its genomic contents despite being isolated from these two distinct epidemiological settings.
in  <dig>  a s. typhi strain was isolated from a sporadic case in kuala lumpur.
this genome has been sequenced earlier  <cit> , and pfge analyses showed that this strain is more distantly related to the outbreak and carrier strains, but the epidemiological link is unknown .
recent whole-genome sequencing of s. typhi has demonstrated that the pathogen shows limited genetic variation with little evidence of purifying selection, antigenic variation or recombination between isolates  <cit> .
this clonal pathogen, however, is associated with varying degrees of disease severity in different regions  <cit> .
previous pfge studies have also demonstrated genome size variations and distinct pfge patterns in relation with fatal and non-fatal typhoid cases  <cit> .
although the health conditions of the host cannot be completely ruled out, various reports have suggested that the gain and loss of genes through mutations and gene transfers that have occurred independently in different lineages have markedly contributed to the varying pathogenic potentials  <cit> .
however, these important factors are poorly understood because there is limited genomic information for s. typhi, particularly involving strains that are associated with diverse epidemiological settings.
its genomic heterogeneity is likely due to the adaptation of the pathogen to the host and its exposure to mobile elements, such as bacteriophages  <cit> .
organisms having common core genomes could differ in their dispensable  genes, reflecting their unique physiological and virulence properties  <cit> .
although not all genetic variations are essential for adaptation, some dispensable genes are believed to be responsible for conferring fitness advantages to the pathogen to thrive in its host.
horizontal gene transfer is also thought to be the predominant force in bacterial evolution, which contributes to novel gene acquisition.
the acquired genes provide new characteristics, which either aid in host adaptation and persistence or enhance virulence capabilities .
a previous study on the pan-genome of salmonella enterica revealed that the pan-genome  of all strains will continue to increase as new genomes are sequenced  <cit> .
with the availability of robust next-generation sequencing technologies, high quality whole genome sequences can be generated and analysed, which will be especially useful for capturing fine variations among highly conserved s. typhi strains.
multiple whole genome sequence comparisons of closely related strains will not only lead to the better understanding of their relationships but also provide novel insights into the functional roles of strain-specific genes.
in this study, we performed detailed and comprehensive comparative functional analyses of three previously sequenced genomes of s. typhi strains that were isolated from typhoid patients during a large outbreak in  <dig>  a sporadic case in  <dig> and an asymptomatic carrier from malaysia, where typhoid is endemic.
these malaysian s. typhi strains were compared with previously published s. typhi genomes with the following aims: 1) to determine and describe the genomes signatures and conserved and unique regions of the strains that were being studied; 2) to elucidate the phylogeny and genetic relatedness of these malaysian strains compared with  <dig> other published strains using phylogenomic analysis; 3) to compare those strains that have been associated with various epidemiological settings , and particularly the regions of plasticity that may contribute to the varying pathogenic potentials; 4) to identify potential novel pathogenic factors that are harboured by the analysed strains; 5) to provide insight into the possible differential functionalities of the genes, focusing mainly on virulence- and persistence-related genes based on non-synonymous snps of closely related strains and particularly on carrier strains to gain insight into the persistence of the carrier state; and 6) to understand how the potential nssnps affect protein structures and functions.
the data that are generated will be useful for the profiling of strains, marker development and the increased understanding of outbreak, sporadic and less studied asymptomatic typhoid carriage infection.
RESULTS
general genome signatures of s. typhi in association with outbreak, sporadic case and carrier
the genomes of the three previously sequenced malaysian s. typhi strains from different sources were compared to identify potential genomic features that may help to elucidate the different disease outcomes.
one strain was isolated during the largest outbreak in the country, one from a carrier  during typhoid surveillance following the outbreak, representing the carrier state, and one from a sporadic case in the metropolitan area of kuala lumpur, which has a relatively lower incidence of typhoid.
however, because of the limited numbers of strains studied and absence of information regarding the pan-genome of the s.typhi population, the genetic differences observed should be taken with caution.
the aim of the detailed comparative analysis was to provide a better understanding and insights into the unusual genome dynamics of s. typhi, a highly clonal organism.
previous genome sequencing analyses have generated high-quality assemblies with an average genome coverage of 100× for the  <dig> malaysian s. typhi genomes, including bl <dig>   <cit> , cr <dig>   <cit>  and st <dig>   <cit> .
the approximate predicted genome sizes and average guanine-plus-cytosine  contents of all  <dig> genomes ranged from  <dig>  mb to  <dig>  mb and  <dig> % to  <dig> %, respectively .
these genomes form a single and circular chromosome with no plasmids detected.
an in silico multi-locus sequence typing  analysis classified both bl <dig> and cr <dig> as st  <dig> and st <dig> as st  <dig>  which are the main sequence types that have been associated with the worldwide distribution of s. typhi out of the  <dig> sts that have been identified to date  <cit> .
the predicted coding sequences  of the genomes based on rast subsystem-based annotations varied from  <dig>  to  <dig>  with an average coding percentage of  <dig> %.
the average sizes of the cdss were similar , indicating that the size differences among the genomes are largely attributable to a number of cdss and intergenic regions.
approximately 12% of the cdss were annotated as uncharacterised proteins .
some of these genes  were observed to vary from one strain to another.
these “dispensable” genomes carried genes that were present in one or more strains and could even be unique to a single strain  <cit> , indicating a possible open pan-genome of s. typhi.
in general, the chromosomes of the three assembled genomes exhibited overall structural conservation and colinearity with each other as evidenced by the homologous and conserved regions that were shared and the very small strain-specific regions , which may harbour genes that are relevant to the specific adaptations and fitness advantages of each of the strains.
these regions most likely represent dna that was acquired during events of hgt that may provide the strains with greater metabolic versatility or even virulence capabilities, as will be further discussed in the plasticity section.figure  <dig>
circular genomes representation map and genome comparison of
s
.
typhi .
the circle is divided into arcs representing the sequence of the chromosomes of all six genomes as labeled.
the inner ring shows the coordinates in scale and total genomes size of the reference sequence, ct <dig> in scale , with black histogram bar representing gc content and purple-green histogram bar representing gc deviation.
orthologs for each genomes with respect to ct <dig>  are showed in order  with the percentage of similarity based on nucleotide sequences  as indicated on the legend located at the right of the figure.
the outermost arch in red represents the location of spi <dig> to spi <dig> in relative to ct <dig> and labeled in red.
the position of the genes related to pathogenicity and host persistent are shown and labeled in blue at the edge of the rings.
comparative genomics of s. typhi
genomic comparisons were performed on the three malaysian s. typhi strains and the three completed s. typhi genomes , ct <dig>   <cit> , ty <dig>   <cit>  and p-stx- <dig>   <cit> .
the comparisons of the studied strains with the reference genomes allowed for the elucidation of novel and additional genes that are carried by the malaysian s. typhi genomes that may be of significance.
we have analysed the shared and unique genes of all six genomes to determine their distinct virulence and pathogenic features.
as expected, the six genomes exhibited high similarities and syntenies with each other with limited evidence of genomic rearrangements, which collectively indicate stable genomic structures.
the majority of the orfs from the compared genomes were part of a conserved genomic core, in which  <dig> orfs were shared among all of the genomes.
these shared orfs provide clear evidence of conservation among the genomes of the s. typhi strains.
the rest of the unshared orfs or accessory genes are present in one or more strains, which represent the salient differences in the genomes.
most of these orfs  that were harboured by each genome were annotated as hypothetical proteins .
our extended homology analysis has shown that the remaining unshared orfs  were likely to encode for functional proteins from diverse categories, including virulence-related proteins, secretory proteins, conserved domain proteins, transporter proteins and phage proteins among others.
these considerable portions of the genomes could provide important functional clues for understanding the virulence and persistence of the pathogen more clearly, anticipating the need for extensive future studies focusing on their possible roles in bacterial pathogenesis.
however, the numbers of shared and unshared orfs may have been underestimated because the genomes were incomplete.
among the shared genes that were found between strains bl <dig> and cr <dig>  uncommon orfs that encoded for the vi icm-f secretion protein, icm-f-related protein and type vi secretion protein evpb were identified whose products are related to the type vi secretion system.
the genes shared 99% similarity with the type vi secretion protein of salmonella typhimurium strain d <dig>  <cit>  and were only found in bl <dig> and cr <dig>  this protein was recently recognised as one of the main virulence determinants in burkholderia pseudomallei, legionella pneumophila and vibrio cholerae, but its function in s. typhi remains to be elucidated  <cit> .
t6ss genes are believed to be involved in either structural components of the secretory apparatus, secretory products or assisting with protein translocation; for example, providing the energy to push substrates through the channel of the apparatus  <cit> .
these genes are also proposed to be involved in surface reorganisation, enhancing adherence to epithelial cells, intracellular multiplication and human macrophage killing .
other t6ss clusters were found intact as in reference genomes.
the high similarities of the genetic contents of bl <dig> and cr <dig> with minor variations, and particularly the presence of unique accessory genes , are in agreement with the pfge pulsotype data, which revealed that both strains are genetically similar, showing a difference of only one band .
the chromosome of salmonella enterica is commonly integrated with a large portion of horizontally acquired dna apart from its core, which are termed the salmonella pathogenicity islands   <cit> .
these acquired spis have led to divergence and host restriction similar to those in s. typhi.
the identification of conserved spis and their variations have important implications in a wide range of microbiological applications, such as antigen and marker discovery and the identification of essential genes and their respective traits.
in this study, we have annotated all spis  and its genetic variant of s. typhi in the genome, which is characterised by its deviated gc content, flanking by trna genes and the presence of phages, integrases, recombinases and genes that are related to dna integration.
although all the spis   were detected in the genomes, there were marked variations.
the presence of a large number of transposition-related genes in these spis suggests that the sites may be actively involved in the integration and transposition of genetic elements, which drive genetic variation.
interestingly, our comparative analysis revealed that the carrier strain p-stx- <dig> lacks a ~10 kb prpz cluster and adjacent gene clusters harbouring  <dig> orf with a deviated gc skew of  <dig>  % at spi- <dig> but remains fully intact in our carrier strain .
previously, a prpz cluster deletion study showed that the mutant has a significantly lower survival rate compared with the parental strain, which may be due to a signalling pathway that controls the long-term survival of s. typhi in host cells, and particularly, the survival in human macrophages  <cit> .
in fact, our results support this study with the fine-tuned postulation that the deletion led to reduced virulence that enabled the carrier strain to coexist with the host; for example, in the tissue of gall bladder.
this possibly explains why long-term survival in the macrophage is no longer necessary, which is presumably because the pathogens have colonised and persisted in other cells of the host during adaptation.
however, the deletion was not detected in our carrier strain, suggesting that the genes may not be the only factors that are relevant to a carrier state.
furthermore, the region is flanked by multiple transposases, integrases, ligases and uncharacterised proteins, which are known to be involved in transposition.
additionally, genes coding for the dna mismatch repair protein mutc and transposase were identified both upstream and downstream of the cluster, suggesting that the region could have possibly been acquired earlier during horizontal gene transfer.
the presence of a dna mismatch protein gene has been previously implicated to be involved in modulating recombination events by incorporating or inhibiting the transfer of mobile genetic elements  <cit> .
the deletion of the gene cluster and the presence of a large number of genes that are related to transposition indicate that spi- <dig> may be unstable and prone to excision similar to the precise excision of the crucial spi- <dig>  which has been recently reported in s. typhi  <cit> , suggesting that gene deletion may be important to the host adaptation of this organism, although independent acquisition or gene gain by other strains cannot be completely ruled out.
these findings expand upon previous studies, which have reported that other spis are relatively stable in the genome  <cit> , highlighting the importance of a future evaluation of the stability of the other spis.
apart from these remarkable differences, the comparison of the two carrier strains, cr <dig> and p-stx- <dig>  revealed that cr <dig> carries several additional genes that were not identified in p-stx- <dig> that encode for unknown functions and phages that are present in the phage region.
almost all of the potential major virulence and persistence-associated genes have homologues in p-stx- <dig>  suggesting that they are not specifically associated with the unique persistence of the carrier strains but are common to s. typhi.
the genomic structure of the sporadic strain st <dig> is more conserved in comparison and has relatively fewer dispensable orfs, which are mainly genes that code for hypothetical proteins and phages, indicating the conservation of large numbers of genes, which is essential for strict host adaptation and virulence optimisation.
phylogenomics of s. typhi revealed shared common ancestry
we determined a core genome-based phylogeny by mapping  <dig> query genomic sequences against ct <dig> into a single non-redundant alignment of  <dig> ,681 bp  .
the phylogenomic tree showed that the outbreak strain bl <dig> and carrier strain cr <dig> were closely related and could be differentiated by only  <dig> snps.
these data are in agreement with our pfge results that showed that both strains are highly related .
the observed close phylogenetic relationship between these two strains is consistent with our earlier speculation that the large outbreak that occurred in  <dig> shared a common ancestor strain with the typhoid carrier that may have been circulating for a long period in the country.
it is challenging to determine how these two strains are related, considering their short-term evolutionary relationship.
however, three evolutionary postulations are possible.
first, the carrier strain may have been derived from the outbreak strain.
second, the carrier strain may have long existed in a carrier who served as a reservoir and the source of the outbreak.
finally, both of these strains may have diverged independently from a common ancestor, which was possibly harboured in a long existed carrier to give rise to two independent cases, considering the geographical proximities.
notably, these two highly related strains with unique gene repertoires clustered with  <dig> epidemiologically and geographically unrelated strains from india , russia , vietnam  and senegal in west africa  .
interestingly, all of these strains, including bl <dig>  cr <dig>  p-stx- <dig>  ty <dig> and e01- <dig>  were subtyped as st  <dig> .
conversely, the malaysian sporadic s. typhi strain st <dig> clustered closely with the multidrug-resistant strain ct <dig> from vietnam together with other geographically related strains from indonesia  and bangladesh  , which were all subtyped as st  <dig>  suggesting the possible movement of clonally related strains among the southeast asian countries.
such close genetic relatedness that is based on macrorestriction and snp typing analyses has been previously reported  <cit> .
the rest of the genomes were clustered as st  <dig>  from the analysis, we found no temporal or geographical signals, but the sequence types was highly correlated with the phylogenomic clustering.
apparently, the cluster subtyped as st  <dig> could be further differentiated into two clusters according to the phylogenomic analysis but was limited in the current mlst scheme for s. typhi.
these results indicate that the phylogenomic analysis has much better resolution power compared with mlst in separating highly clonal strains in s. typhi.
this information is essential for devising a better set of mlst alleles for improved molecular typing.
these data further support the widespread distribution of st  <dig> and st  <dig> as the major genotypes that occur worldwide, although rarely, st  <dig> and st  <dig> have also been isolated in a previous study  <cit> .
it is important to note that the two carriers, cr <dig> and p-stx- <dig>  belonged to two different subtypes, reflect the non-universality of the genotypes relevant for carrier state transformation.figure  <dig>
phylogenomic tree inferred by approximately-maximum-likelihood method from the aligned core genomes.
multiple genomes alignments were generated by mapping genome sequences of the  <dig> global s. typhi strains against ct <dig> at all sites relevant for phylogenomic analysis using realphy  <cit> .
phylogenetic tree  was inferred via approximately-maximum-likelihood method using fasttreemp  <cit> .
strains studied were labelled in blue.
bootstrap support values shown at each node.
genome plasticity
our analyses of the relationships among the strains and their genomic variations are further supported and extended to the study of the genome plasticity of s. typhi.
in many organisms, genomic plasticity is commonly observed, but s. typhi has generally demonstrated few variations compared with other salmonella spp.
<cit> .
however, in our study, we detected considerable genetic variations, predominantly involving is elements, phages and crisprs.
marked variations in the numbers and types of is elements were detected.
is <dig>  and is <dig>  were both abundant in the three malaysian s. typhi genomes.
the is elements, such as is <dig>  have been widely used as molecular markers for subtyping due to their genome-wide distributions and high levels of diversity, but their roles in modulating the gene expression in s. typhi have yet to be clarified  <cit> .
recent studies on enterohaemorrhagic escherichia coli o <dig> have shown that the presence of is could play a role in the gene inactivation and immobilisation of incoming phages and plasmids, leading to the diversification and evolution of the bacterial genome  <cit> .
is elements have also been shown to affect the expression of neighbouring genes and induce genomic rearrangements   <cit> .
however, little is understood about their roles in modulating virulence and gene expression in s. typhi.
the variations that were detected may provide clues on how these differences affect the virulence and fitness strategies of the pathogens.
the s. typhi strains bl <dig>  cr <dig> and st <dig> carry eight, seven and eight phages, respectively.
we have identified substantial phage variations among the s. typhi genomes.
one of the differentiating features was a distinct set of prophages that were harboured by both bl <dig> and cr <dig>  which rendered them less unique compared with st <dig> with the exception of a few orfs that encoded for phages and hypothetical proteins.
interestingly, the phages that were carried by st <dig> had relatively shorter in lengths  compared with the phages that were identified in the other strains.
as expected, both of the closely related strains  had highly similar phage contents and carried an additional intact salmonella phage re- <dig> with uncommon orfs, which mainly contained genes coding for hypothetical proteins, phage proteins, prophage-like proteins, repressor proteins, excisionases, terminases and integrases.
the variations that were detected in the numbers of predicted prophages and prophage-like regions illustrated the dynamics of phage gain and loss that distinguished one strain from the others, indicating that phages may play important roles in genomic diversity.
by correlating phylogenomic analysis and whole genome sequence alignments, this region appears to show the typical gain and loss of sequences during the course of genomic evolution.
this evolutionary relationship is consistent with the phage variations, providing a useful framework for investigating the relationships of strains and their respective phenotypes.
the phage was likely acquired prior to the divergence of the common ancestor of both bl <dig> and cr <dig>  through horizontal gene transfer rather than phage loss.
alternatively, due to the advantageous roles of the hgt events, the phage proteins that were acquired by the strains could promote their in vivo survival and pathogenesis  <cit> .
among the detected phages, the two typical s. typhi phages, gifsy- <dig> and fels- <dig>  were both found to be intact and conserved in all six of the s. typhi genomes, suggesting that these regions may play essential roles and provide fitness advantages to the pathogen.
apart from the intact phages,  <dig> incomplete phages, including burkholderia_phage_bcepmu, enterobacteria_phage_cdti and lactococcus_phage_bil <dig>  were also observed in the six genomes, suggesting that these common phages were likely present in their common ancestors and may be relevant for host adaptation and survival.
clustered regularly interspaced short palindromic repeats  were recently found to be important in bacteria as a primary defence strategy against foreign nucleic acids, including phages and conjugative plasmids  <cit> .
in fact, crispr regions have been found to be integrated in response to infecting phages.
these regions are known to have hypervariable genetic loci due to the high diversities of the interspaced regions between the palindromic repeats and frequently match to phage and other extrachromosomal elements.
the presence of crisprs in genomes may affect short-term phenotypic changes and mediate long-term sublineage divergences  <cit> .
crispr regions were identified in all three of the malaysian s. typhi genomes using crisprs finder online .
the identified regions were found to be located around the crispr-associated protein cas <dig> and flanked by genes coding for alkaline phosphatase isoenzyme conversion aminopeptidase and clusters of crispr-associated genes, including cas and cse.
crispr_ <dig> was identified with the palindromic repeats, showing strikingly high similarities among all of the analysed s. typhi genomes with minor variations in spacers, indicating the important evolutionary conservation of the s. typhi strains.
apparently, the gene orders of crisprs are also conserved among the genomes.
the sporadic s. typhi strain st <dig> harbours a shorter crispr region  of 333 bp in length compared with ct <dig> , ty <dig>  and p-stx- <dig> .
the crispr region of st <dig> has identical palindromic repeats and a repeat length compared with the other genomes but lacks one spacer.
similar observations were also observed in ct <dig>  which had shorter spacers compared with the rest, which is in agreement with phylogenomic analysis, that ct <dig> is genetically related to st <dig> and cluster together with ty <dig> and p-stx- <dig>  which were found to possess crispr regions of identical lengths.
a previous study showed that the addition or deletion of spacers is able to modify the phenotype of phage resistance  <cit> .
alternatively, non-identity spacers have been suggested to mediate the interactions between crisprs and phages  <cit> .
the diversity of spacers in crisprs of s. typhi may be relevant to other interesting roles that are yet to be understood.
interestingly, additional crispr-like regions  were observed to be identical in the two closely related strains, bl <dig> and cr <dig>  we found that all of the s. typhi spacers that were analysed showed homology to many eukaryotic sequences, extrachromosomal sequences and phages, supporting the immunity roles of crisprs against phages and other incoming dna.
variations in crispr regions were identified, including those in bl <dig> and cr <dig>  providing evidence of evolutionary relevance that is useful for distinguishing between closely related strains and inferring ancestral relationships.
putative pathogenomic island harbouring zot, a potential novel virulence-determining factor
a total of eight  and  <dig> orfs  were identified in bl <dig> and cr <dig>  respectively, with eight commonly shared  orfs being identified  compared with the other genomes.
these clusters were predicted to be genomic islands  using islandviewer   <cit> .
gis are non-self-mobilising elements that code for proteins with diverse functions that may be integrated or excised, thus playing important roles in bacterial diversification and adaptation.
the predicted features of the identified gis resemble those of the previously reported pathogenicity islands, such as the presence of iss, integrases, transposases and deviated gc contents  <cit> .
the gis were found to harbour a novel gene coding for the zonular occluden toxin family protein , that is convergently oriented with respect to its flanking genes coding for a conserved domain protein and bacterial type ii and iii secretion proteins respectively, suggesting their involvement in the transportation of zot toxin into the host intestine.
intriguingly, the genomic elements together with the hypothetical proteins shared remarkable sequence similarities of >90% with yersinia pestis a <dig> and yersinia pestis co <dig>  <cit> , the causative agents of the black death .
we further screened for the prevalence of this gene in 41 s. typhi strains from diverse locations that were collected over a span of 25 years from 1983– <dig> and remarkably, zot is only present in strains bl <dig> and cr <dig> .
this is in agreement with our speculations that both the clonally related strains are responsible for the outbreak and carrier cases.
future studies could be carried out to fully characterise the gis and their relevant roles in the pathogenesis of s. typhi.figure  <dig>
schematic representative novel putative pathogenomic island harboring
zot
of cr <dig> and bl <dig>   <dig> -kb and  <dig> -kb genomic island fragment of strain bl <dig> and cr <dig>  respectively.
the schematic diagram shows the presence of various virulence associated genes detected, particularly zot and tblastx comparison  with yersinia pestis a <dig> and yersinia pestis co <dig>  the position of the regions in the genomes is labeled .
the arrow bars denode annotated genes as in the legend based on blast classification .
the arrow direction showed transcription direction of the gene.
the green and pink blocks above the arrow bars denode gc deviation.
salmonella spp.
including s. typhi, that carry zot gene are very limited  in our phylogenetic analysis.
the s. typhi zot homologue formed a monophyletic group with various zot homologues from enterobacteriacea, including e. coli and yersinia pestis but grouped more distantly with the well characterised zot protein in vibrio cholerae and neisseria meningitidis <cit>  .
zot was previously characterised as uropathogenic-specific protein in uropathogenic e. coli and a potentially important toxin in many pathogens that have received minimal attention  <cit> .
recently, campylobacter concisus, isolated from the patients with inflammatory bowel disease, was the only species among the campylobacteriales to harbour zot, suggesting its importance in gastrointestinal pathogenesis  <cit> .
however, little is understood on the mechanism of action of this protein, although earlier study of zot in v. cholera suggested that the protein act as a toxin that disrupts the integrity of the intestinal barrier by targeting the tight junction to increase tissue permeability  <cit> .figure  <dig>
phylogenetic tree of
zot
.
phylogenetic analysis of s. typhi zot with zot genes of  <dig> closely and distantly related bacteria strains using the approximate-maximum-likelihood method.
branch support was assessed with  <dig> bootstrap replicates.
bootstrap support values shown at each node.
our phylogenetic tree analysis  has further revealed that the zot homologue is surprisingly diverse, being present in the genomes of plant pathogens, such as ralstonia solanacearum and xanthomonas campestris, and opportunistic pathogens that are found in soil and water, such as providencia stuartii and pandoraea spp.
interestingly, this zot homologue was also found in lithoautotrophic bacteria, such as caminibacter mediatlanticus, sulphurimonas spp.
that clustered closely with the functionally validated zot in neisseria spp.
and campylobacter spp.
but more distantly related to enterobacteriacea, including our s. typhi.
other species such as pseudogulbenkiana ferrooxidans and beggiatoa alba, which were isolated from extreme environments, such as deep hydrothermal vents and sulphur springs were found to be more closely related and shared a common ancestry with s. typhi .
the importance of the necessity of gene targeting tight junctions in these extremophiles is yet to be understood but may be due to the role of the toxins in facilitating the penetration of the tissue layers that cover the host organisms that they colonise.
despite harboring four highly conserved domains leading to its putative assignment as a zot toxin, the s. typhi zot lacks a previously identified active domain  found in v. cholera.
nevertheless, a previous study has shown that the partial refolding of the denatured binding peptide of this domain did not prevent its specific binding to the zot receptor on caco- <dig> cells, demonstrating that the conformationally varied domain was still able to induce toxin activity  <cit> .
the possible mechanism of zot in s. typhi may be similar to that of v. cholerae, in which tight junction disassembly is induced through the activation of the proteinase-activated receptor  <dig>  <cit> .
however, comparison of the predicted zot protein tertiary structures of s. typhi with that of v. cholerae and n. meningitidis showed that they were highly variable, suggesting that the toxin might have different mechanisms of action despite maintaining some core binding activities.
the high similarities between the gis of the bl <dig> and cr <dig> strains further provide strong evidence of the common ancestry between the strains.
alarmingly, these strains are more diverse than was previously thought, emphasising our concern that such strains that have acquired new genes through hgt are circulating in the country or elsewhere.
microvariations distinguishing closely related strains
single nucleotide polymorphisms  were determined for the comparisons of the two highly related strains, bl <dig> and cr <dig>  we investigated the microvariations of the high quality non-synonymous single nucleotide polymorphism  that were identified.
in this study, only the potential functions altering the nssnp mutations were considered.
despite being highly related and similar, the genomes could be discretely distinguished from each other by  <dig> nssnps .
interestingly, the  <dig> nssnps that were identified, two were found on two highly related virulence determinant genes, rna polymerase sigma factor rpos and the vi polysaccharide biosynthesis protein tvie.
these two mutations are strain-specific and were not detected in any other s. typhi genomes that were analysed.
bl <dig> carried a mutant rpos , which is 100% similar to salmonella pullorum s <dig>  the causative agent of fowl typhoid  <cit> .
interestingly, unusual rpos gene was also observed in p-stx- <dig>  the rpos of p-stx- <dig> contains an additional 57 bp of a response regulator gaca fragment that was fused at the 5’ end of the gene, resulting in a longer rpos.
the gaca fragment was highly similar to the transposons tn10d teta and tetr, which are associated with regulatory and transcription signals, suggesting the occurrence of transposition events and possibly gene regulation.
the rpos mutant and its effects on s. typhi were first reported in the ty <dig> genome and attenuated strain ty21a.
these mutants, which were partially derived from natural mutations in ty <dig>  were found to affect the stress response and other related functions significantly  <cit> .
mutations in rpos are apparently advantageous to the strain for survival in the host during prolonged stress, allowing for the selection of more efficient transcription factors for survival and fitness during unfavourable conditions  <cit> .
rpos is commonly associated with the virulence regulation of pathogens because it regulates over  <dig> genes that are related to the stress response, such as the spv protein, which is involved in host cell survival, and vi-polysaccharide biosynthesis proteins in different osmolarities  <cit> .
in contrast, cr <dig> carried a mutant tvie .
tvie, which is encoded by spi- <dig>  is required for virulence capsular formation and acts as a protective antigen.
the antibody that responds to the vi-positive strain has been shown to be more virulent than that targeting the vi-negative strain  <cit> .
a role of the rpos protein in fine-tuning the synthesis of the vi polysaccharide in s. typhi has also been reported  <cit> , suggesting the possible close regulation of these two genes in modulating adaptation and virulence capabilities in different host environments.
the nssnps that were detected in these two closely related genes may indicate that their adaptive selection is important for host survival.
to address false-positive results, the nssnps were validated using a high-resolution melt  analysis and direct sequencing .
the unique hrm profiles of the wild-type strains and those containing the snp transition mutations in the normalised graph are shown in .
both snps showed unique melting profiles for the strains that were tested.
for the rpos snp, the transition mutation from c to t occurred in strain bl <dig>  and the separation of the melting profile began at ~81°c and ended at 84°c.
for the tvie snp, the transition mutation from c to t occurred in strain cr <dig>  and the separation of the melting profile began at ~74°c and ended at 78°c.
the results of this analysis were further confirmed by the direct sequencing of targeted loci .
the primers and hrm profiles that were developed may be useful for distinguishing between these two strains and as important markers for future surveillance.
another  <dig> high quality nssnps were also identified that mainly encoded for non-virulence factors .
twenty-three genes had well defined functions, including four that were involved in metabolism and  <dig> that played roles in cellular processes, signalling and transport.
the remainder were poorly characterised or had unknown functions.
additionally, out of  <dig> nssnps,  <dig>  mutants were detected in the carrier strain cr <dig>  suggesting the possible functional adaptation of this carrier strain in the host cell relative to its closely related strain.
however, most of the mutant snps that were carried by cr <dig> was not detected in p-stx- <dig> with the exception of the gene encoding the trehalose permease iic component.
in fact, these two carrier strains are genetically different, with the presence of  <dig> snps, including  <dig> nssnp.
notably, the nssnps were mainly found in the genes encoding proteins associated with virulence, metabolism, outer membrane, and other regulatory proteins.
however, nssnps were also observed in a large number of uncharacterised proteins, suggesting that independent genomic factors may have contributed to the carrier state.
it is important to note that the independent mutations in rpos and tvie as well as other genes in strains bl <dig> and cr <dig> support the postulation that these strains possibly diverged from a common ancestor.
we speculate that the source of this ancestral strain may be still circulating in the country.
therefore, national surveillance program and epidemiologic study are pivotal for effective microbial source tracking and dissemination control.
analysing the molecular effects of nssnps, protein structure modelling and molecular dynamics  simulation
point mutations that cause alterations in amino acids can have profound effects on the structural stability of proteins; hence, a study of their effects is necessary to understand their functionalities.
we have modelled both the native and mutant structures of the proteins.
out of the  <dig> native protein structure models that were generated by i-tasser  <cit> , the best structure with the highest confidence score  was collected and used for further investigations.
the modelled native rpos and mutant rpos structures showed good stereochemical properties, with  <dig> % and  <dig> % of the residues being within the most favourable region of the ramachandran plot, respectively, whereas the native and mutant tvie showed  <dig> % and  <dig> % of the residues in most favourable region of the plot, respectively.
all of the structures passed prosa model quality validation with z-scores  falling well within the range of those that are typically reported for native proteins of similar sizes from different sources .
molecular dynamics  simulation with a realistic aqueous solvent environment was performed to reveal the explicit solvent behaviours of the native and mutant structures, which could elicit the differences in their dynamics and stabilities.
the energy minimisation studies were performed for both the native and mutant structures, and the total energies of the native and mutant rpos achieved were − <dig>  j/mol and − <dig>  j/mol, respectively, whereas those of the native and mutant tvie achieved were − <dig>  j/mol and − <dig>  j/mol, respectively.
energy minimisation assessments provide clues regarding protein stability.
the deviation between two structures can be evaluated by the root mean square deviation .
the higher the rmsd value is, the greater the deviation will be between the native and mutant structures.
the rmsd value between the native and mutant rpos was found to be  <dig>  Å and that between the native and mutant tvie was  <dig>  Å. the native, mutant, superimposed protein structures at their corresponding positions for rpos and tvie are shown .
the rmsd values of the native and mutant structures were significantly similar for both rpos and tvie, suggesting similar levels of protein folding alterations.
we analysed the molecular effects and functional modifications based on several predictive tools  targeting various aspects of protein dynamics with confidence scores .
out of  <dig> predictive tools used, all found that the nssnp in rpos was deleterious , whereas in tvie,  <dig> predictive tools found that the nssnp to be deleterious,  <dig> to be neutral and one undetermined.
the nssnp in rpos showed psic score of  <dig>   with polyphen <dig>  <cit>  in addition to a probability score of  <dig>  with sift  <cit>  and a provean  <cit>  score of − <dig>  , suggesting that the nssnp could affect the protein drastically and result in functional modifications.
it is important to note that the deleterious effects could lead to positive or negative functional modifications.
to further validate the results, a method that was based on the hidden markov model  from panther  <cit>  was used.
the nssnp in rpos was found to be deleterious, but undetermined for tvie.figure  <dig>
modelled protein structures of rpos and tvie.
a <dig> showed native rpos with proline at position  <dig>  a <dig> showed mutant rpos with amino acid leucine at position  <dig>  a <dig> showed superimposed structure of rpos native structure  with mutant structure .
b <dig> showed native tvie with histidine at position  <dig>  b <dig> showed mutant tvie with amino acid tyrosine at position  <dig>  b <dig> showed superimposed structure of tvie native structure  with mutant structure .
the predictions of the nssnp in rpos by these tools are in agreement and show strong correlations among various methodologies.
the analysis with i-mutant  <dig>   <cit> , which is based on the support vector machine  and dgg stability changes, revealed that the nssnp in rpos may have led to its greatly increased stability , suggesting that the deleterious effects of the nssnp is favourable to the protein folding and structure and possibly lead to enhanced functions.
the ddg value is calculated from the unfolding gibbs free energy value of the mutated protein minus the unfolding gibbs free energy value of the wild-type strain , which is based on a trained cross-validation procedure using a comprehensive experimental database of protein mutations  <cit> .
these predicted results were consistent with the differences in the energy minimisation of md simulation, suggesting the favourable folding of the mutant structure.
native and mutant structures differ due to the specific properties of the residues that could disrupt the structure and function of the protein.
the mutant residue  of rpos is molecularly larger in size than the wild-type residue , which is a highly rigid, small molecule that is required to induce a unique backbone conformation.
however, the alteration from the small-sized secondary amine structure to a larger-sized primary aliphatic amine at this site can disturb its conformation and lead to bumps in the structure, in which the mutant residue is not in the correct position to make the typical hydrogen bond that is formed by the native residue .
the hydrophobicity of the wild type and mutant differ because the mutation introduces a very high hydrophobic residue in place of a less hydrophobic residue.
this can result in a loss of hydrogen bonding and may disturb proper protein folding.
the wild-type proline residue in rpos is well conserved and is located at the discrete compact three-helical domain within region  <dig> of the protein; however, no known mutant residues with similar properties were observed at this position in the other homologous sequences.
this region is the specific binding site of bacterial promoters containing an extended −10-promoter element and is primarily involved in the binding of the core rna polymerase in the holoenzyme.
the mutation in this important site could affect protein functioning pertaining to the transcription efficiency.
however, the predictions regarding the nssnp in tvie were contrasting, which suggest more benign effects, indicating that nssnps may have more modest effects on the functioning of this protein.
the nssnp in tvie may have decreased its stability  as predicted by i-mutant  <dig>   <cit> .
unlike rpos, the mutant residue  in tvie is not conserved at this position of the helical structure, and other non-similar residues were observed at this position with no known protein binding sites in other homologous sequences.
however, the size of the mutant residue  is larger compared with that of the native residue  despite the fact that both are polar and located at the surface of the protein.
this mutation also introduces a more hydrophobic residue at this position, suggesting that it could possibly alter the correct folding of the protein, subsequently affecting hydrogen bonding but with only modest effects.
these novel mutations could have important impacts on the pathogenesis and persistence of strains in the host.
although we could not determine the true extent of the effects of the nssnps on the protein functions, these data suggest that they alter the protein structures  considerably, potentially leading to the enhanced regulation of rpos and stress response in bl <dig> and reduced efficiency of tvie in the virulence capsular formation of cr <dig>  the close regulation of these two genes may be relevant to the virulence and persistence capabilities of the closely related strains that lead to the different clinical outcomes.
these data provide essential insights into the underlying molecular mechanisms upon mutations and serve as caveats for future functional gene knock-out studies.
CONCLUSIONS
we have thoroughly dissected the genome of s. typhi in association with three important epidemiological settings.
comparative genomics and phylogeny analyses have revealed that the strain that was associated with the large outbreak was highly related and shared common ancestry with the carrier strain.
these findings are supported by their common genomic features and uncommon gene repertoires, including dispensable genes, phages and an additional putative pathogenomic island harbouring virulence-related genes, and zot in particular.
apart from these, variations were also identified in t6ss and spi-related genes, insertion sequences, crisprs and nssnps among the studied genomes, which may be novel factors that contribute to the varied host adaptations and pathogenicities.
despite being highly similar, bl <dig> and cr <dig> may be distinguished by microvariations in their nssnps.
interestingly, the protein modelling and md simulation of the wild-type and mutant rpos and tvie suggest that the potential protein structure and functional modification was more stable in rpos, which plausibly leads to enhanced regulation and stress response.
on the other hand, the mutation in tvie was less stable than that of the wild type, which could potentially lead to lower capsular formation efficiency.
the close association of these virulence-related genes are relevant for long-term host persistence and adaptation, which serve as important caveats for further functional studies.
the analysis also revealed that spi- <dig>  which was previously thought to be relatively stable, is possibly prone to excision.
moreover, multiple regions of genomic plasticity were detected.
in particular, the discovery of new gis in the outbreak strain and the highly related carrier strain are of great concern epidemiologically.
these results suggest the plasticity and open pan-genome of s. typhi, indicating that the pathogen is more diverse than previously thought and that genes may have been acquired or transferred from one another through hgt, posing higher risk for effective disease control.
the genomic information that was obtained in this study provides novel insights into the pathogenesis and control of s. typhi, essentially, gene targets for vaccine development.

 typhoid fever is a human systemic infection that is caused by salmonella enterica subsp. enterica serovar typhi. this human-restricted and highly adapted pathogen is transmitted via the oral-faecal route. s. typhi is responsible for million infections and results in approximately deaths worldwide annually. the disease primarily causes acute systemic infection with life-threatening complications, and the recovering patient may develop into a chronic carrier state. typhoid is endemic, with periodic outbreaks and sporadic cases occurring in developing countries, particularly in southeast asia, south central asia, latin america and southern africa where sanitary conditions are poor. among the states of malaysia, kelantan has a significantly higher incidence of typhoid fever. a large typhoid outbreak occurred in kelantan state, which resulted in cases of infection and two deaths in a short period of 3 months. previously, pulsed-field gel electrophoresis revealed close genetic relatedness between a s. typhi strain that was isolated from an asymptomatic carrier in and a strain that originated from a patient during the outbreak in kelantan. human carriers are the main reservoirs of s. typhi transmission, but the genetic basis, and the underlying mechanisms in particular, are unclear. it has been suggested that a carrier strain will likely lack gene acquisition capabilities and have little fitness advantages compared with those strains causing symptomatic infections because the human reservoir is small and physiologically isolated. therefore, it is of great interest to know as to what extent these closely related strains differed or shared in its genomic contents despite being isolated from these two distinct epidemiological settings. in a s. typhi strain was isolated from a sporadic case in kuala lumpur. this genome has been sequenced earlier, and pfge analyses showed that this strain is more distantly related to the outbreak and carrier strains, but the epidemiological link is unknown. recent whole-genome sequencing of s. typhi has demonstrated that the pathogen shows limited genetic variation with little evidence of purifying selection, antigenic variation or recombination between isolates. this clonal pathogen, however, is associated with varying degrees of disease severity in different regions. previous pfge studies have also demonstrated genome size variations and distinct pfge patterns in relation with fatal and non-fatal typhoid cases. although the health conditions of the host cannot be completely ruled out, various reports have suggested that the gain and loss of genes through mutations and gene transfers that have occurred independently in different lineages have markedly contributed to the varying pathogenic potentials. however, these important factors are poorly understood because there is limited genomic information for s. typhi, particularly involving strains that are associated with diverse epidemiological settings. its genomic heterogeneity is likely due to the adaptation of the pathogen to the host and its exposure to mobile elements, such as bacteriophages. organisms having common core genomes could differ in their dispensable genes, reflecting their unique physiological and virulence properties. although not all genetic variations are essential for adaptation, some dispensable genes are believed to be responsible for conferring fitness advantages to the pathogen to thrive in its host. horizontal gene transfer is also thought to be the predominant force in bacterial evolution, which contributes to novel gene acquisition. the acquired genes provide new characteristics, which either aid in host adaptation and persistence or enhance virulence capabilities. a previous study on the pan-genome of salmonella enterica revealed that the pan-genome of all strains will continue to increase as new genomes are sequenced. with the availability of robust next-generation sequencing technologies, high quality whole genome sequences can be generated and analysed, which will be especially useful for capturing fine variations among highly conserved s. typhi strains. multiple whole genome sequence comparisons of closely related strains will not only lead to the better understanding of their relationships but also provide novel insights into the functional roles of strain-specific genes. in this study, we performed detailed and comprehensive comparative functional analyses of three previously sequenced genomes of s. typhi strains that were isolated from typhoid patients during a large outbreak in a sporadic case in and an asymptomatic carrier from malaysia, where typhoid is endemic. these malaysian s. typhi strains were compared with previously published s. typhi genomes with the following aims: 1) to determine and describe the genomes signatures and conserved and unique regions of the strains that were being studied; 2) to elucidate the phylogeny and genetic relatedness of these malaysian strains compared with other published strains using phylogenomic analysis; 3) to compare those strains that have been associated with various epidemiological settings, and particularly the regions of plasticity that may contribute to the varying pathogenic potentials; 4) to identify potential novel pathogenic factors that are harboured by the analysed strains; 5) to provide insight into the possible differential functionalities of the genes, focusing mainly on virulence- and persistence-related genes based on non-synonymous snps of closely related strains and particularly on carrier strains to gain insight into the persistence of the carrier state; and 6) to understand how the potential nssnps affect protein structures and functions. the data that are generated will be useful for the profiling of strains, marker development and the increased understanding of outbreak, sporadic and less studied asymptomatic typhoid carriage infection.  general genome signatures of s. typhi in association with outbreak, sporadic case and carrier the genomes of the three previously sequenced malaysian s. typhi strains from different sources were compared to identify potential genomic features that may help to elucidate the different disease outcomes. one strain was isolated during the largest outbreak in the country, one from a carrier during typhoid surveillance following the outbreak, representing the carrier state, and one from a sporadic case in the metropolitan area of kuala lumpur, which has a relatively lower incidence of typhoid. however, because of the limited numbers of strains studied and absence of information regarding the pan-genome of the s.typhi population, the genetic differences observed should be taken with caution. the aim of the detailed comparative analysis was to provide a better understanding and insights into the unusual genome dynamics of s. typhi, a highly clonal organism. previous genome sequencing analyses have generated high-quality assemblies with an average genome coverage of 100× for the malaysian s. typhi genomes, including bl, cr and st. the approximate predicted genome sizes and average guanine-plus-cytosine contents of all genomes ranged from  mb to  mb and % to %, respectively. these genomes form a single and circular chromosome with no plasmids detected. an in silico multi-locus sequence typing analysis classified both bl and cr as st and st as st which are the main sequence types that have been associated with the worldwide distribution of s. typhi out of the sts that have been identified to date. the predicted coding sequences of the genomes based on rast subsystem-based annotations varied from to with an average coding percentage of %. the average sizes of the cdss were similar, indicating that the size differences among the genomes are largely attributable to a number of cdss and intergenic regions. approximately 12% of the cdss were annotated as uncharacterised proteins. some of these genes were observed to vary from one strain to another. these “dispensable” genomes carried genes that were present in one or more strains and could even be unique to a single strain, indicating a possible open pan-genome of s. typhi. in general, the chromosomes of the three assembled genomes exhibited overall structural conservation and colinearity with each other as evidenced by the homologous and conserved regions that were shared and the very small strain-specific regions, which may harbour genes that are relevant to the specific adaptations and fitness advantages of each of the strains. these regions most likely represent dna that was acquired during events of hgt that may provide the strains with greater metabolic versatility or even virulence capabilities, as will be further discussed in the plasticity section.figure circular genomes representation map and genome comparison of s. typhi. the circle is divided into arcs representing the sequence of the chromosomes of all six genomes as labeled. the inner ring shows the coordinates in scale and total genomes size of the reference sequence, ct in scale, with black histogram bar representing gc content and purple-green histogram bar representing gc deviation. orthologs for each genomes with respect to ct are showed in order with the percentage of similarity based on nucleotide sequences as indicated on the legend located at the right of the figure. the outermost arch in red represents the location of spi to spi in relative to ct and labeled in red. the position of the genes related to pathogenicity and host persistent are shown and labeled in blue at the edge of the rings. comparative genomics of s. typhi genomic comparisons were performed on the three malaysian s. typhi strains and the three completed s. typhi genomes, ct, ty and p-stx-. the comparisons of the studied strains with the reference genomes allowed for the elucidation of novel and additional genes that are carried by the malaysian s. typhi genomes that may be of significance. we have analysed the shared and unique genes of all six genomes to determine their distinct virulence and pathogenic features. as expected, the six genomes exhibited high similarities and syntenies with each other with limited evidence of genomic rearrangements, which collectively indicate stable genomic structures. the majority of the orfs from the compared genomes were part of a conserved genomic core, in which orfs were shared among all of the genomes. these shared orfs provide clear evidence of conservation among the genomes of the s. typhi strains. the rest of the unshared orfs or accessory genes are present in one or more strains, which represent the salient differences in the genomes. most of these orfs that were harboured by each genome were annotated as hypothetical proteins. our extended homology analysis has shown that the remaining unshared orfs were likely to encode for functional proteins from diverse categories, including virulence-related proteins, secretory proteins, conserved domain proteins, transporter proteins and phage proteins among others. these considerable portions of the genomes could provide important functional clues for understanding the virulence and persistence of the pathogen more clearly, anticipating the need for extensive future studies focusing on their possible roles in bacterial pathogenesis. however, the numbers of shared and unshared orfs may have been underestimated because the genomes were incomplete. among the shared genes that were found between strains bl and cr uncommon orfs that encoded for the vi icm-f secretion protein, icm-f-related protein and type vi secretion protein evpb were identified whose products are related to the type vi secretion system. the genes shared 99% similarity with the type vi secretion protein of salmonella typhimurium strain d and were only found in bl and cr this protein was recently recognised as one of the main virulence determinants in burkholderia pseudomallei, legionella pneumophila and vibrio cholerae, but its function in s. typhi remains to be elucidated. t6ss genes are believed to be involved in either structural components of the secretory apparatus, secretory products or assisting with protein translocation; for example, providing the energy to push substrates through the channel of the apparatus. these genes are also proposed to be involved in surface reorganisation, enhancing adherence to epithelial cells, intracellular multiplication and human macrophage killing. other t6ss clusters were found intact as in reference genomes. the high similarities of the genetic contents of bl and cr with minor variations, and particularly the presence of unique accessory genes, are in agreement with the pfge pulsotype data, which revealed that both strains are genetically similar, showing a difference of only one band. the chromosome of salmonella enterica is commonly integrated with a large portion of horizontally acquired dna apart from its core, which are termed the salmonella pathogenicity islands. these acquired spis have led to divergence and host restriction similar to those in s. typhi. the identification of conserved spis and their variations have important implications in a wide range of microbiological applications, such as antigen and marker discovery and the identification of essential genes and their respective traits. in this study, we have annotated all spis and its genetic variant of s. typhi in the genome, which is characterised by its deviated gc content, flanking by trna genes and the presence of phages, integrases, recombinases and genes that are related to dna integration. although all the spis were detected in the genomes, there were marked variations. the presence of a large number of transposition-related genes in these spis suggests that the sites may be actively involved in the integration and transposition of genetic elements, which drive genetic variation. interestingly, our comparative analysis revealed that the carrier strain p-stx- lacks a ~10 kb prpz cluster and adjacent gene clusters harbouring orf with a deviated gc skew of % at spi- but remains fully intact in our carrier strain. previously, a prpz cluster deletion study showed that the mutant has a significantly lower survival rate compared with the parental strain, which may be due to a signalling pathway that controls the long-term survival of s. typhi in host cells, and particularly, the survival in human macrophages. in fact, our results support this study with the fine-tuned postulation that the deletion led to reduced virulence that enabled the carrier strain to coexist with the host; for example, in the tissue of gall bladder. this possibly explains why long-term survival in the macrophage is no longer necessary, which is presumably because the pathogens have colonised and persisted in other cells of the host during adaptation. however, the deletion was not detected in our carrier strain, suggesting that the genes may not be the only factors that are relevant to a carrier state. furthermore, the region is flanked by multiple transposases, integrases, ligases and uncharacterised proteins, which are known to be involved in transposition. additionally, genes coding for the dna mismatch repair protein mutc and transposase were identified both upstream and downstream of the cluster, suggesting that the region could have possibly been acquired earlier during horizontal gene transfer. the presence of a dna mismatch protein gene has been previously implicated to be involved in modulating recombination events by incorporating or inhibiting the transfer of mobile genetic elements. the deletion of the gene cluster and the presence of a large number of genes that are related to transposition indicate that spi- may be unstable and prone to excision similar to the precise excision of the crucial spi- which has been recently reported in s. typhi, suggesting that gene deletion may be important to the host adaptation of this organism, although independent acquisition or gene gain by other strains cannot be completely ruled out. these findings expand upon previous studies, which have reported that other spis are relatively stable in the genome, highlighting the importance of a future evaluation of the stability of the other spis. apart from these remarkable differences, the comparison of the two carrier strains, cr and p-stx- revealed that cr carries several additional genes that were not identified in p-stx- that encode for unknown functions and phages that are present in the phage region. almost all of the potential major virulence and persistence-associated genes have homologues in p-stx- suggesting that they are not specifically associated with the unique persistence of the carrier strains but are common to s. typhi. the genomic structure of the sporadic strain st is more conserved in comparison and has relatively fewer dispensable orfs, which are mainly genes that code for hypothetical proteins and phages, indicating the conservation of large numbers of genes, which is essential for strict host adaptation and virulence optimisation. phylogenomics of s. typhi revealed shared common ancestry we determined a core genome-based phylogeny by mapping query genomic sequences against ct into a single non-redundant alignment of,681 bp. the phylogenomic tree showed that the outbreak strain bl and carrier strain cr were closely related and could be differentiated by only snps. these data are in agreement with our pfge results that showed that both strains are highly related. the observed close phylogenetic relationship between these two strains is consistent with our earlier speculation that the large outbreak that occurred in shared a common ancestor strain with the typhoid carrier that may have been circulating for a long period in the country. it is challenging to determine how these two strains are related, considering their short-term evolutionary relationship. however, three evolutionary postulations are possible. first, the carrier strain may have been derived from the outbreak strain. second, the carrier strain may have long existed in a carrier who served as a reservoir and the source of the outbreak. finally, both of these strains may have diverged independently from a common ancestor, which was possibly harboured in a long existed carrier to give rise to two independent cases, considering the geographical proximities. notably, these two highly related strains with unique gene repertoires clustered with epidemiologically and geographically unrelated strains from india, russia, vietnam and senegal in west africa. interestingly, all of these strains, including bl cr p-stx- ty and e01- were subtyped as st. conversely, the malaysian sporadic s. typhi strain st clustered closely with the multidrug-resistant strain ct from vietnam together with other geographically related strains from indonesia and bangladesh, which were all subtyped as st suggesting the possible movement of clonally related strains among the southeast asian countries. such close genetic relatedness that is based on macrorestriction and snp typing analyses has been previously reported. the rest of the genomes were clustered as st from the analysis, we found no temporal or geographical signals, but the sequence types was highly correlated with the phylogenomic clustering. apparently, the cluster subtyped as st could be further differentiated into two clusters according to the phylogenomic analysis but was limited in the current mlst scheme for s. typhi. these results indicate that the phylogenomic analysis has much better resolution power compared with mlst in separating highly clonal strains in s. typhi. this information is essential for devising a better set of mlst alleles for improved molecular typing. these data further support the widespread distribution of st and st as the major genotypes that occur worldwide, although rarely, st and st have also been isolated in a previous study. it is important to note that the two carriers, cr and p-stx- belonged to two different subtypes, reflect the non-universality of the genotypes relevant for carrier state transformation.figure phylogenomic tree inferred by approximately-maximum-likelihood method from the aligned core genomes. multiple genomes alignments were generated by mapping genome sequences of the global s. typhi strains against ct at all sites relevant for phylogenomic analysis using realphy. phylogenetic tree was inferred via approximately-maximum-likelihood method using fasttreemp. strains studied were labelled in blue. bootstrap support values shown at each node. genome plasticity our analyses of the relationships among the strains and their genomic variations are further supported and extended to the study of the genome plasticity of s. typhi. in many organisms, genomic plasticity is commonly observed, but s. typhi has generally demonstrated few variations compared with other salmonella spp.. however, in our study, we detected considerable genetic variations, predominantly involving is elements, phages and crisprs. marked variations in the numbers and types of is elements were detected. is and is were both abundant in the three malaysian s. typhi genomes. the is elements, such as is have been widely used as molecular markers for subtyping due to their genome-wide distributions and high levels of diversity, but their roles in modulating the gene expression in s. typhi have yet to be clarified. recent studies on enterohaemorrhagic escherichia coli o have shown that the presence of is could play a role in the gene inactivation and immobilisation of incoming phages and plasmids, leading to the diversification and evolution of the bacterial genome. is elements have also been shown to affect the expression of neighbouring genes and induce genomic rearrangements. however, little is understood about their roles in modulating virulence and gene expression in s. typhi. the variations that were detected may provide clues on how these differences affect the virulence and fitness strategies of the pathogens. the s. typhi strains bl cr and st carry eight, seven and eight phages, respectively. we have identified substantial phage variations among the s. typhi genomes. one of the differentiating features was a distinct set of prophages that were harboured by both bl and cr which rendered them less unique compared with st with the exception of a few orfs that encoded for phages and hypothetical proteins. interestingly, the phages that were carried by st had relatively shorter in lengths compared with the phages that were identified in the other strains. as expected, both of the closely related strains had highly similar phage contents and carried an additional intact salmonella phage re- with uncommon orfs, which mainly contained genes coding for hypothetical proteins, phage proteins, prophage-like proteins, repressor proteins, excisionases, terminases and integrases. the variations that were detected in the numbers of predicted prophages and prophage-like regions illustrated the dynamics of phage gain and loss that distinguished one strain from the others, indicating that phages may play important roles in genomic diversity. by correlating phylogenomic analysis and whole genome sequence alignments, this region appears to show the typical gain and loss of sequences during the course of genomic evolution. this evolutionary relationship is consistent with the phage variations, providing a useful framework for investigating the relationships of strains and their respective phenotypes. the phage was likely acquired prior to the divergence of the common ancestor of both bl and cr through horizontal gene transfer rather than phage loss. alternatively, due to the advantageous roles of the hgt events, the phage proteins that were acquired by the strains could promote their in vivo survival and pathogenesis. among the detected phages, the two typical s. typhi phages, gifsy- and fels- were both found to be intact and conserved in all six of the s. typhi genomes, suggesting that these regions may play essential roles and provide fitness advantages to the pathogen. apart from the intact phages, incomplete phages, including burkholderia_phage_bcepmu, enterobacteria_phage_cdti and lactococcus_phage_bil were also observed in the six genomes, suggesting that these common phages were likely present in their common ancestors and may be relevant for host adaptation and survival. clustered regularly interspaced short palindromic repeats were recently found to be important in bacteria as a primary defence strategy against foreign nucleic acids, including phages and conjugative plasmids. in fact, crispr regions have been found to be integrated in response to infecting phages. these regions are known to have hypervariable genetic loci due to the high diversities of the interspaced regions between the palindromic repeats and frequently match to phage and other extrachromosomal elements. the presence of crisprs in genomes may affect short-term phenotypic changes and mediate long-term sublineage divergences. crispr regions were identified in all three of the malaysian s. typhi genomes using crisprs finder online. the identified regions were found to be located around the crispr-associated protein cas and flanked by genes coding for alkaline phosphatase isoenzyme conversion aminopeptidase and clusters of crispr-associated genes, including cas and cse. crispr_ was identified with the palindromic repeats, showing strikingly high similarities among all of the analysed s. typhi genomes with minor variations in spacers, indicating the important evolutionary conservation of the s. typhi strains. apparently, the gene orders of crisprs are also conserved among the genomes. the sporadic s. typhi strain st harbours a shorter crispr region of 333 bp in length compared with ct, ty and p-stx-. the crispr region of st has identical palindromic repeats and a repeat length compared with the other genomes but lacks one spacer. similar observations were also observed in ct which had shorter spacers compared with the rest, which is in agreement with phylogenomic analysis, that ct is genetically related to st and cluster together with ty and p-stx- which were found to possess crispr regions of identical lengths. a previous study showed that the addition or deletion of spacers is able to modify the phenotype of phage resistance. alternatively, non-identity spacers have been suggested to mediate the interactions between crisprs and phages. the diversity of spacers in crisprs of s. typhi may be relevant to other interesting roles that are yet to be understood. interestingly, additional crispr-like regions were observed to be identical in the two closely related strains, bl and cr we found that all of the s. typhi spacers that were analysed showed homology to many eukaryotic sequences, extrachromosomal sequences and phages, supporting the immunity roles of crisprs against phages and other incoming dna. variations in crispr regions were identified, including those in bl and cr providing evidence of evolutionary relevance that is useful for distinguishing between closely related strains and inferring ancestral relationships. putative pathogenomic island harbouring zot, a potential novel virulence-determining factor a total of eight and orfs were identified in bl and cr respectively, with eight commonly shared orfs being identified compared with the other genomes. these clusters were predicted to be genomic islands using islandviewer. gis are non-self-mobilising elements that code for proteins with diverse functions that may be integrated or excised, thus playing important roles in bacterial diversification and adaptation. the predicted features of the identified gis resemble those of the previously reported pathogenicity islands, such as the presence of iss, integrases, transposases and deviated gc contents. the gis were found to harbour a novel gene coding for the zonular occluden toxin family protein, that is convergently oriented with respect to its flanking genes coding for a conserved domain protein and bacterial type ii and iii secretion proteins respectively, suggesting their involvement in the transportation of zot toxin into the host intestine. intriguingly, the genomic elements together with the hypothetical proteins shared remarkable sequence similarities of >90% with yersinia pestis a and yersinia pestis co, the causative agents of the black death. we further screened for the prevalence of this gene in 41 s. typhi strains from diverse locations that were collected over a span of 25 years from 1983– and remarkably, zot is only present in strains bl and cr. this is in agreement with our speculations that both the clonally related strains are responsible for the outbreak and carrier cases. future studies could be carried out to fully characterise the gis and their relevant roles in the pathogenesis of s. typhi.figure schematic representative novel putative pathogenomic island harboring zot of cr and bl -kb and -kb genomic island fragment of strain bl and cr respectively. the schematic diagram shows the presence of various virulence associated genes detected, particularly zot and tblastx comparison with yersinia pestis a and yersinia pestis co the position of the regions in the genomes is labeled. the arrow bars denode annotated genes as in the legend based on blast classification. the arrow direction showed transcription direction of the gene. the green and pink blocks above the arrow bars denode gc deviation. salmonella spp. including s. typhi, that carry zot gene are very limited in our phylogenetic analysis. the s. typhi zot homologue formed a monophyletic group with various zot homologues from enterobacteriacea, including e. coli and yersinia pestis but grouped more distantly with the well characterised zot protein in vibrio cholerae and neisseria meningitidis. zot was previously characterised as uropathogenic-specific protein in uropathogenic e. coli and a potentially important toxin in many pathogens that have received minimal attention. recently, campylobacter concisus, isolated from the patients with inflammatory bowel disease, was the only species among the campylobacteriales to harbour zot, suggesting its importance in gastrointestinal pathogenesis. however, little is understood on the mechanism of action of this protein, although earlier study of zot in v. cholera suggested that the protein act as a toxin that disrupts the integrity of the intestinal barrier by targeting the tight junction to increase tissue permeability.figure phylogenetic tree of zot. phylogenetic analysis of s. typhi zot with zot genes of closely and distantly related bacteria strains using the approximate-maximum-likelihood method. branch support was assessed with bootstrap replicates. bootstrap support values shown at each node. our phylogenetic tree analysis has further revealed that the zot homologue is surprisingly diverse, being present in the genomes of plant pathogens, such as ralstonia solanacearum and xanthomonas campestris, and opportunistic pathogens that are found in soil and water, such as providencia stuartii and pandoraea spp. interestingly, this zot homologue was also found in lithoautotrophic bacteria, such as caminibacter mediatlanticus, sulphurimonas spp. that clustered closely with the functionally validated zot in neisseria spp. and campylobacter spp. but more distantly related to enterobacteriacea, including our s. typhi. other species such as pseudogulbenkiana ferrooxidans and beggiatoa alba, which were isolated from extreme environments, such as deep hydrothermal vents and sulphur springs were found to be more closely related and shared a common ancestry with s. typhi. the importance of the necessity of gene targeting tight junctions in these extremophiles is yet to be understood but may be due to the role of the toxins in facilitating the penetration of the tissue layers that cover the host organisms that they colonise. despite harboring four highly conserved domains leading to its putative assignment as a zot toxin, the s. typhi zot lacks a previously identified active domain found in v. cholera. nevertheless, a previous study has shown that the partial refolding of the denatured binding peptide of this domain did not prevent its specific binding to the zot receptor on caco- cells, demonstrating that the conformationally varied domain was still able to induce toxin activity. the possible mechanism of zot in s. typhi may be similar to that of v. cholerae, in which tight junction disassembly is induced through the activation of the proteinase-activated receptor. however, comparison of the predicted zot protein tertiary structures of s. typhi with that of v. cholerae and n. meningitidis showed that they were highly variable, suggesting that the toxin might have different mechanisms of action despite maintaining some core binding activities. the high similarities between the gis of the bl and cr strains further provide strong evidence of the common ancestry between the strains. alarmingly, these strains are more diverse than was previously thought, emphasising our concern that such strains that have acquired new genes through hgt are circulating in the country or elsewhere. microvariations distinguishing closely related strains single nucleotide polymorphisms were determined for the comparisons of the two highly related strains, bl and cr we investigated the microvariations of the high quality non-synonymous single nucleotide polymorphism that were identified. in this study, only the potential functions altering the nssnp mutations were considered. despite being highly related and similar, the genomes could be discretely distinguished from each other by nssnps. interestingly, the nssnps that were identified, two were found on two highly related virulence determinant genes, rna polymerase sigma factor rpos and the vi polysaccharide biosynthesis protein tvie. these two mutations are strain-specific and were not detected in any other s. typhi genomes that were analysed. bl carried a mutant rpos, which is 100% similar to salmonella pullorum s the causative agent of fowl typhoid. interestingly, unusual rpos gene was also observed in p-stx- the rpos of p-stx- contains an additional 57 bp of a response regulator gaca fragment that was fused at the 5’ end of the gene, resulting in a longer rpos. the gaca fragment was highly similar to the transposons tn10d teta and tetr, which are associated with regulatory and transcription signals, suggesting the occurrence of transposition events and possibly gene regulation. the rpos mutant and its effects on s. typhi were first reported in the ty genome and attenuated strain ty21a. these mutants, which were partially derived from natural mutations in ty were found to affect the stress response and other related functions significantly. mutations in rpos are apparently advantageous to the strain for survival in the host during prolonged stress, allowing for the selection of more efficient transcription factors for survival and fitness during unfavourable conditions. rpos is commonly associated with the virulence regulation of pathogens because it regulates over genes that are related to the stress response, such as the spv protein, which is involved in host cell survival, and vi-polysaccharide biosynthesis proteins in different osmolarities. in contrast, cr carried a mutant tvie. tvie, which is encoded by spi- is required for virulence capsular formation and acts as a protective antigen. the antibody that responds to the vi-positive strain has been shown to be more virulent than that targeting the vi-negative strain. a role of the rpos protein in fine-tuning the synthesis of the vi polysaccharide in s. typhi has also been reported, suggesting the possible close regulation of these two genes in modulating adaptation and virulence capabilities in different host environments. the nssnps that were detected in these two closely related genes may indicate that their adaptive selection is important for host survival. to address false-positive results, the nssnps were validated using a high-resolution melt analysis and direct sequencing. the unique hrm profiles of the wild-type strains and those containing the snp transition mutations in the normalised graph are shown in. both snps showed unique melting profiles for the strains that were tested. for the rpos snp, the transition mutation from c to t occurred in strain bl and the separation of the melting profile began at ~81°c and ended at 84°c. for the tvie snp, the transition mutation from c to t occurred in strain cr and the separation of the melting profile began at ~74°c and ended at 78°c. the results of this analysis were further confirmed by the direct sequencing of targeted loci. the primers and hrm profiles that were developed may be useful for distinguishing between these two strains and as important markers for future surveillance. another high quality nssnps were also identified that mainly encoded for non-virulence factors. twenty-three genes had well defined functions, including four that were involved in metabolism and that played roles in cellular processes, signalling and transport. the remainder were poorly characterised or had unknown functions. additionally, out of nssnps, mutants were detected in the carrier strain cr suggesting the possible functional adaptation of this carrier strain in the host cell relative to its closely related strain. however, most of the mutant snps that were carried by cr was not detected in p-stx- with the exception of the gene encoding the trehalose permease iic component. in fact, these two carrier strains are genetically different, with the presence of snps, including nssnp. notably, the nssnps were mainly found in the genes encoding proteins associated with virulence, metabolism, outer membrane, and other regulatory proteins. however, nssnps were also observed in a large number of uncharacterised proteins, suggesting that independent genomic factors may have contributed to the carrier state. it is important to note that the independent mutations in rpos and tvie as well as other genes in strains bl and cr support the postulation that these strains possibly diverged from a common ancestor. we speculate that the source of this ancestral strain may be still circulating in the country. therefore, national surveillance program and epidemiologic study are pivotal for effective microbial source tracking and dissemination control. analysing the molecular effects of nssnps, protein structure modelling and molecular dynamics simulation point mutations that cause alterations in amino acids can have profound effects on the structural stability of proteins; hence, a study of their effects is necessary to understand their functionalities. we have modelled both the native and mutant structures of the proteins. out of the native protein structure models that were generated by i-tasser, the best structure with the highest confidence score was collected and used for further investigations. the modelled native rpos and mutant rpos structures showed good stereochemical properties, with % and % of the residues being within the most favourable region of the ramachandran plot, respectively, whereas the native and mutant tvie showed % and % of the residues in most favourable region of the plot, respectively. all of the structures passed prosa model quality validation with z-scores falling well within the range of those that are typically reported for native proteins of similar sizes from different sources. molecular dynamics simulation with a realistic aqueous solvent environment was performed to reveal the explicit solvent behaviours of the native and mutant structures, which could elicit the differences in their dynamics and stabilities. the energy minimisation studies were performed for both the native and mutant structures, and the total energies of the native and mutant rpos achieved were −  j/mol and −  j/mol, respectively, whereas those of the native and mutant tvie achieved were −  j/mol and −  j/mol, respectively. energy minimisation assessments provide clues regarding protein stability. the deviation between two structures can be evaluated by the root mean square deviation. the higher the rmsd value is, the greater the deviation will be between the native and mutant structures. the rmsd value between the native and mutant rpos was found to be  Å and that between the native and mutant tvie was  Å. the native, mutant, superimposed protein structures at their corresponding positions for rpos and tvie are shown. the rmsd values of the native and mutant structures were significantly similar for both rpos and tvie, suggesting similar levels of protein folding alterations. we analysed the molecular effects and functional modifications based on several predictive tools targeting various aspects of protein dynamics with confidence scores. out of predictive tools used, all found that the nssnp in rpos was deleterious, whereas in tvie, predictive tools found that the nssnp to be deleterious, to be neutral and one undetermined. the nssnp in rpos showed psic score of with polyphen in addition to a probability score of with sift and a provean score of −, suggesting that the nssnp could affect the protein drastically and result in functional modifications. it is important to note that the deleterious effects could lead to positive or negative functional modifications. to further validate the results, a method that was based on the hidden markov model from panther was used. the nssnp in rpos was found to be deleterious, but undetermined for tvie.figure modelled protein structures of rpos and tvie. a showed native rpos with proline at position a showed mutant rpos with amino acid leucine at position a showed superimposed structure of rpos native structure with mutant structure. b showed native tvie with histidine at position b showed mutant tvie with amino acid tyrosine at position b showed superimposed structure of tvie native structure with mutant structure. the predictions of the nssnp in rpos by these tools are in agreement and show strong correlations among various methodologies. the analysis with i-mutant, which is based on the support vector machine and dgg stability changes, revealed that the nssnp in rpos may have led to its greatly increased stability, suggesting that the deleterious effects of the nssnp is favourable to the protein folding and structure and possibly lead to enhanced functions. the ddg value is calculated from the unfolding gibbs free energy value of the mutated protein minus the unfolding gibbs free energy value of the wild-type strain, which is based on a trained cross-validation procedure using a comprehensive experimental database of protein mutations. these predicted results were consistent with the differences in the energy minimisation of md simulation, suggesting the favourable folding of the mutant structure. native and mutant structures differ due to the specific properties of the residues that could disrupt the structure and function of the protein. the mutant residue of rpos is molecularly larger in size than the wild-type residue, which is a highly rigid, small molecule that is required to induce a unique backbone conformation. however, the alteration from the small-sized secondary amine structure to a larger-sized primary aliphatic amine at this site can disturb its conformation and lead to bumps in the structure, in which the mutant residue is not in the correct position to make the typical hydrogen bond that is formed by the native residue. the hydrophobicity of the wild type and mutant differ because the mutation introduces a very high hydrophobic residue in place of a less hydrophobic residue. this can result in a loss of hydrogen bonding and may disturb proper protein folding. the wild-type proline residue in rpos is well conserved and is located at the discrete compact three-helical domain within region of the protein; however, no known mutant residues with similar properties were observed at this position in the other homologous sequences. this region is the specific binding site of bacterial promoters containing an extended −10-promoter element and is primarily involved in the binding of the core rna polymerase in the holoenzyme. the mutation in this important site could affect protein functioning pertaining to the transcription efficiency. however, the predictions regarding the nssnp in tvie were contrasting, which suggest more benign effects, indicating that nssnps may have more modest effects on the functioning of this protein. the nssnp in tvie may have decreased its stability as predicted by i-mutant. unlike rpos, the mutant residue in tvie is not conserved at this position of the helical structure, and other non-similar residues were observed at this position with no known protein binding sites in other homologous sequences. however, the size of the mutant residue is larger compared with that of the native residue despite the fact that both are polar and located at the surface of the protein. this mutation also introduces a more hydrophobic residue at this position, suggesting that it could possibly alter the correct folding of the protein, subsequently affecting hydrogen bonding but with only modest effects. these novel mutations could have important impacts on the pathogenesis and persistence of strains in the host. although we could not determine the true extent of the effects of the nssnps on the protein functions, these data suggest that they alter the protein structures considerably, potentially leading to the enhanced regulation of rpos and stress response in bl and reduced efficiency of tvie in the virulence capsular formation of cr the close regulation of these two genes may be relevant to the virulence and persistence capabilities of the closely related strains that lead to the different clinical outcomes. these data provide essential insights into the underlying molecular mechanisms upon mutations and serve as caveats for future functional gene knock-out studies.  we have thoroughly dissected the genome of s. typhi in association with three important epidemiological settings. comparative genomics and phylogeny analyses have revealed that the strain that was associated with the large outbreak was highly related and shared common ancestry with the carrier strain. these findings are supported by their common genomic features and uncommon gene repertoires, including dispensable genes, phages and an additional putative pathogenomic island harbouring virulence-related genes, and zot in particular. apart from these, variations were also identified in t6ss and spi-related genes, insertion sequences, crisprs and nssnps among the studied genomes, which may be novel factors that contribute to the varied host adaptations and pathogenicities. despite being highly similar, bl and cr may be distinguished by microvariations in their nssnps. interestingly, the protein modelling and md simulation of the wild-type and mutant rpos and tvie suggest that the potential protein structure and functional modification was more stable in rpos, which plausibly leads to enhanced regulation and stress response. on the other hand, the mutation in tvie was less stable than that of the wild type, which could potentially lead to lower capsular formation efficiency. the close association of these virulence-related genes are relevant for long-term host persistence and adaptation, which serve as important caveats for further functional studies. the analysis also revealed that spi- which was previously thought to be relatively stable, is possibly prone to excision. moreover, multiple regions of genomic plasticity were detected. in particular, the discovery of new gis in the outbreak strain and the highly related carrier strain are of great concern epidemiologically. these results suggest the plasticity and open pan-genome of s. typhi, indicating that the pathogen is more diverse than previously thought and that genes may have been acquired or transferred from one another through hgt, posing higher risk for effective disease control. the genomic information that was obtained in this study provides novel insights into the pathogenesis and control of s. typhi, essentially, gene targets for vaccine development. 
our study has demonstrated that the pathogen is continually acquiring new genes through horizontal gene transfer in the process of host adaptation, providing novel insight into its unusual genomic dynamics. typhoid fever is an infectious disease of global importance that is caused by salmonella enterica subsp. s. typhi is strictly a human pathogen, and some recovered individuals become long-term carriers who continue to shed the bacteria in their faeces, thus becoming main reservoirs of infection. the understanding of these strains and virulence factors, and particularly the strain that is associated with the large outbreak and the less studied asymptomatic typhi carrier in the population, will have important impact on disease control. additionally, the comparative genomics analysis with all of the other completely sequenced s. typhi genomes revealed that strains bl and cr exhibit unusual genomic variations despite s. typhi being generally regarded as highly clonal. the uncommon genomic variations in the two closely related bl and cr strains suggests that s. typhi is more diverse than previously thought. s. typhityphoidgenomescomparative genomicspathogenzotspit6ssphylogeneticenterobacteriaceaesequenceevolutionproteinstrainvariationvirulenceinfectionprotein modellingissue-copyright-statement© the author 2014 a comparative genomics analysis combined with a phylogenomic analysis revealed that the strains from the outbreak and carrier were closely related with microvariations and possibly derived from a common ancestor.
 typhoid fever is an infectious disease of global importance that is caused by salmonella enterica subsp. enterica serovar typhi. this disease causes an estimated deaths per year and remains a serious global health threat. s. typhi is strictly a human pathogen, and some recovered individuals become long-term carriers who continue to shed the bacteria in their faeces, thus becoming main reservoirs of infection.  a comparative genomics analysis combined with a phylogenomic analysis revealed that the strains from the outbreak and carrier were closely related with microvariations and possibly derived from a common ancestor. additionally, the comparative genomics analysis with all of the other completely sequenced s. typhi genomes revealed that strains bl and cr exhibit unusual genomic variations despite s. typhi being generally regarded as highly clonal. the two genomes shared distinct chromosomal architectures and uncommon genome features; notably, the presence of a ~10 kb novel genomic island containing uncharacterised virulence-related genes, and zot in particular. variations were also detected in the t6ss system and genes that were related to spi- insertion sequences, crisprs and nssnps among the studied genomes. interestingly, the carrier strain cr harboured far more genetic polymorphisms compared with the closely related bl outbreak strain. notably, the two highly related virulence-determinant genes, rpos and tvie, were mutated in strains bl and cr respectively, which revealed that the mutation in rpos is stabilising, while that in tvie is destabilising. these microvariations provide novel insight into the optimisation of genes by the pathogens. however, the sporadic strain was found to be far more conserved compared with the others.  the uncommon genomic variations in the two closely related bl and cr strains suggests that s. typhi is more diverse than previously thought. our study has demonstrated that the pathogen is continually acquiring new genes through horizontal gene transfer in the process of host adaptation, providing novel insight into its unusual genomic dynamics. the understanding of these strains and virulence factors, and particularly the strain that is associated with the large outbreak and the less studied asymptomatic typhi carrier in the population, will have important impact on disease control. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords s. typhityphoidgenomescomparative genomicspathogenzotspit6ssphylogeneticenterobacteriaceaesequenceevolutionproteinstrainvariationvirulenceinfectionprotein modellingissue-copyright-statement© the author 2014 