Molecular typing of Bacillus thuringiensis serovars by RAPD-PCR

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.

Identification of Bacillus Probiotics Isolated from Soil Rhizosphere Using 16S rRNA, recA, rpoB Gene Sequencing and RAPD-PCR

Some Bacillus species, especially Bacillus subtilis and Bacillus pumilus groups, have highly similar 16S rRNA gene sequences, which are hard to identify based on 16S rDNA sequence analysis. To conquer this drawback, rpoB, recA sequence analysis along with randomly amplified polymorphic (RAPD) fingerprinting was examined as an alternative method for differentiating Bacillus species. The 16S rRNA, rpoB and recA genes were amplified via a polymerase chain reaction using their specific primers. The resulted PCR amplicons were sequenced, and phylogenetic analysis was employed by MEGA 6 software. Identification based on 16S rRNA gene sequencing was underpinned by rpoB and recA gene sequencing as well as RAPD-PCR technique. Subsequently, concatenation and phylogenetic analysis showed that extent of diversity and similarity were better obtained by rpoB and recA primers, which are also reinforced by RAPD-PCR methods. However, in one case, these approaches failed to identify one isolate, which in combination with the phenotypical method offsets this issue. Overall, RAPD fingerprinting, rpoB and recA along with concatenated genes sequence analysis discriminated closely related Bacillus species, which highlights the significance of the multigenic method in more precisely distinguishing Bacillus strains. This research emphasizes the benefit of RAPD fingerprinting, rpoB and recA sequence analysis superior to 16S rRNA gene sequence analysis for suitable and effective identification of Bacillus species as recommended for probiotic products.

Molecular typing of native Bacillus thuringiensis isolates from diverse habitats in India using REP-PCR and ERIC-PCR analysis

Bacillus thuringiensis is a bacterium of great agronomic and scientific interest. The subspecies of this bacterium colonize and kill a large variety of host insects and even nematodes, but each strain does so with a high degree of specificity. Therefore molecular typing and diversity analysis of B. thuringiensis has enormous importance for discrimination of strains isolated from different sources. In this study, 113 native B. thuringiensis isolates collected from diverse habitats and locations in India and 27 B. thuringiensis type strains obtained from the Bacillus Genetic Stock Centre (BGSC), Ohio State University, USA and used as reference, were analyzed for molecular typing. Genotypic data of 140 B. thuringiensis isolates and type strains was generated by using REP-PCR and ERIC-PCR primers and unweighted pair group method with arithmetic mean (UPGMA) analysis using NTSYSpc2.2 and grouped into 4 main clusters. All the groups have isolates from diverse origins. No group was found to represent any specific origin or location. The observed patterns of REP-PCR and ERIC-PCR pattern were discriminatory enough to reveal differences in the B. thuringiensis isolates and reference strains. The resolution power and marker index of the ERIC-PCR (RP 9.39, MI 6.34) was found to be higher than that of the REP-PCR (RP 6.20, MI 4.48). The REP-PCR and ERIC-PCR markers have been found to be useful for discrimination of B. thuringiensis isolates and reference strains. ERIC-PCR was the more informative of the two techniques. This study showed that the B. thuringiensis isolates collected from diverse habitats in India had a high degree of genetic diversity.

Evidence of Bacillus thuringiensis intra-serovar diversity revealed by Bacillus cereus group-specific repetitive extragenic palindromic sequence-based PCR genomic fingerprinting

Bacillus thuringiensis is classified into serovars on the basis of H-flagellar antigens. Several alternative typing methods have been described. Among them, a B. cereus group-specific repetitive extragenic palindromic (Rep)-PCR fingerprinting technique was shown to be discriminative and able to identify B. thuringiensis serovars. The aim of this study was to investigate the genomic diversity and relationship among B. thuringiensis strains collected from different Argentinean ecosystems. Thirty-seven B. thuringiensis reference strains and 131 Argentinean isolates were analyzed using a B. cereus group-specific Rep-PCR. Fourteen different patterns were identified among the Argentinean isolates. Eight could not be associated to any pattern obtained from a reference strain. The pattern identical to the serovar kurstaki HD-1 strain was the most frequently identified in 68 native isolates. The profiles allowed tracing a single dendrogram with two groups and eight main lineages. Some strains showed distinctive patterns despite belonging to the same serovar. An intraspecific diversity resulted from this analysis that was highlighted by this technique since strains from a given serovar showed distinct profiles. This study may help to establish a system of B. thuringiensis classification with a higher discrimination level than established by the H antigen serotyping.

