Search

Abstract Background The bacterium Vibrio cholerae causes diarrheal illness and can acquire genetic material leading to multiple drug resistance (MDR). Rapid detection of resistance-conferring mobile genetic elements helps avoid the prescription of ineffective antibiotics for specific strains. Colorimetric loop-mediated isothermal amplification (LAMP) assays provide a rapid and cost-effective means for detection at point-of-care since they do not require specialized equipment, require limited expertise to perform, and can take less than 30 min to perform in resource limited regions. LAMP output is a color change that can be viewed by eye, but it can be difficult to design primer sets, determine target specificity, and interpret subjective color changes. Methods We developed an algorithm for the in silico design and evaluation of LAMP assays within the open-source PCR Signature Erosion Tool (PSET) and a computer vision application for the quantitative analysis of colorimetric outputs. First, Primer3 calculates LAMP primer sequence candidates with settings based on GC-content optimization. Next, PSET aligns the primer sequences of each assay against large sequence databases to calculate sufficient sequence similarity, coverage, and primer arrangement to the intended taxa, ultimately generating a confusion matrix. Finally, we tested assay candidates in the laboratory against synthetic constructs. Results As an example, we generated new LAMP assays targeting drug resistance in V. cholerae and evaluated existing ones from the literature based on in silico target specificity and in vitro testing. Improvements in the design and testing of LAMP assays, with heightened target specificity and a simple analysis platform, increase utility for in-field applications. Overall, 9 of the 16 tested LAMP assays had positive signal through visual and computer vision-based detection methods developed here. Here we show LAMP assays tested on synthetic AMR gene targets for aph(6), varG, floR, qnrVC5, and almG, which allow for resistance to aminoglycosides, penicillins, carbapenems, phenicols, fluoroquinolones, and polymyxins respectively.

ABSTRACT Antigen-based rapid diagnostic tests (Ag-RDTs) provide timely results, are simple to use, and are less expensive than molecular assays. Recent studies suggest that antigen-based testing aligns with virus culture-based results (a proxy of contagiousness at the peak viral phase of illness); however, the performance of Ag-RDTs for newer SARS-CoV-2 variants is unclear. In this study, we (i) assessed the performance of Ag-RDTs and diagnostic antibodies to detect a range of SARS-CoV-2 variants and (ii) determined whether Ag-RDT results correlated with culture positivity. We noted only minor differences in the limit of detection by variant for all assays, and we demonstrated consistent antibody affinity to the N protein among the different variants. We observed moderate to high sensitivity (46.8%–83.9%) for Ag-RDTs when compared to PCR positivity (100%), and all variants were assessed on each assay. Ag-RDT sensitivity and PCR Ct showed an inverse correlation with the detection of viable virus. Collectively, our results demonstrate that commercially available Ag-RDTs offer variable sensitivity compared to PCR, show similar diagnostic validity across variants, and may predict the risk of transmissibility. These findings may be used to support more tailored SARS-CoV-2 isolation strategies, particularly if other studies clarify the direct association between Ag-RDT positivity and transmission risk. The apparent trade-off between sensitivity in the detection of any PCR-positive infection and concordance with infectious virus positivity may also inform new RDT diagnostic development strategies for SARS-CoV-2 and other epidemic respiratory pathogens.IMPORTANCEDespite the availability of vaccines, COVID-19 continues to be a major health concern, and antigen-based rapid diagnostic tests (Ag-RDTs) are commonly used as point-of-care or at-home diagnostic tests. In this study, we evaluated the performance of two commercially available Ag-RDTs and a research Ag-RDT to detect multiple SARS-CoV-2 variants using upper respiratory tract swab samples from clinical COVID-19 cases. Furthermore, we determined whether Ag-RDT results correlated with culture positivity, a potential proxy of viral transmissibility. Our results have important implications to inform future testing and response strategies during periods of high COVID-19 transmission with new variants.

