The phylogeographic pattern of Francisella tularensis in Sweden indicates a Scandinavian origin of Eurosiberian tularaemia.

Previous studies of the causative agent of tularaemia, Francisella tularensis have identified phylogeographic patterns suggestive of environmental maintenance reservoirs. To investigate the phylogeography of tularaemia in Sweden, we selected 163 clinical isolates obtained during 1995-2009 in 10 counties and sequenced one isolate's genome to identify new genetic markers. An improved typing scheme based on two indels and nine SNPs was developed using hydrolysis or TaqMan MGB probe assays. The results showed that much of the known global genetic diversity of F. tularensis subsp. holarctica is present in Sweden. Thirteen of the 163 isolates belonged to a new genetic group that is basal to all other known members of the major genetic clade B.I, which is spread across the Eurosiberian region. One hundred and twenty-five of the 163 Swedish isolates belonged to B.I, but individual clades' frequencies differed from county to county (P < 0.001). Subsequent analyses revealed a correlation between genotype variation over time and recurrent outbreaks at specific places, supporting the 'maintenance reservoir' environmental maintenance hypothesis. Most importantly, the findings reveal the presence of diverse source populations of F. tularensis subsp. holarctica in Sweden and suggest a historical spread of the disease from Scandinavia to other parts of Eurosiberia.

Genome sequence of Francisella tularensis subspecies holarctica strain FSC200, isolated from a child with tularemia.

Here we report the complete, accurate 1.89-Mb genome sequence of Francisella tularensis subsp. holarctica strain FSC200, isolated in 1998 in the Swedish municipality Ljusdal, which is in an area where tularemia is highly endemic. This genome is important because strain FSC200 has been extensively used for functional and genetic studies of Francisella and is well-characterized.

Increased knowledge of Francisella genus diversity highlights the benefits of optimised DNA-based assays.

BACKGROUND: Recent advances in sequencing technologies offer promising tools for generating large numbers of genomes, larger typing databases and improved mapping of environmental bacterial diversity. However, DNA-based methods for the detection of Francisella were developed with limited knowledge about genetic diversity. This, together with the high sequence identity between several Francisella species, means there is a high risk of false identification and detection of the highly virulent pathogen Francisella tularensis. Moreover, phylogenetic reconstructions using single or limited numbers of marker sequences often result in incorrect tree topologies and inferred evolutionary distances. The recent growth in publicly accessible whole-genome sequences now allows evaluation of published genetic markers to determine optimal combinations of markers that minimise both time and laboratory costs. RESULTS: In the present study, we evaluated 38 previously published DNA markers and the corresponding PCR primers against 42 genomes representing the currently known diversity of the genus Francisella. The results highlight that PCR assays for Francisella tularensis are often complicated by low specificity, resulting in a high probability of false positives. A method to select a set of one to seven markers for obtaining optimal phylogenetic resolution or diagnostic accuracy is presented. CONCLUSIONS: Current multiple-locus sequence-typing systems and detection assays of Francisella, could be improved by redesigning some of the primers and reselecting typing markers. The use of only a few optimally selected sequence-typing markers allows construction of phylogenetic topologies with almost the same accuracy as topologies based on whole-genome sequences.

Genome characterisation of the genus Francisella reveals insight into similar evolutionary paths in pathogens of mammals and fish.

BACKGROUND: Prior to this study, relatively few strains of Francisella had been genome-sequenced. Previously published Francisella genome sequences were largely restricted to the zoonotic agent F. tularensis. Only limited data were available for other members of the Francisella genus, including F. philomiragia, an opportunistic pathogen of humans, F. noatunensis, a serious pathogen of farmed fish, and other less well described endosymbiotic species. RESULTS: We determined the phylogenetic relationships of all known Francisella species, including some for which the phylogenetic positions were previously uncertain. The genus Francisella could be divided into two main genetic clades: one included F. tularensis, F. novicida, F. hispaniensis and Wolbachia persica, and another included F. philomiragia and F. noatunensis.Some Francisella species were found to have significant recombination frequencies. However, the fish pathogen F. noatunensis subsp. noatunensis was an exception due to it exhibiting a highly clonal population structure similar to the human pathogen F. tularensis. CONCLUSIONS: The genus Francisella can be divided into two main genetic clades occupying both terrestrial and marine habitats. However, our analyses suggest that the ancestral Francisella species originated in a marine habitat. The observed genome to genome variation in gene content and IS elements of different species supports the view that similar evolutionary paths of host adaptation developed independently in F. tularensis (infecting mammals) and F. noatunensis subsp. noatunensis (infecting fish).

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.

