Abstract

Chickpeas (Cicer arietinum L.) are native to the Middle East (ME) and must be inoculated with symbiotic bacteria (Mesorhizobium ciceri) to fix nitrogen (N) in North American soils. Whether commercial M. ciceri strains are more or less effective than wild strains from ME soils when paired with various chickpea hosts must be elucidated. Wild N-fixing bacterial strains were isolated from ME soils, and their effectiveness was compared against commercial strains on US and ME chickpea varieties. Chickpeas were inoculated with individual strains and grown in chambers for 8 weeks. Plants received 2 mM (15NH4)2 SO4 (5% atom excess) to measure N fixation by isotope dilution. Plant below- and above-ground biomass and proportion of N fixed (PNF) were determined. Commercial and wild ME strains were examined for genetic diversity by sequencing their 16 S rDNA region. The PNF was significantly influenced by inoculant strains and chickpea varieties. Among varieties, Sierra, Troy, and Almaz had the highest PNF of 86.7%, 85.3%, and 85.2%, respectively. Among strains, Jord-M1 contributed to greater PNF (84.7%) compared to Syr-M1 (81.4%). Overall, chickpea varieties had greater effect on PNF than strain selection. These findings support efforts focusing on varietal breeding and strain selection to increase agricultural N fixation.

1. Introduction

Global demand for nitrogen fertilizer is predicted to increase 1.4% annually, which includes a 23.1 million tonnes increase by 2012 [1]. In the US, N fertilizer costs approximately tripled from 2002–2011 [2]. Although prices have decreased somewhat from the 2008 peak, the loss of billions of dollars in farm profit has drawn attention to the need for alternative sources of N. The primary feedstock for synthetic nitrogen fertilizer is natural gas, both as a source of hydrogen and energy [3]. The US is estimated to have only 3.7% of the world natural gas reserves [4], which suggests that American agriculture would benefit from less expensive and more sustainable N sources.

Chickpea (Cicer arietinum L.), a drought-resistant legume from the Middle East (ME), has been cultivated for over 10,000 years [5]. Chickpeas belong to two major types, desi and kabuli. Desi varieties have a small dark seed with a rough coat while kabuli seeds are larger with lighter color and smoother coat [6]. In the US, 51,435 ha were planted with chickpeas in 2011 with a total production of 103,179 tonnes [7]. Chickpeas form root nodules that support biological N fixation (BNF) and host symbiotic N-fixing bacteria. On average the estimated amounts of N fixed by chickpeas under regular precipitation and drought stress conditions are 60 kg/ha [8] and 19–24 kg/ha [9], respectively.

There is increasing evidence that suggests that more N can be fixed by existing legume grain crops if they are inoculated more often or with more effective strains of rhizobia [10]. Significant increases in BNF can result from the use of different commercial or wild strains [11]. The efficiency of biological N fixation is also influenced by cultivars’ selection; significant differences were observed among pea cultivars for the percentage of plant N supplied by bacterial N fixation and also for the number of root nodules formed per plant in a greenhouse experiment [12]. The known N-fixing symbiont of chickpea, M. ciceri, is not native to North American soils, and the number of commercial strains is limited. In addition, there has been little work to enhance the ability of chickpeas to support N fixation [13]. Plant breeders may be able to improve chickpea varieties through the use of breeding lines that are better N-fixing symbiont hosts and could increase the BNF process.

The objectives of this study were to isolate and identify new N-fixing symbiont strains for chickpea from ME soil and to compare the BNF effectiveness of the new strains and current commercial strains when paired with US and ME varieties of chickpeas.

2. Materials and Methods

2.1. N Fixation Study

