Genome sequencing shows that European isolates of Francisella tularensis subspecies tularensis are almost identical to US laboratory strain Schu S4

Francisella tularensis PCR detection in Cape hares (Lepus capensis) and wild rabbits (Oryctolagus cuniculus) in Algeria

Tularemia is a zoonosis caused by the bacterium Francisella tularensis. Leporids are primary sources of human infections in the northern hemisphere. Africa is classically considered free of tularemia, but recent data indicate that this dogma might be wrong. We assessed the presence of this disease in wild leporids in Algeria. Between 2014 and 2018, we collected 74 leporids carcasses from spontaneously dead or hunted animals. Francisella tularensis DNA was detected by specific real-time PCR tests in 7/36 (19.44%) Cape hares (Lepus capensis) and 5/38 (13.15%) wild rabbits (Oryctolagus cuniculus). Known tularemia arthropod vectors infested half of the PCR-positive animals. At necropsy, F. tularensis-infected animals presented with an enlarged spleen (n = 12), enlarged adrenal glands (12), liver discoloration (12), hemorrhages (11), and pneumonia (11). Immunohistological examination of liver tissue from one animal was compatible with the presence of F. tularensis. Our study demonstrates the existence of tularemia in lagomorphs in Algeria. It should encourage investigations to detect this disease among the human population of this country.

Francisella tularensis caused cervical lymphadenopathy in little children after a tick bite: Two case reports and a short literature review

Although Francisella (F.) tularensis is a well-described and understood zoonotic pathogen, its importance in Central Europe is relatively minor and, as such, tularaemia may be missed in the differential diagnosis. The annual incidence of tularaemia in the Czech Republic is relatively stable with up to 100 reported cases per year, except in the epidemic years 1998 and 1999 with 225 and 222 reported cases, respectively. It is, however, higher in comparison with the neighbouring countries. The common route of transmission in Central Europe is handling infected animals. Tularaemia is not commonly recognized as a tick-borne disease. Here we report two rare cases of a tick bite-associated ulceroglandular form of tularaemia in 2.5-year-old and 6.5-year-old children presenting with cervical lymphadenopathy. The unusual and interesting features of those cases are the young age and relatively uncommon route of transmission suggesting possible changes in the epidemiology of tularaemia in the Czech Republic. Therefore, the infection with F. tularensis should be considered in the differential diagnosis after a tick bite even in infants.

Cross-Species Proteomic Comparison of Outer Membrane Vesicles and Membranes of Francisella tularensis subsp. tularensis versus subsp. holarctica

Release of outer membrane vesicles (OMV) is an important phenomenon in Gram-negative bacteria playing multiple roles in their lifestyle, including in relation to virulence and host-pathogen interaction. Francisella tularensis, unlike other bacteria, releases unusually shaped, tubular OMV. We present a proteomic comparison of OMV and membrane fractions from two F. tularensis strains: moderately virulent subsp. holarctica strain FSC200 and highly virulent subsp. tularensis strain SchuS4. Proteomic comparison studies routinely evaluate samples from the same proteome, but sometimes we must compare samples from closely related organisms. This raises quantification issues. We propose a novel approach to cross-species proteomic comparison based on an intersection protein database from the individual single-species databases. This is less prone to quantification errors arising from differences in the sequences. Consecutively comparing subproteomes of OMV and membranes of the two strains allows distinguishing differences in relative protein amounts caused by global expression changes from those caused by preferential protein packing to OMV or membranes. Among the proteins most differently packed into OMV between the two strains, we detected proteins involved in biosynthesis and metabolism of bacterial envelope components like O-antigen, lipid A, phospholipids, and fatty acids, as well as some major structural outer membrane proteins. The data are available via ProteomeXchange with identifier PXD022406.

Tularemia Goes West: Epidemiology of an Emerging Infection in Austria

The zoonotic disease tularemia is caused by the Gram-negative bacterium Francisella tularensis, with the two major subspecies tularensis and holarctica being responsible for infections in humans and animals. The F. tularensis subspecies holarctica is less virulent and prevalent in Europe and Asia. Over the last few centuries, few epidemic outbreaks and low numbers of infections have been registered in the eastern part of Austria, specifically in the provinces of Lower Austria, Burgenland, and Styria. The reported infections were mostly associated with hunting hares and the skinning of carcasses. Within the last decade, ticks have been identified as important vectors in Tyrol and served as first evidence for the spread of F. tularensis to Western Austria. In 2018, the pathogen was detected in hares in the provinces of Tyrol, Vorarlberg, and Salzburg. We presume that F. tularensis is now established in most regions of Austria, and that the investigation of potential host and vector animals should be spotlighted by public institutions. Tularemia in humans presents with various clinical manifestations. As glandular, ulceroglandular, and typhoidal forms occur in Austria, this infectious disease should be considered as a differential diagnosis of unknown fever.

Using affinity propagation clustering for identifying bacterial clades and subclades with whole-genome sequences of Francisella tularensis

