Whole genome single nucleotide polymorphism based phylogeny of Francisella tularensis and its application to the development of a strain typing assay

Evidente-a visual analytics tool for data enrichment in SNP-based phylogenetic trees

Motivation

A common practice in the analysis of pathogens and their strains is using single-nucleotide polymorphisms (SNPs) to reconstruct their evolutionary history. However, genome-wide SNP-based phylogenetic trees are rarely analyzed without any further information. Including the underlying SNP data together with further metadata on the respective samples in the exploration process can facilitate linking the genomic and phenotypic properties of the samples.

Results

We introduce Efficient VIsual analytics tool for Data ENrichment in phylogenetic TreEs (Evidente), a web-application that provides an interactive visual analysis interface for the simultaneous interrogation of phylogenetic relationships, genome-wide SNP data and metadata for samples of an organism. Besides visualizing the phylogenetic tree, Evidente classifies SNPs as supporting or non-supporting of the tree structures and shows the distribution of both types of SNPs among samples and clades of interest. Furthermore, additional metadata can be included in the visualization. Lastly, Evidente includes an enrichment analysis to identify over-represented genomic features encoded by GO-terms within the clades of the tree. We demonstrate the usability of Evidente with the data of the pathogens Treponema pallidum and Mycobacterium leprae.

Availability and implementation

Evidente is available at the TueVis visualization web server at https://evidente-tuevis.cs.uni-tuebingen.de/, it can also be run locally.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Development, Phenotypic Characterization and Genomic Analysis of a Francisella tularensis Panel for Tularemia Vaccine Testing

Francisella tularensis is one of several biothreat agents for which a licensed vaccine is needed to protect against this pathogen. To aid in the development of a vaccine protective against pneumonic tularemia, we generated and characterized a panel of F. tularensis isolates that can be used as challenge strains to assess vaccine efficacy. Our panel consists of both historical and contemporary isolates derived from clinical and environmental sources, including human, tick, and rabbit isolates. Whole genome sequencing was performed to assess the genetic diversity in comparison to the reference genome F. tularensis Schu S4. Average nucleotide identity analysis showed >99% genomic similarity across the strains in our panel, and pan-genome analysis revealed a core genome of 1,707 genes, and an accessory genome of 233 genes. Three of the strains in our panel, FRAN254 (tick-derived), FRAN255 (a type B strain), and FRAN256 (a human isolate) exhibited variation from the other strains. Moreover, we identified several unique mutations within the Francisella Pathogenicity Island across multiple strains in our panel, revealing unexpected diversity in this region. Notably, FRAN031 (Scherm) completely lacked the second pathogenicity island but retained virulence in mice. In contrast, FRAN037 (Coll) was attenuated in a murine pneumonic tularemia model and had mutations in pdpB and iglA which likely led to attenuation. All of the strains, except FRAN037, retained full virulence, indicating their effectiveness as challenge strains for future vaccine testing. Overall, we provide a well-characterized panel of virulent F. tularensis strains that can be utilized in ongoing efforts to develop an effective vaccine against pneumonic tularemia to ensure protection is achieved across a range F. tularensis strains.

The Francisella tularensis Polysaccharides: What Is the Real Capsule?

Francisella tularensis is a tier 1 select agent responsible for tularemia in humans and a wide variety of animal species. Extensive research into understanding the virulence factors of the bacterium has been ongoing to develop an efficacious vaccine. At least two such virulence factors are described as capsules of F. tularensis: the O-antigen capsule and the capsule-like complex (CLC). These two separate entities aid in avoiding host immune defenses but have not been clearly differentiated. These components are distinct and differ in composition and genetic basis. The O-antigen capsule consists of a polysaccharide nearly identical to the lipopolysaccharide (LPS) O antigen, whereas the CLC is a heterogeneous complex of glycoproteins, proteins, and possibly outer membrane vesicles and tubes (OMV/Ts). In this review, the current understanding of these two capsules is summarized, and the historical references to "capsules" of F. tularensis are clarified. A significant amount of research has been invested into the composition of each capsule and the genes involved in synthesis of the polysaccharide portion of each capsule. Areas of future research include further exploration into the molecular regulation and pathways responsible for expression of each capsule and further elucidating the role that each capsule plays in virulence.

Russian isolates enlarge the known geographic diversity of Francisella tularensis subsp. mediasiatica

Francisella tularensis, a small Gram-negative bacterium, is capable of infecting a wide range of animals, including humans, and causes a plague-like disease called tularemia—a highly contagious disease with a high mortality rate. Because of these characteristics, F. tularensis is considered a potential agent of biological terrorism. Currently, F. tularensis is divided into four subspecies, which differ in their virulence and geographic distribution. Two of them, subsp. tularensis (primarily found in North America) and subsp. holarctica (widespread across the Northern Hemisphere), are responsible for tularemia in humans. Subsp. novicida is almost avirulent in humans. The fourth subspecies, subsp. mediasiatica, is the least studied because of its limited distribution and impact in human health. It is found only in sparsely populated regions of Central Asia. In this report, we describe the first focus of naturally circulating F. tularensis subsp. mediasiatica in Russia. We isolated and characterized 18 strains of this subspecies in the Altai region. All strains were highly virulent in mice. The virulence of subsp. mediasiatica in a vaccinated mouse model is intermediate between that of subsp. tularensis and subsp. holarctica. Based on a multiple-locus variable number tandem repeat analysis (MLVA), we show that the Altaic population of F. tularensis subsp. mediasiatica is genetically distinct from the classical Central Asian population, and probably is endemic to Southern Siberia. We propose to subdivide the mediasiatica subspecies into three phylogeographic groups, M.I, M.II and M.III.

