Invasion dynamics of white-nose syndrome fungus, midwestern United States, 2012-2014

White-nose syndrome has devastated bat populations in eastern North America. In Midwestern United States, prevalence increased quickly in the first year of invasion (2012-13) but with low population declines. In the second year (2013-14), environmental contamination led to earlier infection and high population declines. Interventions must be implemented before or soon after fungal invasion to prevent population collapse.

Investigation of Yersinia pestis Laboratory Adaptation through a Combined Genomics and Proteomics Approach

The bacterial pathogen Yersinia pestis, the cause of plague in humans and animals, normally has a sylvatic lifestyle, cycling between fleas and mammals. In contrast, laboratory-grown Y. pestis experiences a more constant environment and conditions that it would not normally encounter. The transition from the natural environment to the laboratory results in a vastly different set of selective pressures, and represents what could be considered domestication. Understanding the kinds of adaptations Y. pestis undergoes as it becomes domesticated will contribute to understanding the basic biology of this important pathogen.

In this study, we performed a parallel serial passage experiment (PSPE) to explore the mechanisms by which Y. pestis adapts to laboratory conditions, hypothesizing that cells would undergo significant changes in virulence and nutrient acquisition systems. Two wild strains were serially passaged in 12 independent populations each for ~750 generations, after which each population was analyzed using whole-genome sequencing, LC-MS/MS proteomic analysis, and GC/MS metabolomics. We observed considerable parallel evolution in the endpoint populations, detecting multiple independent mutations in ail, pepA, and zwf, suggesting that specific selective pressures are shaping evolutionary responses. Complementary LC-MS/MS proteomic data provide physiological context to the observed mutations, and reveal regulatory changes not necessarily associated with specific mutations, including changes in amino acid metabolism and cell envelope biogenesis. Proteomic data support hypotheses generated by genomic data in addition to suggesting future mechanistic studies, indicating that future whole-genome sequencing studies be designed to leverage proteomics as a critical complement.

Efficacy of Visual Surveys for White-Nose Syndrome at Bat Hibernacula

White-Nose Syndrome (WNS) is an epizootic disease in hibernating bats caused by the fungus Pseudogymnoascus destructans. Surveillance for P. destructans at bat hibernacula consists primarily of visual surveys of bats, collection of potentially infected bats, and submission of these bats for laboratory testing. Cryptic infections (bats that are infected but display no visual signs of fungus) could lead to the mischaracterization of the infection status of a site and the inadvertent spread of P. destructans. We determined the efficacy of visual detection of P. destructans by examining visual signs and molecular detection of P. destructans on 928 bats of six species at 27 sites during surveys conducted from January through March in 2012-2014 in the southeastern USA on the leading edge of the disease invasion. Cryptic infections were widespread with 77% of bats that tested positive by qPCR showing no visible signs of infection. The probability of exhibiting visual signs of infection increased with sampling date and pathogen load, the latter of which was substantially higher in three species (Myotis lucifugus, M. septentrionalis, and Perimyotis subflavus). In addition, M. lucifugus was more likely to show visual signs of infection than other species given the same pathogen load. Nearly all infections were cryptic in three species (Eptesicus fuscus, M. grisescens, and M. sodalis), which had much lower fungal loads. The presence of M. lucifugus or M. septentrionalis at a site increased the probability that P. destructans was visually detected on bats. Our results suggest that cryptic infections of P. destructans are common in all bat species, and visible infections rarely occur in some species. However, due to very high infection prevalence and loads in some species, we estimate that visual surveys examining at least 17 individuals of M. lucifugus and M. septentrionalis, or 29 individuals of P. subflavus are still effective to determine whether a site has bats infected with P. destructans. In addition, because the probability of visually detecting the fungus was higher later in winter, surveys should be done as close to the end of the hibernation period as possible.

Best practices for evaluating single nucleotide variant calling methods for microbial genomics

