Genes under positive selection in the core genome of pathogenic Bacillus cereus group members

A Comparative Analysis of the Core Proteomes within and among the Bacillus subtilis and Bacillus cereus Evolutionary Groups Reveals the Patterns of Lineage- and Species-Specific Adaptations

By integrating phylogenomic and comparative analyses of 1104 high-quality genome sequences, we identify the core proteins and the lineage-specific fingerprint proteins of the various evolutionary clusters (clades/groups/species) of the Bacillus genus. As fingerprints, we denote those core proteins of a certain lineage that are present only in that particular lineage and absent in any other Bacillus lineage. Thus, these lineage-specific fingerprints are expected to be involved in particular adaptations of that lineage. Intriguingly, with a few notable exceptions, the majority of the Bacillus species demonstrate a rather low number of species-specific fingerprints, with the majority of them being of unknown function. Therefore, species-specific adaptations are mostly attributed to highly unstable (in evolutionary terms) accessory proteomes and possibly to changes at the gene regulation level. A series of comparative analyses consistently demonstrated that the progenitor of the Cereus Clade underwent an extensive genomic expansion of chromosomal protein-coding genes. In addition, the majority (76–82%) of the B. subtilis proteins that are essential or play a significant role in sporulation have close homologs in most species of both the Subtilis and the Cereus Clades. Finally, the identification of lineage-specific fingerprints by this study may allow for the future development of highly specific vaccines, therapeutic molecules, or rapid and low-cost molecular tests for species identification.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan‐genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well‐studied Bacillus cereus sensu lato group, we identified clade‐specific core genes (present in all clade members) and then used clade‐specific allelic diversity to identify genes under purifying and diversifying selection. Clade‐specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade‐specific core or accessory genomes. Genes under purifying selection were often clade‐specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade‐specific selection, allowing ready identification of strains and their ecological niche.

In-depth genome analysis of Bacillus sp. BH32, a salt stress-tolerant endophyte obtained from a halophyte in a semiarid region

Endophytic strains belonging to the Bacillus cereus group were isolated from the halophytes Atriplex halimus L. (Amaranthaceae) and Tamarix aphylla L. (Tamaricaceae) from costal and continental regions in Algeria. Based on their salt tolerance (up to 5%), the strains were tested for their ability to alleviate salt stress in tomato and wheat. Bacillus sp. strain BH32 showed the highest potential to reduce salinity stress (up to + 50% and + 58% of dry weight improvement, in tomato and wheat, respectively, compared to the control). To determine putative mechanisms involved in salt tolerance and plant growth promotion, the whole genome of Bacillus sp. BH32 was sequenced, annotated, and used for comparative genomics against the genomes of closely related strains. The pangenome of Bacillus sp. BH32 and its closest relative was further analyzed. The phylogenomic analyses confirmed its taxonomic position, a member of the Bacillus cereus group, with intergenomic distances (GBDP analysis) pinpointing to a new taxon (digital DNA-DNA hybridization, dDDH < 70%). Genome mining unveiled several genes involved in stress tolerance, production of anti-oxidants and genes involved in plant growth promotion as well as in the production of secondary metabolites. KEY POINTS : • Bacillus sp. BH32 and other bacterial endophytes were isolated from halophytes, to be tested on tomato and wheat and to limit salt stress adverse effects. • The strain with the highest potential was then studied at the genomic level to highlight numerous genes linked to plant growth promotion and stress tolerance. • Pangenome approaches suggest that the strain belongs to a new taxon within the Bacillus cereus group.

The Distribution of Several Genomic Virulence Determinants Does Not Corroborate the Established Serotyping Classification of Bacillus thuringiensis

Bacillus thuringiensis, commonly referred to as Bt, is an object of the lasting interest of microbiologists due to its highly effective insecticidal properties, which make Bt a prominent source of biologicals. To categorize the exuberance of Bt strains discovered, serotyping assays are utilized in which flagellin serves as a primary seroreactive molecule. Despite its convenience, this approach is not indicative of Bt strains' phenotypes, neither it reflects actual phylogenetic relationships within the species. In this respect, comparative genomic and proteomic techniques appear more informative, but their use in Bt strain classification remains limited. In the present work, we used a bottom-up proteomic approach based on fluorescent two-dimensional difference gel electrophoresis (2D-DIGE) coupled with liquid chromatography/tandem mass spectrometry(LC-MS/MS) protein identification to assess which stage of Bt culture, vegetative or spore, would be more informative for strain characterization. To this end, the proteomic differences for the israelensis-attributed strains were assessed to compare sporulating cultures of the virulent derivative to the avirulent one as well as to the vegetative stage virulent bacteria. Using the same approach, virulent spores of the israelensis strain were also compared to the spores of strains belonging to two other major Bt serovars, namely darmstadiensis and thuringiensis. The identified proteins were analyzed regarding the presence of the respective genes in the 104 Bt genome assemblies available at open access with serovar attributions specified. Of 21 proteins identified, 15 were found to be encoded in all the present assemblies at 67% identity threshold, including several virulence factors. Notable, individual phylogenies of these core genes conferred neither the serotyping nor the flagellin-based phylogeny but corroborated the reconstruction based on phylogenomics approaches in terms of tree topology similarity. In its turn, the distribution of accessory protein genes was not confined to the existing serovars. The obtained results indicate that neither gene presence nor the core gene sequence may serve as distinctive bases for the serovar attribution, undermining the notion that the serotyping system reflects strains' phenotypic or genetic similarity. We also provide a set of loci, which fit in with the phylogenomics data plausibly and thus may serve for draft phylogeny estimation of the novel strains.

