Genomic comparison between a virulent type A1 strain of Francisella tularensis and its attenuated O-antigen mutant

Intraspecies characterization of bacteria via evolutionary modeling of protein domains

The ability to detect and characterize bacteria within a biological sample is crucial for the monitoring of infections and epidemics, as well as for the study of human health and its relationship with commensal microorganisms. To this aim, a commonly used technique is the 16S rRNA gene targeted sequencing. PCR-amplified 16S sequences derived from the sample of interest are usually clustered into the so-called Operational Taxonomic Units (OTUs) based on pairwise similarities. Then, representative OTU sequences are compared with reference (human-made) databases to derive their phylogeny and taxonomic classification. Here, we propose a new reference-free approach to define the phylogenetic distance between bacteria based on protein domains, which are the evolving units of proteins. We extract the protein domain profiles of 3368 bacterial genomes and we use an ecological approach to model their Relative Species Abundance distribution. Based on the model parameters, we then derive a new measurement of phylogenetic distance. Finally, we show that such model-based distance is capable of detecting differences between bacteria in cases in which the 16S rRNA-based method fails, providing a possibly complementary approach , which is particularly promising for the analysis of bacterial populations measured by shotgun sequencing.

HMD-ARG: hierarchical multi-task deep learning for annotating antibiotic resistance genes

BACKGROUND

The spread of antibiotic resistance has become one of the most urgent threats to global health, which is estimated to cause 700,000 deaths each year globally. Its surrogates, antibiotic resistance genes (ARGs), are highly transmittable between food, water, animal, and human to mitigate the efficacy of antibiotics. Accurately identifying ARGs is thus an indispensable step to understanding the ecology, and transmission of ARGs between environmental and human-associated reservoirs. Unfortunately, the previous computational methods for identifying ARGs are mostly based on sequence alignment, which cannot identify novel ARGs, and their applications are limited by currently incomplete knowledge about ARGs.

RESULTS

Here, we propose an end-to-end Hierarchical Multi-task Deep learning framework for ARG annotation (HMD-ARG). Taking raw sequence encoding as input, HMD-ARG can identify, without querying against existing sequence databases, multiple ARG properties simultaneously, including if the input protein sequence is an ARG, and if so, what antibiotic family it is resistant to, what resistant mechanism the ARG takes, and if the ARG is an intrinsic one or acquired one. In addition, if the predicted antibiotic family is beta-lactamase, HMD-ARG further predicts the subclass of beta-lactamase that the ARG is resistant to. Comprehensive experiments, including cross-fold validation, third-party dataset validation in human gut microbiota, wet-experimental functional validation, and structural investigation of predicted conserved sites, demonstrate not only the superior performance of our method over the state-of-art methods, but also the effectiveness and robustness of the proposed method.

CONCLUSIONS

We propose a hierarchical multi-task method, HMD-ARG, which is based on deep learning and can provide detailed annotations of ARGs from three important aspects: resistant antibiotic class, resistant mechanism, and gene mobility. We believe that HMD-ARG can serve as a powerful tool to identify antibiotic resistance genes and, therefore mitigate their global threat. Our method and the constructed database are available at http://www.cbrc.kaust.edu.sa/HMDARG/ . Video abstract (MP4 50984 kb).

Formation of the Francisella tularensis Biofilm is Affected by Cell Surface Glycosylation, Growth Medium, and a Glucan Exopolysaccharide

Biofilms are matrix-associated communities that enable bacteria to colonise environments unsuitable for free-living bacteria. The facultative intracellular pathogen Francisella tularensis can persist in water, amoebae, and arthropods, as well as within mammalian macrophages. F. tularensis Types A and B form poor biofilms, but F. tularensis mutants lacking lipopolysaccharide O-antigen, O-antigen capsule, and capsule-like complex formed up to 15-fold more biofilm than fully glycosylated cells. The Type B live vaccine strain was also 50% less capable of initiating surface attachment than mutants deficient in O-antigen and capsule-like complex. However, the growth medium of all strains tested also influenced the formation of biofilm, which contained a novel exopolysaccharide consisting of an amylose-like glucan. In addition, the surface polysaccharide composition of the bacterium affected the protein:DNA:polysaccharide composition of the biofilm matrix. In contrast, F. novicida attached to surfaces more efficiently and made a more robust biofilm than Type A or B strains, but loss of O-antigen or capsule-like complex did not significantly affect F. novicida biofilm formation. These results indicated that suppression of surface polysaccharides may promote biofilm formation by F. tularensis Types A and B. Whether biofilm formation enhances survival of F. tularensis in aquatic or other environmental niches has yet to be determined.

