Multiple independent origins of Shigella clones of Escherichia coli and convergent evolution of many of their characteristics (original) (raw)

Abstract

The evolutionary relationships of 46 Shigella strains representing each of the serotypes belonging to the four traditional Shigella species (subgroups), Dysenteriae, Flexneri, Boydii, and Sonnei, were determined by sequencing of eight housekeeping genes in four regions of the chromosome. Analysis revealed a very similar evolutionary pattern for each region. Three clusters of strains were identified, each including strains from different subgroups. Cluster 1 contains the majority of Boydii and Dysenteriae strains (B1–4, B6, B8, B10, B14, and B18; and D3–7, D9, and D11–13) plus Flexneri 6 and 6A. Cluster 2 contains seven Boydii strains (B5, B7, B9, B11, B15, B16, and B17) and Dysenteriae 2. Cluster 3 contains one Boydii strain (B12) and the Flexneri serotypes 1–5 strains. Sonnei and three Dysenteriae strains (D1, D8, and D10) are outside of the three main clusters but, nonetheless, are clearly within Escherichia coli. Boydii 13 was found to be distantly related to E. coli. Shigella strains, like the other pathogenic forms of E. coli, do not have a single evolutionary origin, indicating convergent evolution of Shigella phenotypic properties. We estimate the three main Shigella clusters to have evolved within the last 35,000 to 270,000 years, suggesting that shigellosis was one of the early infectious diseases of humans.

Shigella is a well known human pathogen that is prevalent in less developed countries where conditions of poor sanitation and personal hygiene increase the incidence of disease. The low infectious dose (10 cells) (1), allows the disease to be spread effectively by infected food or water, and also by person-to-person contact.

Historically Shigella was first described as Bacillus dysenteriae, the cause of bacilliary dysentery. It was clearly related to Bacillus (now Escherichia)coli but was given a different name because B. coli was known as a commensal organism and, at that time, there was no concept of a species including strains with such diverse characteristics. In the 1940s, four species of the new genus_Shigella_ were recognized. The traditional classification of_Shigella_ is a product of this history, which has been well documented by Ewing (2) and Bensted (3) and follows the recommendations of Ewing (4) with four species: S. boydii, S. dysenteriae, S. flexneri, and S. sonnei, also known as Shigella subgroups A, B, C, and D respectively. Nevertheless, Shigella and E. coli have always been considered to be very closely related, and for some time it has been clear that they are sufficiently similar to be placed in the same species (2, 5).

However, it was necessary to find characteristics that would enable Shigella strains to be distinguished from E. coli to confirm a diagnosis, and the major characteristics used were inability to ferment lactose and nonmotility. The matter was complicated in the 1940s when pathogenic strains were found that had some of the sugar fermentation and other characteristics of E. coli rather than Shigella, and these were treated as pathogenic forms of_E. coli_.

The two genera have been retained largely because of the serious nature of the disease shigellosis and its association with the name of the genus. However it now seems clear that Shigella strains are not a discrete group within E. coli and so do not even constitute a subspecies. For this reason, we refer to all Shigella strains as forms of E. coli (6). We consider it more accurate and less confusing to use a nomenclature that reflects well established natural relationships.

Subdivision into serotypes is based on the O antigen only, because Shigella strains lack the flagella H antigen and capsular K antigen also used in typing other E. coli. Forty-six Shigella serotypes are recognized. Flexneri serotypes 1 through 5 (12 including “subtypes”) are closely related and have a common basic structure (7) encoded by a common O antigen gene cluster (8). The differences are conferred by addition to the basic structure of glucosyl and/or_O_-acetyl residues (9). Flexneri serotypes 6 and 6A are not closely related genetically to the other Flexneri forms but crossreact with them serologically because of similarity in part of the O antigen structure (8).

Boydii and Dysenteriae have 18 and 13 well differentiated serotypes, respectively, with only occasional serological crossreaction between serotypes (2). Their repeat unit structures, where known, are unique among Shigella, although many Boydii and Dysenteriae O antigens are either identical to or related to a conventional E. coli O antigen (2). Sonnei has only one serotype, which has been shown to comprise a single clone (10). The Sonnei and Flexneri O antigens are not found otherwise in E. coli.