Soil was collected at 10–15 cm depth from non-inoculated chickpea fields in the following Middle Eastern cities with the help of the following collaborating entities: (1) Terbol, Lebanon (American University of Beirut (AUB)), (2) Tal Hadya, Syria (International Center for AgriculturaI Research in the Dry Areas (ICARDA)), (3) Erbid, Jordan (National Center for Agricultural Research and Technology Transfer of the Ministry of Agriculture (NCARTT)), and (4) Yozgat, Turkey (International Maize and Wheat Improvement Center (CYMMIT)). To isolate symbiotic BNF bacteria from these soils, chickpea seeds (Arman and Sierra) were sterilized with 0.5% sodium hypochlorite for three minutes and germinated for 5 days on sterile filter paper. Germinated seeds were planted in 164 mL seedling cones (SC10 container, Stuewe & Sons, Inc.) filled with sterile sand and perlite (1 : 1 volume) and inoculated with 10 mL of 10−1 soil solution. The pots (4 cm × 20 cm) were placed for 8 wks in growth chambers at 20°C and 16 hrs light. At harvest, large red nodules were selected, surface sterilized for 1 min with 0.5% sodium hypochlorite, and a loop of inoculum was streaked on Yeast Extract Mannitol Agar (YEMA). Two nodules were collected from two different plants per soil, to produce a total of eight wild strains isolated from single colonies (Leb-M1, Leb-M2, Syr-M1, Syr-M2, Jord-M1, Jord-M2, Turk-M1, Turk-M2). Commercial strains of mesorhizobia were also provided from Novozyme Biologicals (NB-M4, NB-M5, NB-M6, NB-M7, NB-M8), Becker Underwood (BU-M3), and EMD (EMD-M1). Three varieties of US chickpeas (Myles, Troy, and Sierra) were supplied by the USDA-ARS Grain Legume Genetics-Physiology Research Unit, Pullman, WA, and three ME varieties (Arman, Almaz, and Hashem) were supplied from the U.S. National Plant Germplasm System, Regional Plant Introduction Station in Pullman, WA. Five days before seed inoculation, strains were grown in 1 X yeast extract broth at 22°C on an orbital shaker. Cells were collected by centrifugation at 4000 rpm for 4 min. An optical density ( A = 0 . 3 ) corresponding to a bacterial count of 106–107 was determined using a spectrophotometer at 640 nm wavelength. Inocula were diluted to 105–106 cfu per mL. One mL of wild (Leb-M1, Syr-M1, Jord-M1, Turk-M1) or commercial diluted inocula (NB-M4, NB-M5, BU-M3, EMD-M1) was added to each cone, which contained one chickpea seed. Arman was also inoculated with an additional four wild (Leb-M2, Syr-M2, Jord-M2, Turk-M2) and four commercial strains (NB-M6, NB-M7, NB-M8, EMD-M6).

Treatments were a factorial design of two main factors, rhizobial strain and crop variety. Varieties by strain pairings were grown in 164 mL seedling cones filled with sand-perlite mixture (1 : 1 volume) using five replicates. For precise measurement of biological N fixation, an isotopic 15-N method was used. Two conditions must exist for the effective use of the isotope dilution method for determining BNF: (1) the addition of N should be low and (2) the N addition should be thoroughly mixed in the medium [1416].

Plants were irrigated with 10 mL 2 mM (15NH4)2 SO4 (5% atom excess) three times a week for the first two weeks, and with H2O thereafter [15]. Sterile N-free Hoagland solution was added once a week to provide all necessary nutrients except N. After 8 weeks, the whole plants were harvested and dried at 65°C, after which above-ground biomass (AGB) and below-ground biomass (BGB) were separated at the root crown, and weighed. Atom percent 15N/14N (A%) of combined AGB and BGB was measured by an isotope ratio mass spectrometer (Thermo-Finigan, Delta Plus Advantage). The atom % excess (A% − 0.366 A%) was used to calculate the proportion of N fixed (PNF): PNF = 1 − [% atom excess 15N in fixing crop/% atom excess 15N in nonfixing crop] × 100 [18]. Wheat was used as a control in the experiment as described by Knowles and Blackburn [19] and Smith et al. [20].

Five replicates were distributed in a randomized complete block design blocked by time. Analysis of variance (ANOVA) and correlations between various parameters were determined using SAS system for Windows version 9.1. Means were calculated with LSM in Proc Mixed. Data were considered significantly different at 𝑃 0 . 0 5 according to Tukey’s method in Proc Mixed.

2.2. PCR, Sequencing and Sequence Analysis

