A pangenomic study of Bacillus thuringiensis

Bacillus thuringiensis Bt_UNVM-84, a Novel Strain Showing Insecticidal Activity against Anthonomus grandis Boheman (Coleoptera: Curculionidae)

Bacillus thuringiensis is a Gram-positive bacterium known for its insecticidal proteins effective against various insect pests. However, limited strains and proteins target coleopteran pests like Anthonomous grandis Boheman, causing substantial economic losses in the cotton industry. This study focuses on characterizing a Bacillus sp. strain, isolated from Oncativo (Argentina), which exhibits ovoid to amorphous parasporal crystals and was designated Bt_UNVM-84. Its genome encodes genes for the production of two pairs of binary Vpb1/Vpa2 proteins and three Cry-like proteins showing similarity with different Cry8 proteins. Interestingly, this gene content was found to be conserved in a previously characterized Argentine isolate of B. thuringiensis designated INTA Fr7-4. SDS-PAGE analysis revealed a major band of 130 kDa that is proteolytically processed to an approximately 66-kDa protein fragment by trypsin. Bioassays performed with spore-crystal mixtures demonstrated an interesting insecticidal activity against the cotton boll weevil A. grandis neonate larvae, resulting in 91% mortality. Strain Bt_UNVM-84 is, therefore, an interesting candidate for the efficient biological control of this species, causing significant economic losses in the cotton industry in the Americas.

More than just an insect killer: The non-insecticidal activities of Bacillus thuringiensis with biotechnological potential

Bacillus thuringiensis (Bt) is known for the biological control of important insect pests, but scientific advances have revealed several interesting characteristics, in addition to this classical function as a bioinsecticide. To investigate the current knowledge about these non-insecticidal activities, a systematic research on primary data in the scientific literature was conducted on alternative functions of Bt with biotechnological potential. Out of a total of 140 articles selected, 15 non-insecticidal Bt activities were found. Publications related to this topic are available since 1971, and different metadata were reported, such as biomolecules and genes involved in Bt performances in non-insecticidal bioactivities. A total of 11 Bt activities with different effect measures (response variables) were identified, with an average of 48 distinct Bt strains evaluated per activity. Approximately 81.2% of all identified experiments/tests deal with the direct effects of Bt on target cells/organisms, with 36.3% of the strains within these studies tested for antibacterial action; of all microbial targets tested, 92.8% are bacteria, which led to 75.2% of the experimental conditions for all direct activities being performed in vitro. Regarding indirect Bt activities, 67.6% of these studies reported tritrophic Bt-plant-pathogen interactions. Bioremediation also appears as a relevant Bt activity being investigated in-depth. Alternative Bt activities offer innovative ways of developing biotechnology for different areas of anthropic interest; hence, we also focus on the possibility of finding multifunctional strains of Bt, as this may be advantageous from a bioeconomic point of view. Our findings are discussed in terms of research trends, aspects, details and depth of the current knowledge on alternative non-insecticidal Bt traits. We also discuss the potential application of this science for useful technological developments, aiming at solving issues related to human health, sustainable agriculture and environmental preservation/restoration.

mOTUpan: a robust Bayesian approach to leverage metagenome-assembled genomes for core-genome estimation

Recent advances in sequencing and bioinformatics have expanded the tree of life by providing genomes for uncultured environmentally relevant clades, either through metagenome-assembled genomes or through single-cell genomes. While this expanded diversity can provide novel insights into microbial population structure, most tools available for core-genome estimation are sensitive to genome completeness. Consequently, a major portion of the huge phylogenetic diversity uncovered by environmental genomic approaches remains excluded from such analyses. We present mOTUpan, a novel iterative Bayesian method for computing the core genome for sets of genomes of highly diverse completeness range. The likelihood for each gene cluster to belong to core or accessory genome is estimated by computing the probability of its presence/absence pattern in the target genome set. The core-genome prediction is computationally efficient and can be scaled up to thousands of genomes. It has shown comparable estimates to state-of-the-art tools Roary and PPanGGOLiN for high-quality genomes and is capable of using genomes at lower completeness thresholds. mOTUpan wraps a bootstrapping procedure to estimate the quality of a specific core-genome prediction, as the accuracy of each run will depend on the specific completeness distribution and the number of genomes in the dataset under scrutiny. mOTUpan is implemented in the mOTUlizer software package, and available at github.com/moritzbuck/mOTUlizer, under GPL 3.0 license.

