Molecular diversity in Bacillus anthracis

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.     

vrrB, a hypervariable open reading frame in Bacillus anthracis

Bacillus anthracis appears to be the most molecularly homogeneous bacterial species known. Extensive surveys of worldwide isolates have revealed vanishingly small amounts of genomic variation. The biological importance of the resting-stage spore may lead to very low evolutionary rates and, perhaps, to the lack of potentially adaptive genetic variation. In contrast to the overall homogeneity, some gene coding regions contain hypervariability that is translated into protein variation. During marker analysis of diverse strains, we have discovered a novel ca. 750-nucleotide open reading frame (ORF) that contains in-frame, variable-number tandem-repeat sequences. Four distinct variable regions exist within vrrB, giving rise to 11 distinct alleles in eight different length categories among B. anthracis strains. This ORF putatively codes for a 241- to 265-amino-acid protein, rich in glutamine (13.2%), glycine (23.4%), and histidine (23.0%). The variable-region amino acids of the vrrB ORF are strongly hydrophilic. Coupled with putative transmembrane domains flanking the variable regions, this suggests a membrane-anchored cytosolic or extracellular location for the putative protein. Sequence analysis of the complete ORFs from three Bacillus cereus strains shows maintenance of the ORF across species boundaries, including strong conservation of the amino acid sequence and the capacity to vary among strains. The presence of 11 different alleles of the vrrB locus is in stark contrast to the near homogeneity of B. anthracis. Evolution of hypervariable genes can negate the lack of genetic variability in species such as B. anthracis and provide select rapid evolution in other more variable species.     

Comparison of minisatellite polymorphisms in the Bacillus cereus complex: a simple assay for large-scale screening and identification of strains most closely related to Bacillus anthracis

Polymorphism of five tandem repeats that are monomorphic in Bacillus anthracis was investigated in 230 isolates of the B. cereus group and in 5 sequenced B. cereus genomes in search for markers allowing identification of B. cereus and B. thuringiensis strains most closely related to B. anthracis. Using this multiple-locus variable number of tandem repeat analysis (MLVA), a cluster of 30 strains was selected for further characterization. Eventually, six of these were characterized by multilocus sequence type analysis. One of the strains is only six point mutations (of almost 3,000 bp) away from B. anthracis and was also proposed to be closest to B. anthracis by MLVA analysis. However, this strain remains separated from B. anthracis by a number of significant genetic events observed in B. anthracis, including the loss of the hemolysin activity, the presence of four prophages, and the presence of the two virulence plasmids, pXO1 and pXO2. One particular minisatellite marker provides an efficient assay to identify the subset of B. cereus and B. thuringiensis strains closely related to B. anthracis. Based on these results, a very simple assay is proposed that allows the screening of hundreds of strains from the B. cereus complex, with modest equipment and at a low cost, to eventually fill the gap with B. anthracis and better understand the origin and making of this dangerous pathogen.     

Molecular evolution and diversity in Bacillus anthracis as detected by amplified fragment length polymorphism markers.

Bacillus anthracis causes anthrax and represents one of the most molecularly monomorphic bacteria known. We have used AFLP (amplified fragment length polymorphism) DNA markers to analyze 78 B. anthracis isolates and six related Bacillus species for molecular variation. AFLP markers are extremely sensitive to even small sequence variation, using PCR and high-resolution electrophoresis to examine restriction fragments. Using this approach, we examined ca. 6.3% of the Bacillus genome for length mutations and ca. 0.36% for point mutations. Extensive variation was observed among taxa, and both cladistic and phenetic analyses were used to construct a phylogeny of B. anthracis and its closest relatives. This genome-wide analysis of 357 AFLP characters (polymorphic fragments) indicates that B. cereus and B. thuringiensis are the closest taxa to B. anthracis, with B. mycoides slightly more distant. B. subtilis, B. polymyxa, and B. stearothermophilus shared few AFLP markers with B. anthracis and were used as outgroups to root the analysis. In contrast to the variation among taxa, only rare AFLP marker variation was observed within B. anthracis, which may be the most genetically uniform bacterial species known. However, AFLP markers did establish the presence or absence of the pXO1 and pXO2 plasmids and detected 31 polymorphic chromosomal regions among the 79 B. anthracis isolates. Cluster analysis identified two very distinct genetic lineages among the B. anthracis isolates. The level of variation and its geographic distribution are consistent with a historically recent African origin for this pathogenic organism. Based on AFLP marker similarity, the ongoing anthrax epidemic in Canada and the northern United States is due to a single strain introduction that has remained stable over at least 30 years and a 1,000-mile distribution.     