The Food Poisoning Toxins of Bacillus cereus

Bacillus cereus is a ubiquitous soil bacterium responsible for two types of food-associated gastrointestinal diseases. While the emetic type, a food intoxication, manifests in nausea and vomiting, food infections with enteropathogenic strains cause diarrhea and abdominal pain. Causative toxins are the cyclic dodecadepsipeptide cereulide, and the proteinaceous enterotoxins hemolysin BL (Hbl), nonhemolytic enterotoxin (Nhe) and cytotoxin K (CytK), respectively. This review covers the current knowledge on distribution and genetic organization of the toxin genes, as well as mechanisms of enterotoxin gene regulation and toxin secretion. In this context, the exceptionally high variability of toxin production between single strains is highlighted. In addition, the mode of action of the pore-forming enterotoxins and their effect on target cells is described in detail. The main focus of this review are the two tripartite enterotoxin complexes Hbl and Nhe, but the latest findings on cereulide and CytK are also presented, as well as methods for toxin detection, and the contribution of further putative virulence factors to the diarrheal disease.

Identification of bacterial endospores and targeted detection of foodborne viruses in industrially reared insects for food

With edible insects being increasingly produced, food safety authorities have called for the determination of microbiological challenges posed to human health. Here, we find that the bacterial endospore fraction in industrially reared mealworm and cricket samples is largely comprised of Bacillus cereus group members that can pose insect or human health risks. Hepatitis A virus, hepatitis E virus and norovirus genogroup II were not detected in the sample collection, indicating a low food safety risk from these viral pathogens.

Genetic diversity estimation of Yunnan indigenous goat breeds using microsatellite markers

BACKGROUND

To assess the genetic diversity of seven Yunnan indigenous goat populations (Fengqing hornless goat, Mile red-bone goat, Longling goat, Ninglang black goat, Black-bone goat, Yunling black goat, and Zhaotong goat), their population structures were investigated using 20 microsatellite markers.

RESULTS

The results indicated that the genetic diversity of these goats was rich. The observed heterozygosity ranged from 0.4667 ± 0.0243 to 0.5793 ± 0.0230, and the mean number of alleles ranged from 4.80 ± 1.61 and 4.80 ± 1.64 to 6.20 ± 2.93. The population structure analysis showed that these seven goat populations were separated into two clusters, consistent with the results from phylogenetic networks, pairwise differences, and STRUCTURE analyses. We speculate that this may have been caused by natural geographical isolation, human migration and economic and cultural exchanges. We suggest removing CSRD247 and ILSTS005, two loci identified to be under positive selection in the present study, from the microsatellite evaluation system of goats.

CONCLUSIONS

The present study may provide a scientific basis for the conservation and utilization of Yunnan indigenous goats.

Genome-wide analyses reveal genes subject to positive selection in Toxoplasma gondii

Toxoplasma gondii is an important protozoan pathogen that infects many wild and domestic animals and causes infections in immunocompromised humans. However, there has been little investigation of the molecular evolutionary trajectories of this pathogenic protozoa using comparative genomics data. Here, we employed a comparative evolutionary genomics approach to identify genes that are under site- and lineage-specific positive selection in nine strains of T. gondii, including two closely related species, Neospora caninum and Hammondia hammondi. Based on the analyses of five coccidian core genomes, 4.5% of the 5788 core genome genes showed strong signals for positive selection in the site model. In addition, the branch-site model analyses in the nine T. gondii core genomes indicated that 2 to 20 genes underwent significant positive selection along each lineage leading to T. gondii strains. Many of the protein products encoded by the positively selected genes are secretory or surface proteins that have previously been implicated in host pathogenesis. The adaptive changes in these positively selected genes might be related to dynamic interactions between the host immune systems and might play a crucial role in the infection and pathogenic processes of T. gondii.

Whole-genome-based phylogeny of Bacillus cytotoxicus reveals different clades within the species and provides clues on ecology and evolution

Bacillus cytotoxicus is a member of the Bacillus cereus group linked to fatal cases of diarrheal disease. Information on B. cytotoxicus is very limited; in particular comprehensive genomic data is lacking. Thus, we applied a genomic approach to characterize B. cytotoxicus and decipher its population structure. To this end, complete genomes of ten B. cytotoxicus were sequenced and compared to the four publicly available full B. cytotoxicus genomes and genomes of other B. cereus group members. Average nucleotide identity, core genome, and pan genome clustering resulted in clear distinction of B. cytotoxicus strains from other strains of the B. cereus group. Genomic content analyses showed that a hydroxyphenylalanine operon is present in B. cytotoxicus, but absent in all other members of the B. cereus group. It enables degradation of aromatic compounds to succinate and pyruvate and was likely acquired from another Bacillus species. It allows for utilization of tyrosine and might have given a B. cytotoxicus ancestor an evolutionary advantage resulting in species differentiation. Plasmid content showed that B. cytotoxicus is flexible in exchanging genes, allowing for quick adaptation to the environment. Genome-based phylogenetic analyses divided the B. cytotoxicus strains into four clades that also differed in virulence gene content.