By combining a reference-independent SNP analysis and average nucleotide identity (ANI) with affinity propagation clustering (APC), we developed a significantly improved methodology allowing resolving phylogenetic relationships, based on objective criteria. These bioinformatics tools can be used as a general ruler to determine phylogenetic relationships and clustering of bacteria, exemplary done with Francisella (F.) tularensis. Molecular epidemiology of F. tularensis is currently assessed mostly based on laboratory methods and molecular analysis. The high evolutionary stability and the clonal nature makes Francisella ideal for subtyping with single nucleotide polymorphisms (SNPs). Sequencing and real-time PCR can be used to validate the SNP analysis. We investigate whole-genome sequences of 155 F. tularensis subsp. holarctica isolates. Phylogenetic testing was based on SNPs and average nucleotide identity (ANI) as reference independent, alignment-free methods taking small-scale and large-scale differences within the genomes into account. Especially the whole genome SNP analysis with kSNP3.0 allowed deciphering quite subtle signals of systematic differences in molecular variation. Affinity propagation clustering (APC) resulted in three clusters showing the known clades B.4, B.6, and B.12. These data correlated with the results of real-time PCR assays targeting canSNPs loci. Additionally, we detected two subtle sub-clusters. SplitsTree was used with standard-setting using the aligned SNPs from Parsnps. Together APC, HierBAPS, and SplitsTree enabled us to generate hypotheses about epidemiologic relationships between bacterial clusters and describing the distribution of isolates. Our data indicate that the choice of the typing technique can increase our understanding of the pathogenesis and transmission of diseases with the eventual for prevention. This is opening perspectives to be applied to other bacterial species. The data provide evidence that Germany might be the collision zone where the clade B.12, also known as the East European clade, overlaps with the clade B.6, also known as the Iberian clade. Described methods allow generating a new, more detailed perspective for F. tularensis subsp. holarctica phylogeny. These results may encourage to determine phylogenetic relationships and clustering of other bacteria the same way.

Tularemia as a waterborne disease: a review

Francisella tularensis is a Gram-negative, intracellular bacterium causing the zoonosis tularemia. This highly infectious microorganism is considered a potential biological threat agent. Humans are usually infected through direct contact with the animal reservoir and tick bites. However, tularemia cases also occur after contact with a contaminated hydro-telluric environment. Water-borne tularemia outbreaks and sporadic cases have occurred worldwide in the last decades, with specific clinical and epidemiological traits. These infections represent a major public health and military challenge. Human contaminations have occurred through consumption or use of F. tularensis-contaminated water, and various aquatic activities such as swimming, canyoning and fishing. In addition, in Sweden and Finland, mosquitoes are primary vectors of tularemia due to infection of mosquito larvae in contaminated aquatic environments. The mechanisms of F. tularensis survival in water may include the formation of biofilms, interactions with free-living amoebae, and the transition to a 'viable but nonculturable' state, but the relative contribution of these possible mechanisms remains unknown. Many new aquatic species of Francisella have been characterized in recent years. F. tularensis likely shares with these species an ability of long-term survival in the aquatic environment, which has to be considered in terms of tularemia surveillance and control.

Revisiting Francisella tularensis subsp. holarctica, Causative Agent of Tularemia in Germany With Bioinformatics: New Insights in Genome Structure, DNA Methylation and Comparative Phylogenetic Analysis

Francisella (F.) tularensis is a highly virulent, Gram-negative bacterial pathogen and the causative agent of the zoonotic disease tularemia. Here, we generated, analyzed and characterized a high quality circular genome sequence of the F. tularensis subsp. holarctica strain 12T0050 that caused fatal tularemia in a hare. Besides the genomic structure, we focused on the analysis of oriC, unique to the Francisella genus and regulating replication in and outside hosts and the first report on genomic DNA methylation of a Francisella strain. The high quality genome was used to establish and evaluate a diagnostic whole genome sequencing pipeline. A genotyping strategy for F. tularensis was developed using various bioinformatics tools for genotyping. Additionally, whole genome sequences of F. tularensis subsp. holarctica isolates isolated in the years 2008–2015 in Germany were generated. A phylogenetic analysis allowed to determine the genetic relatedness of these isolates and confirmed the highly conserved nature of F. tularensis subsp. holarctica.

Forensic Approaches to Detect Possible Agents of Bioterror

Many biological agents have been strategic pathogenic agents throughout history. Some have even changed history as a consequence of early discoveries of their use as weapons of war. Many of these bioagents can be easily isolated from the environment, and some have recently been genetically manipulated to become more pathogenic for biowarfare. However, it is difficult to determine accidental outbreaks of disease from intentional exposures. In this review, we examine how molecular tools have been used in combination with forensic research to resolve cases of unusual outbreaks and trace the source of the biocrime. New technologies are also discussed in terms of their crucial role impacting forensic science. The anthrax event of 2001 serves as an example of the real threat of bioterrorism and the employment of bioagents as weapons against a population. The Amerithrax investigation has given us lessons of the highest resolution possible with new technologies capable of distinguishing isolates at the base-pair level of sensitivity. In addition, we discuss the implications of proper sanitation to avoid waterborne diseases. The use of new methods in forensic science and health-related surveillance will be invaluable in determining the source of any new disease outbreak, and these data will allow for a quick response to any type of public health threat, whether accidental or purposely initiated.

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.