Innovations in sequencing technologies have allowed biologists to make incredible advances in understanding biological systems. As experience grows, researchers increasingly recognize that analyzing the wealth of data provided by these new sequencing platforms requires careful attention to detail for robust results. Thus far, much of the scientific Communit’s focus for use in bacterial genomics has been on evaluating genome assembly algorithms and rigorously validating assembly program performance. Missing, however, is a focus on critical evaluation of variant callers for these genomes. Variant calling is essential for comparative genomics as it yields insights into nucleotide-level organismal differences. Variant calling is a multistep process with a host of potential error sources that may lead to incorrect variant calls. Identifying and resolving these incorrect calls is critical for bacterial genomics to advance. The goal of this review is to provide guidance on validating algorithms and pipelines used in variant calling for bacterial genomics. First, we will provide an overview of the variant calling procedures and the potential sources of error associated with the methods. We will then identify appropriate datasets for use in evaluating algorithms and describe statistical methods for evaluating algorithm performance. As variant calling moves from basic research to the applied setting, standardized methods for performance evaluation and reporting are required; it is our hope that this review provides the groundwork for the development of these standards.

Phylogenetically typing bacterial strains from partial SNP genotypes observed from direct sequencing of clinical specimen metagenomic data

We describe an approach for genotyping bacterial strains from low coverage genome datasets, including metagenomic data from complex samples. Sequence reads from unknown samples are aligned to a reference genome where the allele states of known SNPs are determined. The Whole Genome Focused Array SNP Typing (WG-FAST) pipeline can identify unknown strains with much less read data than is needed for genome assembly. To test WG-FAST, we resampled SNPs from real samples to understand the relationship between low coverage metagenomic data and accurate phylogenetic placement. WG-FAST can be downloaded from https://github.com/jasonsahl/wgfast.

First confirmation of Pseudogymnoascus destructans in British bats and hibernacula

White-nose syndrome (WNS) is a fatal fungal infection of bats in North America caused by Pseudogymnoascus destructans. P. destructans has been confirmed in Continental Europe but not associated with mass mortality. Its presence in Great Britain was unknown. Opportunistic sampling of bats in GB began during the winter of 2009. Any dead bats or samples from live bats with visible fungal growths were submitted to the Animal Health and Veterinary Laboratories Agency for culture. Active surveillance by targeted environmental sampling of hibernacula was carried out during the winter of 2012/2013. Six hibernacula were selected by their proximity to Continental Europe. Five samples, a combination of surface swabs or sediment samples, were collected. These were sent to the Center for Microbial Genetics and Genomics, Northern Arizona University, for P. destructans PCR. Forty-eight incidents were investigated between March 2009 and July 2013. They consisted of 46 bat carcases and 31 other samples. A suspected P. destructans isolate was cultured from a live Daubenton's bat (Myotis daubentonii) sampled in February 2013. This isolate was confirmed by the Mycology Reference Laboratory, Bristol (Public Health England), as P. destructans. A variety of fungi were isolated from the rest but all were considered to be saprophytic or incidental. P. destructans was also confirmed by the Center for Microbial Genetics and Genomics in five of the six sites surveyed.

Host and pathogen ecology drive the seasonal dynamics of a fungal disease, white-nose syndrome

Seasonal patterns in pathogen transmission can influence the impact of disease on populations and the speed of spatial spread. Increases in host contact rates or births drive seasonal epidemics in some systems, but other factors may occasionally override these influences. White-nose syndrome, caused by the emerging fungal pathogen Pseudogymnoascus destructans, is spreading across North America and threatens several bat species with extinction. We examined patterns and drivers of seasonal transmission of P. destructans by measuring infection prevalence and pathogen loads in six bat species at 30 sites across the eastern United States. Bats became transiently infected in autumn, and transmission spiked in early winter when bats began hibernating. Nearly all bats in six species became infected by late winter when infection intensity peaked. In summer, despite high contact rates and a birth pulse, most bats cleared infections and prevalence dropped to zero. These data suggest the dominant driver of seasonal transmission dynamics was a change in host physiology, specifically hibernation. Our study is the first, to the best of our knowledge, to describe the seasonality of transmission in this emerging wildlife disease. The timing of infection and fungal growth resulted in maximal population impacts, but only moderate rates of spatial spread.