Abstract BackgroundThe COVID-19 pandemic highlighted the need for pathogen surveillance systems to augment both early warning and outbreak monitoring/control efforts. Community wastewater samples provide a rapid and accurate source of environmental surveillance data to complement direct patient sampling. Due to its global presence and critical missions, the US military is a leader in global pandemic preparedness efforts. Clinical testing for COVID-19 on US Air Force (USAF) bases (AFBs) was effective but costly with respect to direct monetary costs and indirect costs due to lost time. To remain operating at peak capacity, such bases sought a more passive surveillance option and piloted wastewater surveillance (WWS) at 17 AFBs to demonstrate feasibility, safety, utility, and cost-effectiveness from May 2021 to January 2022. ObjectiveWe model the costs of a wastewater program for pathogens of public health concern within the specific context of US military installations using assumptions based on the results of the USAF and Joint Program Executive Office for Chemical, Biological, Radiological and Nuclear Defense pilot program. The objective was to determine the cost of deploying WWS to all AFBs relative to clinical swab testing surveillance regimes. MethodsA WWS cost projection model was built based on subject matter expert input and actual costs incurred during the WWS pilot program at USAF AFBs. Several SARS-CoV-2 circulation scenarios were considered, and the costs of both WWS and clinical swab testing were projected. Analysis was conducted to determine the break-even point and how a reduction in swab testing could unlock funds to enable WWS to occur in parallel. ResultsOur model confirmed that WWS is complementary and highly cost-effective when compared to existing alternative forms of biosurveillance. We found that the cost of WWS was between US $10.5-$18.5 million less expensive annually in direct costs as compared to clinical swab testing surveillance. When the indirect cost of lost work was incorporated, including lost work associated with required clinical swab testing, we estimated that over two-thirds of clinical swab testing could be maintained with no additional costs upon implementation of WWS. ConclusionsOur results support the adoption of WWS across US military installations as part of a more comprehensive and early warning system that will enable adaptive monitoring during disease outbreaks in a more cost-effective manner than swab testing alone.

Abstract Objective In this study, we sought to determine the types and prevalence of antimicrobial resistance determinants (ARDs) in Burkholderia spp. strains using the Antimicrobial Resistance Determinant Microarray (ARDM). Results Whole genome amplicons from 22 B. mallei (BM) and 37 B. pseudomallei (BP) isolates were tested for > 500 ARDs using ARDM v.3.1. ARDM detected the following Burkholderia spp.-derived genes, aac(6), bla BP/MBL-3, blaA BPS, penA-BP, and qacE, in both BM and BP while bla BP/MBL-1, macB, bla OXA-42/43 and penA-BC were observed in BP only. The method of denaturing template for whole genome amplification greatly affected the numbers and types of genes detected by the ARDM. Bla TEM was detected in nearly a third of BM and BP amplicons derived from thermally, but not chemically denatured templates. Bla TEM results were confirmed by PCR, with 81% concordance between methods. Sequences from 414-nt PCR amplicons (13 preparations) were 100% identical to the Klebsiella pneumoniae reference gene. Although bla TEM sequences have been observed in B. glumae, B. cepacia, and other undefined Burkholderia strains, this is the first report of such sequences in BM/BP/B. thailandensis (BT) clade. These results highlight the importance of sample preparation in achieving adequate genome coverage in methods requiring untargeted amplification before analysis.

Bacteriophages such as γ and AP50c have been shown to infect strains of Bacillus anthracis with high specificity, and this feature has been exploited in the development of bacterial detection assays. To better understand the emergence of phage resistance, and thus the potential failure of such assays, it is important to identify the host and phage receptors necessary for attachment and entry. Using genetic approaches, the bacterial receptors of AP50c and γ have been identified as sap and GamR, respectively. A second AP50c-like phage, Wip1, also appears to use sap as a receptor. In parallel with this work, the cognate phage-encoded receptor binding proteins (RBPs) have also been identified (Gp14 for γ, P28 for AP50c, and P23 for Wip1); however, the strength of evidence supporting these protein–protein interactions varies, necessitating additional investigation. Here, we present genetic evidence further supporting the interaction between sap and the RBPs of AP50c and Wip1 using fluorescently tagged proteins and a panel of B. anthracis mutants. These results showed that the deletion of the sap gene, as well as the deletion of csaB, whose encoded protein anchors sap to the bacterial S-layer, resulted in the loss of RBP binding. Binding could then be rescued by expressing these genes in trans. We also found that the RBP of the γ-like prophage λBa03 relied on csaB activity for binding, possibly by a different mechanism. RBPλBa03 binding to B. anthracis cells was also unique in that it was not ablated by heat inactivation of vegetative cells, suggesting that its receptor is still functional following incubation at 98°C. These results extend our understanding of the diverse attachment and entry strategies used by B. anthracis phages, enabling future assay development.

