Search

Brucella abortus is the main causative agent of brucellosis in cattle, leading to severe economic consequences in agriculture and affecting public health. The zoonotic nature of the infection increases the need to control the spread and dynamics of outbreaks in animals with the incorporation of high resolution genotyping techniques. Based on such methods, B. abortus is currently divided into three clades, A, B, and C. The latter includes subclades C1 and C2. This study presents the results of whole-genome sequencing of 49 B. abortus strains isolated in Kazakhstan between 1947 and 2015 and of 36 B. abortus strains of various geographic origins isolated from 1940 to 2004. In silico Multiple Locus Sequence Typing (MLST) allowed to assign strains from Kazakhstan to subclades C1 and to a much lower extend C2. Whole-genome Single-Nucleotide Polymorphism (wgSNP) analysis of the 46 strains of subclade C1 with strains of worldwide origins showed clustering with strains from neighboring countries, mostly North Caucasia, Western Russia, but also Siberia, China, and Mongolia. One of the three Kazakhstan strains assigned to subclade C2 matched the B. abortus S19 vaccine strain used in cattle, the other two were genetically close to the 104 M vaccine strain. Bayesian phylodynamic analysis dated the introduction of B. abortus subclade C1 into Kazakhstan to the 19th and early 20th centuries. We discuss this observation in view of the history of population migrations from Russia to the Kazakhstan steppes.

This article describes Bacillus anthracis strains isolated in Kazakhstan since the 1950s until year 2016 from sixty-one independent events associated with anthrax in humans and animals. One hundred and fifty-four strains were first genotyped by Multiple Locus VNTR (variable number of tandem repeats) Analysis (MLVA) using 31 VNTR loci. Thirty-five MLVA31 genotypes were resolved, 28 belong to the A1/TEA group, five to A3/Sterne-Ames group, one to A4/Vollum and one to the B clade. This is the first report of the presence of the B-clade in Kazakhstan. The MLVA31 results and epidemiological data were combined to select a subset of seventy-nine representative strains for draft whole genome sequencing (WGS). Strains from Kazakhstan significantly enrich the known phylogeny of the Ames group polytomy, including the description of a new branch closest to the Texas, United States A.Br.Ames sublineage stricto sensu. Three among the seven currently defined branches in the TEA polytomy are present in Kazakhstan, “Tsiankovskii”, “Heroin”, and “Sanitary Technical Institute (STI)”. In particular, strains from the STI lineage are largely predominant in Kazakhstan and introduce numerous deep branching STI sublineages, demonstrating a high geographic correspondence between “STI” and Kazakhstan, Central Asia. This observation is a strong indication that the TEA polytomy emerged after the last political unification of Asian steppes in the fourteenth century of the Common Era. The phylogenetic analysis of the Kazakhstan data and of currently available WGS data of worldwide origin strengthens our understanding of B. anthracis geographic expansions in the past seven centuries.

