Comparative genomic characterization of Francisella tularensis strains belonging to low and high virulence subspecies

Avirulence of a spontaneous Francisella tularensis subsp. mediasiatica prmA mutant

Francisella tularensis, the causative agent of tularemia, is divided into three subspecies. Two of these, subspecies holarctica and tularensis, are highly pathogenic to humans and consequently relatively well studied. The third subspecies, mediasiatica, is rarely isolated and remains poorly studied. It is distributed in the sparsely populated regions of Central Asia and Siberia. Curently this subspecies is not known to have been responsible for human infections in spite of its high virulence in laboratory animals. Subspecies mediasiatica is currently divided into three subgroups-MI, present in Central Asia, MII, present in southern Siberia, and MIII represented by a unique strain, 60(B)57, isolated in Uzbekistan in 1960. We describe here the unexpected observation that MIII strain 60(B)57 is avirulent and immunogenic. We observed that infection with this strain protected mice from challenge 21 days later with a virulent subsp. mediasiatica strain. With an increase of this interval, the protection for mice was significantly reduced. In contrast, guinea pigs were protected from challenge with strains of the subspecies holarctica and mediasiatica (but not subsp. tularensis) 90 days after infection with 60(B)57. We performed genome assembly based on whole genome sequencing data obtained using the Nanopore MinION for strain 60(B)57 and two subsp. mediasiatica strains representing the Central Asian MI and Siberian MII phylogenetic subgroups. The prmA gene is truncated due to a nonsense mutation in strain 60(B)57. The deletion of gene prmA has previously been shown to induce a loss of virulence in Francisella novicida the closest model organism suggesting that the observed mutation might the cause of the avirulence of strain 60(B)57.

Tularemia - a re-emerging disease with growing concern

Tularemia caused by Gram-negative, coccobacillus bacterium, Francisella tularensis, is a highly infectious zoonotic disease. Human cases have been reported mainly from the United States, Nordic countries like Sweden and Finland, and some European and Asian countries. Naturally, the disease occurs in several vertebrates, particularly lagomorphs. Type A (subspecies tularensis) is more virulent and causes disease mainly in North America; type B (subspecies holarctica) is widespread, while subspecies mediasiatica is present in central Asia. F. tularensis is a possible bioweapon due to its lethality, low infectious dosage, and aerosol transmission. Small mammals like rabbits, hares, and muskrats are primary sources of human infections, but true reservoir of F. tularensis is unknown. Vector-borne tularemia primarily involves ticks and mosquitoes. The bacterial subspecies involved and mode of transmission determine the clinical picture. Early signs are flu-like illnesses that may evolve into different clinical forms of tularemia that may or may not include lymphadenopathy. Ulcero-glandular and glandular forms are acquired by arthropod bite or handling of infected animals, oculo-glandular form as a result of conjunctival infection, and oro-pharyngeal form by intake of contaminated food or water. Pulmonary form appears after inhalation of bacteria. Typhoidal form may occur after infection via different routes. Human-to-human transmission has not been known. Diagnosis can be achieved by serology, bacterial culture, and molecular methods. Treatment for tularemia typically entails use of quinolones, tetracyclines, or aminoglycosides. Preventive measures are necessary to avoid infection although difficult to implement. Research is underway for the development of effective live attenuated and subunit vaccines.

Phenotypic and genotypic discrimination of Francisella tularensis ssp. holarctica clades

Francisella tularensis is the causative agent of tularemia, a zoonotic disease with a wide host range. F. tularensis ssp. holarctica (Fth) is of clinical relevance for European countries, including Germany. Whole genome sequencing methods, including canonical Single Nucleotide Polymorphism (canSNP) typing and whole genome SNP typing, have revealed that European Fth strains belong to a few monophyletic populations. The majority of German Fth isolates belong to two basal phylogenetic clades B.6 (biovar I) and B.12 (biovar II). Strains of B.6 and B.12 seem to differ in their pathogenicity, and it has been shown that strains of biovar II are resistant against erythromycin. In this study, we present data corroborating our previous data demonstrating that basal clade B.12 can be divided into clades B.71 and B.72. By applying phylogenetic whole genome analysis as well as proteome analysis, we could verify that strains of these two clades are distinct from one another. This was confirmed by measuring the intensity of backscatter light on bacteria grown in liquid media. Strains belonging to clades B.6, B.71 or B.72 showed clade-specific backscatter growth curves. Furthermore, we present the whole genome sequence of strain A-1341, as a reference genome of clade B.71, and whole proteomes comparison of Fth strains belonging to clades B.6, B.71 and B.72. Further research is necessary to investigate phenotypes and putative differences in pathogenicity of the investigated different clades of Fth to better understand the relationship between observed phenotypes, pathogenicity and distribution of Fth strains.

Identification of Francisella tularensis Subspecies in a Clinical Setting Using MALDI-TOF MS: An In-House Francisella Library and Biomarkers

Francisella tularensis is a zoonotic bacterium that is endemic in large parts of the world. It is absent in the standard library of the most applied matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) systems: the Vitek MS and the Bruker Biotyper system. The additional Bruker MALDI Biotyper Security library contains F. tularensis without subspecies differentiation. The virulence of F. tularensis differs between the subspecies. The F. tularensis subspecies (ssp.) tularensis is highly pathogenic, whereas the subspecies holarctica displays lower virulence and subspecies novicida and F. tularensis ssp. mediasiatica are hardly virulent. To differentiate the Francisellaceae and the F. tularensis-subspecies, an in-house Francisella library was built with the Bruker Biotyper system and validated together with the existing Bruker databases. In addition, specific biomarkers were defined based on the main spectra of the Francisella strains supplemented with in silico genome data. Our in-house Francisella library accurately differentiates the F. tularensis subspecies and the other Francisellaceae. The biomarkers correctly differentiate the various species within the genus Francisella and the F. tularensis subspecies. These MALDI-TOF MS strategies can successfully be applied in a clinical laboratory setting as a fast and specific method to identify F. tularensis to subspecies level.

Heteroresistance to Amikacin in Carbapenem-Resistant Klebsiella pneumoniae Strains

Heteroresistance can lead to treatment failure and is difficult to detect by the methods currently employed by clinical laboratories. The aim of this study was to investigate the prevalence of the amikacin-heteroresistant Klebsiella pneumoniae strains and explore potential amikacin heteroresistance mechanism through whole-genome sequencing (WGS) and quantitative reverse-transcription PCR (qRT-PCR). In this study, 13 isolates (8.39%) were considered as amikacin-heteroresistant K. pneumoniae strains among a total of 155 K. pneumoniae strains. The majority of the heterogeneous phenotypes (11/13, 84.61%) was unstable and the minimal inhibitory concentrations (MICs) fully or partially reverted back to the level of susceptibility of the parental isolate. The frequency of heteroresistant subpopulation ranged from 2.94×10^-7 to 5.59×10^-6. Whole-genome sequencing and single-nucleotide variants (SNVs) analysis showed that there were different nucleotide and resultant amino acid alterations among an amikacin-heteroresistant strain S38 and the resistant subpopulation S38L in several genes. Quantitative reverse-transcription PCR analysis revealed that the increased expression of aminoglycoside resistance genes detected in amikacin-heteroresistant K. pneumoniae strains might be associated with amikacin heteroresistance. The findings raise concerns for the emergence of amikacin-heteroresistant K. pneumoniae strains and the use of amikacin as therapy for the treatment of multidrug-resistant K. pneumoniae strains.

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.