Whole-Genome Relationships among Francisella Bacteria of Diverse Origins Define New Species and Provide Specific Regions for Detection

Francisella tularensis is a highly virulent zoonotic pathogen that causes tularemia and, because of weaponization efforts in past world wars, is considered a tier 1 biothreat agent. Detection and surveillance of F. tularensis may be confounded by the presence of uncharacterized, closely related organisms. Through DNA-based diagnostics and environmental surveys, novel clinical and environmental Francisella isolates have been obtained in recent years. Here we present 7 new Francisella genomes and a comparison of their characteristics to each other and to 24 publicly available genomes as well as a comparative analysis of 16S rRNA and sdhA genes from over 90 Francisella strains. Delineation of new species in bacteria is challenging, especially when isolates having very close genomic characteristics exhibit different physiological features-for example, when some are virulent pathogens in humans and animals while others are nonpathogenic or are opportunistic pathogens. Species resolution within Francisella varies with analyses of single genes, multiple gene or protein sets, or whole-genome comparisons of nucleic acid and amino acid sequences. Analyses focusing on single genes (16S rRNA, sdhA), multiple gene sets (virulence genes, lipopolysaccharide [LPS] biosynthesis genes, pathogenicity island), and whole-genome comparisons (nucleotide and protein) gave congruent results, but with different levels of discrimination confidence. We designate four new species within the genus; Francisella opportunistica sp. nov. (MA06-7296), Francisella salina sp. nov. (TX07-7308), Francisella uliginis sp. nov. (TX07-7310), and Francisella frigiditurris sp. nov. (CA97-1460). This study provides a robust comparative framework to discern species and virulence features of newly detected Francisella bacteria.

IMPORTANCE

DNA-based detection and sequencing methods have identified thousands of new bacteria in the human body and the environment. In most cases, there are no cultured isolates that correspond to these sequences. While DNA-based approaches are highly sensitive, accurately assigning species is difficult without known near relatives for comparison. This ambiguity poses challenges for clinical cases, disease epidemics, and environmental surveillance, for which response times must be short. Many new Francisella isolates have been identified globally. However, their species designations and potential for causing human disease remain ambiguous. Through detailed genome comparisons, we identified features that differentiate F. tularensis from clinical and environmental Francisella isolates and provide a knowledge base for future comparison of Francisella organisms identified in clinical samples or environmental surveys.

Rapid Identification and Characterization of Francisella by Molecular Biology and Other Techniques

Francisella tularensis is the causative pathogen of tularemia and a Tier 1 bioterror agent on the CDC list. Considering the fact that some subpopulation of the F. tularensis strains is more virulent, more significantly associated with mortality, and therefore poses more threat to humans, rapid identification and characterization of this subpopulation strains is of invaluable importance. This review summarizes the up-to-date developments of assays for mainly detecting and characterizing F. tularensis and a touch of caveats of some of the assays.

Multiple-locus variable-number tandem-repeat analysis of Francisella tularensis from Quebec, Canada

Francisella tularensis is ubiquitous in the Northern Hemisphere. Yet, little is known about the disease and its ecology within Canada as few serological studies have shown exposure to the disease and fewer case studies have been reported. This report is the first to describe the molecular subtyping of F. tularensis isolates within eastern Canada using multiple-locus variable-number tandem-repeat analysis. From 1998 to 2011, a total of 73 specimens were isolated from unique human and animal sources. As expected, F. tularensis subsp. tularensis AI and F. tularensis subsp. holarctica subtypes were observed, corresponding to the known geographical division within this species. The majority of human isolates (78%) and all animal (hare) isolates were of the more virulent, AI type. Half of the B isolates were isolated from patients living in a region of Quebec where muskrat densities are known to be high. A relatively high level of marker diversity was found, suggestive of multiple introductions of the organism to the region, or more likely ongoing endemicity. There was no evidence of ongoing outbreaks or transmission, and the bulk of cases were likely due to interaction between human activity and the environment (e.g. hunting/trapping activities).

SIGNIFICANCE AND IMPACT OF THE STUDY

This study reveals the diversity of Francisella tularensis in eastern Canada using multiple-locus variable-number tandem-repeat analysis. It was initiated to further the understanding of the species within North America as previous studies elucidating the diversity and phylogeography of the species have consisted mostly of specimens from the United States. Type A tularaemia, the most life-threatening subtype of the species and a Category A biothreat agent, is restricted to North America, and this study serves to broaden the knowledge of the epidemiology and diversity of the organism.

Francisella tularensis subsp. tularensis group A.I, United States

We used whole-genome analysis and subsequent characterization of geographically diverse strains using new genetic signatures to identify distinct subgroups within Francisella tularensis subsp. tularensis group A.I: A.I.3, A.I.8, and A.I.12. These subgroups exhibit complex phylogeographic patterns within North America. The widest distribution was observed for A.I.12, which suggests an adaptive advantage.

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.

TaqMan real-time PCR assays for single-nucleotide polymorphisms which identify Francisella tularensis and its subspecies and subpopulations

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

Population analysis of Streptococcus suis isolates from slaughtered swine by use of minimum core genome sequence typing

