Genomic analysis of three African strains of Bacillus anthracis demonstrates that they are part of the clonal expansion of an exclusively pathogenic bacterium

Molecular characterization of an outbreak-involved Bacillus anthracis strain confirms the spillover of anthrax from West Africa

BACKGROUND

Anthrax, a zoonotic disease caused by the spore-forming bacterium Bacillus anthracis, remains a major global public health concern, especially in countries with limited resources. Sierra Leone, a West African country historically plagued by anthrax, has almost been out of report on this disease in recent decades. In this study, we described a large-scale anthrax outbreak affecting both animals and humans and attempted to characterize the pathogen using molecular techniques.

METHODS

The causative agent of the animal outbreak in Port Loko District, Sierra Leone, between March and May 2022 was identified using the nanopore sequencing technique. A nationwide active surveillance was implemented from May 2022 to June 2023 to monitor the occurrence of anthrax-specific symptoms in humans. Suspected cases were subsequently verified using quantitative polymerase chain reaction. Full-genome sequencing was accomplished by combining long-read and short-read sequencing methods. Subsequent phylogenetic analysis was performed based on the full-chromosome single nucleotide polymorphisms.

RESULTS

The outbreak in Port Loko District, Sierra Leone, led to the death of 233 animals between March 26th and May 16th, 2022. We ruled out the initial suspicion of Anaplasma species and successfully identified B. anthracis as the causative agent of the outbreak. As a result of the government's prompt response, out of the 49 suspected human cases identified during the one-year active surveillance, only 6 human cases tested positive, all within the first month after the official declaration of the outbreak. The phylogenetic analysis indicated that the BaSL2022 isolate responsible for the outbreak was positioned in the A.Br.153 clade within the TransEuroAsian group of B. anthracis.

CONCLUSIONS

We successfully identified a large-scale anthrax outbreak in Sierra Leone. The causative isolate of B. anthracis, BaSL2022, phylogenetically bridged other lineages in A.Br.153 clade and neighboring genetic groups, A.Br.144 and A.Br.148, eventually confirming the spillover of anthrax from West Africa. Given the wide dissemination of B. anthracis spores, it is highly advisable to effectively monitor the potential reoccurrence of anthrax outbreaks and to launch campaigns to improve public awareness regarding anthrax in Sierra Leone.

Effects of Bacillus coagulans as an adjunct starter culture on yogurt quality and storage

Bacillus coagulans has been widely studied for its probiotic properties. Therefore, identifying a strain that can be used as an adjunct starter culture for yogurt production would have commercial value. In this study, 30 B. coagulans strains were isolated from vegetable samples from 11 provinces or autonomous regions in China, and their pan-genomic and phylogenetic characteristics were analyzed. Phylogenetic analysis categorized 30 strains into 4 different subphylotypes, including subtype I (11 isolates), subtype II (7 isolates), subtype III (11 isolates), and subtype IV (1 isolate). Four B. coagulans strains (B. coagulans-70, B. coagulans-78, B. coagulans-79, and B. coagulans-100) were randomly selected from each subphylotype of the phylogenetic tree as adjunct starter cultures. Compared with the other tested strains, B. coagulans-70 showed the highest count in yogurt at the end of the manufacturing period. Comparative genome analysis indicated that the different bacterial levels of B. coagulans strains in yogurt may be associated with the abundance of genes related to carbohydrate transport and metabolism (e.g., sucrose utilization). Finally, differences in texture and volatile flavor compound profiles were observed between the yogurt samples. Compared with the other groups, the addition of B. coagulans-70 exerted a positive effect on the appearance and texture of yogurt products. Volatile analysis showed increased quantities of 2-heptanone, 2-nonanone, amyl alcohol, and 2-hydroxy-3-pentanone in the B. coagulans-70 group compared with control yogurts. These results above combined with the results of a sensory evaluation indicated that B. coagulans-70 is the most suitable strain for further use in functional dairy product development.