Structural and functional analysis of the Francisella lysine decarboxylase as a key actor in oxidative stress resistance

Francisella tularensis is one of the most virulent pathogenic bacteria causing the acute human respiratory disease tularemia. While the mechanisms underlying F. tularensis pathogenesis are largely unknown, previous studies have shown that a F. novicida transposon mutant with insertions in a gene coding for a putative lysine decarboxylase was attenuated in mouse spleen, suggesting a possible role of its protein product as a virulence factor. Therefore, we set out to structurally and functionally characterize the F. novicida lysine decarboxylase, which we termed LdcF. Here, we investigate the genetic environment of ldcF as well as its evolutionary relationships with other basic AAT-fold amino acid decarboxylase superfamily members, known as key actors in bacterial adaptative stress response and polyamine biosynthesis. We determine the crystal structure of LdcF and compare it with the most thoroughly studied lysine decarboxylase, E. coli LdcI. We analyze the influence of ldcF deletion on bacterial growth under different stress conditions in dedicated growth media, as well as in infected macrophages, and demonstrate its involvement in oxidative stress resistance. Finally, our mass spectrometry-based quantitative proteomic analysis enables identification of 80 proteins with expression levels significantly affected by ldcF deletion, including several DNA repair proteins potentially involved in the diminished capacity of the F. novicida mutant to deal with oxidative stress. Taken together, we uncover an important role of LdcF in F. novicida survival in host cells through participation in oxidative stress response, thereby singling out this previously uncharacterized protein as a potential drug target.

Reorganized Genomic Taxonomy of Francisellaceae Enables Design of Robust Environmental PCR Assays for Detection of Francisella tularensis

In recent years, an increasing diversity of species has been recognized within the family Francisellaceae. Unfortunately, novel isolates are sometimes misnamed in initial publications or multiple sources propose different nomenclature for genetically highly similar isolates. Thus, unstructured and occasionally incorrect information can lead to confusion in the scientific community. Historically, detection of Francisella tularensis in environmental samples has been challenging due to the considerable and unknown genetic diversity within the family, which can result in false positive results. We have assembled a comprehensive collection of genome sequences representing most known Francisellaceae species/strains and restructured them according to a taxonomy that is based on phylogenetic structure. From this structured dataset, we identified a small number of genomic regions unique to F. tularensis that are putatively suitable for specific detection of this pathogen in environmental samples. We designed and validated specific PCR assays based on these genetic regions that can be used for the detection of F. tularensis in environmental samples, such as water and air filters.

Tularemia as a Mosquito-Borne Disease

Francisella tularensis (Ft) is the etiological agent of tularemia, a disease known for over 100 years in the northern hemisphere. Ft includes four subspecies, of which two are the etiologic agents of tularemia: Ft subsp. tularensis (Ftt) and Ft subsp. holarctica (Fth), mainly distributed in North America and the whole northern hemisphere, respectively. Several routes of human infection with these bacteria exist, notably through bites of Ixodidae ticks. However, mosquitoes represent the main vectors of Fth in Scandinavia, where large tularemia outbreaks have occurred, usually during the warm season. The mechanisms making mosquitoes vectors of Fth are still unclear. This review covers the inventory of research work and epidemiological data linking tularemia to mosquitoes in Scandinavia and highlights the gaps in understanding mosquitoes and Ft interactions.

Differences in Blood-derived Francisella tularensis Type B Strains from Clinical Cases of Tularemia

Francisella tularensis can cause the zoonotic disease tularemia and is partitioned into subspecies due to differences in chromosomal organization and virulence. The subspecies holarctica (type B) is generally considered more clonal than the other subpopulations with moderate virulence compared to the hypervirulent A.I clade. We performed whole genome sequencing (WGS) on six type B strains isolated from the blood of patients with tularemia within a one-year period from the same United States region, to better understand the associated pathogenicity. The WGS data were compared to the prototype strain for this subspecies, specifically FSC200, which was isolated from a patient with tularemia in Europe. These findings revealed 520–528 single nucleotide polymorphisms (SNPs) between the six United States type B strains compared to FSC200, with slightly higher A+T content in the latter strain. In contrast, comparisons between the six type B isolates showed that five of the six type B isolates had only 4–22 SNPs, while one of the strains had 47–53 SNPs. Analysis of SNPs in the core genome for the six United States type B isolates and the FSC200 strain gave similar results, suggesting that some of these mutations may have been nonsynonymous, resulting in altered protein function and pathogenicity.

The Comparative Virulence of Francisella tularensis Subsp. mediasiatica for Vaccinated Laboratory Animals

Tularemia is a severe infectious disease caused by the Gram-negative bacteria Fracisella tularensis. There are four subspecies of F.tularensis: holarctica, tularensis, mediasiatica, and novicida, which differ in their virulence and geographic distribution. One of them, subsp. mediasiatica remains extremely poorly studied, primarily due to the fact that it is found only in the sparsely populated regions of Central Asia and Russia. In particular there is little information in the literature on the virulence and pathogenicity of subsp. mediasiatica. In the present article, we evaluated the comparative virulence of subsp. mediasiatica in vaccinated laboratory animals which we infected with virulent strains: subsp. mediasiatica 678, subsp. holarctica 503, and subsp. tularensis SCHU within 60 to 180 days after vaccination. We found that subsp. mediasiatica is comparable in pathogenicity in mice with subsp. tularensis and in guinea pigs with subsp. holarctica. We also found that the live vaccine does not fully protect mice from subsp. mediasiatica but completely protects guinea pigs for at least six months. In general, our data suggest that subsp. mediasiatica occupies an intermediate position in virulence between spp. tularensis and holarctica.

Francisella novicida CRISPR-Cas Systems Can Functionally Complement Each Other in DNA Defense while Providing Target Flexibility

CRISPR-Cas systems are prokaryotic adaptive immune systems that facilitate protection of bacteria and archaea against infection by external mobile genetic elements. The model pathogen Francisella novicida encodes a CRISPR-Cas12a (FnoCas12a) system and a CRISPR-Cas9 (FnoCas9) system, the latter of which has an additional and noncanonical function in bacterial virulence. Here, we investigated and compared the functional roles of the FnoCas12a and FnoCas9 systems in transformation inhibition and bacterial virulence. Unlike FnoCas9, FnoCas12a was not required for F. novicida virulence. However, both systems were highly effective at plasmid restriction and acted independently of each other. We further identified a critical protospacer-adjacent motif (PAM) necessary for transformation inhibition by FnoCas12a, demonstrating a greater flexibility for target identification by FnoCas12a than previously appreciated and a specificity that is distinct from that of FnoCas9. The effectors of the two systems exhibited different patterns of expression at the mRNA level, suggesting that they may confer distinct benefits to the bacterium in diverse environments. These data suggest that due to the differences between the two CRISPR-Cas systems, together they may provide F. novicida with a more comprehensive defense against foreign nucleic acids. Finally, we demonstrated that the FnoCas12a and FnoCas9 machineries could be simultaneously engineered to restrict the same nonnative target, thereby expanding the toolset for prokaryotic genome manipulation.IMPORTANCE CRISPR-Cas9 and CRISPR-Cas12a systems have been widely commandeered for genome engineering. However, they originate in prokaryotes, where they function as adaptive immune systems. The details of this activity and relationship between these systems within native host organisms have been minimally explored. The human pathogen Francisella novicida contains both of these systems, with the Cas9 system also exhibiting a second activity, modulating virulence through transcriptional regulation. We compared and contrasted the ability of these two systems to control virulence and restrict DNA within their native host bacterium, highlighting differences and similarities in these two functions. Collectively, our results indicate that these two distinct and reprogrammable endogenous systems provide F. novicida with a more comprehensive defense against mobile genetic elements.