Landscape epidemiology of tularemia outbreaks in Sweden.

Summer outbreaks of tularemia that occurred from 1995 through 2005 in 2 locations in Sweden affected 441 persons. We performed an epidemiologic investigation of these outbreaks using a novel strategy, involving high-resolution genotyping of Francisella tularensis isolates obtained from 136 patients (using 18 genetic markers developed from 6 F. tularensis genome sequences) and interviews with the patients. Strong spatial associations were found between F. tularensis subpopulations and the places of disease transmission; infection by some subpopulations occurred within areas as small as 2 km(2), indicating unidentified environmental point sources of tularemia. In both locations, disease clusters were associated with recreational areas beside water, and genetic subpopulations were present throughout the tularemia season and persisted over years. High-resolution genotyping in combination with patients' statements about geographic places of disease transmission provided valuable indications of likely sources of infection and the causal genotypes during these tularemia outbreaks.

A real-time PCR array for hierarchical identification of Francisella isolates.

A robust, rapid and flexible real-time PCR assay for hierarchical genetic typing of clinical and environmental isolates of Francisella is presented. Typing markers were found by multiple genome and gene comparisons, from which 23 canonical single nucleotide polymorphisms (canSNPs) and 11 canonical insertion-deletion mutations (canINDELs) were selected to provide phylogenetic guidelines for classification from genus to isolate level. The specificity of the developed assay, which uses 68 wells of a 96-well real-time PCR format with a detection limit of 100 pg DNA, was assessed using 62 Francisella isolates of diverse genetic and geographical origins. It was then successfully used for typing 14 F. tularensis subsp. holarctica isolates obtained from tularemia patients in Sweden in 2008 and five more genetically diverse Francisella isolates of global origins. When applied to human ulcer specimens for direct pathogen detection the results were incomplete due to scarcity of DNA, but sufficient markers were identified to detect fine-resolution differences among F. tularensis subsp. holarctica isolates causing infection in the patients. In contrast to other real-time PCR assays for Francisella, which are typically designed for specific detection of a species, subspecies, or strain, this type of assay can be easily tailored to provide appropriate phylogenetic and/or geographical resolution to meet the objectives of the analysis.

Molecular evolutionary consequences of niche restriction in Francisella tularensis, a facultative intracellular pathogen.

Francisella tularensis is a potent mammalian pathogen well adapted to intracellular habitats, whereas F. novicida and F. philomiragia are less virulent in mammals and appear to have less specialized lifecycles. We explored adaptations within the genus that may be linked to increased host association, as follows. First, we determined the genome sequence of F. tularensis subsp. mediasiatica, the only subspecies that had not been previously sequenced. This genome, and those of 12 other F. tularensis isolates, were then compared to the genomes of F. novicida (three isolates) and F. philomiragia (one isolate). Signs of homologous recombination were found in approximately 19.2% of F. novicida and F. philomiragia genes, but none among F. tularensis genomes. In addition, random insertions of insertion sequence elements appear to have provided raw materials for secondary adaptive mutations in F. tularensis, e.g. for duplication of the Francisella Pathogenicity Island and multiplication of a putative glycosyl transferase gene. Further, the five major genetic branches of F. tularensis seem to have converged along independent routes towards a common gene set via independent losses of gene functions. Our observations suggest that despite an average nucleotide identity of >97%, F. tularensis and F. novicida have evolved as two distinct population lineages, the former characterized by clonal structure with weak purifying selection, the latter by more frequent recombination and strong purifying selection. F. tularensis and F. novicida could be considered the same bacterial species, given their high similarity, but based on the evolutionary analyses described in this work we propose retaining separate species names.

A 55 kDa hypothetical membrane protein is an iron-regulated virulence factor of Francisella tularensis subsp. novicida U112.

