Evan W. Skowronski

Iridovirus and Microsporidian Linked to Honey Bee Colony Decline

Background In 2010 Colony Collapse Disorder (CCD), again devastated honey bee colonies in the USA, indicating that the problem is neither diminishing nor has it been resolved. Many CCD investigations, using sensitive genome-based methods, have found small RNA bee viruses and the microsporidia, Nosema apis and N. ceranae in healthy and collapsing colonies alike with no single pathogen firmly linked to honey bee losses. Methodology/Principal Findings We used Mass spectrometry-based proteomics (MSP) to identify and quantify thousands of proteins from healthy and collapsing bee colonies. MSP revealed two unreported RNA viruses in North American honey bees, Varroa destructor-1 virus and Kakugo virus, and identified an invertebrate iridescent virus (IIV) (Iridoviridae) associated with CCD colonies. Prevalence of IIV significantly discriminated among strong, failing, and collapsed colonies. In addition, bees in failing colonies contained not only IIV, but also Nosema. Co-occurrence of these microbes consistently marked CCD in (1) bees from commercial apiaries sampled across the U.S. in 2006–2007, (2) bees sequentially sampled as the disorder progressed in an observation hive colony in 2008, and (3) bees from a recurrence of CCD in Florida in 2009. The pathogen pairing was not observed in samples from colonies with no history of CCD, namely bees from Australia and a large, non-migratory beekeeping business in Montana. Laboratory cage trials with a strain of IIV type 6 and Nosema ceranae confirmed that co-infection with these two pathogens was more lethal to bees than either pathogen alone. Conclusions/Significance These findings implicate co-infection by IIV and Nosema with honey bee colony decline, giving credence to older research pointing to IIV, interacting with Nosema and mites, as probable cause of bee losses in the USA, Europe, and Asia. We next need to characterize the IIV and Nosema that we detected and develop management practices to reduce honey bee losses.

Surface sampling of spores in dry-deposition aerosols.

ABSTRACT The ability to reliably and reproducibly sample surfaces contaminated with a biological agent is a critical step in measuring the extent of contamination and determining if decontamination steps have been successful. The recovery operations following the 2001 attacks with Bacillus anthracis spores were complicated by the fact that no standard sample collection format or decontamination procedures were established. Recovery efficiencies traditionally have been calculated based upon biological agents which were applied to test surfaces in a liquid format and then allowed to dry prior to sampling tests, which may not be best suited for a real-world event with aerosolized biological agents. In order to ascertain if differences existed between air-dried liquid deposition and biological spores which were allowed to settle on a surface in a dried format, a study was undertaken to determine if differences existed in surface sampling recovery efficiencies for four representative surfaces. Studies were then undertaken to compare sampling efficiencies between liquid spore deposition and aerosolized spores which were allowed to gradually settle under gravity on four different test coupon types. Tests with both types of deposition compared efficiencies of four unique swabbing materials applied to four surfaces with various surface properties. Our studies demonstrate that recovery of liquid-deposited spores differs significantly from recovery of dry aerosol-deposited spores in most instances. Whether the recovery of liquid-deposited spores is overexaggerated or underrepresented with respect to that of aerosol-deposited spores depends upon the surface material being tested.

Double-Blind Characterization of Non-Genome-Sequenced Bacteria by Mass Spectrometry-Based Proteomics

ABSTRACT Due to the possibility of a biothreat attack on civilian or military installations, a need exists for technologies that can detect and accurately identify pathogens in a near-real-time approach. One technology potentially capable of meeting these needs is a high-throughput mass spectrometry (MS)-based proteomic approach. This approach utilizes the knowledge of amino acid sequences of peptides derived from the proteolysis of proteins as a basis for reliable bacterial identification. To evaluate this approach, the tryptic digest peptides generated from double-blind biological samples containing either a single bacterium or a mixture of bacteria were analyzed using liquid chromatography-tandem mass spectrometry. Bioinformatic tools that provide bacterial classification were used to evaluate the proteomic approach. Results showed that bacteria in all of the double-blind samples were accurately identified with no false-positive assignment. The MS proteomic approach showed strain-level discrimination for the various bacteria employed. The approach also characterized double-blind bacterial samples to the respective genus, species, and strain levels when the experimental organism was not in the database due to its genome not having been sequenced. One experimental sample did not have its genome sequenced, and the peptide experimental record was added to the virtual bacterial proteome database. A replicate analysis identified the sample to the peptide experimental record stored in the database. The MS proteomic approach proved capable of identifying and classifying organisms within a microbial mixture.

Genomic Signatures of Strain Selection and Enhancement in Bacillus atrophaeus var. globigii ,a Historical Biowarfare Simulant

