Rapid, accurate, and sensitive detection of biothreat agents requires a broad-spectrum assay capable of discriminating between closely related microbial or viral pathogens. Moreover, in cases where a biological agent release has been identified, forensic analysis demands detailed genetic signature data for accurate strain identification and attribution. To date, nucleic acid sequences have provided the most robust and phylogentically illuminating signature information. Nucleic acid signature sequences are not often linked to genomic or extrachromosomal determinants of virulence, a link that would further facilitate discrimination between pathogens and closely related species. Inextricably coupling genetic determinants of virulence with highly informative nucleic acid signatures would provide a robust means of identifying human, livestock, and agricultural pathogens. By means of example, we present here an overview of two general applications of microarray-based methods for: (1) the identification of candidate virulence factors; and (2) the analysis of genetic polymorphisms that are coupled to Bacillus anthracis virulence factors using an accurate, low cost solid-phase mini-sequencing assay. We show that microarray-based analysis of gene expression can identify potential virulence associated genes for use as candidate signature targets, and, further, that microarray-based single nucleotide polymorphism assays provide a robust platform for the detection and identification of signature sequences in a manner independent of the genetic background in which the signature is embedded. We discuss the strategy as a general approach or pipeline for the discovery of virulence-linked nucleic acid signatures for biothreat agents.

The potential of microarray data is enormous. It allows us to monitor the expression of thousands of genes simultaneously. A common task with microarray is to determine which genes are differentially expressed between two samples obtained under two different conditions. Recently, several statistical methods have been proposed to perform such a task when there are replicate samples under each condition. Two major problems arise with microarray data. The first one is that the number of replicates is very small (usually 2-10), leading to noisy point estimates. As a consequence, traditional statistics that are based on the means and standard deviations, e.g. t-statistic, are not suitable. The second problem is that the number of genes is usually very large (approximately 10,000), and one is faced with an extreme multiple testing problem. Most multiple testing adjustments are relatively conservative, especially when the number of replicates is small. In this paper we present an empirical Bayes analysis that handles both problems very well. Using different parametrizations, we develop four statistics that can be used to test hypotheses about the means and/or variances of the gene expression levels in both one- and two-sample problems. The methods are illustrated using experimental data with prior knowledge. In addition, we present the result of a simulation comparing our methods to well-known statistics and multiple testing adjustments.

BACKGROUND:Complete sequencing and annotation of the 96.2 kb Bacillus anthracis plasmid, pXO2, predicted 85 open reading frames (ORFs). Bacillus cereus and Bacillus thuringiensis isolates that ranged in genomic similarity to B. anthracis, as determined by amplified fragment length polymorphism (AFLP) analysis, were examined by PCR for the presence of sequences similar to 47 pXO2 ORFs.RESULTS:The two most distantly related isolates examined, B. thuringiensis 33679 and B. thuringiensis AWO6, produced the greatest number of ORF sequences similar to pXO2; 10 detected in 33679 and 16 in AWO6. No more than two of the pXO2 ORFs were detected in any one of the remaining isolates. Dot-blot DNA hybridizations between pXO2 ORF fragments and total genomic DNA from AWO6 were consistent with the PCR assay results for this isolate and also revealed nine additional ORFs shared between these two bacteria. Sequences similar to the B. anthracis cap genes or their regulator, acpA, were not detected among any of the examined isolates.CONCLUSIONS:The presence of pXO2 sequences in the other Bacillus isolates did not correlate with genomic relatedness established by AFLP analysis. The presence of pXO2 ORF sequences in other Bacillus species suggests the possibility that certain pXO2 plasmid gene functions may also be present in other closely related bacteria.

ABSTRACT The complete sequencing and annotation of the 181.7-kb Bacillus anthracis virulence plasmid pXO1 predicted 143 genes but could only assign putative functions to 45. Hybridization assays, PCR amplification, and DNA sequencing were used to determine whether pXO1 open reading frame (ORF) sequences were present in other bacilli and more distantly related bacterial genera. Eighteen Bacillus species isolates and four other bacterial species were tested for the presence of 106 pXO1 ORFs. Three ORFs were conserved in most of the bacteria tested. Many of the pXO1 ORFs were detected in closely related Bacillus species, and some were detected only in B. anthracis isolates. Three isolates, Bacillus cereus D-17, B. cereus 43881, and Bacillus thuringiensis 33679, contained sequences that were similar to more than one-half of the pXO1 ORF sequences examined. The majority of the DNA fragments that were amplified by PCR from these organisms had DNA sequences between 80 and 98% similar to that of pXO1. Pulsed-field gel electrophoresis revealed large potential plasmids present in both B. cereus 43881 (341 kb) and B. thuringiensis ATCC 33679 (327 kb) that hybridized with a DNA probe composed of six pXO1 ORFs.