Screen for fitness and virulence factors of Francisella sp. strain W12-1067 using amoebae

Francisella tularensis is the causative agent of the human disease referred to as tularemia. Other Francisella species are known but less is understood about their virulence factors. The role of environmental amoebae in the life-cycle of Francisella is still under discussion. Francisella sp. strain W12-1067 (F-W12) is an environmental Francisella isolate recently identified in Germany which is negative for the Francisella pathogenicity island, but exhibits a putative alternative type VI secretion system. Putative virulence factors have been identified in silico in the genome of F-W12. In this work, we established a "scatter screen", used earlier for pathogenic Legionella, to verify experimentally and identify candidate fitness factors using a transposon mutant bank of F-W12 and Acanthamoeba lenticulata as host organism. In these experiments, we identified 79 scatter clones (amoeba sensitive), which were further analyzed by an infection assay identifying 9 known virulence factors, but also candidate fitness factors of F-W12 not yet described as fitness factors in Francisella. The majority of the identified genes encoded proteins involved in the synthesis or maintenance of the cell envelope (LPS, outer membrane, capsule) or in the metabolism (glycolysis, gluconeogenesis, pentose phosphate pathway). Further ¹³C-flux analysis of the Tn5 glucokinase mutant strain revealed that the identified gene indeed encodes the sole active glucokinase in F-W12. In conclusion, candidate fitness factors of the new Francisella species F-W12 were identified using the scatter screen method which might also be usable for other Francisella species.

Two-Component Systems in Francisella Species

Bacteria alter gene expression in response to changes in their environment through various mechanisms that include signal transduction systems. These signal transduction systems use membrane histidine kinase with sensing domains to mediate phosphotransfer to DNA-binding proteins that alter the level of gene expression. Such regulators are called two-component systems (TCSs). TCSs integrate external signals and information from stress pathways, central metabolism and other global regulators, thus playing an important role as part of the overall regulatory network. This review will focus on the knowledge of TCSs in the Gram-negative bacterium, Francisella tularensis, a biothreat agent with a wide range of potential hosts and a significant ability to cause disease. While TCSs have been well-studied in several bacterial pathogens, they have not been well-studied in non-model organisms, such as F. tularensis and its subspecies, whose canonical TCS content surprisingly ranges from few to none. Additionally, of those TCS genes present, many are orphan components, including KdpDE, QseC, QseB/PmrA, and an unnamed two-component system (FTN_1452/FTN_1453). We discuss recent advances in this field related to the role of TCSs in Francisella physiology and pathogenesis and compare the TCS genes present in human virulent versus. environmental species and subspecies of Francisella.

Effects of temperature changes in the transcriptional profile of the emerging fish pathogen Francisella noatunensis subsp. orientalis

One of the major challenges in Nile tilapia (Oreochromis niloticus L.) farming is the occurrence of bacterial infections, and the Francisella noatunensis subsp. orientalis (FNO) is an important pathogen that has emerged in last decades. Francisellosis outbreaks have been reported in the literature as occurring seasonally when water temperature is below 24 °C. The aim of this study was to quantify the median lethal doses (LD₅₀) of FNO in experimental challenges at 28 °C and 22 °C, and to investigate the impact of temperature changes in whole genome expression using microarray technology. The LD₅₀ for Nile tilapia at 28 °C was ∼10^5.7, whereas at 22 °C, the LD₅₀ was ∼10^2.2, showing that the decrease in temperature enhanced disease outcome. Out of 1917 genes screened, a total of 31 and 19 genes were down- and up-regulated at 22 °C, respectively. These genes were grouped by orthology into functional categories of: amino acid, inorganic ion, and carbohydrate transport and metabolism; transcription; and posttranslational modification, protein turnover, and chaperones. Expression of genes related to metabolism, oxidative stress, and thermal shock were regulated by temperature changes, reflecting an ability of FNO to adapt to the environment. Expression of virulence genes usually required for the Francisella genus was not changed between tested temperatures, including that of genes located on the Francisella Pathogenicity Island.

Cross Sectional Study and Risk Factors Analysis of Francisella tularensis in Soil Samples in Punjab Province of Pakistan

Tularemia is an endemic zoonotic disease in many parts of the world including Asia. A cross-sectional study was conducted to determine genome-based prevalence of Francisella tularensis (Ft) in soil, assess an association between its occurrence in soil and likely predictors i.e., macro and micro-nutrients and several categorical variables, and determine seroconversion in small and large ruminants. The study included a total of 2,280 soil samples representing 456 villages in eight districts of the Punjab Province of Pakistan followed by an analysis of serum antibodies in 707 ruminants. The genome of Ft was detected in 3.25% (n = 74, 95% CI: 2.60-4.06) of soil samples. Soluble salts (OR: 1.276, 95% CI: 1.043-1.562, p = 0.015), Ni (OR: 2.910, 95%CI: 0.795-10.644, p = 0.106), Mn (OR:0.733, 95% CI:0.565-0.951, p = 0.019), Zn (OR: 4.922, 95% CI:0.929-26.064, p = 0.061) and nutrients clustered together as PC-1 (OR: 4.76, 95% CI: 2.37-9.54, p = 0.000) and PC-3 (OR: 0.357, 95% CI: 0.640, p = 0.001) were found to have a positive association for the presence of Ft in soil. The odds of occurrence of Ft DNA in soil were higher at locations close to a water source, including canals, streams or drains, [χ2 = 6.7, OR = 1.19, 95% CI:1.05-3.09, p = 0.004] as well as places where animals were present [χ2 = 4.09, OR = 2.06, 95% CI: 1.05-4.05, p = 0.02]. The seroconversion was detected in 6.22% (n = 44, 95% CI: 4.67-8.25) of domestic animals. An occurrence of Ft over a wide geographical region indicates its expansion to enzootic range, and demonstrates the need for further investigation among potential disease reservoirs and at-risk populations, such as farmers and veterinarians.

The prevalence of Francisella spp. in different natural surface water samples collected from northwest of Iran

BACKGROUND AND OBJECTIVES

Francisella tularensis has a wide distribution in northern hemisphere of the world. Up to now, there was little information about the Francisella spp. situation in the environmental samples in Iran. In this study we aimed to determine the prevalence of Francisella spp. in the environmental samples in northwest of Iran.

MATERIALS AND METHODS

A total of 237 natural water samples from ponds, rivers, lakes, springs and other surface waters from north western provinces of Iran (Kurdistan and Western Azerbaijan) were collected from September to November 2015. All samples were cultured for Francisella and other bacterial species and Real Time TaqMan PCR was performed on the concentrated and DNA extracted samples. For detection of the presence of bacterial DNA in the samples, two different targets in the genome of Francisella, ISFtu2 and fopA were used.

RESULTS

Among the tested surface water samples, 40 (17.09%; 95% CI: 12.67-22.33%) and 12 (5.13%; 95%CI: 2.81-8.56%) samples were positive for ISFtu2 and fopA respectively. None of them was positive in culture.

CONCLUSION

The prevalence of Francisella spp. in the environmental samples in the west of Iran is high and it is comparable with Turkey, Iran's neighboring country. Use of higher copy number genes or IS like ISFtu2 could improve the detection of this organism in the environmental samples.

