Anthrax, caused by the bacterium Bacillus anthracis, is a disease of historical and current importance that is found throughout the world. The basis of its historical transmission is anecdotal and its true global population structure has remained largely cryptic. Seven diverse B. anthracis strains were whole-genome sequenced to identify rare single nucleotide polymorphisms (SNPs), followed by phylogenetic reconstruction of these characters onto an evolutionary model. This analysis identified SNPs that define the major clonal lineages within the species. These SNPs, in concert with 15 variable number tandem repeat (VNTR) markers, were used to subtype a collection of 1,033 B. anthracis isolates from 42 countries to create an extensive genotype data set. These analyses subdivided the isolates into three previously recognized major lineages (A, B, and C), with further subdivision into 12 clonal sub-lineages or sub-groups and, finally, 221 unique MLVA15 genotypes. This rare genomic variation was used to document the evolutionary progression of B. anthracis and to establish global patterns of diversity. Isolates in the A lineage are widely dispersed globally, whereas the B and C lineages occur on more restricted spatial scales. Molecular clock models based upon genome-wide synonymous substitutions indicate there was a massive radiation of the A lineage that occurred in the mid-Holocene (3,064-6,127 ybp). On more recent temporal scales, the global population structure of B. anthracis reflects colonial-era importation of specific genotypes from the Old World into the New World, as well as the repeated industrial importation of diverse genotypes into developed countries via spore-contaminated animal products. These findings indicate humans have played an important role in the evolution of anthrax by increasing the proliferation and dispersal of this now global disease. Finally, the value of global genotypic analysis for investigating bioterrorist-mediated outbreaks of anthrax is demonstrated.

Background The facultative, intracellular bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a serious infectious disease of humans and animals. We identified and categorized tandem repeat arrays and their distribution throughout the genome of B . pseudomallei strain K96243 in order to develop a genetic typing method for B. pseudomallei . We then screened 104 of the potentially polymorphic loci across a diverse panel of 31 isolates including B. pseudomallei , B. mallei and B. thailandensis in order to identify loci with varying degrees of polymorphism. A subset of these tandem repeat arrays were subsequently developed into a multiple-locus VNTR analysis to examine 66 B. pseudomallei and 21 B. mallei isolates from around the world, as well as 95 lineages from a serial transfer experiment encompassing ~18,000 generations. Results B. pseudomallei contains a preponderance of tandem repeat loci throughout its genome, many of which are duplicated elsewhere in the genome. The majority of these loci are composed of repeat motif lengths of 6 to 9 bp with 4 to 10 repeat units and are predominately located in intergenic regions of the genome. Across geographically diverse B . pseudomallei and B . mallei isolates, the 32 VNTR loci displayed between 7 and 28 alleles, with Nei's diversity values ranging from 0.47 and 0.94. Mutation rates for these loci are comparable (>10 -5 per locus per generation) to that of the most diverse tandemly repeated regions found in other less diverse bacteria. Conclusion The frequency, location and duplicate nature of tandemly repeated regions within the B. pseudomallei genome indicate that these tandem repeat regions may play a role in generating and maintaining adaptive genomic variation. Multiple-locus VNTR analysis revealed extensive diversity within the global isolate set containing B. pseudomallei and B. mallei , and it detected genotypic differences within clonal lineages of both species that were identical using previous typing methods. Given the health threat to humans and livestock and the potential for B. pseudomallei to be released intentionally, MLVA could prove to be an important tool for fine-scale epidemiological or forensic tracking of this increasingly important environmental pathogen.

VNTRs are able to discriminate among closely related isolates of recently emerged clonal pathogens, including Yersinia pestis the etiologic agent of plague, because of their great diversity. Diversity is driven largely by mutation but little is known about VNTR mutation rates, factors affecting mutation rates, or the mutational mechanisms. The molecular epidemiological utility of VNTRs will be greatly enhanced when this foundational knowledge is available. Here, we measure mutation rates for 43 VNTR loci in Y. pestis using an in vitro generated population encompassing approximately 96,000 generations. We estimate the combined 43-locus rate and individual rates for 14 loci. A comparison of Y. pestis and Escherichia coli O157:H7 VNTR mutation rates and products revealed a similar relationship between diversity and mutation rate in these two species. Likewise, the relationship between repeat copy number and mutation rate is nearly identical between these species, suggesting a generalized relationship that may be applicable to other species. The single- versus multiple-repeat mutation ratios and the insertion versus deletion mutation ratios were also similar, providing support for a general model for the mutations associated with VNTRs. Finally, we use two small sets of Y. pestis isolates to show how this general model and our estimated mutation rates can be used to compare alternate phylogenies, and to evaluate the significance of genotype matches, near-matches, and mismatches found in empirical comparisons with a reference database.