Streptococcus suis, an important zoonotic pathogen, is a highly diverse species with only a subset of strains that cause disease in humans. Our previous study proposed a minimum core genome (MCG) sequence typing method and defined seven MCG groups, with MCG group 1 as the prevalent group causing human infections. In this study, we identified a set of 10 single nucleotide polymorphisms (SNPs) distributed in six genes that were used to identify the seven MCG groups. The 10 SNPs were typed for 179 S. suis isolates collected from slaughtered pigs. The most prevalent groups among the tested isolates were MCG groups 6 and 7. Most of the isolates (147/179) were genotyped as mrp negative, epf negative, sly negative, and CDS2157 positive. The 179 isolates were also typed by multilocus sequence typing (MLST) and divided into 115 sequence types (STs), 111 of which were new. The 6 serotypes (29, 11, 5, 12, 30, and 2) represented 72.3% of the serotyped isolates. Our data show that the typing assay facilitates the application of genome data to the surveillance of S. suis.

Comparative analyses of a putative Francisella conjugative element

A large circular plasmid detected in Francisella novicida-like strain PA10-7858, designated pFNPA10, was sequenced completely and analyzed. This 41,013-bp plasmid showed no homology to any of the previously sequenced Francisella plasmids and was 8-10 times larger in size than them. A total of 57 ORFs were identified within pFNPA10 and at least 9 of them encoded putative proteins with homology to different conjugal transfer proteins. The presence of iteron-like direct repeats and an ORF encoding a putative replication protein within pFNPA10 suggested that it replicated by the theta mode. Phylogenetic analyses indicated that pFNPA10 had no near neighbors in the databases and that it may have originated within an environmental Francisella lineage. Based on its features, pFNPA10 appears to be a novel extra-chromosomal genetic element within the genus Francisella. The suitability of pFNPA10 as a vector for transformation of species of Francisella by conjugation remains to be explored.

Role of pathogenicity determinant protein C (PdpC) in determining the virulence of the Francisella tularensis subspecies tularensis SCHU

Francisella tularensis subspecies tularensis, the etiological agent of tularemia, is highly pathogenic to humans and animals. However, the SCHU strain of F. tularensis SCHU P0 maintained by passaging in artificial media has been found to be attenuated. To better understand the molecular mechanisms behind the pathogenicity of F. tularensis SCHU, we attempted to isolate virulent bacteria by serial passages in mice. SCHU P5 obtained after 5th passages in mice remained avirulent, while SCHU P9 obtained after 9th passages was completely virulent in mice. Moreover, SCHU P9 grew more efficiently in J774.1 murine macrophages compared with that in the less pathogenic SCHU P0 and P5. Comparison of the nucleotide sequences of the whole genomes of SCHU P0, P5, and P9 revealed only 1 nucleotide difference among P0, P5 and P9 in 1 of the 2 copies of pathogenicity determinant protein C (pdpC) gene. An adenine residue deletion was observed in the pdpC1 gene of SCHU P0, P5, and P9 and in the pdpC2 gene of SCHU P0, and P5, while P9 was characterized by the wild type pdpC2 gene. Thus, SCHU P0 and P5 expressed only truncated forms of PdpC protein, while SCHU P9 expressed both wild type and truncated versions. To validate the pathogenicity of PdpC, both copies of the pdpC gene in SCHU P9 have been inactivated by Targetron mutagenesis. SCHU P9 mutants with inactivated pdpC gene showed low intracellular growth in J774.1 cells and did not induce severe disease in experimentally infected mice, while virulence of the mutants was restored by complementation with expression of the intact PdpC. These results demonstrate that PdpC is crucial in determining the virulence of F. tularensis SCHU.

Virulence difference between the prototypic Schu S4 strain (A1a) and Francisella tularensis A1a, A1b, A2 and type B strains in a murine model of infection

BACKGROUND

The use of prototypic strains is common among laboratories studying infectious agents as it promotes consistency for data comparability among and between laboratories. Schu S4 is the prototypic virulent strain of Francisella tularensis and has been used extensively as such over the past six decades. Studies have demonstrated virulence differences among the two clinically relevant subspecies of F. tularensis, tularensis (type A) and holarctica (type B) and more recently between type A subpopulations (A1a, A1b and A2). Schu S4 belongs to the most virulent subspecies of F. tularensis, subspecies tularensis.

METHODS

In this study, we investigated the relative virulence of Schu S4 in comparison to A1a, A1b, A2 and type B strains using a temperature-based murine model of infection. Mice were inoculated intradermally and a hypothermic drop point was used as a surrogate for death. Survival curves and the length of temperature phases were compared for all infections. Bacterial burdens were also compared between the most virulent type A subpopulation, A1b, and Schu S4 at drop point.

RESULTS

Survival curve comparisons demonstrate that the Schu S4 strain used in this study resembles the virulence of type B strains, and is significantly less virulent than all other type A (A1a, A1b and A2) strains tested. Additionally, when bacterial burdens were compared between mice infected with Schu S4 or MA00-2987 (A1b) significantly higher burdens were present in the blood and spleen of mice infected with MA00-2987.

CONCLUSIONS

The knowledge gained from using Schu S4 as a prototypic virulent strain has unquestionably advanced the field of tularemia research. The findings of this study, however, indicate that careful consideration of F. tularensis strain selection must occur when the overall virulence of the strain used could impact the outcome and interpretation of results.

Non-synonymous variations in cancer and their effects on the human proteome: workflow for NGS data biocuration and proteome-wide analysis of TCGA data

BACKGROUND