Iron is an important nutritional requirement for bacteria due to its conserved role in many essential metabolic processes. As a consequence of the lack of freely available iron in the mammalian host, bacteria upregulate a range of virulence factors during infection. Transcriptional analysis of Francisella tularensis subsp. novicida U112 grown in iron-deficient medium identified 21 genes upregulated in response to this condition, four of which were attributed to a siderophore operon. In addition, a novel iron-regulated gene, FTT0025, was identified which is part of this operon and encodes a 55 kDa hypothetical membrane protein. When grown on chrome azurol S agar, the F. tularensis subsp. novicida U112deltaFTT0025 mutant produced an increased reaction zone compared with the wild-type, suggesting that siderophore production was unaffected but that the bacteria may have a deficiency in their ability to re-sequester this iron-binding molecule. Furthermore, the deltaFTT0025 mutant was attenuated in a BALB/c mouse model of infection relative to wild-type F. tularensis subsp. novicida U112.

Comparison of Francisella tularensis genomes reveals evolutionary events associated with the emergence of human pathogenic strains.

BACKGROUND: Francisella tularensis subspecies tularensis and holarctica are pathogenic to humans, whereas the two other subspecies, novicida and mediasiatica, rarely cause disease. To uncover the factors that allow subspecies tularensis and holarctica to be pathogenic to humans, we compared their genome sequences with the genome sequence of Francisella tularensis subspecies novicida U112, which is nonpathogenic to humans. RESULTS: Comparison of the genomes of human pathogenic Francisella strains with the genome of U112 identifies genes specific to the human pathogenic strains and reveals pseudogenes that previously were unidentified. In addition, this analysis provides a coarse chronology of the evolutionary events that took place during the emergence of the human pathogenic strains. Genomic rearrangements at the level of insertion sequences (IS elements), point mutations, and small indels took place in the human pathogenic strains during and after differentiation from the nonpathogenic strain, resulting in gene inactivation. CONCLUSION: The chronology of events suggests a substantial role for genetic drift in the formation of pseudogenes in Francisella genomes. Mutations that occurred early in the evolution, however, might have been fixed in the population either because of evolutionary bottlenecks or because they were pathoadaptive (beneficial in the context of infection). Because the structure of Francisella genomes is similar to that of the genomes of other emerging or highly pathogenic bacteria, this evolutionary scenario may be shared by pathogens from other species.

Canonical insertion-deletion markers for rapid DNA typing of Francisella tularensis.

To develop effective and accurate typing of strains of Francisella tularensis, a potent human pathogen and a putative bioterrorist agent, we combined analysis of insertion-deletion (indel) markers with multiple-locus variable-number tandem repeat analysis (MLVA). From 5 representative F. tularensis genome sequences, 38 indel markers with canonical properties, i.e., capable of sorting strains into major genetic groups, were selected. To avoid markers with a propensity for homoplasy, we used only those indels with 2 allelic variants and devoid of substantial sequence repeats. MLVA included sequences with much diversity in copy number of tandem repeats. The combined procedure allowed subspecies division, delineation of clades A.I and A.II of subspecies tularensis, differentiation of Japanese strains from other strains of subspecies holarctica, and high-resolution strain typing. The procedure uses limited amounts of killed bacterial preparations and, because only 1 single analytic method is needed, is time- and cost-effective.

Potential source of Francisella tularensis live vaccine strain attenuation determined by genome comparison.

Francisella tularensis is a bacterial pathogen that causes the zoonotic disease tularemia and is important to biodefense. Currently, the only vaccine known to confer protection against tularemia is a specific live vaccine strain (designated LVS) derived from a virulent isolate of Francisella tularensis subsp. holarctica. The origin and source of attenuation of this strain are not known. To assist with the design of a defined live vaccine strain, we sought to determine the genetic basis of the attenuation of LVS. This analysis relied primarily on the comparison between the genome of LVS and Francisella tularensis holarctica strain FSC200, which differ by only 0.08% of their nucleotide sequences. Under the assumption that the attenuation was due to a loss of function(s), only coding regions were examined in this comparison. To complement this analysis, the coding regions of two slightly more distantly related Francisella tularensis strains were also compared against the LVS coding regions. Thirty-five genes show unique sequence variations predicted to alter the protein sequence in LVS compared to the other Francisella tularensis strains. Due to these polymorphisms, the functions of 15 of these genes are very likely lost or impaired. Seven of these genes were demonstrated to be under stronger selective constraints, suggesting that they are the most probable to be the source of LVS attenuation and useful for a newly defined vaccine.

Direct repeat-mediated deletion of a type IV pilin gene results in major virulence attenuation of Francisella tularensis.