Background Despite 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). Results Archival 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. Conclusions Our 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.

Genetic Barcodes for Improved Environmental Tracking of an Anthrax Simulant

ABSTRACT The development of realistic risk models that predict the dissemination, dispersion and persistence of potential biothreat agents have utilized nonpathogenic surrogate organisms such as Bacillus atrophaeus subsp. globigii or commercial products such as Bacillus thuringiensis subsp. kurstaki. Comparison of results from outdoor tests under different conditions requires the use of genetically identical strains; however, the requirement for isogenic strains limits the ability to compare other desirable properties, such as the behavior in the environment of the same strain prepared using different methods. Finally, current methods do not allow long-term studies of persistence or reaerosolization in test sites where simulants are heavily used or in areas where B. thuringiensis subsp. kurstaki is applied as a biopesticide. To create a set of genetically heterogeneous yet phenotypically indistinguishable strains so that variables intrinsic to simulations (e.g., sample preparation) can be varied and the strains can be tested under otherwise identical conditions, we have developed a strategy of introducing small genetic signatures (“barcodes”) into neutral regions of the genome. The barcodes are stable over 300 generations and do not impact in vitro growth or sporulation. Each barcode contains common and specific tags that allow differentiation of marked strains from wild-type strains and from each other. Each tag is paired with specific real-time PCR assays that facilitate discrimination of barcoded strains from wild-type strains and from each other. These uniquely barcoded strains will be valuable tools for research into the environmental fate of released organisms by providing specific artificial detection signatures.

Comparative genomics of 2009 seasonal plague (Yersinia pestis) in New Mexico.

Plague disease caused by the Gram-negative bacterium Yersinia pestis routinely affects animals and occasionally humans, in the western United States. The strains native to the North American continent are thought to be derived from a single introduction in the late 19th century. The degree to which these isolates have diverged genetically since their introduction is not clear, and new genomic markers to assay the diversity of North American plague are highly desired. To assay genetic diversity of plague isolates within confined geographic areas, draft genome sequences were generated by 454 pyrosequencing from nine environmental and clinical plague isolates. In silico assemblies of Variable Number Tandem Repeat (VNTR) loci were compared to laboratory-generated profiles for seven markers. High-confidence SNPs and small Insertion/Deletions (Indels) were compared to previously sequenced Y. pestis isolates. The resulting panel of mutations allowed clustering of the strains and tracing of the most likely evolutionary trajectory of the plague strains. The sequences also allowed the identification of new putative SNPs that differentiate the 2009 isolates from previously sequenced plague strains and from each other. In addition, new insertion points for the abundant insertion sequences (IS) of Y. pestis are present that allow additional discrimination of strains; several of these new insertions potentially inactivate genes implicated in virulence. These sequences enable whole-genome phylogenetic analysis and allow the unbiased comparison of closely related isolates of a genetically monomorphic pathogen.

Pathoadaptive Mutations in Salmonella enterica Isolated after Serial Passage in Mice

How pathogenic bacteria adapt and evolve in the complex and variable environment of the host remains a largely unresolved question. Here we have used whole genome sequencing of Salmonella enterica serovar Typhimurium LT2 populations serially passaged in mice to identify mutations that adapt bacteria to systemic growth in mice. We found unique pathoadaptive mutations in two global regulators, phoQ and stpA, which increase the competitive indexes of the bacteria 3- to 5-fold. Also, all mouse-adapted lineages had changed the orientation of the hin invertable element, resulting in production of a FliC type of flagellum. Competition experiments in mice with locked flagellum mutants showed that strains expressing the FliC type of flagellum had a 5-fold increase in competitive index as compared to those expressing FljB type flagellum. Combination of the flagellum cassette inversion with the stpA mutation increased competitive indexes up to 20-fold. These experiments show that Salmonella can rapidly adapt to a mouse environment by acquiring a few mutations of moderate individual effect that when combined confer substantial increases in growth.

Population-level variation of the preproricin gene contradicts expectation of neutral equilibrium for generalist plant defense toxins.

The preproricin gene encodes ricin, the highly toxic, type II ribosome-inactivating protein of castor bean (Ricinus communis L.). As a generalist plant defense gene, preproricin is expected to exhibit population-level variation consistent with the neutral equilibrium model and to comprise few functionally different alleles. We first test the hypothesis that the preproricin gene family should comprise six to eight members by searching the publicly available draft genome sequence of R. communis and analyzing its ricin-like loci. We then test the neutral equilibrium expectation for the preproricin gene by characterizing its allelic variation among 25 geographically diverse castor bean plants. We confirm the presence of six ricin-like loci that share with the preproricin gene 62.9-96.3% nucleotide identity and intact A-chains. DNA sequence variation among the preproricin haplotypes significantly rejects tests of the neutral equilibrium model. Replacement mutations preserve the 12 amino acids known to affect catalytic and electrostatic interactions of the native protein toxin, which suggests functional divergence among alleles has been minimal. Nucleotide polymorphism is maintained by purifying selection (omega < 0.3) yet includes an excess of rare silent mutations greater than predicted by the neutral equilibrium model. Development of robust detection methods for ricin contamination must account for the presence of these other ricin-like molecules and should leverage the specificity provided by the numerous single nucleotide polymorphisms in the preproricin gene.