Comparative Transcriptional Analyses of Francisella tularensis and Francisella novicida

Francisella tularensis is composed of a number of subspecies with varied geographic distribution, host ranges, and virulence. In view of these marked differences, comparative functional genomics may elucidate some of the molecular mechanism(s) behind these differences. In this study a shared probe microarray was designed that could be used to compare the transcriptomes of Francisella tularensis subsp. tularensis Schu S4 (Ftt), Francisella tularensis subsp. holarctica OR960246 (Fth), Francisella tularensis subsp. holarctica LVS (LVS), and Francisella novicida U112 (Fn). To gain insight into expression differences that may be related to the differences in virulence of these subspecies, transcriptomes were measured from each strain grown in vitro under identical conditions, utilizing a shared probe microarray. The human avirulent Fn strain exhibited high levels of transcription of genes involved in general metabolism, which are pseudogenes in the human virulent Ftt and Fth strains, consistent with the process of genome decay in the virulent strains. Genes encoding an efflux system (emrA2 cluster of genes), siderophore (fsl operon), acid phosphatase, LPS synthesis, polyamine synthesis, and citrulline ureidase were all highly expressed in Ftt when compared to Fn, suggesting that some of these may contribute to the relative high virulence of Ftt. Genes expressed at a higher level in Ftt when compared to the relatively less virulent Fth included genes encoding isochorismatases, cholylglycine hydrolase, polyamine synthesis, citrulline ureidase, Type IV pilus subunit, and the Francisella Pathogenicity Island protein PdpD. Fth and LVS had very few expression differences, consistent with the derivation of LVS from Fth. This study demonstrated that a shared probe microarray designed to detect transcripts in multiple species/subspecies of Francisella enabled comparative transcriptional analyses that may highlight critical differences that underlie the relative pathogenesis of these strains for humans. This strategy could be extended to other closely-related bacterial species for inter-strain and inter-species analyses.

Tularaemia: clinical aspects in Europe

Tularaemia is a zoonotic disease caused by Francisella tularensis, a Gram-negative, facultative intracellular bacterium. Typically, human and animal infections are caused by F tularensis subspecies tularensis (type A) strains mainly in Canada and USA, and F tularensis subspecies holarctica (type B) strains throughout the northern hemisphere, including Europe. In the past, the epidemiological, clinical, therapeutic, and prognostic aspects of tularaemia reported in the English medical literature were mainly those that had been reported in the USA, where the disease was first described. Tularaemia has markedly changed in the past decade, and a large number of studies have provided novel data for the disease characteristics in Europe. In this Review we aim to emphasise the specific and variable aspects of tularaemia in different European countries. In particular, two natural lifecycles of F tularensis have been described in this continent, although not fully characterised, which are associated with different modes of transmission, clinical features, and public health burdens of tularaemia.

Natural Selection in Virulence Genes of Francisella tularensis

A fundamental tenet of evolution is that alleles that are under negative selection are often deleterious and confer no evolutionary advantage. Negatively selected alleles are removed from the gene pool and are eventually extinguished from the population. Conversely, alleles under positive selection do confer an evolutionary advantage and lead to an increase in the overall fitness of the organism. These alleles increase in frequency until they eventually become fixed in the population. Francisella tularensis is a zoonotic pathogen and a potential biothreat agent. The most virulent type of F. tularensis, Type A, is distributed across North America with Type A.I occurring mainly in the east and Type A.II appearing mainly in the west. F. tularensis is thought to be a genome in decay (losing genes) because of the relatively large number of pseudogenes present in its genome. We hypothesized that the observed frequency of gene loss/pseudogenes may be an artifact of evolution in response to a changing environment, and that genes involved in virulence should be under strong positive selection. To test this hypothesis, we sequenced and compared whole genomes of Type A.I and A.II isolates. We analyzed a subset of virulence and housekeeping genes from several F. tularensis subspecies genomes to ascertain the presence and extent of positive selection. Eleven previously identified virulence genes were screened for positive selection along with 10 housekeeping genes. Analyses of selection yielded one housekeeping gene and 7 virulence genes which showed significant evidence of positive selection at loci implicated in cell surface structures and membrane proteins, metabolism and biosynthesis, transcription, translation and cell separation, and substrate binding and transport. Our results suggest that while the loss of functional genes through disuse could be accelerated by negative selection, the genome decay in Francisella could also be the byproduct of adaptive evolution driven by complex interactions between host, pathogen, and thier environment, as evidenced by several of its virulence genes which are undergoing strong, positive selection.