Earlier population genetic studies using multilocus enzyme electrophoresis (MLEE) showed that Shigella strains fall within_E. coli_ (11). We used MLEE and sequencing of the_mdh_ gene to demonstrate that Shigella strains do not form a single discrete set of strains within E. coli (6). In this study, we look at the relationships of all Shigella serotypes by sequencing four regions of the chromosome.

Materials and Methods

We used 46 Shigella strains to represent the known serotypes. Details are given in supplementary Table 1, which is published as supplemental data on the PNAS web site, www.pnas.org. For Dysenteriae, Flexneri, and Boydii, for which there are multiple serotypes, we refer to strains as D1, F1A, and B1 for Dysenteriae 1, Flexneri 1A, and Boydii 1, respectively, and so forth. For Sonnei, we use the abbreviation SS in tables and figures. We also used ECOR set strains ECOR7, ECOR28, ECOR30, ECOR33, ECOR37, ECOR50, ECOR54, and ECOR59, selected to represent nonpathogenic E. coli strains, and_Salmonella enterica_ LT2, which was used as an outgroup in the analysis.

Four segments of the chromosome were sequenced directly from PCR product. Primers for PCR and sequencing (supplementary Table 2) were based on the E. coli K-12 sequence and were selected taking into account conservation in sequence from related species. The primers were chosen to give two overlapping PCR amplicons for sequencing each of the four segments.

DNA sequences were assembled and edited by using programsphred, phrap, and consed (12). Further analysis was undertaken by using programs available from the Australian National Genomic Information Service (ANGIS) at The University of Sydney. Sequence comparisons were analyzed by using themulticomp package (13), which gives pairwise comparisons of DNA and amino acid sequences. multicomp calculates nucleotide diversity (π) by the method described by Nei and Miller (14) and average pairwise percentage difference. Calculation of synonymous and nonsynonymous substitution rates was done by using the program kindly provided by W.-H. Li (15). Molecular evolutionary relationships among each of the genes studied were examined by the neighbor-joining method of tree construction (16, 17), based on distance estimated by using the two-parameter method of Kimura (18). Phylogenetic trees and bootstrap analysis to determine the statistical stability of each node were done by using phylip (version 3.4 written by Joseph Felsenstein, Department of Genetics, University of Washington, Seattle, WA).

Results

A total of 7,160 bp from four regions was sequenced for each strain, except for D8 and B13, where one or two segments (600–950 bp), respectively, could not be amplified (see supplementary Table 3). The four regions are approximately equidistant on the chromosome and cover eight housekeeping genes. The 7,160 bp comprises 2,032 bp of_thr_B-_thr_C region; 1,486 bp of trpB-_trp_C region; 2,101 bp of_pur_M-_pur_N region; and 1,541 of_mdh_-_arg_R region. Details of the genes are presented in supplementary Table 3.

The sequence alignment for informative sites is shown in Fig.1. It is immediately clear that all but five of the Shigella strains fall into one of three major clusters, and that the same strains are in each cluster for all four regions of the chromosome. It is also clear that the strains within each cluster are very similar in sequence, but that there is little in common among the three clusters. Most of the variation involves base substitutions, but there are exceptions. B11 has a 52-bp deletion in the intergenic region between mdh and argR. B13 is particularly divergent in the purM-purN region, having a codon insertion (cysteine) at position 623 of the purN gene, and deletions of bases 13–14, 28–33, and 45–46 in the intragenic region downstream of the purN gene, resulting in a total of 2,094 bp for this strain.

Figure 1.

Figure 1

Alignment of the informative sites for the four regions sequenced. Note that for B13, D1, D8, D10, and Sonnei, which have very divergent sequences, the use of informative sites leads only to omission of many polymorphic sites.