Defining the Metabolic Pathways and Host-Derived Carbon Substrates Required for Francisella tularensis Intracellular Growth

Francisella tularensis is a Gram-negative, facultative, intracellular bacterial pathogen and one of the most virulent organisms known. A hallmark of F. tularensis pathogenesis is the bacterium's ability to replicate to high densities within the cytoplasm of infected cells in over 250 known host species, including humans. This demonstrates that F. tularensis is adept at modulating its metabolism to fluctuating concentrations of host-derived nutrients. The precise metabolic pathways and nutrients utilized by F. tularensis during intracellular growth, however, are poorly understood. Here, we use systematic mutational analysis to identify the carbon catabolic pathways and host-derived nutrients required for F. tularensis intracellular replication. We demonstrate that the glycolytic enzyme phosphofructokinase (PfkA), and thus glycolysis, is dispensable for F. tularensis SchuS4 virulence, and we highlight the importance of the gluconeogenic enzyme fructose 1,6-bisphosphatase (GlpX). We found that the specific gluconeogenic enzymes that function upstream of GlpX varied based on infection model, indicating that F. tularensis alters its metabolic flux according to the nutrients available within its replicative niche. Despite this flexibility, we found that glutamate dehydrogenase (GdhA) and glycerol 3-phosphate (G3P) dehydrogenase (GlpA) are essential for F. tularensis intracellular replication in all infection models tested. Finally, we demonstrate that host cell lipolysis is required for F. tularensis intracellular proliferation, suggesting that host triglyceride stores represent a primary source of glycerol during intracellular replication. Altogether, the data presented here reveal common nutritional requirements for a bacterium that exhibits characteristic metabolic flexibility during infection.IMPORTANCE The widespread onset of antibiotic resistance prioritizes the need for novel antimicrobial strategies to prevent the spread of disease. With its low infectious dose, broad host range, and high rate of mortality, F. tularensis poses a severe risk to public health and is considered a potential agent for bioterrorism. F. tularensis reaches extreme densities within the host cell cytosol, often replicating 1,000-fold in a single cell within 24 hours. This remarkable rate of growth demonstrates that F. tularensis is adept at harvesting and utilizing host cell nutrients. However, like most intracellular pathogens, the types of nutrients utilized by F. tularensis and how they are acquired is not fully understood. Identifying the essential pathways for F. tularensis replication may reveal new therapeutic strategies for targeting this highly infectious pathogen and may provide insight for improved targeting of intracellular pathogens in general.

Molecular Survey of Tularemia and Plague in Small Mammals From Iran

Introduction: Plague and tularemia are zoonoses and their causative bacteria are circulating in certain regions of Iran. This study was conducted to investigate potential disease reservoirs amongst small wildlife species in different regions of Iran.

Methods: Rodents, insectivores and hares from 17 different provinces of the country were collected in 2014 and 2015. Samples were taken from the spleens of the animals and Real-time PCR was applied to detect nucleic acid sequences that are specific to Francisella tularensis and Yersinia pestis, respectively.

Results: Among 140 collected rodents, 25 distinct species were identified out of which five were the most common: Microtus paradoxus (21% out of 140 rodents), Apodemus witherbyi (12%), Microtus irani (11%), Mus musculus (11%) and Microtus socialis (10%). Seventeen insectivores were collected and identified as Crocidura suaveolens (82%) and C. leucodon (18%). Fifty-one hares were collected and identified as Lepus europaeus (57%), Lepus tolai (14%) and Lepus sp. (29%). Three out of 140 explored rodents (1.91%) were positive for F. tularensis, an A. witherbyi, a Mus musculus domesticus, and a Chionomys nivalis collected from Golestan, Khuzestan and Razavi Khorasan provinces, respectively. Two hares (3.92%) were F. tularensis-positive, a L. europaeus from Khuzestan and a Lepus sp. from the Sistan and Baluchistan province. None of the tested animals were positive for Y. pestis.

Conclusion: This is the first report of direct detection of F. tularensis in mammals of Iran and the first-time observation of the agent in a snow vole, C. nivalis worldwide. The results indicate that tularemia is more widespread in Iran than previously reported including the Northeast and Southwestern parts of the country. Future studies should address genetic characterization of F. tularensis positive DNA samples from Iran to achieve molecular subtyping and rule out assay cross-reactivity with near neighbor Francisella species.

Characterization of a Unique Outer Membrane Protein Required for Oxidative Stress Resistance and Virulence of Francisella tularensis

Francisella tularensis, the causative agent of tularemia, lacks typical bacterial virulence factors and toxins but still exhibits extreme virulence. The bacterial multidrug efflux systems consist of an inner membrane, a transmembrane membrane fusion protein, and an outer membrane (OM) component that form a contiguous channel for the secretion of a multitude of bacterial products. Francisella contains three orthologs of the OM proteins; two of these, termed TolC and FtlC, are important for tularemia pathogenesis. The third OM protein, SilC, is homologous to the silver cation efflux protein of other bacterial pathogens. The silC gene (FTL_0686) is located on an operon encoding an Emr-type multidrug efflux pump of F. tularensis The role of SilC in tularemia pathogenesis is not known. In this study, we investigated the role of SilC in secretion and virulence of F. tularensis by generating a silC gene deletion (ΔsilC) mutant and its transcomplemented strain. Our results demonstrate that the ΔsilC mutant exhibits increased sensitivity to antibiotics, oxidants, silver, diminished intramacrophage growth, and attenuated virulence in mice compared to wild-type F. tularensis However, the secretion of antioxidant enzymes of F. tularensis is not impaired in the ΔsilC mutant. The virulence of the ΔsilC mutant is restored in NADPH oxidase-deficient mice, indicating that SilC resists oxidative stress in vivo Collectively, this study demonstrates that the OM component SilC serves a specialized role in virulence of F. tularensis by conferring resistance against oxidative stress and silver.IMPORTANCEFrancisella tularensis, the causative agent of a fatal human disease known as tularemia, is a category A select agent and a potential bioterror agent. The virulence mechanisms of Francisella are not completely understood. This study investigated the role of a unique outer membrane protein, SilC, of a multidrug efflux pump in the virulence of F. tularensis This is the first report demonstrating that the OM component SilC plays an important role in efflux of silver and contributes to the virulence of F. tularensis primarily by providing resistance against oxidative stress. Characterization of these unique virulence mechanisms will provide an understanding of the pathogenesis of tularemia and identification of potential targets for the development of effective therapeutics and prophylactics for protection from this lethal disease.

Zinc Acquisition Mechanisms Differ between Environmental and Virulent Francisella Species

Zinc is an essential nutrient for bacterial growth. Because host cells can restrict pathogen access to zinc as an antimicrobial defense mechanism, intracellular pathogens such as Francisella must sense their environment and acquire zinc in response. In many bacteria, the conserved transcription factor Zur is a key regulator of zinc acquisition. To identify mechanisms of zinc uptake in Francisella novicida U112, transcriptome sequencing of wild-type and putative zur mutant bacteria was performed. Only three genes were confirmed as directly regulated by Zur and zinc limitation by quantitative reverse transcription-PCR. One of these genes, FTN_0879, is predicted to encode a protein with similarity to the zupT family of zinc transporters, which are not typically regulated by Zur. While a putative znuACB operon encoding a high-affinity zinc transporter was identified in U112, expression of this operon was not controlled by Zur or zinc concentration. Disruption of zupT but not znuA in U112 impaired growth under zinc limitation, suggesting that ZupT is the primary mechanism of zinc acquisition under these conditions. In the virulent Francisella tularensis subsp. tularensis Schu S4 strain, zupT is a pseudogene, and attempts to delete znuA were unsuccessful, suggesting that it is essential in this strain. A reverse TetR repression system was used to knock down the expression of znuA in Schu S4, revealing that znuA is required for growth under zinc limitation and contributes to intracellular growth within macrophages. Overall, this work identifies genes necessary for adaptation to zinc limitation and highlights nutritional differences between environmental and virulent Francisella strains.