Molecular characterization of Korean Bacillus anthracis isolates by amplified fragment length polymorphism analysis and multilocus variable-number tandem repeat analysis

We analyzed the genetic relationships and molecular characteristics of 34 Bacillus anthracis isolates from soil and clinical samples in various regions of Korea and 17 related Bacillus species, using the amplified fragment length polymorphism (AFLP) and multilocus variable-number tandem repeat (MLVA) approaches. Triplicate AFLP profiles of these strains showed high reproducibility and identified 376 polymorphisms. AFLP phylogenetic analysis of B. anthracis isolates showed a high level of similarity, 0.93, and this monomorphic fragment profile proved to be useful to differentiate B. anthracis strains from other Bacillus species. The B. cereus group was separated from other Bacillus species at a level of similarity of 0.68. Among them, some B. cereus strains showed genetic interspersion with B. thuringiensis strains. The evolutionary pattern of nucleotide differences among B. anthracis strains with the eight MLVA markers showed nine MLVA types. Three MLVA types, M1 to M3, were pathogenic B. anthracis isolates and were assigned as new genotypes belonging to the A4 and B3 clusters, compared with 89 genotypes deduced from previous data. This indicates that differences in cluster prevalence and distribution may be influenced more by MLVA markers on two plasmids loci and human activity. Consequently, we suggest that the novel MLVA type may represent significant evidence for historic adaptation to environmental conditions of the Asian continent, particularly Korea. Therefore, MLVA techniques may be available for molecular monitoring on anthrax-release-related bioterrorism and further study is required for the continuous epidemiological study of variable anthrax collections.     

Identification and characterization of Bacillus anthracis by multiplex PCR on DNA chip

Bacillus anthracis can be identified by detecting virulence factor genes located on two plasmids, pXO1 and pXO2. Combining multiplex PCR with arrayed anchored primer PCR and biotin-avidin alkaline phosphatase indicator system, we developed a qualitative DNA chip method for characterization of B. anthracis, and simultaneous confirmation of the species identity independent of plasmid contents. The assay amplifies pag gene (in pXO1), cap gene (in pXO2) and Ba813 gene (a B. anthracis specific chromosomal marker), and the results were indicated by an easy-to-read profile based on the color reaction of alkaline phosphatase. About 1 pg of specific DNA fragments on the chip wells could be detected after PCR. With the proposed method, the avirulent (pXO1+/2-, pXO1-/2+ and pXO1-/2-) strains of B. anthracis and distinguished 'anthrax-like' strains from other B. cereus group bacteria were unambiguously identified, while the genera other than Bacillus gave no positive signal.     

Use of multiple‐locus variable number tandem repeat analysis and phage typing for subtyping ofSalmonella Enteritidis from sporadic human cases in the United States