Genomic DNA was extracted with a Fast DNA SPIN KIT (MP Biomedicals Inc, OH) from seven commercial (EMD-M1, BU-M3, NB-M4, NB-M5, NB-M6, NB-M7, and NB-M8) and four wild (Leb-M1, Turk-M2, Syr-M2, and Jord-M1) isolates. In addition, DNA was extracted from two commercial strains of Rhizobium leguminoserum bv. viciae (NB-R2, NB-R7). Using the primer pair fD1 (5′-AGAGTTTGATCCTGGCTCAG-3′) and rD1 (5′-AAGGAGGTGATCCAGCC-3′), an approximately 1340 bp product specific to nearly full length of 16S rRNA gene was amplified by PCR [21]. Cycling conditions were as follows: one cycle of 95°C for 2 min followed by 30 cycles of 94°C for 40 s, 52°C for 40 s, and 72°C for 1 min. This was followed by a final extension at 72°C for 7 min. Amplicons were resolved on a 1.5% agarose (1 X TBE, 90 mM Tris pH 8.0, 90 mM boric acid, 2 mM EDTA) gel, stained with ethidium bromide, and visualized with UV light. The amplicons were purified using a QIAquick PCR Purification kit (QIAGEN Inc, CA) following manufacturer’s instructions, and the purified DNA samples were sequenced in both orientations using the ABI PRISM 377 DNA Sequencer (Applied Biosystems Inc, CA). Sequences were edited, and a consensus sequence for each amplicon was derived using Vector NTI Advance 11 (Invitrogen Inc, CA). Consensus sequences were compared with GenBank accessions by a BLAST search [22]. Multiple sequence alignments were performed by ClustalW [23]. Phylogenetic analyses by the neighbor-joining method and a consensus tree were generated using Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0 [24].

3. Results

3.1. Fixed N

Varieties and strains both had a major influence on the PNF. Averaged across all strains, Sierra had greater PNF than Arman, Hashem, and Myles (Table 1). Myles had a lower PNF than all other varieties. The selection of Sierra could increase PNF in the plant 13.2% over Myles. Varieties from the US and ME were statistically similar in PNF.

Averaged across all varieties, strain Jord-M1 had greater PNF than Syr-M1 (Table 2), but the difference was small. The Jord-M1 strain provided an average 84.7% PNF compared to 81.4% using Syr-M1. Strains from the US and ME provided statistically similar levels of PNF.

3.2. Above- and Below-Ground Biomass

Varieties substantially influenced above and below-ground biomass ( 𝑃 < 0 . 0 0 0 1 ) (Table 1). Sierra and Troy had higher AGB than Arman, Hashem, and Myles. The AGB in Almaz was statistically similar to Sierra, Troy, and Arman and greater than Hashem and Myles (Table 1). This could be due to the larger seed size of these varieties. The average weight per seed in Sierra, Troy, and Almaz was, respectively 0.57 g, 0.60 g, and 0.50 g while Myles weighed 0.16 g (Table 1). Almaz, Arman, and Troy had greater BGB than Myles, Hashem, and Sierra. The BGB in Troy was statistically similar to Sierra, Almaz, and Arman but greater than Hashem, and Myles (Table 1). The ME varieties had significantly greater BGB than the US varieties, which could be directly or indirectly related to the breeding traits selected in the US and ME. A greater BGB may provide the ME varieties a greater or wider tolerance to adverse environmental factors such as water stress. The PNF was highly correlated to AGB ( 𝑟 = 0 . 6 9 ; 𝑃 < 0 . 0 0 0 1 ) and BGB ( 𝑟 = 0 . 3 2 ; 𝑃 < 0 . 0 0 0 1 ), indicating that generally the larger plants had a lower proportion of plant N supplied by the fertilizer.

Arman tests with four additional ME and four additional commercial strains did not reveal a significant difference among these strains on the measured parameters (data not shown). All commercial and ME strains used in this trial were infective and effective on Arman.

3.3. Sequence Analysis