IMPORTANCEFrancisella tularensis is a tier 1 select agent with a high potential for lethality and no approved vaccine. A better understanding of Francisella virulence factors is required for the development of therapeutics. While acquisition of zinc has been shown to be required for the virulence of numerous intracellular pathogens, zinc uptake has not been characterized in Francisella. This work characterizes the Zur regulon in F. novicida and identifies two transporters that contribute to bacterial growth under zinc limitation. In addition, these data identify differences in mechanisms of zinc uptake and tolerance to zinc limitation between F. tularensis and F. novicida, highlighting the role of znuA in the growth of Schu S4 under zinc limitation.

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.

Francisella novicida inhibits spontaneous apoptosis and extends human neutrophil lifespan

Francisella novicida is a Gram-negative bacterium that is closely related to the highly virulent facultative intracellular pathogen, Francisella tularensis Data published by us and others demonstrate that F. tularensis virulence correlates directly with its ability to impair constitutive apoptosis and extend human neutrophil lifespan. In contrast, F. novicida is attenuated in humans, and the mechanisms that account for this are incompletely defined. Our published data demonstrate that F. novicida binds natural IgG that is present in normal human serum, which in turn, elicits NADPH oxidase activation that does not occur in response to F. tularensis As it is established that phagocytosis and oxidant production markedly accelerate neutrophil death, we predicted that F. novicida may influence the neutrophil lifespan in an opsonin-dependent manner. To test this hypothesis, we quantified bacterial uptake, phosphatidylserine (PS) externalization, and changes in nuclear morphology, as well as the kinetics of procaspase-3, -8, and -9 processing and activation. To our surprise, we discovered that F. novicida not only failed to accelerate neutrophil death but also diminished and delayed apoptosis in a dose-dependent, but opsonin-independent, manner. In keeping with this, studies of conditioned media (CM) showed that neutrophil longevity could be uncoupled from phagocytosis and that F. novicida stimulated neutrophil secretion of CXCL8. Taken together, the results of this study reveal shared and unique aspects of the mechanisms used by Francisella species to manipulate neutrophil lifespan and as such, advance understanding of cell death regulation during infection.

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.

Virulence of representative Japanese Francisella tularensis and immunologic consequences of infection in mice

Francisella tularensis, which causes tularemia, is widely distributed in the Northern hemisphere. F. tularensis strains isolated in Japan are genetically unique from non-Japanese strains; however, their phenotypic properties have not been well studied. Thus, mice were infected with representative Japanese strains of F. tularensis and their virulence and mouse immune responses to them assessed. Of four representative Japanese strains, the Ebina, Jap and Tsuchiya strains were susceptible to H2 O2 and did not grow well intracellularly. Only Yama strain grew intracellularly and was lethal to mice. Infection with Yama strain resulted in drastic increases in IFN-γ, CD4 and CD8 double-positive T cells and Th1 cells (CD3, CD4 and Tim3-positive cells), and a decrease in the ratio of CD8-positive CD4-negative T cells in mice. C57BL/6J mice that survived infection produced IgM antibodies to LPS and IgG2c antibodies to 43, 19 and 17 kDa proteinase K-sensitive components. These data are valuable for understanding the phenotypic properties of F. tularensis in Japan.

Avoidance and Subversion of Eukaryotic Homeostatic Autophagy Mechanisms by Bacterial Pathogens

Autophagy is a conserved lysosomal recycling process, which maintains cellular homeostasis during stress and starvation conditions by degrading and recycling proteins, lipids, and carbohydrates, ultimately increasing nutrient availability in eukaryotes. An additional function of autophagy, termed xenophagy, is to detect, capture, and destroy invading microorganisms, such as viruses, bacteria, and protozoa, providing autophagy with a role in innate immunity. Many intracellular pathogens have, however, developed mechanisms to avoid xenophagy and have evolved strategies to take advantage of select autophagic processes to undergo their intracellular life cycle. This review article will discuss the molecular mechanisms used by the intracellular bacterial pathogens Francisella spp. and Brucella spp. to manipulate components of the autophagic pathway, promoting cytosolic growth in the case of Francisella spp. and facilitating cellular egress and cell-to-cell spread in the case of Brucella spp. These examples highlight how successful, highly infectious bacterial pathogens avoid or subvert host autophagy mechanisms normally employed to maintain eukaryotic homeostasis.

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.

Identification of Genes Required for Secretion of the Francisella Oxidative Burst-Inhibiting Acid Phosphatase AcpA

Francisella tularensis is a Tier 1 bioterror threat and the intracellular pathogen responsible for tularemia in humans and animals. Upon entry into the host, Francisella uses multiple mechanisms to evade killing. Our previous studies have shown that after entering its primary cellular host, the macrophage, Francisella immediately suppresses the oxidative burst by secreting a series of acid phosphatases including AcpA-B-C and HapA, thereby evading the innate immune response of the macrophage and enhancing survival and further infection. However, the mechanism of acid phosphatase secretion by Francisella is still unknown. In this study, we screened for genes required for AcpA secretion in Francisella. We initially demonstrated that the known secretion systems, the putative Francisella-pathogenicity island (FPI)-encoded Type VI secretion system and the Type IV pili, do not secrete AcpA. Using random transposon mutagenesis in conjunction with ELISA, Western blotting and acid phosphatase enzymatic assays, a transposon library of 5450 mutants was screened for strains with a minimum 1.5-fold decrease in secreted (culture supernatant) AcpA, but no defect in cytosolic AcpA. Three mutants with decreased supernatant AcpA were identified. The transposon insertion sites of these mutants were revealed by direct genomic sequencing or inverse-PCR and sequencing. One of these mutants has a severe defect in AcpA secretion (at least 85% decrease) and is a predicted hypothetical inner membrane protein. Interestingly, this mutant also affected the secretion of the FPI-encoded protein, VgrG. Thus, this screen identified novel protein secretion factors involved in the subversion of host defenses.

Towards Development of Improved Serodiagnostics for Tularemia by Use of Francisella tularensis Proteome Microarrays

Tularemia in humans is caused mainly by two subspecies of the Gram-negative facultative anaerobe Francisella tularensis: F. tularensis subsp. tularensis (type A) and F. tularensis subsp. holarctica (type B). The current serological test for tularemia is based on agglutination of whole organisms, and the reactive antigens are not well understood. Previously, we profiled the antibody responses in type A and B tularemia cases in the United States using a proteome microarray of 1,741 different proteins derived from the type A strain Schu S4. Fifteen dominant antigens able to detect antibodies to both types of infection were identified, although these were not validated in a different immunoassay format. Since type A and B subspecies are closely related, we hypothesized that Schu S4 antigens would also have utility for diagnosing type B tularemia caused by strains from other geographic locations. To test this, we probed the Schu S4 array with sera from 241 type B tularemia cases in Spain. Despite there being no type A strains in Spain, we confirmed the responses against some of the same potential serodiagnostic antigens reported previously, as well as determined the responses against additional potential serodiagnostic antigens. Five potential serodiagnostic antigens were evaluated on immunostrips, and two of these (FTT1696/GroEL and FTT0975/conserved hypothetical protein) discriminated between the Spanish tularemia cases and healthy controls. We conclude that antigens from the type A strain Schu S4 are suitable for detection of antibodies from patients with type B F. tularensis infections and that these can be used for the diagnosis of tularemia in a deployable format, such as the immunostrip.