Further Characterization of the Capsule-Like Complex (CLC) Produced by Francisella tularensis Subspecies tularensis: Protective Efficacy and Similarity to Outer Membrane Vesicles

Francisella tularensis is the etiologic agent of tularemia, and subspecies tularensis (type A) is the most virulent subspecies. The live vaccine strain (LVS) of subspecies holarctica produces a capsule-like complex (CLC) that consists of a large variety of glycoproteins. Expression of the CLC is greatly enhanced when the bacteria are subcultured in and grown on chemically defined medium. Deletion of two genes responsible for CLC glycosylation in LVS results in an attenuated mutant that is protective against respiratory tularemia in a mouse model. We sought to further characterize the CLC composition and to determine if a type A CLC glycosylation mutant would be attenuated in mice. The CLCs isolated from LVS extracted with 0.5% phenol or 1 M urea were similar, as determined by gel electrophoresis and Western blotting, but the CLC extracted with urea was more water-soluble. The CLC extracted with either 0.5% phenol or 1 M urea from type A strains was also similar to the CLC of LVS in antigenic properties, electrophoretic profile, and by transmission electron microscopy (TEM). The solubility of the CLC could be further enhanced by fractionation with Triton X-114 followed by N-Lauroylsarcosine detergents; the largest (>250 kDa) molecular size component appeared to be an aggregate of smaller components. Outer membrane vesicles/tubules (OMV/T) isolated by differential centrifugation and micro-filtration appeared similar to the CLC by TEM, and many of the proteins present in the OMV/T were also identified in soluble and insoluble fractions of the CLC. Further investigation is warranted to assess the relationship between OMV/T and the CLC. The CLC conjugated to keyhole limpet hemocyanin or flagellin was highly protective against high-dose LVS intradermal challenge and partially protective against intranasal challenge. A protective response was associated with a significant rise in cytokines IL-12, IL-10, and IFN-γ. However, a type A CLC glycosylation mutant remained virulent in BALB/c mice, and immunization with the CLC did not protect mice against high dose respiratory challenge with type A strain SCHU S4.

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.

Glycosylation of a Capsule-Like Complex (CLC) by Francisella novicida Is Required for Virulence and Partial Protective Immunity in Mice

Francisella tularensis is a Gram-negative bacterium and the etiologic agent of tularemia. F. tularensis may appear encapsulated when examined by transmission electron microscopy (TEM), which is due to production of an extracellular capsule-like complex (CLC) when the bacterium is grown under specific environmental conditions. Deletion of two glycosylation genes in the live vaccine strain (LVS) results in loss of apparent CLC and attenuation of LVS in mice. In contrast, F. novicida, which is also highly virulent for mice, is reported to be non-encapsulated. However, the F. novicida genome contains a putative polysaccharide locus with homology to the CLC glycosylation locus in F. tularensis. Following daily subculture of F. novicida in Chamberlain's defined medium, an electron dense material surrounding F. novicida, similar to the F. tularensis CLC, was evident. Extraction with urea effectively removed the CLC, and compositional analysis indicated the extract contained galactose, glucose, mannose, and multiple proteins, similar to those found in the F. tularensis CLC. The same glycosylation genes deleted in LVS were targeted for deletion in F. novicida by allelic exchange using the same mutagenesis vector used for mutagenesis of LVS. In contrast, this mutation also resulted in the loss of five additional genes immediately upstream of the targeted mutation (all within the glycosylation locus), resulting in strain F. novicida Δ1212-1218. The subcultured mutant F. novicida Δ1212-1218 was CLC-deficient and the CLC contained significantly less carbohydrate than the subcultured parent strain. The mutant was severely attenuated in BALB/c mice inoculated intranasally, as determined by the lower number of F. novicida Δ1212-1218 recovered in tissues compared to the parent, and by clearance of the mutant by 10-14 days post-challenge. Mice immunized intranasally with F. novicida Δ1212-1218 were partially protected against challenge with the parent, produced significantly reduced levels of inflammatory cytokines, and their spleens contained only areas of lymphoid hyperplasia, whereas control mice challenged with the parent exhibited hypercytokinemia and splenic necrosis. Therefore, F. novicida is capable of producing a CLC similar to that of F. tularensis, and glycosylation of the CLC contributed to F. novicida virulence and immunoprotection.

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.

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

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

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.

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.