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.

Insights into DNA-binding motifs and mechanisms of Francisella tularensis novicida two-component system response regulator proteins QseB, KdpE, and BfpR.

Two component system bacterial response regulators are typically DNA-binding proteins which enable the genetic regulation of many adaptive bacterial behaviors. Despite structural similarity across response regulator families, there is a diverse array of DNA-binding mechanisms. Bacteria usually encode several dozen two-component system response regulators, but Francisella tularensis only encodes three. Due to their simplified response regulatory network, Francisella species are a model for studying the role of response regulator proteins in virulence. Here, we show that Francisella response regulators QseB, KdpE, and BfpR all utilize different DNA-binding mechanisms. Our evidence suggests that QseB follows a simple mechanism whereby it binds a single inverted repeat sequence with a higher affinity upon phosphorylation. This behavior is independent of whether QseB is a positive or negative regulator of the gene as demonstrated by qseB and priM promoter sequences, respectively. Similarly, KdpE binds DNA more tightly upon phosphorylation, but also exhibits a cooperative binding isotherm. While we propose a KdpE binding site, it is possible that KdpE has a complex DNA-binding mechanism potentially involving multiple copies of KdpE being recruited to a promoter region. Finally, we show that BfpR appears to bind a region of its own promoter sequence with a lower affinity upon phosphorylation. Further structural and enzymatic work will need to be performed to deconvolute the KdpE and BfpR binding mechanisms.

The geographic distribution, and the biotic and abiotic predictors of select zoonotic pathogen detections in Canadian polar bears.

Increasing Arctic temperatures are facilitating the northward expansion of more southerly hosts, vectors, and pathogens, exposing naïve populations to pathogens not typical at northern latitudes. To understand such rapidly changing host-pathogen dynamics, we need sensitive and robust surveillance tools. Here, we use a novel multiplexed magnetic-capture and droplet digital PCR (ddPCR) tool to assess a sentinel Arctic species, the polar bear (Ursus maritimus; n = 68), for the presence of five zoonotic pathogens (Erysipelothrix rhusiopathiae, Francisella tularensis, Mycobacterium tuberculosis complex, Toxoplasma gondii and Trichinella spp.), and observe associations between pathogen presence and biotic and abiotic predictors. We made two novel detections: the first detection of a Mycobacterium tuberculosis complex member in Arctic wildlife and the first of E. rhusiopathiae in a polar bear. We found a prevalence of 37% for E. rhusiopathiae, 16% for F. tularensis, 29% for Mycobacterium tuberculosis complex, 18% for T. gondii, and 75% for Trichinella spp. We also identify associations with bear age (Trichinella spp.), harvest season (F. tularensis and MTBC), and human settlements (E. rhusiopathiae, F. tularensis, MTBC, and Trichinella spp.). We demonstrate that monitoring a sentinel species, the polar bear, could be a powerful tool in disease surveillance and highlight the need to better characterize pathogen distributions and diversity in the Arctic.

Polyaminated, acetylated and stop codon readthrough of recombinant Francisella tularensis universal stress protein in Escherichia coli.

Recombinant Francisella tularensis universal stress protein with a C-terminal histidine-tag (rUsp/His6) was expressed in Escherichia coli. Endogenous F. tularensis Usp has a predicted molecular mass of 30 kDa, but rUsp/His6 had an apparent molecular weight of 33 kDa based on Western blot analyses. To determine the source of the higher molecular weight for rUsp/His6, post translational modifications were examined. Tryptic peptides of purified rUsp/His6 were subjected to liquid chromatography tandem mass spectrometry (LC-MS/MS) and fragmentation spectra were searched for acetylated lysines and polyaminated glutamines. Of the 24 lysines in rUsp/His6, 10 were acetylated (K63, K68, K72, K129, K175, K201, K208, K212, K233, and K238) and three of the four glutamines had putrescine, spermidine and spermine adducts (Q55, Q60 and Q267). The level of post-translational modification was substoichiometric, eliminating the possibility that these modifications were the sole contributor to the 3 kDa extra mass of rUsp/His6. LC-MS/MS revealed that stop codon readthrough had occurred resulting in the unexpected addition of 20 extra amino acids at the C-terminus of rUsp/His6, after the histidine tag. Further, the finding of polyaminated glutamines in rUsp/His6 indicated that E. coli is capable of transglutaminase activity.