Aims: To investigate the genetic diversity among S. Enteritidis isolates from different geographic regions to evaluate the relationship between phage types (PTs) and variable number tandem repeat analysis (VNTR) loci. Methods and Results: We performed multiple‐locus variable number tandem repeat analysis (MLVA) and phage typing on 245 S. Enteritidis isolates collected from sporadic human clinical cases in Michigan, Minnesota, New York, and Washington states between 2000 and 2007. Ninety‐four MLVA types and 22 different PTs were identified. Specific PTs were associated with a predominant allele for certain VNTR loci. Cluster analysis using a minimum‐spanning tree demonstrated two major clusters (I, II) and one minor cluster of isolates. PTs 8, 13a, 13 and 34 were significantly associated with MLVA cluster I. Phage types 1, 4, 6a, and 18 were significantly associated with MLVA cluster II. Conclusions: We found significant association between MLVA‐based clusters and PTs. Certain VNTR loci were associated with specific PTs and could serve as useful molecular markers for S. Enteritidis in epidemiological investigations. Significance and Impact of the Study: MLVA genotyping in combination with phage typing can be used for effective characterization of S. Enteritidis isolates. It can also be useful for tracing possible sources during investigations of sporadic and outbreak cases of S. Enteritidis.     

Virulotypes of Bacillus anthracis strains isolated in Poland

Bacillus anthracis--the causative agent of anthrax--possesses several virulence genes located in the chromosome as well as in two B. anthracis virulence plasmids: pXO1 and pXO2. In the presented study, we determined occurrence of six virulence markers located in the virulence plasmids (capA, capB, capC, pagA, lef and cya) for capsule and toxin production together with virulence-associated gene gerXA and chromosomal gene sap, which are responsible for germination and S-layer biosynthesis respectively. Fourteen strains of B. anthracis isolated in Poland and belonging to five different MLVA genotypes were analyzed by PCR for presence of the aforementioned genes. Two virulotypes were found in tested strains. The only variation was absence of capA, capB and capC due to a lack of pXO2.     

Establishment of a genotyping scheme for Coxiella burnetii

Coxiella burnetii is the causative agent of Q fever. The bacterium is highly infectious and is classified as a category B biological weapon. The tools of molecular biology are of utmost importance in a rapid and unambiguous identification of C. burnetii in naturally occurring Q fever outbreaks, or in cases of a deliberate release of the infectious agent. In this work, development of a multiple locus variable number tandem repeats (VNTR) analysis (MLVA) for the characterization of C. burnetii is described. Sixteen C. burnetii isolates and five passage history/laboratory variants were characterized. The VNTR markers revealed many polymorphisms resulting in nine unique MLVA types that cluster into five different clusters. This proves that the MLVA system is highly discriminatory. The selected VNTR markers were stable. The MLVA method developed in this report is a promising tool for the characterization of C. burnetii isolates and their epidemiological study.     

The use of MLVA for Corynebacterium diphtheriae genotyping--preliminary studies

The complete genome sequence of strain NCTC 13129 C. diphtheriae were investigated in order to identify tandem repeats (VNTR). From 75 VNTR loci identified in the genome 14 were selected. Primers were designed and PCR conditions were optimized for amplification of the selected VNTR markers. Preliminary studies of usefulness of selected VNTR markers were conducted using a group of 28 C. diphtheriae strains. From 14 markers 8 were regarded as potentially useful. The diversity of individual markers ranged from 1 to 6 alleles (Simpson index from 0 to 0,746). No diversity were observed for 3 VNTR markers but it could be a results of too small group of strains analyzed in the tests. Simpson diversity index calculated for all the markers tested on 28 strains was 0,87. Results of the preliminary studies showed usefulness of MLVA for C. diphtheriae genotyping. Nevertheless, confirmation of reliability of the method should be done using a large group of strains. Moreover, the method should be compared with other genotyping methods.     

Highly diverse variable number tandem repeat loci in the E. coli O157:H7 and O55:H7 genomes for high-resolution molecular typing