ABSTRACT Q (query) fever is an infectious zoonotic disease caused by the Gram-negative bacterium Coxiella burnetii. Although the disease has been studied for decades, it still represents a threat due to sporadic outbreaks across farms in Europe. The absence of a central platform for Coxiella typing data management is an important epidemiological gap that is relevant in the case of an outbreak. To fill this gap, we have designed and implemented an online, open-source, web-based platform called CoxBase (https://coxbase.q-gaps.de). This platform includes a database that holds genotyping information on more than 400 Coxiella isolates alongside metadata that annotate them. We have also implemented features for in silico genotyping of completely or minimally assembled Coxiella sequences using five different typing methods, querying of existing isolates, visualization of isolate geodata via aggregation on a world map, and submission of new isolates. We tested our in silico typing method on 50 Coxiella genomes downloaded from the RefSeq database, and we successfully genotyped all genomes except for cases where the sequence quality was poor. We identified new spacer sequences using our implementation of the multispacer sequence typing (MST) in silico typing method and established adaA gene phenotypes for all 50 genomes as well as their plasmid types. IMPORTANCE Q fever is a zoonotic disease that is a source of active epidemiological concern due to its persistent threat to public health. In this project, we have identified areas in the field of Coxiella research, especially regarding public health and genomic analysis, where there is an inadequacy of resources to monitor, organize, and analyze genomic data from C. burnetii. Subsequently, we have created an open, web-based platform that contains epidemiological information, genome typing functions comprising all the available Coxiella typing methods, and tools for isolate data discovery and visualization that could help address the above-mentioned challenges. This is the first platform to combine all disparate genotyping systems for Coxiella burnetii as well as metadata assets with tools for genomic comparison and analyses. This platform is a valuable resource for laboratory researchers as well as research epidemiologists interested in investigating the relatedness or dissimilarity among C. burnetii strains.

The genus Brucella comprises major pathogenic species causing disease in livestock and humans, e.g. B. melitensis. In the past few years, the genus has been significantly expanded by the discovery of phylogenetically more distant lineages comprising strains from diverse wildlife animal species, including amphibians and fish. The strains represent several potential new species, with B. inopinata as solely named representative. Being genetically more distant between each other, relative to the “classical” Brucella species, they present distinct atypical phenotypes and surface antigens. Among surface protein antigens, the Omp2a and Omp2b porins display the highest diversity in the classical Brucella species. The genes coding for these proteins are closely linked in the Brucella genome and oriented in opposite directions. They share between 85 and 100% sequence identity depending on the Brucella species, biovar, or genotype. Only the omp2b gene copy has been shown to be expressed and genetic variation is extensively generated by gene conversion between the two copies. In this study, we analyzed the omp2 loci of the non-classical Brucella spp. Starting from two distinct ancestral genes, represented by Australian rodent strains and B. inopinata, a stepwise nucleotide reduction was observed in the omp2b gene copy. It consisted of a first reduction affecting the region encoding the surface L5 loop of the porin, previously shown to be critical in sugar permeability, followed by a nucleotide reduction in the surface L8 loop-encoding region. It resulted in a final omp2b gene size shared between two distinct clades of non-classical Brucella spp. (African bullfrog isolates) and the group of classical Brucella species. Further evolution led to complete homogenization of both omp2 gene copies in some Brucella species such as B. vulpis or B. papionis. The stepwise omp2b deletions seemed to be generated through recombination with the respective omp2a gene copy, presenting a conserved size among Brucella spp., and may involve short direct DNA repeats. Successive Omp2b porin alteration correlated with increasing porin permeability in the course of evolution of Brucella spp. They possibly have adapted their porin to survive environmental conditions encountered and to reach their final status as intracellular pathogen.

We describe the genetic diversity of 1327 Brucella strains from human patients in Kazakhstan using multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA). All strains were assigned to the Brucella melitensis East Mediterranean group and clustered into 16 MLVA11 genotypes, nine of which are reported for the first time. MLVA11 genotype 116 predominates (86.8%) and is present all over Kazakhstan indicating existence and temporary preservation of a “founder effect” among B. melitensis strains circulating in Central Eurasia. The diversity pattern observed in humans is highly similar to the pattern previously reported in animals. The diversity observed by MLVA suggested that the epidemiological status of brucellosis in Kazakhstan is the result of the introduction of a few lineages, which have subsequently diversified at the most unstable tandem repeat loci. This investigation will allow to select the most relevant strains for testing these hypotheses via whole genome sequencing and to subsequently adjust the genotyping scheme to the Kazakhstan epidemiological situation.