Dual proteomics of infected macrophages reveal bacterial and host players involved in the Francisella intracellular life cycle and cell to cell dissemination by merocytophagy.

Bacterial pathogens adapt and replicate within host cells, while host cells develop mechanisms to eliminate them. Using a dual proteomic approach, we characterized the intra-macrophage proteome of the facultative intracellular pathogen, Francisella novicida. More than 900 Francisella proteins were identified in infected macrophages after a 10-h infection. Biotin biosynthesis-related proteins were upregulated, emphasizing the role of biotin-associated genes in Francisella replication. Conversely, proteins encoded by the Francisella pathogenicity island (FPI) were downregulated, supporting the importance of the F. tularensis Type VI Secretion System for vacuole escape, not cytosolic replication. In the host cell, over 300 proteins showed differential expression among the 6200 identified during infection. The most upregulated host protein was cis-aconitate decarboxylase IRG1, known for itaconate production with antimicrobial properties in Francisella. Surprisingly, disrupting IRG1 expression did not impact Francisella's intracellular life cycle, suggesting redundancy with other immune proteins or inclusion in larger complexes. Over-representation analysis highlighted cell-cell contact and actin polymerization in macrophage deregulated proteins. Using flow cytometry and live cell imaging, we demonstrated that merocytophagy involves diverse cell-to-cell contacts and actin polymerization-dependent processes. These findings lay the groundwork for further exploration of merocytophagy and its molecular mechanisms in future research.Data are available via ProteomeXchange with identifier PXD035145.

Identification of an N-terminal tag (580N) that improves the biosynthesis of fluorescent proteins in Francisella tularensis and other Gram-negative bacteria.

Utilization of fluorescent proteins is widespread for the study of microbial pathogenesis and host-pathogen interactions. Here, we discovered that linkage of the 36 N-terminal amino acids of FTL_0580 (a hypothetical protein of Francisella tularensis) to fluorescent proteins increases the fluorescence emission of bacteria that express these recombinant fusions. This N-terminal peptide will be referred to as 580N. Western blotting revealed that the linkage of 580N to Emerald Green Fluorescent Protein (EmGFP) in F. tularensis markedly improved detection of this protein. We therefore hypothesized that transcripts containing 580N may be translated more efficiently than those lacking the coding sequence for this leader peptide. In support, expression of emGFPFt that had been codon-optimized for F. tularensis, yielded significantly enhanced fluorescence than its non-optimized counterpart. Furthermore, fusing emGFP with coding sequence for a small N-terminal peptide (Serine-Lysine-Isoleucine-Lysine), which had previously been shown to inhibit ribosomal stalling, produced robust fluorescence when expressed in F. tularensis. These findings support the interpretation that 580N enhances the translation efficiency of fluorescent proteins in F. tularensis. Interestingly, expression of non-optimized 580N-emGFP produced greater fluorescence intensity than any other construct. Structural predictions suggested that RNA secondary structure also may be influencing translation efficiency. When expressed in Escherichia coli and Klebsiella pneumoniae bacteria, 580N-emGFP produced increased green fluorescence compared to untagged emGFP (neither allele was codon optimized for these bacteria). In conclusion, fusing the coding sequence for the 580N leader peptide to recombinant genes might serve as an economical alternative to codon optimization for enhancing protein expression in bacteria.

Mutation of wbtJ, a N-formyltransferase involved in O-antigen synthesis, results in biofilm formation, phase variation and attenuation in Francisella tularensis.

Two clinically important subspecies, Francisella tularensis subsp. tularensis (type A) and F. tularensis subsp. holarctica (type B) are responsible for most tularaemia cases, but these isolates typically form a weak biofilm under in vitro conditions. Phase variation of the F. tularensis lipopolysaccharide (LPS) has been reported in these subspecies, but the role of variation is unclear as LPS is crucial for virulence. We previously demonstrated that a subpopulation of LPS variants can constitutively form a robust biofilm in vitro, but it is unclear whether virulence was affected. In this study, we show that biofilm-forming variants of both fully virulent F. tularensis subspecies were highly attenuated in the murine tularaemia model by multiple challenge routes. Genomic sequencing was performed on these strains, which revealed that all biofilm-forming variants contained a lesion within the wbtJ gene, a formyltransferase involved in O-antigen synthesis. A ΔwbtJ deletion mutant recapitulated the biofilm, O-antigen and virulence phenotypes observed in natural variants and could be rescued through complementation with a functional wbtJ gene. Since the spontaneously derived biofilm-forming isolates in this study were a subpopulation of natural variants, reversion events to the wbtJ gene were detected that eliminated the phenotypes associated with biofilm variants and restored virulence. These results demonstrate a role for WbtJ in biofilm formation, LPS variation and virulence of F. tularensis.