Trees for the four regions are shown in Fig.2, and a tree for the combined data is shown in Fig. 3. Cluster 1 contains the majority of Boydii and Dysenteriae strains, plus F6 and F6A. Cluster 2 consists of seven Boydii strains (B5, B7, B9, B11, B15, B16, and B17) and D2. Cluster 3 consists of the Flexneri strains except F 6 and F6A, and B12. Sonnei D1, D8, and D10 are outside of the three main Shigella clusters but, nonetheless, are clearly within E. coli. B13 is only distantly related to other E. coli strains (Figs. 1 and 2, and supplementary Table 4). We applied bootstrap analysis using 1,000 replicas, and the high bootstrap values (100% for all three Clusters) confirm that the three clusters identified are robust.

Figure 2.

Figure 2

Phylogenetic trees generated by the neighbor-joining method for the four regions sequenced, showing different placement of the three clusters among ECOR set strains. The three major clusters are identified by numbers 1, 2, and 3 respectively. For detailed branching within each cluster, see Supplementary Fig. 4. Bootstrap values are percentages of 1000 replications and are indicated at the nodes. LT2 is used as the outgroup.

Figure 3.

Figure 3

Combined phylogenetic tree generated by the neighbor-joining method. Includes the sequences from the four regions sequenced to give a total of 7160 bp. F, Flexneri; B, Boydii; D, Dysenteriae; SS, Sonnei, followed by the serotype number. Bootstrap values are percentages of 1000 replications and are indicated at the nodes. LT2 is used as the outgroup.

Discussion

Historically, what is now Dysenteriae 1 (D1), the type form of Shigella, was discovered soon after E. coli and, because of its human pathogenesis, put in a separate genus. As more pathogenic strains were found, the absence of motility and lactose fermentation were found to be useful diagnostic criteria for dysentery bacteria placed in Shigella, although their relationship to_E. coli_ was always evident. In our previous work (6), we showed by MLEE and sequence of the mdh gene that some Shigella strains at least fall into discrete clusters within E. coli. The sequence data we now present from four regions around the chromosome consistently give three clusters of strains, each including strains from more than one traditional Shigella species. In addition, Sonnei falls outside of these clusters, as do D1, D8, D10, and B13, the last being very divergent. The consistency of the data from four regions gives great confidence in the new groupings, showing that the four existing species are not valid taxa even as groups within_E. coli_, raising questions as to their origins and convergent evolution.

Multiple Origins of Shigella Strains.

The presence of three major clusters and five forms not closely related to any other suggests that the Shigella phenotype has arisen eight times, or seven if one disregards B13 for this purpose, it being so divergent (supplementary Table 4) that it does not fall into E. coli, and is best treated as the only known representative of an unnamed species and left out of current discussion of the origins of the Shigella phenotype. We note that Manolov (19) considered “Shigella 13” strains not to be members of the genus_Shigella_, and that the classification of B13 as a_Shigella_ or E. coli strain has also been observed to be incorrect by Ewing (20) and Brenner (21, 22).

It is interesting to note that D1, the type form for Shigella and one of the prominent strains in outbreaks, and Sonnei, the most commonly isolated in industrialized countries, are both outside of the three main clusters. D1 is also unusual in production of Shiga toxin.

The current study included eight ECOR set E. coli strains, and, for three of the four regions, some fall between major Shigella clusters (Fig. 2), indicating that the three clusters are not derived from a single ancestral Shigella form. The distribution of Shigella strains among the ECOR set of 72 strains is more clearly seen in the earlier work using MLEE (6). There were cluster 1 and cluster 3 strains in that study, and it can be seen they are not closely related. Unfortunately, the limited number of Shigella strains used did not include any in cluster 2. The relationships among the ECOR set strains vary, as do the relationships of the three Shigella clusters and five nonclustered Shigella strains among themselves and with the ECOR set strains. This variation is not unexpected because E. coli has a quite high level of recombination (23). The mdh gene has been used in several studies, and a tree was constructed of all available E. coli sequences, including 31 ECOR set strains (6, 24), 32 pathogenic strains (6), and 21 isolates sampled from native rats in Australia (unpublished data) and all strains in this study. No non-Shigella strains fell within any of the three Shigella clusters, confirming that each is a cluster of closely related Shigella-only strains presumably derived from a single parental Shigella form.

