Previous Article | Next Article ![]()
Journal of Clinical Microbiology, May 2008, p. 1741-1746, Vol. 46, No. 5
0095-1137/08/$08.00+0 doi:10.1128/JCM.02249-07
Copyright © 2008, American Society for Microbiology. All Rights Reserved.
,
The Wellcome Trust Sanger Institute, Hinxton, Cambridgeshire, United Kingdom,1 Oxford University Clinical Research Unit, Ho Chi Minh City, Vietnam,2 Department of Infectious Diseases, Leiden University Medical Centre, Leiden, The Netherlands,3 Centre National de la Recherche Scientifique, Saint Paul Lez Durance, France,4 Centre National de Référence des Salmonella, Institut Pasteur, Paris, France,5 Microbiological Diagnostic Unit-Public Health Laboratory, University of Melbourne, Melbourne, Victoria, Australia,6 Department of Biological Sciences, University of Calgary, Calgary, Alberta, Canada,7 Environmental Research Institute, University College Cork, Cork, Ireland8
Received 20 November 2007/ Returned for modification 8 February 2008/ Accepted 17 February 2008
|
|
|---|
|
|
|---|
Approaches involving direct interrogation of DNA sequence are extremely reproducible and can also supply phylogenetic information. For example, multilocus sequence typing (MLST) involves the comparative sequencing of DNA fragments from different bacterial genes and facilitates the calculation of evolutionary relationships based on sequence diversity (5). MLST has proved highly discriminative at the serovar level within S. enterica (33), but this approach identified few sequence types in serovar Typhi, indicating that this is a monophyletic, recently emerged human pathogen, which is perhaps only 10,000 to 50,000 years old (15, 24).
Recently, Roumagnac et al. (24) used denaturing high-performance liquid chromatography and DNA sequencing of 199 gene fragments to identify rare single nucleotide polymorphisms (SNPs) that defined specific haplotypes illustrative of the evolutionary history of serovar Typhi. SNP typing unequivocally assigns individual serovar Typhi strains to a particular haplotype and can simultaneously gather information about traits such as fluoroquinolone resistance (24, 36). This is because SNP analysis, which has been used extensively to analyze human DNA, can be performed using equipment that interrogates thousands of nucleotide changes simultaneously in a single DNA sample (10).
SNP typing can also potentially be exploited to analyze serovar Typhi populations circulating in a particular geographical region (36). Typhoid fever has a particularly high incidence (810/100,000) and mortality rate in Indonesia (27), where cases are sometimes associated with unusual clinical manifestations, including neurological complications (11). In addition, serovar Typhi isolated in Indonesia can express uncommon flagellum antigens that are not routinely found elsewhere, such as H:j (j) and H:z66 (z66). The j antigen, like the d antigen, is encoded by the chromosomal gene fliC, but the j fliC allele harbors a 251-bp deletion in the variable region of the gene, thus changing the dominant antigenic epitope on the flagellum (6). Expression of j antigen has been associated with a milder clinical presentation of typhoid fever (8). The z66 antigen was first described in 1981, and the genetic determinants for expression were recently shown to be located on a novel linear plasmid pBSSB1 (1, 9, 13). Thus, z66+ serovar Typhi harbor two different flagellin genes, although z66 expression is dominant due to the action of a repressor protein encoded on pBSSB1 that targets the chromosomal fliC gene (2).
It has been hypothesized that serovar Typhi strains from Indonesia have a higher level of genetic diversity than isolates from many other typhoid endemic regions (19). Furthermore, because z66+ strains harbor two flagellin genes they were postulated to be ancestral to global serovar Typhi (7). However, this conclusion was not supported by SNP typing of 105 strains, according to which, all seven z66+ strains were restricted to a local branch (24). In the present study we have used a high-throughput platform to assign 140 Indonesian serovar Typhi isolates to particular haplotypes. Furthermore, we have mapped the origin of some of these haplotyped isolates within a particular geographical area and have associated haplotypes with a particular flagellum antigen type.
|
|
|---|
Serovar Typhi genotyping assay. Primers for the study were specifically designed to target the previously described single nucleotide base changes (24) using MassARRAY software and were supplied by Sigma Genosys (Haverhill, United Kingdom). Genotyping was performed by using Sequenom homogenous mass extend (Sequenom, Inc., San Diego, CA). Multiplexed assays of 1 to 10 SNPs per multiplex were designed by using Sequenom Assay Design v3.0.2.0. Genomic DNA was prepared predominantly using the Wizard Genomic DNA purification kit (Promega) and was diluted to a final concentration of 4 ng/µl. The genomic DNA underwent a locus-specific PCR before being treated with shrimp alkaline phosphatase to remove unincorporated nucleotides. The Sequenom homogenous mass extend reaction was then performed according to the manufacturer's specifications. Oligonucleotide primers anneal adjacent to the SNP of interest, and a mixture of terminator nucleotides allow extension through the SNP site, thus creating allele specific extension products each with unique mass. Extension products were then desalted through the addition of an anion-exchange resin, before their weights were analyzed by matrix-assisted laser desorption ionization-time of flight mass spectrometry. In total, 84 DNA fragments were analyzed for calculation of the genotype for each individual bacterial isolate. Genotypes were assigned and viewed by using the SpectroCALLER and SpectroACQUIRE softwares (Sequenom), respectively. The specific haplotypes for all strains were assigned on the basis of the previously described haplotypes (24). A control strain E00-7866 (superscript c in Fig. 1) originating from Morocco was originally designated H45 but was redesignated H46 on the basis of our data. The specific polymorphisms in E00-7866 were verified by DNA sequencing, and this strain was found to be correctly assigned to H46. Furthermore, two Indonesian serovar Typhi controls (403Ty and 404Ty) were assigned to H59 by our analysis, rather than to H5 and H2, respectively, as previously described (24). It was subsequently shown that the SNPs differentiating 403Ty H5 and 404Ty H2 from H59 were introduced during laboratory manipulations and were therefore not investigated.
![]() View larger version (29K): [in a new window] |
FIG. 1. Minimal spanning tree based upon SNPs, demonstrating the relationship of the isolates used in the present study. This illustration is modified from Fig. 1 of Roumagnac et al. (24) and demonstrates the relationship of the 161 strains used here. The 18 serovar Typhi control strains with known haplotypes are labeled green and include CT18 (superscript a), Ty2 (superscript b), E00-7866 (superscript c), and AG3 (superscript d). The 140 experimental serovar Typhi strains are distinguished by red clusters, and the three serovar Paratyphi A isolates are colored yellow. The specific haplotype for each cluster is labeled (e.g., H45), and undetected haplotypes are gray. The size of each cluster corresponds to the number of strains identified in each group; this is demonstrated by the key in the bottom right corner.
|
|
|
|---|
The minimal spanning tree assembled from the raw data is shown in Fig. 1 and is adapted from Roumagnac et al. (24). Sequenom analysis of the 161 serovar Typhi isolates correctly assigned all control strains to the predicted haplotype (green circles; Fig. 1). The novel Indonesian serovar Typhi isolates were assigned to nine different haplotypes (red circles in Fig. 1), five of which also contained control isolates (green circles within red circles in Fig. 1). Significantly, the remaining 41 haplotypes reported previously in a global serovar Typhi collection, including three haplotypes described in Indonesia—H14, H58, and H70 (24), were undetected (gray circles in Fig. 1). Therefore, only a limited selection of the known haplotypes in the serovar Typhi population are circulating in this region.
The assay was additionally validated by assigning three experimental isolates to the ancestral haplotype, H45 (yellow circle in Fig. 1). These isolates were originally designated serovar Typhi, but when reserotyped they were found to be serovar Paratyphi A (588, 618, and 2585; see Table S1 in the supplemental material). All of the SNPs in the assay are specific for serovar Typhi; therefore, analysis of any other Salmonella serovars with this SNP assay would assign them to the ancestral haplotype H45.
Haplotypes H59 and H8 accounted for 53% (73) and 24% (34), respectively, of the serovar Typhi isolates recovered, and all haplotypes, except H84, included more than one isolate. Serovar Typhi haplotypes H42, H45, H50, H52, and H85 have been described previously as originating from geographically dispersed locations including Asia, Africa, and South America, providing evidence for the global dissemination of some strains (24). In contrast, H59 appears to be restricted to Indonesia and may be a recently evolved variant. H59 strains were identified in all five collections (see Materials and Methods) of experimental strains and were consistently isolated over a 30-year period (see Table S1 in the supplemental material). In addition, H42, H50, and H85 were also identified in different collections, with H85 representatives being isolated as far apart as 1987 and 2003 (see Table S1 in the supplemental material). This supports our previous observations that serovar Typhi strains can persist in a particular geographical region for periods extending several decades (24).
Serovar Typhi strains of haplotype H58 are frequently isolated in many parts of Southeast Asia and are associated with resistance to nalidixic acid, which may be a direct consequence of the common use of fluoroquinolones (17). We were unable to identify any H58 strains within these experimental Indonesian strains (superscript d in Fig. 1), although one fluoroquinolone-resistant strain belonging to H58 has been detected from a French traveler returning from Indonesia in 2003 (24). Importantly, our assay also interrogated the DNA for mutations in DNA gyrase, which leads to fluoroquinolone resistance (34). No such mutations in DNA gyrase were detected in this experimental population of serovar Typhi. This result indicates that there has been no recent clonal expansion of H58 in Indonesia, regardless of the fact that fluoroquinolone-resistant H58 strains may have been introduced in this country from near neighbors, such as Vietnam, where such strains are common. These data concur with those of other studies demonstrating the antimicrobial susceptibility of serovar Typhi and other enteric bacteria in Indonesia (14, 32) and may indicate a different epidemiological dynamic of the disease in this country.
Multiple serovar Typhi haplotypes are circulating in a local region.
Serovar Typhi described in collection 4 (see Table S1 in the supplemental material) were isolated between June 2001 and February 2003 in the Jatinegara district in Eastern Jakarta in a study by Vollaard et al. (35). To understand the circulation of serovar Typhi strains in this local region, we attempted to retrospectively correlate data, locating the residence of the individual typhoid patients with a specific serovar Typhi haplotype. A map demonstrating the distribution of individual serovar Typhi strains isolated within the
10.5-km2 Jatinegara district is shown in Fig. 2. We were able to pinpoint the serovar Typhi haplotype from all of the residence of the patients for which we had residential mapping data (54 of 84 cases in total). We identified eight circulating haplotypes in the Jatinegara district, providing unequivocally evidence that there are independent clusters of typhoid transmission involving distinct hapotypes in this small area. Four serovar Typhi isolates of identical haplotype (H59) were isolated in the same household (red circle within a black circle), suggesting direct horizontal transmission. The methodology presented here indicates that haplotyping data may prove to be a powerful tool to assist the tracking of serovar Typhi within a community by stratifying the bacterial population in a reliable and reproducible manner. A prospective study combining mapping information and haplotyping data would give insight into the mechanisms of transmission and long-term carriage within a population where typhoid is endemic.
![]() View larger version (34K): [in a new window] |
FIG. 2. Specific locations of serovar Typhi haplotypes isolated in Jatinegara, Jakarta. A map of the Jatinegara area in Eastern Jakarta highlighted to demonstrate the association of the household of the typhoid patient with the haplotype of the serovar Typhi bacteria isolated from the patient. Each specific haplotype corresponds to a colored circle as depicted on the right. H45* demonstrates the location of a household from which serovar Paratyphi A was isolated. The red circle enclosed by a black circle depicts the location of four individual typhoid patients from which H59 serovar Typhi was isolated.
|
![]() View larger version (21K): [in a new window] |
FIG. 3. Association of the fliC and fljB(z66) flagellin genes with specific haplotypes. (A) Diagram displaying the association of fliC allele with haplotype of 140 serovar Typhi strains. The strains with harbor the fliC gene which directs expression of the d antigen are colored gray the strains which harbor the fliC gene which directs expression of the j antigen are colored white. The specific haplotypes are labeled in association with Fig. 1, and the sizes of the clusters correspond to the number of strains located within each group; this is demonstrated by the key in the bottom right corner. (B) Diagram displaying the association of fljB(z66) gene with haplotype H59 in 140 serovar Typhi strains. The groups of strains lacking the fljB(z66) gene are white, and the strains harboring the fljB(z66) gene are shaded gray. The strains harboring the fljB(z66) gene are further subdivided in to groups carrying the fliC j allele (solid gray) or the fliC d allele (broken gray). The specific haplotypes are labeled in association with Fig. 1, and the sizes of the clusters correspond to the number of strains located within each group, as Fig. 3A.
|
|
|
|---|
10.5 km2) there are several different haplotypes of serovar Typhi circulating within the local population, suggesting the existence of independent transmission clusters. These data are in agreement with other studies, which report multiple strain types circulating within a specific location (28-30). Furthermore, we identified bacteria of the same haplotype isolated over a 30-year period, providing further support for a vital role of persistent carriage and/or prolonged dissemination of serovar Typhi in the environment. Typhoid fever remains a significant public health problem in many parts of Southeast Asia. Dissemination of serovar Typhi is of particular importance considering the wave of multiple drug-resistant strains that is currently spreading across Asia and the movement of people associated with current rapid economic growth. There is a pressing need for more effective epidemiological scrutiny of this organism, enabling better understanding of the nature of spread, which will ultimately facilitate control policies. Molecular methods appear to be the most robust, but such methods require standardization for the communication and analysis of the resulting data. The SNP-based approach is reproducible, results are applicable in any setting, and the data can be integrated into datasets produced in multiple laboratories using a variety of SNP detection technologies. This is in contrast to other schemes that detect strain differences but produce variable results in different laboratory settings. In addition, such schemes typically provide little or no information on phylogeny or phenotype and are unable to identify clonal expansion in a region.
The high-throughput system utilized here has been used to perform MLST on a number of Neisseria meningitidis reference strains (12). We have demonstrated the effectiveness of the same platform for detecting haplotypes in a genetically monomorphic gram-negative bacterial pathogen. Further flexibility can be anticipated when additional sensitivity is incorporated after more SNP variation is discovered.
Although the system has benefits, high-throughput SNP typing is currently not suitable for routine use in laboratories without specific infrastructure or financial support to run such assays. However, SNP detection technology has been driven in recent years by the need for massive throughput (currently >1 million SNPs for human genotyping studies). Microbial genotyping is much less demanding, with studies such as this providing epidemiologically informative data by assaying fewer than 100 SNP loci. With current high-throughput SNP detection technology, running costs for bacterial genotyping assays should be feasible for clinical laboratories, although setup costs may be prohibitive in some settings. The development of novel medium-throughput SNP detection technologies will be an important step toward the adoption of SNP-based microbial genotyping in clinical laboratories worldwide.
A practical solution for short-term development would be to house SNP typing facilities in regional reference laboratories linked via the internet through a global database. Simple "in the field" assays can also be established to detect a subset of discriminatory SNPs using a targeted approach. For example, MLST is performed routinely in many laboratories, while PCR analysis of nine gene fragments would have detected all of the haplotypes circulating in the present study. The resulting data would provide information about population structure in a particular region. Indeed, we are currently developing such a method based on PCR that can distinguish many of the major serovar Typhi haplotypes, which could potentially be universally adopted in clinical and research laboratories alike. However, it is important to note that the validity of this kind of simplified genotyping scheme depends on prior knowledge of the underlying phylogeny provided by comprehensive SNP typing studies.
It was originally hypothesized that z66+ serovar Typhi isolates from Indonesia were the precursor to global serovar Typhi, and it was assumed that the possession of a second phase of the flagella antigen was an ancestral state which has subsequently been lost (19). Conversely, an "out of Africa" hypothesis was proposed as the global source of Typhi (25), which would predict that the z66 gene was acquired later, possibly by horizontal gene transfer. We now know that z66 is present in only a single haplotype, and that the z66 flagellin gene [fljB(z66)] and fliC repressor [fljA(z66)] are located on a plasmid, indicating a relatively recent origin (1). The present study demonstrates that the acquisition of the pBSSB1 linear plasmid permitting the expression of z66 antigen most likely occurred only once and does not readily transfer to other genotypes. Furthermore, we suggest that pBSSB1 was horizontally transferred into an serovar Typhi strain carrying the d fliC allele, and over time the gene became truncated in some strains to form the j flagellin epitope. Expression of the z66 antigen is the default in serovar Typhi stains harboring the fljB(z66) locus, and fliC is effectively silenced (2). The truncated j allele can be found in three haplotypes, i.e., H50, H59, and H85, suggesting either that deletion has occurred spontaneously on several independent occasions or that homologous recombination has taken place between different haplotypes. However, it seems that strains harboring fljB(z66) preferentially possess the j fliC chromosomal allele.
In conclusion, our results demonstrate the adaptation of a global serovar Typhi genotyping scheme for a single country and ultimately a localized typhoid endemic region. The methodology offers a high level of sensitivity, thus allowing interrogation of phenotypes or pinpointing the locality of specific strains. The technique can in principle be used to define serovar Typhi circulating globally and is also potentially applicable to other bacterial pathogens. Our data prove that a number of distinct haplotypes of serovar Typhi can circulate in relatively small geographical area. However, somewhat paradoxically, we found that H58 strains that are currently circulating in other parts of South Asia and are associated with multiple drug resistance were not detected in Indonesia. Furthermore, strains harboring the z66 antigen are associated with the Indonesian archipelago, there is no evidence for the spread of these organisms to other countries where typhoid is endemic. It is likely that understanding the "Indonesian exception" would aid our understanding of the global epidemiology of typhoid fever. By combining haplotyping with phenotyping we were able to gain considerable insight into the population structure of serovar Typhi circulating in this region. Surveillance of strains using these methods combined with assessment of social and medical practices over a prolonged period will add vital information on how serovar Typhi is evolving and spreading in the human population.
This study was funded by The Wellcome Trust, London, United Kingdom. K.E.S. was supported by a Discovery Grant from the Natural Sciences and Engineering Research Council of Canada. M.A. is a Principal Investigator of the Scientific Foundation of Ireland (grant 05/FE1/B882) and the Max Planck Society for the Advancement of Science.
Published ahead of print on 5 March 2008. ![]()
Supplemental material for this article may be found at http://jcm.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»