The rLVS ΔcapB/iglABC vaccine provides potent protection in Fischer rats against inhalational tularemia caused by various virulent Francisella tularensis strains.

Francisella tularensis is one of the several biothreat agents for which a licensed vaccine is needed. To ensure vaccine protection is achieved across a range of virulent F. tularensis strains, we assembled and characterized a panel of F. tularensis isolates to be utilized as challenge strains. A promising tularemia vaccine candidate is rLVS ΔcapB/iglABC (rLVS), in which the vector is the LVS strain with a deletion in the capB gene and which additionally expresses a fusion protein comprising immunodominant epitopes of proteins IglA, IglB, and IglC. Fischer rats were immunized subcutaneously 1-3 times at 3-week intervals with rLVS at various doses. The rats were exposed to a high dose of aerosolized Type A strain Schu S4 (FRAN244), a Type B strain (FRAN255), or a tick derived Type A strain (FRAN254) and monitored for survival. All rLVS vaccination regimens including a single dose of 107 CFU rLVS provided 100% protection against both Type A strains. Against the Type B strain, two doses of 107 CFU rLVS provided 100% protection, and a single dose of 107 CFU provided 87.5% protection. In contrast, all unvaccinated rats succumbed to aerosol challenge with all of the F. tularensis strains. A robust Th1-biased antibody response was induced in all vaccinated rats against all F. tularensis strains. These results demonstrate that rLVS ΔcapB/iglABC provides potent protection against inhalational challenge with either Type A or Type B F. tularensis strains and should be considered for further analysis as a future tularemia vaccine.

Pathogenicity and virulence of Francisella tularensis.

Tularaemia is a zoonotic disease caused by the Gram-negative bacterium, Francisella tularensis. Depending on its entry route into the organism, F. tularensis causes different diseases, ranging from life-threatening pneumonia to less severe ulceroglandular tularaemia. Various strains with different geographical distributions exhibit different levels of virulence. F. tularensis is an intracellular bacterium that replicates primarily in the cytosol of the phagocytes. The main virulence attribute of F. tularensis is the type 6 secretion system (T6SS) and its effectors that promote escape from the phagosome. In addition, F. tularensis has evolved a peculiar envelope that allows it to escape detection by the immune system. In this review, we cover tularaemia, different Francisella strains, and their pathogenicity. We particularly emphasize the intracellular life cycle, associated virulence factors, and metabolic adaptations. Finally, we present how F. tularensis largely escapes immune detection to be one of the most infectious and lethal bacterial pathogens.

Positive Francisella tularensis meningitis outcome despite delayed identification: a case report.

Francisella tularensis is a Gram-negative bacteria, that may cause a zoonotic disease, tularemia. Here, we describe a patient case, where a previously healthy young woman in Northern Finland contacted health care because of fever and headache. Due to the symptoms and lack of further diagnostic tools in primary health care, she was transferred to University Hospital (UH) where ampicillin and ceftriaxone was given empirically. A cerebrospinal fluid sample (CSF) was drawn showing small Gram-negative rods that grew on chocolate agar after 2 days of incubation. Matrix-assisted laser-desorption-ionization time of-flight (Maldi-tof) did not provide identification, but the bacteria was interpreted as sensitive to ciprofloxacin and the treatment was changed to ciprofloxacin. During the time the patient was infected, there were several positive tularemia samples found in the area. Therefore, an in house tularemia nucleic acid method (PCR) was used on the bacterial culture. Additionally, 16S rDNA sequencing was performed and these methods identified the bacteria as F. tularensis. Fortunately, the patient recovered completely with ciprofloxacin and was discharged without any complications. Our case underlines the need to understand the limits of specific diagnostic methods, such as Maldi-tof, used in clinical laboratory settings. It also highlights the need of both clinicians and laboratory staff to be aware of the many clinical presentations of tularemia when working in an endemic area.

Three Cases of Tickborne Francisella tularensis Infection, Austria, 2022.