Convergent Evolution of the Shigella Phenotype.

Shigella strains have a characteristic form of pathogenesis involving invasion of mucosal epithelium cells of the large intestine. The genes for the invasive property reside on a plasmid present in all Shigella strains tested thus far (25). However, Shigella strains have traditionally been defined and identified by characteristics either known or thought to be determined by chromosomal genes. In particular, Shigella strains lack catabolic pathways otherwise widely present in_E. coli_ (2). With few or no exceptions, lactose and mucate are not used, and lysine is not decarboxylated. Shigella strains are also nonmotile. Other E. coli strains are commonly positive for these properties. The independent evolution of the Shigella phenotype many times raises questions on the loss of motility and catabolic pathways by convergent evolution.

The studies by Al Mamun et al. (26) show that the basis for loss of motility in Shigella strains varies and can be due to deletion in the fliF operon (flagellar coding region) or an IS_I_ insertion mutation in the flhD operon. The lack of lactose fermentation was studied in several Shigella strains by Ito (27). Southern hybridization of lacY, lacA, and_lacZ_ genes showed that Flexneri 1 and 3 (cluster 3 in this study) and Bodyii 2 (cluster 2) and 4 (cluster 1) do not contain any of the lac genes and that Dysenteriae 1 has lacY and_lacA_ but not lacZ, whereas Sonnei has all three genes. It was further shown that, in Sonnei, the defect lies in the permase activity. These observations support the concept of multiple origins of the Shigella phenotype by convergent evolution.

It seems most likely that acquisition of the plasmid is a prerequisite for adoption of Shigella properties, and that the other characteristics, originally used together with symptoms of dysentery to define Shigella strains, are acquired by mutation later.

Convergent evolution of phenotypes is also seen for characters not universally present in Shigella strains. Shigella strains commonly lack ability to metabolize substrates otherwise widely utilized by_E. coli_. For example, mannitol utilization is present in 98% of typical E. coli, and its absence in many Shigella strains is presumably due to loss of function. There appear to have been seven independent events: three in outlier strains D1, D8, and D10, D2 in cluster 2, D5 and D7 in cluster 1, and the ancestor of the seven Dysenteriae strains (D3, D4, D6, D9, D11, D12, and D13) grouped together within cluster 1. Likewise, 96% of E. coli are positive for indole production, and loss of function is presumed to have occurred in Shigella cluster 3 and also in three of the outlier strains. A similar situation applies for utilization of xylose, rhamnose, and glycerol, present in 82%, 83%, and 89% of E. coli isolates, respectively [biochemical data from Ewing (2)].

The Shigella subgroups (2) are generally distinguished by differences in their biochemical profiles (Supplementary Table 5) although not unambiguously and often only in combination with serotyping. It is interesting that, when the strains are rearranged according to their respective sequence-based clusters, the profile of biochemical reactions is no less consistent than that of strains arranged by subgroup (compare supplementary Tables 5 and 6). In the case of indole production, there is a much better fit with the sequence-based clusters than the traditional subgroups. Mannitol fermentation, used to differentiate Boydii and Dysenteriae, naturally gives a very clear picture when strains are arranged in traditional subgroups. However, the far better correlation of indole production with genetic relatedness suggests that indole production would have been a better criterion for subdivision of Shigella strains. Note that the Dysenteriae strains in cluster 1 are grouped together, except D5 and D7, effectively allowing the use of mannitol to subdivide cluster 1.

It seems that Shigella strains have a general tendency to lose catabolic pathways, and this tendency gives rise to convergent evolution. Some of these characteristics have been used to define Shigella strains whereas others are not so widespread. One can speculate that Shigella strains, which live inside epithelial cells, do not need the range of catabolic pathways that generally characterize_E. coli_ as a species. In the case of lysine decarboxylase, it has been shown that this property is deleterious for Shigella such that there would be strong selection for its loss (28).

