http://orcid.org/0000-0001-9816-7674
Figures
Abstract
The pathogenesis of Salmonella enterica serovar Typhi (S. Typhi), the cause of typhoid fever in humans, is mainly attributed to the acquisition of horizontally acquired DNA elements. Salmonella pathogenicity islands (SPIs) are indubitably the most important form of horizontally acquired DNA with respect to pathogenesis of this bacterium. The insertion or deletion of any of these transferrable SPIs may have impact on the virulence potential of S. Typhi. In this study, the virulence potential and genetic relatedness of 35 S. Typhi isolates, collected from 2004 to 2013 was determined by identification of SPI and non-SPI virulence factors through a combination of techniques including virulotyping, Whole Genome Sequencing (WGS), and Variable Number of Tandem Repeats (VNTR) profiling. In order to determine the virulence potential of local S. Typhi isolates, 56 virulence related genes were studied by PCR. These genes are located in the core as well as accessory genome (SPIs and plasmid). Major variations among studied virulence determinants were found in case of SPI-7 and SPI-10 associated genes. On the basis of presence of virulence related genes, the studied S. Typhi isolates from Pakistan were clustered into two virulotypes Vi-positive and Vi-negative. Interestingly, SPI-7 and SPI-10 were collectively absent or present in Vi-negative and Vi-positive strains, respectively. Two Vi-negative and 11 Vi-positive S. Typhi strains were also analyzed by whole genome sequencing (WGS) and their results supported the PCR results. Genetic diversity was tested by VNTR-based molecular typing. All 35 isolates were clustered into five groups. Overall, all Vi-negative isolates were placed in a single group (T5) whereas Vi-positive isolates were grouped into four types. Vi-negative and Vi-positive isolates were mutually exclusive. This is the first report on the comparative distribution of SPI and non-SPI related virulence genes in Vi-negative and Vi-positive S. Typhi isolates with an important finding that SPI-10 is absent in all Vi-negative isolates.
Author summary
The distribution of virulence factors in S. Typhi can vary in isolates from different geographical regions and can have significant effect on the disease control. In this study, we have checked the distribution of 56 reported virulence associated factors in 35 local isolates of S. Typhi to identify any variations that can help in designing effective control strategies for typhoid. We have identified four naturally occurring variants which are simultaneously lacking SPI-7 and SPI-10, two adjacently located pathogenicity islands on S. Typhi chromosome. These isolates are not producing Vi capsular antigen hence the Vi based vaccines will not be effective against them. These findings highlight the need to develop typhoid vaccines specifically effective in Pakistan.
Citation: Liaquat S, Sarwar Y, Ali A, Haque A, Farooq M, Martinez-Ballesteros I, et al. (2018) Virulotyping of Salmonella enterica serovar Typhi isolates from Pakistan: Absence of complete SPI-10 in Vi negative isolates. PLoS Negl Trop Dis 12(11): e0006839. https://doi.org/10.1371/journal.pntd.0006839
Editor: Florian Marks, International Vaccine Institute, REPUBLIC OF KOREA
Received: February 22, 2018; Accepted: September 13, 2018; Published: November 30, 2018
Copyright: © 2018 Liaquat et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: The Whole Genome Sequencing (WGS) was funded by the Consolider-Ingenio 2010 program of the Spanish Government (CDS 2009-00006), and Innuendo project co-funded by the European Food Safety Agency (EFSA) (GP/EFSA/AFSCO/2015/01/CT2).The rest of the research work was funded by Higher Education Commission of Pakistan (HEC) and Vital Foundation of Vitoria-Gastiez. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Pathogenicity islands are distinct genetic components located on the pathogenic bacterial chromosomes. The pathogenesis of Salmonella enterica serovar Typhi (S. Typhi) is mainly accredited to the possession of horizontally acquired large DNA elements that transcribe in a coordinated manner to produce an array of symptoms for the onset of disease. S. Typhi is a human adapted pathogen. It causes a severe systemic infection, the typhoid fever, which is a serious worldwide public health problem. According to the World Health Organization (WHO) the annual global burden of typhoid fever is about 11–20 million new cases per year and 1% of which are fatal. More than 90% of typhoid fever cases occurred in Asia [1,2]. It is highly prevalent in Asia and Africa due to shortage of hygienic water and poor sanitation. It is also a significant travel-associated disease [3]. Therefore, S. Typhi infection poses substantial burden on healthcare system throughout the world especially in Southeast Asia (including Pakistan) and other endemic countries. Typhoid fever is clinically manifested by prolonged fever, abdominal discomfort, headache, and general lethargy. Early diagnosis and treatment using an appropriate antibiotic are essential for optimal management of typhoid fever, especially in children. Unfortunately, the emergence of multidrug-resistant S. Typhi strains causes difficulty in its treatment and poses a serious threat to future treatment options [4,5,6,7,8].
The complex pathogenesis of systemic Salmonella infections is associated with the presence of various defensive as well as offensive virulence factors. These factors contribute for its success as an intracellular human pathogen and participate at various stages of invasion, intracellular replication and survival within the host. Many of the Salmonella virulence genes are distributed on large genomic regions of 10–134 kb known as Salmonella pathogenicity islands (SPIs) [9]. The SPIs are characterized by a base composition different from the core genome and are often associated with tRNA genes and mobile genetic elements, like IS elements, transposons or phage genes [10]. Virulence factors encoded by SPI genes tamper with host cellular mechanisms and are thought to dictate the host specificity of different S. enterica serovars. Some virulence genes not located on SPIs such as the chromosomally-encoded phoP/Q (two component global regulator), rpoS (global stationary phase regulator) and fliA (RNA polymerase sigma factor for flagellar operon) also play important roles in the virulence of Salmonella [11,12].
Twenty one SPIs are known to date in Salmonella [13]. Out of these, 17 SPIs, 1 to 13 and 15 to 18 have been reported in S. Typhi. The largest of these islands, SPI-7 contains 134 kb of S. Typhi-specific DNA and carries biosynthesis genes (viaB locus) for the production of the Vi capsular antigen [14]. The Vi capsular antigen is a significant virulence factor for typhoid fever, as isolates positive for Vi production have higher rates of infection [15,16], and it continues to be the focus for prophylaxis for this disease. Volunteer studies have indicated that Vi-positive strains of serovar Typhi are more virulent in humans than Vi-negative isolates, although Vi production is not essential for the infection process in humans [17].
Vi-negative isolates (lacking Vi capsular polysaccharide antigen) of serovar Typhi have been reported in regions where typhoid fever is endemic. Previously, we have reported the existence of two types of naturally occurring Vi negative S. Typhi in Faisalabad region of Pakistan with partial (viaB operon only) or total absence of SPI-7 [18].
S. Typhi isolates which are genetically diverse with clonal expansion and genome variations have been reported in Malaysia and Southeast Asia. These genetic variations may be important for virulence [19,20,21,22]. The severity of the illness varies in different areas and this may be due to genetic diversity among the endemic strains [23]. Recent genome sequence projects demonstrated that S. Typhi strains showed limited genetic variation [24].
Genome sequence data on S. Typhi strains from different countries are required to clearly recognize their virulence potential. A well-known quick typing method is virulotyping that is used for detection and profiling in pathogenic bacteria. It increases our understanding of possible risk for human and animal infections. Virulotyping is a valuable tool for the characterization of Salmonella isolates [25]. In most of the studies involved in virulotyping, virulence factors with reported contributions to virulence were screened by PCR using gene specific primers. To find the presence of virulence genes, monoplex and multiplex PCR is routinely carried out. The distribution of SPIs has already been investigated in a reference strain CT18 [26] but such studies on naturally occurring Vi-positive and Vi-negative strains of S. Typhi are infrequent. This study was designed to find differences, if any, in the distribution of SPIs and related virulence factors as wells as non-SPI virulence factors of Vi-positive and Vi-negative isolates from Punjab, Pakistan. For this purpose, the comparative distribution of a significant number of virulence factors in clinical isolates of Vi-negative and Vi-positive S. Typhi collected from local sources was investigated.
Results
Distribution of virulence genes in Vi-positive and Vi-negative S. Typhi isolates
In this study 56 virulence related genes involved in mobility, secretion systems, metabolic regulation and toxin production were screened by PCR. The distribution of the virulence related genes among local isolates is presented in Table 1. In this study, S. Typhi isolates showed clearly distinct virulence-gene profiles: Vi antigen-positive and Vi antigen-negative according to the association of the virulence genes with SPI-7 and SPI-10.
Virulence determinant associated to core genome.
Ten virulence related genes associated with core genome were included in the present work. These genes encode outer membrane protein (tolC), adhesion factors (agfA, staA, stgA, sifA and fimA), sigma factor (rpoS, fliA), virulence transcriptional regulatory gene (phoP) and peptidase (STY1450). Each of the studied S. Typhi isolate harbored all these core genome associated virulence genes. No variation was observed in the genes located on core genome of S. Typhi.
Virulence determinant associated to accessory genome.
Virulence genes associated to fifteen SPIs (SPI-1-12, and SPI-16-18) previously reported to be related to virulence were investigated in this study. All of the major virulence markers located on SPI-1 to 6, 8, 9, 11, 12 and 16–18 (Table 1) were identified in 100% of local isolates. On the contrary, some of the virulence determinants located on SPI-7 (pilS, tviA and tviB) and all virulence determinants of SPI-10 (sefBCR, prpZ, prkY and prkX) were absent in 11% (n = 4) of S. Typhi isolates. Both tviA and tviB are required for Vi capsule synthesis. As these four isolates (11% of total isolates) were expected to be incapable of Vi expression, they were designated as ‘Vi-negative’. The SPI-7 associated phage-related sopE gene was detected in all S. Typhi isolates. However, partial deletion in sopE prophage was found in case of Vi-negative isolates.
Virulence related genes located on plasmids.
Three virulence genes (spvB, spvR and thrW) located on plasmids were also studied and 71% of strains did not possess any plasmid-associated virulence gene. Only one plasmid-associated gene trhW (located on pHCM1; encodes plasmid transfer protein involved in fimbrial regulation) was detected in 29% of isolates. This virulence gene was only detected in Vi-positive isolates.
Evaluating the presence or absence of SPI-7 and SPI-10 in Vi-negative isolates
Polymerase chain reaction.
Presence or absence of complete SPI-7 and SPI-10 was searched using primers specific to their flanking sides. SPI-7 is flanked by two partially duplicated tRNApheU loci. Primers DE0032-F and DE0083-R have been previously used to demonstrate the lack of an insertion at the tRNApheU locus [27,28]. These primers generate a PCR amplicon of 1275bp if the island is absent. SPI-7 is 134kb in length; therefore, the presence of the island is outside the constraints of the PCR. All Vi-negative and Vi-positive isolates failed to give any amplification with these primers, suggesting the presence of SPI-7. It was found that in our all Vi-negative isolates SPI-7 was present but it lacked viaB and pil operons. The primer pair SPI10up-F, SPI10up-R, was designed to amplify the upstream flanking side whereas SPI10dn-F, SPI10dn-R for downstream flanking side of SPI-10. Fig 1 shows the position of these primers in the genome of S. Typhi CT18.
Genes are depicted by arrows with gene designation indicated within arrows. Virulence genes are dotted and in black are phage-related genes. The schematic locations of primer pairs specific for SPI-10up, SPI-10dn, prkX, prkY, prpZ, sefB, sefC, and sefR are shown below genes. The position and length of expected PCR product (including region and genes investigated in this study) are mentioned according to the CT18 genome coordinates [26].
Interestingly, no PCR amplification was obtained with SPI-10 flanking side primers in case of Vi-negative isolates whereas Vi-positive isolates showed the respective amplifications for SPI-10 flanking side primers. The forward primer of upstream insertion site, SPI10up-F and reverse primer of downstream insertion site, SPI10dn-R were also used to amplify any genome sequences flanking the two insertion sites to verify the absence of SPI-10 but no amplification was observed. Therefore, it was concluded that SPI-10 was completely absent including the insertion sites from all our Vi-negative S. Typhi isolates. As expected, SPI-10 associated virulence genes (sefBCR, prpZ, prkY and prkX) were detected only in Vi-positive S. Typhi. In order to confirm the specificity of amplicons of SPI-10, associated virulence genes as well as amplified fragments of flanking sides of SPI-10 were sequenced. The sequences were then verified by BLASTn [29] against Salmonella enterica serovar Typhi sequences. All these sequences were confirmed as related to corresponding genes and region of SPI-10 of the S. Typhi genome.
Whole genome sequencing results.
Complete absence of SPI-10 in Vi-negative S. Typhi isolates was confirmed by the whole genome sequence analysis [29]. Altogether, 76.79 million read pairs of 100bp were obtained for the genome of 13 S. Typhi strains (2 Vi-negative and 11 Vi-positive) with an average of 5.9 million read pairs per sample, comprising average throughput of 1.18Gb per sample, as shown in supplementary table (S1 Table). All of the sequenced strains were having a minimum of 230x genome coverage, which is fairly good enough for a reliable comparison of bacterial genomes. The Vi-negative strains have 4062 complete and 49 partial genes (pseudogenes) as compared to the reference genes set (4395) of S. Typhi CT18 which is fairly low than the other Vi-positive strains. The NG50 statistics for the de novo assemblies is in the range of 145kb to 173kb. The results show that both Vi-negative strains (ST5 and ST25) have a higher level of divergence from the S. Typhi reference genome as the number of contig mismatches (>9300) and indels (>200) are higher than other strains, as shown in supplementary table (S1 Table). The comparative analysis of Nx, NGx, cumulative length, GC content and coverage histograms has been provided in supplementary figures S1 Fig–S31 Fig. Blastn [29] similarity search of assembled contigs with the S. Typhi reference genome confirmed the PCR findings that both Vi-negative strains (ST5 and ST25) did not have SPI-7 and SPI-10 except for some of the genes of SPI-7, as shown in Fig 2.
The outer circle corresponds to the reference genome of S. Typhi (Acc: AL513382) whereas inner circles depict the strains with clear depiction of absence of SPI-7 and SPI-10 in S. Typhi strain 5 and 25. The gene labels are displayed at 1500bp distance window size for figure adjustment.
These data validate our hypothesis of simultaneous absence of SPI-7 and SPI-10 from Vi-negative isolates. The WGS data of the 13 isolates of S. Typhi have been submitted in GenBank and the accession numbers were assigned as SAMN08195664, SAMN08195665, SAMN08195666, SAMN08195667, SAMN08195668, SAMN08195669, SAMN08195670, SAMN08195671, SAMN08195672, SAMN08195673, SAMN08195674, SAMN08195675, SAMN08195676.
Variable Number of Tandem Repeats (VNTR) based molecular typing of S. Typhi Local Isolates
On the basis of multiplex PCR results, all 35 S. Typhi isolates were grouped into 5 VNTR types (Fig 3). Overall, 3, 4 and 2 alleles were observed for TR1, TR2 and TR3, respectively. Therefore, all these nine type of amplicons from representative isolates were sequenced and number of repeats in each case was calculated by using tandem repeat finder (70). Table 2 briefly describes each observed VNTR type.
Lane 1 and 8:100bp DNA ladder (Invitrogen 15628–019; showing bands of 2072, 1500, 1400, 1300, 1200, 1100, 1000, 900, 800, 700, 600, 500, 400, 300, 200 and 100 bp), Lane 2–6: Vi-positive S. Typhi isolates, Lane 7: Vi- negative S. Typhi isolates.
Discussion
This study was focused on assessing the distribution of Salmonella pathogenicity islands (SPIs) in Vi-positive and Vi-negative variants of Salmonella enterica serovar Typhi and comparative virulotyping from local clinical samples. The virulence factors of the Salmonella serovars are mostly encoded by the Salmonella pathogenicity islands (SPIs). Two hallmarks of Salmonella pathogenesis, the host invasion and intracellular proliferation, are directly linked to the genes on SPIs. SPI-1 contains invasion genes whereas SPI-2 is required for intracellular pathogenesis. These SPIs encode two Type-III secretion systems (T3SS) and have crucial role for systemic S. enterica infections [30]. The invA, hilA, sipA, prgl and prgH genes are located on SPI-1 whereas spiC, sifA and sseB genes belong to SPI-2. In this study, we detected these genes in all S. Typhi isolates. SPI-3 is involved in intracellular survival and encodes a magnesium transporter (mgtB and mgtC) [31]. These genes were also detected in both types of isolates. Similar genetic homogeneity among Vi positive and Vi negative S. Typhi isolates was observed in virulence genes located on SPI-4 to 6, SPI-8, 9, 11, 12, and SPI-16-18. These findings are in accordance with previous studies [26,32].
Virulence factors on SPI-7 and SPI-10 are less stable than those associated with other SPIs [9]. In S. Typhi, SPI-7 is also referred as major SPI and represents the largest SPI with a size of 134 kb. It comprises of four parts: type IVB pilus operon, sopE prophage, Vi biosynthetic operon, and a 15 kb phage-like segment [26,27]. It codifies for the surface Vi polysaccharide antigen that contributes to virulence. However, lack of Vi expression can also be beneficial to some key steps of S. Typhi infectivity, for example, invasion, as Vi is the target of protective immune responses [33]. Vi-positive S. Typhi can benefit by inhibiting complement deposition at the bacterial surface and the post-phagocytic oxidative burst, thus resulting in reduced bacterial internalization and killing by phagocytes [34]. Serovar Typhi lacking Vi capsular polysaccharide antigen has been known and reported worldwide for several decades. Molecular evidence of the loss of Vi antigen has suggested that Vi-negative strains can be derived by the excision of SPI-7 or by a spontaneous base change in the viaB operon [28,35,36].
In the present study, all local isolate were carrying SPI-7 but some of them (11%) were deficient in SPI-7 associated genes (tviA, tviB and pilS). Results of this study indicated that 89% of analyzed S. Typhi isolates possessed viaB operon which is involved in biosynthesis of Vi-antigen. This virulence factor is also observed in S. Typhi isolates from typhoid patient’s blood, bone marrow as well as stool. It has been reported that most isolates were Vi-positive and only a few S. Typhi were Vi-negative [18,37,38]. The Vi-negative S. Typhi isolates that only lose expression of viaB operon have also been reported from India [37], and Nepal [38], but the absence of SPI-7 has been first time demonstrated in this study supported by whole genome sequencing (WGS).
SPI-10 is a pathogenicity island found next to the tRNAleuX gene at centisome 93. In the S. Typhi genome, this island corresponds to a 33 kb fragment carrying a full P4-related prophage, termed ST46 and the sefA-R chaperone-usher fimbrial operon [26,39]. The tRNAleuX region is a hypervariable hot spot for horizontal gene transfer in the Salmonella genus and different isolates from the same S. enterica serovar can exhibit significant variation in this region. Presence of mobile genetic elements and P4 phage play a major role in driving the variability of this region [39]. Genome sequence analysis identified three open reading frames carried by P4 like phage which have integrated within SPI-10 and are termed as prpZ gene cluster [39,40]. These are present in Ty2 and multi-drug resistant CT18 genomes of S. Typhi but absent in all other sequenced serovars of S. enterica. Several lines of evidence indicate that S. Typhi has acquired these three ORFs through Horizontal Gene Transfer (HGT). The prpZ gene cluster consists of three ORFs coding for proteins with homology to eukaryotic-type Ser/Thr protein phosphatases 2C (prpZ) and Ser/Thr protein kinases (prkY and prkX). It has been found to promote survival in macrophages [41,42].
The results in this study related to prpZ gene cluster were apparently distinct from previous reports [42,43,44]. Based on the findings of this study, we report for the first time that all our naturally occurring Vi-negative S. Typhi are lacking SPI-10. None of SPI-10 associated virulence genes (sefB, sefC, sefR, prpZ, prkX, and prkY) were detected in the Vi-negative isolates. In addition PCR targeting flanking regions of SPI-10 also yielded negative result. On the other hand all Vi-positive S. Typhi isolates yielded positive results under same conditions. WGS analysis of our Vi-negative and Vi-positive strains also confirms absence of SPI-7 and SPI-10.
Unlike SPI-7 which can be either partially or fully lost, SPI-10 is observed as completely absent in Vi-negative strains. It might be because its size is much smaller (33kb). It can be inferred that virulence related genes located on SPIs associated with prophages make up the major differences in gene contents among Vi-negative and Vi-positive S. Typhiisolates.
Repetitive DNA including VNTRs in bacterial genomes reflects their genomic diversity. VNTRs are commonly used to differentiate strains of homogenous clones. These are convenient for typing Salmonella enterica serovars and many previous studies have reported the VNTR based variability among S. Typhi isolates [45,46,47,48,49]. In this study, genetic diversity among studied S. Typhi isolates was tested by VNTR based molecular typing. All 35 isolates were clustered into five groups. Overall, all Vi-negative isolates were placed in a single group (T5) whereas Vi-positive isolated were grouped into four types. Both Vi-negative and Vi-positive isolates did not fall into same VNTR type. Unlike our isolates, S. Typhi strains from other countries including Nepal [45,46], China [48], Indonesia, Bangladesh, India, Singapore and Malaysia [45]showed more variety in their VNTR types. Octavia and Lan also searched VNTR profiles among S. Typhi strains collected from many countries except Pakistan. They also found variety of VNTR types among these global isolates [47]. Overall results of our study showed that Vi-negative and Vi-positive isolates were mutually exclusive in VNTR typing. Unlike S. Typhi strains from other counties, isolates from Pakistani showed less variety in their VNTR types.
S. Typhi clinical isolates from Pakistan showed genetic homogeneity in most of the virulence genes. They were found to be different only with regard to that of SPI-7 and SPI-10 regions. Our Vi-negative strains were deficient in both SPI-7 and SPI-10, so are reporting first time their simultaneous absencein S. Typhi Vi-negative strains. Interestingly, the absence of SPI-10 was also reported in an Indian strain (P-stx-12) that has intact SPI-7 [22]. Both SPI-7 and SPI-10 are associated with bacteriophages and mobile genetic elements, so this is probably the reason for their less stability.
It is hypothesized that Vi-negative isolates have evolved alternative ways to survive and colonize even without virulence genes associated with SPI-7 and SPI-10. The absence of both SPIs from the genome of Vi-negative strain could provide important functional clues for understanding the virulence and persistence of the pathogen, anticipating the need for extensive future studies focusing on their possible roles in bacterial pathogenesis. Our unique finding of SPI-10 deficient Vi-negative variants has further highlighted the importance of naturally occurring Vi-negative S. Typhi and thus the need to focus on universally present somatic antigens rather than Vi antigen for the preparation of successful vaccines effective against all isolates of S. Typhi.
Materials and methods
S. Typhi isolates
A total of 35 of S. Typhi isolates (Vi-negative (n = 4) and Vi-positive (n = 31)) were included in this study. All isolates were revived from National Institute for Biotechnology and Genetic Engineering (NIBGE, Pakistan) stock cultures previously collected different patients suffering from typhoid fever in various hospitals in Faisalabad region of Pakistan between 2002 and 2006 and stored in 20% glycerol at -20°C.
Identification and confirmation of isolates
Identification of isolates was performed by conventional biochemical methods after growing them on TSI medium (Merck, Germany). The strains were tested by agglutination for the presence of Vi antigen using Vi monovalent antisera (Monovalent Vi, Bio-Rad, France). These isolates were confirmed by regular and nested polymerase chain reaction (PCR) targeting fliC-d, tviA, and tviB genes as previously reported [50].
PCR virulotyping
For virulotyping analysis, all 35 S. Typhi isolates were screened for the presence/absence of 56 reported virulence genes (Table 3). Virulence determinants were categorized according to their locations and function. Virulence genes located on chromosome and associated with 15 SPIs (SPI-1 to 12, SPI-16 to 18) and other virulence genes not associated with any SPI were studied. These genes are responsible for encoding Type I and III secretion systems, invasions, adhesions, motility, sigma factor, virulence transcriptional regulatory protein, outer membrane protein, peptidase, phosphatase, fimbriae, pili, toxin and purine ribonucleoside efflux pump (Table 3). Most of the genes were detected by PCR using previously reported primer sequences whereas we designed new oligonucleotide primers for some of the virulence genes not studied previously by PCR (Table 3).
PCR primers were designed for this study from S. Typhi CT18 reference genome sequence from the NCBI Genbank (Accession number AL513382) [26]. Primer-BLAST [29] was used for designing the gene specific primers. Each 50μL of reaction mixture for each of the PCR included, 10μL of template DNA, contained 1.5mM MgCl2, 50nmol of each dNTP, 40pM of each (forward and reverse) primer and 2U of Taq DNA polymerase (Thermo Scientific, USA). PCR conditions were as follows: initial denaturation at 94°C for 5 min and 30 cycles of denaturation at 94°C for 1 min, annealing at temperature mentioned in Table 1 for each primer set for 1 min and extension at 72°C for 1 min, with a final extension at 72°C for 5 min using T100™ Thermal Cycler (Bio-Rad, USA). A non-template control was included in each run. The PCR products were analyzed by gel electrophoresis on 1.5% agarose stained with ethidium bromide under UV transilluminator.
Confirmation of absence of SPI-7 and SPI-10 in Vi-negative isolates
Presence of complete SPI-7 was searched by using flanking sides primers [27]. To confirm the complete absence of SPI-10 in Vi-negative isolates, primers specific to flanking sides of SPI-10 were designed (Table 4). The location of these primers on the physical map of the SPI-10 region is shown in Fig 1.
Whole Genome Sequencing (WGS) analysis.
Two Vi-negative and 11 Vi-positive S. Typhi strains were analysed by whole genome sequencing to confirm our findings. DNA was isolated by genomic DNA purification kit (Cat# K0152, Thermo Scientific, EU) and whole genome paired end sequencing was performed at aScidea Computational Biology Solutions, S.L., Barcelona (Spain) using Illumina HiSeq 2000 system. The sequencing adaptors were removed using built-in Illumina pipeline and raw data was quality checked using FASTQC program and quality filtering and trimming was performed using Trimmomatic v0.35. de novo assembly of the short reads was performed using Edena v3 [73] using default parameters. Comparative genome analysis of assembled contigs with S. Typhi reference (Accession number: AL513382)was performed using NCBI Blastn [29] and visualized using Circos program [74].
VNTR based molecular typing
Genetic relatedness of all studied isolates was investigated by Variable Number of Tandem Repeat (VNTR) analysis as described by Liu et al [45]. Three primer pairs, TR1, TR2, and TR3 were used to perform a multiplex PCR (Table 5), and all amplified fragments of different sizes were sequenced. The sequences were then verified by BLASTn[75] against Salmonella enterica serovar Typhi. In order to calculate copy number, Tandem Repeat Finder program [76] was used to analyze these sequences.
Supporting information
S1 Table. Whole genome sequencing and assembly statisitics of SalmonellaTyphi Strains.
https://doi.org/10.1371/journal.pntd.0006839.s001
(XLSX)
S1 Fig. Comparison of de novo assembled contig lengths between all strains.
The x-axis shows the percentage of length of the assembled genome for any strain and y-axis shows the length of contigs in kilobases used for a particular length percentage of assembled genome.
https://doi.org/10.1371/journal.pntd.0006839.s002
(TIF)
S2 Fig. Comparison of de novo assembled contig lengths between all strains with respect to the reference genome size [S. Typhi, Acc: AL513382].
The x-axis shows the percentage of length assembled with respect to the reference genome for any strain and y-axis shows the length of contigs in kilobases used for a particular length percentage.
https://doi.org/10.1371/journal.pntd.0006839.s003
(TIF)
S3 Fig. Comparison of growth of contig lengths.
On the x-axis, contigs are ordered from the largest to smallest. The y-axis gives the size of the x largest contigs in the assembly.
https://doi.org/10.1371/journal.pntd.0006839.s004
(TIF)
S4 Fig. Comparative distribution of # contigs with GC percentage in a certain range.
The x value is the GC percentage intervals. The y value is the number of contigs which GC content lies in the corresponding interval.
https://doi.org/10.1371/journal.pntd.0006839.s005
(TIF)
S5 Fig. %GC content of a ST5 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s006
(TIF)
S6 Fig. %GC content of a ST6 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s007
(TIF)
S7 Fig. %GC content of a ST7 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s008
(TIF)
S8 Fig. %GC content of a ST8 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s009
(TIF)
S9 Fig. %GC content of a ST11 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s010
(TIF)
S10 Fig. %GC content of a ST16 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s011
(TIF)
S11 Fig. %GC content of a ST17 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s012
(TIF)
S12 Fig. %GC content of a ST18 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s013
(TIF)
S13 Fig. %GC content of a ST19 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s014
(TIF)
S14 Fig. %GC content of a ST22 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s015
(TIF)
S15 Fig. %GC content of a ST25 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s016
(TIF)
S16 Fig. %GC content of a ST26 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s017
(TIF)
S17 Fig. %GC content of a ST38 de novo assembled contigs.
https://doi.org/10.1371/journal.pntd.0006839.s018
(TIF)
S18 Fig. Distribution of total contig lengths (y-axis) at different read coverage depths (x-axis, grouped in bins).
https://doi.org/10.1371/journal.pntd.0006839.s019
(TIF)
S19 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST5.
https://doi.org/10.1371/journal.pntd.0006839.s020
(TIF)
S20 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST6.
https://doi.org/10.1371/journal.pntd.0006839.s021
(TIF)
S21 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST7.
https://doi.org/10.1371/journal.pntd.0006839.s022
(TIF)
S22 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST8.
https://doi.org/10.1371/journal.pntd.0006839.s023
(TIF)
S23 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST11.
https://doi.org/10.1371/journal.pntd.0006839.s024
(TIF)
S24 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST16.
https://doi.org/10.1371/journal.pntd.0006839.s025
(TIF)
S25 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST17.
https://doi.org/10.1371/journal.pntd.0006839.s026
(TIF)
S26 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST18.
https://doi.org/10.1371/journal.pntd.0006839.s027
(TIF)
S27 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST19.
https://doi.org/10.1371/journal.pntd.0006839.s028
(TIF)
S28 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST22.
https://doi.org/10.1371/journal.pntd.0006839.s029
(TIF)
S29 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST25.
https://doi.org/10.1371/journal.pntd.0006839.s030
(TIF)
S30 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST26.
https://doi.org/10.1371/journal.pntd.0006839.s031
(TIF)
S31 Fig. Distribution of contig lengths (y-axis) at different read coverage depths (x-axis) of ST38.
https://doi.org/10.1371/journal.pntd.0006839.s032
(TIF)
Acknowledgments
The acknowledgments are due to National Institute for Biotechnology and Genetic Engineering (NIBGE), Faisalabad, Pakistan, and Department of Immunology, Microbiology and Parasitology, University of Basque Country, Spain, for technical facilities to conduct this research project and Consolider-Ingenio and Innuendo projects.
References
- 1. Crump JA, Luby SP, Mintz ED (2004) The global burden of typhoid fever. Bull World Health Organ 82: 346–353. pmid:15298225
- 2.
WHO (2018) Typhoid, Fact sheet.
- 3. Connor BA, Schwartz E (2005) Typhoid and paratyphoid fever in travellers. Lancet Infect Dis 5: 623–628. pmid:16183516
- 4. Bhan MK, Bahl R, Bhatnagar S (2005) Typhoid and paratyphoid fever. Lancet 366: 749–762. pmid:16125594
- 5. Wain J, Kidgell C (2004) The emergence of multidrug resistance to antimicrobial agents for the treatment of typhoid fever. Trans R Soc Trop Med Hyg 98: 423–430. pmid:15138079
- 6. Saleem MZ, Arshad A, Qayyum M, Shabbir MI, Ali A, et al. (2017) Changing trends in antibiogram and molecular analysis of quinolone resistant Salmonella Typhi isolates in Pakistan. J Infec Dis Treat 3.
- 7. Klemm EJ, Shakoor S, Page AJ Q, FN JK, Saeed DK, et al. (2018) Emergence of an extensively drugresistant Salmonella enterica serovar Typhi clone harboring a promiscuous plasmid encoding resistance to fluoroquinolones and third-generation cephalosporins. mBio 9 e00105–00118. pmid:29463654
- 8. Rasul F, Sughra K, Mushtaq A, Zeeshan N, Mehmood S, et al. (2017) Surveillance report on typhoid fever epidemiology and risk factor assessment in district Gujrat, Punjab, Pakistan. Biomedical Research 28: 6921–6926.
- 9. Hensel M (2004) Evolution of pathogenicity islands of Salmonella enterica. Int J Med Microbiol 294: 95–102. pmid:15493819
- 10. Schmidt H, Hensel M (2004) Pathogenicity islands in bacterial pathogenesis. Clin Microbiol Rev 17: 14–56. pmid:14726454
- 11. Haghjoo E (2004) Salmonella Typhi encodes a functional cytolethal distending toxin that is delivered into host cells by a bacterial-internalization pathway. Proc Natl Acad Sci USA 101: 4614–4619. pmid:15070766
- 12. Daigle F (2008) Typhi genes expressed during infection or involved in pathogenesis. J Infect Dev Ctries 2: 431–437. pmid:19745519
- 13. Sabbagh SC, Forest CG, Lepage C, Leclerc JM, Daigle F (2010) So similar, yet so different: uncovering distinctive features in the genomes of Salmonella enterica serovars Typhimurium and Typhi. FEMS Microbiol Lett 305: 1–13. pmid:20146749
- 14. Hashimoto Y, Li N, Yokoyama H, Ezaki T (1993) Complete nucleotide sequence and molecular characterization of viaB region encoding Vi antigen in Salmonella Typhi. J Bacteriol 175: 4456–4465. pmid:8331073
- 15. Hornick RB, Greisman SE, Woodward TE, DuPont HL, Dawkins AT, et al. (1970) Typhoid fever: pathogenesis and immunologic control. 2. N Engl J Med 283: 739–746. pmid:4916916
- 16. Hone DM, Attridge SR, Forrest B, Morona R, Daniels D, et al. (1988) A galE via (Vi antigen-negative) mutant of Salmonella Typhi Ty2 retains virulence in humans. Infect Immun 56: 1326–1333. pmid:3356467
- 17. Hone DM, Harris AM, Chatfield S, Dougan G, Levine MM (1991) Construction of genetically defined double aro mutants of Salmonella Typhi. Vaccine 9: 810–816. pmid:1759503
- 18. Baker S, Sarwar Y, Aziz H, Haque A, Ali A, et al. (2005) Detection of Vi-negative Salmonella enterica serovar Typhi in the peripheral blood of patients with typhoid fever in the Faisalabad region of Pakistan. J Clin Microbiol 43: 4418–4425. pmid:16145086
- 19. Baddam R, Kumar N, Thong KL, Ngoi ST, Teh CS, et al. (2012) Genetic fine structure of a Salmonella enterica serovar Typhi strain associated with the 2005 outbreak of typhoid fever in Kelantan, Malaysia. J Bacteriol 194: 3565–3566. pmid:22689247
- 20. Thong KL, Puthucheary S, Yassin RM, Sudarmono P, Padmidewi M, et al. (1995) Analysis of Salmonella Typhi isolates from Southeast Asia by pulsed-field gel electrophoresis. J Clin Microbiol 33: 1938–1941. pmid:7665677
- 21. Thong KL, Puthucheary SD, Pang T (1997) Genome size variation among recent human isolates of Salmonella Typhi. Res Microbiol 148: 229–235. pmid:9765803
- 22. Yap KP, Gan HM, Teh CS, Chai LC, Thong KL (2014) Comparative genomics of closely related Salmonella enterica serovar Typhi strains reveals genome dynamics and the acquisition of novel pathogenic elements. BMC Genomics 15: 1007. pmid:25412680
- 23. Parry CM, Hien TT, Dougan G, White NJ, Farrar JJ (2002) Typhoid fever. N Engl J Med 347: 1770–1782. pmid:12456854
- 24. Holt KE, Parkhill J, Mazzoni CJ, Roumagnac P, Weill FX, et al. (2008) High-throughput sequencing provides insights into genome variation and evolution in Salmonella Typhi. Nat Genet 40: 987–993. pmid:18660809
- 25. Huehn S, La Ragione RM, Anjum M, Saunders M, Woodward MJ, et al. (2010) Virulotyping and antimicrobial resistance typing of Salmonella enterica serovars relevant to human health in Europe. Foodborne Pathog Dis 7: 523–535. pmid:20039795
- 26. Parkhill J, Dougan G, James KD, Thomson NR, Pickard D, et al. (2001) Complete genome sequence of a multiple drug resistant Salmonella enterica serovar Typhi CT18. Nature 413: 848–852. pmid:11677608
- 27. Pickard D, Wain J, Baker S, Line A, Chohan S, et al. (2003) Composition, acquisition, and distribution of the Vi exopolysaccharide-encoding Salmonella enterica pathogenicity island SPI-7. J Bacteriol 185: 5055–5065. pmid:12923078
- 28. Wain J, House D, Zafar A, Baker S, Nair S, et al. (2005) Vi antigen expression in Salmonella enterica serovar Typhi clinical isolates from Pakistan. J Clin Microbiol 43: 1158–1165. pmid:15750077
- 29. Johnson M, Zaretskaya I, Raytselis Y, Merezhuk Y, McGinnis S, et al. (2008) NCBI BLAST: a better web interface. Nucleic Acids Res 36: W5–9. pmid:18440982
- 30. Hansen-Wester I, Hensel M (2001) Salmonella pathogenicity islands encoding type III secretion systems. Microbes Infect 3: 549–559. pmid:11418329
- 31. Blanc-Potard AB, Solomon F, Kayser J, Groisman EA (1999) The SPI-3 pathogenicity island of Salmonella enterica. J Bacteriol 181: 998–1004. pmid:9922266
- 32. Marcus SL, Brumell JH, Pfeifer CG, Finlay BB (2000) Salmonella pathogenicity islands: big virulence in small packages. Microbes and Infection 2: 145–156. pmid:10742687
- 33. Janis C, Grant AJ, McKinley TJ, Morgan FJ, John VF, et al. (2011) In vivo regulation of the Vi antigen in Salmonella and induction of immune responses with an in vivo-inducible promoter. Infect Immun 79: 2481–2488. pmid:21402763
- 34. Arricau N, Hermant D, Waxin H, Ecobichon C, Duffey PS, et al. (1998) The RcsB-RcsC regulatory system of Salmonella typhi differentially modulates the expression of invasion proteins, flagellin and Vi antigen in response to osmolarity. Mol Microbiol 29: 835–850. pmid:9723922
- 35. Bueno SM, Santiviago CA, Murillo AA, Fuentes JA, Trombert AN, et al. (2004) Precise excision of the large pathogenicity island, SPI7, in Salmonella enterica serovar Typhi. J Bacteriol 186: 3202–3213. pmid:15126483
- 36. Nair S, Alokam S, Kothapalli S, Porwollik S, Proctor E, et al. (2004) Salmonella enterica serovar Typhi strains from which SPI7, a 134-kilobase island with genes for Vi exopolysaccharide and other functions, has been deleted. J Bacteriol 186: 3214–3223. pmid:15126484
- 37. Saha MR, Ramamurthy T, Dutta P, Mitra U (2000) Emergence of Salmonella Typhi Vi antigen-negative strains in an epidemic of multidrug-resistant typhoid fever cases in Calcutta, India. Natl Med J India 13: 164.
- 38. Pulickal AS, Callaghan MJ, Kelly DF, Maskey M, Mahat S, et al. (2013) Prevalence and genetic analysis of phenotypically Vi- negative Salmonella Typhi isolates in children from Kathmandu, Nepal. J Trop Pediatr 59: 317–320. pmid:23598894
- 39. Bishop AL, Baker S, Jenks S, Fookes M, Gaora PO, et al. (2005) Analysis of the hypervariable region of the Salmonella enterica genome associated with tRNA(leuX). J Bacteriol 187: 2469–2482. pmid:15774890
- 40. Ong SY, Ng FL, Badai SS, Yuryev A, Alam M (2010) Analysis and construction of pathogenicity island regulatory pathways in Salmonella enterica serovar Typhi. J Integr Bioinform 7.
- 41. Thomson N, Baker S, Pickard D, Fookes M, Anjum M, et al. (2004) The Role of Prophage-like Elements in the Diversity of Salmonella enterica Serovars. Journal of Molecular Biology 339: 279–300. pmid:15136033
- 42. Faucher SP, Viau C, Gros PP, Daigle F, Le Moual H (2008) The prpZ gene cluster encoding eukaryotic-type Ser/Thr protein kinases and phosphatases is repressed by oxidative stress and involved in Salmonella enterica serovar Typhi survival in human macrophages. FEMS Microbiol Lett 281: 160–166. pmid:18312574
- 43. Deng W, Liou SR, Plunkett G 3rd, Mayhew GF, Rose DJ, et al. (2003) Comparative genomics of Salmonella enterica serovar Typhi strains Ty2 and CT18. J Bacteriol 185: 2330–2337. pmid:12644504
- 44. Lai SM, Le Moual H (2005) PrpZ, a Salmonella enterica serovar Typhi serine/threonine protein phosphatase 2C with dual substrate specificity. Microbiology 151: 1159–1167. pmid:15817783
- 45. Liu Y, Lee M-A, Ooi E-E, Mavis Y, Tan A-L, et al. (2003) Molecular typing of Salmonella enterica serovar Typhi isolates from various countries in Asia by a multiplex PCR assay on variable-number tandem repeats. J Clin Microbiol 41: 4388–4394. pmid:12958274
- 46. Acharya B, Aryal G, Banjade N, Sharma S, Pokharel B, et al. (2012) VNTR molecular typing of Salmonella enterica serovar Typhi isolates in Kathmandu valley. J Patho Nepal 2: 220–223.
- 47. Octavia S, Lan R (2009) Multiple-locus variable-number tandem-repeat analysis of Salmonella enterica serovar Typhi. J Clin Microbiol 47: 2369–2376. pmid:19535521
- 48. Wang H, Diao B, Cui Z, Yan M, Kan B (2016) Genotyping of Salmonella Typhi using 8-loci multi locus VNTR analysis Gut Pathog 8. pmid:27047570
- 49. Ramisse V, Houssu P, Hernandez E, Denoeud F, Hilaire V, et al. (2004) Variable number of tandem repeats in Salmonella enterica subsp. enterica for typing purposes. J Clin Microbiol 42: 5722–5730. pmid:15583305
- 50. Haque A, Ahmed J, Qureshi JA (1999) Early detection of typhoid by polymerase chain reaction. Ann Saudi Med 19: 337–340. pmid:17277535
- 51. Rahn K, De Grandis SA, Clarke RC, McEwen SA, Galán JE, et al. (1992) Amplification of an invA gene sequence of Salmonella Typhimurium by polymerase chain reaction as a specific method of detection of Salmonella. Mol Cell Probes 6: 271–279. pmid:1528198
- 52. Faucher SP, Porwollik S, Dozois CM, McClelland M, Daigle F (2006) Transcriptome of Salmonella enterica serovar Typhi within macrophages revealed through the selective capture of transcribed sequences. Proc Natl Acad Sci U S A 103: 1906–1911. pmid:16443683
- 53. Tran QT, Gomez G, Khare S, Lawhon SD, Raffatellu M, et al. (2010) The Salmonella enterica serotype Typhi Vi capsular antigen is expressed after the bacterium enters the ileal mucosa. Infect Immun 78: 527–535. pmid:19901065
- 54. Hughes LA, Shopland S, Wigley P, Bradon H, Leatherbarrow AH, et al. (2008) Characterisation of Salmonella enterica serotype Typhimurium isolates from wild birds in northern England from 2005–2006. BMC Vet Res 4: 4. pmid:18230128
- 55. Parvathi A, Vijayan J, Murali G, Chandran P (2011) Comparative virulence genotyping and antimicrobial susceptibility profiling of environmental and clinical Salmonella enterica from Cochin, India. Curr Microbiol 62: 21–26. pmid:20490498
- 56. Retamal P, Castillo-Ruiz M, Mora GC (2009) Characterization of MgtC, a virulence factor of Salmonella enterica serovar Typhi. PLoS One 4: e5551. pmid:19436747
- 57. Sanchez-Jimenez MM, Cardona-Castro N, Canu N, Uzzau S, Rubino S (2010) Distribution of pathogenicity islands among Colombian isolates of Salmonella. J Infect Dev Ctries 4: 555–559. pmid:21045367
- 58. Lim S, Jung J, Kim D (2007) The effect of γ radiation on the expression of the virulence genes of Salmonella Typhimurium and Vibrio spp. Radiat Phys Chem 76: 1763–1766.
- 59. Townsend SM, Kramer NE, Edwards R, Baker S, Hamlin N, et al. (2001) Salmonella enterica serovar Typhi possesses a unique repertoire of fimbrial gene sequences. Infect Immun 69: 2894–2901. pmid:11292704
- 60. Hashimoto Y, Itho Y, Fujinaga Y, Khan AQ, Sultana F, et al. (1995) Development of nested PCR based on the viaB sequence to detect Salmonella Typhi. J Clin Microbiol 33: 3082. pmid:8576386
- 61.
Sarwar Y (2008) Detection, molecular analysis and expressional studies of local Vi negative strains of Salmonella Typhi: Quaid-i-Azam University, Islamabad. 129 p.
- 62. Saroj SD, Shashidhar R, Karani M, Bandekar JR (2008) Distribution of Salmonella pathogenicity island (SPI)-8 and SPI-10 among different serotypes of Salmonella. J Med Microbiol 57: 424–427. pmid:18349359
- 63. Faucher SP, Forest C, Beland M, Daigle F (2009) A novel PhoP-regulated locus encoding the cytolysin ClyA and the secreted invasin TaiA of Salmonella enterica serovar Typhi is involved in virulence. Microbiology 155: 477–488. pmid:19202096
- 64. Skyberg JA, Logue CM, Nolan LK (2006) Virulence genotyping of Salmonella spp. with multiplex PCR. Avian Dis 50: 77–81. pmid:16617986
- 65. von Rhein C, Bauer S, Lopez Sanjurjo EJ, Benz R, Goebel W, et al. (2009) ClyA cytolysin from Salmonella: distribution within the genus, regulation of expression by SlyA, and pore-forming characteristics. Int J Med Microbiol 299: 21–35. pmid:18715828
- 66. Baumler AJ, Gilde AJ, Tsolis RM, van der Velden AW, Ahmer BM, et al. (1997) Contribution of horizontal gene transfer and deletion events to development of distinctive patterns of fimbrial operons during evolution of Salmonella serotypes. J Bacteriol 179: 317–322. pmid:8990281
- 67. Bishop A, House D, Perkins T, Baker S, Kingsley RA, et al. (2008) Interaction of Salmonella enterica serovar Typhi with cultured epithelial cells: roles of surface structures in adhesion and invasion. Microbiology 154: 1914–1926. pmid:18599820
- 68. Ngan GJ, Ng LM, Lin RT, Teo JW (2010) Development of a novel multiplex PCR for the detection and differentiation of Salmonella enterica serovars Typhi and Paratyphi A. Res Microbiol 161: 243–248. pmid:20381608
- 69. Cohen HJ, Mechanda SM, Lin W (1996) PCR amplification of the fimA gene sequence of Salmonella Typhimurium, a specific method for detection of Salmonella spp. Appl Environ Microbiol 62: 4303–4308. pmid:8953701
- 70. Way JS, Josephson KL, Pillai SD, Abbaszadegan M, Gerba CP, et al. (1993) Specific detection of Salmonella spp. by multiplex polymerase chain reaction. Appl Environ Microbiol 59: 1473–1479. pmid:8517741
- 71. Wu S, Li Y, Xu Y, Li Q, Chu Y, et al. (2010) A Salmonella enterica serovar Typhi plasmid induces rapid and massive apoptosis in infected macrophages. Cell Mol Immunol 7: 271–278. pmid:20473323
- 72. Sheikh A, Charles RC, Sharmeen N, Rollins SM, Harris JB, et al. (2011) In vivo expression of Salmonella enterica serotype Typhi genes in the blood of patients with typhoid fever in Bangladesh. PLoS Negl Trop Dis 5: e1419. pmid:22180799
- 73. Hernandez D, Francois P, Farinelli L, Osteras M, Schrenzel J (2008) De novo bacterial genome sequencing: millions of very short reads assembled on a desktop computer. Genome Res 18: 802–809. pmid:18332092
- 74. Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, et al. (2009) Circos: an information aesthetic for comparative genomics. Genome Res 19: 1639–1645. pmid:19541911
- 75. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. J Mol Biol 215: 403–410. pmid:2231712
- 76. Benson G (1999) Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res 27: 573–580. pmid:9862982