- Research article
- Open Access
A tandem repeats database for bacterial genomes: application to the genotyping of Yersinia pestis and Bacillus anthracis
BMC Microbiology volume 1, Article number: 2 (2001)
Some pathogenic bacteria are genetically very homogeneous, making strain discrimination difficult. In the last few years, tandem repeats have been increasingly recognized as markers of choice for genotyping a number of pathogens. The rapid evolution of these structures appears to contribute to the phenotypic flexibility of pathogens. The availability of whole-genome sequences has opened the way to the systematic evaluation of tandem repeats diversity and application to epidemiological studies.
This report presents a database (http://minisatellites.u-psud.fr) of tandem repeats from publicly available bacterial genomes which facilitates the identification and selection of tandem repeats. We illustrate the use of this database by the characterization of minisatellites from two important human pathogens, Yersinia pestis and Bacillus anthracis. In order to avoid simple sequence contingency loci which may be of limited value as epidemiological markers, and to provide genotyping tools amenable to ordinary agarose gel electrophoresis, only tandem repeats with repeat units at least 9 bp long were evaluated. Yersinia pestis contains 64 such minisatellites in which the unit is repeated at least 7 times. An additional collection of 12 loci with at least 6 units, and a high internal conservation were also evaluated. Forty-nine are polymorphic among five Yersinia strains (twenty-five among three Y. pestis strains). Bacillus anthracis contains 30 comparable structures in which the unit is repeated at least 10 times. Half of these tandem repeats show polymorphism among the strains tested.
Analysis of the currently available bacterial genome sequences classifies Bacillus anthracis and Yersinia pestis as having an average (approximately 30 per Mb) density of tandem repeat arrays longer than 100 bp when compared to the other bacterial genomes analysed to date. In both cases, testing a fraction of these sequences for polymorphism was sufficient to quickly develop a set of more than fifteen informative markers, some of which show a very high degree of polymorphism. In one instance, the polymorphism information content index reaches 0.82 with allele length covering a wide size range (600-1950 bp), and nine alleles resolved in the small number of independent Bacillus anthracis strains typed here.
The polymorphism associated with tandem repeats has been instrumental in mammalian genetics for the construction of genetic maps and still is the basis of DNA fingerprinting in forensic applications. Tandem repeats are usually classified among satellites (spanning megabases of DNA, associated with heterochromatin), minisatellites (repeat units in the range 6-100 bp, spanning hundreds of base-pairs) and microsatellites (repeat units in the range 1-5 bp, spanning a few tens of nucleotides).
More recently, a number of studies have supported the notion that tandem repeats reminiscent of mini and microsatellites are likely to be a highly significant source of very informative markers for the identification of pathogenic bacteria even when these pathogens are recently emerged, highly monomorphic species [1–5]. This probably reflects the important contribution of tandem repeats to the adaptation of the pathogen to its host. Tandem repeats appear to contribute to phenotypic variation in bacteria in at least two ways. Tandem repeats located within the regulatory region of a gene can constitute an on/off switch of gene expression at the transcriptional level [6,7]. Similarly, tandem repeats within coding regions with repeat units length not a multiple of three can induce a reversible premature end of translation when a mutation changes the number of repeats (reviewed in [8–10]). In other instances, the repeated unit length is a multiple of three, and the tandem repeat contributes to a coding region. In such cases, variations in the number of copies modify the gene product itself .
Mutation mechanisms of micro and minisatellites have been studied in some detail in eukaryotes, essentially human and yeast (reviewed in ). In brief, the data obtained so far suggest that microsatellites mutate by replication slippage processes; mutation rates depend upon the efficiency of mismatch repair mechanisms and an internal heterogeneity within the array strongly stabilizes the tandem repeat. In contrast, minisatellites mutate predominantly as the result of the repair of a double strand break initiated within, or very close to, the tandem repeat. In eukaryotes at least, these events can be of replicative origin , or can be genetically controlled, and specifically induced, during meiosis, at double strand breaks hot-spots. Minisatellite mutation rate in eukaryotes appears to be insensitive to mismatch repair efficiency, and internal heterogeneity is compatible with a high mutation rate [12, 14].
In bacteria, loci containing a tandem repeat from the microsatellite class (repeat unit sizes of 1-8 bp) have been called simple sequence contingency loci . Altered number of repeats allows for reversible on and off states of expression for the corresponding gene. The mutation rate of a tetranucleotide (microsatellite) tract in Haemophilus influenzae is higher than 10-4 and contributes to the adaptation of the pathogen to its hosts as the infection progresses . In such an extreme situation, the microsatellite is of limited value for strain identification, epidemiological and phylogenetic studies. The tandem repeat array is composed of perfect copies of the elementary unit, and different alleles are observed in a single culture. In contrast, the phylogenetic identity of minisatellite alleles of identical size can usually be further checked by DNA sequencing, since the repeated units are often not perfect . The pattern of variants along the array provides an additional level of allele identification and phylogenetic information. In addition, tandem repeats with longer repeat unit length can be relatively easily typed in the size range of a few hundred base-pairs using ordinary horizontal gel electrophoresis.
In this report, we will first describe the use of a tandem repeats database for bacterial genomes (http://minisatellites.u-psud.fr) and briefly compare the general characteristics of tandem repeats in a number of bacterial genomes for which the sequence has been determined and made publicly available. We will then show how this tool can easily be applied to the rapid characterization of new highly polymorphic markers in two pathogens, Y. pestis and B. anthracis.
Both Y. pestis (causative agent of plague) and B. anthracis (causative agent of anthrax) are recently emerged clones of respectively Y. pseudotuberculosis  and B. cereus . In the case of Y. pestis, a high resolution typing tool based on RFLP (Restriction Fragment Length Polymorphism) analysis of IS100 locations has already been developed . However this technology is more demanding than PCR typing, which justifies the development of such an assay. In the case of B. anthracis, polymorphisms were initially identified essentially using AFLP (Amplified Fragment Length Polymorphism) typing . Subsequent analyses demonstrated that the most informative fragments in AFLP patterns resulted from tandem repeat array length variations (five minisatellite loci were characterized in this way ).
Results and discussion
Use of the tandem repeats database
To date, 36 bacterial genome sequences from 32 species have been released in the public domain and are included in the database (Figure 1A; the nine archaebacteria genomes sequenced to date are presented in an other page, which can be accessed from http://minisatellites.u-psud.fr/). As many other sequencing projects are under way (http://www.ncbi.nlm.nih.gov/PMGifs/Genomes/bact.html; http://www.tigr.org/tdb/mdb/mdbinprogress.html; http://www.sanger.ac.uk/Projects/Microbes/), the database will be regularly updated. The collection of tandem repeats present in a given genome can be queried according to a combination of criteria, total tandem repeat array length (L), repeat unit length (U), number of repeats (N), percentage of conservation of the repeats along the array (V), position on the genome (Pos), average GC percent of the repeats (%GC), strand bias in nucleotide composition (B) (these values have been precomputed using the Tandem Repeats Finder software described in ). The results shown on Figure 1B use the "Tandem Repeats Distribution according to repeat unit length" option (Figure 1A). Three genomes were searched for tandem repeat arrays longer than 100 base-pairs (L ≥ 100). The genomes selected illustrate three different behaviors. On the right panel, Pseudomonas aeruginosa shows a very striking bias towards minisatellites with a motif length multiple of three. On the left and middle panels of Figure 1B, Buchnera sp and Y. pestis, show no such bias. The overall density of tandem repeat arrays longer than 100 base-pairs varies in the different genomes. Buchnera sp. contains 103 such loci, for a total genome size of 641 kb, which corresponds to a density per megabase of 161. Pseudomonas aeruginosa, with a total genome length of 6.3 Mb, has a density of 48. Y. pestis has an intermediate value of 30. Figure 2 summarizes the values observed in the 32 species. Ten non pathogenic species are presented in the upper part, 22 pathogenic species on the lower part. The species are ordered from top to bottom according to increasing genome size. The dark bars indicate for each genome the density per megabase of tandem repeat arrays longer than 100 bp. The clear bars reflect the excess of tandem repeats with unit length a multiple of three. A wide range of situations is observed, with a remarkable excess of tandem repeats multiples of three in Mycobacterium tuberculosis and Pseudomonas aeruginosa, presumably reflecting a significant contribution of tandem repeats to coding regions in these two bacteria.
As a quick illustration of the use of this database to facilitate the development of genotyping tools for bacterial genomes, we have evaluated the polymorphism associated with tandem repeats from Y. pestis on one hand and B. anthracis on the other (in this second instance, the genome sequence has not been completed yet and does not appear on the publicly accessible Tandem Repeats Database page, Figure 1A).
Application to Y. pestis
Figure 3A presents the result of a query run on Y. pestis, to identify tandem repeats with repeat units longer than 9 base-pairs repeated at least 7 times in the strain which has been sequenced (CO-92 biovar Orientalis).Sixty-four tandem repeats fulfill these criteria (an additional group of forty-nine have 6 copies of the motif; the twelve loci with the highest internal conservation were also included in this study). The output includes links to individual alignment files, as produced by the Tandem Repeat Finder software . The alignment file also includes 200 base-pairs of flanking sequence from each side of the tandem repeat, from which primers can be selected for PCR amplification. Figure 3B shows an annotated extract of one alignment file. The positions of the primers selected for subsequent PCR amplification are underlined. Three Y. pestis (representing the Antiqua, Medievalis, and Orientalis biovars ) and two Y. pseudotuberculosis strains were used for the initial identification of minisatellites sufficiently polymorphic to be of interest for further studies. Table 1 summarizes the PCR conditions used for each polymorphic locus and the results obtained. A total of 76 tandem repeats were tested. PCR amplification failed in 6 cases. Twenty one loci are monomorphic in the five Yersinia strains typed here. Forty-nine of the loci are polymorphic (Table 1). Twenty-five of these are polymorphic among the Y. pestis strains. Seven present a different allele in each of the five Yersinia strains, thirteen have a different allele in each of the three Y. pestis strains. Gel images for the 25 loci polymorphic among Y. pestis are shown in Figure 4. As can be seen, the repeat unit size and the overall length of the PCR products are such that tandem repeats differing by a single repeat unit can be distinguished by simple agarose gel electrophoresis.
Application to B. anthracis
Given the relatively low overall size of most bacterial tandem repeats, tandem repeat search can be run even on unfinished sequences. Tandem Repeats Finder was applied to B. anthracis sequence obtained from The Institute for Genomic Research through the website at http://www.tigr.org. The sequence was recovered as approximately 1000 contigs, for a total amount of slightly more than 5 Mb. Thirty tandem repeats have at least 10 copies of a repeat unit longer than 9 base-pairs. Fourteen of them are polymorphic among the 31 B. anthracis strains typed here (Table 2). Twenty-seven different genotypes are identified. Polymorphism information content (PIC) indexes based on the 27 genotypes vary from 0.07 to 0.82. Nine PIC values are above 0.5. Eight alleles are identified for CEB-Bams30, in a size range 270-900 base-pairs (Figure 5). In this case, the resolution of the largest alleles would probably be improved by using an automated DNA sequencer, and more alleles might be resolved. There are clear gaps in the size range coverage shown in Figure 5, and it is likely that the typing of additional strains would uncover new alleles. The genotyping data obtained was used to construct a phylogenetic tree based upon the Neighbor-Joining method (http://www.infobiogen.fr). In order to be able to correlate the tree obtained here with earlier studies , 5 minisatellites and one microsatellite reported previously were also typed. Figure 6 presents the data obtained and the resulting tree, using the nomenclature previously proposed . Six Bacillus cereus strains have also been included and used as an outgroup in the analysis. Occasionally B. cereus strains will not amplify (scored as 0 in Figure 6) or will give weak amplification signals (Figure 5, last six lanes on the right). The proposed tree is in good agreement with earlier results. In particular, the A and B clusters are well defined. We have apparently no representatives for the A1b and A3a group, whereas strains 9533 and 9502 to 9505 appear to define a new branch. The correspondence between allele numbering and allele size is indicated in Table 3.
Correlations between polymorphism and structural characteristics of minisatellites
We have looked for correlations between on one hand the number of alleles and polymorphism of the minisatellites, and on the other, simple structural characteristics of the tandem repeats in the sequenced strain : motif size, number of motifs, total length, conservation of the motifs along the array (percent identity), GC content, strand bias. In the case of B. anthracis, a highly significant correlation (0.01 level) is observed between polymorphism and both total length and GC content. This is not true for Y. pestis in which a strong correlation is seen between the number of alleles and the conservation of the motifs (Figure 7).
We limited here our investigation of tandem repeats to minisatellites, i.e. repeat units longer than 9 base-pairs, so as to avoid simple sequence contingency loci  of limited epidemiological value, and to facilitate the typing of alleles with agarose gel electrophoresis. However, simple sequence contingency loci are also represented in the database and are of great interest for molecular pathogenicity studies [6–8]. The use of the tandem repeats database was demonstrated here on two of the most genetically homogeneous human pathogenes, Y. pestis and B. anthracis. There is consequently a possibility that a common database format for identification and epidemiological analyses of pathogens amenable to minisatellite typing be developed. As more data becomes available on polymorphism associated with tandem repeats, it will be added to the database presented here in order to avoid duplication of work and nomenclature.
Bacterial species differ very significantly in the density of tandem repeats within their genome, and also in their use of tandem repeats. Some species have a very strong excess of tandem repeats with repeat units length which are multiple of three, the most striking examples being M. tuberculosis and P. aeruginosa. Polymorphism in such tandem repeats is likely to modulate the protein structure rather than gene activity. In M. tuberculosis, all tandem repeats with total length (L) higher than 100 bp and 9 or 15 base-pairs long units are located with ORFs . An important proportion of these tandem repeats correspond to the so-called PE and PPE multigene families .
In the two species studied here, tandem repeat polymorphism is strongly correlated with one or more of the sequenced allele characteristics, as illustrated in Figure 7. In Yersinia pestis a strong correlation is observed between number of alleles observed and homogeneity of the tandem array. In Bacillus anthracis, the strongest correlations are with total array length and GC content. It appears that the correlations are not the same in the two species, so that at present at least, the polymorphism associated with a tandem repeat cannot be inferred from its primary sequence. In particular, and in contrast to what is known for microsatellites (1-5 bp repeat units), some of the minisatellites are highly polymorphic in spite of a poor internal homogeneity of the sequenced allele, as is also the case for minisatellites in the human genome . However, more systematic allele sequencing will be required to demonstrate that polymorphism is not associated with a subclass of alleles showing a higher internal homogeneity. Similarly, allele sequencing will be required to formally establish that the allele size variations observed are indeed (as is likely) the consequence of variations in the number of repeats.
Five among the B. anthracis markers described here (Ceb-Bams1, 3, 7, 13 and 30) are highly polymorphic with PIC values (or Nei's index) above 0.7. In this respect, it is important to observe that the length of the allele observed for Ceb-Bams1 in the Ames strain is not of the size expected from the sequence data (Table 2). This may result either from a high mutation rate at Ceb-Bams1 or from a sequencing error. The expected allele size corresponds to allele 4 (Table 3), which is unlikely for the Ames strain because Ceb-Bams1 allele 4 is observed only in cluster B strains (Figure 6) and Ames is well apart of cluster B . A similar situation is observed for Ceb-Bams28, for which the expected product does not correspond to any existing allele in the collection of strains typed. In this case however, the locus is moderately polymorphic, with a PIC value of 0.26 and only three alleles observed (Table 2), so that a sequencing error is the most likely interpretation. This issue could be easily solved by typing with Ceb-Bams1 and Ceb-Bams28 the very strain which has been used for the sequencing project.
It is interesting to observe that, although the magnitude of allele size difference has not been taken into account when building the distance matrix, the resulting phylogenetic tree proposed in Figure 6 tends to group together strains with alleles of similar size at these most variable loci. This is reminiscent of observations made in H. influenzae  and suggest that mutation events are predominantly small size changes. Here again, more detailed studies involving full allele sequencing should now help understand the succession of events producing a population of alleles.
Materials and methods
Bacterial genomes DNA sequences
Finished sequences in the public domain were recovered by ftp from the NCBI or the Sanger center sites (http://www.ncbi.nlm.nih.gov/PMGifs/Genomes/bact.html; http://www.sanger.ac.uk/Projects/Microbes/). Preliminary sequence data for B. anthracis was obtained from The Institute for Genomic Research through the website at http://www.tigr.org.
All strains used here are part of the collection maintained by the Centre d'Etudes du Bouchet (CEB). They originate either from the CIP (Collection Institut Pasteur, http://www.pasteur.fr/) or from AFSSA (Agence Française de Sécurité Sanitaire des Aliments, http://www.afssa.fr/, Dr Josée Vaissaire). DNA from each isolate was obtained by large-batch procedures or by the simplified procedure as described in . In addition, 15 μg of DNA from the B. anthracis Ames strain were kindly provided by Dr Mats Forsman, FOA, Sweden.
Minisatellite PCR amplification and genotyping
PCR reactions were performed in 15 μl containing 1 ng of DNA, 1x Long Range Reaction Buffer 3 (Roche-Boehringer), 1 unit of Taq DNA polymerase, 200 μM of each dNTP, 0.3 μM of each flanking primer. The Taq DNA polymerase was either prepared essentially as described in  or purchased from Qbiogen or Roche-Boehringer. The 1x LongRange Buffer 3 is 1.75 mM MgCl2, 50 mM Tris-HCl pH9.2, 16 mM (NH4)2SO4.
PCR reactions were run on a Perkin-Elmer 9600 or a MJResearch PTC200 thermocycler. An initial denaturation at 96°C for five minutes was followed by 34 cycles of denaturation at 96°C for 20 seconds, annealing at 60°C for 30 seconds, elongation at 65°C for 1 minute, followed by a final extension step of 5 minutes at 65°C. In few cases, other annealing temperatures and/or elongation times were used (see tables 1 and 2). Five microliters of the PCR products where run on standard 1% or 2% agarose gel (Qbiogen) in 0.5 x TBE buffer at a voltage of 10 V/cm as indicated in Tables 1 and 2. Gel length of 10 to 40 cm were used according to PCR product size and motif length. Gels were stained with ethidium bromide and visualized under UV light. Allele sizes were estimated using as size markers the 1 kb ladder plus (Gibco-BRL which also includes a 100 bp ladder between 100 bp and 500 bp, plus 650, 850 and 1000 bp bands) or the 50 bp ladder (Euromedex) which provides a 50 bp ladder between 50 and 300 bp and a 100 bp ladder from 300 bp to 1000 bp.
Tandem Repeats Finder analysis
Sequences were processed using the Tandem Repeats Finder software (http://c3.biomath.mssm.edu/trf.html). The output was processed to eliminate duplicates before being imported in a database (running under Access2000, Microsoft Corp.) as described previously . The B. anthracis preliminary sequence data file uses FASTA type of headers (i.e. >sequenceId) to separate the independent contigs. The headers were replaced by runs of 10 Ns before running Tandem Repeats Finder.
Blast queries against the M. tuberculosis genome
The identifications of the open reading frames containing a given tandem repeat from M. tuberculosis were done by running a BLAST search on the dedicated web page at http://www.sanger.ac.uk/Projects/M_tuberculosis/blast_server.shtml.
Estimation of the excess of tandem repeats with motif length multiple of three
A χ2 test was calculated for the difference between the observed number of tandem repeats with motif length multiple of 3 and the expected number of tandem repeats with motif length multiple of 3 (expected value in the absence of bias being the total number of tandem repeats divided by 3). The χ2 values vary from 0.01 to 253.5. There is a significant excess (χ2 > 3.841) for all species but 6 (Buchnera sp, T. maritima, H. influenzae, M. genitalium, R. prowazekii, Y. pestis).
Polymorphism Information Index (PIC) or Nei's diversity index is calculated as 1 - Σ (allele frequency)2 based upon the unique genotypes.
A phenetic approach, based on a distance matrix was used. Distance matrix between strains was obtained by counting the number of differences between the corresponding genotypes. Then, Neighbor Joining cluster analysis was performed with Phylip  accessed at http://www.infobiogen.fr/. An outgroup was arbitrary chose among Bacillus cereus strains (9785) and input order of species was randomised.
Data (genotypes, distance matrix, phylogenetic tree) are available at http://minisatellites.u-psud.fr/ASPSamp/Phylogenie/data.htm
Correlations were calculated with the statistical program SPSS: Pearson correlation, and non-parametric correlations (Kendall's tau and Spearman's rho) show similar results.
van Belkum A, Scherer S, van Leeuwen W, Willemse D, van Alphen L, Verbrugh H: Variable number of tandem repeats in clinical strains of Haemophilus influenzae. Infect Immun. 1997, 65: 5017-27.
Keim P, Price LB, Klevytska AM, Smith KL, Schupp JM, Okinaka R, Jackson PJ, Hugh-Jones ME: Multiple-Locus Variable-Number Tandem Repeat Analysis Reveals Genetic Relationships within Bacillus anthracis. J Bacteriol. 2000, 182: 2928-2936. 10.1128/JB.182.10.2928-2936.2000.
Frothingham R, Meeker-O'Connell WA: Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. Microbiology. 1998, 144: 1189-96.
Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, Locht C: Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol. 2000, 36: 762-71. 10.1046/j.1365-2958.2000.01905.x.
Adair DM, Worsham PL, Hill KK, Klevytska AM, Jackson PJ, Friedlander AM, Keim P: Diversity in a variable-number tandem repeat from Yersinia pestis. J Clin Microbiol. 2000, 38: 1516-9.
van Ham SM, van Alphen L, Mooi FR, van Putten JP: Phase variation of H. influenzae fimbriae: transcriptional control of two divergent genes through a variable combined promoter region. Cell. 1993, 73: 1187-96.
Weiser JN, Love JM, Moxon ER: The molecular mechanism of phase variation of H. influenzae lipopolysaccharide. Cell. 1989, 59: 657-65.
Bayliss CD, Field D, Moxon ER: The simple sequence contingency loci of Haemophilus influenzae and Neisseria meningitidis. J Clin Invest. 2001, 107: 657-66.
Henderson IR, Owen P, Nataro JP: Molecular switches - the ON and OFF of bacterial phase variation. Mol Microbiol. 1999, 33: 919-32. 10.1046/j.1365-2958.1999.01555.x.
Wang G, Ge Z, Rasko DA, Taylor DE: Lewis antigens in Helicobacter pylori: biosynthesis and phase variation. Mol Microbiol. 2000, 36: 1187-96. 10.1046/j.1365-2958.2000.01934.x.
Wilton JL, Scarman AL, Walker MJ, Djordjevic SP: Reiterated repeat region variability in the ciliary adhesin gene of Mycoplasma hyopneumoniae. Microbiology. 1998, 144: 1931-43.
Vergnaud G, Denoeud F: Minisatellites: Mutability and Genome Architecture. Genome Res. 2000, 10: 899-907. 10.1101/gr.10.7.899.
Kokoska RJ, Stefanovic L, Tran HT, Resnick MA, Gordenin DA, Petes TD: Destabilization of yeast micro- and minisatellite DNA sequences by mutations affecting a nuclease involved in Okazaki fragment processing (rad27) and DNA polymerase delta (pol3-t). Mol Cell Biol. 1998, 18: 2779-88.
Debrauwère H, Buard J, Tessier J, Aubert D, Vergnaud G, Nicolas A: Meiotic instability of human minisatellite CEB1 in yeast requires DNA double-strand breaks. Nat Genet. 1999, 23: 367-71. 10.1038/15557.
De Bolle X, Bayliss CD, Field D, van de Ven T, Saunders NJ, Hood DW, Moxon ER: The length of a tetranucleotide repeat tract in Haemophilus influenzae determines the phase variation rate of a gene with homology to type III DNA methyltransferases. Mol Microbiol. 2000, 35: 211-22. 10.1046/j.1365-2958.2000.01701.x.
van Belkum A, Scherer S, van Alphen L, Verbrugh H: Short-sequence DNA repeats in prokaryotic genomes. Microbiol Mol Biol Rev. 1998, 62: 275-93.
Achtman M, Zurth K, Morelli G, Torrea G, Guiyoule A, Carniel E: Yersinia pestis, the cause of plague, is a recently emerged clone of Yersinia pseudotuberculosis [published erratum appears in Proc Natl Acad Sci U S A 2000 Jul 5;97(14):8192]. Proc Natl Acad Sci U S A. 1999, 96: 14043-8. 10.1073/pnas.96.24.14043.
Helgason E, Okstad OA, Caugant DA, Johansen HA, Fouet A, Mock M, Hegna I, Kolsto : Bacillus anthracis, Bacillus cereus, and Bacillus thuringiensis - one species on the basis of genetic evidence. Appl Environ Microbiol. 2000, 66: 2627-30. 10.1128/AEM.66.6.2627-2630.2000.
Keim P, Kalif A, Schupp J, Hill K, Travis SE, Richmond K, Adair DM, Hugh-Jones M, Kuske CR, Jackson P: Molecular evolution and diversity in Bacillus anthracis as detected by amplified fragment length polymorphism markers. J Bacteriol. 1997, 179: 818-24.
Benson G: Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 1999, 27: 573-80. 10.1093/nar/27.2.573.
Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, Gordon SV, Eiglmeier K, Gas S, Barry CE, Tekaia F, Badcock K, Basham D, Brown D, Chillingworth T, Connor R, Davies R, Devlin K, Feltwell T, Gentles S, Hamlin N, Holroyd S, Hornsby T, Jagels K, Barrell BG: Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998, 393: 537-44. 10.1038/31159.
Engelke DR, Krikos A, Bruck ME, Ginsburg D: Purification of Thermus aquaticus DNA polymerase expressed in Escherichia coli. Anal. Biochem. 1990, 191: 396-400.
Felsenstein J: PHYLIP - Phylogeny Inference Package (Version 3.2). Cladistics. 1989, 5: 164-166.
Minisatellite investigations in the laboratory are supported by grants from Délégation Générale de l'Armement (DGA/DSA/STTC and DGA/DSA/SPNuc). Preliminary sequence data for B. anthracis was obtained from The Institute for Genomic Research through the website at http://www.tigr.org. Sequencing of B. anthracis was accomplished with support from Office of Naval Research, Department of Energy, and National Institute of Allergy and Infectious diseases. The Yersinia pestis sequence data were produced by the Y. pestis Sequencing Group at the Sanger Centre and can be obtained from ftp://ftp.sanger.ac.uk/pub/pathogens/yp/. We wish to thank the referees for the significant improvements they have suggested.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.