Development and assessment of multiplex high resolution melting assay as a tool for rapid single-tube identification of five Brucella species

BACKGROUND

The zoonosis brucellosis causes economically significant reproductive problems in livestock and potentially debilitating disease of humans. Although the causative agent, organisms from the genus Brucella, can be differentiated into a number of species based on phenotypic characteristics, there are also significant differences in genotype that are concordant with individual species. This paper describes the development of a five target multiplex assay to identify five terrestrial Brucella species using real-time polymerase chain reaction (PCR) and subsequent high resolution melt curve analysis. This technology offers a robust and cost effective alternative to previously described hydrolysis-probe Single Nucleotide Polymorphism (SNP)-based species defining assays.

RESULTS

Through the use of Brucella whole genome sequencing five species defining SNPs were identified. Individual HRM assays were developed to these target these changes and, following optimisation of primer concentrations, it was possible to multiplex all five assays in a single tube. In a validation exercise using a panel of 135 Brucella strains of terrestrial and marine origin, it was possible to distinguish the five target species from the other species within this panel.

CONCLUSION

The HRM multiplex offers a number of diagnostic advantages over previously described SNP-based typing approaches. Further, and uniquely for HRM, the successful multiplexing of five assays in a single tube allowing differentiation of five Brucella species in the diagnostic laboratory in a cost-effective and timely manner is described. However there are possible limitations to using this platform on DNA extractions direct from clinical material.

MetaGeniE: characterizing human clinical samples using deep metagenomic sequencing

With the decreasing cost of next-generation sequencing, deep sequencing of clinical samples provides unique opportunities to understand host-associated microbial communities. Among the primary challenges of clinical metagenomic sequencing is the rapid filtering of human reads to survey for pathogens with high specificity and sensitivity. Metagenomes are inherently variable due to different microbes in the samples and their relative abundance, the size and architecture of genomes, and factors such as target DNA amounts in tissue samples (i.e. human DNA versus pathogen DNA concentration). This variation in metagenomes typically manifests in sequencing datasets as low pathogen abundance, a high number of host reads, and the presence of close relatives and complex microbial communities. In addition to these challenges posed by the composition of metagenomes, high numbers of reads generated from high-throughput deep sequencing pose immense computational challenges. Accurate identification of pathogens is confounded by individual reads mapping to multiple different reference genomes due to gene similarity in different taxa present in the community or close relatives in the reference database. Available global and local sequence aligners also vary in sensitivity, specificity, and speed of detection. The efficiency of detection of pathogens in clinical samples is largely dependent on the desired taxonomic resolution of the organisms. We have developed an efficient strategy that identifies “all against all” relationships between sequencing reads and reference genomes. Our approach allows for scaling to large reference databases and then genome reconstruction by aggregating global and local alignments, thus allowing genetic characterization of pathogens at higher taxonomic resolution. These results were consistent with strain level SNP genotyping and bacterial identification from laboratory culture.

Changing climate and the altitudinal range of avian malaria in the Hawaiian Islands - an ongoing conservation crisis on the island of Kaua'i