Previous studies have shown the usefulness of MLVA16 as a rapid molecular identification and classification method for Brucella species and biovars including recently described novel Brucella species from wildlife. Most studies were conducted on a limited number of strains from limited geographic/host origins. The objective of this study was to assess genetic diversity of Brucella spp. by MLVA16 on a larger scale. Thus, 1404 animal or human isolates collected from all parts of the world over a period of 32 years (1974-2006) were investigated. Selection of the 1404 strains was done among the approximately 4000 strains collection of the BCCN (Brucella Culture Collection Nouzilly), based on classical biotyping and on the animal/human/geographic origin over the time period considered. MLVA16 was performed on extracted DNAs using high throughput capillary electrophoresis. The 16 loci were amplified in four multiplex PCR reactions. This large scale study firstly confirmed the accuracy of MLVA16 typing for Brucella species and biovar identification and its congruence with the recently described Extended Multilocus Sequence Analysis. In addition, it allowed identifying novel MLVA11 (based upon 11 slowly evolving VNTRs) genotypes representing an increase of 15% relative to the previously known Brucella MLVA11 genotypes. Cluster analysis showed that among the MLVA16 genotypes some were genetically more distant from the major classical clades. For example new major clusters of B. abortus biovar 3 isolated from cattle in Sub-Saharan Africa were identified. For other classical species and biovars this study indicated also genotypic expansion within the population structure of classical Brucella species. MLVA proves to be a powerful tool to rapidly assess genetic diversity of bacterial populations on a large scale, as here on a large collection of strains of the genomically homogeneous genus Brucella. The highly discriminatory power of MLVA appears of particular interest as a first step for selection of Brucella strains for whole-genome sequencing. The MLVA data of this study were added to the public Brucella MLVA database at http://microbesgenotyping.i2bc.paris-saclay.fr. Current version Brucella_4_3 comprises typing data from more than 5000 strains including in silico data analysis of public whole genome sequence datasets.

Transposable phages, also called saltoviruses, of which the Escherichia coli phage Mu is the reference, are temperate phages that multiply their genome through replicative transposition at multiple sites in their host chromosome. The viral genome is packaged together with host DNA at both ends. In the present work, genome sequencing of three Pseudomonas aeruginosa transposable phages, HW12, 2P1, and Ab30, incidentally gave us access to the location of thousands of replicative integration sites and revealed the existence of a variable number of hotspots. Taking advantage of deep sequencing, we then designed an experiment to study 13,000,000 transposon integration sites of bacteriophage Ab30. The investigation revealed the presence of 42 transposition hotspots adjacent to bacterial interspersed mosaic elements (BIME) accounting for 5% of all transposition sites. The rest of the sites appeared widely distributed with the exception of coldspots associated with low G-C content segments, including the putative O-antigen biosynthesis cluster. Surprisingly, 0.4% of the transposition events occurred in a copy of the phage genome itself, indicating that the previously described immunity against such events is slightly leaky. This observation allowed drawing an image of the phage chromosome supercoiling into four loops.

Abstract Background Staphylococcus aureus infection in patients with cystic fibrosis (CF) is frequent and may be due to colonization by a few pathogenic lineages. Systematic genotyping of all isolates, methicillin-susceptible S. aureus (MSSA) as well as methicillin-resistant S. aureus (MRSA) is necessary to identify such lineages and follow their evolution in patients. Multiple-locus variable-number tandem repeat analysis (MLVA/VNTR) was used to survey S. aureus clinical isolates in a French paediatric CF centre. Results During a 30 months period, 108 patients, aged 2 to 21 years, regularly followed up at the centre, provided sputum for culture. From 79 patients, a total of 278 isolates were genotyped by MLVA, resolving into 110 genotypes and 19 clonal complexes (CC) composed of similar or closely related isolates. 71% of the strains were distributed into four main CCs, in term of number of isolates and number of genotypes. Spa (Staphylococcus protein A) typing was performed on representative samples, showing an excellent concordance with MLVA. In 17 patients, strains from two to four different CCs were recovered over time. On six occasions, S. aureus isolates with the same genotype were shared by 2 different patients and they belonged to one of the four main clusters. Methicillin-resistance was observed in 60% of the isolates, 90% of which belonged to the main clonal complexes CC8, CC45 and CC5. In 5 patients, methicillin-resistance of S. aureus isolates was not associated with the mecA gene: for four patients, it was due to overproduction of β-lactamase, leading to BOR-SA (borderline S. aureus) isolates, while a strain showing probably a new modified penicillin-binding capacity (MOD-SA) was observed from one patient. Conclusion Systematic genotyping of S. aureus isolates recovered from sputum of CF children allows a thorough analysis of the strains responsible for sporadic as well as chronic colonization and the follow up of their evolution over time. We show here that more than 70% of these strains belong to 4 major CCs. MSSA as well as MRSA, BOR-SA and MOD-SA isolates can persist over several years, despite antibiotic treatments.

