Global genetic population structure of Bacillus anthracis.

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.

Tandem repeat regions within the Burkholderia pseudomallei genome and their application for high resolution genotyping.

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.

Mutations, mutation rates, and evolution at the hypervariable VNTR loci of Yersinia pestis.

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.

Principles of 3' splice site selection and alternative splicing for an unusual group II intron from Bacillus anthracis.

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.

Molecular typing of Mycobacterium tuberculosis by using nine novel variable-number tandem repeats across the Beijing family and low-copy-number IS6110 isolates.

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 IS6110 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 IS6110 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.

Bacillus anthracis virulence in Guinea pigs vaccinated with anthrax vaccine adsorbed is linked to plasmid quantities and clonality.

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.

Strain typing of Borrelia burgdorferi, Borrelia afzelii, and Borrelia garinii by using multiple-locus variable-number tandem repeat analysis.

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.

Diversity among French Bacillus anthracis isolates.

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.

Antimicrobial susceptibilities of diverse Bacillus anthracis isolates.

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 beta-lactamase production. Continued surveillance of B. anthracis field isolates is recommended to monitor antimicrobial susceptibility.

Comparative genome sequencing for discovery of novel polymorphisms in Bacillus anthracis.

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.

Molecular analysis of rifampin resistance in Bacillus anthracis and Bacillus cereus.

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 x 10(-9) mutations/generation by a Luria-Delbrück fluctuation test.