Complete Genome Sequence of a Francisella tularensis subsp. holarctica Strain from Germany Causing Lethal Infection in Common Marmosets

Marmosets as models of infectious diseases

Animal models of infectious disease often serve a crucial purpose in obtaining licensure of therapeutics and medical countermeasures, particularly in situations where human trials are not feasible, i.e., for those diseases that occur infrequently in the human population. The common marmoset (Callithrix jacchus), a Neotropical new-world (platyrrhines) non-human primate, has gained increasing attention as an animal model for a number of diseases given its small size, availability and evolutionary proximity to humans. This review aims to (i) discuss the pros and cons of the common marmoset as an animal model by providing a brief snapshot of how marmosets are currently utilized in biomedical research, (ii) summarize and evaluate relevant aspects of the marmoset immune system to the study of infectious diseases, (iii) provide a historical backdrop, outlining the significance of infectious diseases and the importance of developing reliable animal models to test novel therapeutics, and (iv) provide a summary of infectious diseases for which a marmoset model exists, followed by an in-depth discussion of the marmoset models of two studied bacterial infectious diseases (tularemia and melioidosis) and one viral infectious disease (viral hepatitis C).

Strengthening the genomic surveillance of Francisella tularensis by using culture-free whole-genome sequencing from biological samples

Introduction

Francisella tularensis is a highly infectious bacterium that causes the zoonotic disease tularemia. The development of genotyping methods, especially those based on whole-genome sequencing (WGS), has recently increased the knowledge on the epidemiology of this disease. However, due to the difficulties associated with the growth and isolation of this fastidious pathogen in culture, the availability of strains and subsequently WGS data is still limited.

Methods

To surpass these constraints, we aimed to implement a culture-free approach to capture and sequence F. tularensis genomes directly from complex samples. Biological samples obtained from 50 common voles and 13 Iberian hares collected in Spain were confirmed as positive for F. tularensis subsp. holarctica and subjected to a WGS target capture and enrichment protocol, using RNA oligonucleotide baits designed to cover F. tularensis genomic diversity.

Results

We obtained full genome sequences of F. tularensis from 13 animals (20.6%), two of which had mixed infections with distinct genotypes, and achieved a higher success rate when compared with culture-dependent WGS (only successful for two animals). The new genomes belonged to different clades commonly identified in Europe (B.49, B.51 and B.262) and subclades. Despite being phylogenetically closely related to other genomes from Spain, the detected clusters were often found in other countries. A comprehensive phylogenetic analysis, integrating 599 F. tularensis subsp. holarctica genomes, showed that most (sub)clades are found in both humans and animals and that closely related strains are found in different, and often geographically distant, countries.

Discussion

Overall, we show that the implemented culture-free WGS methodology yields timely, complete and high-quality genomic data of F. tularensis, being a highly valuable approach to promote and potentiate the genomic surveillance of F. tularensis and ultimately increase the knowledge on the genomics, ecology and epidemiology of this highly infectious pathogen.

Genotyping of Francisella tularensis subsp. holarctica from Hares in Germany

Francisella tularensis is the causative agent of the zoonotic disease tularemia. In Germany, most human infections are caused by contact with infected hares. The aim of this study was to characterize Francisella tularensis subsp. holarctica strains isolated from hares in Germany and to develop bioinformatics tools to analyze their genetic relatedness. In total, 257 German isolates-obtained mainly from hares (n = 233), other vertebrate animals, and ticks, but also from humans (n = 3)-were analyzed within this study. Publically available sequence data from 49 isolates were used to put our isolates into an epidemiological context and to compare isolates from natural foci and humans. Whole-genome sequences were analyzed using core-genome Multi-Locus-Sequence-Typing, canonical Single Nucleotide Polymorphism (SNP) typing and whole-genome SNP typing. An overall conformity of genotype clustering between the typing methods was found, albeit with a lower resolution for canonical single SNP typing. The subclade distribution, both on local and national levels, among strains from humans and hares was similar, suggesting circulation of the same genotypes both in animals and humans. Whilst close to identical isolates of the same subclade were found distributed over large areas, small geographical foci often harbored members of different subclades. In conclusion, although genomic high-resolution typing was shown to be robust, reproducible and allowed the identification of highly closely related strains, genetic profiling alone is not always conclusive for epidemiological linkage of F. tularensis strains.

Multiple domains of bacterial and human Lon proteases define substrate selectivity