Francisella tularensis, the causative agent of tularaemia, is a highly infectious and virulent intracellular pathogen. There are two main human pathogenic subspecies, Francisella tularensis ssp. tularensis (type A), and Francisella tularensis ssp. holarctica (type B). So far, knowledge regarding key virulence determinants is limited but it is clear that intracellular survival and multiplication is one major virulence strategy of Francisella. In addition, genome sequencing has revealed the presence of genes encoding type IV pili (Tfp). One genomic region encoding three proteins with signatures typical for type IV pilins contained two 120 bp direct repeats. Here we establish that repeat-mediated loss of one of the putative pilin genes in a type B strain results in severe virulence attenuation in mice infected by subcutaneous route. Complementation of the mutant by introduction of the pilin gene in cis resulted in complete restoration of virulence. The level of attenuation was similar to that of the live vaccine strain and this strain was also found to lack the pilin gene as result of a similar deletion event mediated by the direct repeats. Presence of the pilin had no major effect on the ability to interact, survive and multiply inside macrophage-like cell lines. Importantly, the pilin-negative strain was impaired in its ability to spread from the initial site of infection to the spleen. Our findings indicate that this putative pilin is critical for Francisella infections that occur via peripheral routes.

A mutant of Francisella tularensis strain SCHU S4 lacking the ability to express a 58-kilodalton protein is attenuated for virulence and is an effective live vaccine.

Francisella tularensis subsp. tularensis (type A) strain SCHU S4 is a prototypic strain of the pathogen that is highly virulent for humans and other mammals. Its intradermal (i.d.) 50% lethal dose (LD50) for mice is <10 CFU. We discovered a spontaneous mutant, designated FSC043, of SCHU S4 with an i.d. LD50 of >10(8) CFU. FSC043 effectively vaccinated mice against challenge with a highly virulent type A strain, and the protective efficacy was at least as good as that of F. tularensis LVS, an empirically attenuated strain which has been used as an efficacious human vaccine. Comparative proteomics was used to identify two proteins of unknown function that were identified as defective in LVS and FSC043, and deletion mutants of SCHU S4 were created for each of the two encoding genes. One mutant, the DeltaFTT0918 strain, failed to express a 58-kDa protein, had an i.d. LD50 of approximately 10(5) CFU, and was found to be less capable than SCHU S4 of growing in peritoneal mouse macrophages. Mice that recovered from sublethal infection with the DeltaFTT0918 mutant survived when challenged 2 months later with >100 LD50s of the highly virulent type A strain FSC033. This is the first report of the generation of defined mutants of F. tularensis subsp. tularensis and their use as live vaccines.

Evolution of subspecies of Francisella tularensis.

Analysis of unidirectional genomic deletion events and single nucleotide variations suggested that the four subspecies of Francisella tularensis have evolved by vertical descent. The analysis indicated an evolutionary scenario where the highly virulent F. tularensis subsp. tularensis (type A) appeared before the less virulent F. tularensis subsp. holarctica (type B). Compared to their virulent progenitors, attenuated strains of F. tularensis exhibited specific unidirectional gene losses.

The complete genome sequence of Francisella tularensis, the causative agent of tularemia.

Francisella tularensis is one of the most infectious human pathogens known. In the past, both the former Soviet Union and the US had programs to develop weapons containing the bacterium. We report the complete genome sequence of a highly virulent isolate of F. tularensis (1,892,819 bp). The sequence uncovers previously uncharacterized genes encoding type IV pili, a surface polysaccharide and iron-acquisition systems. Several virulence-associated genes were located in a putative pathogenicity island, which was duplicated in the genome. More than 10% of the putative coding sequences contained insertion-deletion or substitution mutations and seemed to be deteriorating. The genome is rich in IS elements, including IS630 Tc-1 mariner family transposons, which are not expected in a prokaryote. We used a computational method for predicting metabolic pathways and found an unexpectedly high proportion of disrupted pathways, explaining the fastidious nutritional requirements of the bacterium. The loss of biosynthetic pathways indicates that F. tularensis is an obligate host-dependent bacterium in its natural life cycle. Our results have implications for our understanding of how highly virulent human pathogens evolve and will expedite strategies to combat them.

The crystal structure of the carbohydrate-recognition domain of the glycoprotein sorting receptor p58/ERGIC-53 reveals an unpredicted metal-binding site and conformational changes associated with calcium ion binding.