Abstract Background Since 1994, Brucella strains have been isolated from a wide range of marine mammals. They are currently recognized as two new Brucella species, B. pinnipedialis for the pinniped isolates and B. ceti for the cetacean isolates in agreement with host preference and specific phenotypic and molecular markers. In order to investigate the genetic relationships within the marine mammal Brucella isolates and with reference to terrestrial mammal Brucella isolates, we applied in this study the Multiple Loci VNTR (Variable Number of Tandem Repeats) Analysis (MLVA) approach. A previously published assay comprising 16 loci (MLVA-16) that has been shown to be highly relevant and efficient for typing and clustering Brucella strains from animal and human origin was used. Results 294 marine mammal Brucella strains collected in European waters from 173 animals and a human isolate from New Zealand presumably from marine origin were investigated by MLVA-16. Marine mammal Brucella isolates were shown to be different from the recognized terrestrial mammal Brucella species and biovars and corresponded to 3 major related groups, one specific of the B. ceti strains, one of the B. pinnipedialis strains and the last composed of the human isolate. In the B. ceti group, 3 subclusters were identified, distinguishing a cluster of dolphin, minke whale and porpoise isolates and two clusters mostly composed of dolphin isolates. These results were in accordance with published analyses using other phenotypic or molecular approaches, or different panels of VNTR loci. The B. pinnipedialis group could be similarly subdivided in 3 subclusters, one composed exclusively of isolates from hooded seals (Cystophora cristata) and the two others comprising other seal species isolates. Conclusion The clustering analysis of a large collection of marine mammal Brucella isolates from European waters significantly strengthens the current view of the population structure of these two species, and their relative position with respect to the rest of the Brucella genus. MLVA-16 is confirmed as being a rapid, highly discriminatory and reproducible method to classify Brucella strains including the marine mammal isolates. The Brucella2009 MLVA-16 genotyping database available at http://mlva.u-psud.fr/ is providing a detailed coverage of all 9 currently recognized Brucella species.

Abstract Background Anthrax and plague are diseases caused by Bacillus anthracis and Yersinia pestis respectively. These bacteria are etiological agents for worldwide zoonotic diseases and are considered among the most feared potential bioterror agents. Strain differentiation is difficult for these microorganisms because of their high intraspecies genome homogeneity. Moreover, fast strain identification and comparison with known genotypes may be crucial for naturally occurring outbreaks versus bioterrorist events discrimination. Results Thirty-nine B. anthracis and ten Y. pestis strains, representative of the species genetic diversity, were genotyped by Agilent 2100 Bioanalyzer using previously described Multiple Locus VNTR Analysis assays (MLVA). Results were compared to previous data obtained by standard genotyping system (capillary electrophoresis on automatic sequencer) and, when necessary, direct amplicon sequencing. A reference comparison table containing actual fragment sizes, sequencer sizes and Agilent sizes was produced. Conclusion In this report an automated DNA electrophoresis apparatus which provides a cheaper alternative compared to capillary electrophoresis approaches was applied for genotyping of B. anthracis and Y. pestis. This equipment, uses pre-cast gels and provides easy transportation, low maintenance and overall general logistic requirements and costs, is easy to set up and provides rapid analysis. This platform is a candidate for on-site MLVA genotyping of biothreat agents as well as other bacterial pathogens. It is an alternative to the more expensive and demanding capillary electrophoresis methods, and to the less expensive but more time-consuming classical gel electrophoresis approach.