Bacillus pumilus Group Comparative Genomics: Toward Pangenome Features, Diversity, and Marine Environmental Adaptation

Background

Members of the Bacillus pumilus group (abbreviated as the Bp group) are quite diverse and ubiquitous in marine environments, but little is known about correlation with their terrestrial counterparts. In this study, 16 marine strains that we had isolated before were sequenced and comparative genome analyses were performed with a total of 52 Bp group strains. The analyses included 20 marine isolates (which included the 16 new strains) and 32 terrestrial isolates, and their evolutionary relationships, differentiation, and environmental adaptation.

Results

Phylogenomic analysis revealed that the marine Bp group strains were grouped into three species: B. pumilus, B. altitudinis and B. safensis. All the three share a common ancestor. However, members of B. altitudinis were observed to cluster independently, separating from the other two, thus diverging from the others. Consistent with the universal nature of genes involved in the functioning of the translational machinery, the genes related to translation were enriched in the core genome. Functional genomic analyses revealed that the marine-derived and the terrestrial strains showed differences in certain hypothetical proteins, transcriptional regulators, K⁺ transporter (TrK) and ABC transporters. However, species differences showed the precedence of environmental adaptation discrepancies. In each species, land specific genes were found with possible functions that likely facilitate survival in diverse terrestrial niches, while marine bacteria were enriched with genes of unknown functions and those related to transcription, phage defense, DNA recombination and repair.

Conclusion

Our results indicated that the Bp isolates show distinct genomic features even as they share a common core. The marine and land isolates did not evolve independently; the transition between marine and non-marine habitats might have occurred multiple times. The lineage exhibited a priority effect over the niche in driving their dispersal. Certain intra-species niche specific genes could be related to a strains adaptation to its respective marine or terrestrial environment(s). In summary, this report describes the systematic evolution of 52 Bp group strains and will facilitate future studies toward understanding their ecological role and adaptation to marine and/or terrestrial environments.

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.

Comparative genomic analysis and mosquito larvicidal activity of four Bacillus thuringiensis serovar israelensis strains

Bacillus thuringiensis serovar israelensis (Bti) is used to control insect vectors of human and animal diseases. In the present study, the toxicity of four strains of Bti, named T0124, T0131, T0137, and T0139, toward Aedes aegypti and Culex quinquefasciatus larvae was analyzed. The T0131 strain showed the highest larvicidal activity against A. aegypti (LC₅₀ = 0.015 µg/ml) and C. quinquefasciatus larvae (LC₅₀ = 0.035 µg/ml) when compared to the other strains. Furthermore, the genomic sequences of the four strains were obtained and compared. These Bti strains had chromosomes sizes of approximately 5.4 Mb with GC contents of ~35% and 5472–5477 putative coding regions. Three small plasmids (5.4, 6.8, and 7.6 kb) and three large plasmids (127, 235, and 359 kb) were found in the extrachromosomal content of all four strains. The SNP-based phylogeny revealed close relationship among isolates from this study and other Bti isolates, and SNPs analysis of the plasmids 127 kb did not reveal any mutations in δ-endotoxins genes. This newly acquired sequence data for these Bti strains may be useful in the search for novel insecticidal toxins to improve existing ones or develop new strategies for the biological control of important insect vectors of human and animal diseases.

Pan- and core- gene association networks: Integrative approaches to understanding biological regulation