p58/ERGIC-53 is a calcium-dependent animal lectin that acts as a cargo receptor, binding to a set of glycoproteins in the endoplasmic reticulum (ER) and transporting them to the Golgi complex. It is similar in structure to calcium-dependent leguminous lectins. We have determined the structure of the carbohydrate-recognition domain of p58/ERGIC-53 in its calcium-bound form. The structure reveals localized but large conformational changes in relation to the previously determined metal ion-free structure, mapping mostly to the ligand-binding site. It reveals the presence of two calcium ion-binding sites located 6A apart, one of which has no equivalent in the plant lectins. The second metal ion-binding site present in that class of lectins, binding Mn(2+), is absent from p58/ERGIC-53. The absence of a short loop in the ligand-binding site in this protein suggests that it has adapted to optimally bind the high-mannose Man(8)(GlcNAc)(2) glycan common to glycoproteins at the ER exit stage.

VIPL, a VIP36-like membrane protein with a putative function in the export of glycoproteins from the endoplasmic reticulum.

Subsets of glycoproteins are thought to require lectin-like membrane receptors for efficient export out of the endoplasmic reticulum (ER). To identify new members related to two previously characterized intracellular lectins ERGIC-53/p58 and VIP36, we carried out an extensive database search using the conserved carbohydrate recognition domain (CRD) as a search string. A gene, more closely related to VIP36 than to ERGIC-53/p58, and hence called VIPL (VIP36-Like), was identified. VIPL has been conserved through evolution from zebra fish to man. The 2.4-kb VIPL mRNA was widely expressed to varying levels in different tissues. Using an antiserum prepared against the CRD, the 32-kDa VIPL protein was detected in various cell lines. The single N-linked glycan of VIPL remained endoglycosidase H-sensitive during a 2-h pulse-chase, even when the protein was overexpressed or mutated to allow export to the plasma membrane. VIPL localized primarily to the ER and partly to the Golgi complex. Like VIP36, the cytoplasmic tail of VIPL terminates in the sequence KRFY, a motif characteristic for proteins recycling between the ER and ERGIC/cis-Golgi. Mutating the retrograde transport signal KR to AA resulted in transport of VIPL to the cell surface. Finally, knock-down of VIPL mRNA using siRNA significantly slowed down the secretion of two glycoproteins (M(R) 35 and 250 kDa) to the medium, suggesting that VIPL may also function as an ER export receptor.

Expression, purification, refolding and crystallization of the carbohydrate-recognition domain of p58/ERGIC-53, an animal C-type lectin involved in export of glycoproteins from the endoplasmic reticulum.

p58/ERGIC-53 is a mammalian calcium-dependent lectin that serves as a glycoprotein-sorting receptor between the endoplasmic reticulum (ER) and the Golgi complex. It is a type I transmembrane protein with two lumenal domains, one of which is a carbohydrate-recognition domain (CRD) and homologous to leguminous lectins. The CRD of p58, the rat homologue of human ERGIC-53, was overexpressed in insect cells and Escherichia coli, purified and crystallized using Li(2)SO(4) as a precipitant. The crystals belong to space group I222, with unit-cell parameters a = 49.6, b = 86.1, c = 128.1 A, and contain one molecule per asymmetric unit, corresponding to a packing density of 2.4 A(3)Da(-1). Knowledge of the structure of p58/ERGIC-53 will provide a starting model for understanding receptor-mediated glycoprotein sorting between the ER and the Golgi.

Crystal structure of the carbohydrate recognition domain of p58/ERGIC-53, a protein involved in glycoprotein export from the endoplasmic reticulum.

p58/ERGIC-53 is an animal calcium-dependent lectin that cycles between the endoplasmic reticulum (ER) and the Golgi complex and appears to act as a cargo receptor for a subset of soluble glycoproteins exported from the ER. We have determined the crystal structure of the carbohydrate recognition domain (CRD) of p58, the rat homologue of human ERGIC-53, to 1.46 A resolution. The fold and ligand binding site are most similar to those of leguminous lectins. The structure also resembles that of the CRD of the ER folding chaperone calnexin and the neurexins, a family of non-lectin proteins expressed on neurons. The CRD comprises one concave and one convex beta-sheet packed into a beta-sandwich. The ligand binding site resides in a negatively charged cleft formed by conserved residues. A large surface patch of conserved residues with a putative role in protein-protein interactions and oligomerization lies on the opposite side of the ligand binding site. Together with previous functional data, the structure defines a new and expanding class of calcium-dependent animal lectins and provides a starting point for the understanding of glycoprotein sorting between the ER and the Golgi.