Many SARS-CoV-2 variants have emerged during the course of the COVID-19 pandemic. These variants have acquired mutations conferring phenotypes such as increased transmissibility or virulence, or causing diagnostic, therapeutic, or immune escape. Detection of Alpha and the majority of Omicron sublineages by PCR relied on the so-called S gene target failure due to the deletion of six nucleotides coding for amino acids 69–70 in the spike (S) protein. Detection of hallmark mutations in other variants present in samples relied on whole genome sequencing. However, whole genome sequencing as a diagnostic tool is still in its infancy due to geographic inequities in sequencing capabilities, higher cost compared to other molecular assays, longer turnaround time from sample to result, and technical challenges associated with producing complete genome sequences from samples that have low viral load and/or high background. Hence, there is a need for rapid genotyping assays. In order to rapidly generate information on the presence of a variant in a given sample, we have created a panel of four triplex RT-qPCR assays targeting 12 mutations to detect and differentiate all five variants of concern: Alpha, Beta, Gamma, Delta, and Omicron. We also developed an expanded pentaplex assay that can reliably distinguish among the major sublineages (BA.1–BA.5) of Omicron. In silico, analytical and clinical testing of the variant panel indicate that the assays exhibit high sensitivity and specificity. This panel can help fulfill the need for rapid identification of variants in samples, leading to quick decision making with respect to public health measures, as well as treatment options for individuals. Compared to sequencing, these genotyping PCR assays allow much faster turn-around time from sample to results—just a couple hours instead of days or weeks.

Polymerase chain reaction (PCR) remains the gold standard in disease diagnostics due to its extreme sensitivity and specificity. However, PCR tests are expensive and complex, require skilled personnel and specialized equipment to conduct the tests, and have long turnaround times. On the other hand, lateral flow immunoassay-based antigen tests are rapid, relatively inexpensive, and can be performed by untrained personnel at the point of care or even in the home. However, rapid antigen tests are less sensitive than PCR since they lack the inherent target amplification of PCR. It has been argued that rapid antigen tests are better indicators of infection in public health decision-making processes to test, trace, and isolate infected people to curtail further transmission. Hence, there is a critical need to increase the sensitivity of rapid antigen tests and create innovative solutions to achieve that goal. Herein, we report the development of a low-cost diagnostic platform, enabling rapid detection of SARS-CoV-2 under field or at-home conditions. This platform (Halo™) is a small, highly accurate, consumer-friendly diagnostic reader paired with fluorescently labeled lateral flow assays and custom software for collection and reporting of results. The focus of this study is to compare the analytical performance of HaloTM against comparable tests that use either colloidal gold nanoparticles or fluorescence-based reporters in simulated nasal matrix and not in clinical samples. Live virus data has demonstrated limit of detection performance of 1.9 TCID50/test in simulated nasal matrix for the delta variant, suggesting that single-assay detection of asymptomatic SARS-CoV-2 infections may be feasible. Performance of the system against all tested SARS CoV-2 virus variants showed comparable sensitivities indicating mutations in SARS-CoV-2 variants do not negatively impact the assay.

