Received 22 January 2016; accepted 17 February 2016; published 22 February 2016
To establish sustainable agricultural system, by which crops and vegetables had been produced stably, maintaining soil fertility was primary important. In order to utilize the soil microorganisms effectively for soil management, property of soil microorganisms had to be evaluated as similarly as those of physical and chemical properties; while analyses of microbial properties in soil and manure had always included the problem that there was no available standard method to evaluate microbial property. The one of the major problems was the vast diversity and the enormous population of soil microorganisms  , the other was an existence of numerically dominant unculturable microorganisms which comprise 99% of soil habitat  .
Although denaturing gradient gel electrophoresis (DGGE), by which microbial flora could be analyzed without any culture steps  , had widely been used within this decade for environmental microbiological researches, the method was found to have the following weak points: 1) PCR bias always disturbed to provide the exact information not only of microbial number but also ratio of each microbial groups by changing a ratio of the amplified DNAs from that of the original genomes  - . 2) Microorganisms, affiliated by this method, were not a representative of whole microorganism but one part, which produced discrete bands in the DGGE gel  - , and the large portion in background smear remained unclear as “interband” region  . All of the PCR products could only been affiliated when the selective primer for the specific groups were used  -  , which resulted in another selection bias caused by PCR primers. 3) The migration rates of each bands varied depending on the prepared gel condition, which inhibited to use the past data as reference database. 4) It required a lot of time and labor not only to remove completely PCR inhibiting substances, such as humic substance contained in the sample  , but also to prepare and run DGGE gel, and affiliation of each band. These weak points disturbed the method to come into wide use. Especially as evaluation system, there was no method available by which numbers and contained microorganisms could systemically be analyzed without any preliminary information of microorganisms included in the sample.
Until now, we had found a new affiliation method of microorganisms based on restriction fragment polymorphism analysis, and developed a system and method by which bacterial affiliations could be completed systematically  . By using isolated environmental bacteria, precision of bacterial affiliation have been evaluated   . Its combined use of the most probable number method (MPN) was found useful to provide numbers and taxonomies of each bacterial group without isolation in the former papers   .
As the method also seemed suitable as simple evaluation method of microbial properties in soils and manures, we presented our evaluation results whether the method could be used as systematic analyses method of soils and manure in this manuscript.
2. Materials and Methods
2.1. Samples of Soils and Manure
Soil samples were obtained from surface of upland field (U; Gleysol) at Itoshima, Fukuoka, Japan, where vegetable had been cultivated under conventional field management, and from the surface horizon of forest soil (F; brown forest soil) at Wajiro-Hgashi, Fukuoka, Japan. Commercial manure (M) made from rice straw and cattle feces was used in this study. To test vials (5 replicates) including Biolog Universal Growth Medium (BUGM; BIOLOG Hayworth, CA, USA) broth  , serial 10-fold dilutions (10−4 to 10−10) prepared from samples (1 g fresh wt.) were inoculated. After 3 days incubation at 30˚C, bacterial DNA in each vial was extracted described previously   and purified by the conventional methods.
2.2. PCR Amplification and Restriction Digestion
Using the V2 forward primer (41f), and the V6 reverse primer (1066r)  , 16S rDNA was amplified according to the former study   . After restriction digestion of the PCR product (10 μl) by each of 10 units of the restriction enzyme, Hae III or Hha I or Rsa I or Scr F1 (Takara Bio Co. Ltd., Shiga, Japan) in Low salt buffer solution (10xLow salt buffer, Takara Bio Co. Ltd.) and 5 folds dilution by de-ionized water (for Low salt buffer), restriction fragment lengths were measured by microchip electrophoresis system (MCE-202 MultiNA; Shimadzu Co., Ltd., Kyoto, Japan).
2.3. Used Reference MERFL Database for Phylogenetic Estimation
The newly constructed database was used for this research, which was edited using the method of Watanabe and Okuda  described previously  . For 41f/1066r primers, 30,844 post-amplification sequence files, which were consisted from 1379 bacterial genera, including uncultured and unidentified bacteria, were mainly re- edited using small subunit rRNA files in RDP II release 9_61  under 5-bases mismatches in the both in primer annealing sites.
2.4. Selection of the Measured MERFLP Originated from the Homogeneous 16S rDNA and Phylogenetic Estimation
As the reference MERFL database was edited from the homogeneous 16S rDNA sequences, the measured MERFL digested from the homogeneous 16S rDNA had to be used for phylogenetic estimation.
The major RFs, which had the highest relative mole concentration (ratio of fluorescent intensity to fragment size) and represented as H in Table 1, were selected among the mixed heterogeneous FRs as described previously. The 2nd major RFs, represented as M in Table 1, were similarly selected among the remained mixed heterogeneous FRs after subtraction of the major RFs. The 3rd major gene, represented as L in Table 1, were similarly selected using the remained mixed heterogeneous FRs after subtraction of the 2nd major RFs.
The pairwise distance (DAB) between the measured RFLP (A) and the theoretical RFLP (B) was calculated according to Nei and Li  . For similarity search, the theoretical MERFLP (B) having the smallest pairwise distance (DABME), which was an average of all the DABs for used restriction enzymes, to the measured MERFLP (A) were searched in the reference database as described previously (Watanabe et al., 2008). Similarity (%) in Table 1 was calculated as the following equation; (1-DABME) × 100.
In phylogenetic estimation, identical theoretical MERFL (100%) was searched preferentially by using all the 4 measured MERFL data at first. When the completely identical theoretical MERFL was not found, combinations of 3 restriction enzymes were used for the next searches (Table 1). When the completely identical theoretical MERFL (100%) was not found, combinations of 2 restriction enzymes were used for the next searches (Ta- ble 1). When the completely identical theoretical MERFL (100%) was not found using 2 restriction enzymes, the theoretical MERFL having the highest similarity to the measured MERFL was indicated in Table 1   .
2.5. Estimation of Numbers of Each Bacterial Group by MPN
After differentiation of the measured MERFLs into 8 groups (A~J) based on the phylogenetic estimation. Numbers of each group were estimated by MPN for five-tube, three-decimal-dilution experiment (Table 2). Confidence limits shown in Table 2 were obtained using FDA’s Bacterial Analytical Manual  .
3. Results and Discussion
3.1. Grouping and Affiliation of Bacteria Based on MERFLP
Affiliations of fifty MERFLs were summarized in Table 1. The MERFLs in this study was found to have a higher bacterial diversity than those in the former studies as the followings; all of the 50 MERFLs were divided into 47 OTUs, then ratio of total number of the OTUs to that of MERFLs was 94%, which was higher than that of upland field using selective medium (62.2%)  , that of manures during composting (60.4%)  and that of commercial food products (34.6%; unpublished results). The higher diversity of MERFLs was caused from higher bacterial diversity of samples and non- selectivity of the used incubation medium. They were divided into 8 groups for the MPN calculation as the followings; Actinobacteria (Group A, 10 MERFLPs), Bacillus spp. (Group B, 5 MERFLs), the other Firmicutes (Group C, 11 MERFLs), α-Proteobacteria (Group D, 4 MERFLs), β-Proteobacteria (Group E, 8 MERFLs), γ-Proteobacteria (Group F, 7 MERFLs), δ-Proteobacteria (Group G, 3 MERFLs), and the other gram negative bacterial group (Group H, 2 MERFLs) (Table 1).
Table 1. Affiliation of bacteria grown in serially diluted BUGM medium by MERFLa.
aGrouping was based on affiliation by MERFL; Actinobacteria (Group A), Bacillus spp. (Group B), the other Firmicutes (Group C), α-Proteobacteria (Group D), β-Proteobacteria (Group E), γ-Proteobacteria (Group F), δ-Proteobacteria (Group G), and the other gram negative bacterial group (Group H); bThe 1st letter in vial indicates samples; “F” stands for the sample from forest soil, “M” stands for the sample from commercial manure, and “U” stands for the sample from upland field soil . Exponential of vial number represents the decimal dilution of the vial. The 2nd number of vial number (1 - 5) represents number in 5 replicates for the each decimal dilution. H of last letter represents MERFL originating from the major 16S rDNA, M represents from the 2nd major 16S rDNA, and L represents from the 3rd major 16S rDNA; cRestriction enzymes used for similarity search; “Ha”, “R”, “Sc”, and “Hh” stand for Hae III, Rsa I, Scr F1, and Hha I. For the measured MERFLP which had no completely identical theoretical MERFLP, the theoretical MERFLP having the highest similarity using all the RFLPs was presented with the similarity as described in the materials and method; dSpecies name (accession number) of the theoretical MERFL having the highest similarity with the measured MERFL; eAdditional name (accession number) of the theoretical MERFL using the different restriction enzymes; fDifferent accession number of the theoretical MERFL in the same group using the different restriction enzymes.
Table 2. Most probable numbers of each groups (A~H) and 5% confidence limits obtained using FDA’s Bacterial Analytical Manual  .
aGroups: A: B. cereus, B: Bacillus spp., C: Clostridium, D: The other Fumicutes, E: Actinobacteria, F: Proteobacteria, G: Prevotella, H: Cytophagales, I: Gram negative bacteria; bTotal number of bacteria.
The precision of the affiliations of each MERFLs was lower than that of the former studies. With respect to the major MERFL, ratio of the MERFLs with 100% similarity to the corresponding theoretical MERFLs (43.3%) was lower than that of field soils using selective medium (90.5%)  , that of the manures during composting (62.9%)   , and that of the commercial food products (59.6%). The lower precision of the affiliations was caused from higher bacterial diversity of sample. As the diversity of sample became higher, the each MPN vial included many kinds of bacteria, which made it difficult to select MERFLs originated from homogenous 16S rDNA.
3.2. Estimation of Numbers of Each Bacterial Groups by MPN
There was a large difference in microbial properties among the three samples as the followings. In the forest soil (F), Group E, which included Burkholderia sp., Ralstonia sp., and Alcaligenes sp., was numerically dominant bacterial group (3.64 × 106 MPN g−1 dry soil), followed by Group F (1.32 × 106 MPN g−1), Group G (0.006 × 106 MPN g−1), and Group H (0.006 × 106 MPN g−1) (Table 2, Figure 1). The some bacterial group detected here, e.g., Burkholderia spp., was reported to be detected using clone library sequencing  except for phyla Acidobacteria which could not be detected by the culture based method  .
In the commercial manure (M), Group A, which included Streptoverticillium salmonis, Mycrococcus sp., Streptomyces bikiniensis, and Microbacterium ulmi, was numerically dominant bacterial group (30.8 × 106 MPN g−1), followed by Group D (26.0 × 106 MPN g−1), which included Agrobacterim sp., Sphingomonas sp., Erythrobacter citreus and Erlichia sp., Group E (17.1 × 106 MPN g−1), which included Alcaligenes sp., and Ralstonia sp., Group G (11.2 × 106 MPN g−1), which included various sulfate reducing bacteria and Chondromyces robustus, Group C (1.71 × 106 MPN g−1), which included Paenibacillus gluconolyticus, Eubacterium cylindoides, and Staphylococcus sp., Group F (0.5 × 106 MPN g−1), and the Group H (0.05 × 106 MPN g−1) (Table 2, Figure 1).
The microbial property of M was different from those of the manures during composting in the former paper as the followings  : Total bacterial number (3.64 × 108 MPN g−1) was lower than that of the manure after thermophilic phase (7.89 × 1010 MPN g−1), and that after maturing phase (14.83 × 1010 MPN g−1). The reason of the lower number was attributed to an absence of Bacillus spp., which was the dominant bacteria in thermophilic phase, and the decrease of the other Firmicutes, sulfate reducing bacteria and the other gram negative bacterial group, which were dominant bacteria in maturing phase  . While α and β-Proteobacteria, which once disappeared after maturing phase, recovered in considerable number, and Actinobacteria, which increased during maturing phase, was remained  .
In the upland field (U), Group C, which included Paenibacillus sp., and Weissella paramesenteroides, was numerically dominant bacterial group (4.41 × 106 MPN g−1), followed by Group A (2.14 × 106 MPN g−1), which
Figure 1. Numbers of bacterial groups estimated by MPN and MERFLP in forest soil (F), commercial manure (M), and upland field (U). Number of Actinobacteria (Group A; ), the other Bacillus spp. (Group B; ), the other Firmicutes (Group C; ), α-Proteobacteria (Group D; ), β-Proteobacteria (Group E; ), γ-Proteobacteria (Group F; ), δ-Proteobacteria (Group G; ), t and the other gram negative bacterial group (Group H; ) were presented.
included Mycobacterium sp., Corynebacterium genitalium, and Arthrobacter citeus, Group B (2.14 × 106 MPN g−1), which included B.cereus, B.fusiformis/B.spaericus and B.firmus/B.smithii/B.azotoformans, and Group F (0.35 × 106 MPN g−1), which included Pseudomonas sp. (Table 2, Figure 1). Number of Bacillus spp. was similar to that of the upland Andosol field in the former study using selective medium for Bacillus spp. (1.58 × 106 MPN g−1), and those estimated by dilution plate method (2.30 × 106 CFU g−1)  . The most bacterial groups detected here were reported to be detected using culture-based method  , or DGGE  , or clone library sequencing  except for phyla Acidobacteria  .
There was a difference between the total bacterial number estimated by MPN using all the amplified vials and those of the sum of the each bacterial MPN (Table 2). In the forest soil (F), the sum of the each MPN (4.7 x106 MPN g−1) was higher than that of the total bacterial MPN (2.9 × 106 MPN g−1). The over estimation was caused from the reason that some positive vial was repeatedly counted not only as the major MERFLs, but also as the 2nd major, and the 3rd major MERFLs. In the forest soil, two 2nd major MERFLs and one 3rd major MERFL were additionally counted as the positive results to the 6 major MERFL. In the manure (M), the sum of the each MPN (74.9 × 106 MPN g−1) was lower than the total bacterial MPN (312 × 106 MPN g−1). Differentiation of the whole MERFLs into the 6 sub-groups was the major factor of the underestimation (Table 2). Because not all bacteria in each vial were detected by the method due to the PCR bias, the MPN scores of each group were lower than the true MPN scores. This under estimation could be decreased by conjugating the small sub-groups into the larger group with higher MPN score. In upland field soil (U), the sum of the each MPN (9.0 × 106 MPN g−1) was much lower than that of the total bacterial MPN (4088 × 106 MPN g−1). The under estimation was also caused from low PCR amplification rate of the numerically dominant bacteria, for which bands detected in the highest dilutions (10−8 and 10−9) were too weak to afford visible fragments after restriction digestion. The other our research indicated that amplification rates of some Antinobacteria and Firmicutes were low, and the other Actinobacteria was not amplified by the used PCR condition (unpublished results). The analysis by a new PCR condition including newly designed PCR primer for these bacteria will be presented in the following manuscripts.
In this method PCR inhibiting substances included in manure and soils had no serious effect on the results in spite of the used DNA extraction method, which included no extra purification step, because the effect of humic substances was decreased by using DNA extracted after proliferation in the growth medium, especially in higher decimal dilution vials of MPN, where the numerically dominant microorganisms were detected, concentration of the inhibiting substance was minimized. Only in forest soil, the inhibiting substance might cause the under estimation of some microbial group, because amplification band was observed until under 10−7 dilution vials and PCR inhibition was observed until 10−5 ~ 10−6 dilution vails, which afforded 6 positive vials. The under estimation might be avoided by using the conventional extraction method for environmental DNA, which included purification step.
Although the present method was culture based method, which eliminated unculturable microorganisms, we thought that the method was suitable as evaluation system of soil and manure in aim to maintain soil fertility. Because one of the unculturable microorganism, which couldn’t proliferate without the other microorganisms, was detected by this method   , and the other unculturable microorganism, which lost proliferation ability and had no effect on soil fertility, was eliminated. As to another type of unculturable microorganism, which could proliferate, but no-one knew how to proliferate them such as phyla Acidobacteria  , we already developed an unculture-based new method (unpublished result), which was found useful for unculturable microbial analysis in activated sludge. The difference of the results between culture based and unculture based method will be presented in the following manuscripts.
Classification and affiliation in species or genus level was possible by this method   . However, enumeration of each differentiated sub-groups by MPN depended on a diversity of samples; e.g., in order to enumerate in species or genus level, three-tube, three-decimal-dilution experiment was sufficient for a sample having lower microbial diversity such as food (unpublished results), while five-tube, three-decimal-dilution experiment was insufficient for a sample having higher microbial diversity such as soil or manure and ten-tube, four- decimal-dilution experiment would improve accuracy of enumeration. Because differentiation of whole MERFLs into detailed sub-groups with lower number of MERFLs caused underestimation due to PCR bias which preferentially amplified specific DNA in MPN vials   .
As the system required lower cost for instrument and running and RFLP data was automatically obtained by MultiNA, the method was suitable as evaluation system of soil and manure. Although some data processing was manually processed at this moment, the method was the versatile system used not only as evaluation system of environmental microorganisms, but also inspection method of food microorganisms (unpublished results). Compared to the next-generation method such as pyro-sequencing, reliable affiliations of all the bacteria might be difficult by our method, our method might not be suitable for pure research purpose, but suitable as inspection method due to its lower running cost and simplicity. A difference of the results obtained by this culture- based technique and by the unculture-based technique, such as DGGE, will be presented in the next manuscripts.
We thank Mr. Y. Sogabe, Global Application Center, Shimadzu Co., for variable suggestion and support for MultiNA. We thank Dr. H. Yosikawa, the former Prof. of Fukuoka Institute of Technology, Dr. A. Hosoda, and Prof. H. Tamura, Meijyo University, for their encouragements during this work.
 Muyzer, G., Waal, E.C.D. and Uitterlinden, A.G. (1993) Profiling of Complex Microbial Populations by Denaturing Gradient Gel Electrophoresis Analysis of Polymerase Chain Reaction-Amplified Genes Coding for 16S rRNA. Applied and Environmental Microbiology, 59, 695-700.
 Ishii, K. and Fukui, M. (2001) Optimization of Annealing Temperature to Reduce Bias Caused by a Primer Mismatch in Multitemplate PCR. Applied and Environmental Microbiology, 67, 3753-3755.
 Neilson, J.W., Jordan, F.L. and Maier, R.M. (2013) Analysis of Artifacts Suggests DGGE Should Not Be Used for Quantitative Diversity Analysis. Journal of Microbiological Method, 92, 256-263.
 Direito, S.O.L., Zaura, E., Little, M., Ehrenfreund, P. and Roling, W.F.M. (2014) Systematic Evaluation of Bias in Microbial Community Profiles Induced by Whole Genome Amplification. Environmental Microbiology, 16, 643-657.
 Griffiths, R.I., Whiteley, A.S., O’Donnell, A.G. and Bailey, M.J. (2000) Rapid Method for Coextraction of DNA and RNA from Natural Environments for Analysis of Ribosomal DNA- and rRNA-Based Microbial Community Composition. Applied and Environmental Microbiology, 66, 5488-5491.
 Duineveld, B.M., Kowalchuk, G.A., Keijzer, A., van Elsas, J.D. and van Veen, J.A. (2001) Analysis of Bacterial Communities in the Rhizosphere of Chrysanthemum via Denaturing Gradient Gel Electrophoresis of PCR-Amplified 16S rRNA as Well as DNA Fragments Coding for 16S rRNA. Applied and Environmental Microbiology, 67, 172-178.
 Sun, H.Y., Deng, S.P. and Raun, W.R. (2004) Bacterial Community Structure and Diversity in a Century-Old Manure-Treated Agroecosystem. Applied and Environmental Microbiology, 70, 5868-5874.
 Huang, J., Sheng, X., He, L., Huang, Z., Wang, Q. and Zhang, Z. (2013) Characterization of Depth-Related Changes in Bacterial Community Compositions and Functions of a Paddy Soil Profile. FEMS Microbiology Letters, 347, 33-42.
 Nikolausz, M., Sipos, R., Revesz, S., Szekely, A. and Marialigeti, K. (2005) Observation of Bias Associated with Re-Amplification of DNA Isolated from Denaturing Gradient Gels. FEMS Microbiological Letters, 244, 385-390.
 Santegoeds, C.M., Ferdelman, T.G., Muyzer, G. and Beer, D.D. (1998) Structural and Functional Dynamics of Sulfate-Reducing Populations in Bacterial Biofilms. Applied and Environmental Microbiology, 64, 3731-3739.
 Watanabe, T., Asakawa, S., Nakamura, A., Nagaoka, K. and Kimura, M. (2004) DGGE Method for Analyzing 16S rDNA of Methanogenic Archaeal Community in Paddy Field Soil. FEMS Microbiology Letters, 232, 153-163.
 Chu, H., Fujii, T., Morimoto, S., Lin, X., Yagi, K., Hu, J. and Jiabao, Z. (2007) Community Structure of Ammonia-Oxidizing Bacteria under Long-Term Application of Mineral Fertilizer and Organic Manure in a Sandy Loam Soil. Applied and Environmental Microbiology, 73, 485-491.
 Watanabe, K. and Koga, N. (2009) Use of a Microchip Electrophoresis System for Estimation of Bacterial Phylogeny and Analysis of NO3- Reducing Bacterial Flora in Field Soils. Bioscience Biotechnology and Biochemistry, 73, 479-488.
 Watanabe, K., Okuda, M. and Koga, N. (2008) A Newly Developed System Based on Multiple Enzyme Restriction Fragment Length Polymorphism—An Application to Proteolytic Bacterial Flora Analysis. Soil Science and Plant Nutrition, 54, 204-215.
 Watanabe, K., Horinishi, N. and Matumoto, K. (2015) Antibiotic-Resistant Bacterial Group in Field Soil Evaluated by a Newly Developed Method Based on Restriction Fragment Length Polymorphism Analysis. Advances in Microbiology, 5, 807-816. http://dx.doi.org/10.4236/aim.2015.512085
 Watanabe, K., Horinishi, N., Matumoto, K., Tanaka, A. and Yakushido, K. (2015) Bacterial Groups Concerned with Maturing Process in Manure Production Analyzed by a Method Based on Restriction Fragment Length Polymorphism Analysis. Advances in Microbiology, 5, 832-841. http://dx.doi.org/10.4236/aim.2015.513088
 Weidner, S., Arnold, W. and Puhler, A. (1996) Diversity of Uncultured Microorganisms Associated with the Seagrass Halophila stipulacea Estimated by Restriction Fragment Length Polymorphism Analysis of PCR-Amplified 16S rRNA Genes. Applied and Environmental Microbiology, 62, 766-771.
 Cole, J.R., Chai, B., Farris, R., Wang, Q., Kulam-Syed-Mohideen, A.S., McGarrell, D.M., Bandela, A.M., Cardenas, E., Garrity, G.M. and Tiedje, J.M. (2007) The Ribosomal Database Project (RDP-II): Introducing myRDP Space and Quality Controlled Public Data. Nucleic Acids Research, 35, D169-D172.
 Nei, M. and Li, W.H. (1979) Mathematical Model for Studying Genetic Variation in Terms of Restriction Endonucleases. Proceedings of the National Academy of Sciences of the United States of America, 76, 5269-5273.
 Tsai, S.H., Selvam, A., Chang, Y.P. and Yang, S.S. (2009) Soil Bacterial Community Composition across Different Topographic Sites Characterized by 16S rRNA Gene Clones in the Fushan Forest of Taiwan. Botanical Studies, 50, 57-68.
 Naether, A., Foesel, B.U., Naegele, V., Wust, P.K., Weinert, J., Bonkowski, M., Alt, F., Oelmann, Y., Polle, A., Lohaus, G., Gockel, S., Hemp, A., Kalko, E.K.V., Linsenmair, K.E., Pfeiffer, S., Renner, S., Schoning, I., Weisser, W.W., Wells, K., Fischer, M., Overmann, J. and Friedricha, M.W. (2012) Environmental Factors Affect Acidobacterial Communities below the Subgroup Level in Grassland and Forest Soils. Applied and Environmental Microbiology, 78, 7398-7406.
 Tang, H., Xiao, C., Ma, J., Yu, M., Li, Y., Wang, G. and Zhang, L. (2009) Prokaryotic Diversity in Continuous Cropping and Rotational Cropping Soybean Soil. FEMS Microbiology Letters, 298, 267-273.
 Bae, J.W., Kim, J.J., Jeon, C.O., Kim, K., Song, J.J., Lee, S.G., Poo, H., Jung, C.M., Park, Y.H. and Sung, M.H. (2003) Application of Denaturing Gradient Gel Electrophoresis to Estimate the Diversity of Commensal Thermophiles. Journal of Microbiology and Biotechnology, 13, 1008-1012.