Coalescence modeling of intrainfection Bacillus anthracis populations allows estimation of infection parameters in wild populations.

Bacillus anthracis, the etiological agent of anthrax, is a well-established model organism. For B. anthracis and most other infectious diseases, knowledge regarding transmission and infection parameters in natural systems, in large part, comprises data gathered from closely controlled laboratory experiments. Fatal, natural anthrax infections transmit the bacterium through new host-pathogen contacts at carcass sites, which can occur years after death of the previous host. For the period between contact and death, all of our knowledge is based upon experimental data from domestic livestock and laboratory animals. Here we use a noninvasive method to explore the dynamics of anthrax infections, by evaluating the terminal diversity of B. anthracis in anthrax carcasses. We present an application of population genetics theory, specifically, coalescence modeling, to intrainfection populations of B. anthracis to derive estimates for the duration of the acute phase of the infection and effective population size converted to the number of colony-forming units establishing infection in wild plains zebra (Equus quagga). Founding populations are small, a few colony-forming units, and infections are rapid, lasting roughly between 1 d and 3 d in the wild. Our results closely reflect experimental data, showing that small founding populations progress acutely, killing the host within days. We believe this method is amendable to other bacterial diseases from wild, domestic, and human systems.

Lethal exposure: An integrated approach to pathogen transmission via environmental reservoirs.

To mitigate the effects of zoonotic diseases on human and animal populations, it is critical to understand what factors alter transmission dynamics. Here we assess the risk of exposure to lethal concentrations of the anthrax bacterium, Bacillus anthracis, for grazing animals in a natural system over time through different transmission mechanisms. We follow pathogen concentrations at anthrax carcass sites and waterholes for five years and estimate infection risk as a function of grass, soil or water intake, age of carcass sites, and the exposure required for a lethal infection. Grazing, not drinking, seems the dominant transmission route, and transmission is more probable from grazing at carcass sites 1-2 years of age. Unlike most studies of virulent pathogens that are conducted under controlled conditions for extrapolation to real situations, we evaluate exposure risk under field conditions to estimate the probability of a lethal dose, showing that not all reservoirs with detectable pathogens are significant transmission pathways.

Bacillus anthracis comparative genome analysis in support of the Amerithrax investigation.

Before the anthrax letter attacks of 2001, the developing field of microbial forensics relied on microbial genotyping schemes based on a small portion of a genome sequence. Amerithrax, the investigation into the anthrax letter attacks, applied high-resolution whole-genome sequencing and comparative genomics to identify key genetic features of the letters' Bacillus anthracis Ames strain. During systematic microbiological analysis of the spore material from the letters, we identified a number of morphological variants based on phenotypic characteristics and the ability to sporulate. The genomes of these morphological variants were sequenced and compared with that of the B. anthracis Ames ancestor, the progenitor of all B. anthracis Ames strains. Through comparative genomics, we identified four distinct loci with verifiable genetic mutations. Three of the four mutations could be directly linked to sporulation pathways in B. anthracis and more specifically to the regulation of the phosphorylation state of Spo0F, a key regulatory protein in the initiation of the sporulation cascade, thus linking phenotype to genotype. None of these variant genotypes were identified in single-colony environmental B. anthracis Ames isolates associated with the investigation. These genotypes were identified only in B. anthracis morphotypes isolated from the letters, indicating that the variants were not prevalent in the environment, not even the environments associated with the investigation. This study demonstrates the forensic value of systematic microbiological analysis combined with whole-genome sequencing and comparative genomics.

Historical distribution and molecular diversity of Bacillus anthracis, Kazakhstan.

To map the distribution of anthrax outbreaks and strain subtypes in Kazakhstan during 1937-2005, we combined geographic information system technology and genetic analysis by using archived cultures and data. Biochemical and genetic tests confirmed the identity of 93 archived cultures in the Kazakhstan National Culture Collection as Bacillus anthracis. Multilocus variable number tandem repeat analysis genotyping identified 12 genotypes. Cluster analysis comparing these genotypes with previously published genotypes indicated that most (n = 78) isolates belonged to the previously described A1.a genetic cluster, 6 isolates belonged to the A3.b cluster, and 2 belonged to the A4 cluster. Two genotypes in the collection appeared to represent novel genetic sublineages; 1 of these isolates was from Krygystan. Our data provide a description of the historical, geographic, and genetic diversity of B. anthracis in this Central Asian region.

Comparison of two multiple-locus variable-number tandem-repeat analysis methods for molecular strain typing of human Brucella melitensis isolates from the Middle East.