The rapid increase in transcriptome data provides an opportunity to access the complex regulatory mechanisms in cellular systems through gene association network (GAN). Nonetheless, GANs derived from single datasets generally allow us to envisage only one side of the regulatory network, even under the particular condition of study. The circumstance is well demonstrated by inconsistent GANs of individual datasets proposed for similar experimental conditions, which always leads to ambiguous interpretation. Here, pan- and core-gene association networks (pan- and core-GANs), analogous to the pan- and core-genome concepts, are proposed to increase the power of inference through the integration of multiple, diverse datasets. The core-GAN represents the consensus associations of genes that were inferred from all individual networks. On the other hand, the pan-GAN represents the extensive gene-gene associations that occurred in each individual network. The pan- and core-GANs prospects were demonstrated based on three time series microarray datasets in leaves of Arabidopsis thaliana grown under diurnal conditions. We showed the overall performance of pan- and core-GANs was more robust to the number of data points in gene expression data compared to the GANs inferred from individual datasets. In addition, the incorporation of multiple data broadened our understanding of the biological regulatory system. While the pan-GAN enabled us to observe the landscape of gene association system, core-GAN highlighted the basic gene-associations in essence of the regulation regulating starch metabolism in leaves of Arabidopsis.

Cloning, Expression, Purification, and Insecticidal Activity of a Novel Cry1Na3 Toxin From Bacillus thuringiensis BRC-ZYR2

Bacillus thuringiensis produces a variety of insecticidal crystal proteins (ICPs). Genome sequencing is a promising strategy for detecting and identifying B. thuringiensis ICPs, which are of great interest to the biocontrol field. In this study, a novel ICP gene was cloned from B. thuringiensis BRC-ZYR2 based on genomic data from 454 GS-FLX Titanium sequencing and an analysis of the results using the B. thuringiensis Toxin_Scanner ( http://bcam.hzaubmb.net/BtToxin_scanner/index.php ). cry1Na3 designated by the B. thuringiensis Toxin Nomenclature Committee, encoded a 601-amino acid, 68.0-kDa protein that exhibited 95% identity with Cry1Na1 and 99% identity with Cry1Na2. Cry1Na3 contained three conserved domains commonly found in three-domain ICPs. Cry1Na3 was toxic to Plutella xylostella (L.) and Ostrinia furnacalis (Guenée), with LC 50 values of 3.69 μg/ml and 31.30 μg/ml, respectively. However, Laodelphax striatellus (Fallén) nymphs were unaffected when fed purified Cry1Na3 (250 μg/ml) in their diet. Spodoptera exigua (Hübner) and Colaphellus bowringi (Baly) larvae survived even when the concentration of Cry1Na3 protein reached 500 μg/ml. Cry1Na3 is a promising agent for the control of lepidopteran insect pests.

A brief review of software tools for pangenomics

Since the proposal for pangenomic study, there have been a dozen software tools actively in use for pangenomic analysis. By the end of 2014, Panseq and the pan-genomes analysis pipeline (PGAP) ranked as the top two most popular packages according to cumulative citations of peer-reviewed scientific publications. The functions of the software packages and tools, albeit variable among them, include categorizing orthologous genes, calculating pangenomic profiles, integrating gene annotations, and constructing phylogenies. As epigenomic elements are being gradually revealed in prokaryotes, it is expected that pangenomic databases and toolkits have to be extended to handle information of detailed functional annotations for genes and non-protein-coding sequences including non-coding RNAs, insertion elements, and conserved structural elements. To develop better bioinformatic tools, user feedback and integration of novel features are both of essence.

Analysis of the Pantoea ananatis pan-genome reveals factors underlying its ability to colonize and interact with plant, insect and vertebrate hosts

Background

Pantoea ananatis is found in a wide range of natural environments, including water, soil, as part of the epi- and endophytic flora of various plant hosts, and in the insect gut. Some strains have proven effective as biological control agents and plant-growth promoters, while other strains have been implicated in diseases of a broad range of plant hosts and humans. By analysing the pan-genome of eight sequenced P. ananatis strains isolated from different sources we identified factors potentially underlying its ability to colonize and interact with hosts in both the plant and animal Kingdoms.

Results

The pan-genome of the eight compared P. ananatis strains consisted of a core genome comprised of 3,876 protein coding sequences (CDSs) and a sizeable accessory genome consisting of 1,690 CDSs. We estimate that ~106 unique CDSs would be added to the pan-genome with each additional P. ananatis genome sequenced in the future. The accessory fraction is derived mainly from integrated prophages and codes mostly for proteins of unknown function. Comparison of the translated CDSs on the P. ananatis pan-genome with the proteins encoded on all sequenced bacterial genomes currently available revealed that P. ananatis carries a number of CDSs with orthologs restricted to bacteria associated with distinct hosts, namely plant-, animal- and insect-associated bacteria. These CDSs encode proteins with putative roles in transport and metabolism of carbohydrate and amino acid substrates, adherence to host tissues, protection against plant and animal defense mechanisms and the biosynthesis of potential pathogenicity determinants including insecticidal peptides, phytotoxins and type VI secretion system effectors.

Conclusions

P. ananatis has an ‘open’ pan-genome typical of bacterial species that colonize several different environments. The pan-genome incorporates a large number of genes encoding proteins that may enable P. ananatis to colonize, persist in and potentially cause disease symptoms in a wide range of plant and animal hosts.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-404) contains supplementary material, which is available to authorized users.