Twenty years of bacterial genome sequencing

Twenty years ago, the publication of the first bacterial genome sequence, from Haemophilus influenzae, shook the world of bacteriology. In this Timeline, we review the first two decades of bacterial genome sequencing, which have been marked by three revolutions: whole-genome shotgun sequencing, high-throughput sequencing and single-molecule long-read sequencing. We summarize the social history of sequencing and its impact on our understanding of the biology, diversity and evolution of bacteria, while also highlighting spin-offs and translational impact in the clinic. We look forward to a 'sequencing singularity', where sequencing becomes the method of choice for as-yet unthinkable applications in bacteriology and beyond.

Francisella tularensis Subtype A.II Genomic Plasticity in Comparison with Subtype A.I

Although Francisella tularensis is considered a monomorphic intracellular pathogen, molecular genotyping and virulence studies have demonstrated important differences within the tularensis subspecies (type A). To evaluate genetic variation within type A strains, sequencing and assembly of a new subtype A.II genome was achieved for comparison to other completed F. tularensis type A genomes. In contrast with the F. tularensis A.I strains (SCHU S4, FSC198, NE061598, and TI0902), substantial genomic variation was observed between the newly sequenced F. tularensis A.II strain (WY-00W4114) and the only other publically available A.II strain (WY96-3418). Genome differences between WY-00W4114 and WY96-3418 included three major chromosomal translocations, 1580 indels, and 286 nucleotide substitutions of which 159 were observed in predicted open reading frames and 127 were located in intergenic regions. The majority of WY-00W4114 nucleotide deletions occurred in intergenic regions, whereas most of the insertions and substitutions occurred in predicted genes. Of the nucleotide substitutions, 48 (30%) were synonymous and 111 (70%) were nonsynonymous. WY-00W4114 and WY96-3418 nucleotide polymorphisms were predominantly G/C to A/T allelic mutations, with WY-00W4114 having more A+T enrichment. In addition, the A.II genomes contained a considerably higher number of intact genes and longer repetitive sequences, including transposon remnants than the A.I genomes. Together these findings support the premise that F. tularensis A.II may have a fitness advantage compared to the A.I subtype due to the higher abundance of functional genes and repeated chromosomal sequences. A better understanding of the selective forces driving F. tularensis genetic diversity and plasticity is needed.

A combined enrichment and aptamer pulldown assay for Francisella tularensis detection in food and environmental matrices

Francisella tularensis, a Gram-negative bacterium and causative agent of tularemia, is categorized as a Class A select agent by the Centers for Disease Control and Prevention due to its ease of dissemination and ability to cause disease. Oropharyngeal and gastrointestinal tularemia may occur due to ingestion of contaminated food and water. Despite the concern to public health, little research is focused on F. tularensis detection in food and environmental matrices. Current diagnostics rely on host responses and amplification of F. tularensis genetic elements via Polymerase Chain Reaction; however, both tools are limited by development of an antibody response and limit of detection, respectively. During our investigation to develop an improved culture medium to aid F. tularensis diagnostics, we found enhanced F. tularensis growth using the spent culture filtrate. Addition of the spent culture filtrate allowed for increased detection of F. tularensis in mixed cultures of food and environmental matrices. Ultraperformance liquid chromatography (UPLC)/MS analysis identified several unique chemicals within the spent culture supernatant of which carnosine had a matching m/z ratio. Addition of 0.625 mg/mL of carnosine to conventional F. tularensis medium increased the growth of F. tularensis at low inoculums. In order to further enrich F. tularensis cells, we developed a DNA aptamer cocktail to physically separate F. tularensis from other bacteria present in food and environmental matrices. The combined enrichment steps resulted in a detection range of 1-106 CFU/mL (starting inoculums) in both soil and lettuce backgrounds. We propose that the two-step enrichment process may be utilized for easy field diagnostics and subtyping of suspected F. tularensis contamination as well as a tool to aid in basic research of F. tularensis ecology.

Genome sequence and phenotypic analysis of a first German Francisella sp. isolate (W12-1067) not belonging to the species Francisella tularensis

BACKGROUND

Francisella isolates from patients suffering from tularemia in Germany are generally strains of the species F. tularensis subsp. holarctica. To our knowledge, no other Francisella species are known for Germany. Recently, a new Francisella species could be isolated from a water reservoir of a cooling tower in Germany.

RESULTS

We identified a Francisella sp. (isolate W12-1067) whose 16S rDNA is 99% identical to the respective nucleotide sequence of the recently published strain F. guangzhouensis. The overall sequence identity of the fopA, gyrA, rpoA, groEL, sdhA and dnaK genes is only 89%, indicating that strain W12-1067 is not identical to F. guangzhouensis. W12-1067 was isolated from a water reservoir of a cooling tower of a hospital in Germany. The growth optimum of the isolate is approximately 30°C, it can grow in the presence of 4-5% NaCl (halotolerant) and is able to grow without additional cysteine within the medium. The strain was able to replicate within a mouse-derived macrophage-like cell line. The whole genome of the strain was sequenced (~1.7 mbp, 32.2% G + C content) and the draft genome was annotated. Various virulence genes common to the genus Francisella are present, but the Francisella pathogenicity island (FPI) is missing. However, another putative type-VI secretion system is present within the genome of strain W12-1067.