Transmission of avian malaria in the Hawaiian Islands varies across altitudinal gradients and is greatest at elevations below 1500 m where both temperature and moisture are favorable for the sole mosquito vector, Culex quinquefasciatus, and extrinsic sporogonic development of the parasite, Plasmodium relictum. Potential consequences of global warming on this system have been recognized for over a decade with concerns that increases in mean temperatures could lead to expansion of malaria into habitats where cool temperatures currently limit transmission to highly susceptible endemic forest birds. Recent declines in two endangered species on the island of Kaua'i, the 'Akikiki (Oreomystis bairdi) and 'Akeke'e (Loxops caeruleirostris), and retreat of more common native honeycreepers to the last remaining high elevation habitat on the Alaka'i Plateau suggest that predicted changes in disease transmission may be occurring. We compared prevalence of malarial infections in forest birds that were sampled at three locations on the Plateau during 1994-1997 and again during 2007-2013, and also evaluated changes in the occurrence of mosquito larvae in available aquatic habitats during the same time periods. Prevalence of infection increased significantly at the lower (1100 m, 10.3% to 28.2%), middle (1250 m, 8.4% to 12.2%), and upper ends of the Plateau (1350 m, 2.0% to 19.3%). A concurrent increase in detections of Culex larvae in aquatic habitats associated with stream margins indicates that populations of the vector are also increasing. These increases are at least in part due to local transmission because overall prevalence in Kaua'i 'Elepaio (Chasiempis sclateri), a sedentary native species, has increased from 17.2% to 27.0%. Increasing mean air temperatures, declining precipitation, and changes in streamflow that have taken place over the past 20 years are creating environmental conditions throughout major portions of the Alaka'i Plateau that support increased transmission of avian malaria.

No selection on immunological markers in response to a highly virulent pathogen in an Arctic breeding bird

In natural populations, epidemics provide opportunities to look for intense natural selection on genes coding for life history and immune or other physiological traits. If the populations being considered are of management or conservation concern, then identifying the traits under selection (or 'markers') might provide insights into possible intervention strategies during epidemics. We assessed potential for selection on multiple immune and life history traits of Arctic breeding common eiders (Somateria mollissima) during annual avian cholera outbreaks (summers of 2006, 2007 & 2008). We measured prelaying body condition, immune traits, and subsequent reproductive investment (i.e., clutch size) and survival of female common eiders and whether they were infected with Pasteurella multocida, the causative agent of avian cholera. We found no clear and consistent evidence of directional selection on immune traits; however, infected birds had higher levels of haptoglobin than uninfected birds. Also, females that laid larger clutches had slightly lower immune responses during the prelaying period reflecting possible downregulation of the immune system to support higher costs of reproduction. This supports a recent study indicating that birds investing in larger clutches were more likely to die from avian cholera and points to a possible management option to maximize female survival during outbreaks.

Genetic relatedness of Brucella suis biovar 2 isolates from hares, wild boars and domestic pigs

Porcine brucellosis generally manifests as disorders in reproductive organs potentially leading to serious losses in the swine industry. Brucella suis biovar 2 is endemic in European wild boar (Sus scrofa) and hare (Lepus europeus, Lepus capensis) populations, thus these species may play a significant role in disease spread and serve as potential sources of infection for domestic pigs. The aim of this study was an epidemiologic analysis of porcine brucellosis in Hungary and a comparative analysis of B. suis bv. 2 strains from Europe using multiple-locus variable-number tandem repeat analysis (MLVA). MLVA-16 and its MLVA-11 subset were used to determine the genotypes of 68 B. suis bv. 2 isolates from Hungary and results were then compared to European MLVA genotypes. The analyses indicated relatively high genetic diversity of B. suis bv. 2 in Hungary. Strains isolated from hares and wild boars from Hungary showed substantial genetic divergence, suggesting separate lineages in each host and no instances of cross species infections. The closest relatives of strains from Hungarian wild boars and domestic pigs were mainly in the isolates from German and Croatian boars and pigs. The assessment of the European MLVA genotypes of wild boar isolates generally showed clustering based on geographic origin. The hare strains were relatively closely related to one another and did not cluster based on geographic origin. The limited relationships between geographic origin and genotype in isolates from hares might be the result of cross-border live animal translocation. The results could also suggest that certain B. suis strains are more adapted to hares. Across Europe, isolates from domestic pigs were closely related to isolates originating from both hares and wild boars, supporting the idea that wild animals are a source of brucellosis in domestic pigs.

