Genetic comparison of Bacillus anthracis and its close relatives using amplified fragment length polymorphism and polymerase chain reaction analysis

Amplified fragment length polymorphism (AFLP) analysis allows a rapid, relatively simple analysis of a large portion of a microbial genome, providing information about the species and its phylogenetic relationship to other microbes (Vos et al. 1995). The method simply surveys the genome for length and sequence polymorphisms. The AFLP pattern identified can be used for comparison to the genomes of other species. Unlike other methods, it does not rely on analysis of a single genetic locus that may bias the interpretation of results and does not require any prior knowledge of the targeted organism. Moreover, a standard set of reagents can be applied to any species without using species-specific information or molecular probes. We are using AFLP analysis to rapidly identify different bacterial species. A comparison of AFLP profiles generated from a large battery of Bacillus anthracis strains shows very little variability among different isolates (Keim et al. 1997). By contrast, there is a significant difference between AFLP profiles generated for any B. anthracis strain and even the most closely related Bacillus species. Sufficient variability is apparent among all known microbial species to allow phylogenetic analysis based on large numbers of genetically unlinked loci. These striking differences among AFLP profiles allow unambiguous identification of previously identified species and phylogenetic placement of newly characterized isolates relative to known species based on a large number of independent genetic loci. Data generated thus far show that the method provides phylogenetic analyses that are consistent with other widely accepted phylogenetic methods. However, AFLP analysis provides a more detailed analysis of the targets and samples a much larger portion of the genome. Consequently, it provides an inexpensive, rapid means of characterizing microbial isolates to further differentiate among strains and closely related microbial species. Such information cannot be rapidly generated by other means. AFLP sample analysis quickly generates a very large amount of molecular information about microbial genomes. However, this information cannot be analysed rapidly using manual methods. We are developing a large archive of electronic AFLP signatures that is being used to identify isolates collected from medical, veterinary, forensic and environmental samples. We are also developing the computational packages necessary to rapidly and unambiguously analyse the AFLP profiles and conduct a phylogenetic comparison of these data relative to information already in our database. We will use this archive and the associated algorithms to determine the species identity of previously uncharacterized isolates and place them phylogenetically relative to other microbes based on their AFLP signatures. This study provides significant new information about microbes with environmental, veterinary and medical significance. This information can be used in further studies to understand the relationships among these species and the factors that distinguish them from one another. It should also allow the identification of unique factors that contribute to important microbial traits, including pathogenicity and virulence. We are also using AFLP data to identify, isolate and sequence DNA fragments that are unique to particular microbial species and strains. The fragment patterns and sequence information provide insights into the complexity and organization of bacterial genomes relative to one another. They also provide the information necessary for the development of species-specific polymerase chain reaction primers that can be used to interrogate complex samples for the presence of B. anthracis, other microbial pathogens or their remnants.     

Estimating plant abundance using inflated beta distributions: Applied learnings from a lichen–caribou ecosystem

Quantifying abundance and distribution of plant species can be difficult because data are often inflated with zero values due to rarity or absence from many ecosystems. Terrestrial fruticose lichens (Cladonia and Cetraria spp.) occupy a narrow ecological niche and have been linked to the diets of declining caribou and reindeer populations (Rangifer tarandus) across their global distribution, and conditions related to their abundance and distribution are not well understood. We attempted to measure effects related to the occupancy and abundance of terrestrial fruticose lichens by sampling and simultaneously modeling two discrete conditions: absence and abundance. We sampled the proportion cover of terrestrial lichens at 438 vegetation plots, including 98 plots having zero lichens. A zero‐inflated beta regression model was employed to simultaneously estimate both the absence and the proportion cover of terrestrial fruticose lichens using fine resolution satellite imagery and light detection and ranging (LiDAR) derived covariates. The probability of lichen absence significantly increased with shallower groundwater, taller vegetation, and increased Sphagnum moss cover. Vegetation productivity, Sphagnum moss cover, and seasonal changes in photosynthetic capacity were negatively related to the abundances of terrestrial lichens. Inflated beta regression reliably estimated the abundance of terrestrial lichens (R² = .74) which was interpolated on a map at fine resolution across a caribou range to support ecological conservation and reclamation. Results demonstrate that sampling for and simultaneously estimating both occupancy and abundance offer a powerful approach to improve statistical estimation and expand ecological inference in an applied setting. Learnings are broadly applicable to studying species that are rare, occupy narrow niches, or where the response variable is a proportion value containing zero or one, which is typical of vegetation cover data.     

The genome and variation of Bacillus anthracis

The Bacillus anthracis genome reflects its close genetic ties to Bacillus cereus and Bacillus thuringiensis but has been shaped by its own unique biology and evolutionary forces. The genome is comprised of a chromosome and two large virulence plasmids, pXO1 and pXO2. The chromosome is mostly co-linear among B. anthracis strains and even with the closest near neighbor strains. An exception to this pattern has been observed in a large inversion in an attenuated strain suggesting that chromosome co-linearity is important to the natural biology of this pathogen. In general, there are few polymorphic nucleotides among B. anthracis strains reflecting the short evolutionary time since its derivation from a B. cereus-like ancestor. The exceptions to this lack of diversity are the variable number tandem repeat (VNTR) loci that exist in genic and non genic regions of the chromosome and both plasmids. Their variation is associated with high mutability that is driven by rapid insertion and deletion of the repeats within an array. A notable example is found in the vrrC locus which is homologous to known DNA translocase genes from other bacteria.     