CONCLUSIONS

Isolate W12-1067 is closely related to the recently described F. guangzhouensis species and it replicates within eukaryotic host cells. Since W12-1067 exhibits a putative new type-VI secretion system and F. tularensis subsp. holarctica was found not to be the sole species in Germany, the new isolate is an interesting species to be analyzed in more detail. Further research is needed to investigate the epidemiology, ecology and pathogenicity of Francisella species present in Germany.

Comparative review of Francisella tularensis and Francisella novicida

Francisella tularensis is the causative agent of the acute disease tularemia. Due to its extreme infectivity and ability to cause disease upon inhalation, F. tularensis has been classified as a biothreat agent. Two subspecies of F. tularensis, tularensis and holarctica, are responsible for tularemia in humans. In comparison, the closely related species F. novicida very rarely causes human illness and cases that do occur are associated with patients who are immune compromised or have other underlying health problems. Virulence between F. tularensis and F. novicida also differs in laboratory animals. Despite this varying capacity to cause disease, the two species share ~97% nucleotide identity, with F. novicida commonly used as a laboratory surrogate for F. tularensis. As the F. novicida U112 strain is exempt from U.S. select agent regulations, research studies can be carried out in non-registered laboratories lacking specialized containment facilities required for work with virulent F. tularensis strains. This review is designed to highlight phenotypic (clinical, ecological, virulence, and pathogenic) and genomic differences between F. tularensis and F. novicida that warrant maintaining F. novicida and F. tularensis as separate species. Standardized nomenclature for F. novicida is critical for accurate interpretation of experimental results, limiting clinical confusion between F. novicida and F. tularensis and ensuring treatment efficacy studies utilize virulent F. tularensis strains.

Metabolic network analysis-based identification of antimicrobial drug targets in category A bioterrorism agents

The 2001 anthrax mail attacks in the United States demonstrated the potential threat of bioterrorism, hence driving the need to develop sophisticated treatment and diagnostic protocols to counter biological warfare. Here, by performing flux balance analyses on the fully-annotated metabolic networks of multiple, whole genome-sequenced bacterial strains, we have identified a large number of metabolic enzymes as potential drug targets for each of the three Category A-designated bioterrorism agents including Bacillus anthracis, Francisella tularensis and Yersinia pestis. Nine metabolic enzymes- belonging to the coenzyme A, folate, phosphatidyl-ethanolamine and nucleic acid pathways common to all strains across the three distinct genera were identified as targets. Antimicrobial agents against some of these enzymes are available. Thus, a combination of cross species-specific antibiotics and common antimicrobials against shared targets may represent a useful combinatorial therapeutic approach against all Category A bioterrorism agents.

First indication for a functional CRISPR/Cas system in Francisella tularensis

Francisella tularensis is a zoonotic agent and the subspecies novicida is proposed to be a water-associated bacterium. The intracellular pathogen F. tularensis causes tularemia in humans and is known for its potential to be used as a biological threat. We analyzed the genome sequence of F. tularensis subsp. novicida U112 in silico for the presence of a putative functional CRISPR/Cas (clustered regularly interspaced short palindromic repeats/CRISPR-associated) system. CRISPR/Cas systems are known to encode an RNA-guided adaptive immunity-like system to protect bacteria against invading genetic elements like bacteriophages and plasmids. In this work, we present a first indication that F. tularensis subsp. novicida encodes a functional CRISPR/Cas defence system. Additionally, we identified various spacer DNAs homologous to a putative phage present within the genome of F. tularensis subsp. novicida-like strain 3523. CRISPR/Cas is also present in F. tularensis subsp. tularensis, holarctica, and mediasiatica, but these systems seem to be non-functional.

The class A β-lactamase FTU-1 is native to Francisella tularensis

The class A β-lactamase FTU-1 produces resistance to penicillins and ceftazidime but not to any other β-lactam antibiotics tested. FTU-1 hydrolyzes penicillin antibiotics with catalytic efficiencies of 10(5) to 10(6) M(-1) s(-1) and cephalosporins and carbapenems with catalytic efficiencies of 10(2) to 10(3) M(-1) s(-1), but the monobactam aztreonam and the cephamycin cefoxitin are not substrates for the enzyme. FTU-1 shares 21 to 34% amino acid sequence identity with other class A β-lactamases and harbors two cysteine residues conserved in all class A carbapenemases. FTU-1 is the first weak class A carbapenemase that is native to Francisella tularensis.

Common ancestry and novel genetic traits of Francisella novicida-like isolates from North America and Australia as revealed by comparative genomic analyses