AIM: Evaluation of the Escherichia coli genome for variable number tandem repeat (VNTR) loci in order to provide a subtyping tool with greater discrimination and more efficient capacity. METHODS AND RESULTS: Twenty-nine putative VNTR loci were identified from the E. coli genomic sequence. Their variability was validated by characterizing the number of repeats at each locus in a set of 56 E. coli O157:H7/HN and O55:H7 isolates. An optimized multiplex assay system was developed to facility high capacity analysis. Locus diversity values ranged from 0.23 to 0.95 while the number of alleles ranged from two to 29. This multiple-locus VNTR analysis (MLVA) data was used to describe genetic relationships among these isolates and was compared with PFGE (pulse field gel electrophoresis) data from a subset of the same strains. Genetic similarity values were highly correlated between the two approaches, through MLVA was capable of discrimination amongst closely related isolates when PFGE similar values were equal to 1.0. CONCLUSIONS: Highly variable VNTR loci exist in the E. coli O157:H7 genome and are excellent estimators of genetic relationships, in particular for closely related isolates. SIGNIFICANCE AND IMPACT OF THE STUDY: Escherichia coli O157:H7 MLVA offers a complimentary analysis to the more traditional PFGE approach. Application of MLVA to an outbreak cluster could generate superior molecular epidemiology and result in a more effective public health response.     

Mass spectrometry provides accurate characterization of two genetic marker types in Bacillus anthracis

Epidemiological and forensic analyses of bioterrorism events involving Bacillus anthracis could be improved if both variable number tandem repeats (VNTRs) and single nucleotide polymorphisms (SNPs) could be combined on a single analysis platform. Here we present the use of electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICR-MS) to characterize 24 alleles from 6 VNTR loci and 11 alleles from 7 SNP loci in B. anthracis. The results obtained with ESI-FTICR-MS were consistent with independent results obtained from traditional approaches using electrophoretic detection of fluorescent products. However, ESI-FTICR-MS improves on the traditional approaches because it does not require fluorescent labeling of PCR products, minimizes post-PCR processing, obviates electrophoresis, and provides unambiguous base composition of both SNP and VNTR PCR products. In addition, ESI-FTICR-MS allows both marker types to be examined simultaneously and at a rate of approximately 1 sample per min. This technology represents a significant advance in our ability to rapidly characterize B. anthracis isolates using VNTR and SNP loci.     

Identification and characterization of Bacillus anthracis by multiplex PCR analysis of sequences on plasmids pXO1 and pXO2 and chromosomal DNA

Abstract Bacillus anthracis can be identified on the basis of the detection of virulence factor genes located on two plasmids, pXO1 and pXO2. Thus isolates lacking both pXO1 and pXO2 are indistinguishable from closely related B. cereus group bacteria. We developed a multiplex PCR assay for characterization of B. anthracis isolates, and simultaneous confirmation of the species identity independent of plasmid content. The assay amplifies lef, cya, pag (pXO1) and cap (pXO2) genes, and a B. anthracis specific chromosomal marker, giving an easy-to-read profile. This system unambiguously identified virulent (pXO1⁺/2⁺) and avirulent (pXO1⁺/2⁻, pXO1⁻/2⁺ and pXO1⁻/2⁻) strains of B. anthracis and distinguished 'anthrax-like' strains from other B. cereus group bacteria.     

Multiple-locus variable-number tandem repeat analysis of Salmonella Enteritidis isolates from human and non-human sources using a single multiplex PCR

Simplified multiple-locus variable-number tandem repeat analysis (MLVA) was developed using one-shot multiplex PCR for seven variable-number tandem repeats (VNTR) markers with high diversity capacity. MLVA, phage typing, and PFGE methods were applied on 34 diverse Salmonella Enteritidis isolates from human and non-human sources. MLVA detected allelic variations that helped to classify the S. Enteritidis isolates into more evenly distributed subtypes than other methods. MLVA-based S. Enteritidis clonal groups were largely associated with sources of the isolates. Nei's diversity indices for polymorphism ranged from 0.25 to 0.70 for seven VNTR loci markers. Based on Simpson's and Shannon's diversity indices, MLVA had a higher discriminatory power than pulsed field gel electrophoresis (PFGE), phage typing, or multilocus enzyme electrophoresis. Therefore, MLVA may be used along with PFGE to enhance the effectiveness of the molecular epidemiologic investigation of S. Enteritidis infections.     