Characterization of native Bacillus thuringiensis strains by PCR-RAPD based fingerprinting

Seventy isolates of Bacillus thuringiensis were isolated from soil samples collected from cotton fields. These isolates were characterized by randomly amplified poylmorphic DNA (RAPD) markers to determine their genetic diversity pattern based on their source of origin. Different random decamer primers were used for RAPD amplification, which generated a total of 1935 fragments; of these 1865 were polymorphic and 68 monomorphic. The primers OPA03, OPA08, OPD14, OPD19, OPD20, OPE17 and OPD19 produced 100% polymorphic fragments, whereas primers OPC06, OPC20 and OPD17 produced 20, 31 and 17 monomorphic fragments, respectively. When the RAPD banding pattern data was subjected to dendrogram construction, the 70 isolates fell into two separate clusters, cluster I and cluster II, which includes 26 and 44 B. thuringiensis isolates, respectively. These two main clusters were further divided into four subclusters at Eucledian distance of 150 and 80% similarity index. All primers showed amplification and indicated the good diversity of B. thuringiensis isolates. The RAPD pattern showed 4-10 bands per isolate, with MWt in the range of 0.4-3.5 Kb and an average of 193.5 fragments were produced per primer. The primer OPE17 was found to be the most discriminatory as it produced 286 polymorphic bands.

Phenotypic and genotypic features of new autoagglutinating Bacillus thuringiensis strains

A total of 28 autoagglutinating strains of Bacillus thuringiensis were isolated from different ecologic niches and distinct sites. Twenty-six strains demonstrated toxicity to mosquito larvae of Aedes aegypti and Culex quinquefasciatus. The electrophoretic protein profiles of the crystal components were studied. Twenty-three out of the 28 strains showed the same larvicidal activity and the same protein profiles as B. thuringiensis serovar israelensis. Using isoenzyme analysis (MLEE), it was observed the presence of three electrophoretic types (ETs). The mosquitocidal strains grouped into one ET. The random amplified polymorphic DNA analysis (RAPD) was evaluated using six primers, which demonstrated three different patterns for the 28 autoagglutinating strains, allowing correlation of the profiles obtained with the toxicity observed in the bioassays. The RAPD patterns for mosquitocidal strains were identical to the one of serovar israelensis. However, to strains of low toxicity, each primer generated distinctive RAPD patterns, which demonstrated that these strains belong to different serovars. Although the antigenic classification the 26 autoagglutinating strains of B. thuringiensis could not be determined by classical flagellar serotyping, MLEE and RAPD profiles proved these strains to be compatible with B. thuringiensis serovar israelensis.

The Bacillus cereus group: novel aspects of population structure and genome dynamics

AIMS

To provide new insights into the population and genomic structure of the Bacillus cereus group of bacteria.

METHODS AND RESULTS

The genetic relatedness among B. cereus group strains was assessed by multilocus sequence typing (MLST) using an optimized scheme based on seven chromosomal housekeeping genes. A set of 48 strains from different clinical sources was included, and six clonal complexes containing several genetically similar isolates from unrelated patients were identified. Interestingly, several clonal groups contained strains that were isolated from similar human sources. Furthermore, comparative whole genome sequence analysis of 16 strains led to the discovery of novel ubiquitous genome features of the B. cereus group, such as atypical group II introns, IStrons, and hitherto uncharacterized repeated elements.

CONCLUSIONS

The B. cereus group constitutes a coherent population unified by the presence of ubiquitous and specific genetic elements which do not show any pattern, either in their sequences or genomic locations, which allows to differentiate between the member species of the group. Nevertheless, the population is very dynamic, as particular lineages of clinical origin can evolve to form clonal complexes. At the genome level, the dynamic behaviour is indicated by the presence of numerous mobile and repeated elements.