Shared strategies for β-lactam catabolism in the soil microbiome

The soil microbiome can produce, resist, or degrade antibiotics and even catabolize them. While resistance genes are widely distributed in the soil, there is a dearth of knowledge concerning antibiotic catabolism. Here we describe a pathway for penicillin catabolism in four isolates. Genomic and transcriptomic sequencing revealed β-lactamase, amidase, and phenylacetic acid catabolon upregulation. Knocking out part of the phenylacetic acid catabolon or an apparent penicillin utilization operon (put) resulted in loss of penicillin catabolism in one isolate. A hydrolase from the put operon was found to degrade in vitro benzylpenicilloic acid, the β-lactamase penicillin product. To test the generality of this strategy, an Escherichia coli strain was engineered to co-express a β-lactamase and a penicillin amidase or the put operon, enabling it to grow using penicillin or benzylpenicilloic acid, respectively. Elucidation of additional pathways may allow bioremediation of antibiotic-contaminated soils and discovery of antibiotic-remodeling enzymes with industrial utility.A β-lactamase, a novel type of amidase, and the phenylacetic acid catabolon comprise a catabolic pathway, revealed by genomic and transcriptomic analysis, that enables multiple soil bacteria to use β-lactam antibiotics as a carbon source.

Hybrid Vibrio cholerae El Tor Lacking SXT Identified as the Cause of a Cholera Outbreak in the Philippines

ABSTRACT Cholera continues to be a global threat, with high rates of morbidity and mortality. In 2011, a cholera outbreak occurred in Palawan, Philippines, affecting more than 500 people, and 20 individuals died. Vibrio cholerae O1 was confirmed as the etiological agent. Source attribution is critical in cholera outbreaks for proper management of the disease, as well as to control spread. In this study, three V. cholerae O1 isolates from a Philippines cholera outbreak were sequenced and their genomes analyzed to determine phylogenetic relatedness to V. cholerae O1 isolates from recent outbreaks of cholera elsewhere. The Philippines V. cholerae O1 isolates were determined to be V. cholerae O1 hybrid El Tor belonging to the seventh-pandemic clade. They clustered tightly, forming a monophyletic clade closely related to V. cholerae O1 hybrid El Tor from Asia and Africa. The isolates possess a unique multilocus variable-number tandem repeat analysis (MLVA) genotype (12-7-9-18-25 and 12-7-10-14-21) and lack SXT. In addition, they possess a novel 15-kb genomic island (GI-119) containing a predicted type I restriction-modification system. The CTXΦ-RS1 array of the Philippines isolates was similar to that of V. cholerae O1 MG116926, a hybrid El Tor strain isolated in Bangladesh in 1991. Overall, the data indicate that the Philippines V. cholerae O1 isolates are unique, differing from recent V. cholerae O1 isolates from Asia, Africa, and Haiti. Furthermore, the results of this study support the hypothesis that the Philippines isolates of V. cholerae O1 are indigenous and exist locally in the aquatic ecosystem of the Philippines. IMPORTANCE Genetic characterization and phylogenomics analysis of outbreak strains have proven to be critical for probing clonal relatedness to strains isolated in different geographical regions and over time. Recently, extensive genetic analyses of V. cholerae O1 strains isolated in different countries have been done. However, genome sequences of V. cholerae O1 isolates from the Philippines have not been available for epidemiological investigation. In this study, molecular typing and phylogenetic analysis of Vibrio cholerae isolated from both clinical and environmental samples in 2011 confirmed unique genetic features of the Philippines isolates, which are helpful to understand the global epidemiology of cholera. Genetic characterization and phylogenomics analysis of outbreak strains have proven to be critical for probing clonal relatedness to strains isolated in different geographical regions and over time. Recently, extensive genetic analyses of V. cholerae O1 strains isolated in different countries have been done. However, genome sequences of V. cholerae O1 isolates from the Philippines have not been available for epidemiological investigation. In this study, molecular typing and phylogenetic analysis of Vibrio cholerae isolated from both clinical and environmental samples in 2011 confirmed unique genetic features of the Philippines isolates, which are helpful to understand the global epidemiology of cholera.