Francisella novicida is a close relative of Francisella tularensis, the causative agent of tularemia. The genomes of F. novicida-like clinical isolates 3523 (Australian strain) and Fx1 (Texas strain) were sequenced and compared to F. novicida strain U112 and F. tularensis strain Schu S4. The strain 3523 chromosome is 1,945,310 bp and contains 1,854 protein-coding genes. The strain Fx1 chromosome is 1,913,619 bp and contains 1,819 protein-coding genes. NUCmer analyses revealed that the genomes of strains Fx1 and U112 are mostly colinear, whereas the genome of strain 3523 has gaps, translocations, and/or inversions compared to genomes of strains Fx1 and U112. Using the genome sequence data and comparative analyses with other members of the genus Francisella, several strain-specific genes that encode putative proteins involved in RTX toxin production, polysaccharide biosynthesis/modification, thiamine biosynthesis, glucuronate utilization, and polyamine biosynthesis were identified. The RTX toxin synthesis and secretion operon of strain 3523 contains four open reading frames (ORFs) and was named rtxCABD. Based on the alignment of conserved sequences upstream of operons involved in thiamine biosynthesis from various bacteria, a putative THI box was identified in strain 3523. The glucuronate catabolism loci of strains 3523 and Fx1 contain a cluster of nine ORFs oriented in the same direction that appear to constitute an operon. Strains U112 and Schu S4 appeared to have lost the loci for RTX toxin production, thiamine biosynthesis, and glucuronate utilization as a consequence of host adaptation and reductive evolution. In conclusion, comparative analyses provided insights into the common ancestry and novel genetic traits of these strains.

Whole-genome sequencing reveals distinct mutational patterns in closely related laboratory and naturally propagated Francisella tularensis strains

The F. tularensis type A strain FSC198 from Slovakia and a second strain FSC043, which has attenuated virulence, are both considered to be derivatives of the North American F. tularensis type A strain SCHU S4. These strains have been propagated under different conditions: the FSC198 has undergone natural propagation in the environment, while the strain FSC043 has been cultivated on artificial media in laboratories. Here, we have compared the genome sequences of FSC198, FSC043, and SCHU S4 to explore the possibility that the contrasting propagation conditions may have resulted in different mutational patterns. We found four insertion/deletion events (INDELs) in the strain FSC043, as compared to the SCHU S4, while no single nucleotide polymorphisms (SNPs) or variable number of tandem repeats (VNTRs) were identified. This result contrasts with previously reported findings for the strain FSC198, where eight SNPs and three VNTR differences, but no INDELs exist as compared to the SCHU S4 strain. The mutations detected in the laboratory and naturally propagated type A strains, respectively, demonstrate distinct patterns supporting that analysis of mutational spectra might be a useful tool to reveal differences in past growth conditions. Such information may be useful to identify leads in a microbial forensic investigation.

Activation of B cell apoptotic pathways in the course of Francisella tularensis infection

Francisella tularensis is a facultative intracellular, gram-negative bacterium that induces apoptosis in macrophages and B cells. Here we show apoptotic pathways that are activated in the Ramos human B cell line in the course of F. tularensis infection. Live bacteria F. tularensis FSC200 activate caspases 8, 9 and 3, as well as Bid; release cytochrome c and apoptosis-inducing factor from mitochondria; and induce depolarization of mitochondrial membrane potential in the Ramos cell line, thus leading these cells to apoptosis. Unlike live bacteria, killed F. tularensis FSC200 bacteria activated only caspase 3, and did not cause apoptosis of Ramos cells as measured by annexin V. Killed bacteria also caused accumulation of anti-apoptotic protein Bclx(L) in mitochondrial membranes. Thus, live F. tularensis activates both caspase pathways (receptor-mediated and intrinsic) as well as caspase-independent mitochondrial death.

Molecular complexity orchestrates modulation of phagosome biogenesis and escape to the cytosol of macrophages by Francisella tularensis

Upon entry of Francisella tularensis to macrophages, the Francisella-containing phagosome (FCP) is trafficked into an acidified late endosome-like phagosome with limited fusion to the lysosomes followed by rapid escape into the cytosol where the organism replicates. Although the Francisella Pathogenicity Island (FPI), which encodes a type VI-like secretion apparatus, is required for modulation of phagosome biogenesis and escape into the cytosol, the mechanisms involved are not known. To decipher the molecular bases of modulation of biogenesis of the FCP and bacterial escape into the macrophage cytosol, we have screened a comprehensive mutant library of F. tularensis ssp. novicida for their defect in proliferation within human macrophages, followed by characterization of modulation of phagosome biogenesis and bacterial escape into the cytosol. Our data show that at least 202 genes are required for intracellular proliferation within macrophages. Among the 125 most defective mutants in intracellular proliferation, we show that the FCP of at least 91 mutants colocalize persistently with the late endosomal/lysosomal marker LAMP-1 and fail to escape into the cytosol, as determined by fluorescence-based phagosome integrity assays and transmission electron microscopy. At least 34 genes are required for proliferation within the cytosol but do not play a detectable role in modulation of phagosome biogenesis and bacterial escape into the cytosol. Our data indicate a tremendous adaptation and metabolic reprogramming by F. tularensis to adjust to the micro-environmental and nutritional cues within the FCP, and these adjustments play essential roles in modulation of phagosome biogenesis and escape into the cytosol of macrophages as well as proliferation in the cytosol. The plethora of the networks of genes that orchestrate F. tularensis-mediated modulation of phagosome biogenesis, phagosomal escape and bacterial proliferation within the cytosol is novel, complex and involves an unusually large portion of the genome of an intracellular pathogen.