Multiple-locus variable-number of tandem-repeats analysis distinguishes Vibrio parahaemolyticus pandemic O3:K6 strains

A specific serotype of Vibrio parahaemolyticus, O3:K6, has recently been linked to epidemics of gastroenteritis in Southeast Asia, Japan, and North America. These pandemic O3:K6 strains appear to have recently spread across continents from a single origin to reach global coverage, based on profiling of strains by several molecular typing methods. In this study, variable-number tandem repeats (VNTR)-based fingerprinting was applied to clinical and environmental V. parahaemolyticus O3:K6 strains in an attempt to develop a molecular method with increased sensitivity for discriminating strains; the relative discriminatory powers were compared with ribotyping and pulsed-field gel electrophoresis (PFGE). All clinical strains tested were independent human isolates obtained from different outbreaks or from sporadic cases in Tokyo during the period from 1996 to 2003. Multiple-locus VNTR analysis (MLVA) was shown to have high resolution and reproducibility for typing of V. parahaemolyticus clones. MLVA analysis of 28 pandemic V. parahaemolyticus O3:K6 strains isolated from human cases produced 28 distinct VNTR patterns. The VNTR loci displayed between 2 and 15 alleles at each of eight loci with Nei's diversity index ranging from 0.35 and 0.91. These data demonstrated that MLVA is useful for individual strain typing of new O3:K6 strains, which appear to be closely related by other molecular methods.     

Molecular characterization of Bacillus anthracis using multiplex PCR, ERIC-PCR and RAPD

AIMS: To investigate the molecular characterization of Bacillus anthracis strains by multiplex PCR, enterobacterial repetitive intergenic consensus-PCR (ERIC-PCR) and random amplification of polymorphic DNA (RAPD). METHODS AND RESULTS: Three primers were used to amplify the cya, cap and cereolysinAB genes in the multiplex PCR. Two distinct ERIC-PCR and RAPD fragments, which separated B. anthracis into two groups, were used as probes in Southern hybridization experiments. The probes hybridized only to the cya+ B. anthracis strains identified by the multiplex PCR. Nucleotide sequence analysis of the two cloned fragments showed they were from the pXO1 plasmid of B. anthracis. CONCLUSION: Multiplex PCR simultaneously identified isolates of the Bacillus cereus group and the B. anthracis virulence factors. ERIC-PCR and RAPD, combined with the Southern hybridization analyses, differentiated B. anthracis strains and separated them from the closely related B. cereus group bacteria. SIGNIFICANCE AND IMPACT OF THE STUDY: ERIC-PCR and RAPD assay could be effective in differentiating virulent from avirulent B. anthracis. Our results also show that the amplification of the large plasmids was allowed in the ERIC-PCR and RAPD assay.     

Markers for species-specific detection of bacteria from Bacillus cereus group