Comparative genomics and pangenome-oriented studies reveal high homogeneity of the agronomically relevant enterobacterial plant pathogen Dickeya solani

Background

Dickeya solani is an important plant pathogenic bacterium causing severe losses in European potato production. This species draws a lot of attention due to its remarkable virulence, great devastating potential and easier spread in contrast to other Dickeya spp. In view of a high need for extensive studies on economically important soft rot Pectobacteriaceae, we performed a comparative genomics analysis on D. solani strains to search for genetic foundations that would explain the differences in the observed virulence levels within the D. solani population.

Results

High quality assemblies of 8 de novo sequenced D. solani genomes have been obtained. Whole-sequence comparison, ANIb, ANIm, Tetra and pangenome-oriented analyses performed on these genomes and the sequences of 14 additional strains revealed an exceptionally high level of homogeneity among the studied genetic material of D. solani strains. With the use of 22 genomes, the pangenome of D. solani, comprising 84.7% core, 7.2% accessory and 8.1% unique genes, has been almost completely determined, suggesting the presence of a nearly closed pangenome structure. Attribution of the genes included in the D. solani pangenome fractions to functional COG categories showed that higher percentages of accessory and unique pangenome parts in contrast to the core section are encountered in phage/mobile elements- and transcription- associated groups with the genome of RNS 05.1.2A strain having the most significant impact. Also, the first D. solani large-scale genome-wide phylogeny computed on concatenated core gene alignments is herein reported.

Conclusions

The almost closed status of D. solani pangenome achieved in this work points to the fact that the unique gene pool of this species should no longer expand. Such a feature is characteristic of taxa whose representatives either occupy isolated ecological niches or lack efficient mechanisms for gene exchange and recombination, which seems rational concerning a strictly pathogenic species with clonal population structure. Finally, no obvious correlations between the geographical origin of D. solani strains and their phylogeny were found, which might reflect the specificity of the international seed potato market.

Interest of bacterial pangenome analyses in clinical microbiology

Thanks to the progress and decreasing costs in genome sequencing technologies, more than 250,000 bacterial genomes are currently available in public databases, covering most, if not all, of the major human-associated phylogenetic groups of these microorganisms, pathogenic or not. In addition, for many of them, sequences from several strains of a given species are available, thus enabling to evaluate their genetic diversity and study their evolution. In addition, the significant cost reduction of bacterial whole genome sequencing as well as the rapid increase in the number of available bacterial genomes have prompted the development of pangenomic software tools. The study of bacterial pangenome has many applications in clinical microbiology. It can unveil the pathogenic potential and ability of bacteria to resist antimicrobials as well identify specific sequences and predict antigenic epitopes that allow molecular or serologic assays and vaccines to be designed. Bacterial pangenome constitutes a powerful method for understanding the history of human bacteria and relating these findings to diagnosis in clinical microbiology laboratories in order to optimize patient management.

Bioinformatics for Marine Products: An Overview of Resources, Bottlenecks, and Perspectives

The sea represents a major source of biodiversity. It exhibits many different ecosystems in a huge variety of environmental conditions where marine organisms have evolved with extensive diversification of structures and functions, making the marine environment a treasure trove of molecules with potential for biotechnological applications and innovation in many different areas. Rapid progress of the omics sciences has revealed novel opportunities to advance the knowledge of biological systems, paving the way for an unprecedented revolution in the field and expanding marine research from model organisms to an increasing number of marine species. Multi-level approaches based on molecular investigations at genomic, metagenomic, transcriptomic, metatranscriptomic, proteomic, and metabolomic levels are essential to discover marine resources and further explore key molecular processes involved in their production and action. As a consequence, omics approaches, accompanied by the associated bioinformatic resources and computational tools for molecular analyses and modeling, are boosting the rapid advancement of biotechnologies. In this review, we provide an overview of the most relevant bioinformatic resources and major approaches, highlighting perspectives and bottlenecks for an appropriate exploitation of these opportunities for biotechnology applications from marine resources.

Virulence Plasmids of the Pathogenic Clostridia