Phylogeography of Francisella tularensis from Tibet, China: Evidence for an asian origin and radiation of holarctica-type Tularemia

The geographical origin and radiation of holarctica-type tularemia, which has spread across the northern hemisphere, is open to scientific debate. Here, through phylogenetics, we show that five Tibetan Francisella tularensis isolates subsp. holarctica cluster between basal-positioned Japanese isolates and all other subspecies strains in the world, providing evidence for a previously unknown intermediate lineage next to the Japanese isolates. Importantly, identification of this new intermediate lineage complements current knowledge of tularemia epidemiology, supporting a geographical origin and radiation of the subsp. holarctica in Asia. In addition, thirteen Tibetan isolates belonging to a clade previously found only in North America and Scandinavia, further increases the diversity of holarctica strains in Asia. In summary, this study provides evidence for an Asian origin and radiation of holarctica-type tularemia.

Morphological analysis of Francisella novicida epithelial cell infections in the absence of functional FipA

Francisella novicida is a surrogate pathogen commonly used to study infections by the potential bioterrorism agent, Francisella tularensis. One of the primary sites of Francisella infections is the liver where >90% of infected cells are hepatocytes. It is known that once Francisella enter cells it occupies a membrane-bound compartment, the Francisella-containing vacuole (FCV), from which it rapidly escapes to replicate in the cytosol. Recent work examining the Francisella disulfide bond formation (Dsb) proteins, FipA and FipB, have demonstrated that these proteins are important during the Francisella infection process; however, details as to how the infections are altered in epithelial cells have remained elusive. To identify the stage of the infections where these Dsbs might act during epithelial infections, we exploited a hepatocyte F. novicida infection model that we recently developed. We found that F. novicida ΔfipA-infected hepatocytes contained bacteria clustered within lysosome-associated membrane protein 1-positive FCVs, suggesting that FipA is involved in the escape of F. novicida from its vacuole. Our morphological evidence provides a tangible link as to how Dsb FipA can influence Francisella infections.

From the Outside-In: The Francisella tularensis Envelope and Virulence

Francisella tularensis is a highly-infectious bacterium that causes the rapid, and often lethal disease, tularemia. Many studies have been performed to identify and characterize the virulence factors that F. tularensis uses to infect a wide variety of hosts and host cell types, evade immune defenses, and induce severe disease and death. This review focuses on the virulence factors that are present in the F. tularensis envelope, including capsule, LPS, outer membrane, periplasm, inner membrane, secretion systems, and various molecules in each of aforementioned sub-compartments. Whereas, no single bacterial molecule or molecular complex single-handedly controls F. tularensis virulence, we review here how diverse bacterial systems work in conjunction to subvert the immune system, attach to and invade host cells, alter phagosome/lysosome maturation pathways, replicate in host cells without being detected, inhibit apoptosis, and induce host cell death for bacterial release and infection of adjacent cells. Given that the F. tularensis envelope is the outermost layer of the bacterium, we highlight herein how many of these molecules directly interact with the host to promote infection and disease. These and future envelope studies are important to advance our collective understanding of F. tularensis virulence mechanisms and offer targets for future vaccine development efforts.

Draft Genome Sequence of Francisella tularensis Strain 410108 from Tibet, China

Francisella tularensis is the etiological agent of the potentially fatal disease tularemia. Here, we report the draft genome sequence of a virulent human isolate from Tibet, China in 1962, F. tularensis strain 410108, an intermediate-genotype strain of F. tularensis subsp. holarctica between biovar japonica and non-japonica strains in the world.

Phylogeographical pattern of Francisella tularensis in a nationwide outbreak of tularaemia in Norway, 2011

In 2011, a nationwide outbreak of tularaemia occurred in Norway with 180 recorded cases. It was associated with the largest peak in lemming density seen in 40 years. Francisella tularensis was isolated from 18 patients. To study the geographical distribution of F. tularensis genotypes in Norway and correlate genotype with epidemiology and clinical presentation, we performed whole genome sequencing of patient isolates. All 18 genomes from the outbreak carried genetic signatures of F. tularensis subsp. holarctica and were assigned to genetic clades using canonical single nucleotide polymorphisms. Ten isolates were assigned to major genetic clade B.6 (subclade B.7), seven to clade B.12, and one to clade B.4. The B.6 subclade B.7 was most common in southern and central Norway, while clade B.12 was evenly distributed between the southern, central and northern parts of the country. There was no association between genotype and clinical presentation of tularaemia, time of year or specimen type. We found extensive sequence similarity with F. tularensis subsp. holarctica genomes from high-endemic tularaemia areas in Sweden. Finding nearly identical genomes across large geographical distances in Norway and Sweden imply a life cycle of the bacterium without replication between the outbreaks and raise new questions about long-range migration mechanisms.

Tularemia, a re-emerging infectious disease in Iran and neighboring countrie

OBJECTIVES

Tularemia is a zoonotic disease transmitted by direct contact with infected animals and through arthropod bites, inhalation of contaminated aerosols, ingestion of contaminated meat or water, and skin contact with any infected material. It is widespread throughout the northern hemisphere, including Iran and its neighbors to the north, northeast, and northwest.

METHODS

In this paper, the epidemiology of tularemia as a re-emerging infectious disease in the world with a focus on Iran and the neighboring countries is reviewed.

RESULTS

In Iran, positive serological tests were first reported in 1973, in wildlife and domestic livestock in the northwestern and southeastern parts of the country. The first human case was reported in 1980 in the southwest of Iran, and recent studies conducted among at-risk populations in the western, southeastern, and southwestern parts of Iran revealed seroprevalences of 14.4, 6.52, and 6%, respectively.

CONCLUSIONS

Several factors may explain the absence of reported tularemia cases in Iran since 1980. Tularemia may be underdiagnosed in Iran because Francisella tularensis subspecies holarctica is likely to be the major etiological agent and usually causes mild to moderately severe disease. Furthermore, tularemia is not a disease extensively studied in the medical educational system in Iran, and empirical therapy may be effective in many cases. Finally, it should be noted that laboratories capable of diagnosing tularemia have only been established in the last few years. Since both recent and older studies have consistently found tularemia antibodies in humans and animals, the surveillance of this disease should receive more attention. In particular, it would be worthwhile for clinical researchers to confirm tularemia cases more often by isolating F. tularensis from infected humans and animals.

Restriction of Francisella novicida genetic diversity during infection of the vector midgut