The complete genome sequence of Bacillus anthracis Ames "Ancestor"

The pathogenic bacterium Bacillus anthracis has become the subject of intense study as a result of its use in a bioterrorism attack in the United States in September and October 2001. Previous studies suggested that B. anthracis Ames Ancestor, the original Ames fully virulent plasmid-containing isolate, was the ideal reference. This study describes the complete genome sequence of that original isolate, derived from a sample kept in cold storage since 1981.     

Structure and function of the abasic site specificity pocket of an AP endonuclease from Archaeoglobus fulgidus

The major AP endonuclease in Escherichia coli Exonuclease III (ExoIII) is frequently used in gene technology due to its strong exonucleolytic activity. A thermostabilized variant of ExoIII or a homologous enzyme from thermophilic organisms could be most useful for further applications. For this purpose we characterized a nuclease from the hyperthermophilic archaeon Archaeoglobus fulgidus (Af_Exo), which shares 33% overall sequence identity and 55% similarity to ExoIII. The gene coding for this thermostable enzyme was cloned and expressed in E. coli. The purified protein shows a strong Mg²⁺-dependent nicking activity at AP-sites, nicking of undamaged double-stranded (ds) DNA and a weak exonucleolytic activity. A V217G variant of the enzyme was crystallized with decamer ds-DNA molecule, and the three-dimensional structure was determined to 1.7 Å resolution. Besides our goal to find or produce a thermostable exonuclease, the structural and catalytic data of Af_Exo and a series of mutant proteins, based on the crystal structure, provide new insight into the mechanism of abasic site recognition and repair. Each of the hydrophobic residues Phe 200, Trp 215 and Val 217, forming a binding pocket for the abasic deoxyribose in Af_Exo, were mutated to glycine or serine. By expanding the size of the binding pocket the unspecific endonucleolytic activity is increased. Thus, size and flexibility of the mostly hydrophobic binding pocket have a significant influence on AP-site specificity. We suggest that its tight fitting to the flipped-out deoxyribose allows for a preferred competent binding of abasic sites. In a larger or more flexible pocket however, intact nucleotides more easily bind in a catalytically competent conformation, resulting in loss of specificity. Moreover, with mutations of Phe 200 and Trp 215 we induced a strong exonucleolytic activity on undamaged DNA.     

Population differences in host immune factors may influence survival of Gunnison's prairie dogs (Cynomys gunnisoni) during plague outbreaks

Over the past 40 yr, epizootics of plague (Yersinia pestis) in northern Arizona have reduced populations of the Gunnison's prairie dog (Cynomys gunnisoni), with the exception of a large population found in the Aubrey Valley (AV). To examine potential mechanisms accounting for their survival, we collected prairie dog serum samples in 2005-2006 from AV and a neighboring population near Seligman (SE), Arizona. We quantified gene expression at 58 diverse immune proteins using a multiplexed enzyme-linked immunosorbent assay panel. We found a subset of proteins important in coagulation and inflammation (tissue factor [TF], calbindin [Cal], and thrombopoietin [TPO]) and T-cell responses (CD40L and CD40) that were present in AV at levels two to eight times greater than SE. These results suggest that AV and SE animals might differ in their ability to mount an immune response.     

TaqMan Real-Time PCR Assays for Single-Nucleotide Polymorphisms Which Identify Francisella tularensis and Its Subspecies and Subpopulations

Francisella tularensis, the etiologic agent of tularemia and a Class A Select Agent, is divided into three subspecies and multiple subpopulations that differ in virulence and geographic distribution. Given these differences, there is a need to rapidly and accurately determine if a strain is F. tularensis and, if it is, assign it to subspecies and subpopulation. We designed TaqMan real-time PCR genotyping assays using eleven single nucleotide polymorphisms (SNPs) that were potentially specific to closely related groups within the genus Francisella, including numerous subpopulations within F. tularensis species. We performed extensive validation studies to test the specificity of these SNPs to particular populations by screening the assays across a set of 565 genetically and geographically diverse F. tularensis isolates and an additional 21 genetic near-neighbor (outgroup) isolates. All eleven assays correctly determined the genetic groups of all 565 F. tularensis isolates. One assay differentiates F. tularensis, F. novicida, and F. hispaniensis from the more genetically distant F. philomiragia and Francisella-like endosymbionts. Another assay differentiates F. tularensis isolates from near neighbors. The remaining nine assays classify F. tularensis-confirmed isolates into F. tularensis subspecies and subpopulations. The genotyping accuracy of these nine assays diminished when tested on outgroup isolates (i.e. non F. tularensis), therefore a hierarchical approach of assay usage is recommended wherein the F. tularensis-specific assay is used before the nine downstream assays. Among F. tularensis isolates, all eleven assays were highly sensitive, consistently amplifying very low concentrations of DNA. Altogether, these eleven TaqMan real-time PCR assays represent a highly accurate, rapid, and sensitive means of identifying the species, subspecies, and subpopulation of any F. tularensis isolate if used in a step-wise hierarchical scheme. These assays would be very useful in clinical, epidemiological, and/or forensic investigations involving F. tularensis.