The clostridia cause a spectrum of diseases in humans and animals ranging from life-threatening tetanus and botulism, uterine infections, histotoxic infections and enteric diseases, including antibiotic-associated diarrhea, and food poisoning. The symptoms of all these diseases are the result of potent protein toxins produced by these organisms. These toxins are diverse, ranging from a multitude of pore-forming toxins to phospholipases, metalloproteases, ADP-ribosyltransferases and large glycosyltransferases. The location of the toxin genes is the unifying theme of this review because with one or two exceptions they are all located on plasmids or on bacteriophage that replicate using a plasmid-like intermediate. Some of these plasmids are distantly related whilst others share little or no similarity. Many of these toxin plasmids have been shown to be conjugative. The mobile nature of these toxin genes gives a ready explanation of how clostridial toxin genes have been so widely disseminated both within the clostridial genera as well as in the wider bacterial community.

Clustering of Pan- and Core-genome of Lactobacillus provides Novel Evolutionary Insights for Differentiation

Background

Bacterial taxonomy aims to classify bacteria based on true evolutionary events and relies on a polyphasic approach that includes phenotypic, genotypic and chemotaxonomic analyses. Until now, complete genomes are largely ignored in taxonomy. The genus Lactobacillus consists of 173 species and many genomes are available to study taxonomy and evolutionary events.

Results

We analyzed and clustered 98 completely sequenced genomes of the genus Lactobacillus and 234 draft genomes of 5 different Lactobacillus species, i.e. L. reuteri, L. delbrueckii, L. plantarum, L. rhamnosus and L. helveticus. The core-genome of the genus Lactobacillus contains 266 genes and the pan-genome 20′800 genes. Clustering of the Lactobacillus pan- and core-genome resulted in two highly similar trees. This shows that evolutionary history is traceable in the core-genome and that clustering of the core-genome is sufficient to explore relationships. Clustering of core- and pan-genomes at species’ level resulted in similar trees as well. Detailed analyses of the core-genomes showed that the functional class “genetic information processing” is conserved in the core-genome but that “signaling and cellular processes” is not. The latter class encodes functions that are involved in environmental interactions. Evolution of lactobacilli seems therefore directed by the environment. The type species L. delbrueckii was analyzed in detail and its pan-genome based tree contained two major clades whose members contained different genes yet identical functions. In addition, evidence for horizontal gene transfer between strains of L. delbrueckii, L. plantarum, and L. rhamnosus, and between species of the genus Lactobacillus is presented. Our data provide evidence for evolution of some lactobacilli according to a parapatric-like model for species differentiation.

Conclusions

Core-genome trees are useful to detect evolutionary relationships in lactobacilli and might be useful in taxonomic analyses. Lactobacillus’ evolution is directed by the environment and HGT.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4601-5) contains supplementary material, which is available to authorized users.

Comparative Genomics for Prokaryotes

Bacteria and archaea, collectively known as prokaryotes, have in general genomes that are much smaller than those of eukaryotes. As a result, thousands of these genomes have been sequenced. In prokaryotes, gene architecture lacks the intron-exon structure of eukaryotic genes (with an occasional exception). These two facts mean that there is an abundance of data for prokaryotic genomes, and that they are easier to study than the more complex eukaryotic genomes. In this chapter, we provide an overview of genome comparison tools that have been developed primarily (sometimes exclusively) for prokaryotic genomes. We cover methods that use only the DNA sequences, methods that use only the gene content, and methods that use both data types.

Comparison of French and Worldwide Bacillus anthracis Strains Favors a Recent, Post-Columbian Origin of the Predominant North-American Clade

Background

Bacillus anthracis, the highly dangerous zoonotic bacterial pathogen species is currently composed of three genetic groups, called A, B and C. Group A is represented worldwide whereas group B is present essentially in Western Europe and Southern Africa. Only three strains from group C have been reported. This knowledge is derived from the genotyping of more than 2000 strains collected worldwide. Strains from both group A and group B are present in France. Previous investigations showed that the majority of sporadic French strains belong to the so-called A.Br.011/009 group A clade and define a very remarkable polytomy with six branches. Here we explore the significance of this polytomy by comparing the French B. anthracis lineages to worldwide lineages. We take advantage of whole genome sequence data previously determined for 122 French strains and 45 strains of various origins.