The genetic diversity of pathogens, and interactions between genotypes, can strongly influence pathogen phenotypes such as transmissibility and virulence. For vector-borne pathogens, both mammalian hosts and arthropod vectors may limit pathogen genotypic diversity (number of unique genotypes circulating in an area) by preventing infection or transmission of particular genotypes. Mammalian hosts often act as “ecological filters” for pathogen diversity, where novel variants are frequently eliminated because of stochastic events or fitness costs. However, whether vectors can serve a similar role in limiting pathogen diversity is less clear. Here we show using Francisella novicida and a natural tick vector of Francisella spp. (Dermacentor andersoni), that the tick vector acted as a stronger ecological filter for pathogen diversity compared to the mammalian host. When both mice and ticks were exposed to mixtures of F. novicida genotypes, significantly fewer genotypes co-colonized ticks compared to mice. In both ticks and mice, increased genotypic diversity negatively affected the recovery of available genotypes. Competition among genotypes contributed to the reduction of diversity during infection of the tick midgut, as genotypes not recovered from tick midguts during mixed genotype infections were recovered from tick midguts during individual genotype infection. Mediated by stochastic and selective forces, pathogen genotype diversity was markedly reduced in the tick. We incorporated our experimental results into a model to demonstrate how vector population dynamics, especially vector-to-host ratio, strongly affected pathogen genotypic diversity in a population over time. Understanding pathogen genotypic population dynamics will aid in identification of the variables that most strongly affect pathogen transmission and disease ecology.

Co-infection, the presence of multiple genotypes of the same pathogen species within an infected individual, is common. Genotype diversity, defined as the number of unique genotypes, and the interaction between genotypes, can strongly influence virulence and pathogen transmission. Understanding how genotypic diversity affects transmission of pathogens that naturally cycle among disparate hosts, such as vector-borne pathogens, is especially important as the capacity of the host and vector to sustain genotypic diversity may differ. To address this, we exposed Dermacentor andersoni ticks, via infected mice, to variably diverse populations of Francisella novicida genotypes. Interestingly, we found that ticks served as greater ecological filters for genotypic diversity compared to mice. This loss in genotypic diversity was due to both stochastic and selective forces. Based on these data and a model, we determined that high numbers of ticks in an environment support high genotypic diversity, while genotypic diversity will be lost rapidly in environments with low tick numbers. Together, these results provide evidence that vector population dynamics, vector-to-host ratios, and competition among pathogen genotypes play critical roles in the maintenance of pathogen genotypic diversity.

An overview: tularemia and travel medicine

Tularemia is a bacterial zoonotic infection. The disease is endemic in most parts of the world, has been reported through the northern hemisphere between 30 and 71° N latitude. Francisella tularensis causes infection in a wide range of vertebrates (rodents, lagomorphs) and invertebrates (ticks, mosquitoes and other arthropods). Humans can acquire this infection through several routes including; a bite from an infected tick, deerfly or mosquito, contact with an infected animal or its dead body. It can also be spread to human by drinking contaminated water or breathing contaminated dirt or aerosol. Clinical manifestation of this disease varies depending on the biotype, inoculum and port of entry. Infection is potentially life threatening, but can effectively be treated with antibiotics. Travelers visiting rural and agricultural areas in endemic countries may be at greater risk. Appropriate clothing and use of insect repellants is essential to prevent tick borne illness. Travelers also should be aware of food and waterborne disease; avoid consuming potentially contaminated water and uncooked meat. Physicians should be aware of any clinical presentation of tularemia in the patients returning from endemic areas.

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.

Comparative genomic analysis reveals distinct genotypic features of the emerging pathogen Haemophilus influenzae type f

BACKGROUND

The incidence of invasive disease caused by encapsulated Haemophilus influenzae type f (Hif) has increased in the post-H. influenzae type b (Hib) vaccine era. We previously annotated the first complete Hif genome from a clinical isolate (KR494) that caused septic shock and necrotizing myositis. Here, the full genome of Hif KR494 was compared to sequenced reference strains Hib 10810, capsule type d (Hid) Rd Kw20, and finally nontypeable H. influenzae 3655. The goal was to identify possible genomic characteristics that may shed light upon the pathogenesis of Hif.

RESULTS

The Hif KR494 genome exhibited large regions of synteny with other H. influenzae, but also distinct genome rearrangements. A predicted Hif core genome of 1390 genes was shared with the reference strains, and 6 unique genomic regions comprising half of the 191 unique coding sequences were revealed. The majority of these regions were inserted genetic fragments, most likely derived from the closely-related Haemophilus spp. including H. aegyptius, H. haemolyticus and H. parainfluenzae. Importantly, the KR494 genome possessed several putative virulence genes that were distinct from non-type f strains. These included the sap2 operon, aef3 fimbriae, and genes for kanamycin nucleotidyltranserase, iron-utilization proteins, and putative YadA-like trimeric autotransporters that may increase the bacterial virulence. Furthermore, Hif KR494 lacked a hisABCDEFGH operon for de novo histidine biosynthesis, hmg locus for lipooligosaccharide biosynthesis and biofilm formation, the Haemophilus antibiotic resistance island and a Haemophilus secondary molybdate transport system. We confirmed the histidine auxotrophy and kanamycin resistance in Hif by functional experiments. Moreover, the pattern of unique or missing genes of Hif KR494 was similar in 20 Hif clinical isolates obtained from different years and geographical areas. A cross-species comparison revealed that the Hif genome shared more characteristics with H. aegyptius than Hid and NTHi.

CONCLUSIONS

The genomic comparative analyses facilitated identification of genotypic characteristics that may be related to the specific virulence of Hif. In relation to non-type f H. influenzae strains, the Hif genome contains differences in components involved in metabolism and survival that may contribute to its invasiveness.

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.

Francisella tularensis intracellular survival: to eat or to die

Francisella tularensis is a highly infectious facultative intracellular bacterium causing the zoonotic disease tularemia. Numerous attributes required for F. tularensis intracellular multiplication have been identified recently. However, the mechanisms by which the majority of them interfere with the infected host are still poorly understood. The following review summarizes our current knowledge on the different steps of Francisella intramacrophagic life cycle and expands on the importance of nutrient acquisition.

Synthetic promoters functional in Francisella novicida and Escherichia coli

In this work, we describe the identification of synthetic, controllable promoters that function in the bacterial pathogen Francisella novicida, a model facultative intracellular pathogen. Synthetic DNA fragments consisting of the tetracycline operator (tetO) flanked by a random nucleotide sequence were inserted into a Francisella novicida shuttle plasmid upstream of a promoterless artificial operon containing the reporter genes cat and lacZ. Fragments able to promote transcription were selected for based on their ability to drive expression of the cat gene, conferring chloramphenicol resistance. Promoters of various strengths were found, many of which were repressed in the presence of the tetracycline repressor (TetR) and promoted transcription only in the presence of the TetR inducer anhydrotetracycline. A subset of both constitutive and inducible synthetic promoters were characterized to find their induction ratios and to identify their transcription start sites. In cases where tetO was located between or downstream of the -10 and -35 regions of the promoter, control by TetR was observed. If the tetO region was upstream of the -35 region by more than 9 bp, it did not confer TetR control. We found that three of three promoters isolated in F. novicida functioned at a comparable level in E. coli; however, none of the 10 promoters isolated in E. coli functioned at a significant level in F. novicida. Our results allowed us to isolate minimal F. novicida promoters of 47 and 48 bp in length.

Degeneration of a CRISPR/Cas system and its regulatory target during the evolution of a pathogen

CRISPR/Cas systems are bacterial RNA-guided endonuclease machineries that target foreign nucleic acids. Recently, we demonstrated that the Cas protein Cas9 controls gene expression and virulence in Francisella novicida by altering the stability of the mRNA for an immunostimulatory bacterial lipoprotein (BLP). Genomic analyses, however, revealed that Francisella species with increased virulence harbor degenerated CRISPR/Cas systems. We hypothesize that CRISPR/Cas degeneration removed a barrier against genome alterations, which resulted in enhanced virulence. Importantly, the BLP locus was also lost; likely a necessary adaptation in the absence of Cas9-mediated repression. CRISPR/Cas systems likely play regulatory roles in numerous bacteria, and these data suggest additional genomic changes may be required to maintain fitness after CRISPR/Cas loss in such bacteria, having important evolutionary implications.