High-throughput whole-genome sequencing to dissect the epidemiology of Acinetobacter baumannii isolates from a hospital outbreak

Shared care of military and civilian patients has resulted in transmission of multidrug-resistant Acinetobacter baumannii (MDR-Aci) from military casualties to civilians. Current typing technologies have been useful in revealing relationships between isolates of A. baumannii but they are unable to resolve differences between closely related isolates from small-scale outbreaks, where chains of transmission are often unclear. In a recent hospital outbreak in Birmingham, six patients were colonised with MDR-Aci isolates indistinguishable using standard techniques. We used whole-genome sequencing to identify single nucleotide polymorphisms in these isolates, allowing us to discriminate between alternative epidemiological hypotheses in this setting.

[Mission oriented diagnostic real-time PCR]

In out of area military missions soldiers are potentially exposed to bacteria that are endemic in tropical areas and can be used as biological agents. It can be difficult to culture these bacteria due to sample contamination, low number of bacteria or pretreatment with antibiotics. Commercial biochemical identification systems are not optimized for these agents which can result in misidentification. Immunological assays are often not commercially available or not specific. Real-time PCR assays are very specific and sensitive and can shorten the time required to establish a diagnosis markedly. Therefore, real-time PCRs for the identification of Bacillus anthracis, Brucella spp., Burkholderia mallei und Burkholderia pseudomallei, Francisella tularensis und Yersinia pestis have been developed. PCR results can be false negative due to inadequate clinical samples, low number of bacteria in samples, DNA degradation, inhibitory substances and inappropriate DNA preparation. Hence, it is crucial to cultivate the organisms as a prerequisite for adequate antibiotic therapy and typing of the agent. In a bioterrorist scenario samples have to be treated according to rules applied in forensic medicine and documentation has to be flawless.

Large direct repeats flank genomic rearrangements between a new clinical isolate of Francisella tularensis subsp. tularensis A1 and Schu S4

Francisella tularensis subspecies tularensis consists of two separate populations A1 and A2. This report describes the complete genome sequence of NE061598, an F. tularensis subspecies tularensis A1 isolated in 1998 from a human with clinical disease in Nebraska, United States of America. The genome sequence was compared to Schu S4, an F. tularensis subspecies tularensis A1a strain originally isolated in Ohio in 1941. It was determined that there were 25 nucleotide polymorphisms (22 SNPs and 3 indels) between Schu S4 and NE061598; two of these polymorphisms were in potential virulence loci. Pulsed-field gel electrophoresis analysis demonstrated that NE061598 was an A1a genotype. Other differences included repeat sequences (n = 11 separate loci), four of which were contained in coding sequences, and an inversion and rearrangement probably mediated by insertion sequences and the previously identified direct repeats I, II, and III. Five new variable-number tandem repeats were identified; three of these five were unique in NE061598 compared to Schu S4. Importantly, there was no gene loss or gain identified between NE061598 and Schu S4. Interpretation of these data suggests there is significant sequence conservation and chromosomal synteny within the A1 population. Further studies are needed to determine the biological properties driving the selective pressure that maintains the chromosomal structure of this monomorphic pathogen.

GroEL and lipopolysaccharide from Francisella tularensis live vaccine strain synergistically activate human macrophages

Francisella tularensis, the causative agent of tularemia, interacts with host cells of innate immunity in an atypical manner. For most Gram-negative bacteria, the release of lipopolysaccharide (LPS) from their outer membranes stimulates an inflammatory response. When LPS from the attenuated live vaccine strain (LVS) or the highly virulent Schu S4 strain of F. tularensis was incubated with human umbilical vein endothelial cells, neither species of LPS induced expression of the adhesion molecule E-selectin or secretion of the chemokine CCL2. Moreover, a high concentration (10 microg/ml) of LVS or Schu S4 LPS was required to stimulate production of CCL2 by human monocyte-derived macrophages (huMDM). A screen for alternative proinflammatory factors of F. tularensis LVS identified the heat shock protein GroEL as a potential candidate. Recombinant LVS GroEL at a concentration of 10 microg/ml elicited secretion of CXCL8 and CCL2 by huMDM through a TLR4-dependent mechanism. When 1 microg of LVS GroEL/ml was added to an equivalent amount of LVS LPS, the two components synergistically activated the huMDM to produce CXCL8. Schu S4 GroEL was less stimulatory than LVS GroEL and showed a lesser degree of synergy when combined with Schu S4 LPS. These findings suggest that the intrinsically low proinflammatory activity of F. tularensis LPS may be increased in the infected human host through interactions with other components of the bacterium.