In July 1993, a liquid suspension of Bacillus anthracis was aerosolized from the roof of an eight-story building in Kameido, Tokyo, Japan, by the religious group Aum Shinrikyo. During 1999 to 2001, microbiologic tests were conducted on a liquid environmental sample originally collected during the 1993 incident. Nonencapsulated isolates of B. anthracis were cultured from the liquid. Multiple-locus, variable-number tandem repeat analysis found all isolates to be identical to a strain used in Japan to vaccinate animals against anthrax, which was consistent with the Aum Shinrikyo members’ testimony about the strain source. In 1999, a retrospective case-detection survey was conducted to identify potential human anthrax cases associated with the incident, but none were found. The use of an attenuated B. anthracis strain, low spore concentrations, ineffective dispersal, a clogged spray device, and inactivation of the spores by sunlight are all likely contributing factors to the lack of human cases.

We investigated the self-splicing properties of two introns from the bacterium Bacillus anthracis. One intron (B.a.I1) splices poorly in vitro despite having typical structural motifs, while the second (B.a.I2) splices well while having apparently degenerated features. The spliced exons of B.a.I2 were sequenced, and splicing was found to occur at a 3' site shifted one nucleotide from the expected position, thus restoring missing gamma-gamma' and IBS3-EBS3 pairings, but leaving the two conserved exonic ORFs out of frame. Because of the unexpected splice site, the principles for 3' intron definition were examined, which showed that the 3' splice site is flexible but contingent on gamma-gamma' and IBS3-EBS3 pairings, and can be as far away as four nucleotides from the wild-type site. Surprisingly, alternative splicing occurs at position +4 for wild-type B.a.I2 intron, both in vitro and in vivo, and the alternative event fuses the two conserved exon ORFs, presumably leading to translation of the downstream ORF. The finding suggests that the structural irregularities of B.a.I2 may be an adaptation to facilitate gene expression in vivo.

ABSTRACT Molecular epidemiological tools for genotyping clinical isolates of Mycobacterium tuberculosis have been developed and used to help track and contain transmission of tuberculosis. We identified 87 short sequence repeat loci within the genome of the M. tuberculosis H37Rv strain. Nine tandem repeats were found to be variable (variable-number tandem repeats [VNTRs]) in a set of 91 isolates. Fifty-seven of the isolates had only four IS 6110 bands. The other 34 isolates were members of the Beijing strain family. The number of alleles of each these nine VNTRs was determined by examining each isolate. Six of the loci (Mtb-v1, -v4, -v10, -v15, -v18, and -v20) were able to differentiate the Beijing spoligotype identical isolates into seven distinct genotypes. Five of the loci (Mtb-v3, -v5, -v6, -v10, and -v15) were informative in discriminating the four-band IS 6110 restriction fragment length polymorphism isolates from each other. The Nei's diversity values of each marker ranged from 0.02 to 0.59, with the number of alleles ranging from two to eight across the entire strain set. These nine loci provide a useful, discriminatory extension of VNTR typing methods for application to molecular epidemiologic studies of M. tuberculosis .

ABSTRACT Bacillus anthracis is a bacterial pathogen of great importance, both historically and in the present. This study presents data collected from several investigations and indicates that B. anthracis virulence is associated with the clonality and virulence of plasmids pXO1 and pXO2. Guinea pigs vaccinated with Anthrax Vaccine Adsorbed were challenged with 20 B. anthracis isolates representative of worldwide genetic diversity. These same isolates were characterized with respect to plasmid copy number by using a novel method of quantitative PCR developed for rapid and efficient detection of B. anthracis from environmental samples. We found that the copy numbers for both pXO1 and pXO2 differed from those in previously published reports. By combining the data on survival, plasmid copy numbers, and clonality, we developed a model predicting virulence. This model was validated by using a randomly chosen set of 12 additional B. anthracis isolates. Results from this study will be helpful in future efforts to elucidate the basis for variation in the virulence of this important pathogen.

ABSTRACT While outbreaks of animal anthrax zoonoses still regularly occur in France, little is known about the epidemiology links between them. We have used the eight-locus multilocus variable-number tandem repeat analysis typing technique against a collection of 50 Bacillus anthracis isolates from France. There were eight distinct genotypes belonging to two dissimilar genetic clusters. Regional strain patterns were observed, with the B2 genotypes prevalent in southern France and the A1a genotypes found only in northern France.

