Variable Region of 16s rRNA is Essential For the Identification of Group 1 Mosquito-Pathogenic Strains of Lysinibacillus
Ajeet Kumar Mohanty1, Sandeep Garg2, Kulvir Dhindsa3 and Ashwani Kumar1*
1National Institute of Malaria Research (ICMR), India
2Department of Microbiology, Goa University, India
3Department of Zoology, D. A. VCollege, India
Submission: November 27, 2016; Published: January 31, 2017
*Corresponding author: Ashwani Kumar, M.Sc., Ph D, Scientist 'F & Officer-in-Charge, National Institute of Malaria Research (ICMR), Field Unit, DHS Building, Campal, Panaji, Goa- 403001, India, Tel: +91 832 2222444; Email: ashwani7@gmail.com
How to cite this article: Ajeet K M, Sandeep G, Kulvir D, Ashwani K. Variable Region of 16s rRna is Essential For the Identification of Group 1 Mosquito-Pathogenic Strains of Lysini bacillus. Adv Biotech & Micro. 2017; 2(2): 555583. DOI: 10.19080/AIBM.2017.02.555583
Abstract
A study was carried out to characterize and identify the potential mosquito-pathogenic bacilli strains isolated from the soil samples obtained from various locations in Goa, India. Six isolates were characterized morphologically, biochemically and phylogentically. Phylogenetic relationships of these strains were determined based on comparison of the 16S rRNA sequences with that of the closest Lysinibacillus species obtained from nucleotide database. Their mosquitocidal abilities were confirmed by bioassay and presence of genes responsible for pathogenicity and toxin production. Molecular characterization by partial sequencing of 16S rRNA gene has confirmed their identity as Lysinibacillus sp. Group 1.
Comparison of the 16S rRNA sequences with that of the five groups of "Bacillus sphaericus like organisms”, re-designated off lately as Lysinibacillus, revealed that a variable region of 695 bp and hypervariable region of 22 bp were highly conserved among these five groups. These isolates were found effective against Culex quinquefasciatus larvae with LC50 values ranging from 0.018 to 0.58 ppm. The PCR amplification also showed the presence of binA, binB binary genes and mtx1, mtx2, mtx3 mosquitocidal genes. Larvicidal activity against Culex quinquefasciatus appeared related with expression of binary toxins. These strains were found more potent than the commercial strain and hence could be formulated for controlling Culex species transmitters of Japanese encephalitis and filariasis.
Keywords: Biocontrol agent; Binary toxins; Hypervariable; Lysinibacillus strains; 16S rRNA gene
Abbreviations: Btx: Binary Toxins; Mtx: Mosquitocidal Toxins; HV: Hyper Variable; NJ: Neighbor Joining; NYSM: Nutrient Yeast Sporulating Medium
Introduction
Some strains of "Bacillus sphaericus like organisms” are toxic to dipteran larvae and used operationally for control of mosquito-borne diseases. These micro-organisms are aerobic, spore-forming bacteria producing spherical or oval endospores and contain endotoxins responsible for the killing of mosquito larvae [1] The larvicidal strains contain two types of insecticidal toxins viz. binary and mosquitocidal. The binary toxins (Btx) are spore associated toxins produced during sporulation phase and form crystals in the mother cell. The crystal comprises of two toxic proteins with molecular weights of about 42 kDa (binA) and 51 kDa (binB). Both proteins are required in equal amounts for larval toxicity [2,3].
In contrast, mosquitocidal toxins (Mtx) are produced during vegetative stage. There are three Mtx toxins (Mtx1, Mtx2 and Mtx3) that have been identified in some mosquitocidal strains [4-6]. Interestingly, Mtx proteins are degraded by extracellular proteases which are produced when the cell growth reaches sporulation phase [7 ]. Highly toxic strains are known to produce both Btx and Mtx toxins, while other strains synthesize either of them [8]. Krych et al. [9] on the basis of DNA-DNA homology demonstrated that Bacillus sphaericus is a conglomerate of 5 distinct homology groups. Group I represents Bacillus sphaericus sensu stricto, while Group II consists of two subgroups, IIA & IIB. All mosquito pathogenic strains are placed in subgroup IIA. Group IIB consists of non-pathogenic Bacillus fusiformis.
Alexander & Priest [10] attempted to use phenotypic features to separate these groups. The toxic strains of Group IIA produce ribosomal RNA gene restriction patterns which are different from Group IIB and non-toxic Bacillussphaericus sensu stricto [11]. De muro & Priest [12] eveloped an oligonucleotide probe based on specific region of 16S rDNA which differentiates Group IIA from other DNA homology groups. However, Jhanz et al. [13] showed that this probe also hybridizes with non-toxic strains. RAPD and repetitive primer amplified polymorphic DNA analysis showed close relatedness among mosquito pathogenic strains belonging to homology group IIA [14,15]. Nakamura [16] redistributed Bacillus sphaericus-like organism based on 16S rRNA gene sequence into seven distinct clusters. In this schema, mosquitocidal and few non-mosquitocidal strains were clustered in Group 1 which corresponded to Krych Group IIA and closely linked to Group 2 representing Krych Group IIB named as Bacillus fusiformis.
Group 1 to 4 of Bacillus sphaericus like organisms were re-designated as new genus Lysinibacillus on the basis of the presence of lysine in the cell wall [17]. In the past, different techniques have been used for grouping and classifying mosquito pathogens. Except the pathogenicity to mosquito larvae, there are no definitive characters that can separate this from other taxa, therefore the taxonomic status of this group is still unclear.
Here, we have studied variable and hypervariable regions of 16S rRNA gene specific to the strains of Lysinibacillus which correspond to group 1 to 5 of Nakamura [16]. This study was aimed at finding whether the distinctive phenotypic characters, the presence of unique nucleotide sequences and mosquitocidal genes altogether could form basis for distinguishing this group from other spore forming Lysinibacillus.
Materials and Methods
Source of strain
The bacterial strains were isolated from soil samples which were collected from various mosquito breeding habitats viz., stagnant ponds, paddy fields and mangrove vegetations at different location in Goa, India [18]. Stock cultures were maintained on NYSM agar slants stored at 4 °C.
Characterization of the isolates
Morphological, biochemical and physiological characteristics of the isolates were studied as per Krych et al. [9] and Alexander & Priest [10].
Amplification of 16S rRNA gene of the isolates
The genomic DNA was extracted from 24hrs old culture grown in J-Broth using the method described in Sambrook et al. [19]. Gene of 16S rRNA was amplified using the universal primers, a forward primer (S-D-Bact-0011-a-5-17: 5'-GTTTGATCCTGGCTCAG-3') and a reverse primer (S-*-Univ- 1392-b-A-15: 5'-ACGGGCGGTGTGTNC-3'). The conditions for PCR amplification were used as described by Prabhu et al. [20]. PCR products were purified using purification kit (Quiagen). Purified PCR product obtained was outsourced to Bangalore Genei for sequencing. The sequences obtained were submitted in the Gene Bank under the accession numbers FJ473362 (KSD- 1), FJ473363 (KSD-2), FJ473364 (KSD-3), FJ473365 (KSD-4), FJ473368 (KSD-8) and FJ473369 (KSD-7).
The sequences obtained were compared with those available in NCBI database using BLASTn. The sequences from phylogenetic groups described by Nakamura [16] for Bacillus sphaericus like organisms were compared. Clustal X version 2.0.7 was used to generate multiple sequence alignment between these selected sequences. Consensus sequences obtained from these multiple alignment within groups were compared with sequences available in NCBI databases. Consensus sequences for these groups were then analyzed in relation to variable and hypervariable (HV) regions to develop phylogenetic tree using neighbor joining (NJ) method. Tree was obtained with 1000 seeds and 10000 bootstraps. Final tree was rooted and drawn using MEGA4 [21].
PCR amplification of genes encoding the mosquitocidal toxins
Genomic DNA of isolates was subjected to PCR amplification to detect the presence of genes encoding toxins. Primers used for amplification of toxin genes binA, binB, mtx1, mtx2 and mtx3 are listed in Table 1. Conditions used for amplification were as described in [6,22,23].
Source of immature of test mosquito species
Culex quinquefasciatus larvae were obtained from the insectary of National Institute of Malaria Research, Field Unit, Goa, where this mosquito species is maintained at a temperature of 27±2 °C, relative humidity of 70±5% and a photoperiod: scotoperiod of 12:12 h (light: dark). From this laboratory bred mosquitoes, 3rd instars larvae were selected and used for larval bioassays to determine larvicidal activity of bacterial strains.
Preparation of Lyophilized powder for Bioassay
Isolates were grown in a 250ml conical flask containing 100ml of NYSM broth and incubated at 30 °C on a rotary shaker for 96 hrs. Spore crystals were harvested by centrifugation at 5000rpm for 10 min, cell pellet was washed twice with sterilized distilled water and the final pellet was kept for lyophilization which was carried out in Delvac Model (Lyodel) at temperature -40 °C. Vacuum was applied for 4-6 hours for complete drying. Lyophilized culture was stored at 4 °C until use.
Larval Bioassays
Stock dilutions were prepared by dissolving lyophilized powder in sterile distilled water. Dilutions were made to obtain appropriate range of dosage [24]. For each of the doses, 3 replicates of 25 healthy III instar larvae of Culex quinquefasciatus were introduced in 250ml plastic bowl containing 100ml of sterile distilled water. Concurrent control was maintained under similar conditions without the addition of spore suspension. The mortality was recorded by counting live larvae after 24 & 48 hours. LC50 were determined by Probit analysis with SPSS PASW 18.0 indicating mean and standard errors [25].
Results and Discussion
Phenotypic characterisation of isolates
Colonies of all six mosquito pathogenic isolates were white- creamish in colour. The vegetative cells were Gram's positive, long and rod-shaped. Sporulating cells consisted of racket shaped swollen sporangia were containing round endospores associated with the crystal proteins. These isolates did not produce nitrate reductase, amylases, tyrosinase and showed the production of catalase, urease and acid from glucose. All isolates grew in broth with 2 to 7% NaCl, sensitive to erythromycin and resistant to streptomycin. These isolates cleared adenine and were unable to utilize citrate and L-isoleucine in the medium.
Earlier, mosquito-pathogenic bacterial strains were classified on the basis of limited biochemical tests as Bacillus sphaericus Krych et al. [9] demonstrated using DNA-DNA homology that Bacillus sphaericus actually a conglomerate of 5 homology groups. The mosquito-pathogenic strains belong to subgroup IIA. Non-pathogenic strains belong to subgroup IIB were identified as Bacillus fusiformis. Alexander & Priest [10] used several phenotypic tests in order to give features specific to mosquito-pathogenic strains belong to subgroup IIA. The above isolates showed all the features in agreement to the described feature of subgroup IIA. It was noted that subgroup IIA pathogens separated from the subgroup IIB non-pathogens only by the ability of the former to clear adenine from nutrient agar and sensitivity to erythromycin. An additional feature is mosquito-pathogenicity, which is based on the production of toxins. Recently, members of group I to IV of Krych et al. [9] were re-designated as members of genus Lysinibacillus due to presence of lysine in the cell wall [17].
Phylogenetic characterization
Sequences of 16S rRNA genes obtained from the isolates were used to search NCBI databases using BLASTn which showed high similarity with Lysinibacillus fusiformis and Lysinibacillus sphaericus. (Figure 1) presents phylogenetic tree based on comparison of the variable region of 16S rRNA sequences generated with that of sequences of species belonging to Lysinibacillus genus, including strains of Bacillus sphaericus like organisms earlier reported as pathogenic to mosquitoes in the GenBank database. Phylogenetic tree separates these organisms into five groups. All the six mosquito-pathogenic isolates were clustered in Group 1 in accordance with Nakamura [16] that correspond to sub-group IIA as described above.
Consensus sequences derived in this study for five groups of Lysinibacillus are shown in (Figure 2). A variable region (V) of 695 bp nucleotide position stretching from 62 to 756 bps and a hypervariable (HV) region of 22 bps at nucleotide position from 178 to 199 bps were found among all these five groups. The HV regions were conserved and specific to respective groups. Interestingly, the highly conserved HV region in Group 1 and Group 2 were different only by three nucleotides at position 197199 bps. Additionally, we found six single nucleotide differences in the variable regions of Group 1 and 2 which were used to differentiate members of these groups. These unique nucleotide differences among these two groups are at nucleotide position 85, 221, 571, 586, 665 and 720.
Group 2 and Group 3 of 16S rRNA gene alignments correspond to sub-group IIB and Group I of DNA-DNA homology and were named as Lysinibacillus fusiformis and Lysinibacillus sphaericus, respectively. Group 4 has been named as Lysinibacillus boronitolerans which corresponds to group III of DNA-DNA homology [26]. So far Group 1 containing mosquito pathogenic strains and Group 5 (DNA-DNA homology Group IV) has not been recognised as a separate species.
The phylogenetic tree shown in (Figure 1) included the consensus sequences drawn for five groups. Sequence similarities of the above mentioned variable region between these five groups are shown in (Table 2). This analysis reveals that Group 1 and Group 2 have maximum commonalities of only 98%. Our findings showed that variable region was sufficient in distinguishing the five groups as the nodes were showing very high bootstrap values (Figure 1). These findings match with the phylogenetic tree drawn using 16S rDNA sequences more than 1400 bases [16]. Although, HV regions showed the signature sequences of Group 1 and Group 2 but were incapable of separating these groups on its own during phylogenetic analysis. This region was earlier used by de Muro et al. [12] for designing probe for identifying the mosquito-pathogenic isolates.
Mosquito-pathogenicity of isolates
+, positive; -, negative; V, variable;
*Data from other study as follows: Group2, Group3, Group4 and Group5 organism's biochemical characteristics as described in Alexander& Priest [10].
All 6 isolates were tested against laboratory reared III instar larvae of Culex quinquefasciatus. These isolates were found highly toxic to Culex quinquefasciatus larvae, with LC50 values ranging from 0.05to 0.87 ppm at 24 hrs post treatment and 0.018 to 0.58 ppm at 48 hrs (Table 3). Isolate KSD-2 had lowest LC50 value of 0.05 & 0.018 ppm at 24 and 48 hrs respectively and hence was most lethal. Lethal conc. of KSD-2 was 32 folds lower as compared with commercial Bacillus sphaericus 101 (H5a5b)
PCR amplification of binary and mosquitocidal genes confirmed the presence of binA, binB, mtx1, mtx2 and mtx3 in all 6 isolates as shown in (Figure 3). Larvicidal activity was found to be primarily associated with binary toxins. This could be explained as only sporulating stage was tested for determination of larvicidal activity. Previous studies also showed that highly toxic strains are known to produce both binary and Mtx toxin proteins, while other strains synthesize either of these toxins [8]. The DNA sequences of the genes encoding these proteins are highly homologous between the strains, although they were isolated from different continents [27]. Previously, pathogenicity of Bacilli was considered important for distinguishing taxa as there were no reliable phenotypic characters [10], but now probe based on 16S rRNA and toxin genes are available to segregate mosquito-pathogenic strains morphologically similar to Bacillus sphaericus like organisms [12,28]. In the present study, the presence of both bin and mtx genes in the mosquito- pathogenic strains has been confirmed.
Proposal for species novo
Classical characterization methods for describing Bacillus sphaericus have failed to provide enough taxonomic information [29]. Krych et al. [9] were the first to identify 5 distinct DNA-DNA homology groups among the 50 strains, which include 7 strains showing pathogenicity to mosquito larvae, which were placed in the sub-Group IIA. There was about 62% DNA-DNA sequence homology, with ATm of about 7oC between the strains of subgroup IIA and IIB. It is generally argued that strains within a species show greater than 70% DNA sequence homology and thermal hybrid stability (ATm) of less than 6oC [30]. This amply suggested that sub-group IIA and IIB cannot be allocated to the same species. The creation of a separate species for the mosquito pathogens was contemplated earlier but due to lack of positive phenotypic tests, other than mosquito pathogenicity, its identification as a taxon has been delayed. Alexander &Priest [10] used several phenotypic tests and tried to develop numerical taxonomic group in order to identify features specific to mosquito-pathogenic strains belonging to subgroup IIA. Further support for distinction of sub-group IIA pathogens from non-pathogenic sub-group IIB became possible by rRNA gene pattern [11], RAPD finger print and isoenzyme analysis [31].
These studies revealed a very low similarity between the two subgroups and led to the conclusion that the pathogenic strains are an independent lineage distinct from the other round spore forming bacilli of sub-group IIB. Nakamura [16] attempted to raise the status of sub-group IIA and sub-group IIB to independent groups namely Group 1 and Group 2, respectively using 16S rRNA sequences. Now with re-designation of the genus due to presence of Lysine in cell wall the sub-group IIB is known as Lysinibacillus fusiformis.
PROBIT model: PROBIT(p) = Intercept + BX (Covariates X (dose) are transformed using the base 10.000 logarithm)
Acknowledgement
We thank the Indian Council of Medical Research, National Institute of Malaria Research, New Delhi for institutional support. We are thankful to the field and laboratory staff of National Institute of Malaria Research (ICMR), Field Unit, Goa, India. The authors also thank Dr. Barbuddhe of ICAR, Research Complex, Goa for providing the Lyophiliser facility. This paper bears the NIMR publication screening committee approval no. 049/2016.
The present investigation clearly identified the differences of 16S rRNA sequences between various members of Lysininbacillus including Lysinibacillus fusiformis and group IIA containing mosquito pathogenic isolates. The phylogenetic tree constructed by taking only variable region was also stable at its nodes with high value of bootstrap. The distinct features specific to each of the species of Lysinibacillus are listed in (Table 4). These taxonomic signatures could be reliably utilized for the screening and identification of mosquito pathogenic strains. Due to the taxonomic value of the unique sequences observed based on variable region in the mosquito-pathogenic strains and phenotypic characters, we propose to designate mosquito- pathogenic strains as a new species Lysinibacillus kychi sp. novo that are hitherto clustered in Group IIA by Krych et al. [9] and Group 1 by Nakamura [16].
References
- Yousten AA (1984) Bacillus sphaericus: microbiological factors related to its potential as a mosquito larvicide. Adv Biotechnol Processes 3: 315-343.
- Baumann P, Clark MA, Baumann L, Broadwell AH (1991) Bacillus sphaericus as a Mosquito Pathogen: Properties of the Organism and Its Toxins. Microbiol Reviews 55(3): 425-436.
- Charles JF, LeRoux CN, Delecluse A (1996) Bacillus sphaericus Toxins: Molecular Biology and Mode of Action. Annual Review Entomol 41: 451-472.
- Thanabalu T, Hindley J, Jackson Yap J, Berry C (1991) Cloning, sequencing, and expression of a gene encoding a 100-kilodalton mosquitocidal toxin from Bacillus sphaericus SSII-1. J Bacteriol 173(9): 2776-278S.
- Thanabalu T, Porter AG (1996) A Bacillus sphaericus gene encoding a novel type of mosquitocidal toxin of 31.8 kDa. Gene 170(1): 8S-89.
- Liu JW, Porter AG, Wee BY, Thanabalu T (1996) New gene from nine Bacillus sphaericus strains encoding highly conserved 3S.8-kilodalton mosquitocidal toxin. Appl Environ Microbiol 62(6): 2174- 2176.
- Yang YL, Wang A, Gaviria Z, Berry C (2007) Proteolytic stability of insecticidal toxins expressed in recombinant BaciIIi. Appl Environ Microbiol 73(1): 218-22S.
- Liu JW, Hindley WJ, Porter AG, Priest FG (1993) New high-toxicity mosquitocidal strains of Bacillus sphaericus lacking a 100-kilodalton toxin gene. Appl Environ Microbiol S9(10): 3470-3473.
- Krych VK, Johnson JL, Yousten AA (1980) Deoxyribonucleic Acid Homologies among Strains of BaciIIus sphaericus. Int J Syst Bacteriol 30(2): 476-484.
- Alexander B, Priest FG (1990) Numerical classification and identification of Bacillus sphaericus including some strains pathogenic for mosquito larvae. J Gen Microbiol 136(2): 367-376.
- Aquino de Muro M, Mitchell WJ, Priest FG (1992) Differentiation of mosquito-pathogenic strains of Bacillus sphaericus from non-toxic varieties by ribosomal RNA gene restriction patterns. J Gen Microbiol 138(6): 11S9-1166.
- Aquino de Muro M, Priest FG (1993) Phylogenetic analysis of Bacillus sphaericus and development of an ologonucleotide probe specific for mosquito-pathogenic strains. FEMS Microbiol lett 112(2): 20S-210.
- Jahnz U , Fitch A, Priest F G (1996) Evaluation of an rRNA-targeted oligonucleotide probe for the detection of mosquitocidal strains of BaciIIus sphaericus in soils: characterization of novel strains lacking toxin genes. FEMS Microbiol Eco 20(2): 91-99.
- Woodburn M, Yousten A, Hilu K (199S) Random amplified polymorphic DNA fingerprinting of mosquito pathogenic and non-pathogenic strains of Bacillus sphaericus. Int J Syst Bacteriol 4S(2): 212-217.
- Miteva V, Pobell SS (1999) Random and repetitive primer amplified polymorphic DNA analysis of BaciIIus sphaericus. J Appl Microbiol 86(6): 928-936.
- Nakamura LK (2000) Phylogeny of Bacillussphaericus-like organisms. Int J syst and Evol Micribiol S0(Pt S): 171S-1722.
- Ahmed I, Yokota A, Yamazoe A, Fujiwara T (2007) Proposal of Lysinibacillus boronitolerans gen.nov.sp. nov., and transfer of BaciIIus fusiformis to LysinibaciIIus fusiformis comb. nov. and Bacillus sphaericus to LysinibaciIIus sphaericus comb. nov. Int J syst and Evol Microbiol S7(pt S): 1117-112S.
- Dhindsa KS, Sangodkar UMX, Kumar A (2002) Novel Cost-effective method of screening soils for the presence of mosquito-pathogenic bacilli. Lett Appl Microbiol 35(6): 457-461.
- Sambrook J, Maniatis T, Fritsch EF (1989) Molecular cloning: a laboratory manual. (2nd edn), Cold Spring Harbor Laboratory Press, New York, USA, pp. 1546.
- Prabhu NN, Santimano MC, Mavinkurve S, Bhosle SN, Garg S (2010) Native granule associated short chain length polyhydroxyalkanoate synthase from a marine derived Bacillus sp. NQ-11/A2. Antonie van Leeuwenhoek 97(1): 41-50.
- Tamura K, Dudley J, Nei M, Kumar S (2007) MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 24(8): 1596-1599.
- Chan SW, Thanabalu T, Wee BY, Porter AG (1996) Unusual Amino Acid Determinants of Host Range in the Mtx2 Family of Mosquitocidal Toxins. J Biol Chem 271(24): 14183-14187.
- Otsuki K, Guaycurus TV, Vicente ACP (1997) Bacillus sphaericus entomocidal potential determined by polymerase chain reaction. Mem Inst Oswaldo Cruz 92(1): 107-108.
- Sieo CC, Wong BT, Abdullah N, Ho YW ( 2007) Effects of extraction methods and age of cells on the whole-cell protein patterns of Lactobacillus. Res J Microbiol 2(10): 727-734.
- Finney DJ (1971) Probit analysis-a statistical treatment of the sigmoid response curve, (3rd edn) Cambridge: Cambridge University Press, UK, pp. 333.
- Carboulec N, Priest FG (1989) Phenotypic characterization of some strains of Bacillus sphaericus. Lett Appl Microbia 9(3): 113-116.
- Berry C, Jackson-Tap J, Oei C, Hindley J (1989) Nucleotide sequence of two toxin genes from Bacillus sphaericus MAB59: sequence comparisons between five highly toxinogenic strains. Nucleic Acids Research 17(18): 7516.
- Aquino de Muro M, Priest F G. (1994) A colony hybridization procedure for the identification of mosquitocida1 strains of Bacillus sphaericus on isolation plates. J Invert Pathol 63(3): 310-313.
- Gordon RE, Haynes WC, Pang C H-N (1973) "The Genus Bacillus.” United States Department of Agriculture, Washington, DC, pp. 427
- Grimont F, Grimont PAD (1991) DNA fingerprinting. In Nucleic Acid Techniques in Bacterial Systematics, Edited by E. Stackebrdndt & M Goodfellow (Eds.), Chichester: John Wiley. pp. 249-219.
- Zahner V, Rabinovitch L, Cavados CFG, Momen H (1994) Multilocus enzyme electrophoresis on agarose gel as an aid to the identification of entomopathogenic Bacillus sphaericus. J Appl Bacterial 76(4): 327-335.