Tubagus Haeru.pmd Volume 3, Number 2, August 2009 p 56-60 ISSN 1978-3477 *Corresponding author, Phone/Fax: +62-21-7890962, E-mail: haeru_tb@yahoo.com The growing economic importance of the aquaculture industry worldwide and the need to find a new approach for a probiotic method has led to increasing interest in rapid and reliable methods for detection and identification of bacteria (Nilsson and Strom 2002). Recently, nucleic acid-based assays have been studied intensively. One of methods being developed is that of 16S ribosomal RNA genes sequencing. These assays are generally considered to be more accurate, more sensitive, inexpensive and more rapid to perform than those of more conventional variety, such as morphological, biochemical/physiological and serological assays, which led to various results of species resemblance in physiology and susceptibility to environmental exchanges (Macrae 2000). Broad-range PCR amplification of 16S rDNA with universal bacterial primers and subsequent sequence analysis of cloned products is a widely used method in molecular biology. The method can be used for identification of microorganism identity to the genus or species level (Tiirola et al. 2002) and construct phylogenetic analysis which plays on important role in taxonomic positioning. It can be undertaken by comparing the gene sequences with other bacteria sequences in the GenBank database using basic local alignment search tool (BLAST) search and continue with method of neighbor-joining to construct a phylogenetic tree (Saitou and Nei 1987; Holmes 2003). In this paper, we identified the selected bacterial isolates from L. vannamei culture system and gut environment using the 16S ribosomal RNA gene sequence and constructed and analyzed the phylogenetic trees using the program Clustal X. It is important to understand that genetic diversity of selected isolates of probiotic bacteria and their taxonomic position provides basic information for developing disease control strategies in the shrimp culture industry. MATERIALS AND METHODS Sampling of Bacteria and Isolation. Bacterial isolates were obtained from the water, sediments and shrimp’s gut using healthy wild shrimps (15 ± 1 g) collected from four districts, i.e., Pandeglang, Serang, Tangerang and Karawang. Sampling stations were determined randomly with 5 to 7 spots per pond. Substrates were taken from 3 substrates, i.e. pond water, mud and shrimp’ gut following the procedures of Hjelm et al. (2004). Water aliquots of 500 mL and mud aliquot of 50 to100 g were taken using sterile plastic bags and stored in refrigerator at 4 °C. Samples from shrimp’ guts were conducted by taking living shrimps to the laboratorium using closed transportation with air and water at a ratio 3:1. Sample of 5% of obtained material were taken and then isolated and enriched on TSB medium with 2% NaCl and incubated for 24 h at 30 °C before streaking bacteria on marine agar. Samples from mud of 50% (w/v) were diluted using sterile distilled water, whereas samples from shrimp’s guts were taken by gut-scraping. All subsequent isolates were streaked on marine agar and incubated for 24 h at 30 °C. Representative colonies were selected based on phenotype, Gram staining and acidity of the medium. The selected isolates were cultured on slash- marine-agar and stored at 4-10 °C after separation into three groups, i.e. original culture, stock culture and working culture. Development of Pure Cultures of Isolates. Nineteen selected isolates from the slash-agar stock were streaked in Identification and Phylogenetic Analysis of Bacterial Isolates from Litopenaeus vannamei Shrimp Culture System and Gut Environment Based on 16SrRNA Gene Sequence Data TUBAGUS HAERU RAHAYU 1,2*, INDRAWATI GANDJAR1, ETTY RIANI3, IIN SITI DJUNAIDAH4, AND WELLYZAR SJAMSURIDZAL1 1Department of Biology, Faculty of Mathematics and Natural Sciences, Universitas Indonesia, Depok 16424, Indonesia; 2Department of Aquaculture, Sekolah Tinggi Perikanan Jakarta; Jalan AUP Pasar Minggu, Jakarta 12520, Indonesia; 3Faculty of Fisheries and Marine Science, Institut Pertanian Bogor, Kampus Darmaga, Bogor 16680, Indonesia; 4Agency for Human Resources Development of Marine and Fisheries, Departemen Kelautan dan Perikanan, Jalan MT. Haryono Kav. 52-53, Jakarta 12770, Indonesia Selected bacterial isolates from a Litopenaeus vannamei shrimp culture system and gut environment were assessed using 16S rRNA gene sequencing method to identify their identity and to construct their phylogenetic relationship. In a preliminary study, a total of 19 isolates were selected as probiotics. These isolates were prepared using freeze and heat-shock method to obtain the DNA template. PCR amplification of 16S ribosomal RNA gene of isolates was carried out using bacterial universal primers 9F and 1510R and was sequenced using an automated DNA sequencer. These gene sequences were compared with other gene sequences in the GenBank database (NCBI) using a BLAST search to find closely related sequences. Alignment of these sequences with sequences available from GenBank database was carried out to construct a phylogenetic tree for these bacteria. Most of the isolates obtained, i.e. 17 out of the 19 isolates, belonged to different species of Bacillus, sharing 95 to 99% 16S ribosomal RNA identity with the respective type-strain, whereas the remaining 2 isolates belonged to Micrococcus sp. and Micrococcus luteus, with 97 to 99% 16S rRNA homology, consecutively. Key words: Bacillus, Micrococcus, PCR amplification of 16S rRNA gene, phylogeny Microbiol Indones 57 duplicate on the respective media (marine agar, DIFCO, 55.1 g.l-1) plates and incubated at 30 ºC for 48 h. The isolates were then stored in a refrigerator for short periods until further use. Preliminary Characterization of Isolated Bacteria. Characterization of isolates was done following the method of Akhtar et al. (2008). Colonies on respective media plates were examined using the light microscope and their characteristics i.e. color, form, elevation and margin of colonies were recorded. Isolates were also characterized on the basis of Gram staining. DNA Template Preparation. DNA templates were prepared by using the freeze and heat-shock method (Sjamsuridzal and Oetari 2003). Incubated isolates (48 h) were taken using sterile microtips and put into 1.5 mL microtube containing 500 µL of nuclease-free water. The cell suspension was vortexed for 1 min and then frozen for 24 h. The suspension was then boiled using a water bath (>95 ºC) for 20 min. The cells suspension was centrifuged at 13 000 rpm for 15 min. The supernatant containing the DNA was stored in a refrigerator until further use. Phylogenetic Analysis of Isolates. Polymerase chain reaction amplification of 16S ribosomal RNA of isolates were carried out using following primers and set of conditions: forward primer (9F; 5'-GAGTTTGATCCTGGCTCAG-3' and reverse primer (1510R; 5’GGCTACCTTGTTACGA-3') (Nilsson and Strom 2002). Each vial contained 1.25 µL of each primer, 15 µL RTG suspended with nuclease free water and 7.5 µL DNA template in a total 25 µL reaction volume. The reaction mixture was heated for 3 min at 95 °C and then amplification was carried out in 35 cycles. Each cycle was comprised of 30 sec at 95 °C, 15 sec at 55 °C and 1 min at 72 °C. The final extension was for 5 min at 72 °C followed by storage at 25 °C. The PCR product was purified using sodium acetate/ethanol precipitation to remove excess primers and free nucleotides. The next step was cycle-sequencing reaction, in which a similar process to amplification of 16S ribosomal RNA was performed. Each vial contained 1 µL Big Dye Terminator Ready Reaction Mix V.3.1, 0.5 µL primer 9F, 7.0 µL sequence buffer, 0.5 µL nuclease free water and 1 µL DNA template in a total of 10 µL reaction volume. The reaction mixture was heated for 1 min at 96 °C and then amplification was carried out over 30 cycles. Each cycle was comprised of 10 sec at 96 °C, 5 sec at 55 °C and 90 sec at 60 °C. The final step was storage at 25 °C. The PCR product was re-purified using sodium acetate/ethanol precipitation to remove excess dye- terminators, followed by denaturation using 15 µL Hi-DI Formamide and heat denaturation at 95 °C for 5 min and then immediately transferred on ice. Amplified 16S ribosomal RNA fragments of about 1500 bp size were obtained from each isolate and were partially sequenced (around 500 bp) using an automated DNA sequencer (ABI 310 PRISM), following the procedure given by the manufacturer. The gene sequences were compared with others in the GenBank databases using the NCBI BLAST (www.ncbi.nlm.nih.gov). Gene sequences of 16S ribosomal RNA of selected organisms were obtained from GenBank and aligned with gene sequence of our isolates using a multiple sequence alignment program CLUSTAL X software (Thompson 1997). The aligned sequences were used to construct a distance matrix, after the generation of 100 bootstrap sets, which was subsequently used to construct a phylogenetic tree using the neighbor-joining method (Holmes 2003). RESULTS Colony and Cell Morphologies of Isolates. Nineteen isolates of selected putative probiotic bacteria were examined for studies on the basis of the colony color and bacterial morphology (Table 1). Colony pigmentation of the isolates included white, light pink, pink-orange, light yellow, yellow and pinkish red. These colonies were dominated by circular forms and only 4 isolates had circular-tending to a irregular form. Colony elevation was convex in all cases. All isolates had entire margins and gave Gram positive staining. PCR Amplification of 16S rRNA Genes. Electrophoretic visualization of amplified 16S rRNA genes from nineteen selected probiotics using primer 9F and 1510R is shown in Table 1 Some characteristics of different isolates of selected putative probiotic bacteria Isolates code S 18 T 1 T 17 T 28 P 18 T 9 P 43 T 21 T 26 K 48 P 10 P 16 P 11 S 9 T 23 S 23 T 18 K 52 P 7 Substrate water water gut gut gut gut water gut gut sediment sediment water water gut gut sediment water water sediment Color White White White White White Light Pink Light Pink Pink Orange Pink Orange Pink Orange Pink Orange Pink Orange Pink Orange Light yellow Yellow Pinkish red Pinkish red Pinkish Red Pinkish Red Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Glossy Surface Forms Circular Circular Circular Circular Circular Circular Circular Circular Circular Circular tend to Irregular Circular tend to Irregular Circular tend to Irregular Circular tend to Irregular Circular Circular Circular Circular Circular Circular Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Entire Margins Elevation Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Convex Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Positive Gram staining Macroscopic observation of colony S18 until P1 shown the isolates code of selected indigenous bacteria Volume 3, 2009 Microbiol Indones Fig 1. The characteristic bacterial 16S gene amplicon of about 1500 bp was clearly visible in all cases. Identity of Isolates Based on 16S-rRNA-Sequence Method. Nucleotides BLAST search of partial sequences (500 bp) of the 16S rRNA for various isolates showed that most of the isolates had variable percentage identity with Bacillus (95-99%), whereas the remainder, i.e. isolate S 23 and S 9, had 97% and 99% homology to Micrococcus sp. and Micrococcus luteus, respectively (Table 2). Construction of Phylogenetic Tree. A phylogenetic tree was developed by aligning 16S rRNA sequences of different bacteria taken from GenBank, NCBI and sequences of the isolates from this study (Fig 2). Isolate S 9 was related to Micrococcus luteus clone B14. It shared 99% similarity in its 16S ribosomal RNA sequence, whereas isolate T 23 was related to Micrococcus sp. Mali2 and Kocuria sp. with 97% similarity. Based on the morphology of colony (Table 1), both isolates had similarity in patters of surface, form, margin, elevation and Gram staining i.e., glossy, circular, entire, convex and Gram positive, consecutively. The defining feature of both isolates was that isolate S 9 had a light yellow color, whereas T 23 was yellow. Meanwhile, the cluster of the Bacillus group was divided into six clusters. The first cluster consisted of isolates P 18 , S 18 , T 28 , T 1 and T 17 , which were closely related to Bacillus pumilus (DURCK14, S6-05, AU39 and NUC-F) and Bacillus antrachis (me-12), with similarity of 98, 98, 99, 99, 95%, consecutively, and bootstrap value 99.9%. The second cluster consisted of two isolates, i.e. T 21 and T 26 . Both isolates related to Bacillus sp. WRB-4, with similarity 95 and 97% consecutively with bootstrap value of 100%. The third cluster was divided into two sub-clusters. The first was isolate P 7 , which was related to Bacillus sp. Mali10 with similarity of both 97% while the second sub-cluster was isolates S 23 and K 52 , which were close to Bacillus sp. By231Ydz-fq with a similarity of 96-97%. The fourth cluster was isolate T 18 . This was close to Bacillus sp. K38T with a 99% similarity. The fifth cluster was represented by T 9 and P 43 which were related to Bacillus subtilis S8-04 and Bacillus megaterium TK1 with a similarity of 96 and 98% consecutively. The sixth cluster was again divided into sub-clusters. The first sub-cluster consisted of isolate P 10 with a 99% similarity and related to Bacillus flexus LLH. The second sub-cluster was represented by isolates K 48 , P 16 and P 11 , and related to Bacillus flexus GS11 and Bacillus sp. Bsi20565 with similarities of 95, 97, 98 and bootstrap values higher than 60%. DISCUSSION This is a preliminary study on the characteristics of nineteen isolated colonies. Colonies were grown on agar media and were visible after 12 h incubation at 30 °C. This is due to the presence of sufficient nutrient, appropriateness of the agar used and appropriateness of the culture period. Marine agar (DIFCO; 55.1 g l-1) is a specific agar for heterotrophic marine bacteria. It contains 2% (w/v) NaCl. Such conditions were appropriated for isolates which were obtained from shrimp ponds with a salinity of about 25 ppt. Thus the physiological activity from the bacteria was not disrupted. The 48 h-incubation-period was adequate enough for bacterial growth. Our study showed the colonies grew Table 2. Nearest relative of various selected isolates of probiotic bacteria SubstrateIsolate S 18 T 1 T 17 T 28 P 18 T 9 P 43 T 21 T 26 K 48 P 10 P 16 P 11 S 9 T 23 S 23 T 18 K 52 P 7 water water gut gut gut gut water gut gut sediment sediment water water gut gut sediment water water sediment Nearest relative Bacillus pumilus isolate NUC-F Bacillus pumilus strain S6-05 Bacillus pumilus strain AU39 Bacillus pumilus strain DURCK14 Bacillus pumilus isolate EGU275 Bacillus subtilis strain S8-04 Bacillus megaterium strain TK1 Bacillus sp. WRB-4 Bacillus sp. WRB-4 Bacillus flexus strain GS11 Bacillus flexus ISOLAT LLH Bacillus flexus strain GS11 Bacillus sp. BSi20565 Micrococcus luteus clone B14 Micrococcus sp. “Mali2” Bacillus sp. By231Ydz-fq Bacillus sp. K38T Bacillus sp. ‘Mali10’ Bacillus sp. ‘Mali10’ DQ833752 E U 6 2 4 4 2 9 EF032679 A M 7 7 8 1 9 1 EF633222 E U 6 2 0 4 1 2 E U 5 8 6 0 3 4 EF636891 EF636891 DQ365587 DQ333292 DQ365587 E U 3 3 0 3 4 1 E U 1 9 6 5 3 1 AY 2 11 0 9 6 E U 0 7 0 3 7 2 A M 9 8 3 5 2 5 AY 2 111 0 4 AY 2 111 0 4 Accession Number Homology (%) 9 8 9 9 9 5 9 9 9 8 9 6 9 8 9 5 9 7 9 5 9 9 9 7 9 8 9 9 9 7 9 6 9 9 9 7 9 7 S18 until P1 shown the isolates code of selected indigenous bacteria Fig 1 Visualization of amplification of 16S rRNA genes from 19 selected using primer 9F and 1510R. M, Marker; bp, base pair; S18 until P1 shown the code of selected indigenous bacteria isolates. 500 bp M S18 T 1 T 17 T 28 P18 T 9 P43 T 21 T 26 K48 P10 P16 P11 S9 T 23 S23 T 18 K52 P 7 ¢ 58 RAHAYU ET AL. Microbiol Indones 59 Fig 2 The inferred relationship based on partial 16S ribosomal RNA sequence, of selected putative probiotic bacteria (boldface) to other bacteria. The tree was rooted with Escherichia coli as an out-group. Scale bar represents the number of inferred nucleotide substitution per site. Bootstrap values (1000 replicates) are shown at the nodes. well, with the typical color of bacteria. The results were in line with Akhtar et al. (2008), in that the growth of bacteria will also be determined by the period of incubation by giving enough time to reach the optimum stage. On a molecular basis, all selected isolates belonged to the Gram positive group. Most of the isolates (17 isolates) were bacilli (89.7%), with the remaining two isolates (10.53%) being Micrococci (Table 1). It indicated that bacteria belonging to Bacilli group were more prevalent in healthy L. vannamei shrimp ponds. This finding was in line with those of Verschuere et al. (2000), Balcázar et al. (2006) and Geovanny et al. (2007) in which several bacteria were identified and can be considered as probiotic bacteria in aquaculture such as: Bacillus, Vibrio, Pseudomonas, Micrococcus, Bacteroides and Clostridium. Gram negative bacteria were not found in the 19 culture isolates, after identification based on the sequence of 16S ribosomal RNA. It was probably due to the fact Gram negative bacteria which were successful obtained during sampling of bacteria and isolation were not showing the characteristics of probiotic bacteria i.e. do not improve shrimps’ survival rate during challenge test (Rengpipat et al. 2000). Isolates related to Bacilli were represented by Bacillus sp. (isolates code: T 21 , T 26 , T 18 , K 52 , P 7 , P 11 and S 23 ), B. pumilus (isolates code: S 18 , T 1 , T 17 , T 28 and P 18 ), B. subtilis (isolate code: T 9 ), B. megaterium (isolate code: P 43 ), and B. flexus (isolates codes: K 48 , P 10 and P 16 ). Meanwhile, Micrococcus, was represented by Micrococcus sp. (isolate code: T 23 ) and M. luteus (isolate code: S 9 ). As the major species found, the dominance of Bacilli can be comprehended since they have a wide range of diversity in physiological ability with respect to heat, acidity and salinity tolerance. bacilli can be found in water, soil and the digestive tracts of several organisms, including shrimps (L. vannamei) (Rengpipat et al. 2000). Although the genus Micrococcus has not a wide range of tolerance compared to Bacillus, it is however considered Volume 3, 2009 0.05Knuc 1 0 0 0 1 0 0 0 9 7 3 9 9 9 S 9 T23 9 9 9 Micrococcus luteus clone B14 (EU196531) Kocuria sp G3DM 46-EU037282 Micrococcus sp Mali2 (AY211096)7 3 3 P 1 8 S 1 8 T1 Bacillus pumilus DURCK14 (AM778191) T28 Bacillus pumilus strain AU39 (EF032679) Bacillus anthracis-me-12 (EU6520) Bacillus pumilus isolate NUC-FDQ833752 Bacillus pumilus strain S6-05 (EU624429) Bacillus pumilus isolate EGU275 (EF633222) T17 9 9 9 5 5 4 1 0 0 0 T9 Bacillus sp.WRB-4 (EF636891) T21 T26 9 5 0 9 9 7 8 8 0 5 9 2 9 9 9 P 7 Bacillus sp Mali10 (AY211104) Bacillus sp By231Ydz-fq (EU070372) S 2 3 K52 Bacillus sp K38T (AM983525) T18 9 8 0 9 9 1 6 2 3 7 3 8 P 4 3 Bacillus subtilis (EU620412) Bacillus megaterium (EU586034) 6 5 3 6 5 2 Bacillus flexus LLH (DQ333292) P 1 0 6 4 6 Bacillus sp BSi20565 (EU330341) 9 1 3 Bacillus flexus strain GS11(DQ365587) 6 5 4 K48 P 1 6 P11 OUT-GROUP Escherichia coli strain SVUB1(AM982782) Microbiol Indones60 RAHAYU ET AL. as a halotolerant group. Micrococcus can grow at 5% (w/v) NaCl, and manage to survive in the digestive tract as well. Geovany et al. (2007) reported that M. luteus was the major component of the population of the gut of Rainbow Trout (Oncorhynchus mykiss). Phylogenetic analysis indicated that nineteen strains isolated from selected local probiotic bacteria can be identified as Gram positive bacteria of the Bacillus and Micrococcus species (Fig 2). When the type of Escherichia coli used for an out-group, it was clearly seen that the distribution of the nineteen of the isolates fell into two big clusters, i.e. Bacillus cluster and Micrococcus cluster. One was composed of seventeen Bacillus genera consisted of isolate P 18 , S 18 , T 28 , T 1 , T 17 , T 21 , T 26 , P 7 , S 23 , K 52 , T 18 , T 9 , P 43 , P 10 , K 48 , P 16 and P 11 , indicated 95% in bootstrap value. The other was composed of two micrococcal strains, S 9 and T 23 , indicated 100% in bootstrap value. The higher bootstrap value (≥95%) indicated that both groups were clearly distinct (Meerak et al. 2007). Based on the findings above, it may be concluded that the 19 isolates of bacteria can be used as probiotic bacteria and are useful for developing shrimp culture systems. Bacillus and Micrococcus are frequently used in shrimp aquaculture to improve the health status toward disease outbreaks through diet or giving out into water culture (Gatesoupe 1999; Verschuere et al. 2000; Balcázar et al. 2006). ACKNOWLEDGEMENT This research was supported partly by a Local Graduate Scholarship, awarded by the Ministry of Marine Affairs and Fisheries, and partly by the Research and Development for Public Affairs Unit of Sekolah Tinggi Perikanan. REFERENCES Akhtar N, Ghauri MA, Iqbal A, Anwar MA, Akhtar K. 2008. Biodiversities and phylogenetic analysis of culturable bacteria indigenous to Khewra salt mine of Pakistan and their industrial importance. Braz J Microbiol 39:143-50. Balcázar JL, Blas I de, Zarzuela IR, Cunningham D, Vendrell D, Múszquiz JL. 2006. The role of probiotics in aquaculture. Vet Microbiol 114:173-86. Hjelm M, Bergh Ø, Riaza A, Nielse J, Melchiorsen J, Jensen S, Duncan H, Ahrens P, Birkbeck H, Gram L. 2004. Selection and identification of autochthonous potential probiotic bacteria from Turbot Larvae (Scopthalmus maximus) rearing units. Syst Appl Microbial 27:360-71. Gatesoupe FJ. 1999. The use of probiotics in aquaculture. A review. Aquaculture 180:147-65. Geovany GR, Luis BJ, Shen MA. 2007. Probiotics as control agents in aquaculture. A review. J Ocean Univ China 6:76-9. Holmes S. 2003. Bootstrapping phylogenetic trees: theory and methods. Stat Sci 18:241-55. Macrae A. 2000. The use of 16S rDNA methods in soil microbial ecology. Braz J Microbiol 31:77-82. Meerak JHL, Watanabe Y, Miyashita M, Sato H, Nakagawa Y, Tahara Y. 2007. Phylogeny of γ-polyglutamic acid-producing Bacillus strains isolated from fermented soybean foods manufactured in Asian countries. J Gen Appl Microbiol 53:315-23. Nilsson WB, Strom M S. 2002. Detection and identification of bacterial pathogens of fish in kidney tissue using terminal restriction fragment length polymorphism (T-RFLP) analysis of 16S RNA genes. Dis Aqua Org 48:175-85. Rengpipat S, Rukpratanporn S, Piyatiratitivorakul S, Menasaveta P. 2000. Immunity enhancement in black tiger shrimp (Penaeus monodon) by a probiont bacterium (Bacillus S11). Aquaculture 191:271-88. Saitou N, Nei M. 1987. The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406-25. Sjamsuridzal W, Oetari A. 2003. Rapid preparation of fungal and bacterial genomic DNA for PCR. Hayati 10:122-4. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG. 1997. The Clustal X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 24:4876-82. Tiirola M, Valtonen ET, Kinnunen PR, Kulomaa MS. 2002. Diagnosis of flavobacteriosis by direct amplification of rRNA genes. Dis Aqua Org 51:93-100. Verschuere L, Rombaut G, Sorgeloos P, Verstraete W. 2000. Probiotic bacteria as biological control agents in aquaculture. Microbiol Mol Biol Rev 64:655-71.