ABSTRACT Human Lyme borreliosis (LB) is the most prevalent arthropod-borne infection in temperate climate zones around the world and is caused by Borrelia spirochetes. We have identified 10 variable-number tandem repeat (VNTR) loci present within the genome of Borrelia burgdorferi and subsequently developed a multiple-locus VNTR analysis (MLVA) typing system for this disease agent. We report here the successful application of MLVA for strain discrimination among a group of 41 globally diverse Borrelia isolates including B. burgdorferi, B. afzelii, and B. garinii. PCR assays displayed diversity at these loci, with total allele numbers ranging from two to nine and Nei's diversity (D) values ranging from 0.10 to 0.87. The average D value was 0.53 across all VNTR loci. A clear correlation exists between the repeat copy number and the D value (r = 0.62) or the number of alleles (r = 0.93) observed across diverse strains. Cluster analysis by the unweighted pair-group method with arithmetic means resolved the 30 observed unique Borrelia genotypes into five distinct groups. B. burgdorferi, B. afzelii, and B. garinii clustered into distinct affiliations, consistent with current 16S rRNA phylogeny studies. Genetic similarity and diversity suggest that B. afzelii and B. garinii are close relatives and were perhaps recently derived from B. burgdorferi. MLVA provides both phylogenetic relationships and additional resolution to discriminate among strains of Borrelia species. This new level of strain identification and discrimination will allow more detailed epidemiological and phylogenetic analysis in future studies.

Comparison of the whole-genome sequence of Bacillus anthracis isolated from a victim of a recent bioterrorist anthrax attack with a reference reveals 60 new markers that include single nucleotide polymorphisms (SNPs), inserted or deleted sequences, and tandem repeats. Genome comparison detected four high-quality SNPs between the two sequenced B. anthracis chromosomes and seven differences among different preparations of the reference genome. These markers have been tested on a collection of anthrax isolates and were found to divide these samples into distinct families. These results demonstrate that genome-based analysis of microbial pathogens will provide a powerful new tool for investigation of infectious disease outbreaks.

ABSTRACT A test of 25 genetically diverse isolates of Bacillus anthracis was conducted to determine their susceptibility to seven clinically relevant antimicrobial agents. Etest strips (AB BIODISK, Solna, Sweden) were used to measure the MICs for the isolates. Using the National Committee for Clinical Laboratory Standards MIC breakpoints for staphylococci, three isolates were found to be resistant to penicillin and five were found to be resistant to cefuroxime. The penicillin-resistant isolates were negative for β-lactamase production. Continued surveillance of B. anthracis field isolates is recommended to monitor antimicrobial susceptibility.

Bacillus anthracis is a pathogen that is widely distributed around the globe. However, this great distribution is not accompanied by great genetic diversity. Although subtle morphological and biochemical differences exist, the underlying genetic basis for this plasticity is not known. Indeed, very few single nucleotide differences have been detected among isolates and the only documented high variable sequences are associated with variable number tandem repeated (VNTR) sequences. The differences among the VNTRs has been used to suggest phylogenetic relationships among the worldwide isolates. There is one major clonal split in B. anthracis (A and B), with two minor clusters established within the B branch and four or more minor branches within the A group. The A branch is the most common worldwide, though the B branch is locally important in certain areas. The ecology and the evolution of B. anthracis have been greatly influenced by the spore phase of the lifecycle.

ABSTRACT Rifampin-resistant mutants were selected from UV-light-treated Bacillus cereus (20 mutants) and attenuated B. anthracis (23 mutants). In addition, spontaneous rifampin-resistant mutants were also isolated in attenuated B. anthracis (22 mutants). The rifampin resistance clusters of the rpoB gene were sequenced for all 65 mutants. Mutations associated with resistance were consistent with those from other bacteria, though two novel changes were observed. The spontaneous rate of resistance was estimated at 1.57 × 10 −9 mutations/generation by a Luria-Delbrück fluctuation test.

ABSTRACT In 1993, the Aum Shinrikyo cult aerosolized Bacillus anthracis spores over Kameido, Japan. Spore samples were obtained from the release site, cultured, and characterized by molecular genetic typing. The isolates were consistent with strain Sterne 34F2, which is used in Japan for animal prophylaxis against anthrax.

ABSTRACT Francisella tularensis , the etiological agent of tularemia, is found throughout the Northern hemisphere. After analyzing the F. tularensis genomic sequence for potential variable-number tandem repeats (VNTRs), we developed a multilocus VNTR analysis (MLVA) typing system for this pathogen. Variation was detected at six VNTR loci in a set of 56 isolates from California, Oklahoma, Arizona, and Oregon and the F. tularensis live vaccine strain. PCR assays revealed diversity at these loci with total allele numbers ranging from 2 to 20, and Nei's diversity index values ranging from 0.36 to 0.93. Cluster analysis identified two genetically distinct groups consistent with the current biovar classification system of F. tularensis . These findings suggest that these VNTR markers are useful for identifying F. tularensis isolates at this taxonomic level. In this study, biovar B isolates were less diverse than those in biovar A, possibly reflecting the history of tularemia in North America. Seven isolates from a recent epizootic in Maricopa County, Ariz., were identical at all VNTR marker loci. Their identity, even at a hypervariable VNTR locus, indicates a common source of infection. This demonstrates the applicability of MLVA for rapid characterization and identification of outbreak isolates. Future construction of reference databases will allow faster outbreak tracking as well as providing a foundation for deciphering global genetic relationships.

