Abundance of culturable endophytic bacteria in wheat seeds
The mean population densities of culturable, bacterial endophytes ranged from 2.79+ 0.05 to 5.19 + 0.06 log CFU g-1 d.w. seed, with variations which could be attributed to the type of medium tested and the wheat cultivation zones. The total endophytic bacterial population densities were the least in the CZ with the use of nutrient agar (NA) medium while the highest was in the NHZ with the use of newly designed, Wheat Flour-based Agar (WFA) medium (SaiPrasad, et al. 2020; Table 2). The zonal influences were apparently lesser on the total population densities, when cultured using the Trypticase Soy Agar (TSA) medium. In general, the culturable endophytic bacterial populations were lesser in the seeds sampled from the CZ, PZ, and SHZ than those from the NEPZ, NHZ, and NWPZ. The culturability of bacterial endophytes differed among the media used and additionally, the abundance of these bacteria was typically more in the seeds from the NEPZ, followed by the NWPZ and NHZ. On the basis of morphometric analyses including morphology, size, and pigmentation of colonies, a total of 220 endophytes (43, 40, 38, 28, 34 and 37 from the NWPZ, NEPZ, NHZ, CZ, PZ and SHZ, respectively) were selected for biochemical and phylogenetic investigations.
Phenotyping and Phylotyping of Wheat Seed Endophytic Bacteria
The phenotyping of Wheat Seed Endophytic Bacteria (WSEB) which was performed by examining their growth at different temperatures (4-40 oC), pH (3 to 10), salt concentrations (5-15%), and polyethylene glycol concentrations (PEG 6000, 10-15%) showed considerable variations in their potentials (Supplementary Fig. 1).
The most predominant, distinctly different isolates from each zone were further subjected to the phylotyping by sequencing the partial 16S rRNA genes. The accession numbers of the NCBI GenBank for the phylotypes (46) were MT184815–MT184857(Supplementary Fig. 2). The diversity indices of these phylotypes in the wheat cultivation zones showed marginal differences (Table 2). The Shannon diversity (H) value was the highest in the PZ while those of Chao1 and Simpson’s reciprocal indices were the lowest in the NEPZ. The values of species evenness were higher in the NHZ, followed by the NEPZ. Irrespective of these minor differences in the diversity indices, seeds were enriched with three bacterial phyla (i.e., Firmicutes, Actinobacteria and Proteobacteria with the distribution ratios of 87.0%, 6.5% and 6.5%, respectively) (Supplementary Fig. 3). In all the six wheat cultivation zones, the members of Bacillus belonging to the Phylum Firmicutes were predominant and the numbers of other identified species were about three each in Actinobacteria and Proteobacteria.The niche specific species were Bacillus cereus in the NWPZ, Saccharibacillussacchari in the NEPZ, B. paranthracis in the NHZ, B. australimaris in the PZ and B. aerius in the SHZ, respectively (Fig. 1). The phylogenetic relationships among the identified bacterial species showed that the endophytes were more diverse in the NWPZ than in other zones.
Table 1
Major wheat cultivation zones with the characteristic features of area, productivity, number of varieties released and constraints
Wheat cultivation zone
|
Area (Million ha)
|
Productivity (t/ha)
|
No. of varieties notified
|
Major varieties
|
Major zonal constraints
|
Number of samples
|
|
Norther Hills Zone - NHZ
(31° 06' 12'' N: 77° 10' 20'' E)
|
0.8
|
16.6
|
36
|
VL 616, HS 240, HPW 251, HS 507, VL 907, VL 892, HPW 349
|
Low temperature, Moisture stress, Low soil fertility, frost damage
|
15
|
North Western Plains Zone - NWPZ
(28° 36' 50'' N: 77° 12' 32'' E)
|
11.6
|
39.4
|
80
|
HD 2009, HD 2285, HD 2329, PBW 343, DBW 17, HD 2967, HD 3086
|
High fertilizer dose, Decline of water table, intensive tillage, terminal heat stress
|
18
|
North Eastern Plains Zone – NEPZ
(25° 51' 39'' N: 85° 46' 56'' E)
|
10.5
|
25.1
|
55
|
UP 262, HUW 234, HD 2967, DBW 14, DBW 39, LBW 38, K 0307, K 1006, HD 2733, NW 5054
|
High temperature, Moisture stress, Micronutrient deficiencies
|
12
|
Central Zone – CZ
(22° 43' 31'' N: 75° 51' 55''E)
|
5.2
|
24.1
|
56
|
Jairaj, Sujata, Swati, HI 1544, HI 8498, MP 4010, GW 322, GW 366, GW 273
|
Moisture stress, High temperature
|
6
|
Peninsular Zone- PZ
(20° 04' 59'' N: 74° 07' 00'' E)
|
1.6
|
29.8
|
57
|
HD 2189, DWR 162, NI 5439, NIAW 917, MACS 6222
|
High temperature, Imbalanced fertilizer usage, Light soils, Moisture stress
|
12
|
Southern Hill Zone-SHZ
(11° 22' 12'' N: 06° 48' 00'' E)
|
0.1
|
10.1
|
9
|
NP 200, HW 5216, HW 1085, HUW 318, HW 517, HD 2135
|
Moisture stress, Acidic soils
|
9
|
Table 2
Population densities of culturable endophytic bacteria and the diversity indices of phylotypes of predominant cultural isolates
Wheat cultivation zone
|
Repertoire of culturesa
|
Diversity indices of phylotypes
|
Medium used
|
Shannon-H
|
Simpson-D
|
Chao1
|
Evenness J
|
NA*
|
TSA*
|
WFA*
|
Norther Hills Zone
NHZ
|
4.87±0.07ab
|
3.26±0.05bc
|
5.19 ±0.06b
|
2.72
|
0.95
|
16
|
0.95
|
North Western Plains Zone NWPZ
|
4.74±0.02b
|
3.77 ±0.04d
|
4.32 ±0.04a
|
2.66
|
0.95
|
15
|
0.95
|
North Eastern Plains Zone
NEPZ
|
5.02±0.11a
|
4.25±0.08a
|
4.40 ±0.05a
|
2.41
|
0.94
|
13
|
0.94
|
Central Zone
CZ
|
2.79 ±0.05e
|
3.12±0.03c
|
3.32 ±0.05e
|
2.08
|
0.96
|
14
|
0.96
|
Peninsular Zone
PZ
|
3.73±0.06c
|
3.38±0.05b
|
2.99±0.04c
|
2.74
|
0.96
|
16
|
0.96
|
Southern Hill Zone
SHZ
|
3.14 ±0.04d
|
3.26±0.03bc
|
4.06±0.07d
|
2.73
|
0.96
|
16
|
0.96
|
a Colony forming units are given as the mean Log CFU+ Standard Error. NA represents ‘Nutrient Agar,’ TSA-‘Trypticase Soy Agar,’ and WFA-‘Wheat Flour-based Agar.’ |
3.2.1. Qualitative and quantitative analyses of phenotypic traits
Of all the culturable WSEB isolates, about 83%, 72%, 72%, 46%, 41%, 57%, 67%, and 50% tested positive for the activities of amylase, cellulase, xylanase, protease, phytase, lipase, pectinase and esterase lytic enzymes, respectively. The endophytic isolates such as NH-4, PZ-24, CZ-39, and SHZ-35 were positive for all the lytic enzymes tested. Both NWP-9 and NWP-10 were positive for all the enzymes except lipase while the isolates such as NEP-19, PZ-27 and NWPZ-60 were positive for all the enzymes except esterase (Supplementary Table 1).
The potentials for hydrolytic enzyme production varied quantitatively among the representative isolates, randomly selected eight each from the wheat cultivation zones (Table 3). The amylase activities of the selected endophytic isolates were between 8.15 and 15.17 nkatal with the minimum by the isolates from NEPZ but the maximum by the isolates from the NWPZ. Highest activities of phytase were by the WSEB of the NHZ but those from the SHZ has the lowest. The pectinase activities were lesser in the representative WSEB from the NWPZ, followed by those from the NHZ. The representative WSEB from all other zones had comparable activities. The xylanase activities were between 2.09 and 5.69 nkatal while the higher activities were in both NHZ and PZ while the lower were in the NEPZ and NWPZ. The cellulase activities were the highest in the PZ (6.64 nkatal), followed by the SHZ and NWPZ. The protease activities ranged from 2.78 to 4.79 nkatal, maximum in the PZ, followed by the NHZ. The activities of esterase were more in the culturable bacteria of the NHZ while those in the NEPZ and CZ had lesser activities. The lipase activities were the highest in the culturable bacteria from the NHZ, followed by the NWPZ while those from all other zones had comparable activities.
Plant growth promoting traits and antagonism against fungal pathogens by WSEB
Traits such as N2 fixation, and the solubilization of phosphorous and potassium were qualitatively observed in about 46%, 78% and 9% of all the endophytic bacterial isolates from the wheat seeds, respectively. Likewise, the production of IAA, siderophores, ammonia and HCN production were observed in 89%, 11%, 72% and 59% of the isolates tested, respectively. The NWP-11 isolate possessed all the traits tested qualitatively. The isolate NWP-10 was positive for all traits except the ammonia production while the isolate NEP-22 had an exception of siderophore production(Supplementary Table 2).
The production of indole acetic acid was considerably higher in the culturable bacteria from the NWPZ than those from other zones. While the WSEB from the PZ had the least potentials, those from the NHZ and SHZ, and from the NEPZ and CZ had comparable potentials for the IAA production. The acetylene reduction activity (ARA), a quantitative estimate of N2 fixation of putative endophytes were between 0.81 and 32.06 nmol ethylene h−1 mg−1 protein(Table 4). The isolates from the NWPZ had the highest potential for the ARA. On the contrary, the endophytic isolates from the NHZ, followed by those from the CZ were the poorest in their potentials for the ARA. The potentials for the phosphate solubilization were higher in the isolates from the NHZ and NWPZ than those from other zones. The phosphate solubilization was the least in the WSEB from the SHZ (Table 4).
The endophytes showed variable antagonistic reactions against three potent fungal pathogens tested (F. graminearum, B. sorokiniana and T. indica). The endophytic isolates from each of these zones which reacted antagonistically to individual pathogens were tested for per cent inhibition. In general, the endophytic isolates from the NWPZ had higher levels of per cent inhibition against three pathogens tested(Supplementary Table 3).
Table 3
Hydrolytic enzyme activities of representative culturable endophytic bacteria isolated from different wheat cultivation zones
Wheat cultivation Zone
|
Amylase
|
Phytase
|
Pectinase
|
Xylanase
|
Cellulase
|
Protease
|
Esterase
|
Lipase
|
Norther Hills Zone- NHZ
|
12.61 ± 0.16c
|
5.65 ± 0.13d
|
2.98 ± 0.04b
|
5.69±0.09d
|
3.78±0.06a
|
4.33±0.07c
|
3.10±0.04c
|
6.55± 0.07c
|
North Western Plains Zone- NWPZ
|
15.17 ± 0.16d
|
3.79 ± 0.04c
|
1.01 ± 0.01a
|
2.54±0.03a
|
5.78±0.10c
|
3.30±0.05ab
|
2.54±0.04b
|
4.29±0.06b
|
North Eastern Plains Zone-NEPZ
|
8.15 ± 0.14a
|
3.00 ± 0.04b
|
3.25 ± 0.02c
|
2.09±0.01a
|
4.80±0.10b
|
3.55±0.03b
|
1.20±0.01a
|
2.34±0.03a
|
Central Zone-CZ
|
9.27 ± 0.17b
|
3.72 ± 0.04c
|
3.59 ± 0.05c
|
3.78±0.04b
|
4.69±0.08b
|
2.78±0.05a
|
1.88±0.02a
|
2.61±0.04a
|
Peninsular Zone-PZ
|
13.93 ± 0.24c
|
3.93 ± 0.08c
|
3.26 ± 0.03c
|
5.40±0.05d
|
6.64±0.09d
|
4.79±0.08c
|
2.42±0.04b
|
2.43±0.03a
|
Southern Hill Zone-SHZ
|
8.45 ± 0.09a
|
2.11 ± 0.03a
|
3.97 ± 0.06c
|
4.86±0.04c
|
5.86 ±0.05c
|
3.03±0.04a
|
4.00±0.07d
|
2.65±0.05a
|
The enzyme activity measurements of the endophytic isolates were expressed in nkatal. The mean values (+ S.E.) followed by the same letter(s) are not significantly different from each other (p = 0.05 ANOVA followed by the DMRT test). |
Table 4
Plant growth promoting traits and antagonism against fungal pathogens by the representative culturable endophytic bacteria isolated from different wheat cultivation zones
Wheat cultivation Zone
|
IAA1
|
ARA2
|
Phosphate solubilization3
|
Antagonism against4
|
F. graminearum
|
B. sorokiniana
|
T. indica
|
Norther Hills Zone- NHZ
|
91.75±1.25b
|
0.81 ±0.02a
|
415.03±5.25d
|
39.71 ±0.48c(7)
|
41.47 ±0.37c(6)
|
35.88±0.41ac(6)
|
North Western Plains Zone-NWPZ
|
203.02±3.38d
|
32.06 ±0.31e
|
440.25±5.63d
|
52.35 ±0.89d(7)
|
35.88 ±0.46ac (5)
|
43.38±0.53c(6)
|
North Eastern Plains Zone-NEPZ
|
152.25 ±1.75c
|
17.94 ±0.18d
|
203.25±3.38b
|
30.81 ±0.27ab (6)
|
29.85 ±0.38ab (6)
|
24.63±0.41ab (4)
|
Central Zone-CZ
|
158.36±1.63c
|
4.92 ±0.13b
|
219.25±2.01bc
|
26.25 ±0.31b (6)
|
34.78 ±0.54ac (5)
|
19.12±0.22b (6)
|
Peninsular Zone-PZ
|
68.88±0.63a
|
20.36 ±0.23d
|
226.25±1.50bc
|
38.97 ±0.42ac (6)
|
29.78 ±0.46ab (4)
|
33.75±0.41ac (6)
|
Southern Hill Zone-SHZ
|
111.25±1.63b
|
9.06 ±0.06c
|
146.75±2.50a
|
36.54 ±0.55ac (4)
|
24.56 ±0.49b (5)
|
35.15±0.32ac (3)
|
1IAA-Indole Acetic Acid in µg ml-1; 2ARA-Acetylene Reduction Assay- nmol ethylene h-1 mg-1 protein; 3P solubilisation µg ml-1; 4Antagonism in per cent (%). Values in parentheses are the number of isolates showed antagonism. The mean values (+ S.E.) followed by the same letter(s) are not significantly different from each other (p = 0.05 ANOVA followed by the DMRT test). |
The principal component analysis was performed for all the cultural phenotypes tested for the representative cultures from the six wheat cultivation zones. The analysis showed that the contributions of the first principal component (PC1) was 40.0% while that of the second component (PC2) was 26.25%, together accounting for 66.25% of variations(Fig.2). The principal component (PC1) explained the maximum variation with significant negative correlations in the activities of cellulase and pectinase while those related to the activities of esterase, xylanase, protease, lipase, and phytase, the production of IAA and phosphate solubilization, and the antagonistic activities against T. indica, F. graminearum, and B. sorokiniana were found positively correlated. The phenotypic trait accounting for PC2 is xylanase showing maximum variation over cellulase.
The enzyme activities of phytase, lipase, and xylanase had positive correlation whereas amylase was negatively correlated with PC2. The majority of traits (i.e., xylanase, esterase, protease, and lipase activities and antagonism against B. sorokiniana) that contributed to the separation of samples had higher coefficients with the NHZ; the activities of esterase and phytase had with the PZ. The representative endophytic bacterial cultures from the NHZ and NWPZ showed maximum variance through PC1 axis and are considered to be highly efficient with most of the enzyme activities and antagonistic activities against three pathogens tested while those in the NEPZ had less variance, on both sides of the axis, considered in having lesser efficiencies. The analysis also showed that the CZ, SHZ, and NEPZ are closely related and the NWPZ and NHZ as the most diverged from all other cultivation zones.
Colonization abilities of selected WSEB in wheat seedling roots
The root colonization abilities of three endophytic isolates of Bacillus and one each of Pantoea and Pseudomonas, and of the non-colonizing bacterial isolate of E. coli as the experimental control were assayed using the TTC staining. The wheat seedlings treated with the selected endophytic isolates showed differential responses as compared to the control and the untreated seedlings. The visual observations after 10 h of treating the seedlings with bacterial cells showed the motility, gathering around root tip first and stained as pink tips of root hairs (Fig. 3).On the extension of incubation for 24 to 48 h led to the increased colonization of root hairs and the whole root system. The inoculated bacteria were alive, colonizing the root tips first, and then entering into root hairs.The control treatment using E. coli had no pigmentation, neither stained after the TTC or H2O2 staining. The TTC-stained roots, root hairs, and root sections clearly showed that the bacterial cells of B. megaterium NEP-22, P. agglomeransNWP-9, and P. putida NWP-10 adhered closely, on the surface of epidermal parenchyma, entered intracellularly and subsequently colonized the roots (Fig. 3a, 3b, 3c).
The root cells showed the intracellular presence of bacteria, stained dark brown due to the H2O2 staining, in tissues under the bright field microscopy (Fig.3e, 3f). Through microscopic examinations of seedlings, we observed intracellular bacteria in the roots of selected species examined. Several bacteria were observed in various locations in seedling tissues, including root hairs, root epidermal cells, and root tips. In all treated seedling roots, bacteria were seen to be located within cells, located intercellularly in seedling roots and its hairs. Bacteria were observed in root tissues but could not be visualized in shoot tissues.
GFP-tagged bacterial endophytes in wheat seedlings
The GFP plasmid (pCambia1301) vector DNA isolated from E. coli was about 11 kb. The five predominant isolates of T. aestivum such as Bacillus cereus-NWPZ-5, Pantoea agglomerans-NWPZ-9, Pseudomonas putida-NWPZ-10, Bacillus megaterium-NEP-22, and Bacillus subtilis-PZ-23 were selected for the preparation of competent cells, and tagged with GFP plasmid vector through transformation process at different temperatures and timings. The bacterial colonies from the Kan50 NA plate were prepared for the confocal microscopy and the tagged cells were observed as in Fig.4. The colonization patterns of the tagged isolates differed in the tissues. Out of the isolates tested, two of them belonging to the most prevalent bacterial groups (Pantoeaspp. and Pseudomonas spp.) among all the wheat seeds were easier to transform than the other prevalent group of Bacillus spp. When the GFP-tagged endophytes were treated with wheat seeds, more light-green autofluorescence was observed in leaves, compared to root and stem of seedlings. In general, the GFP-tagged cells were more on 7 d after inoculation, than on 5 or 15 d after inoculation. After 15 d of inoculation, the endophytic colonization was more in inter- and intra-cellular spaces in roots, leaves, and the xylem vessels of the stem(Fig.4a, 4b, 4c, 4d, 4e). The matrix such as the soft agar medium did not affect the colonization pattern while the rapid spread in the vascular system suggesting the systemic colonization. The seedlings under the hydroponic condition led to greater colonization of GFP-tagged cells in stem, with lesser in root apices, root hairs and leaves. But the GFP-tagged bacterial cells colonized the zone of lateral root emergence, root tips and in the intercellular spaces of root epidermis.