Abstract Background The genome of Bacillus anthracis, the etiological agent of anthrax, is highly monomorphic which makes differentiation between strains difficult. A Multiple Locus Variable-number tandem repeats (VNTR) Analysis (MLVA) assay based on 20 markers was previously described. It has considerable discrimination power, reproducibility, and low cost, especially since the markers proposed can be typed by agarose-gel electrophoresis. However in an emergency situation, faster genotyping and access to representative databases is necessary. Results Genotyping of B. anthracis reference strains and isolates from France and Italy was done using a 25 loci MLVA assay combining 21 previously described loci and 4 new ones. DNA was amplified in 4 multiplex PCR reactions and the length of the resulting 25 amplicons was estimated by automated capillary electrophoresis. The results were reproducible and the data were consistent with other gel based methods once differences in mobility patterns were taken into account. Some alleles previously unresolved by agarose gel electrophoresis could be resolved by capillary electrophoresis, thus further increasing the assay resolution. One particular locus, Bams30, is the result of a recombination between a 27 bp tandem repeat and a 9 bp tandem repeat. The analysis of the array illustrates the evolution process of tandem repeats. Conclusion In a crisis situation of suspected bioterrorism, standardization, speed and accuracy, together with the availability of reference typing data are important issues, as illustrated by the 2001 anthrax letters event. In this report we describe an upgrade of the previously published MLVA method for genotyping of B. anthracis and apply the method to the typing of French and Italian B. anthracis strain collections. The increased number of markers studied compared to reports using only 8 loci greatly improves the discrimination power of the technique. An Italian strain belonging to the B branch was described, and two new branches, D and E, are proposed. Owing to the upgrading achieved here, precise genotyping can now be produced either by automated capillary electrophoresis, or by the more accessible but slower and for some markers slightly less accurate agarose gel methodology.

Abstract Background Coxiella burnetii, the causative agent of Q fever, has a wide host range. Few epidemiological tools are available, and they are often expensive or not easily standardized across laboratories. In this work, C. burnetii isolates from livestock and ticks were typed using infrequent restriction site-PCR (IRS-PCR) and multiple loci variable number of tandem repeats (VNTR) analysis (MLVA). Results By applying IRS-PCR, 14 C. burnetii isolates could be divided into six groups containing up to five different isolates. Clustering as deduced from MLVA typing with 17 markers provided an increased resolution with an excellent agreement to IRS-PCR, and with the plasmid type of each strain. MLVA was then applied to 28 additional C. burnetii isolates of different origin and 36 different genotypes were identified among the 42 isolates investigated. The clustering obtained is in agreement with published Multiple Locus Sequence Typing (MLST) data. Two panels of markers are proposed, panel 1 which can be confidently typed on agarose gel at a lower cost and in any laboratory setting (10 minisatellite markers with a repeat unit larger than 9 bp), and panel 2 which comprises 7 microsatellites and provides a higher discriminatory power. Conclusion Our analyses demonstrate that MLVA is a powerful and promising molecular typing tool with a high resolution and of low costs. The consistency of the results with independent methods suggests that MLVA can be applied for epidemiological studies. The resulting data can be queried on a dedicated MLVA genotyping Web service.

Abstract Background The classification of Brucella into species and biovars relies on phenotypic characteristics and sometimes raises difficulties in the interpretation of the results due to an absence of standardization of the typing reagents. In addition, the resolution of this biotyping is moderate and requires the manipulation of the living agent. More efficient DNA-based methods are needed, and this work explores the suitability of multiple locus variable number tandem repeats analysis (MLVA) for both typing and species identification. Results Eighty tandem repeat loci predicted to be polymorphic by genome sequence analysis of three available Brucella genome sequences were tested for polymorphism by genotyping 21 Brucella strains (18 reference strains representing the six 'classical' species and all biovars as well as 3 marine mammal strains currently recognized as members of two new species). The MLVA data efficiently cluster the strains as expected according to their species and biovar. For practical use, a subset of 15 loci preserving this clustering was selected and applied to the typing of 236 isolates. Using this MLVA-15 assay, the clusters generated correspond to the classical biotyping scheme of Brucella spp. The 15 markers have been divided into two groups, one comprising 8 user-friendly minisatellite markers with a good species identification capability (panel 1) and another complementary group of 7 microsatellite markers with higher discriminatory power (panel 2). Conclusion The MLVA-15 assay can be applied to large collections of Brucella strains with automated or manual procedures, and can be proposed as a complement, or even a substitute, of classical biotyping methods. This is facilitated by the fact that MLVA is based on non-infectious material (DNA) whereas the biotyping procedure itself requires the manipulation of the living agent. The data produced can be queried on a dedicated MLVA web service site.