ABSTRACT Bacillus anthracis is one of the most genetically homogeneous pathogens described, making strain discrimination particularly difficult. In this paper, we present a novel molecular typing system based on rapidly evolving variable-number tandem repeat (VNTR) loci. Multiple-locus VNTR analysis (MLVA) uses the combined power of multiple alleles at several marker loci. In our system, fluorescently labeled PCR primers are used to produce PCR amplification products from eight VNTR regions in the B. anthracis genome. These are detected and their sizes are determined using an ABI377 automated DNA sequencer. Five of these eight loci were discovered by sequence characterization of molecular markers ( vrrC 1 , vrrC 2 , vrrB 1 , vrrB 2 , and CG3), two were discovered by searching complete plasmid nucleotide sequences (pXO1-aat and pXO2-at), and one was known previously ( vrrA ). MLVA characterization of 426 B. anthracis isolates identified 89 distinct genotypes. VNTR markers frequently identified multiple alleles (from two to nine), with Nei's diversity values between 0.3 and 0.8. Unweighted pair-group method arithmetic average cluster analysis identified six genetically distinct groups that appear to be derived from clones. Some of these clones show worldwide distribution, while others are restricted to particular geographic regions. Human commerce doubtlessly has contributed to the dispersal of particular clones in ancient and modern times.

ABSTRACT The Kruger National Park (KNP), South Africa, has a recorded history of periodic anthrax epidemics causing widespread disease among wild animals. Bacillus anthracis is the causative agent of anthrax, a disease primarily affecting ungulate herbivores. Worldwide there is little diversity among B. anthracis isolates, but examination of variable-number tandem repeat (VNTR) loci has identified six major clones, with the most dissimilar types split into the A and B branches. Both the A and B types are found in southern Africa, giving this region the greatest genetic diversity of B. anthracis worldwide. Consequently, southern Africa has been hypothesized to be the geographic origin of B. anthracis . In this study, we identify the genotypic types of 98 KNP B. anthracis isolates using multiple-locus VNTR analysis. Two major types are evident, the A branch and the B branch. The spatial and temporal distribution of the different genotypes indicates that anthrax epidemic foci are independent, though correlated through environmental cues. Kruger B isolates were found on significantly higher-calcium and higher-pH soils than were Kruger type A. This relationship between genotype and soil chemistry may be due to adaptive differences among divergent anthrax strains. While this association may be simply fortuitous, adaptation of A types to diverse environmental conditions is consistent with their greater geographic dispersal and genetic dissimilarity.

Molecular typing of Bacillus anthracis has been extremely difficult due to the lack of polymorphic DNA markers. We have identified nine novel variable number tandemly repeated loci from previously known amplified fragment length polymorphism markers or from the DNA sequence. In combination with the previously known vrrA locus, these markers provide discrimination power to genetically characterize B. anthracis isolates. The variable number tandem repeat (VNTR) loci are found in both gene coding (genic) and non-coding (non-genic) regions. The genic differences are 'in frame' and result in additions or deletion of amino acids to the predicted proteins. Due the rarity of molecular differences, the VNTR changes represent a significant portion of the genetic variation found within B. anthracis. This variation could represent an important adaptive mechanism. Marker similarity and differences among diverse isolates have identified seven major diversity groups that may represent the only world-wide B. anthracis clones. The lineages reconstructed using these data may reflect the dispersal and evolution of this pathogen.

It has only recently been possible to detect sufficient genetic diversity among anthrax isolates to allow genotype grouping (Keim et al. 1997). Early results of such grouping suggest that the southern African subcontinent may be the geographical origin of Bacillus anthracis. This report describes a pilot investigation of the genetic diversity of a study group of isolates from the Kruger National Park, South Africa, and efforts to detect spatio-temporal clustering within the study group. This study has also served as further validation for the newly developed Multi-Locus VNTR Analysis (MLVA), designed to simplify genotyping of B. anthracis isolates. The results reveal a diverse range of genotypes within the park allied with three genotype reference groups, and show that the MLVA procedure is robust for rapid analysis of B. anthracis genotypes. We also observed multiple genotype groups within epidemics and between geographically and temporally close epidemic episodes. This is in contrast to earlier characterizations of anthrax epidemics. The result of a Mantel test for time-space clustering indicates clustering of the anthrax isolates selected for the study.