Grain legumes play an important nutritional role in the diet of millions of people in the developing countries . Black gram (Vigna mungo L. Hepper) is a short duration crop that belongs to grain legumes family and rich in protein . It is an important legume crop in Asia . In Myanmar, it is one of the major exportable crops and it is the second largest cultivated Legume crop as well . During 2017-2018, black gram growing areas in Myanmar was about 9.77 million hectares with a production of 1.37 metric tons, and the average yield was 1.41 metric ton∙ha−1 .
Nitrogen (N) fixation through legume-Rhizobium symbiosis is important for enhancing agricultural productivity and is therefore of great economic interest . Leguminous crops have a reputation for maintaining soil fertility since it can assimilate nitrogen from the atmosphere through symbiotic biological N2 fixation (BNF) with Rhizobia . Biological nitrogen fixation (BNF) is an important component of sustainable agriculture , and rhizobial inoculants have been applied frequently as biofertilizers.
Rhizobial inoculant can be used to substitute the nitrogenous fertilizers in food legume crops. Recently, peat-based root nodule bacterial inoculants containing TAL strains are using as biofertilizer in seven legumes distributed by Ministry of Agriculture and Irrigation, Myanmar . However, the effectiveness of symbiotic N2 fixation may be an important factor to take into consideration through successful management of symbiosis between black gram cultivars and effective bradyrhizobial strains.
Commercial production of Rhizobium inoculant has been developed for several decades. The main objective of using Rhizobium inoculant is to substitute the nitrogenous fertilizers in food legume production. It is cheaper and lighter in weight than urea and easier to use for the farmers . To improve high quality inoculant production, the evaluation of highly effective rhizobial strains for specific legume is one of the principle obligations  . Indigenous rhizobial strains also play an important role since they have adapted to local environmental conditions. Therefore, the investigation of effective indigenous rhizobial strains should be considered for current inoculant production in Myanmar.
Myanmar farmers have used, and continue to use, rhizobial inoculants when sowing legumes, but the practice is currently not extensive. The Department of Agricultural Research, Ministry of Agriculture, Livestock and Irrigation is responsible for producing inoculants in Myanmar for their distribution to farmers. Production by DAR peaked during the 1980s at 600 - 700,000 packets annually. Current production is < 100,000 packets, due to limitations in the whole supply chain from production and quality assurance to distribution to demand. Myanmar farmers use nitrogenous (N) fertilizers sparingly, particularly on legume crops. Thus, low-nodulation induced N deficiencies of the legumes are not remedied by inputs of fertilizer N and the value of lost production could exceed $100 million annually .
Several studies also reported significant increase growth parameters and yield due to the inoculation of rhizobial isolates in chickpea  in mung bean and in soybean     in Myanmar to identify the strains, plant growth and nitrogen fixation. However, there was nobody investigated the native bradyrhizobial isolates from Myanmar black gram cultivars on diversity, plant growth and N2 fixation. So, this research is very first report for phylogenetic diversity and evaluation the effectiveness of indigenous Bradyrhizobium strains for Myanmar black gram (Vigna mungo L.) cultivars and this research finding is needed in order to guarantee the potential of these indigenous Bradyrhizobium strains for black gram production. Not only bradyrhizobial inoculation is important for black gram cultivation by poor Myanmar farmers as valuable cheap source of nitrogen enough N2 fixation by biological means are also necessary and attractive topic for investigation.
For this reason, we aimed to isolate indigenous root nodule bacteria from collected soil samples of major black gram growing areas of Myanmar, to identify the phylogenetic diversity of indigenous black gram-nodulating bradyrhizobia in Myanmar based on sequence analysis of the 16S rRNA region of the isolates, and to evaluate the effectiveness of indigenous Myanmar Bradyrhizobium strains for plant growth and nitrogen fixation of Myanmar black gram varieties are required for investigation.
2. Materials and Methods
2.1. Origin of Soil Sample Collection Sites to Obtain Black Gram Bradyrhizboial Isolates
Soil samples were collected from eight major black gram growing regions of Myanmar (Figure 1). Meadow soil occurs in Delta Plains Region of Nyaunglebin Bago Region (17˚57'N 96˚43'E) with a pH of 5.04, Costal Strips Region of Chaungzon Mon State (16˚23'N 97˚32'E) with a pH of 4.69 and Sittwe Rakhine State (20˚30'N 93˚20'E) with a pH of 5.91. The Costal Strips Region of Launglon Tanintharyi Region (14˚05'N 98˚12'E) with a pH of 4.79 and Mountain and Hill
Figure 1. Map of Myanmar showing soil sampling sites from black gram major growing areas in Myanmar. Numbers in parentheses represent the total number of isolates collected from each site.
Region of Hpa-an Kayin State (16˚53'N 97˚38'E) with a pH 5.13 was Latritic Soils . The Central Regions of Pyinmanar Nay Pyi Taw Region (19˚45'N 96˚12'E) was a Meadow Alluvial Soil with a pH of 6.01 and Salingyi Sagaing Region (21˚58'N 95˚05'E) with a pH 6.97 was Meadow and Meadow Alluvial Soil are performing a tropical wet and dry climate . Alluvial Soil occurs in Delta Plains Region of Danubyu Ayeyarwady Region (17˚22'N 95˚27'E) with a pH of 6.70 (Table 1, Table 2 and Figure 1).
Table 1. Soil classification, location and climate of soil samples collected from major black gram growing areas of Myanmar.
Sources: 19)Shein, H.A. The soil types and characteristics of Myanmar. Department of Agriculture, Ministry of Agriculture, Livestocks and Irrigation: Nay Pyi Taw, Myanmar, 2015. 20)Aung, L.L.; Zin, E.E.; Theingi, P.; Elvera, N.; Aung, P.P.; Han, T.T.; Oo, Y.; Skaland, R.G. Myanmar Climate Report published by Department of Meteorology and Hydrology Myanmar, Ministry of Transport and Communications, Government of the Republic of the Union of Myanmar, 2017.
Table 2. Characteristics of soils used to obtain black gram rhizobia isolate from major black gram growing regions of Myanmar.
Soil sample analyses were performed in the plant nutrition laboratory, faculty of agriculture, Kyushu University, Japan.
2.2. Analysis of Collected Soil Samples
The collected soil samples were analyzed at Plant Nutrition Laboratory, Faculty of Agriculture, Kyushu University, Fukuoka, Japan. For each collected soil sample, soil pH H2O (1:2.5 soil: H2O) was measured using a pH meter (HM-10P; DKK-TOA Corp., Tokyo, Japan). Soil samples were also digested using the salicylic acid-sulfuric acid-hydrogen peroxide method  ; then, total N was examined using the indophenol method , and total P was tested using the ascorbic acid method . Total K was analyzed using an atomic absorption spectrophotometer (Z-5300; Hitachi, Tokyo, Japan) afterward samples were digested. To analyze mineralizable N, we firstly used the soil incubation method  followed by the indophenol method .
2.3. Isolation of Indigenous Root Nodule Bacteria from Soil Samples of Major Black Gram (Vigna mungo L.) Growing Areas of Myanmar
One gram of each composite soil sample was diluted with 99 ml of sterilized one-half strength modified Hoagland nutrient solution (MHN) in a 200 mL conical flask. The flasks were shaken on a rotary shaker at 120 rpm for one hour to prepare a well-mixed soil suspension. The culture pots (1 L volume) were filled with 1 L of vermiculite and 0.6 L of MHN nutrient solution. The pots were covered with aluminum foil and autoclaved at 121˚C for 20 min. For surface sterilization, the seeds were soaked in a 2.5% sodium hypochlorite solution for 5 min, rinsed five times with 10 ml of 99.5% ethanol and washed five times with sterilized MHN nutrient solution to remove traces of sodium hypochlorite and ethanol. Myanmar Yezin-7 black gram variety, common use in Myanmar was used as trap host for all soil samples.
The culture pots using black gram Yezin-7 variety with eight soil suspensions and control pot were prepared. The seeds were surface-sterilized and planted in the sterilized vermiculite pots. 5 ml aliquot of soil suspension was inoculated per seed. The control was planted without inoculation to assess the possibility of contamination. The plants were cultivated in the incubator (25˚C and 16 hours light) for four weeks. Autoclaved deionized water was poured when the original weight of the pots decreased by around 300 g.
After carefully uprooting, the five nodules (≥2 mm in diameter) were collected per pot. For surface sterilization, the nodules were soaked in 70% ethanol for 3 min, 2.5% sodium hypochlorite (NaClO) solution for 15 min and washed five times with 0.9% autoclaved sodium chloride (NaCl) solution. The surface sterilized nodules were transferred separately into the autoclaved small test tubes and crushed. For every sample, a loopful of the suspension was streaked on Yeast Extract Mannitol Agar (YMA) plates containing 25-µg Congo red . The plates were incubated at 30˚C for 1 - 3 days for fast growing bacteria and 5 - 7 days for slow growing bacteria. Finally, purified 40 indigenous root nodule bacteria isolates were isolated from soil samples of major Black gram (Vigna mungo L.) growing areas of Myanmar. The collected isolates were nominated as HpaBG1 to LauBG40.
2.4. DNA Extraction, PCR Analysis and Phylogenetic Analysis
For DNA extraction, the collected isolated were streaked onto A1E agar plates and incubated at 30˚C for 7 days. A single pure colony of each isolate from A1E plates was cultured in AIE liquid medium at 30˚C for 5 days to obtain the required optimum density (0.4 < OD600 nm < 0.6). Total DNA was extracted using ISOPLANT (Nippon gene, Tokyo, Japan), following instructions from the manufacturer. The DNA concentrations were calculated using NIH Image 1.62 (National Institutes of Health, Bethesda, MD, USA) after agarose gel electrophoresis (0.3% agarose gel in 1 TAE buffer), staining with ethidium bromide (Toyobo, Tokyo, Japan), and destaining in 1 TAE buffer.
The primers 16S-F (5’-AGAGTTTGATCCTGGCTCAG-3’) and 16S-R2 (5’-CGGCTACCTTGTTACGACTT-3’) were used to amplify the 16S rRNA region of mesorhizobia. The PCR reaction consisted of a pre-run at 94˚C for 5 min, denaturation at 94˚C for 30 s, annealing at 60˚C for 30 s, and extension at 72˚C for 1 min . The cycle was repeated for 33 cycles, followed by a final extension at 72˚C for 10 min. PCR products were purified using the Wizard Gel and PCR Clean-up System (Promega, Madison, WI, USA). Purified PCR products (≥50 ng∙µL−1) were subjected to direct sequencing by Macrogen (Tokyo, Japan), using the primer set described above. Raw sequence results were edited using MEGA version 6-software  to create 16S sequence fragments.
For homology searches, sequences were compared with the DNA Data Bank of Japan (DDBJ) using the Basic Local Alignment Search Tool (BLAST) program . To construct the phylogenetic tree, sequences of type strains and closely related strains of Bradyrhizobium genospecies were retrieved from the BLAST database. All selected sequences including type strains and closet strains were aligned using the CLUSTALW function of the MEGA version 6-software . After alignment, a phylogenetic tree was constructed according to the neighbor-joining method . The phylogenetic tree was bootstrapped with 1000 replications of each sequence to evaluate the tree topology for reliability. Genetic distances were calculated using the Kimura two-parameter model .
2.5. Nucleotide Sequence Accession Numbers
The nucleotide sequences of 16S rRNA genes of 40 Bradyrhizobium strains were deposited in the DDBJ under the set of accession numbers LC515811 to LC515850.
2.6. Myanmar Black Gram (Vigna mungo L.) Cultivars and Bradyrhizobium Strains
Myanmar black gram (Vigna mungo L.) cultivars, Yezin-4 and Yezin-7 were collected from Food Legumes Section, Department of Agricultural Research, Yezin, Myanmar. The black gram Yezin-7 variety was the most widely grown cultivar in Myanmar and used for screening of the 40 indigenous Bradyrhizobium strains. For effectiveness of selected Bradyrhizboium strains, Black gram Yezin-4 and Yezin-7 cultivars from Myanmar were investigated. The purified forty Bradyrhizobium strains were cultured in A1E liquid media  on a rotary shaker (100 rpm) at 30˚C for 7 days. The cultures were diluted with sterilized N-free half-strength modified Hoagland nutrient (MHN) solution  to ca 107 cells mL−1 for Bradyrhizobium strains.
2.7. Screening of Effective Bradyrhizobial Strains by Myanmar Black Gram Yezin-4 for Nitrogen Fixation
The purified forty indigenous bradyrhizobial strains were screened for nitrogen fixing effectiveness on Myanmar black gram Yezin-7 in pots with vermiculite and MHN solution were covered with aluminum foil and autoclaved at 121˚C for 20 min. Black gram Yezin-7 seeds were surface sterilized  and germinated on sterile petri dishes with filter paper. The germinated seeds were sown into the sterilized pot with vermiculite. In the non-inoculated treatment, a control treatment was also provided. The weights of original pots were measured. The plants were cultivated in incubator (25˚C and 70% relative humidity). During the growing period, sterilized water was irrigated. The plants in each pot were uprooted and carefully washed with water so as not to detach the nodules. The acetylene reduction assay (ARA) was performed according to Haider et al., 1991  to measure nitrogenase activity. The black gram plants were cut at the cotyledonary nodes. Then, the black gram roots with intact nodules were placed in a 100 mL conical flask and sealed with a serum stopper. A 12 mL aliquot of acetylene (C2H2) gas was injected into the flask to replace the air with acetylene. The flasks containing roots with intact nodules were incubated at room temperature and 1.0 mL subsamples were analyzed at 5 and 65 min, respectively. The ARA value, in terms of ethylene (C2H4) production per plant, was measured using a flame ionization gas chromatograph (GC-14A, Shimadzu, Kyoto, Japan) equipped with a stainless steel column (3 mm diameter, 0.5 m length). The column was filled with Porapak R 60 - 80 mesh (Nicalai Tesque, Inc., Kyoto Japan). Column, injection and detection temperatures were 35˚C, 45˚C and 170˚C, respectively. N2 gas was used as the carrier gas at a flow rate of 45 mL∙min−1. The number of nodules was counted after the assay. Shoots, roots and nodules were collected separately and oven dried for at 70˚C for 48 hours to record the dry weight determination.
2.8. Evaluation the Effectiveness of Selected Bradyrhizboium Strains on Two Myanmar Black Gram Cultivars: Yezin-4 and Yezin-7
The five Bradyrhizobium strains of B. elkanii HpaBG5, B. liaoningense HpaBG6, B. liaoningense HpaBG7, B. liaoningense HpaBG12, B. elkanii LauBG38 (Table 3) were selected based on the results of the above screening experiment in their
Table 3. The morphological characteristics of black gram bradyrhizobial isolate from Myanmar.
Bl1 = Bradyrhizobium liaoningense cluster 1; Bl2 = Bradyrhizobium liaoningense cluster 2; Bo1 = Bradyrhizobium ottawaense cluster 1; Bs1 = Bradyrhizobium sp. cluster 1; Bs2 = Bradyrhizobium sp. cluster 2; By1 = Bradyrhizobium yunamingense cluster 1; Be1 = Bradyrhizobium elkanii cluster 1; UP = Undulated-pulvinate, EP = Entirely-pulvinate.
potential efficiency on ARA per plant. The experiment was performed in completely randomized design with three replicates. The inoculation and growing condition of pot experiment were also conducted as the above experiment. ARA per plant, nodule, root and shoot dry weight were determined after four weeks. Data were analyzed using the STATISTIX 8 software (Analytical Software, Tallahassee, FL, USA), and treatment means were compared by Tukey’s HSD test (P < 0.05) for the collected parameters.
3.1. Diversity of Indigenous Bradyrhizobium Strains for Myanmar Black Gram (Vigna mungo L.) Cultivars
The forty root nodule bacteria were isolated from eight different soil samples from major black gram growing areas in Myanmar (Table 1, Table 2 and Figure 1). According to Somasegaran and Hoben, 1994 , these strains were proved as pure Bradyrhizobium strains. In YMA plates, the bradyrhizobial colonies reached 1 - 3 mm diameter with undulated pulvinate and entirely pulvinate shapes after 5 - 7 days incubation (Table 3).
Neighbor-joining trees for each gene had similar overall tree topologies. Groups were selected on the basis of the minimum standard changes between named species in the 16S rRNA phylogram (Figure 2), and all groups were well supported in neighbor-joining analyses which had less than 50% bootstrap support in the neighbor-joining tree. The results of the phylogenetic analysis based on the 16S rRNA sequence, indicated that all the 40 isolates belonged to the genus Bradyrhizobium (Table 3 and Figure 2).
The seven clusters were identified in the phylogenetic tree including four clusters of B. liaoningense (Bl1 and Bl2), two clusters of B. sp. (Bs1 and Bs2), one cluster of each B. japonicum (Bj1), B. yunamingense (By1) and B. elkanii (Be1) (Figure 2). Among these clusters, two clusters of Bl1 and Bl2 were 98% sequence similarity with B. liaoningense CCNWSX 0360T and B. liaoningense LMG 18230T while two clusters of Bs1 and Bs2 were 98% sequence similar to B. sp. CI 110T and B. sp. JNVU DC11T. Among the clusters belonged to B. japonicum, Bj1 was showing 97% sequence similarity with B. japonicum LMG 6138T, B. japonicum N2 225T B. japonicum CCBAU 83623T and B. japonicum R33T. The last two clusters of By1 and Be1 groups were related to B. yunamingense 65010T and B. elkanii Cte 503T with at least 96% and 98% sequence similarity.
Almost all the collected isolates from Nyaunglebin Bago Region, Chaungzon Mon State and Sittwe Rakhine State black gram growing regions were identified as B. liaoningense. At Danubyu Ayeyarwady Region and Pyinmanar Nay Pyi Taw Region, most of the strains were identified as B. japonicum. On the other hand, more or less all the isolates from Launglon Tanintharyi Region and Hpa-an Kayin State were related to B. elkanii. However, all B. sp. strains were found in Salingyi Sagaing Region black gram growing region (Table 3 and Figure 2).
Figure 2. Position of the 40 strains in the phylogenetic tree based on the 16S rRNA sequences of related Bradyrhizobium strains (in italics) retrieved from GenBank. The tree was constructed by the neighbor-joining method with the Kimura 2-parameter (K2P) distance correlation model and 1000 bootstrap replications. Bootstrap values above 50% are indicated at the nodes. Bar, 0.02 Knuc in nucleotide sequences. Accession numbers of the reference strains, including all type strains of Bradyrhizobium, are shown in parentheses. B: Bradyrhizobium and S: Sinorhizobium. The clustering of isolates and their distribution throughout the total studied area is noted into the tree: Bl1 = Bradyrhizobium liaoningense cluster 1; Bl2 = Bradyrhizobium liaoningense cluster 2; Bj1 = Bradyrhizobium japonicum cluster 1; Bs1 = Bradyrhizobium sp. cluster 1; Bs2 = Bradyrhizobium sp. cluster 2; By1 = Bradyrhizobium yunamingense cluster 1; Be1 = Bradyrhizobium elkanii cluster 1.
In this study, cluster Bl1 and Bl2 were observed in major black gram growing areas Nyaunglebin Bago Regio, Chaungzon Mon State, Sittwe Rakhine State, Danubyu Ayeyarwady Region and Launglon Tanintharyi Region. In Hpa-an Kayin State, Salingyi Sagaing Region and Pyinmanar Nay Pyi Taw Region, cluster Bs1 and Bs1 were found but in Chaungzon Mon State, Danubyu Ayeyarwady Region, Sittwe Rakhine State and Pyinmanar Nay Pyi Taw Region, cluster Bj1 was distributed. However, By1 cluster was only distributed in Salingyi Sagaing Region and Be1 cluster was scattered in both Hpa-an Kayin State and Launglon Tanintharyi Region of black growing areas of Myanmar (Figure 2).
3.2. Screening of Effective Bacterial Strains by Yezin-7 for Nitrogen Fixation
In the screening experiment, the effective strains were determined their potential ability in the N fixation analyzed by means of ARA per plant. Each bacterial strain that responded on black gram Yezin-7 was expressed in Figure 3. The higher ARA per plant was found in plants inoculated by B. elkanii HpaBG5, B. liaoningense ChaBG6, B. liaoningense ChaBG7 and B. elkanii LauBG38 (Figure 3). From this experiment, five indigenous Bradyrhizibium strains were screened out for further studies based on their effectiveness in ARA per plant. There were four Bradyrhizobial strains from the highest ARA per plant value of B. elkanii HpaBG5, B. liaoningense ChaBG6, B. liaoningense ChaBG7, B. elkanii LauBG38 and one strain from the middle ARA per plant value of B. liaoningense NyaBG12.
Figure 3. Acetylene reduction activity per plant of Yezin-7 Black gram (Vigna mungo) variety inoculated with 40 purified Bradyrhizobium strains after four weeks, Bars mean standard deviation levels. The middle line indicates average acetylene reduction activity per plant of all strains; the upper line means the higher acetylene reduction activity per plant and the lower line refers the lower acetylene reduction activity per plant of all strain.
3.3. Effectivity of Selected Bradyrhizobium Strains on Yezin-4 and Yezin-7 of Two Myanmar Black Gram Varieties
Table 4 shows the effectiveness of selected Bradyrhizobial strains on ARA per plant, nodule and shoot dry weight of Yezin-4 and Yezin-7 black gram varieties. The ARA per plant of selected five Bradyrhizobial strains was significant affected on Yezin-4 and Yezin-7 black gram varieties (Table 4). It was found that significant higher in ARA per plant (P < 0.05) was observed by B. elkanii LauBG38 than B. liaoningense NyaBG12 in both Yezin-4 and Yezin-7 black gram varieties. In Yezin-4 black gram variety, the B. elkanii LauBG38 was insignificant differences that compared with B. elkanii HpaBG5, B. liaoningense ChaBG6 and B. liaoningense ChaBG7 in ARA per plant in Table 4.
On the other hand, the B. elkanii LauBG38 showed the highest nodule and shoot dry weights of 9.17 mg plant−1, 0.17 g plant−1 in Yezin-4 and 10.97 mg plant−1, 0.18 g plant−1 in Yezin-7 black gram variety, respectively. In Yezin-4 variety, inoculation with B. elkanii LauBG38 gave significant higher in nodule dry weight than B. liaoningense ChaBG6, B. liaoningense ChaBG7, and B. liaoningense NyaBG12. In Yezin-7 variety, B. elkanii LauBG38 gave the highest nodule dry weight among the tested strains and this strain was significant difference from those given by B. liaoningense ChaBG6 and B. liaoningense NyaBG12 (Table 4). However, shoot dry weight of tested B. elkanii LauBG38 was only significant difference than B. liaoningense NyaBG12 in both Myanmar black gram varieties (Table 4).
The symbiotic association between legumes and rhizobia is one of the most important contributors to the world's supply of biologically fixed nitrogen to agriculture. Effective symbiosis can only be achieved when the nodules are formed
Table 4. Response of selected Myanmar bradyrhizobial strains on acetylene reduction activity, nodule and shoot dry weight of Yezin-4 and Yezin-7 black gram varieties after four weeks.
Means in each column followed by different letters differed significantly at P < 0.05 (Tukey’s test), NDW, SDW means nodule and shoot dry weight per plant and ARA means C2H4 produced per hour per plant, NDW and ARA of non-inoculated treatment is zero for both black gram varieties, SDW of non-inoculated treatment is 0.05 g plant−1 for both back gram varieties.
by effective rhizobia. The symbiotic relationship between rhizobia and black gram has not been extensively analyzed. Therefore, this is the first report investigation of addressing the genetic diversity and evaluation the effectiveness of indigenous Bradyrhizobium strains for Myanmar black gram cultivars.
Sequence analysis of 16S ribosomal RNA (rRNA) has been developed used as one of the most important methods in taxonomy and phylogenic analysis of bacteria    . The phylogenetic investigates based on 16S rRNA gene sequences showed that the rhizobial strains nodulating black gram in India were dispersed under the genus Bradyrhizobium . These results were dependable with former reports that showed black gram rhizobia were more closely associated with Bradyrhizobium species   . The present study was focused on investigation of 16S rRNA region of 40 Bradyrhizobial strains that were successfully isolated from the different soil samples of major black gram growing areas in Myanmar and proved as pure Bradyrhizobium strains .
The genus Bradyrhizobium was proposed by Jordan (1982)  for the slow-growing root nodule bacteria on yeast extract-mannitol agar (YMA) medium. The Bradyrhizobium genus contained 12 defined species: Bradyrhizobium japonicum , B. elkanii , B. liaoningense , B. yuanmingense , B. betae , B. canariense , B. denitrificans , B. iriomotense , B. jicamae, B. pachyrhizi , B. lablabi  and B. cytisi . With the exception of B. betae, all species are symbiotic nitrogen fixing bacteria associated with different legumes.
Previously, Soe et al., 2013  and Htwe et al., 2015  confirmed that the collected rood nodule bacteria isolates from different soybean growing of Myanmar were categorized and identified as the B. japonicum, B. elkanii, B. yuanmingense and Bradyrhizobium sp. and B. liaoningense from Myanmar soybean growing areas. Several scientists found the B. japonicum, B. elkanii, B. liaoningense and B. yuanmingense     have been isolated from root nodules of soybean grown in different regions of China as well. Clear biogeographic patterns have been revealed in the soybean rhizobia and soil types appeared to be the main factor of the biogeography of these bacteria . Cowpea bradyrhizobia isolated in Africa , in China , and in Brazil  were identified as B. elkanii, B. japonicum, B. liaoningense, B. yuanmingense or as novel Bradyrhizobium lineages. Less information is available for indigenous rhizobia nodulating green gram and black gram. Also, unidentified slow growers with different 16S rRNA gene haplotypes were isolated from black gram plants growing in the south of India , while some strains isolated from green gram and black gram plants growing in Thailand were closely related to B. japonicum . In our present investigation, B. liaoningense, B. japonicum, B. sp., B. yunamingense and B. elkanii were collected from major black gram growing areas of Myanmar with different locations, weathers and soil conditions.
The agro-ecological origin of rhizobial inoculants and thus most possibly edaphic and climatic variation are often not considered sufficiently to make inoculation successful. A variety of biotic and abiotic factors, such as host plant, cultivation history, drought, soil pH, salinity, mineral nutrient availability, soil organic carbon content and texture, are known to affect rhizobial diversity and distribution . This is evidence that the collected 40 Bradyrhizobia isolates from major black gram growing areas of Myanmar with various soil types, soil pH and soil nutrients were affected to bradyrhizobial strains diversity and distribution.
A phylogenetic analysis of present study showed that indigenous B. liaoningense as a dominant strain was distributed throughout the five major black gram growing areas of Myanmar with a pH range of 4.69 - 6.70. B. liaoningense strain was isolated from soybean growing areas in China  and in Myanmar .
Symbiotic N2 fixation can recompense for absent soil nitrogen (N) and thus potentially save costly mineral N fertilizer . Rhizobial inocula for inoculating legumes increasingly account for differences in symbiotic specificity and effectivity, two parameters that are often correlated . Several studies also reported significant increase growth parameters and yield due to the inoculation of rhizobial isolates in chickpea in Myanmar , green gram and black gram in Thailand . In addition, indigenous rhizobial strains also play important role since they have adapted to local environmental conditions.
An effective Rhizobium-legume symbiosis largely depends on the presence of a specific and compatible strain in the soil for a particular legume. Several studies have reported a significant increase in green gram and black gram growth parameters and yield due to the inoculation of bradyrhizobial isolates . In the present study, the effectiveness of indigenous 40 Bradyrhizobium strains in symbiosis association with Yezin-7 Myanmar black gram variety by using sterilized vermiculite with MHN as a growth media. When the nitrogen fixation in terms of ARA per plant with the test Bradyrhizobium strains was compared, four Bradyrhizobium strains designated as B. elkanii HpaBG5, B. liaoningense ChaBG6, B. liaoningense ChaBG7 and B. elkanii LauBG38 were found to be more effective than other tested strains. The effectiveness of a strain of rhizobia is due to the genetic interaction with the host plant, which is known as host-strain specificity. Thus, the selection of strains of rhizobia for cultivated legume varieties is a critical step in the production of the legume inoculants .
In this study, the evaluation the effectiveness of selected Bradyrhizboium strains on two Myanmar black gram varieties, Yezin-4 and Yezin-7 were investigated. It was observed that the inoculation of Bradyrhizobium elkanii LauBG38 gave significantly higher in ARA per plant for nitrogen fixation and nodule dry weight in both black gram varieties. ARA per plant for nitrogen fixation and nodule dry weight of Yezin-7 variety was higher than those of Yeizn-4 black gram. Therefore based on the results of nodulation efficiency of nodule dry weight of the plants Yezin-7 could be used in future experiment as superior host genotypes for high nitrogen fixation. Symbiotic nitrogen fixation depends on interactions among the genotype of the host plant, rhizobial strain genotype and environment. In grain legume species, genotypic variability affected nodule number or mass or nitrogenase activity  pointed out that, through the use of plant genotypes in symbiotic ability, it is possible to identify genes responsible for a particular part of the process, depending on a particular rhizobial strain used. The selection of effective rhizobial strains for cultivated legumes is a critical step in the production of high quality legume inoculant.
The effectiveness of indigenous Bradyrhizobium strains was observed in Myanmar back gram cultivars using the correct varieties, and proper nodulated bacteria. This is the first report of phylogenetic diversity and evaluation the effectiveness of indigenous Bradyrhizobium strains for Myanmar black gram cultivars. The selected Bradyrhizobium elkanii LauBG38 strain might be considered for rhizobial inoculatns to use as Biofertilizer in Myanmar near future.
This is the first report on describing Bradyrhizobium strains that were isolated from soil samples of major black gram growing areas of Myanmar and the effectiveness of those strains for plant growth and nitrogen fixation of Myanmar Black gram varieties. In total, 40 indigenous bradyrhizobia were successfully isolated and their geographic distribution was determined based on the analysis of the 16S rRNA region. Our results indicated that B. liaoningense strain was widely distributed in the major black gram growing regions of Myanmar whereas B. japonicum was found more abundant in Danubyu Ayeyarwady Region and Pyinmanar Nay Pyi Taw Region. However, B. elkanii weas found more abundant in Launglon Tanintharyi Region and Hpa-an Kayin State. The effectiveness of those strains for plant growth and nitrogen fixation of Myanmar black gram varieties was investigated in the present study. The 40 Bradyrhizobium strains were screened for their effectiveness on Yezin-7 black gram variety and five Bradyrhizobium strains were selected. These selected strains were tested for their effectiveness on Yezin-4 and Yezin-7 black gram varieties. The Bradyrhizobium elkanii LauBG38 was significantly superior in both black gram varieties. All of these experiments were conducted under the control conditions by growing the plants in the sterilized vermiculite with MHN solution. So, the bradyrhizobial strains selected in the control room trials then must be evaluated in the field. Although it is a preliminary study, it can help for the future study in the inoculants production. The further investigation of symbiotic effectiveness of selected indigenous bradyrhizobia strains on Myanmar black gram cultivars will be examined in the field condition. We do hope that Myanmar Bradyrhizobium strains will be able to use as Biofertilizer for black gram cultivars that enhance crop production through nitrogen fixation and yield.
The authors are thankful to Japan Society for the Promotion of Science (JP19F18079) for their financial support of the present study. We are also very grateful to the members of Land Use Division, Department of Agriculture and members of Department of Agricultural Research, Ministry of Agriculture, Livestock and Irrigation, Myanmar who helped for collecting soil samples and black gram seeds for these experiments.
 Nedumaran, S., Abinaya, P., Jyosthnaa, P., Shraavya, B., Rao, P. and Bantilan, C. (2015) Grain Legumes Production, Consumption and Trade Trends in Developing Countries. Working Paper Series No. 60. ICRISAT Research Program, Telangana, India.
 Jat, B.L., Choudhary. P., Singh, G. and Reddy, G.L. (2017) Effect of Bio-Fertilizer on Different Varieties of Black Gram (Vigna mungo L.). International Journal of Current Microbiology and Applied Sciences, 6, 302-316.
 Somta, P., Chen, J., Yundaeng, C., Yuan, X., Yimram, T., Tomooka N. and Chen, X. (2019) Development of an SNP-Based High-Density Linkage Map and QTL Analysis for Bruchid (Callosobruchus maculatus F.) Resistance in Black Gram (Vigna mungo (L.) Hepper). Scientific Reports, 9, 3930.
 Wai, P.P. (2019) Comparative Study of Profitability of Black Gram Production before and after India’s Import Suspension in Kyauktaga Township, Bago Region (Ph.D. Thesis). Yezin Agriculture University, Myanmar, 118.
 Korir, H., Mungai, N.W., Thuita, M., Hamba Y. and Masso C. (2017) Co-inoculation Effect of Rhizobia and Plant Growth Promoting Rhizobacteria on Common Bean Growth in a Low Phosphorus Soil. Frontiers in Plant Science, 8, 141.
 Sessitsch, A., Howieson, J.G., Perret, X., Antoun, H. and Martínez-Romero, E. (2002) Advances in Rhizobium Research. Critical Reviews in Plant Sciences, 21, 323-378.
 Than, M.M., San, K.K. and Thein, M.M. (2006) Evaluation of Effective Rhizobial Strains for Commercial Legume Inoculants. Journal of Agriculture, Forestry Livestock and Fishery Science, 6, 264-280.
 Than, M.M., Myint, H., Lwin, T. and Myint, Y.Y. (2003) Comparison of Nodulation Efficiency between Native Rhizobium Strains and Exotic TAL Strains on Green Gram (Vigna radiata). In: Proceedings of the 3rd Agricultural Conference, Yezin Agriculture University, Yezin, Myanmar, 50-58.
 Brockwell, J., Dudman, W.F., Gibson, A.H., Hely, F.W. and Erbinsion, A.C. (1968) An Integrated Programme for the Improvement of Legume Inoculant Strains. Transactions 9th International Congress Soil Science, 2, 103-114.
 Herridge, D., Maw, J.B., Thein, M.M., Rupela, O.P., Boonkerd, N., Thao, T.Y., Deaker, R., Hartley, E. and Gemell, G. (2008) Expanding Production and Use of Legume Inoculants in Myanmar and Vietnam.
 Than, M.M. (2010) Evaluation and Selection of Root Nodule Bacteria (Mesorhizobium ciceri) and Chickpea Germplasm for High Nitrogen Fixation. Ph.D. Thesis, Yezin Agriculture University, Myanmar, 118.
 Soe, K.M., Bhromsiri, A. and Karladee, D. (2010) Effect of Selected Endophytic actinomycetes (Steptomyces sp.) and Bradyrhizobia from Myanmar on Growth, Nodulation, Nitrogen Fixation and Yield of Different Soybean Varieties. Chiang Mai University Journal of Natural Sciences, 9, 95-109.
 Soe, K.M. and Yamakawa, T. (2013a) Low-Density Co-Inoculation of Myanmar Bradyrhizobium yuanmingense MAS34 and Streptomyces griseoflavus P4 to Enhance Symbiosis and Seed Yield in Soybean Varieties. American Journal of Plant Sciences, 4, 1879-1892.
 Htwe, A.Z., Moh, S.M., Soe, K.M., Moe, K. and Yamakawa, T. (2019) Effects of Biofertilizer Produced from Bradyrhizobium and Streptomyces griseoflavus on Plant Growth, Nodulation, Nitrogen Fixation, Nutrient Uptake, and Seed Yield of Mung Bean, Cowpea, and Soybean. Agronomy, 9, 77.
 Aung, L.L., Zin, E.E., Theingi, P., Elvera, N., Aung, P.P., Han, T.T., Oo, Y. and Skaland, R.G. (2017) Myanmar Climate Report. Department of Meteorology and Hydrology Myanmar, Ministry of Transport and Communications, Government of the Republic of the Union of Myanmar.
 Ohyama, T., Ito, M., Kobayashi, K., Araki, S., Yasuyoshi, S., Sasaki, O., Yamazaki, T., Soyama, K., Tanemura, R. and Mizuno, Y. (1991) Analytical Procedures of N, P, K Contents in Plant and Manure Materials Using H2SO4-H2O2 Kjeldahl Digestion method. Bulletin Faculty of Agriculture, Niigata University, Japan, 43, 110-120.
 Cataldo, D.A., Schrader, L.E. and Youngs, V.L. (1974) Analysis by Digestion and Colorimetric Assay of Total Nitrogen in Plant Tissues High in Nitrate. Crop Science, 14, 854-856.
 Murphy, J. and Riley, J. (1962) A Modified Single Solution Method for the Determination of Phosphate in Natural Waters. Analytica Chimica Acta, 27, 31-36.
 Muramoto, J., Goto, I. and Ninaki, M. (1992) Rapid Analysis of Exchangeable Cations and Cation Exchange Capacity (CEC) of Soils by Shaking Extraction Method. Journal of Soil Science and Plant Nutrition, 63, 210-215.
 Sarr, P.S., Yamakawa, T., Fujimoto, S., Saeki, Y., Thao, H.T.B. and Myint, A.K. (2009) Phylogenetic Diversity and Symbiotic Effectiveness of Root-Nodulating Bacteria Associated with Cowpea in the South-West Area of Japan. Microbes and Environments, 24, 285-292.
 Tamura, K., Stecher, G., Peterson, D., Filipski, A. and Kumar, S. (2013) MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and Evolution, 30, 2725-2729.
 Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J., Zhang, Z., Miller, W. and Lipman, D.J. (1997) Gapped BLAST and PSI-BLAST: A New Generation of Protein Database Search Programs. Nucleic Acids Research, 25, 3389-3402.
 Kimura, M.A. (1980) Simple Method for Estimating Evolutionary Rates of Base Substitutions through Comparative Studies of Nucleotide Sequences. Journal of Molecular Evolution, 16, 111-120.
 Kuykendall, L.D. (1979) Transfer of R Factors to and between Genetically Marked Sublines of Rhizobium japonicum. Applied and Environmental Microbiology, 37, 862-866.
 Tsurumaru, H., Yamakawa, T., Tanaka, M. and Sakai, M. (2008) Tn5 Mutants of Bradyrhizobium japonicum Is-1 with Altered Compatibility with Rj2-Soybean Cultivars. The Journal of Soil Science and Plant Nutrition, 54, 197-203.
 Yamakawa, T., Shirai, T. and Ishizuka, J. (2000) Effects of Symbiosis with Rhizobium fredii on Transport of Fixed Nitrogen in the Xylem of Soybean Plant. The Journal of Soil Science and Plant Nutrition, 46, 885-892.
 Haider, J., Hussam, A.K.M.A., Ikeda, M., Yamakawa, T. and Ishizuka, J. (1991) Effects of Nitrate Application on Growth, Nodulation and Nitrogen Fixation of Nitrate-Tolerant Mutant Soybean. The Journal of Soil Science and Plant Nutrition, 37, 521-529.
 Weisburg, W.G., Barns, S.M., Pelletier, D.A. and Lane, D.J. (1991) 16S Ribosomal DNA Amplification for Phylogenetic Study. Journal of Bacteriology, 173, 697-703.
 Chen, W.M., Zhu, W.F., Bontemps, C., Young, J.P.W. and Wei, G.H. (2010) Mesorhizobium alhagi sp. nov., Isolated from Wild Alhagi sparsifolia in North-Western China. International Journal of Systematic and Evolutionary Microbiology, 60, 958-962.
 Martens, M., Delaere, M., Coopman, R., De Vos, P., Gillis, M. and Willems, A. (2007) Multilocus Sequence Analysis of Ensifer and Related Taxa. International Journal of Systematic and Evolutionary Microbiology, 57, 489-503.
 Laranjo, M., Young, J.P.W. and Oliveira, S. (2012) Multilocus Sequence Analysis Reveals Multiple Symbiovars within Mesorhizobium Species. Systematic and Applied Microbiology, 35, 359-367.
 Saleena L.M., Loganathan, P., Rangarajan, S. and Nair, S. (2001) Genetic Diversity and Relationship between Bradyrhizobium Strains Isolated from Blackgram and Cowpea. Biol Fertil Soils, 34, 276-281.
 Javaid, A., Bajwa, R. and Anjum, T. (2006) Response of Black Gram (Vigna mungo (L.) Hepper) to Bradyrhizobium japonicum Inoculation under Different Soil Amendment Systems. Pakistan Journal of Botany, 38, 851-857.
 Appunua C., N’Zoued, A., Moulind, L., Depret, G. and Laguerre, G. (2009) Vigna mungo, V. radiata and V. unguiculata Plants Sampled in Different Agronomical-Ecological-Climatic Regions of India Are Nodulated by Bradyrhizobium yuanmingense. Systematic and Applied Microbiology, 32, 460-470.
 Jordan, D. C. (1982) Transfer of Rhizobium japonicum to Bradyrhizobium gen. nov., a Genus of Slow-Growing, Root Nodule Bacteria from Leguminous Plants. International Journal of Systematic and Evolutionary Microbiology, 32, 136-139.
 Kuykendall, L.D., Saxena, B., Devine, T.E. and Udell, S.E. (1992) Genetic Diversity in Bradyrhizobium japonicum Jordan 1982 and a Proposal for Bradyrhizobium elkanii sp. nov. The Canadian Journal of Microbiology, 38, 501-505.
 Xu, L.M., Ge, C., Cui, Z., Li, J. and Fan, H. (1995) Bradyrhizobium liaoningense sp. nov., Isolated from the Root Nodules of Soybeans. International Journal of Systematic and Evolutionary Microbiology, 45, 706-711.
 Yao, Z.Y., Kan, F.L., Wang, E.T., Wei, G.H. and Chen, W.X. (2002) Characterization of Rhizobia that Nodulate Legume Species of the Genus Lespedeza and Description of Bradyrhizobium yuanmingense sp. nov. International Journal of Systematic and Evolutionary Microbiology, 52, 2219-2230.
 Rivas, R., Willems, A., Palomo, J.L., Garcia-Benavides, P., Mateos, P.F., Martinez-Molina, E., Gillis, M. and Velazquez, E. (2004) Bradyrhizobium betae sp. nov., Isolated from Roots of Beta vulgaris Affected by Tumour-Like Deformations. International Journal of Systematic and Evolutionary Microbiology, 54, 1271-1275.
 Vinuesa, P., Leon-Barrios, M., Silva, C., Willems, A., Jarabo-Lorenzo, A., Perez-Galdona, R., Werner, D. and Martinez-Romero, E. (2005) Bradyrhizobium canariense sp. nov., an Acid-Tolerant Endosymbiont that Nodulates Genistoid Legumes (Papilionoideae: Genisteae) from the Canary Islands, along with Bradyrhizobium japonicum bv. Genistearum, Bradyrhizobium genospecies Alpha and Bradyrhizobium genospecies Beta. International Journal of Systematic and Evolutionary Microbiology, 55, 569-575.
 Van Berkum, P., Leibold, J.M. and Eardly, B.D. (2006) Proposal for Combining Bradyrhizobium spp. (Aeschynomene indica) with Blastobacter denitrificans and to Transfer Blastobacter denitrificans (Hirsch and Muller, 1985) to the Genus Bradyrhizobium as Bradyrhizobium denitrificans (comb. nov.). Systematic and Applied Microbiology, 29, 207-215.
 Islam, M.S., Kawasaki, H., Muramatsu, Y., Nakagawa, Y. and Seki, T. (2008) Bradyrhizobium iriomotense sp. nov., Isolated from a Tumorlike Root of the Legume Entada koshunensis from Iriomote Island in Japan. Bioscience, Biotechnology, and Biochemistry, 72, 1416-1429.
 Ramirez-Bahena, M.H., Peix, A., Rivas, R., Camacho, M., Rodriguez-Navarro, D.N., Mateos, P.F., Martinez-Molina, E., Willems, A. and Velazquez, E. (2009) Bradyrhizobium pachyrhizi sp. nov. and Bradyrhizobium jicamae sp. nov., Isolated from Effective Nodules of Pachyrhizus erosus. International Journal of Systematic and Evolutionary Microbiology, 59, 1929-1934.
 Chang, Y.L., Wang, J.Y., Wang, E.T., Liu, H.C., Sui, X.H. and Chen, W.X. (2011) Bradyrhizobium lablabi sp. nov., Isolated from Effective Nodules of Lablab purpureus and Arachis hypogaea. International Journal of Systematic and Evolutionary Microbiology, 61, 2496-2502.
 Soe, K.M., Yamakawa, T., Hashimotoa, S. and Sarr, P.S. (2013) Phylogenetic Diversity of Indigenous Soya Bean Bradyrhizobia from Different Agro-Climatic Regions in Myanmar. Science Asia, 39, 574-583.
 Htwe, A.Z., Yamakawa, T., Sarr. P.S. and Sakata, T. (2015) Diversity and Distribution of Soybean-Nodulating Bradyrhizobia Isolated from Major Soybean-Growing Regions in Myanmar. African Journal of Microbiology Research, 9, 2183-2196.
 Han, L., Wang, E., Han, T., Liu, J., Sui, X., Chen, W. and Chen, W. (2009) Unique Community Structure and Biogeography of Soybean Rhizobia in the Saline-Alkaline Soils of Xinjiang, China. Plant & Soil, 324, 291-305.
 Li, Q.Q., Wang, E.T., Zhang, Y.Z., Zhang, Y.M., Tian, C.F., Sui, X.H., Chen, W.F. and Chen, W.X. (2011) Diversity and Biogeography of Rhizobia Isolated from Root Nodules of Glycine max Grown in Hebei Province, China. Microbial Ecology, 61, 917-931.
 Man, C., Wang, H., Chen, W., Sui, X., Wang, E. and Chen, W. (2008) Diverse Rhizobia Associated with Soybean Grown in the Subtropical and Tropical Regions of China. Plant & Soil, 310, 77-87.
 Wang, H., Man, C., Wang, E. and Chen, W. (2009) Diversity of Rhizobia and Interactions among the Host Legumes and Rhizobial Genotypes in an Agricultural-Forestry Ecosystem. Plant & Soil, 314, 169-182.
 Zhang, Y.M., Li Jr., Y., Chen, W.F., Wang, E.T., Tian, C.F., Li, Q.Q., Zhang, Y.Z., Sui, X.H. and Chen, W.X. (2011) Biodiversity and Biogeography of Rhizobia Associated with Soybean Plants Grown in the North China Plain. Applied and Environmental Microbiology, 77, 6331-6342.
 Steenkamp, E.T., Stepkowski, T., Przymusiak, A., Botha, W.J. and Law, I.J. (2008) Cowpea and Peanut in Southern Africa Are Nodulated by Diverse Bradyrhizobium Strains Harbouring Nodulation Genes that Belong to the Large Pantropical Clade Common in Africa. Molecular Phylogenetics and Evolution, 48, 1131-1144.
 Zhang, Y.F., Wang, E.T., Tian, C.F., Wang, F.Q., Han, L.L., Chen, W.F. and Chen, W.X. (2008) Bradyrhizobium elkanii, Bradyrhizobium yuanmingense and Bradyrhizobium japonicum Are the Main Rhizobia Associated with Vigna unguiculata and Vigna radiata in the Subtropical Region of China. FEMS Microbiology Letters, 285, 146-154.
 Zilli, J.E., Valisheski, R.R., Filho, F.R.F., Neves, M.C.P. and Rumjanek, N.G. (2004) Assessment of Cowpea Rhizobium Diversity in Cerrado Areas of Northeast Brazil. Brazilian Journal of Microbiology, 35, 281-287.
 Saleena, L.M., Loganathan, P., Sunita, R. and Nair, S. (2001) Genetic Diversity and Relationship between Bradyrhizobium Strains Isolated from Black Gram and Cowpea. Biology and Fertility of Soils, 34, 276-281.
 Yokoyama, T.N., Tomooka, N., Okabayashi, M., Kaga, A., Boonkerd, N. and Vaughan, D.A. (2006) Variation in the Nod Gene RFLPs, Nucleotide Sequences of 16S rRNA Genes, Nod Factors, and Nodulation Abilities of Bradyrhizobium Strains Isolated from Thai Vigna Plants. The Canadian Journal of Microbiology, 52, 31-46.
 Wade, T.K., Le Quere, A., Laguerre, G., N’Zoue, A., Ndione, J.-A., Dorego, F., Sadio, O., Ndoye, I. and Neyra, M. (2014) Eco-Geographical Diversity of Cowpea Bradyrhizobia in Senegal Is Marked by Dominance of Two Genetic Types. Systematic and Applied Microbiology, 37, 129-139.
 Guimaraes, A.A., Duque Jaramillo, P.M., Abrahao Nobrega, R.S., Florentino, L.A., Silva, K.B. and de Souza Moreira, F.M. (2012) Genetic and Symbiotic Diversity of Nitrogen Fixing Bacteria Isolated from Agricultural Soils in the Western Amazon by Using Cowpea as the Trap Plant. Applied and Environmental Microbiology, 78, 6726-6733.
 Batista, L., Irisarri, P., Rebuffo, M., Jose Cuitino, M., Sanjuan, J. and Monza, J. (2015) Nodulation Competitiveness as a Requisite for Improved Rhizobial Inoculants of Trifolium pratense. Biology and Fertility of Soils, 51, 11-20.
 Payakapong, W., Tittabutr, P., Teaumroong N. and Boonkerd, N. (2004) Soybean Cultivars Affect Nodulation Competition of Bradyrhizobium japonicum Strains. World Journal of Microbiology and Biotechnology, 20, 311-315.
 Sivasakthi, S., Saranraj, P. and Sivasakthivelan, P. (2019) Beneficial Effect of Bradyrhizobium sp. on Black Gram (Vigna mungo L.) Cultivation: A Review. Life Science Archives (LSA), 3, 1195-1206.
 Wani, S.P., Rupela, O.P. and Lee, K.K. (1995) Sustainable Agriculture in a Semiarid Tropic through Biological Nitrogen Fixation in Grain Legume. Plant & Soil, 174, 29-49.