Incongruence between multi-locus sequence analysis (MLSA) and whole-genome-based phylogenies: Pseudomonas syringae pathovar pisi as a cautionary tale

Previous phylogenies, built using a subset of genomic loci, split Pseudomonas syringae pv. pisi into two well-supported clades and implied convergence in host range for these lineages. The analysis of phenotypic and genotypic data within the context of this phylogenetic relationship implied further convergence at the level of virulence gene loss and acquisition. We generate draft genome assemblies for two additional P. syringae strains, isolated from diseased pea plants, and demonstrate incongruence between phylogenies created from a subset of the data compared with the whole genomes. Our whole-genome analysis demonstrates that strains classified as pv. pisi actually form a coherent monophyletic clade, so that apparent convergence is actually the product of shared ancestry. We use this example to urge caution when making evolutionary inferences across closely related strains of P. syringae.

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

Accurate and rapid identification of the Burkholderia pseudomallei near-neighbour, Burkholderia ubonensis, using real-time PCR

Burkholderia ubonensis is an environmental bacterium belonging to the Burkholderia cepacia complex (Bcc), a group of genetically related organisms that are associated with opportunistic but generally nonfatal infections in healthy individuals. In contrast, the near-neighbour species Burkholderia pseudomallei causes melioidosis, a disease that can be fatal in up to 95% of cases if left untreated. B. ubonensis is frequently misidentified as B. pseudomallei from soil samples using selective culturing on Ashdown’s medium, reflecting both the shared environmental niche and morphological similarities of these species. Additionally, B. ubonensis shows potential as an important biocontrol agent in B. pseudomallei-endemic regions as certain strains possess antagonistic properties towards B. pseudomallei. Current methods for characterising B. ubonensis are laborious, time-consuming and costly, and as such this bacterium remains poorly studied. The aim of our study was to develop a rapid and inexpensive real-time PCR-based assay specific for B. ubonensis. We demonstrate that a novel B. ubonensis-specific assay, Bu550, accurately differentiates B. ubonensis from B. pseudomallei and other species that grow on selective Ashdown’s agar. We anticipate that Bu550 will catalyse research on B. ubonensis by enabling rapid identification of this organism from Ashdown’s-positive colonies that are not B. pseudomallei.

Francisella guangzhouensis sp. nov., isolated from air-conditioning systems

Four strains (08HL01032(T), 09HG994, 10HP82-6 and 10HL1960) were isolated from water of air-conditioning systems of various cooling towers in Guangzhou city, China. Cells were Gram-stain-negative coccobacilli without flagella, catalase-positive and oxidase-negative, showing no reduction of nitrate, no hydrolysis of urea and no production of H2S. Growth was characteristically enhanced in the presence of l-cysteine, which was consistent with the properties of members of the genus Francisella. The quinone system was composed of ubiquinone Q-8 with minor amounts of Q-9. The polar lipid profile consisted of the predominant lipids phosphatidylethanolamine, diphosphatidylglycerol, phosphatidylglycerol, phosphatidylcholine, two unidentified phospholipids (PL2, PL3), an unidentified aminophospholipid and an unidentified glycolipid (GL2). The polyamine pattern consisted of the major compounds spermidine, cadaverine and spermine. The major cellular fatty acids were C10 : 0, C14 : 0, C16 : 0, C18 : 1ω9c and C18 : 1 3-OH. A draft whole-genome sequence of the proposed type strain 08HL01032(T) was generated. Comparative sequence analysis of the complete 16S and 23S rRNA genes confirmed affiliation to the genus Francisella, with 95 % sequence identity to the closest relatives in the database, the type strains of Francisella philomiragia and Francisella noatunensis subsp. orientalis. Full-length deduced amino acid sequences of various housekeeping genes, recA, gyrB, groEL, dnaK, rpoA, rpoB, rpoD, rpoH, fopA and sdhA, exhibited similarities of 67-92 % to strains of other species of the genus Francisella. Strains 08HL01032(T), 09HG994, 10HP82-6 and 10HL1960 exhibited highly similar pan-genome PCR profiles. Both the phenotypic and molecular data support the conclusion that the four strains belong to the genus Francisella but exhibit considerable divergence from all recognized Francisella species. Therefore, we propose the name Francisella guangzhouensis sp. nov., with the type strain 08HL01032(T) ( = CCUG 60119(T) = NCTC 13503(T)).