Humans and evolutionary and ecological forces shaped the phylogeography of recently emerged diseases

The development of human civilizations and global commerce has led to the emergence and worldwide circulation of many infectious diseases. Anthrax, plague and tularaemia are three zoonotic diseases that have been intensely studied through genome characterization of the causative species and phylogeographical analyses. A few highly fit genotypes in each species represent the causative agents for most of the observed disease cases. Together, mutational and selective forces create highly adapted pathogens, but this must be coupled with ecological opportunities for global expansion. This Review describes the distributions of the bacteria that cause anthrax, plague and tularaemia and investigates the forces that created clonal structures in these species.

Francisella tularensis type A strains cause the rapid encystment of Acanthamoeba castellanii and survive in amoebal cysts for three weeks postinfection

Francisella tularensis, the causative agent of the zoonotic disease tularemia, has recently gained increased attention due to the emergence of tularemia in geographical areas where the disease has been previously unknown and to the organism's potential as a bioterrorism agent. Although F. tularensis has an extremely broad host range, the bacterial reservoir in nature has not been conclusively identified. In this study, the ability of virulent F. tularensis strains to survive and replicate in the amoeba Acanthamoeba castellanii was explored. We observe that A. castellanii trophozoites rapidly encyst in response to F. tularensis infection and that this rapid encystment phenotype is caused by factor(s) secreted by amoebae and/or F. tularensis into the coculture medium. Further, our results indicate that in contrast to the live vaccine strain LVS, virulent strains of F. tularensis can survive in A. castellanii cysts for at least 3 weeks postinfection and that the induction of rapid amoeba encystment is essential for survival. In addition, our data indicate that pathogenic F. tularensis strains block lysosomal fusion in A. castellanii. Taken together, these data suggest that interactions between F. tularensis strains and amoebae may play a role in the environmental persistence of F. tularensis.

Phylogeography of Francisella tularensis: global expansion of a highly fit clone

Francisella tularensis contains several highly pathogenic subspecies, including Francisella tularensis subsp. holarctica, whose distribution is circumpolar in the northern hemisphere. The phylogeography of these subspecies and their subclades was examined using whole-genome single nucleotide polymorphism (SNP) analysis, high-density microarray SNP genotyping, and real-time-PCR-based canonical SNP (canSNP) assays. Almost 30,000 SNPs were identified among 13 whole genomes for phylogenetic analysis. We selected 1,655 SNPs to genotype 95 isolates on a high-density microarray platform. Finally, 23 clade- and subclade-specific canSNPs were identified and used to genotype 496 isolates to establish global geographic genetic patterns. We confirm previous findings concerning the four subspecies and two Francisella tularensis subsp. tularensis subpopulations and identify additional structure within these groups. We identify 11 subclades within F. tularensis subsp. holarctica, including a new, genetically distinct subclade that appears intermediate between Japanese F. tularensis subsp. holarctica isolates and the common F. tularensis subsp. holarctica isolates associated with the radiation event (the B radiation) wherein this subspecies spread throughout the northern hemisphere. Phylogenetic analyses suggest a North American origin for this B-radiation clade and multiple dispersal events between North America and Eurasia. These findings indicate a complex transmission history for F. tularensis subsp. holarctica.

Inferring genomic flux in bacteria

Acquisition and loss of genetic material are essential forces in bacterial microevolution. They have been repeatedly linked with adaptation of lineages to new lifestyles, and in particular, pathogenicity. Comparative genomics has the potential to elucidate this genetic flux, but there are many methodological challenges involved in inferring evolutionary events from collections of genome sequences. Here we describe a model-based method for using whole-genome sequences to infer the patterns of genome content evolution. A fundamental property of our model is that it allows the rates at which genetic elements are gained or lost to vary in time and from one lineage to another. Our approach is purely sequence based, and does not rely on gene identification. We show how inference can be performed under our model and illustrate its use on three datasets from Francisella tularensis, Streptococcus pyogenes, and Escherichia coli. In all three examples, we found interesting variations in the rates of genetic material gain and loss, which strongly correlate with their lifestyle. The algorithms we describe are implemented in a computer software named GenoPlast.

Bacterial pathogenomics

Genomes from all of the crucial bacterial pathogens of humans, plants and animals have now been sequenced, as have genomes from many of the important commensal, symbiotic and environmental microorganisms. Analysis of these sequences has revealed the forces that shape pathogen evolution and has brought to light unexpected aspects of pathogen biology. The finding that horizontal gene transfer and genome decay have key roles in the evolution of bacterial pathogens was particularly surprising. It has also become evident that even the definitions for 'pathogen' and 'virulence factor' need to be re-evaluated.