The situation with enteroinvasive E. coli (EIEC) strains is also interesting. Many have Shigella-like characteristics: they may be lac−, nonmotile, low level indole producing, and/or may not produce gas during fermentation. Their invasion properties are due to presence of a plasmid similar to those of Shigella strains. However, the characteristics of these strains are not well studied in comparison to those of Shigella strains, and we have not included them in this study.

The EIEC strains do not have the full set of characters that define Shigella strains, but they may be strains developing the full Shigella phenotype. None of the five EIEC strains included in an earlier study (6) was in any of the three clusters of Shigella strains in the_mdh_ tree referred to above, suggesting that the distinction between Shigella and EIEC strains is not entirely arbitrary.

Age of the Three Main Shigella Clusters.

The mean distance between alleles at synonymous sites (Ks) (supplementary Table 7) was used to estimate the time since descent from a common ancestor for each of the three main clusters of Shigella strains. We used the molecular clock rates of Whittam (29) and of Gutman and Dykhuizen (23) as done by Achtman et al. (30) to obtain a range of estimates. We calculated the Ks values for each gene and used the average of the eight genes as the distance for calculations. The estimated time is between 35,000 to 270,000 years for the three clusters. Clusters 1 and 2 have similar estimated times of origin 50,000–270,000 year ago, whereas cluster 3 apparently arose later, at 35,000 to 170,000 years ago. This is a very short time in comparison with that for diversification of E. coli in general, as illustrated by divergence of group A, one of the major ECOR set groups, estimated to have diverged from other major groups 8–22 million years ago (31).

These conclusions on the origins of Shigella can be compared with those for other human-specific pathogens for which similar data are available.

Yersinia pestis, the cause of plague, lacks variation in any of six gene segments sequenced and is proposed to be a clone that evolved from closely related Y. pseudotuberculosis 1,500 to 20,000 years ago (30). Plasmodium falciparum, a major cause of malaria, likewise lacks variation in five housekeeping genes and is estimated to have an origin as a clone after a selective sweep about 6,000 years ago (32).

The Mycobacterium tuberculosis complex, hich includes_M. bovis_, M. africanum,and M. microtis as well as M. tuberculosis, has virtually no silent substitution even between “species. There are many substitutions associated with drug resistance, but these are likely to have arisen by selection for drug resistance, and it is silent substitutions that give an estimate for age of the clone. On this basis, M. tuberculosis is proposed to have undergone an evolutionary bottleneck about 15,000 to 20,000 years ago (30) (this time frame presumably applies to the complex although there are significant host specificity differences between M. bovis, M. microtis, and M. tuberculosis, with M. africanum resembling M. tuberculosis in high degree of specificity for humans).

The Shigella strains of E. coli are essentially human specific as are M. tuberculosis and P. falciparum.

The Shigella strains and Y. pestis are in effect clones of more diverse species found in a range of animal hosts. However, the three clusters of Shigella strains carry much more neutral variation than the Y. pseudotuberculosis clone known as Y. pestis, or M. tuberculosis or Plasmodium falciparum. The dates for origin of Y. pestis and the sweep/bottlenecks of M. tuberculosis and P. falciparum are within or close to the start of the neolithic, with a change to settlement based population structure and rapid expansion of human population, and for Y. pestis a change in the ecology of the rat. The sequence variation within each of the three major Shigella clusters suggests that they arose much earlier, which is very interesting when looked at in relation to evolution of humans, the host species.

It has been proposed that that organisms like Shigella, which cause short-term infections with no carrier state and no alternative host, could not have survived in the paleolithic because early man lived in small bands with infrequent contact with other closed bands (33, 34). It is thought that such organisms could only survive after the development of agriculture with large settlements, a situation thought to have arisen in the neolithic about 10,000 years ago.

The estimated time of origin of Shigella forms of E. coli correlates with the origin and expansion of Homo sapiens rather than with the development of agriculture. It is not clear what this correlation means. It may be that Shigella strains had a greater capacity to survive in small hunter gatherer bands in the paleolithic than currently assumed. Perhaps we should also consider the possibility that we have underestimated the complexity of human populations 50,000 to 200,000 years ago during the rise to dominance of H. sapiens and the paleolithic expansion.