SIGNIFICANCE AND IMPACT OF THE STUDY

The B. cereus group of bacteria comprises species that are of medical and economic importance. The MLST data, along with the primers and protocols used, will be available in a public, web-accessible database (http://mlstoslo.uio.no).

Fingerprinting of Bacillus thuringiensis type strains and isolates by using Bacillus cereus group-specific repetitive extragenic palindromic sequence-based PCR analysis

A total of 119 Bacillus thuringiensis strains (83 type strains and 26 native isolates), as well as five B. cereus group species, were analyzed by repetitive extragenic palindromic sequence-based PCR analysis (Rep-PCR) fingerprinting. Primers Bc-REP-1 and Bc-REP-2 were specifically designed according to an extragenic 26-bp repeated sequence found in the six B. cereus group genomes reported. A total of 47 polymorphic bands were detected, and the patterns varied from 5 to 13 bands in number and from 0.2 to 3.8 kb in size. Virtually each type strain showed a distinctive B. cereus (Bc)-Rep-PCR pattern, except for B. thuringiensis serovars dakota (H serotype 15 [H15]) and sotto (H4a,4b), as well as serovars amagiensis (H29) and seoulensis (H35), which shared the same patterns. As expected, serovar entomocidus (H6) and its biovar subtoxicus showed an identical pattern; similarly, serovars sumiyoshiensis (H3a,3d) and fukuokaensis (H3a,3d,3e), which share two antigenic determinants, also showed identical Bc-Rep-PCR patterns. Interestingly, serovars israelensis (H14) and malaysiensis (H36), which share several phenotypic attributes, also showed identical Bc-Rep-PCR patterns. Native, coleopteran-active strains, including the self-agglutinated LBIT-74 strain, showed Bc-Rep-PCR patterns identical or very similar to that of the tenebrionis strain. Likewise, native mosquitocidal strains (including some self-agglutinated strains) also showed patterns identical or very similar to that of the serovar israelensis IPS-82 strain. Additionally, native beta-exotoxin-producing strains from serovar thuringiensis showed patterns identical to that of the B. thuringiensis type strain. The B. cereus group-specific Bc-Rep-PCR fingerprinting technique was shown to be highly discriminative, fast, easy, and able to identify B. thuringiensis serotypes, including nonflagellar and self-agglutinated strains.

Population structure and evolution of the Bacillus cereus group

Representative strains of the Bacillus cereus group of bacteria, including Bacillus anthracis (11 isolates), B. cereus (38 isolates), Bacillus mycoides (1 isolate), Bacillus thuringiensis (53 isolates from 17 serovars), and Bacillus weihenstephanensis (2 isolates) were assigned to 59 sequence types (STs) derived from the nucleotide sequences of seven alleles, glpF, gmk, ilvD, pta, pur, pycA, and tpi. Comparisons of the maximum likelihood (ML) tree of the concatenated sequences with individual gene trees showed more congruence than expected by chance, indicating a generally clonal structure to the population. The STs followed two major lines of descent. Clade 1 comprised B. anthracis strains, numerous B. cereus strains, and rare B. thuringiensis strains, while clade 2 included the majority of the B. thuringiensis strains together with some B. cereus strains. Other species were allocated to a third, heterogeneous clade. The ML trees and split decomposition analysis were used to assign STs to eight lineages within clades 1 and 2. These lineages were defined by bootstrap analysis and by a preponderance of fixed differences over shared polymorphisms among the STs. Lineages were named with reference to existing designations: Anthracis, Cereus I, Cereus II, Cereus III, Kurstaki, Sotto, Thuringiensis, and Tolworthi. Strains from some B. thuringiensis serovars were wholly or largely assigned to a single ST, for example, serovar aizawai isolates were assigned to ST-15, serovar kenyae isolates were assigned to ST-13, and serovar tolworthi isolates were assigned to ST-23, while other serovars, such as serovar canadensis, were genetically heterogeneous. We suggest a revision of the nomenclature in which the lineage and clone are recognized through name and ST designations in accordance with the clonal structure of the population.