The Lon protease selectively degrades abnormal proteins or certain normal proteins in response to environmental and cellular conditions in many prokaryotic and eukaryotic organisms. However, the mechanism(s) behind the substrate selection of normal proteins remains largely unknown. In this study, we identified 10 new substrates of F. tularensis Lon from a total of 21 candidate substrates identified in our previous work, the largest number of novel Lon substrates from a single study. Cross-species degradation of these and other known Lon substrates revealed that human Lon is unable to degrade many bacterial Lon substrates, suggestive of a “organism-adapted” substrate selection mechanism for the natural Lon variants. However, individually replacing the N, A, and P domains of human Lon with the counterparts of bacterial Lon did not enable the human protease to degrade the same bacterial Lon substrates. This result showed that the “organism-adapted” substrate selection depends on multiple domains of the Lon proteases. Further in vitro proteolysis and mass spectrometry analysis revealed a similar substrate cleavage pattern between the bacterial and human Lon variants, which was exemplified by predominant representation of leucine, alanine, and other hydrophobic amino acids at the P(−1) site within the substrates. These observations suggest that the Lon proteases select their substrates at least in part by fine structural matching with the proteins in the same organisms.

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

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

Comparison of virulence of Francisella tularensis ssp. holarctica genotypes B.12 and B.FTNF002-00

Background

Two main genetic groups (B.12 and B.FTNF002-00) of Francisella tularensis ssp. holarctica are endemic in Europe. The B.FTNF002-00 group proved to be dominant in Western European countries, while strains of the B.12 group were isolated mainly in Northern, Central and Eastern Europe. The clinical course of tularemia in the European brown hare (Lepus europaeus) also shows distinct patterns according to the geographical area. Acute course of the disease is observed in hares in Western European countries, while signs of sub-acute or chronic infection are more frequently detected in the eastern part of the continent. The aim of the present study was to examine whether there is any difference in the virulence of the strains belonging to the B.FTNF002-00 and B.12 genetic clades.

Results

Experimental infection of Fischer 344 rats was performed by intra-peritoneal injection of three dilutions of a Hungarian (B.12 genotype) and an Italian (B.FTNF002-00 genotype) F. tularensis ssp. holarctica strain. Moderate difference was observed in the virulence of the two genotypes. Significant differences were observed in total weight loss values and scores of clinical signs between the two genotypes with more rats succumbing to tularemia in groups infected with the B.FTNF002-00 genotype.

Conclusions

Results of the experimental infection are consistent with previous clinical observations and pathological studies suggesting that F. tularensis ssp. holarctica genotype B.FTNF002-00 has higher pathogenic potential than the B.12 genotype.

Electronic supplementary material

The online version of this article (doi:10.1186/s12917-017-0968-9) contains supplementary material, which is available to authorized users.

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.

Long-range dispersal moved Francisella tularensis into Western Europe from the East

For many infections transmitting to humans from reservoirs in nature, disease dispersal patterns over space and time are largely unknown. Here, a reversed genomics approach helped us understand disease dispersal and yielded insight into evolution and biological properties of Francisella tularensis, the bacterium causing tularemia. We whole-genome sequenced 67 strains and characterized by single-nucleotide polymorphism assays 138 strains, collected from individuals infected 1947-2012 across Western Europe. We used the data for phylogenetic, population genetic and geographical network analyses. All strains (n=205) belonged to a monophyletic population of recent ancestry not found outside Western Europe. Most strains (n=195) throughout the study area were assigned to a star-like phylogenetic pattern indicating that colonization of Western Europe occurred via clonal expansion. In the East of the study area, strains were more diverse, consistent with a founder population spreading from east to west. The relationship of genetic and geographic distance within the F. tularensis population was complex and indicated multiple long-distance dispersal events. Mutation rate estimates based on year of isolation indicated null rates; in outbreak hotspots only, there was a rate of 0.4 mutations/genome/year. Patterns of nucleotide substitution showed marked AT mutational bias suggestive of genetic drift. These results demonstrate that tularemia has moved from east to west in Europe and that F. tularensis has a biology characterized by long-range geographical dispersal events and mostly slow, but variable, replication rates. The results indicate that mutation-driven evolution, a resting survival phase, genetic drift and long-distance geographical dispersal events have interacted to generate genetic diversity within this species.

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.

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

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

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.

Draft genome sequence of Francisella tularensis subsp. holarctica BD11-00177

Francisella tularensis is a facultative intracellular bacterium in the class Gammaproteobacteria. This strain is of interest because it is the etiologic agent of tularemia and a highly virulent category A biothreat agent. Here we describe the draft genome sequence and annotation of Francisella tularensis subsp. holarctica BD11-00177, isolated from the first case of indigenous tularemia detected in The Netherlands since 1953. Whole genome DNA sequence analysis assigned this isolate to the genomic group B.FTNF002–00, which previously has been exclusively reported from Spain, France, Italy, Switzerland and Germany. Automatic annotation of the 1,813,372 bp draft genome revealed 2,103 protein-coding and 46 RNA genes.