Antigenic Diversity in Shigella Strains.

There is considerable diversity of O antigen in the three clusters of Shigella strains. Those of cluster 3 are mostly flexneri strains with O antigens that differ only in phage-encoded properties, but each of the 28 strains in clusters 1 and 2 has a quite different O antigen, with the exception of F6 and F6A, which have variants of the same O antigen. Of the 27 unique O antigens, 9 are reported to be identical to antigens in other E. coli strains (serological data from Ewing, ref.2). That leaves 18 O antigens found only in one or other of these clusters. This diversity must have arisen by transfer of O antigen gene clusters to strains of these clusters during the 50,000 to 270,000 years since the last common ancestor of the current members of each cluster. The fact that 18 of the 27 O antigens have not been reported in other E. coli strains suggests that a high proportion of them may have come directly from outside of E. coli. D1 has one of the genes for O antigen synthesis on a plasmid (35), and the whole of the O antigen gene cluster of Sonnei is on a plasmid (36). The plasmid location is a good indication of recent lateral transfer and there is direct evidence for recent transfer of the Sonnei O antigen gene cluster from Plesiomonas shigelloides (J. Sheperd, L. Wang, and P.R.R., unpublished data). In any case, this expansion of O antigen diversity over a very short period relative to the time frame for diversification of the extant E. coli (see above) indicates that O antigen specificity is very important for Shigella strains. This study provides the first data that allow us to give a time frame for the generation of antigenic diversity in a closely related group of strains and an indication of the frequency of lateral transfer of such gene clusters.

The Classification of Shigella Strains.

As discussed above with the exception of Sonnei, which is a clone (10), the four generally recognized subgroups of Shigella do not represent natural groupings. Some of the relationships evident in Fig. 2 have been recognized for some time: e.g., that F6 is distinct from the other Flexneri strains, and that B13 is very distinct. Indeed Brenner (37) suggested the transfer of F6 to the Boydii subgroup but was rejected for historical reasons.

It is interesting to compare the results of this study with those of Dodd and Jones (38). They used 102 Shigella strains to determine relationships by using the methods of numerical taxonomy. In that study, Shigella strains formed a taxon not only distinct from E. coli but more closely related to Proteus/Providencia,Citrobacter, and Salmonella than to E. coli. We suggest that the reason for the discrepancy is that, in numerical taxonomy, the characters used are in general those developed by taxonomists because of their use in identification of generally recognized taxa. In a group like Shigella, which has been difficult to distinguish from E. coli, these characters will be strongly biased in favor of those that confirm the preferred taxonomy. This example illustrates the value of sequence data for determining true relationships in bacteria.

A good understanding of the origins of Shigella will be most important as more detail of the molecular basis for pathogenicity emerges. It will also help us to better understand clinical variation among Shigella forms.

Supplementary Material

Supplemental Data

Acknowledgments

We are grateful to Drs. Johanne Lefebvre and Pierre Harbec for providing strains shown in supplementary Table 1, and to Christine Dodd for providing additional information on the work reported in ref. 38. We thank David Ryan for technical assistance. We thank the anonymous referees for suggestions on improvement of the manuscript. The research was supported by a grant from the National Health and Medical Research Council (Australia).

Abbreviations

MLEE

multilocus enzyme electrophoresis

EIEC

enteroinvasive E. coli

Footnotes

This paper was submitted directly (Track II) to the PNAS office.

Data deposition: The sequences reported in this paper have been deposited in the GenBank database (accession nos. AF293105–AF293320).

Article published online before print: Proc. Natl. Acad. Sci. USA, 10.1073/pnas.180094797.

Article and publication date are at www.pnas.org/cgi/doi/10.1073/pnas.180094797

References

Associated Data

This section collects any data citations, data availability statements, or supplementary materials included in this article.

Supplementary Materials

Supplemental Data