Brucella species are highly monomorphic, with minimal genetic variation among species, hindering the development of reliable subtyping tools for epidemiologic and phylogenetic analyses. Our objective was to compare two distinct multiple-locus variable-number tandem-repeat analysis (MLVA) subtyping methods on a collection of 101 Brucella melitensis isolates from sporadic human cases of brucellosis in Egypt (n = 83), Qatar (n = 17), and Libya (n = 1). A gel-based MLVA technique, MLVA-15(IGM), was compared to an automated capillary electrophoresis-based method, MLVA-15(NAU), with each MLVA scheme examining a unique set of variable-number tandem repeats. Both the MLVA(IGM) and MLVA(NAU) methods were highly discriminatory, resolving 99 and 101 distinct genotypes, respectively, and were able to largely separate genotypes from Egypt and Qatar. The MLVA-15(NAU) scheme presented higher strain-to-strain diversity in our test population than that observed with the MLVA-15(IGM) assay. Both schemes were able to genetically correlate some strains originating from the same hospital or region within a country. In addition to comparing the genotyping abilities of these two schemes, we also compared the usability, limitations, and advantages of the two MLVA systems and their applications in the epidemiological genotyping of human B. melitensis strains.

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.

Strain-specific single-nucleotide polymorphism assays for the Bacillus anthracis Ames strain.

Highly precise diagnostics and forensic assays can be developed through a combination of evolutionary analysis and the exhaustive examination of genomic sequences. In Bacillus anthracis, whole-genome sequencing efforts revealed ca. 3,500 single-nucleotide polymorphisms (SNPs) among eight different strains and evolutionary analysis provides the identification of canonical SNPs. We have previously shown that SNPs are highly evolutionarily stable, and the clonal nature of B. anthracis makes them ideal signatures for subtyping this pathogen. Here we identified SNPs that define the lineage of B. anthracis that contains the Ames strain, the strain used in the 2001 bioterrorist attacks in the United States. Sequencing and real-time PCR were used to validate these SNPs across B. anthracis strains, including (i) 88 globally and genetically diverse isolates; (ii) isolates that were shown to be genetic relatives of the Ames strain by multiple-locus variable number tandem repeat analysis (MLVA); and (iii) several different lab stocks of the Ames strain, including a clinical isolate from the 2001 letter attack. Six SNPs were found to be highly specific for the Ames strain; four on the chromosome, one on the pX01 plasmid, and one on the pX02 plasmid. All six SNPs differentiated the B. anthracis Ames strain from the 88 unique B. anthracis strains, while five of the six separated Ames from its close genetic relatives. The use of these SNPs coupled with real-time PCR allows specific and sensitive (<100 fg of template DNA) identification of the Ames strain. This evolutionary and genomics-based approach provides an effective means for the discovery of strain-specific SNPs in B. anthracis.

Use of a real-time PCR TaqMan assay for rapid identification and differentiation of Burkholderia pseudomallei and Burkholderia mallei.

A TaqMan allelic-discrimination assay designed around a synonymous single-nucleotide polymorphism was used to genotype Burkholderia pseudomallei and Burkholderia mallei isolates. The assay rapidly identifies and discriminates between these two highly pathogenic bacteria and does not cross-react with genetic near neighbors, such as Burkholderia thailandensis and Burkholderia cepacia.

Specific detection of bacillus anthracis using a TaqMan mismatch amplification mutation assay.

Single nucleotide polymorphisms (SNPs) are increasingly recognized as important diagnostic markers for the detection and differentiation of Bacillus anthracis. The use of SNP markers for identifying B. anthracis DNA in environmental samples containing genetically similar bacteria requires the ability to amplify and detect DNA with single nucleotide specificity. We designed a TaqMan mismatch amplification mutation assay (TaqMAMA) around a SNP in the plcR gene of B. anthracis. The assay permits specific, low-level detection (25 fg DNA) of this B. anthracis-specific SNP, even in the presence of environmental DNA extracts containing a 20,000-fold excess of the alternate allele. We anticipate that the ability to selectively amplify and detect low copy number DNAs with single nucleotide specificity will represent a valuable tool in the arena of biodefense and microbial forensics.

Characterization of a Bacillus anthracis isolate causing a rare case of fatal anthrax in a 2-year-old boy from Hong Kong.

We used multiple-locus variable-number tandem repeat analysis (MLVA) and pagA sequencing to genotype a Bacillus anthracis isolate from a fatal case of human anthrax in Hong Kong. The isolate has a unique MLVA genotype, is related to the Sterne and Ames strains, and is consistent with genotypes identified in China.