Comparative phosphoproteomics reveals components of host cell invasion and post-transcriptional regulation during Francisella infection

Francisella tularensis is a facultative intracellular bacterium that causes the deadly disease tularemia. Most evidence suggests that Francisella is not well recognized by the innate immune system that normally leads to cytokine expression and cell death. In previous work, we identified new bacterial factors that were hyper-cytotoxic to macrophages. Four of the identified hyper-cytotoxic strains (lpcC, manB, manC, and kdtA) had an impaired lipopolysaccharide (LPS) synthesis and produced an exposed lipid A lacking the O-antigen. These mutants were not only hyper-cytotoxic but also were phagocytosed at much higher rates compared with the wild type parent strain. To elucidate the cellular signaling underlying this enhanced phagocytosis and cell death, we performed a large-scale comparative phosphoproteomic analysis of cells infected with wild-type and delta-lpcC F. novicida. Our data suggest that not only actin but also intermediate filaments and microtubules are important for F. novicida entry into the host cells. In addition, we observed differential phosphorylation of tristetraprolin, a key component of the mRNA-degrading machinery that controls the expression of a variety of genes including many cytokines. Infection with the delta-lpcC mutant induced the hyper-phosphorylation and inhibition of tristetraprolin, leading to the production of cytokines such as IL-1beta and TNF-alpha that may kill the host cells by triggering apoptosis. Together, our data provide new insights for Francisella invasion and a post-transcriptional mechanism that prevents the expression of host immune response factors that control infection by this pathogen.

Generation of protection against Francisella novicida in mice depends on the pathogenicity protein PdpA, but not PdpC or PdpD

Previous results suggest that mutations in most genes in the Francisella pathogenicity island (FPI) attenuate the bacterium. Using a mouse model, here we determined the impact of mutations in pdpA, pdpC, and pdpD in Francisella novicida on in vitro replication in macrophages, and in vivo immunogenicity. In contrast to most FPI genes, deletion of pdpC (FnΔpdpC) and pdpD (FnΔpdpD) from F. novicida did not impact growth in mouse bone-marrow derived macrophages. Nonetheless, both FnΔpdpC and FnΔpdpD were highly attenuated when administered intradermally. Infected mice produced relatively normal anti-F. novicida serum antibodies. Further, splenocytes from infected mice controlled intramacrophage Francisella replication, indicating T cell priming, and mice immunized by infection with FnΔpdpC or FnΔpdpD survived secondary lethal parenteral challenge with either F. novicida or Francisella tularensis LVS. In contrast, deletion of pdpA (FnΔpdpA) ablated growth in macrophages in vitro. FnΔpdpA disseminated and replicated poorly in infected mice, accompanied by development of some anti-F. novicida serum antibodies. However, primed Th1 cells were not detected, and vaccinated mice did not survive even low dose challenge with either F. novicida or LVS. Taken together, these results suggest that successful priming of Th1 cells, and protection against lethal challenge, depends on expression of PdpA.

Proteins involved in Francisella tularensis survival and replication inside macrophages

Francisella tularensis, the etiological agent of tularemia, is a member of the γ-proteobacteria class of Gram-negative bacteria. This highly virulent bacterium can infect a large range of mammalian species and has been recognized as a human pathogen for a century. F. tularensis is able to survive in vitro in a variety of cell types. In vivo, the bacterium replicates mainly in infected macrophages, using the cytoplasmic compartment as a replicative niche. To successfully adapt to this stressful environment, F. tularensis must simultaneously: produce and regulate the expression of a series of dedicated virulence factors; adapt its metabolic needs to the nutritional context of the host cytosol; and control the innate immune cytosolic surveillance pathways to avoid premature cell death. We will focus here on the secretion or release of bacterial proteins in the host, as well as on the envelope proteins, involved in bacterial survival inside macrophages.

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.

PanG, a new ketopantoate reductase involved in pantothenate synthesis

Pantothenate, commonly referred to as vitamin B(5), is an essential molecule in the metabolism of living organisms and forms the core of coenzyme A. Unlike humans, some bacteria and plants are capable of de novo biosynthesis of pantothenate, making this pathway a potential target for drug development. Francisella tularensis subsp. tularensis Schu S4 is a zoonotic bacterial pathogen that is able to synthesize pantothenate but is lacking the known ketopantoate reductase (KPR) genes, panE and ilvC, found in the canonical Escherichia coli pathway. Described herein is a gene encoding a novel KPR, for which we propose the name panG (FTT1388), which is conserved in all sequenced Francisella species and is the sole KPR in Schu S4. Homologs of this KPR are present in other pathogenic bacteria such as Enterococcus faecalis, Coxiella burnetii, and Clostridium difficile. Both the homologous gene from E. faecalis V583 (EF1861) and E. coli panE functionally complemented Francisella novicida lacking any KPR. Furthermore, panG from F. novicida can complement an E. coli KPR double mutant. A Schu S4 ΔpanG strain is a pantothenate auxotroph and was genetically and chemically complemented with panG in trans or with the addition of pantolactone. There was no virulence defect in the Schu S4 ΔpanG strain compared to the wild type in a mouse model of pneumonic tularemia. In summary, we characterized the pantothenate pathway in Francisella novicida and F. tularensis and identified an unknown and previously uncharacterized KPR that can convert 2-dehydropantoate to pantoate, PanG.

Whole genome sequencing of the fish pathogen Francisella noatunensis subsp. orientalis Toba04 gives novel insights into Francisella evolution and pathogenecity

BACKGROUND

Francisella is a genus of gram-negative bacterium highly virulent in fishes and human where F. tularensis is causing the serious disease tularaemia in human. Recently Francisella species have been reported to cause mortality in aquaculture species like Atlantic cod and tilapia. We have completed the sequencing and draft assembly of the Francisella noatunensis subsp. orientalisToba04 strain isolated from farmed Tilapia. Compared to other available Francisella genomes, it is most similar to the genome of Francisella philomiragia subsp. philomiragia, a free-living bacterium not virulent to human.

RESULTS

The genome is rearranged compared to the available Francisella genomes even though we found no IS-elements in the genome. Nearly 16% percent of the predicted ORFs are pseudogenes. Computational pathway analysis indicates that a number of the metabolic pathways are disrupted due to pseudogenes. Comparing the novel genome with other available Francisella genomes, we found around 2.5% of unique genes present in Francisella noatunensis subsp. orientalis Toba04 and a list of genes uniquely present in the human-pathogenic Francisella subspecies. Most of these genes might have transferred from bacterial species through horizontal gene transfer. Comparative analysis between human and fish pathogen also provide insights into genes responsible for pathogenecity. Our analysis of pseudogenes indicates that the evolution of Francisella subspecies's pseudogenes from Tilapia is old with large number of pseudogenes having more than one inactivating mutation.

CONCLUSIONS

The fish pathogen has lost non-essential genes some time ago. Evolutionary analysis of the Francisella genomes, strongly suggests that human and fish pathogenic Francisella species have evolved independently from free-living metabolically competent Francisella species. These findings will contribute to understanding the evolution of Francisella species and pathogenesis.