Within-host evolution of Brucella canis during a canine brucellosis outbreak in a kennel

Background

Little is currently known about Brucella evolution within the host during infection. The current study is the first to employ fine-scale genotyping on an isolate collection derived from a Brucella canis outbreak. Eight isolates of B. canis, cultured from different tissues of three dogs (female, stud dog, puppy of another female) from a single kennel over three months were genetically characterized with a 15-marker multi-locus, variable-number tandem repeat (VNTR) analysis (MLVA) to assess the genetic relatedness of isolates and potential rapid mutational changes.

Results

MLVA discriminated among the otherwise indistinguishable isolates from different animals and from isolates collected at different time points within each host, with different VNTR alleles being detected at multiple dates and tissue sites. We suspect that all isolates cultured from the female, puppy, and stud dogs originated from the same strain, with subsequent rapid in vivo mutations. However, high mutation rates and apparent in several of the loci prevented making definitive epidemiological relationships among isolates.

Conclusions

This investigation highlights the rapid in vivo genetic mutations of several VNTRs of B. canis over a short time period in the host and the emergence of alternate alleles. However, this work also suggests the challenges of using highly mutable VNTRs to infer epidemiological relationships of strains within a short duration outbreak.

Molecular epidemiologic investigation of an anthrax outbreak among heroin users, Europe

In December 2009, two unusual cases of anthrax were diagnosed in heroin users in Scotland. A subsequent anthrax outbreak in heroin users emerged throughout Scotland and expanded into England and Germany, sparking concern of nefarious introduction of anthrax spores into the heroin supply. To better understand the outbreak origin, we used established genetic signatures that provided insights about strain origin. Next, we sequenced the whole genome of a representative Bacillus anthracis strain from a heroin user (Ba4599), developed Ba4599-specific single-nucleotide polymorphism assays, and genotyped all available material from other heroin users with anthrax. Of 34 case-patients with B. anthracis-positive PCR results, all shared the Ba4599 single-nucleotide polymorphism genotype. Phylogeographic analysis demonstrated that Ba4599 was closely related to strains from Turkey and not to previously identified isolates from Scotland or Afghanistan, the presumed origin of the heroin. Our results suggest accidental contamination along the drug trafficking route through a cutting agent or animal hides used to smuggle heroin into Europe.

Genotyping of Brucella species using clade specific SNPs

Background

Brucellosis is a worldwide disease of mammals caused by Alphaproteobacteria in the genus Brucella. The genus is genetically monomorphic, requiring extensive genotyping to differentiate isolates. We utilized two different genotyping strategies to characterize isolates. First, we developed a microarray-based assay based on 1000 single nucleotide polymorphisms (SNPs) that were identified from whole genome comparisons of two B. abortus isolates , one B. melitensis, and one B. suis. We then genotyped a diverse collection of 85 Brucella strains at these SNP loci and generated a phylogenetic tree of relationships. Second, we developed a selective primer-extension assay system using capillary electrophoresis that targeted 17 high value SNPs across 8 major branches of the phylogeny and determined their genotypes in a large collection ( n = 340) of diverse isolates.

Results

Our 1000 SNP microarray readily distinguished B. abortus, B. melitensis, and B. suis, differentiating B. melitensis and B. suis into two clades each. Brucella abortus was divided into four major clades. Our capillary-based SNP genotyping confirmed all major branches from the microarray assay and assigned all samples to defined lineages. Isolates from these lineages and closely related isolates, among the most commonly encountered lineages worldwide, can now be quickly and easily identified and genetically characterized.