Real-time reverse transcription polymerase chain reaction (RT-PCR) assays are the most widely used molecular tests for the detection of SARS-CoV-2 and diagnosis of COVID-19 in clinical samples. PCR assays target unique genomic RNA regions to identify SARS-CoV-2 with high sensitivity and specificity. In general, assay development incorporates the whole genome sequences available at design time to be inclusive of all target species and exclusive of near neighbors. However, rapid accumulation of mutations in viral genomes during sustained growth in the population can result in signature erosion and assay failures, creating situational blind spots during a pandemic. In this study, we analyzed the signatures of 43 PCR assays distributed across the genome against over 1.6 million SARS-CoV-2 sequences. We present evidence of significant signature erosion emerging in just two assays due to mutations, while adequate sequence identity was preserved in the other 41 assays. Failure of more than one assay against a given variant sequence was rare and mostly occurred in the two assays noted to have signature erosion. Assays tended to be designed in regions with statistically higher mutations rates. in silico analyses over time can provide insights into mutation trends and alert users to the emergence of novel variants that are present in the population at low proportions before they become dominant. Such routine assessment can also potentially highlight false negatives in test samples that may be indicative of mutations having functional consequences in the form of vaccine and therapeutic failures. This study highlights the importance of whole genome sequencing and expanded real-time monitoring of diagnostic PCR assays during a pandemic.