Bacillus thuringiensis Is an Environmental Pathogen and Host-Specificity Has Developed as an Adaptation to Human-Generated Ecological Niches

Bacillus thuringiensis (Bt) has been used successfully as a biopesticide for more than 60 years. More recently, genes encoding their toxins have been used to transform plants and other organisms. Despite the large amount of research on this bacterium, its true ecology is still a matter of debate, with two major viewpoints dominating: while some understand Bt as an insect pathogen, others see it as a saprophytic bacteria from soil. In this context, Bt's pathogenicity to other taxa and the possibility that insects may not be the primary targets of Bt are also ideas that further complicate this scenario. The existence of conflicting research results, the difficulty in developing broader ecological and genetics studies, and the great genetic plasticity of this species has cluttered a definitive concept. In this review, we gathered information on the aspects of Bt ecology that are often ignored, in the attempt to clarify the lifestyle, mechanisms of transmission and target host range of this bacterial species. As a result, we propose an integrated view to account for Bt ecology. Although Bt is indeed a pathogenic bacterium that possesses a broad arsenal for virulence and defense mechanisms, as well as a wide range of target hosts, this seems to be an adaptation to specific ecological changes acting on a versatile and cosmopolitan environmental bacterium. Bt pathogenicity and host-specificity was favored evolutionarily by increased populations of certain insect species (or other host animals), whose availability for colonization were mostly caused by anthropogenic activities. These have generated the conditions for ecological imbalances that favored dominance of specific populations of insects, arachnids, nematodes, etc., in certain areas, with narrower genetic backgrounds. These conditions provided the selective pressure for development of new hosts for pathogenic interactions, and so, host specificity of certain strains.

Seasonally variable intestinal metagenomes of the red palm weevil (Rhynchophorus ferrugineus)

The intestinal microbes residing in the red palm weevil (RPW, Rhynchophorus ferrugineus) larva consume tender interior fibrous tissues of date palm trunks. The understanding of such microbiota at molecular level provides vital clues for the biological control of this devastating pest. Using pyrosequencing and shotgun strategy, we first study taxonomic profiles of the microbiota sampled at different months (March, July and November), and then confirm the impact of high-temperature stress on the microbial populations based on data from 16S rRNA amplicons using both field and laboratory samples. We further identify Klebsiella pneumoniae in November and Lactococcus lactis in July as the dominant species of the microbiota. We find that the RPW gut microbiota degrades polysaccharides and sucrose with hydrolases and that different active bacterial species in November and July are responsible for the symbiotic relationship between the microbiota and the host. Our results provide vital information for pest control and cellulolytic bacterial species characterization.