Conclusions

We have identified clade-specific SNPs in Brucella that can be used for rapid assignment into major groups below the species level in the three main Brucella species. Our assays represent SNP genotyping approaches that can reliably determine the evolutionary relationships of bacterial isolates without the need for whole genome sequencing of all isolates.

Development and validation of Burkholderia pseudomallei-specific real-time PCR assays for clinical, environmental or forensic detection applications

The bacterium Burkholderia pseudomallei causes melioidosis, a rare but serious illness that can be fatal if untreated or misdiagnosed. Species-specific PCR assays provide a technically simple method for differentiating B. pseudomallei from near-neighbor species. However, substantial genetic diversity and high levels of recombination within this species reduce the likelihood that molecular signatures will differentiate all B. pseudomallei from other Burkholderiaceae. Currently available molecular assays for B. pseudomallei detection lack rigorous validation across large in silico datasets and isolate collections to test for specificity, and none have been subjected to stringent quality control criteria (accuracy, precision, selectivity, limit of quantitation (LoQ), limit of detection (LoD), linearity, ruggedness and robustness) to determine their suitability for environmental, clinical or forensic investigations. In this study, we developed two novel B. pseudomallei specific assays, 122018 and 266152, using a dual-probe approach to differentiate B. pseudomallei from B. thailandensis, B. oklahomensis and B. thailandensis-like species; other species failed to amplify. Species specificity was validated across a large DNA panel (>2,300 samples) comprising Burkholderia spp. and non-Burkholderia bacterial and fungal species of clinical and environmental relevance. Comparison of assay specificity to two previously published B. pseudomallei-specific assays, BurkDiff and TTS1, demonstrated comparable performance of all assays, providing between 99.7 and 100% specificity against our isolate panel. Last, we subjected 122018 and 266152 to rigorous quality control analyses, thus providing quantitative limits of assay performance. Using B. pseudomallei as a model, our study provides a framework for comprehensive quantitative validation of molecular assays and provides additional, highly validated B. pseudomallei assays for the scientific research community.

Epidemiological tracking and population assignment of the non-clonal bacterium, Burkholderia pseudomallei

Rapid assignment of bacterial pathogens into predefined populations is an important first step for epidemiological tracking. For clonal species, a single allele can theoretically define a population. For non-clonal species such as Burkholderia pseudomallei, however, shared allelic states between distantly related isolates make it more difficult to identify population defining characteristics. Two distinct B. pseudomallei populations have been previously identified using multilocus sequence typing (MLST). These populations correlate with the major foci of endemicity (Australia and Southeast Asia). Here, we use multiple Bayesian approaches to evaluate the compositional robustness of these populations, and provide assignment results for MLST sequence types (STs). Our goal was to provide a reference for assigning STs to an established population without the need for further computational analyses. We also provide allele frequency results for each population to enable estimation of population assignment even when novel STs are discovered. The ability for humans and potentially contaminated goods to move rapidly across the globe complicates the task of identifying the source of an infection or outbreak. Population genetic dynamics of B. pseudomallei are particularly complicated relative to other bacterial pathogens, but the work here provides the ability for broad scale population assignment. As there is currently no independent empirical measure of successful population assignment, we provide comprehensive analytical details of our comparisons to enable the reader to evaluate the robustness of population designations and assignments as they pertain to individual research questions. Finer scale subdivision and verification of current population compositions will likely be possible with genotyping data that more comprehensively samples the genome. The approach used here may be valuable for other non-clonal pathogens that lack simple group-defining genetic characteristics and provides a rapid reference for epidemiologists wishing to track the origin of infection without the need to compile population data and learn population assignment algorithms.