Abstract Background The state-of-the-art in nucleic acid based biodetection continues to be polymerase chain reaction (PCR), and many real-time PCR assays targeting biodefense pathogens for biosurveillance are in widespread use. These assays are predominantly singleplex; i.e. one assay tests for the presence of one target, found in a single organism, one sample at a time. Due to the intrinsic limitations of such tests, there exists a critical need for high-throughput multiplex assays to reduce the time and cost incurred when screening multiple targets, in multiple pathogens, and in multiple samples. Such assays allow users to make an actionable call while maximizing the utility of the small volumes of test samples. Unfortunately, current multiplex real-time PCR assays are limited in the number of targets that can be probed simultaneously due to the availability of fluorescence channels in real-time PCR instruments. Results To address this gap, we developed a pipeline in which the amplicons produced by a 14-plex end-point PCR assay using spiked samples were subsequently sequenced using Nanopore technology. We used bar codes to sequence multiple samples simultaneously, leading to the generation and subsequent analysis of sequence data resulting from a short sequencing run time (

The revelation in May 2015 of the shipment of γ irradiation–inactivated wild-type Bacillus anthracis spore preparations containing a small number of live spores raised concern about the safety and security of these materials. The finding also raised doubts about the validity of the protocols and procedures used to prepare them. Such inactivated reference materials were used as positive controls in assays to detect suspected B. anthracis in samples because live agent cannot be shipped for use in field settings, in improvement of currently deployed detection methods or development of new methods, or for quality assurance and training activities. Hence, risk-mitigated B. anthracis strains are needed to fulfill these requirements. We constructed a genetically inactivated or attenuated strain containing relevant molecular assay targets and tested to compare assay performance using this strain to the historical data obtained using irradiation-inactivated virulent spores.

In recent years, an increasing diversity of species has been recognized within the family Francisellaceae. Unfortunately, novel isolates are sometimes misnamed in initial publications or multiple sources propose different nomenclature for genetically highly similar isolates. Thus, unstructured and occasionally incorrect information can lead to confusion in the scientific community. Historically, detection of Francisella tularensis in environmental samples has been challenging due to the considerable and unknown genetic diversity within the family, which can result in false positive results. We have assembled a comprehensive collection of genome sequences representing most known Francisellaceae species/strains and restructured them according to a taxonomy that is based on phylogenetic structure. From this structured dataset, we identified a small number of genomic regions unique to F. tularensis that are putatively suitable for specific detection of this pathogen in environmental samples. We designed and validated specific PCR assays based on these genetic regions that can be used for the detection of F. tularensis in environmental samples, such as water and air filters.

Accurate pathogen detection and diagnosis is paramount in clinical success of treating patients. There are two general paradigms in pathogen detection: molecular and immuno-based, and phage-based detection is a third emerging paradigm due to its sensitivity and selectivity. Molecular detection methods look for genetic material specific for a given pathogen in a sample usually by polymerase chain reaction (PCR). Immuno-methods look at the pathogen components (antigens) by antibodies raised against that pathogen specific antigens. There are different variations and products based on these two paradigms with advantages and disadvantages. The third paradigm at least for bacterial pathogen detection entails bacteriophages specific for a given bacterium. Sensitivity and specificity are the two key parameters in any pathogen detection system. By their very nature, bacteriophages afford the best sensitivity for bacterial detection. Bacteria and bacteriophages form the predator-prey pair in the evolutionary arms race and has coevolved over time to acquire the exquisite specificity of the pair, in some instances at the strain level. This specificity has been exploited for diagnostic purposes of various pathogens of concern in clinical and other settings. Many recent reviews focus on phage-based detection and sensor technologies. In this review, we focus on a very special group of pathogens that are of concern in biodefense because of their potential misuse in bioterrorism and their extremely virulent nature and as such fall under the Centers for Disease and Prevention (CDC) Category A pathogen list. We describe the currently available phage methods that are based on the usual modalities of detection from culture, to molecular and immuno- and fluorescent methods. We further highlight the gaps and the needs for more modern technologies and sensors drawing from technologies existing for detection and surveillance of other pathogens of clinical relevance.

Chromosomal inheritance in bacteria usually entails bidirectional replication of a single chromosome from a single origin into two copies and subsequent partitioning of one copy each into daughter cells upon cell division. However, the human pathogen Vibrio cholerae and other Vibrionaceae harbor two chromosomes, a large Chr1 and a small Chr2. Chr1 and Chr2 have different origins, an oriC-type origin and a P1 plasmid-type origin, respectively, driving the replication of respective chromosomes. Recently, we described naturally occurring exceptions to the two-chromosome rule of Vibrionaceae: i.e., Chr1 and Chr2 fused single chromosome V. cholerae strains, NSCV1 and NSCV2, in which both origins of replication are present. Using NSCV1 and NSCV2, here we tested whether two types of origins of replication can function simultaneously on the same chromosome or one or the other origin is silenced. We found that in NSCV1, both origins are active whereas in NSCV2 ori2 is silenced despite the fact that it is functional in an isolated context. The ori2 activity appears to be primarily determined by the copy number of the triggering site, crtS, which in turn is determined by its location with respect to ori1 and ori2 on the fused chromosome.

Genome sequence analyses of the 2014 Ebola Virus (EBOV) isolates revealed a potential problem with the diagnostic assays currently in use; i.e., drifting genomic profiles of the virus may affect the sensitivity or even produce false-negative results. We evaluated signature erosion in ebolavirus molecular assays using an in silico approach and found frequent potential false-negative and false-positive results. We further empirically evaluated many EBOV assays, under real time PCR conditions using EBOV Kikwit (1995) and Makona (2014) RNA templates. These results revealed differences in performance between assays but were comparable between the old and new EBOV templates. Using a whole genome approach and a novel algorithm, termed BioVelocity, we identified new signatures that are unique to each of EBOV, Sudan virus (SUDV), and Reston virus (RESTV). Interestingly, many of the current assay signatures do not fall within these regions, indicating a potential drawback in the past assay design strategies. The new signatures identified in this study may be evaluated with real-time reverse transcription PCR (rRT-PCR) assay development and validation. In addition, we discuss regulatory implications and timely availability to impact a rapidly evolving outbreak using existing but perhaps less than optimal assays versus redesign these assays for addressing genomic changes.

The genetic make-up of most bacteria is encoded in a single chromosome while about 10% have more than one chromosome. Among these, Vibrio cholerae, with two chromosomes, has served as a model system to study various aspects of chromosome maintenance, mainly replication, and faithful partitioning of multipartite genomes. Here, we describe the genomic characterization of strains that are an exception to the two chromosome rules: naturally occurring single-chromosome V. cholerae. Whole genome sequence analyses of NSCV1 and NSCV2 (natural single-chromosome vibrio) revealed that the Chr1 and Chr2 fusion junctions contain prophages, IS elements, and direct repeats, in addition to large-scale chromosomal rearrangements such as inversions, insertions, and long tandem repeats elsewhere in the chromosome compared to prototypical two chromosome V. cholerae genomes. Many of the known cholera virulence factors are absent. The two origins of replication and associated genes are generally intact with synonymous mutations in some genes, as are recA and mismatch repair (MMR) genes dam, mutH, and mutL; MutS function is probably impaired in NSCV2. These strains are ideal tools for studying mechanistic aspects of maintenance of chromosomes with multiple origins and other rearrangements and the biological, functional, and evolutionary significance of multipartite genome architecture in general.

Historically, cholera outbreaks have been linked to V. cholerae O1 serogroup strains or its derivatives of the O37 and O139 serogroups. A genomic study on the 2010 Haiti cholera outbreak strains highlighted the putative role of non O1/non-O139 V. cholerae in causing cholera and the lack of genomic sequences of such strains from around the world. Here we address these gaps by scanning a global collection of V. cholerae strains as a first step towards understanding the population genetic diversity and epidemic potential of non O1/non-O139 strains. Whole Genome Mapping (Optical Mapping) based bar coding produces a high resolution, ordered restriction map, depicting a complete view of the unique chromosomal architecture of an organism. To assess the genomic diversity of non-O1/non-O139 V. cholerae, we applied a Whole Genome Mapping strategy on a well-defined and geographically and temporally diverse strain collection, the Sakazaki serogroup type strains. Whole Genome Map data on 91 of the 206 serogroup type strains support the hypothesis that V. cholerae has an unprecedented genetic and genomic structural diversity. Interestingly, we discovered chromosomal fusions in two unusual strains that possess a single chromosome instead of the two chromosomes usually found in V. cholerae. We also found pervasive chromosomal rearrangements such as duplications and indels in many strains. The majority of Vibrio genome sequences currently in public databases are unfinished draft sequences. The Whole Genome Mapping approach presented here enables rapid screening of large strain collections to capture genomic complexities that would not have been otherwise revealed by unfinished draft genome sequencing and thus aids in assembling and finishing draft sequences of complex genomes. Furthermore, Whole Genome Mapping allows for prediction of novel V. cholerae non-O1/non-O139 strains that may have the potential to cause future cholera outbreaks.

ABSTRACT Thanks to high-throughput sequencing technologies, genome sequencing has become a common component in nearly all aspects of viral research; thus, we are experiencing an explosion in both the number of available genome sequences and the number of institutions producing such data. However, there are currently no common standards used to convey the quality, and therefore utility, of these various genome sequences. Here, we propose five “standard” categories that encompass all stages of viral genome finishing, and we define them using simple criteria that are agnostic to the technology used for sequencing. We also provide genome finishing recommendations for various downstream applications, keeping in mind the cost-benefit trade-offs associated with different levels of finishing. Our goal is to define a common vocabulary that will allow comparison of genome quality across different research groups, sequencing platforms, and assembly techniques.

In Staphylococcus aureus, the disease impact of chromosomally integrated prophages on virulence is well described. However, the existence of extra-chromosomal prophages, both plasmidial and episomal, remains obscure. Despite the recent explosion in bacterial and bacteriophage genomic sequencing, studies have failed to specifically focus on extra-chromosomal elements. We selectively enriched and sequenced extra-chromosomal DNA from S. aureus isolates using Roche-454 technology and uncovered evidence for the widespread distribution of multiple extra-chromosomal prophages (ExPΦs) throughout both antibiotic-sensitive and -resistant strains. We completely sequenced one such element comprised of a 43.8 kbp, circular ExPΦ (designated ФBU01) from a vancomycin-intermediate S. aureus (VISA) strain. Assembly and annotation of ФBU01 revealed a number of putative virulence determinants encoded within a bacteriophage immune evasion cluster (IEC). Our identification of several potential ExPΦs and mobile genetic elements (MGEs) also revealed numerous putative virulence factors and antibiotic resistance genes. We describe here a previously unidentified level of genetic diversity of stealth extra-chromosomal elements in S. aureus, including phages with a larger presence outside the chromosome that likely play a prominent role in pathogenesis and strain diversity driven by horizontal gene transfer (HGT).

In May of 2011, an enteroaggregative Escherichia coli O104:H4 strain that had acquired a Shiga toxin 2-converting phage caused a large outbreak of bloody diarrhea in Europe which was notable for its high prevalence of hemolytic uremic syndrome cases. Several studies have described the genomic inventory and phylogenies of strains associated with the outbreak and a collection of historical E. coli O104:H4 isolates using draft genome assemblies. We present the complete, closed genome sequences of an isolate from the 2011 outbreak (2011C-3493) and two isolates from cases of bloody diarrhea that occurred in the Republic of Georgia in 2009 (2009EL-2050 and 2009EL-2071). Comparative genome analysis indicates that, while the Georgian strains are the nearest neighbors to the 2011 outbreak isolates sequenced to date, structural and nucleotide-level differences are evident in the Stx2 phage genomes, the mer/tet antibiotic resistance island, and in the prophage and plasmid profiles of the strains, including a previously undescribed plasmid with homology to the pMT virulence plasmid of Yersinia pestis. In addition, multiphenotype analysis showed that 2009EL-2071 possessed higher resistance to polymyxin and membrane-disrupting agents. Finally, we show evidence by electron microscopy of the presence of a common phage morphotype among the European and Georgian strains and a second phage morphotype among the Georgian strains. The presence of at least two stx2 phage genotypes in host genetic backgrounds that may derive from a recent common ancestor of the 2011 outbreak isolates indicates that the emergence of stx2 phage-containing E. coli O104:H4 strains probably occurred more than once, or that the current outbreak isolates may be the result of a recent transfer of a new stx2 phage element into a pre-existing stx2-positive genetic background.

BackgroundDespite the decades-long use of Bacillus atrophaeus var. globigii (BG) as a simulant for biological warfare (BW) agents, knowledge of its genome composition is limited. Furthermore, the ability to differentiate signatures of deliberate adaptation and selection from natural variation is lacking for most bacterial agents. We characterized a lineage of BGwith a long history of use as a simulant for BW operations, focusing on classical bacteriological markers, metabolic profiling and whole-genome shotgun sequencing (WGS).ResultsArchival strains and two "present day" type strains were compared to simulant strains on different laboratory media. Several of the samples produced multiple colony morphotypes that differed from that of an archival isolate. To trace the microevolutionary history of these isolates, we obtained WGS data for several archival and present-day strains and morphotypes. Bacillus-wide phylogenetic analysis identified B. subtilis as the nearest neighbor to B. atrophaeus. The genome of B. atrophaeus is, on average, 86% identical to B. subtilis on the nucleotide level. WGS of variants revealed that several strains were mixed but highly related populations and uncovered a progressive accumulation of mutations among the "military" isolates. Metabolic profiling and microscopic examination of bacterial cultures revealed enhanced growth of "military" isolates on lactate-containing media, and showed that the "military" strains exhibited a hypersporulating phenotype.ConclusionsOur analysis revealed the genomic and phenotypic signatures of strain adaptation and deliberate selection for traits that were desirable in a simulant organism. Together, these results demonstrate the power of whole-genome and modern systems-level approaches to characterize microbial lineages to develop and validate forensic markers for strain discrimination and reveal signatures of deliberate adaptation.

BACKGROUND: The anthrax letter attacks of 2001 highlighted the need for rapid identification of biothreat agents not only for epidemiological surveillance of the intentional outbreak but also for implementing appropriate countermeasures, such as antibiotic treatment, in a timely manner to prevent further casualties. It is clear from the 2001 cases that survival may be markedly improved by administration of antimicrobial therapy during the early symptomatic phase of the illness; i.e., within 3 days of appearance of symptoms. Microbiological detection methods are feasible only for organisms that can be cultured in vitro and cannot detect all genetic modifications with the exception of antibiotic resistance. Currently available immuno or nucleic acid-based rapid detection assays utilize known, organism-specific proteins or genomic DNA signatures respectively. Hence, these assays lack the ability to detect novel natural variations or intentional genetic modifications that circumvent the targets of the detection assays or in the case of a biological attack using an antibiotic resistant or virulence enhanced Bacillus anthracis, to advise on therapeutic treatments. METHODOLOGY/PRINCIPAL FINDINGS: We show here that the Roche 454-based pyrosequencing can generate whole genome draft sequences of deep and broad enough coverage of a bacterial genome in less than 24 hours. Furthermore, using the unfinished draft sequences, we demonstrate that unbiased identification of known as well as heretofore-unreported genetic modifications that include indels and single nucleotide polymorphisms conferring antibiotic and phage resistances is feasible within the next 12 hours. CONCLUSIONS/SIGNIFICANCE: Second generation sequencing technologies have paved the way for sequence-based rapid identification of both known and previously undocumented genetic modifications in cultured, conventional and newly emerging biothreat agents. Our findings have significant implications in the context of whole genome sequencing-based routine clinical diagnostics as well as epidemiological surveillance of natural disease outbreaks caused by bacterial and viral agents.

BACKGROUND: Despite the global threat caused by arthropod-borne viruses, there is not an efficient method for screening vector populations to detect novel viral sequences. Current viral detection and surveillance methods based on culture can be costly and time consuming and are predicated on prior knowledge of the etiologic agent, as they rely on specific oligonucleotide primers or antibodies. Therefore, these techniques may be unsuitable for situations when the causative agent of an outbreak is unknown. METHODOLOGY/PRINCIPAL FINDINGS: In this study we explored the use of high-throughput pyrosequencing for surveillance of arthropod-borne RNA viruses. Dengue virus, a member of the positive strand RNA Flavivirus family that is transmitted by several members of the Aedes genus of mosquitoes, was used as a model. Aedes aegypti mosquitoes experimentally infected with dengue virus type 1 (DENV-1) were pooled with noninfected mosquitoes to simulate samples derived from ongoing arbovirus surveillance programs. Using random-primed methods, total RNA was reverse-transcribed and resulting cDNA subjected to 454 pyrosequencing. CONCLUSIONS/SIGNIFICANCE: In two types of samples, one with 5 adult mosquitoes infected with DENV-1- and the other with 1 DENV-1 infected mosquito and 4 noninfected mosquitoes, we identified DENV-1 DNA sequences. DENV-1 sequences were not detected in an uninfected control pool of 5 adult mosquitoes. We calculated the proportion of the Ae. aegypti metagenome contributed by each infecting Dengue virus genome (p(IP)), which ranged from 2.75×10(-8) to 1.08×10(-7). DENV-1 RNA was sufficiently concentrated in the mosquito that its detection was feasible using current high-throughput sequencing instrumentation. We also identified some of the components of the mosquito microflora on the basis of the sequence of expressed RNA. This included members of the bacterial genera Pirellula and Asaia, various fungi, and a potentially uncharacterized mycovirus.

Complete sequences of 9.5-kb pPCP1 plasmids in three Yersinia pestis strains from the former Soviet Union (FSU) were determined and compared with those of pPCP1 plasmids in three well-characterized, non-FSU Y. pestis strains (KIM, CO92, and 91001). Two of the FSU plasmids were from strains C2614 and C2944, isolated from plague foci in Russia, and one plasmid was from strain C790 from Kyrgyzstan. Sequence analyses identified four sequence types among the six plasmids. The pPCP1 plasmids in the FSU strains were most genetically related to the pPCP1 plasmid in the KIM strain and least related to the pPCP1 plasmid in Y. pestis 91001. The FSU strains generally had larger pPCP1 plasmid copy numbers compared to strain CO92. Expression of the plasmid's pla gene was significantly (P≤.05) higher in strain C2944 than in strain CO92. Given pla's role in Y. pestis virulence, this difference may have important implications for the strain's virulence.