Wild ME nodulators and commercial inoculant strains separated into distinct phylogenetic clades, which were also distinct from the R. leguminoserum outgroup (Figure 1). One commercial strain, NB-M5, was also distinct from all other isolates. A BLAST search of 16S rDNA sequences identified all commercial strains as M. ciceri and all wild ME symbionts as Burkholderia sp. Syr-M2 had 99% sequence identity with Burkholderia fungorum, Jord-M1 had 99% sequence identity with Burkholderia cepacia, and Turk-M2 and Leb-M1 were identified only as Burkholderia sp. Wild ME symbiont sequences were 97-98% similar among them, but only 79-80% similar to commercial strains of M. ciceri and 80-81% similar to R. leguminoserum (data not shown). This is the first report of Burkholderia that nodulates chickpea.

4. Discussion

In this study we used a nonlegume crop as a control and isotope dilution methods to track the proportion of N fixed [19] via Fried and Middleboe’s equation [18]. The total N considered “fixed” in Fried and Middleboe’s equation actually originates from all nonfertilizer sources. There was a positive correlation ( 𝑟 = 0 . 8 0 ) between the weight of the seeds and PNF in the harvested plant. An increase in the seed weight corresponded to an increase in PNF and suggests that 14N was supplied from the seeds. When seed N was subtracted from PNF as defined in Fried and Middleboe equation, only two varieties demonstrated true N fixation. Myles had 8.8% true PNF and Hashem had 14.5% true PNF. In this case, Fried and Middleboe’s equation for the determination of PNF that the plants derived from the atmosphere needs to be adjusted by subtracting the amount of N in the plant derived from the seeds.

The ability of Burkholderia sp. to symbiotically fix N on other legumes has been recognized recently. Bontemps et al. [25] tested 143 symbionts from 47 Brazilian wild species of Mimosa and found that 98% were Burkholderia sp. Two strains of Burkholderia caribensis had been isolated previously from Mimosa nodules [26]. Burkholderia tuberum isolated from a South African legume has been proved to contain nod and nif genes similar to Rhizobiaceae and showed a genetic response with bean and cowpea [27]. Burkholderia cepacia, which was isolated from Acacia decurrens and Ormosia nitida in Brazil, contained nod C and nif H genes [28]. Several studies have demonstrated Burkholderia sp. isolated from nodules of papilionoid and other mimosoid legumes are able to form nodules on selective variety of legumes [29, 30]. Burkholderia belongs to the betaproteobacteria class and carries similar nod genes to rhizobia, which are alphaproteobacteria. The transfer of nod genes to the new N-fixing bacteria may have occurred via lateral gene transfer [29, 31].

Other research has demonstrated a strong influence of plant variety on N fixation. Bello et al. [32] found 70% and 25% variability in kg N fixed ha−1 among three soybean varieties at two sites. The survival and number of inocululated rhizobia are very critical for successful nodulation and yet limited efforts have been made to improve rhizobial inoculants [11]. Inoculated rhizobia strains may not compete well against resident rhizobia. Inoculation by commercial strains filled 5–20% of nodules on soybeans in studies by Caldwell and Vest [33] and Kuykendall and Weber [34]. In these trials, variety had a stronger influence on PNF than did strain, and all strains tested were effective on all varieties. Consequently, the more efficient strategy for BNF improvement may be in breeding crops that support greater N fixation rather than rhizobia strain selection.

5. Conclusion

Increased BNF by legume-bacterial associations can be used to further improve global crop production or reduce the need for synthetic N fertilizers. As the cost of fertilizers increases, the value of improving BNF also increases. Variety of chickpeas had a greater effect on variation in PNF than did the symbionts strain. This study indicates the possibility to improve the efficiency of N fixation through selection of chickpea varieties that host effective symbionts and have a high demand for N. Four newly isolated N-fixing chickpea symbionts from Middle Eastern soils were identified as Burkholderia sp. Studies are needed to investigate additional N-fixing bacteria beyond those found within the Rhizobiaceae. A new strain from Jordan, Jord-M1 contributed to the highest PNF. These findings suggest that variety selection may potentially improve PNF with a wide variety of symbionts, leading to higher yield, higher protein content, and/or more residual N for following crops.

Abbreviations

PNF:Proportion of plant nitrogen provided by fixation
BNF:Biological nitrogen fixation
N:Nitrogen
ME:Middle East
US:United States
AGB:Above-ground biomass
BGB:Below-ground biomass.