Next-generation sequencing (NGS) technologies have resulted in petabytes of scattered data, decentralized in archives, databases and sometimes in isolated hard-disks which are inaccessible for browsing and analysis. It is expected that curated secondary databases will help organize some of this Big Data thereby allowing users better navigate, search and compute on it.

RESULTS

To address the above challenge, we have implemented a NGS biocuration workflow and are analyzing short read sequences and associated metadata from cancer patients to better understand the human variome. Curation of variation and other related information from control (normal tissue) and case (tumor) samples will provide comprehensive background information that can be used in genomic medicine research and application studies. Our approach includes a CloudBioLinux Virtual Machine which is used upstream of an integrated High-performance Integrated Virtual Environment (HIVE) that encapsulates Curated Short Read archive (CSR) and a proteome-wide variation effect analysis tool (SNVDis). As a proof-of-concept, we have curated and analyzed control and case breast cancer datasets from the NCI cancer genomics program - The Cancer Genome Atlas (TCGA). Our efforts include reviewing and recording in CSR available clinical information on patients, mapping of the reads to the reference followed by identification of non-synonymous Single Nucleotide Variations (nsSNVs) and integrating the data with tools that allow analysis of effect nsSNVs on the human proteome. Furthermore, we have also developed a novel phylogenetic analysis algorithm that uses SNV positions and can be used to classify the patient population. The workflow described here lays the foundation for analysis of short read sequence data to identify rare and novel SNVs that are not present in dbSNP and therefore provides a more comprehensive understanding of the human variome. Variation results for single genes as well as the entire study are available from the CSR website (http://hive.biochemistry.gwu.edu/dna.cgi?cmd=csr).

CONCLUSIONS

Availability of thousands of sequenced samples from patients provides a rich repository of sequence information that can be utilized to identify individual level SNVs and their effect on the human proteome beyond what the dbSNP database provides.

Genome-wide single nucleotide polymorphism-based assay for high-resolution epidemiological analysis of the methicillin-resistant Staphylococcus aureus hospital clone EMRSA-15

The EMRSA-15 clone is a major cause of nosocomial methicillin-resistant Staphylococcus aureus (MRSA) infections in the UK and elsewhere but existing typing methodologies have limited capacity to discriminate closely related strains, and are often poorly reproducible between laboratories. Here, we report the design, development and validation of a genome-wide single nucleotide polymorphism (SNP) typing method and compare it to established methods for typing of EMRSA-15. In order to identify discriminatory SNPs, the genomes of 17 EMRSA-15 strains, selected to represent the breadth of genotypic and phenotypic diversity of EMRSA-15 isolates in Scotland, were determined and phylogenetic reconstruction was carried out. In addition to 17 phylogenetically informative SNPs, five binary markers were included to form the basis of an EMRSA-15 genotyping assay. The SNP-based typing assay was as discriminatory as pulsed-field gel electrophoresis, and significantly more discriminatory than staphylococcal protein A (spa) typing for typing of a representative panel of diverse EMRSA-15 strains, isolates from two EMRSA-15 hospital outbreak investigations, and a panel of bacteraemia isolates obtained in healthcare facilities in the east of Scotland during a 12-month period. The assay is a rapid, and reproducible approach for epidemiological analysis of EMRSA-15 clinical isolates in Scotland. Unlike established methods the DNA sequence-based method is ideally suited for inter-laboratory comparison of identified genotypes, and its flexibility lends itself to supplementation with additional SNPs or markers for the identification of novel S. aureus strains in other regions of the world.

Adaptive mutations and replacements of virulence traits in the Escherichia coli O104:H4 outbreak population

The sequencing of highly virulent Escherichia coli O104:H4 strains isolated during the outbreak of bloody diarrhea and hemolytic uremic syndrome in Europe in 2011 revealed a genome that contained a Shiga toxin encoding prophage and a plasmid encoding enteroaggregative fimbriae. Here, we present the draft genome sequence of a strain isolated in Sweden from a patient who had travelled to Tunisia in 2010 (E112/10) and was found to differ from the outbreak strains by only 38 SNPs in non-repetitive regions, 16 of which were mapped to the branch to the outbreak strain. We identified putatively adaptive mutations in genes for transporters, outer surface proteins and enzymes involved in the metabolism of carbohydrates. A comparative analysis with other historical strains showed that E112/10 contained Shiga toxin prophage genes of the same genotype as the outbreak strain, while these genes have been replaced by a different genotype in two otherwise very closely related strains isolated in the Republic of Georgia in 2009. We also present the genome sequences of two enteroaggregative E. coli strains affiliated with phylogroup A (C43/90 and C48/93) that contain the agg genes for the AAF/I-type fimbriae characteristic of the outbreak population. Interestingly, C43/90 also contained a tet/mer antibiotic resistance island that was nearly identical in sequence to that of the outbreak strain, while the corresponding island in the Georgian strains was most similar to E. coli strains of other serotypes. We conclude that the pan-genome of the outbreak population is shared with strains of the A phylogroup and that its evolutionary history is littered with gene replacement events, including most recently independent acquisitions of antibiotic resistance genes in the outbreak strains and its nearest neighbors. The results are summarized in a refined evolutionary model for the emergence of the O104:H4 outbreak population.

German Francisella tularensis isolates from European brown hares (Lepus europaeus) reveal genetic and phenotypic diversity

BACKGROUND

Tularemia is a zoonotic disease caused by Francisella tularensis that has been found in many different vertebrates. In Germany most human infections are caused by contact with infected European brown hares (Lepus europaeus). The aim of this study was to elucidate the epidemiology of tularemia in hares using phenotypic and genotypic characteristics of F. tularensis.

RESULTS

Cultivation of F. tularensis subsp. holarctica bacteria from organ material was successful in 31 of 52 hares that had a positive PCR result targeting the Ft-M19 locus. 17 isolates were sensitive to erythromycin and 14 were resistant. Analysis of VNTR loci (Ft-M3, Ft-M6 and Ft-M24), INDELs (Ftind33, Ftind38, Ftind49, RD23) and SNPs (B.17, B.18, B.19, and B.20) was shown to be useful to investigate the genetic relatedness of Francisella strains in this set of strains. The 14 erythromycin resistant isolates were assigned to clade B.I, and 16 erythromycin sensitive isolates to clade B.IV and one isolate was found to belong to clade B.II. MALDI-TOF mass spectrometry (MS) was useful to discriminate strains to the subspecies level.

CONCLUSIONS

F. tularensis seems to be a re-emerging pathogen in Germany. The pathogen can easily be identified using PCR assays. Isolates can also be identified within one hour using MALDI-TOF MS in laboratories where specific PCR assays are not established. Further analysis of strains requires genotyping tools. The results from this study indicate a geographical segregation of the phylogenetic clade B.I and B.IV, where B.I strains localize primarily within eastern Germany and B.IV strains within western Germany. This phylogeographical pattern coincides with the distribution of biovar I (erythromycin sensitive) and biovar II (erythromycin resistance) strains. When time and costs are limiting parameters small numbers of isolates can be analysed using PCR assays combined with DNA sequencing with a focus on genetic loci that are most likely discriminatory among strains found in a specific area. In perspective, whole genome data will have to be investigated especially when terrorist attack strains need to be tracked to their genetic and geographical sources.

Rapid focused sequencing: a multiplexed assay for simultaneous detection and strain typing of Bacillus anthracis, Francisella tularensis, and Yersinia pestis

Background

The intentional release of Bacillus anthracis in the United States in 2001 has heightened concern about the use of pathogenic microorganisms in bioterrorism attacks. Many of the deadliest bacteria, including the Class A Select Agents Bacillus anthracis, Francisella tularensis, and Yersinia pestis, are highly infectious via the pulmonary route when released in aerosolized form. Hence, rapid, sensitive, and reliable methods for detection of these biothreats and characterization of their potential impact on the exposed population are of critical importance to initiate and support rapid military, public health, and clinical responses.

Methodology/Principal Findings

We have developed microfluidic multiplexed PCR and sequencing assays based on the simultaneous interrogation of three pathogens per assay and ten loci per pathogen. Microfluidic separation of amplified fluorescently labeled fragments generated characteristic electrophoretic signatures for identification of each agent. The three sets of primers allowed significant strain typing and discrimination from non-pathogenic closely-related species and environmental background strains based on amplicon sizes alone. Furthermore, sequencing of the 10 amplicons per pathogen, termed “Rapid Focused Sequencing,” allowed an even greater degree of strain discrimination and, in some cases, can be used to determine virulence. Both amplification and sequencing assays were performed in microfluidic biochips developed for fast thermal cycling and requiring 7 µL per reaction. The 30-plex sequencing assay resulted in genotypic resolution of 84 representative strains belonging to each of the three biothreat species.

Conclusions/Significance

The microfluidic multiplexed assays allowed identification and strain differentiation of the biothreat agents Bacillus anthracis, Francisella tularensis, and Yersinia pestis and clear discrimination from closely-related species and several environmental background strains. The assays may be extended to detect a large number of pathogens, are applicable to the evaluation of both environmental and clinical samples, and have the potential to be applied in military, public health, and clinical diagnostic settings.

Use of temperature for standardizing the progression of Francisella tularensis in mice

The study of infectious agents, their pathogenesis, the host response and the evaluation of newly developed countermeasures often requires the use of a living system. Murine models are frequently used to undertake such investigations with the caveat that non-biased measurements to assess the progression of infection are underutilized. Instead, murine models predominantly rely on symptomology exhibited by the animal to evaluate the state of the animal's health and to determine when euthanasia should be performed. In this study, we used subcutaneous temperature as a non-subjective measurement to follow and compare infection in mice inoculated with Francisella tularensis, a Gram-negative pathogen that produces an acute and fatal illness in mice. A reproducible temperature pattern defined by three temperature phases (normal, febrile and hypothermic) was identified in all mice infected with F. tularensis, regardless of the infecting strain. More importantly and for the first time a non-subjective, ethical, and easily determined surrogate endpoint for death based on a temperature, termed drop point, was identified and validated with statistical models. In comparative survival curve analyses for F. tularensis strains with differing virulence, the drop point temperature yielded the same results as those obtained using observed time to death. Incorporation of temperature measurements to evaluate F. tularensis was standardized based on statistical models to provide a new level of robustness for comparative analyses in mice. These findings should be generally applicable to other pathogens that produce acute febrile disease in animal models and offers an important tool for understanding and following the infection process.

Genomic comparison between a virulent type A1 strain of Francisella tularensis and its attenuated O-antigen mutant

We report the complete genome sequences of TI0902, a highly virulent type A1 strain, and TIGB03, a related, attenuated chemical mutant strain. Compared to the wild type, the mutant strain had 45 point mutations and a 75.9-kb duplicated region that had not been previously observed in Francisella species.

Melt analysis of mismatch amplification mutation assays (Melt-MAMA): a functional study of a cost-effective SNP genotyping assay in bacterial models

Single nucleotide polymorphisms (SNPs) are abundant in genomes of all species and biologically informative markers extensively used across broad scientific disciplines. Newly identified SNP markers are publicly available at an ever-increasing rate due to advancements in sequencing technologies. Efficient, cost-effective SNP genotyping methods to screen sample populations are in great demand in well-equipped laboratories, but also in developing world situations. Dual Probe TaqMan assays are robust but can be cost-prohibitive and require specialized equipment. The Mismatch Amplification Mutation Assay, coupled with melt analysis (Melt-MAMA), is flexible, efficient and cost-effective. However, Melt-MAMA traditionally suffers from high rates of assay design failures and knowledge gaps on assay robustness and sensitivity. In this study, we identified strategies that improved the success of Melt-MAMA. We examined the performance of 185 Melt-MAMAs across eight different pathogens using various optimization parameters. We evaluated the effects of genome size and %GC content on assay development. When used collectively, specific strategies markedly improved the rate of successful assays at the first design attempt from ∼50% to ∼80%. We observed that Melt-MAMA accurately genotypes across a broad DNA range (∼100 ng to ∼0.1 pg). Genomic size and %GC content influence the rate of successful assay design in an independent manner. Finally, we demonstrated the versatility of these assays by the creation of a duplex Melt-MAMA real-time PCR (two SNPs) and conversion to a size-based genotyping system, which uses agarose gel electrophoresis. Melt-MAMA is comparable to Dual Probe TaqMan assays in terms of design success rate and accuracy. Although sensitivity is less robust than Dual Probe TaqMan assays, Melt-MAMA is superior in terms of cost-effectiveness, speed of development and versatility. We detail the parameters most important for the successful application of Melt-MAMA, which should prove useful to the wider scientific community.

DNA fingerprinting of Mycobacterium tuberculosis: from phage typing to whole-genome sequencing

Current typing methods for Mycobacterium tuberculosis complex evolved from simple phenotypic approaches like phage typing and drug susceptibility profiling to DNA-based strain typing methods, such as IS6110-restriction fragment length polymorphisms (RFLP) and variable number of tandem repeats (VNTR) typing. Examples of the usefulness of molecular typing are source case finding and epidemiological linkage of tuberculosis (TB) cases, international transmission of MDR/XDR-TB, the discrimination between endogenous reactivation and exogenous re-infection as a cause of relapses after curative treatment of tuberculosis, the evidence of multiple M. tuberculosis infections, and the disclosure of laboratory cross-contaminations. Simultaneously, phylogenetic analyses were developed based on single nucleotide polymorphisms (SNPs), genomic deletions usually referred to as regions of difference (RDs) and spoligotyping which served both strain typing and phylogenetic analysis. National and international initiatives that rely on the application of these typing methods have brought significant insight into the molecular epidemiology of tuberculosis. However, current DNA fingerprinting methods have important limitations. They can often not distinguish between genetically closely related strains and the turn-over of these markers is variable. Moreover, the suitability of most DNA typing methods for phylogenetic reconstruction is limited as they show a high propensity of convergent evolution or misinfer genetic distances. In order to fully explore the possibilities of genotyping in the molecular epidemiology of tuberculosis and to study the phylogeny of the causative bacteria reliably, the application of whole-genome sequencing (WGS) analysis for all M. tuberculosis isolates is the optimal, although currently still a costly solution. In the last years WGS for typing of pathogens has been explored and yielded important additional information on strain diversity in comparison to the classical DNA typing methods. With the ongoing cost reduction of DNA sequencing it is possible that WGS will become the sole diagnostic tool in the secondary laboratory diagnosis of tuberculosis for identification, drug susceptibility testing and genetic characterization.

Identification and subtyping of Francisella by pyrosequencing and signature matching of 16S rDNA fragments

AIMS

To analyse the V1 region of the 16S rDNA gene by a universal pyrosequencing protocol to identify and subtype Francisella in 31 strains from a repository collection and 96 patient isolates.

METHODS AND RESULTS

Pyrosequencing was used to determine the nucleotide sequence of PCR amplification products of the variable region (V1) of the 16S rDNA from 31 repository strains and 96 isolates from Swedish patients with ulceroglandular tularaemia. Pyrosequencing resulted in a 37 nucleotide sequence, specific for Francisella sp., for all repository strains and patient samples analysed. In addition, the isolates could be divided into two groups based on the analysis of a single nucleotide polymorphism in the sequence: one group included Francisella tularensis ssp. tularensis, ssp. holarctica and ssp. mediasiatica, whereas the other group included Francisella tularensis ssp. novicida and other species of Francisella. The analysis of samples taken from patients suffering from ulceroglandular tularaemia revealed that all isolates belonged to the first group comprising subspecies of F. tularensis virulent for humans.

CONCLUSIONS

The pyrosequencing analysis of the 16S rDNA V1 is a useful molecular tool for the rapid identification of suspected isolates of Francisella sp. in clinical or environmental samples.

SIGNIFICANCE AND IMPACT OF THE STUDY

Virulent F. tularensis ssp. causing ulceroglandular tularaemia, or those with a potential to be used in a bioterrorism event, could rapidly be discriminated from subspecies less virulent for humans.

An integer programming approach to DNA sequence assembly

De novo sequence assembly is a ubiquitous combinatorial problem in all DNA sequencing technologies. In the presence of errors in the experimental data, the assembly problem is computationally challenging, and its solution may not lead to a unique reconstruct. The enumeration of all alternative solutions is important in drawing a reliable conclusion on the target sequence, and is often overlooked in the heuristic approaches that are currently available. In this paper, we develop an integer programming formulation and global optimization solution strategy to solve the sequence assembly problem with errors in the data. We also propose an efficient technique to identify all alternative reconstructs. When applied to examples of sequencing-by-hybridization, our approach dramatically increases the length of DNA sequences that can be handled with global optimality certificate to over 10,000, which is more than 10 times longer than previously reported. For some problem instances, alternative solutions exhibited a wide range of different ability in reproducing the target DNA sequence. Therefore, it is important to utilize the methodology proposed in this paper in order to obtain all alternative solutions to reliably infer the true reconstruct. These alternative solutions can be used to refine the obtained results and guide the design of further experiments to correctly reconstruct the target DNA sequence.

Recombination and population structure in Salmonella enterica

Salmonella enterica is a bacterial pathogen that causes enteric fever and gastroenteritis in humans and animals. Although its population structure was long described as clonal, based on high linkage disequilibrium between loci typed by enzyme electrophoresis, recent examination of gene sequences has revealed that recombination plays an important evolutionary role. We sequenced around 10% of the core genome of 114 isolates of enterica using a resequencing microarray. Application of two different analysis methods (Structure and ClonalFrame) to our genomic data allowed us to define five clear lineages within S. enterica subspecies enterica, one of which is five times older than the other four and two thirds of the age of the whole subspecies. We show that some of these lineages display more evidence of recombination than others. We also demonstrate that some level of sexual isolation exists between the lineages, so that recombination has occurred predominantly between members of the same lineage. This pattern of recombination is compatible with expectations from the previously described ecological structuring of the enterica population as well as mechanistic barriers to recombination observed in laboratory experiments. In spite of their relatively low level of genetic differentiation, these lineages might therefore represent incipient species.

Everything at once: comparative analysis of the genomes of bacterial pathogens

The sum of unique genes in all genomes of a bacterial species is referred to as the pan-genome and is comprised of variably absent or present accessory genes and universally present core genes. The accessory genome is an important source of genetic variability in bacterial populations, allowing sub-populations of bacteria to better adapt to specific niches. Such subgroups may themselves have a relatively stable core genome that may influence host preference, virulence, or an association with specific disease syndromes. The core genome provides a useful means of phylogenetic reconstruction as well as contributing to phenotypic heterogeneity. Variation within the pan-genome forms the basis of comparative genotyping techniques, which have evolved alongside technology. Current high-throughput sequencing platforms have created an unprecedented opportunity for comparisons among multiple, closely related genomes. The computer algorithms and software for such comparisons continue to evolve and promise exciting advances in the world of bacterial comparative genomics. We review genotyping techniques based upon phenotypic traits, both core and accessory genomes, and look at some of the software programs currently available to perform whole-genome comparative analyses.

Differential chitinase activity and production within Francisella species, subspecies, and subpopulations

Genotyping of Francisella tularensis (A1a, A1b, A2, and type B) and Francisella novicida has identified multiple differences between species and among F. tularensis subspecies and subpopulations. Variations in virulence, geographic distribution, and ecology are also known to exist among this group of bacteria, despite the >95% nucleotide identity in their genomes. This study expands the description of phenotypic differences by evaluating the ability of F. tularensis and F. novicida to degrade chitin analogs and produce active chitinases. Endochitinase activities were observed to vary among F. tularensis and F. novicida strains. The activity observed for F. tularensis strains was predominantly associated with whole-cell lysates, while the chitinase activity of F. novicida localized to the culture supernatant. In addition, the overall level of chitinase activity differed among the subpopulations of F. tularensis and between the species. Bioinformatic analyses identified two new putative chitinase genes (chiC and chiD), as well as the previously described chiA and chiB. However, the presence of these four open reading frames as intact genes or pseudogenes was found to differ between Francisella species and F. tularensis subspecies and subpopulations. Recombinant production of the putative chitinases and enzymatic evaluations revealed ChiA, ChiB, ChiC, and ChiD possessed dissimilar chitinase activities. These biochemical studies coupled with bioinformatic analyses and the evaluation of chiA and chiC knockouts in F. tularensis A1 and A2 strains, respectively, provided a molecular basis to explain the differential chitinase activities observed among the species and subpopulations of Francisella.

Monitoring the long-term molecular epidemiology of the pneumococcus and detection of potential 'vaccine escape' strains

BACKGROUND

While the pneumococcal protein conjugate vaccines reduce the incidence in invasive pneumococcal disease (IPD), serotype replacement remains a major concern. Thus, serotype-independent protection with vaccines targeting virulence genes, such as PspA, have been pursued. PspA is comprised of diverse clades that arose through recombination. Therefore, multi-locus sequence typing (MLST)-defined clones could conceivably include strains from multiple PspA clades. As a result, a method is needed which can both monitor the long-term epidemiology of the pneumococcus among a large number of isolates, and analyze vaccine-candidate genes, such as pspA, for mutations and recombination events that could result in 'vaccine escape' strains.

METHODOLOGY

We developed a resequencing array consisting of five conserved and six variable genes to characterize 72 pneumococcal strains. The phylogenetic analysis of the 11 concatenated genes was performed with the MrBayes program, the single nucleotide polymorphism (SNP) analysis with the DNA Sequence Polymorphism program (DnaSP), and the recombination event analysis with the recombination detection package (RDP).

RESULTS

The phylogenetic analysis correlated with MLST, and identified clonal strains with unique PspA clades. The DnaSP analysis correlated with the serotype-specific diversity detected using MLST. Serotypes associated with more than one ST complex had a larger degree of sequence polymorphism than a serotype associated with one ST complex. The RDP analysis confirmed the high frequency of recombination events in the pspA gene.

CONCLUSIONS

The phylogenetic tree correlated with MLST, and detected multiple PspA clades among clonal strains. The genetic diversity of the strains and the frequency of recombination events in the mosaic gene, pspA were accurately assessed using the DnaSP and RDP programs, respectively. These data provide proof-of-concept that resequencing arrays could play an important role within research and clinical laboratories in both monitoring the molecular epidemiology of the pneumococcus and detecting 'vaccine escape' strains among vaccine-candidate genes.

BIGSdb: Scalable analysis of bacterial genome variation at the population level

Background

The opportunities for bacterial population genomics that are being realised by the application of parallel nucleotide sequencing require novel bioinformatics platforms. These must be capable of the storage, retrieval, and analysis of linked phenotypic and genotypic information in an accessible, scalable and computationally efficient manner.

Results

The Bacterial Isolate Genome Sequence Database (BIGSDB) is a scalable, open source, web-accessible database system that meets these needs, enabling phenotype and sequence data, which can range from a single sequence read to whole genome data, to be efficiently linked for a limitless number of bacterial specimens. The system builds on the widely used mlstdbNet software, developed for the storage and distribution of multilocus sequence typing (MLST) data, and incorporates the capacity to define and identify any number of loci and genetic variants at those loci within the stored nucleotide sequences. These loci can be further organised into 'schemes' for isolate characterisation or for evolutionary or functional analyses. Isolates and loci can be indexed by multiple names and any number of alternative schemes can be accommodated, enabling cross-referencing of different studies and approaches. LIMS functionality of the software enables linkage to and organisation of laboratory samples. The data are easily linked to external databases and fine-grained authentication of access permits multiple users to participate in community annotation by setting up or contributing to different schemes within the database. Some of the applications of BIGSDB are illustrated with the genera Neisseria and Streptococcus.

The BIGSDB source code and documentation are available at http://pubmlst.org/software/database/bigsdb/.

Conclusions

Genomic data can be used to characterise bacterial isolates in many different ways but it can also be efficiently exploited for evolutionary or functional studies. BIGSDB represents a freely available resource that will assist the broader community in the elucidation of the structure and function of bacteria by means of a population genomics approach.

Subpopulations of Francisella tularensis ssp. tularensis and holarctica: identification and associated epidemiology

Tularemia is primarily caused by two subspecies of Francisella tularensis worldwide, ssp. tularensis (type A) and ssp. holarctica (type B), which were originally delineated by phenotypic differences. Application of molecular typing methods to investigate population structure of F. tularensis has confirmed that categorizing the two subspecies via phenotypic characteristics corresponds with genotypic differentiation. In addition, genotyping methods have demonstrated that both subspecies, type A and type B, can be further distinguished into subpopulations and, in some cases, biological relevance has been ascribed to these identified subpopulations. Genetic variation among both type A and type B subpopulations has been shown to correlate with differences in geographic distribution and has also been coupled to distinct ecological niches, animal hosts and replication foci. Among type A subpopulations, strain variation is linked to differing clinical manifestations in humans and virulence in mice. This article will highlight our current understanding of F. tularensis subpopulations, including methods for their detection, their observed epidemiologic differences, implications for public health and basic research programs, as well as future challenges yet to be solved.

Ecological niche modeling of Francisella tularensis subspecies and clades in the United States

Two subspecies of Francisella tularensis are recognized: F. tularensis subsp. tularensis (type A) and F. tularensis subsp. holartica (type B). Type A has been subdivided further into A1a, A1b, and A2, which differ geographically and clinically. The aim of this work was to determine whether or not differences among subspecies and clades translate into distinct ecological niches. We used 223 isolates from humans and wildlife representing all six genotypes (type A, B, A1, A2, A1a, or A1b). Ecological-niche models were built independently for each genotype, using the genetic algorithm for rule-set prediction. The resulting models were compared using a non-parametric multivariate analysis-of-variance method. A1 and A2 are ecologically distinct, supporting the previously observed geographic division, whereas ecological niches for types A and B overlapped notably but A1a and A1b displayed no appreciable differences in their ecological niches.

Virulence differences among Francisella tularensis subsp. tularensis clades in mice

Francisella tularensis subspecies tularensis (type A) and holarctica (type B) are of clinical importance in causing tularemia. Molecular typing methods have further separated type A strains into three genetically distinct clades, A1a, A1b and A2. Epidemiological analyses of human infections in the United States suggest that A1b infections are associated with a significantly higher mortality rate as compared to infections caused by A1a, A2 and type B. To determine if genetic differences as defined by molecular typing directly correlate with differences in virulence, A1a, A1b, A2 and type B strains were compared in C57BL/6 mice. Here we demonstrate significant differences between survival curves for infections caused by A1b versus A1a, A2 and type B, with A1b infected mice dying earlier than mice infected with A1a, A2 or type B; these results were conserved among multiple strains. Differences were also detected among type A clades as well as between type A clades and type B with respect to bacterial burdens, and gross anatomy in infected mice. Our results indicate that clades defined within F. tularensis subsp. tularensis by molecular typing methods correlate with virulence differences, with A1b strains more virulent than A1a, A2 and type B strains. These findings indicate type A strains are not equivalent with respect to virulence and have important implications for public health as well as basic research programs.