Tularemia is increasing in Austria. We report Francisella tularensis subspecies holarctica isolated from 3 patients who had been bitten by arthropods. Next-generation sequencing showed substantial isolate similarity. Clinicians should consider bloodstream F. tularensis infections for patients with signs/symptoms of ulceroglandular tularemia, and surveillance of potential vectors should be intensified.

Ribosome heterogeneity results in leader sequence-mediated regulation of protein synthesis in Francisella tularensis.

Although ribosomes are generally examined in aggregate, ribosomes can be heterogenous in composition. Evidence is accumulating that changes in ribosome composition may result in altered function, such that ribosome heterogeneity may provide a mechanism to regulate protein synthesis. Ribosome heterogeneity in the human pathogen Francisella tularensis results from incorporation of one of three homologs of bS21, a small ribosomal subunit protein demonstrated to regulate protein synthesis in other bacteria. Loss of one homolog, bS21-2, results in genome-wide post-transcriptional changes in protein abundance. This suggests that bS21-2 can, either directly or indirectly, lead to preferential translation of particular mRNAs. Here, we examine the potential of bS21-2 to function in a leader sequence-dependent manner and to function indirectly, via Hfq. We found that the 5´ untranslated region (UTR) of some bS21-2-responsive genes, including key virulence genes, is sufficient to alter translation in cells lacking bS21-2. We further identify features of a 5´ UTR that allow responsiveness to bS21-2. These include an imperfect Shine-Dalgarno sequence and a particular six nucleotide sequence. Our results are consistent with a model in which a bS21 homolog increases the efficiency of translation initiation through interactions with specific leader sequences. With respect to bS21-2 indirectly regulating translation via the RNA-binding protein Hfq, we found that Hfq controls transcript abundance rather than protein synthesis, impacting virulence gene expression via a distinct mechanism. Together, we determined that ribosome composition in F. tularensis regulates translation in a leader sequence-dependent manner, a regulatory mechanism which may be used in other bacteria. IMPORTANCE Ribosome heterogeneity is common in bacteria, and there is mounting evidence that ribosome composition plays a regulatory role in protein synthesis. However, mechanisms of ribosome-driven gene regulation are not well understood. In the human pathogen Francisella tularensis, which encodes multiple homologs for the ribosomal protein bS21, loss of one homolog impacts protein synthesis and virulence. Here, we explore the mechanism behind bS21-mediated changes in protein synthesis, finding that they can be linked to altered translation initiation and are dependent on specific sequences in the leaders of transcripts. Our data support a model in which ribosome composition regulates gene expression through translation, a strategy that may be conserved in diverse organisms with various sources of ribosome heterogeneity.

Human Tularemia Epididymo-Orchitis Caused by Francisella tularensis Subspecies holartica, Austria.

A previously healthy man in Austria had tularemia epididymo-orchitis develop, leading to unilateral orchiectomy. Francisella tularensis subspecies holartica was detected by 16S rRNA gene sequencing analysis of inflamed granulomatous testicular tissue. Clinicians should suspect F. tularensis as a rare etiologic microorganism in epididymo-orchitis patients with relevant risk factors.

Discovery of a glutathione utilization pathway in Francisella that shows functional divergence between environmental and pathogenic species.

Glutathione (GSH) is an abundant metabolite within eukaryotic cells that can act as a signal, a nutrient source, or serve in a redox capacity for intracellular bacterial pathogens. For Francisella, GSH is thought to be a critical in vivo source of cysteine; however, the cellular pathways permitting GSH utilization by Francisella differ between strains and have remained poorly understood. Using genetic screening, we discovered a unique pathway for GSH utilization in Francisella. Whereas prior work suggested GSH catabolism initiates in the periplasm, the pathway we define consists of a major facilitator superfamily (MFS) member that transports intact GSH and a previously unrecognized bacterial cytoplasmic enzyme that catalyzes the first step of GSH degradation. Interestingly, we find that the transporter gene for this pathway is pseudogenized in pathogenic Francisella, explaining phenotypic discrepancies in GSH utilization among Francisella spp. and revealing a critical role for GSH in the environmental niche of these bacteria.

Ex vivo infection model for Francisella using human lung tissue.