Abstract Background Precise identification of bacterial pathogens at the strain level is essential for epidemiological purposes. In Streptococcus pneumoniae, the existence of 90 different serotypes makes the typing particularly difficult and requires the use of highly informative tools. Available methods are relatively expensive and cannot be used for large-scale or routine typing of any new isolate. We explore here the potential of MLVA (Multiple Loci VNTR Analysis; VNTR, Variable Number of Tandem Repeats), a method of growing importance in the field of molecular epidemiology, for genotyping of Streptococcus pneumoniae. Results Available genome sequences were searched for polymorphic tandem repeats. The loci identified were typed across a collection of 56 diverse isolates and including a group of serotype 1 isolates from Africa. Eventually a set of 16 VNTRs was proposed for MLVA-typing of S. pneumoniae. These robust markers were sufficient to discriminate 49 genotypes and to aggregate strains on the basis of the serotype and geographical origin, although some exceptions were found. Such exceptions may reflect serotype switching or horizontal transfer of genetic material. Conclusion We describe a simple PCR-based MLVA genotyping scheme for S. pneumoniae which may prove to be a powerful complement to existing tools for epidemiological studies. Using this technique we uncovered a clonal population of strains, responsible for infections in Burkina Faso. We believe that the proposed MLVA typing scheme can become a standard for epidemiological studies of S. pneumoniae.

Abstract Background Yersinia pestis, the agent of plague, is a young and highly monomorphic species. Three biovars, each one thought to be associated with the last three Y. pestis pandemics, have been defined based on biochemical assays. More recently, DNA based assays, including DNA sequencing, IS typing, DNA arrays, have significantly improved current knowledge on the origin and phylogenetic evolution of Y. pestis. However, these methods suffer either from a lack of resolution or from the difficulty to compare data. Variable number of tandem repeats (VNTRs) provides valuable polymorphic markers for genotyping and performing phylogenetic analyses in a growing number of pathogens and have given promising results for Y. pestis as well. Results In this study we have genotyped 180 Y. pestis isolates by multiple locus VNTR analysis (MLVA) using 25 markers. Sixty-one different genotypes were observed. The three biovars were distributed into three main branches, with some exceptions. In particular, the Medievalis phenotype is clearly heterogeneous, resulting from different mutation events in the napA gene. Antiqua strains from Asia appear to hold a central position compared to Antiqua strains from Africa. A subset of 7 markers is proposed for the quick comparison of a new strain with the collection typed here. This can be easily achieved using a Web-based facility, specifically set-up for running such identifications. Conclusion Tandem-repeat typing may prove to be a powerful complement to the existing phylogenetic tools for Y. pestis. Typing can be achieved quickly at a low cost in terms of consumables, technical expertise and equipment. The resulting data can be easily compared between different laboratories. The number and selection of markers will eventually depend upon the type and aim of investigations.

Abstract Background Currently available reference methods for the molecular epidemiology of the Mycobacterium tuberculosis complex either lack sensitivity or are still too tedious and slow for routine application. Recently, tandem repeat typing has emerged as a potential alternative. This report contributes to the development of tandem repeat typing for M. tuberculosis by summarising the existing data, developing additional markers, and setting up a freely accessible, fast, and easy to use, internet-based service for strain identification. Results A collection of 21 VNTRs incorporating 13 previously described loci and 8 newly evaluated markers was used to genotype 90 strains from the M. tuberculosis complex (M. tuberculosis (64 strains), M. bovis (9 strains including 4 BCG representatives), M. africanum (17 strains)). Eighty-four different genotypes are defined. Clustering analysis shows that the M. africanum strains fall into three main groups, one of which is closer to the M. tuberculosis strains, and an other one is closer to the M. bovis strains. The resulting data has been made freely accessible over the internet http://bacterial-genotyping.igmors.u-psud.fr/bnserver to allow direct strain identification queries. Conclusions Tandem-repeat typing is a PCR-based assay which may prove to be a powerful complement to the existing epidemiological tools for the M. tuberculosis complex. The number of markers to type depends on the identification precision which is required, so that identification can be achieved quickly at low cost in terms of consumables, technical expertise and equipment.

Abstract Background 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. Results 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. Conclusions 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.