Use of single nucleotide polymorphisms in the plcR gene for specific identification of Bacillus anthracis.

A TaqMan-minor groove binding assay designed around a nonsense mutation in the plcR gene was used to genotype Bacillus anthracis, B. cereus, and B. thuringiensis isolates. The assay differentiated B. anthracis from these genetic near-neighbors and determined that the nonsense mutation is ubiquitous across 89 globally and genetically diverse B. anthracis strains.

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.

Phylogenetic discovery bias in Bacillus anthracis using single-nucleotide polymorphisms from whole-genome sequencing.

Phylogenetic reconstruction using molecular data is often subject to homoplasy, leading to inaccurate conclusions about phylogenetic relationships among operational taxonomic units. Compared with other molecular markers, single-nucleotide polymorphisms (SNPs) exhibit extremely low mutation rates, making them rare in recently emerged pathogens, but they are less prone to homoplasy and thus extremely valuable for phylogenetic analyses. Despite their phylogenetic potential, ascertainment bias occurs when SNP characters are discovered through biased taxonomic sampling; by using whole-genome comparisons of five diverse strains of Bacillus anthracis to facilitate SNP discovery, we show that only polymorphisms lying along the evolutionary pathway between reference strains will be observed. We illustrate this in theoretical and simulated data sets in which complex phylogenetic topologies are reduced to linear evolutionary models. Using a set of 990 SNP markers, we also show how divergent branches in our topologies collapse to single points but provide accurate information on internodal distances and points of origin for ancestral clades. These data allowed us to determine the ancestral root of B. anthracis, showing that it lies closer to a newly described "C" branch than to either of two previously described "A" or "B" branches. In addition, subclade rooting of the C branch revealed unequal evolutionary rates that seem to be correlated with ecological parameters and strain attributes. Our use of nonhomoplastic whole-genome SNP characters allows branch points and clade membership to be estimated with great precision, providing greater insight into epidemiological, ecological, and forensic questions.

Anthrax molecular epidemiology and forensics: using the appropriate marker for different evolutionary scales.

Precise identification of Bacillus anthracis isolates has aided forensic and epidemiological analyses of natural anthrax cases, bioterrorism acts and industrial scale accidents by state-sponsored bioweapons programs. Because there is little molecular variation among B. anthracis isolates, identifying and using rare variation is crucial for precise strain identification. We think that mutation is the primary diversifying force in a clonal, recently emerged pathogen, such as B. anthracis, since mutation rate is correlated with diversity on a per locus basis. While single nucleotide polymorphisms (SNPs) are rare, their detection is facilitated by whole genome discovery approaches. As highly stable phylogenetic markers, SNPs are useful for identifying long branches or key phylogenetic positions. Selection of single, diagnostic "Canonical SNPs" (canSNPs) for these phylogenetic positions allows for efficient and defining assays. We have taken a nested hierarchal strategy for subtyping B. anthracis, which is consistent with traditional diagnostics and applicable to a wide range of pathogens. Progressive hierarchical resolving assays using nucleic acids (PHRANA) uses a progression of diagnostic genomic loci that are initially highly stable but with low resolution and, ultimately, very unstable but with high resolution. This approach mitigates the need for data weighting and provides both a deeply rooted phylogenetic hypothesis and high resolution discrimination among closely related isolates.

Cyclin mRNA stability does not vary during the cell cycle.

The cyclins are tightly regulated elements governing eukaryotic cell cycle progression by means of sequential activation-inactivation of cyclin-dependent kinases. In one manifestation of this regulation, the mRNA levels of several cyclin genes oscillate during the cycle in mammalian cells. Such cycle-dependent fluctuations in transcript levels could result from changes not only in rates of transcription, but also in mRNA stability. Here we used a new, minimally-disturbing method for producing multi-cycle synchronous growth of human MOLT-4 cells, in combination with quantitative real-time RT-PCR, to compare cell cycle-dependent transcript levels and half-lives of cyclin A2, B1, D3, E and PCNA mRNAs. While all mRNA levels except cyclin D3 varied in the cycle, there were no apparent variations in message half-lives. This differs from several previous reports of dramatic fluctuations in the stabilities of cyclin mRNAs, and infers that fluctuations in cyclin mRNA transcript levels during the MOLT-4 cell cycle are not due to variations in half-lives. The discrepancy in mRNA stability determinations could be due to differences in cell types or synchronization methods, but our findings may be representative of mRNA processing in the cycle of cells in unstressed steady-state growth.