Burkholderia pseudomallei is a soil-dwelling bacterium that can infect a large range of hosts. In humans, B. pseudomallei causes melioidosis, and typical routes of entry include open wounds, inhalation, or ingestion. Clinical features are diverse, although pneumonia and abscess formation are common. High rates of recombination within the genome of this bacterium have confounded attempts to match clinical samples to geographically defined populations. Here we provide a reference that simplifies source attribution issues. We applied population assignment software to previously generated sequence data from seven B. pseudomallei genes to define the major geographic populations within this species. We evaluated the robustness of our results by comparison with two additional population assignment programs. We present the likelihood that each variant is assigned to a particular geographic population. This information can be used to assign novel B. pseudomallei isolates to a geographic population without needing to learn and run cumbersome population assignment applications. This method can also be used for other bacteria that are difficult to source-attribute due to high levels of genomic variation and recombination.

Castor bean organelle genome sequencing and worldwide genetic diversity analysis

Castor bean is an important oil-producing plant in the Euphorbiaceae family. Its high-quality oil contains up to 90% of the unusual fatty acid ricinoleate, which has many industrial and medical applications. Castor bean seeds also contain ricin, a highly toxic Type 2 ribosome-inactivating protein, which has gained relevance in recent years due to biosafety concerns. In order to gain knowledge on global genetic diversity in castor bean and to ultimately help the development of breeding and forensic tools, we carried out an extensive chloroplast sequence diversity analysis. Taking advantage of the recently published genome sequence of castor bean, we assembled the chloroplast and mitochondrion genomes extracting selected reads from the available whole genome shotgun reads. Using the chloroplast reference genome we used the methylation filtration technique to readily obtain draft genome sequences of 7 geographically and genetically diverse castor bean accessions. These sequence data were used to identify single nucleotide polymorphism markers and phylogenetic analysis resulted in the identification of two major clades that were not apparent in previous population genetic studies using genetic markers derived from nuclear DNA. Two distinct sub-clades could be defined within each major clade and large-scale genotyping of castor bean populations worldwide confirmed previously observed low levels of genetic diversity and showed a broad geographic distribution of each sub-clade.

Phylogeography of Francisella tularensis subspecies holarctica from the country of Georgia

Background

Francisella tularensis, the causative agent of tularemia, displays subspecies-specific differences in virulence, geographic distribution, and genetic diversity. F. tularensis subsp. holarctica is widely distributed throughout the Northern Hemisphere. In Europe, F. tularensis subsp. holarctica isolates have largely been assigned to two phylogenetic groups that have specific geographic distributions. Most isolates from Western Europe are assigned to the B.Br.FTNF002-00 group, whereas most isolates from Eastern Europe are assigned to numerous lineages within the B.Br.013 group. The eastern geographic extent of the B.Br.013 group is currently unknown due to a lack of phylogenetic knowledge about populations at the European/Asian juncture and in Asia. In this study, we address this knowledge gap by describing the phylogenetic structure of F. tularensis subsp. holarctica isolates from the country of Georgia, and by placing these isolates into a global phylogeographic context.

Results

We identified a new genetic lineage of F. tularensis subsp. holarctica from Georgia that belongs to the B.Br.013 group. This new lineage is genetically and geographically distinct from lineages previously described from the B.Br.013 group from Central-Eastern Europe. Importantly, this new lineage is basal within the B.Br.013 group, indicating the Georgian lineage diverged before the diversification of the other known B.Br.013 lineages. Although two isolates from the Georgian lineage were collected nearby in the Ukrainian region of Crimea, all other global isolates assigned to this lineage were collected in Georgia. This restricted geographic distribution, as well as the high levels of genetic diversity within the lineage, is consistent with a relatively older origin and localized differentiation.

Conclusions

We identified a new lineage of F. tularensis subsp. holarctica from Georgia that appears to have an older origin than any other diversified lineages previously described from the B.Br.013 group. This finding suggests that additional phylogenetic studies of F. tularensis subsp. holarctica populations in Eastern Europe and Asia have the potential to yield important new insights into the evolutionary history and phylogeography of this broadly dispersed F. tularensis subspecies.