rpoB and gyr genes (and their fragments) of chromosomal DNA of bacteria from Bacillus cereus group - B. anthracis, B. cereus, and B. thuringiensis - which are the potential markers for their genotyping were sequenced and phylogenetic trees were constructed. Sets of primers for species-specific detection of B. anthracis, B. cereus, and B. thuringiensis by multiplex polymerase chain reaction were designed. Also primers sets, which allow to differentiate strains of B. anthracis with various plasmid profiles (containing both plasmids (pXO1+, pXO2+), and without one (pXO1+, pXO2- or pXO1-, pXO2+) or both plasmids (pXO1-, pXO2-), determining pathogenic characteristics of the strains, were developed. For multiplex PCR primer sets were optimized on the annealing temperature of primers and amplicon length. Itwas shown that phylogenetic tree can be applied as an indicator of reliability and accuracy of taxonomical classification of microorganisms' species and subspecies. Comparison of pXO1 and pXO2 plasmid sequences of B. anthracis showed that these plasmids contain 18 and 4 palindrome sequences respectively which can potentially form thermodynamically stable hairpin-loop structures.     

Occurrence and genetic diversity of Bacillus anthracis strains isolated in an active wool-cleaning factory

Culturable microorganisms from various samples taken at an active factory performing wool and goat hair cleaning were isolated and analyzed. Bacillus anthracis was found in air filter dust, wastewater, and goat hairs, where it accounted for approximately 1% of the total counts of viable bacteria. Consistent with the countries of origin of the processed material (South Caucasian and Middle Eastern), all B. anthracis isolates belonged to the same phylogenetic cluster, as determined by variable-number tandem repeat (VNTR) typing at eight loci. Within this cluster, five closely related VNTR subtypes could be identified, of which two were previously unreported. Additional diversity was observed when more sensitive genetic markers were assayed, demonstrating the multifocal nature of goat hair contamination. Goat hair originating from areas where anthrax is endemic remains a material with high biological risk for modern woolworkers.     

The use of synthetic tandem repeats to isolate new VNTR loci: cloning of a human hypermutable sequence.

Synthetic tandem repeats (STRs) of oligonucleotides have previously been shown to detect polymorphic loci in the human genome. Here, we report results from the use of three such probes to screen a human cosmid library. Nine of the 45 positive clones that were analyzed appear to contain highly polymorphic minisatellite or VNTR loci. The degree of enrichment for minisatellite sequences varied with the choice of STR: one provided a 15- to 20-fold enrichment (4 polymorphic loci among 10 clones), whereas 2 others gave a 3- to 5-fold enrichment (5 polymorphic probes in a total of 35 clones) compared to random screening. The 9 VNTR markers have been localized by linkage analysis in the CEPH panel and/or by in situ hybridization. Eight probes identify new loci, one of which maps to an interstitial region. One of the VNTR loci (identified by probe CEB1) was found to be hypermutable, with 52 mutation events identified among 310 children characterized in 40 CEPH families. The parental origin of the mutation could be identified in all instances, and only one mutation was found to be of maternal origin. The mutation rate in males was estimated to be approximately 15%. Segregation analysis of flanking markers suggests that mutations are not associated with crossing over. As the only previously described hypermutable minisatellite loci in humans have equal rates of male and female mutations, these observations establish that a second type of hypermutable minisatellite exists in the human genome. In neither case does the generation of new alleles appear to be associated with unequal crossing over.     

Molecular characterisation of the haemolytic isolates of Bacillus anthracis originated in Poland

Bacillus anthracis is generally considered non-haemolytic, when cultured on the solid media. However, strains capable to lyse sheep erythrocytes have been reported. Anthrolysin O, an orthologue of cereolysin was proposed as a putative haemolysin of B. anthracis. AIM: to determine whether anthrolysin O, haemolytic enterotoxin HBL and the pleiotropic regulator PlcR that activates antrholysin O production are associated with a haemolytic activity of B. anthracis strains isolated in Poland. MATERIAL: in total 8 B. anthracis strains - the fully virulent BL1 and seven the pXO2 lacking strains including: a vaccine strain Sterne 34F2 together with three haemolytic and three non-haemolytic strains isolated from different samples of the same animal died from anthrax in Poland. METHODS: The haemolytic activity was detected using Columbia agar plates supplemented with 5% of sheep blood. Anthrolvsin O, cereolysin and gene hblA encoding the key subunit of the HBL were detected by PCR. In addition, the plcR gene fragment containing the B. anthracis specific non-sense mutation was analysed by the DNA sequencing. Ten marker loci based MLVA genotyping was performed to distinguish tested strains. RESULTS: The alo gene encoding anthrolysin O was detected in both the haemolytic and non-haemolytic strains while hblA was absent. The B. anthracis specific plcR non-sense mutation was detected in both the groups of tested strains, suggesting that the haemolysis in tested strains may rather be conferred by the PlcR-independent factors. Moreover, haemolytic and non-haemolytic strains were indistinguishable by the MLVA. Obtained results may argue the haemolytic and non-haemolytic strains are isogenic and most probably a single mutational event is responsible for the haemolytic phenotype induction.