Introduction: Tularemia is mainly caused by Francisella tularensis (Ft) subsp. tularensis (Ftt) and Ft subsp. holarctica (Ftt) in humans and in more than 200 animal species including rabbits and hares. Human clinical manifestations depend on the route of infection and range from flu-like symptoms to severe pneumonia with a mortality rate up to 60% without treatment. So far, only 2D cell culture and animal models are used to study Francisella virulence, but the gained results are transferable to human infections only to a certain extent. Method: In this study, we firstly established an ex vivo human lung tissue infection model using different Francisella strains: Ftt Life Vaccine Strain (LVS), Ftt LVS ΔiglC, Ftt human clinical isolate A-660 and a German environmental Francisella species strain W12-1067 (F-W12). Human lung tissue was used to determine the colony forming units and to detect infected cell types by using spectral immunofluorescence and electron microscopy. Chemokine and cytokine levels were measured in culture supernatants. Results: Only LVS and A-660 were able to grow within the human lung explants, whereas LVS ΔiglC and F-W12 did not replicate. Using human lung tissue, we observed a greater increase of bacterial load per explant for patient isolate A-660 compared to LVS, whereas a similar replication of both strains was observed in cell culture models with human macrophages. Alveolar macrophages were mainly infected in human lung tissue, but Ftt was also sporadically detected within white blood cells. Although Ftt replicated within lung tissue, an overall low induction of pro-inflammatory cytokines and chemokines was observed. A-660-infected lung explants secreted slightly less of IL-1ß, MCP-1, IP-10 and IL-6 compared to Ftt LVS-infected explants, suggesting a more repressed immune response for patient isolate A-660. When LVS and A-660 were used for simultaneous co-infections, only the ex vivo model reflected the less virulent phenotype of LVS, as it was outcompeted by A-660. Conclusion: We successfully implemented an ex vivo infection model using human lung tissue for Francisella. The model delivers considerable advantages and is able to discriminate virulent Francisella from less- or non-virulent strains and can be used to investigate the role of specific virulence factors.

Comparative evaluation of protective immunity against Francisella tularensis induced by subunit or adenovirus-vectored vaccines.

Tularemia is a highly contagious disease caused by infection with Francisella tularensis (Ft), a pathogenic intracellular gram-negative bacterium that infects a wide range of animals and causes severe disease and death in people, making it a public health concern. Vaccines are the most effective way to prevent tularemia. However, there are no Food and Drug Administration (FDA)-approved Ft vaccines thus far due to safety concerns. Herein, three membrane proteins of Ft, Tul4, OmpA, and FopA, and a molecular chaperone, DnaK, were identified as potential protective antigens using a multifactor protective antigen platform. Moreover, the recombinant DnaK, FopA, and Tul4 protein vaccines elicited a high level of IgG antibodies but did not protect against challenge. In contrast, protective immunity was elicited by a replication-defective human type 5 adenovirus (Ad5) encoding the Tul4, OmpA, FopA, and DnaK proteins (Ad5-Tul4, Ad5-OmpA, Ad5-FopA, and Ad5-DnaK) after a single immunization, and all Ad5-based vaccines stimulated a Th1-biased immune response. Moreover, intramuscular and intranasal vaccination with Ad5-Tul4 using the prime-boost strategy effectively eliminated Ft lung, spleen and liver colonization and provided nearly 80% protection against intranasal challenge with the Ft live vaccine strain (LVS). Only intramuscular, not intranasal vaccination, with Ad5-Tul4 protected mice from intraperitoneal challenge. This study provides a comprehensive comparison of protective immunity against Ft provided by subunit or adenovirus-vectored vaccines and suggests that mucosal vaccination with Ad5-Tul4 may yield desirable protective efficacy against mucosal infection, while intramuscular vaccination offers greater overall protection against intraperitoneal tularemia.

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.

Case Studies and Literature Review of Francisella tularensis-Related Prosthetic Joint Infection.

Tularemia is a zoonotic infection caused by Francisella tularensis. Its most typical manifestations in humans are ulceroglandular and glandular; infections in prosthetic joints are rare. We report 3 cases of F. tularensis subspecies holarctica-related prosthetic joint infection that occurred in France during 2016-2019. We also reviewed relevant literature and found only 5 other cases of Francisella-related prosthetic joint infections worldwide, which we summarized. Among those 8 patients, clinical symptoms appeared 7 days to 19 years after the joint placement and were nonspecific to tularemia. Although positive cultures are typically obtained in only 10% of tularemia cases, strains grew in all 8 of the patients. F. tularensis was initially identified in 2 patients by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry; molecular methods were used for 6 patients. Surgical treatment in conjunction with long-term antimicrobial treatment resulted in favorable outcomes; no relapses were seen after 6 months of follow-up.