Results

A total of 6690 SNPs was identified among the available dataset and used to draw the phylogeny. The phylogeny of the French B group strains which belongs to B.Br.CNEVA indicates an expansion from the south-east part of France (the Alps) towards the south-west (Massif-Central and Pyrenees). The relatively small group A strains belonging to A.Br.001/002 results from at least two independent introductions. Strikingly, the data clearly demonstrates that the currently predominant B. anthracis lineage in North America, called WNA for Western North American, is derived from one branch of the A.Br.011/009 polytomy predominant in France.

Conclusions/Significance

The present work extends the range of observed substitution rate heterogeneity within B. anthracis, in agreement with its ecology and in contrast with some other pathogens. The population structure of the six branches A.Br.011/009 polytomy identified in France, diversity of branch length, and comparison with the WNA lineage, suggests that WNA is of post-Columbian and west European origin, with France as a likely source. Furthermore, it is tempting to speculate that the polytomy’s most recent common ancestor -MRCA- dates back to the Hundred Years' war between France and England started in the mid-fourteenth century. These events were associated in France with deadly epidemics and major economic and social changes.

Whole-genome sequencing investigation of animal-skin-drum-associated UK anthrax cases reveals evidence of mixed populations and relatedness to a US case

There have been two anthrax cases affecting people that played and/or made animal-skin drums in the UK during the last 10 years, with single fatal occurrences in Scotland in 2006 and London in 2008. Investigations by the Health Protection Agency (now Public Health England) employing multi-locus-variable number tandem repeat analysis had previously linked the clinical cases to spores associated with animal skins and drums the patients had been in contact with. In this study, whole-genome sequencing of 23 Bacillus anthracis isolates harvested during the investigations was performed. High-quality draft assemblies of these genomes provided greater characterization of the B. anthracis strains present and placed them all upon a new branch of the global phylogeny. Although closely related, the clinical isolates from the two events, and another isolated from a drum-skin-associated case in New York in 2006, were distinct from each other. Multiple distinct genotypes were found during both investigations, implying either multiple contamination events or a single heterogeneous contamination. One environmental isolate from the Scottish incident was more closely related to London isolates than to the other Scottish isolates. As B. anthracis of this subgroup was present at both geographically and temporally distinct events, it may be more widespread at the source of contamination. All isolates were distinct from currently characterized West African strains, despite this being the likely origin of the drums and hides, therefore adding to our knowledge of B. anthracis diversity in the region.

Pan-genomic analysis to redefine species and subspecies based on quantum discontinuous variation: the Klebsiella paradigm

Background

Various methods are currently used to define species and are based on the phylogenetic marker 16S ribosomal RNA gene sequence, DNA-DNA hybridization and DNA GC content. However, these are restricted genetic tools and showed significant limitations.

Results

In this work, we describe an alternative method to build taxonomy by analyzing the pan-genome composition of different species of the Klebsiella genus. Klebsiella species are Gram-negative bacilli belonging to the large Enterobacteriaceae family. Interestingly, when comparing the core/pan-genome ratio; we found a clear discontinuous variation that can define a new species.

Conclusions

Using this pan-genomic approach, we showed that Klebsiella pneumoniae subsp. ozaenae and Klebsiella pneumoniae subsp. rhinoscleromatis are species of the Klebsiella genus, rather than subspecies of Klebsiella pneumoniae. This pan-genomic analysis, helped to develop a new tool for defining species introducing a quantic perspective for taxonomy.

Reviewers

This article was reviewed by William Martin, Pierre Pontarotti and Pere Puigbo (nominated by Dr Yuri Wolf).

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-015-0085-2) contains supplementary material, which is available to authorized users.