Complete genome sequence of the emerging pathogen Cysteiniphilum spp. and comparative genomic analysis with genus Francisella: Insights into its genetic diversity and potential virulence traits.

Cysteiniphilum is a newly discovered genus in 2017 and is phylogenetically closely related to highly pathogenic Francisella tularensis. Recently, it has become an emerging pathogen in humans. However, the complete genome sequence of genus Cysteiniphilum is lacking, and the genomic characteristics of genetic diversity, evolutionary dynamics, and pathogenicity have not been characterized. In this study, the complete genome of the first reported clinical isolate QT6929 of genus Cysteiniphilum was sequenced, and comparative genomics analyses to Francisella genus were conducted to unveil the genomic landscape and diversity of the genus Cysteiniphilum. Our results showed that the complete genome of QT6929 consists of one 2.61 Mb chromosome and a 76,819 bp plasmid. The calculated average nucleotide identity and DNA-DNA hybridization values revealed that two clinical isolates QT6929 and JM-1 should be reclassified as two novel species in genus Cysteiniphilum. Pan-genome analysis revealed genomic diversity within the genus Cysteiniphilum and an open pan-genome state. Genomic plasticity analysis exhibited abundant mobile genetic elements including genome islands, insertion sequences, prophages, and plasmids on Cysteiniphilum genomes, which facilitated the broad exchange of genetic material between Cysteiniphilum and other genera like Francisella and Legionella. Several potential virulence genes associated with lipopolysaccharide/lipooligosaccharide, capsule, and haem biosynthesis specific to clinical isolates were predicted and might contribute to their pathogenicity in humans. Incomplete Francisella pathogenicity island was identified in most Cysteiniphilum genomes. Overall, our study provides an updated phylogenomic relationship of members of the genus Cysteiniphilum and comprehensive genomic insights into this rare emerging pathogen.

Comparison of Illumina and Oxford Nanopore Technology for genome analysis of Francisella tularensis, Bacillus anthracis, and Brucella suis.

BACKGROUND: Bacterial epidemiology needs to understand the spread and dissemination of strains in a One Health context. This is important for highly pathogenic bacteria such as Bacillus anthracis, Brucella species, and Francisella tularensis. Whole genome sequencing (WGS) has paved the way for genetic marker detection and high-resolution genotyping. While such tasks are established for Illumina short-read sequencing, Oxford Nanopore Technology (ONT) long-read sequencing has yet to be evaluated for such highly pathogenic bacteria with little genomic variations between strains. In this study, three independent sequencing runs were performed using Illumina, ONT flow cell version 9.4.1, and 10.4 for six strains of each of Ba. anthracis, Br. suis and F. tularensis. Data from ONT sequencing alone, Illumina sequencing alone and two hybrid assembly approaches were compared. RESULTS: As previously shown, ONT produces ultra-long reads, while Illumina produces short reads with higher sequencing accuracy. Flow cell version 10.4 improved sequencing accuracy over version 9.4.1. The correct (sub-)species were inferred from all tested technologies, individually. Moreover, the sets of genetic markers for virulence, were almost identical for the respective species. The long reads of ONT allowed to assemble not only chromosomes of all species to near closure, but also virulence plasmids of Ba. anthracis. Assemblies based on nanopore data alone, Illumina data alone, and both hybrid assemblies correctly detected canonical (sub-)clades for Ba. anthracis and F. tularensis as well as multilocus sequence types for Br. suis. For F. tularensis, high-resolution genotyping using core-genome MLST (cgMLST) and core-genome Single-Nucleotide-Polymorphism (cgSNP) typing produced highly comparable results between data from Illumina and both ONT flow cell versions. For Ba. anthracis, only data from flow cell version 10.4 produced similar results to Illumina for both high-resolution typing methods. However, for Br. suis, high-resolution genotyping yielded larger differences comparing Illumina data to data from both ONT flow cell versions. CONCLUSIONS: In summary, combining data from ONT and Illumina for high-resolution genotyping might be feasible for F. tularensis and Ba. anthracis, but not yet for Br. suis. The ongoing improvement of nanopore technology and subsequent data analysis may facilitate high-resolution genotyping for all bacteria with highly stable genomes in future.