- 1ICAR-KVK Tamenglong, Indian Council of Agricultural Research–Research Complex for North–Eastern Hill Region, Manipur Centre, Imphal, Manipur, India
- 2Department of Soil Science, Assam Agricultural University, Jorhat, Assam, India
- 3Division of Agronomy, Indian Council of Agricultural Research–Indian Agricultural Research Institute, New Delhi, India
- 4Division of Crop Production, Indian Council of Agricultural Research–Central Potato Research Institute, Shimla, India
- 5ICAR-National Bureau of Agriculturally Important Microorganism, Mau, India
- 6Department of Agricultural Chemistry and Soil Science, Bidhan Chandra Krishi Viswavidyalaya, Mohanpur, West Bengal, India
- 7Multi Technology Testing Centre and Vocational Training Centre, College of Agricultural Engineering and Post Harvest Technology, Central Agricultural University, Ranipool, Sikkim, India
- 8College of Agriculture, Central Agricultural University, Imphal, Manipur, India
- 9ICAR RC for NEH Region, Sikkim Centre, Tadong, Sikkim, India
- 10Regional Research Station, Indian Cardamom Research Institute, Spices Board, Tadong, Gangtok, India
- 11Amity Institute of Microbial Technology, Amity University Uttar Pradesh, Noida, India
- 12Central Horticultural Experiment Station, Indian Council of Agricultural Research–Indian Institute of Horticultural Research, Bhubaneswar, Odisha, India
- 13Department of Crop Physiology, Assam Agricultural University, Jorhat, Assam, India
- 14Farm Science Centre, Indian Council of Agricultural Research–Indian Institute of Horticultural Research, Kodagu, Karnataka, India
- 15ICAR-Indian Grassland and Fodder Research Institute, Jhansi, Uttar Pradesh, India
- 16Department of Agronomy, College of Agriculture, Central Agricultural University, Imphal, India
- 17ICAR-KVK Chandel, ICAR RC for NEH Region, Manipur Centre, Imphal, Manipur, India
- 18ICAR RC for NEH Region, Manipur Centre, Imphal, Manipur, India
Information on the role of boron (B) on soil physico-chemical and biological entities is scarce, and the precise mechanism in soil is still obscure. Present field investigation aimed to assessing the implication of direct and residual effect of graded levels of applied-B on soil biological entities and its concomitant impact on crop productivity. The treatments comprised of five graded levels of B with four replications. To assess the direct effect of B-fertilization, cauliflower was grown as a test crop wherein, B-fertilization was done every year. For assessment of succeeding residual effects of B-fertilization, cowpea and okra were grown as test crops and, B-fertilization was phased out in both crops. The 100% recommended dose of NPK (RDF) along with FYM was uniformly applied to all crops under CCOCS. Results indicated that the direct effect of B had the edge over residual effect of B in affecting soil physico-chemical and biological entities under CCOCS. Amongst the graded levels of B, application of the highest B level (2 kg ha–1) was most prominent in augmenting microbiological pools in soil at different crop growth stages. The order of B treatments in respect of MBC, MBN, and soil respiration at different crop growth stages was 2.0 kg B ha–1 > 1.5 kg B ha–1 > 1.0 kg B ha–1 > 0.5 kg B ha–1 > 0 kg B ha–1, respectively. Moreover, maximum recoveries of potentially mineralizable-C (PMC) and potentially mineralizable-N (PMN) were noticed under 2 kg B ha–1. Analogous trend was recorded in soil microbial populations at different crop growth stages. Similarly, escalating B levels up to 2 kg B ha–1 exhibited significantly greater soil enzymatic activities viz., arylsulphatase (AS), dehydrogenase (DH), fluorescein diacetate (FDA) and phosphomonoesterase (PMA), except urease enzyme (UE) which showed an antagonistic effect of applied-B in soil. Greater geometric mean enzyme activity (GMEA) and soil functional diversity index were recorded under 2 kg B ha–1 in CCOCS, at all crop growth stages over control. The inclusive results indicated that different soil physico-chemical and biological properties CCOCS can be invariably improved by the application of graded levels of B up to 2 kg B ha–1 in an acid Inceptisol.
Introduction
Boron (B) is a vital micronutrient that is indispensable for proper crop growth (Bhupenchandra et al., 2021a). Boron is a necessary micro-element for plant cell wall structural integrity and is involved in various plant processes like cell division, calcium utilization, pollen production, and anther development during the reproductive phase (Nadeem et al., 2019). Currently, the B deficits in soils are widespread globally causing B micronutrient malformations that impinge on agricultural production (Shorrocks, 1997; Liu, 2000; Choudhary et al., 2014; Behera et al., 2022). Due to excess deficiency symptoms manifested as an implication of B deficit, the assessment of key functions of B in plants has long been the main concern from a nutritional point of view. For plants, managing B is difficult because the optimum B range is narrow which can fluctuate from soil to soil (Gupta, 1993; Marschner, 1995). Normally, B averaged nearly 30 mg kg–1 soil depending on the main rock wherein its content varies extensively. Satisfactory B content for flora in soils is more or less 25 mg kg–1 (GreenFacts, 2002). Soil microbial biomass holds a vital role in nutrient-cycling, plant-pathogen suppression, the disintegration of debris, and decay of pollutants establishing the vibrant living entity of soil, and thus, attributing to ecological sustainability owing to their diverse existence, enormous effective genetic pools, catabolic adaptability, and stress tolerance ability in a holistic manner (Deluca et al., 2019). The dimension and activity of the microbial biomass determine the nutrient availability and production potential of the agro-ecosystems (Friedel et al., 1996; Singh et al., 2020, 2021). Therefore, it becomes obligatory to determine microbial biomass nitrogen (MBN). Since it becomes vital for the quantification of N-dynamics in agro-ecosystems as it controls the soil inorganic-N accessibility and loss and its contribution to the primary N-sources of potentially mineralizable-N (PMN) in the soil (Bonde et al., 1988). Microbial biomass carbon (MBC), MBN, and microbial respiration, have further garnered added interest owing to their sensitivity to crop management practices than the bulk soil organic matter (Awale et al., 2017). Soil microbes, the existing fraction of soil organic matter (SOM) acts as a transitory nutrient-sink and are accountable for unleashing nutrients from SOM for exploitation by plants (Bollag et al., 2017). Basal respiration (BR) and C-mineralization are ample indicators of microbial activity, depending on the substrate accessibility and the soil edaphic environment (Balota et al., 2003). On the whole, the CO2 respired during a year in terrestrial ecosystems is the consequence of C-mineralization of the minute active fraction pools, which are mainly accountable for unleashing nutrients in the soil (Brdar-Jokanović, 2020).
It is implicit to address that B acts a crucial function in the biological activities of living organisms as proven earlier by establishing the necessity of B for diatoms and cyanobacteria (Bonilla et al., 1997). Soil microorganisms are by far the most important producers of soil enzymes that perform many ecological processes such as bio-geochemical cycling and decomposing pollutants and debris from flora and fauna and the microbes (Goswami et al., 2017; Furtak and Gajda, 2018). Soil microbes are accountable for the transformation of SOM and soil nutrients (Mooshammer et al., 2014; Singh et al., 2020, 2021). The microbes and their enzymatic outputs are indispensable to plants, while plant roots generate organic substances that are vital to the populace expansion of microbes (Jjemba and Alexander, 1999). Soil microbial indices is touted as a substitute for organic carbon cycling and its related nutrients viz., N, P, and S, signifying that elevated microbial action implied increased soil productivity and vice-versa (Pavan et al., 2005). Soil enzymes are touted as an index for examining the activities of microbes, soil productivity and soil quality owing to the symbiosis of microbes and flora (Dick, 1994; Bandick and Dick, 1999). Evaluation of diverse soil extracellular enzymes established it as a potent means for assessment of the soil functions for nutrient-cycling and microbial nutrient requirements (Sinsabaugh et al., 2012). We also hypothesized that in lieu of specific enzyme activity, an index merging diverse enzymes would be a more efficient and appropriate indicator of soil quality, since it could specify an inclusive diversity of soil functions. Relationships between crop yields, nutrient availability, and these enzyme activities are also obscure, as information on the enzymes’ ability to predict soil quality attributes is scarce. With the induction of intensive crop management practices like fertilizer application exhibiting complicated and harmful implications on plants and microbial associations, the studies on impact assessment of B-fertilization on soil microbes and enzymes for recuperating agricultural output in B-deficient soils become highly imperative (Dick, 1994; Tabatabai, 1994). The necessity of B-fertilization for the augmentation and maturity of plants has already been established (Shelp, 1993; Marschner, 1995). Boron acts as a vital function in the translocation and assimilation of complex carbohydrates in the plant, production of plant hormones and nucleic acids, germinating pollen, flower induction and fruiting, and water utilization. The main significant roles of B in plants are its structural role in cell wall growth and stimulating or inhibiting of precise metabolic pathways (Ahmad et al., 2009). In addition, B plays a crucial role in N assimilation, N fixation, and the growth of legume root nodules (Bolaños et al., 2004a; Bellaloui et al., 2014).
Boron accessibility depends on many criteria that exist in the soil–plant system such as SOM, soil texture, cultivation, soil moisture, temperature, soil pH and liming, soil fertility, and microbial activity (Shorrocks, 1997; Kumar et al., 2016; Shireen et al., 2018). Microbes assimilate SOM, which in turn, helps in releasing the B from organic complexes. Boron present in the soil is considered as a vital fraction related to SOM and is unleashed via microbial activities (Berger and Pratt, 1963). Despite the prime role of B on floral growth and functioning, no credibility has been established to explain that B is an enzyme component and possesses a direct role in enzyme actions. It is obscure to claim that these processes are precursor of the direct functioning of B or the indirect role of B. The biological effects of B are better understood in plants, where it has been proven that it can influence physiology and biological activities (Grattan et al., 2015). Though the impact of B on the soil microbial community is little known, there is a paucity of information on the element’s consolidative impacts on several characteristics of the soil microbial community, including activity, biomass, and diversity. Most of the soil fertility experiments stressed the aspects of the changes in soil chemical pools without giving much attention to biological attributes in soil rendering a lacuna on comprehensive fertility evaluation in soil fertility research. To date, the study on the role of B on soil biological attributes is obscure and very limited. Consequently, advanced research is necessitated to have a clear insight into the functioning of B in plant development and soil biology. The study hypothesized B-fertilization would improve the soil physio-chemical and biological properties and productivity of cauliflower-cowpea-okra cropping system. The objective of the study is to assess the impact of graded levels of B-fertilization on the soil physio-chemical and biological properties and their relationships in a cauliflower-cowpea-okra cropping system in North East India.
Materials and methods
Experimental site
A field experiment was conducted on cauliflower-cowpea-okra cropping system at the Horticultural Experimental Farm, Assam Agricultural University, Jorhat, India (26°47′N latitude, 94°12′E longitude, 86.6 m altitude) during 2015–2017. The climate of the experimental site is sub-tropical with hot humid summers and comparatively dry and cool winters. Normal annual rainfall varies between 1,500 and 2,000 mm. Usually, rain commences from June and continues up to September with the pre-monsoon showers commencing from mid-March. The highest temperature of 34°C during summers and the lowest about 7°C during winters is usually prevalent. Agro-meteorological information is presented in Figure 1. The soil of the experimental site is Inceptisol having a sandy clay-loam texture with pH 4.8.
Figure 1. Meteorological data during the period of investigation for the year 2015–2017. Source: Department of Agrometeorology, Assam Agricultural University, Jorhat, Assam. The vertical bars represent the SE(m)±. T, temperature; RH, relative humidity; and BSSH, bright sunshine hour.
Experimental design and treatments
The field experiments were uniformly laid-out for all the three crops CCOCS in a completely randomized block design with 4-replicates. Soil application of B was imposed at the rates of 0, 0.5, 1.0, 1.5 and 2.0 kg ha–1 in cauliflower, wherein its application was restricted only to cauliflower in both years of experimentation to assess the direct effect of B fertilization. However, B-fertilization was exempted in succeeding crops (cowpea and okra) to assess the residual effect of B-fertilization in the sequence. Borax (Na2B4O7.10H2O, analytical reagent grade with 10.5% B) was applied as the source of B for soil application. The 100% recommended dose of fertilizers (RDF) (supplied via urea, single super phosphate, and muriate of potash) and farmyard manure (FYM, well-decomposed cowdung) was uniformly applied to all the crops at the time of sowing (Supplementary Table 1). FYM used in the experiment had the bulk density of 0.24 Mg m–3, pH (7.7), N (1.4%), P (0.34%), K (0.8%), Mg (0.5%), Ca (1.4%), and C: N ratio (28:1) as determined using standard procedures (Rana et al., 2014).
Soil sampling and analyses
For soil physio-chemical analysis, soil samples were collected from a depth of 0–15 cm at different crop growth stages in the sequence. While for soil biological properties, moist soil samples were collected at the initial and different crop growth stages of cauliflower, cowpea and okra for the two years. Soil samples were then stored in the refrigerator at 4°C for analysis of biological parameters. All the analyses were made in triplicate following standard protocols. The soil samples were analyzed for soil physico-chemical and biological properties (Supplementary Table 2).
The microbial quotient was calculated as the ratio of MBC to SOC and expressed in percentage (Anderson and Domsch, 1989).
Crop equivalent yields
After the harvest of each crop, the yield was recorded. Crop equivalent yield (CEY) was computed to evaluate system performance after converting the yield of one crop (assumed as x) into the equivalent yield of another crop (assumed as y) on a pricing basis:
Where, EYx = yield of x crop converted to yield of y crop, Yx is the yield of x crops (kg ha–1), Px is the price of x crops (US$ kg–1), and Py is the price of y crop (US$ kg–1). All the yields of crops were converted to the equivalent yield of cauliflower, which was planted first in the cropping sequence.
Enzyme activity-based index for calculation of soil functional diversity
Soil quality index
Geometric mean of enzymatic activities
To better explicate the impact of B on soil enzyme activities, we computed the geometric mean of enzymatic activities (GMEA), as it can replicate the inclusive enzyme activity levels (Hinojosa et al., 2004). GMEA is a consolidative method to pool the enzyme activities associated to diverse soil functions and nutrients; therefore, possibly it will reflect soil quality index (Paz-Ferreiro et al., 2012).
GMEA of the assayed enzymes was computed for each sample as:
Soil functional diversity
The Shannon index and Simpson-Yule index cater the data about the spread or distribution of C source usage by the microbial community (Kumar et al., 2017). The ensuing indexes can be used for quantification of richness, evenness, and diversity of the soil microbial community.
It was computed using the following equations
Shannon’s diversity index (H)
where, Pi is the ratio of each enzyme activity to the summation of whole enzymes activities for a specific sample. Enzyme activities were deciphered as μg product formed per g of soil per hour.
Simpson-Yule index (SYI)
SYI was calculated as for each sample as:
The diversity of the community is directly proportional to H.
Multivariate analysis
Hierarchical cluster analysis
The data obtained on the different groups of biological entities present in the soil samples were subjected to agglomerative hierarchical cluster analysis (HCA) (Euclidean distances, Ward’s agglomeration rule) to establish homogeneous groupings of data. The nodes depicted clusters retrieved on each step of hierarchical clustering.
Principal component analysis
The principal component analysis (PCA) of all the data was performed (Andrewsi and Carroll, 2001; Andrews et al., 2002) to ascertain the variability and show the relationship among the various soil properties, and to extract the dominant principal components from the whole data set in soil resorting to R studio. PCA is a multivariate statistical dimension reduction tool that resorts to an orthogonal transformation to transform a set of correlated variables to linearly uncorrelated variables known as principal components (PC). The extracted results of a PCA are displayed in terms of component scores, also called factor scores and loadings (Wold et al., 1987). In the PCA algorithm, diminution of the number of components was yielded via the eigenvalue-one criterion i.e., eigenvalue >1 is retained, also called as Kaiser criterion (Kaiser, 1960), and the scree test (Cattell, 1966).
Stepwise multiple linear regression analysis
Stepwise multiple linear regression analyses were performed using R studio with the backward exclusion method to explore the significance of dominant soil biological entities in the prediction of crop yield in the sequence. The relationship between a single response variable (dependent variable) and two or more controlled variables was evaluated using multiple stepwise linear regressions (MLR) (independent variables). MLR used in the research states that the higher R2 generates good results in model fitting (Bowerman et al., 2005). The analytical model, used to develop a model for predicting crop yield from the biological attribute’s relationship is given by the equation:
Where, βs are coefficients, Xi are the predictors, Y is the crop yield (response) and β0 is a constant.
The null hypothesis of sequential uncorrelated errors was tested independently on regression residuals using the Durbin–Watson statistic.
Path analysis (causal modeling)
Path analysis is a standardized partial regression analysis used to determine the significance of the relationship between sets of variables and to provide estimates of the magnitude to make the multiple regressions easier to comprehend. It also helps to figure out the direct, indirect and total impact of predictor variables on the response variable. An evaluation of correlation does not specify the precise influence of the attributes to crop yield and this correlation can be segregated into direct and indirect effects via path coefficient analysis. It permits the separation of the direct and their indirect effects via other traits through allocating the correlations (Wright, 1921) for clarity of explanation of cause and effect.
Data analysis and visualization
The experimental data obtained from different observations were analyzed statistically by using Fisher’s method of ANOVA in randomized block design (Panse and Sukhatme, 1985). Significance or non-significance of the variance due to different treatment effects was estimated by computing concerned ‘F’ values. At a 95% confidence level, the experimental means were compared. To compare treatment means, the Duncan Multiple-Range-Test (DMRT) was employed. Univariate Pearson’s correlation analysis was executed to determine the interrelationship between biological entities in the soil samples and crop yield. A correlograms was built using the “corrgram package” in R Studio.
Results
Crop yield and cauliflower-equivalent-yield
In general, the crop yields in the cauliflower-cowpea-okra cropping system were significantly impacted by the imposition of a graded level of B, as evident by the significant augmentation in yield (Figure 2A). A satisfactory cauliflower curd yield (highest) of 23.25 Mg ha–1 was obtained as a ramification of the direct effect of 2 kg B ha–1 imposition in cauliflower with drastic yield augmentation up to 21.8% over the control (19.09 Mg ha–1). While the residual implication of 2 kg B ha–1 in cowpea and okra, also leveraged the pod yield to the tune of 7.15 Mg ha–1 and fruit yield of 20.61 Mg ha–1, thereby, improving the crop yield by 25.7% over the control (5.69 Mg ha–1) in cowpea and 21.2% over the control in okra (17.0 Mg ha–1). Likewise, the CEY of the cauliflower-cowpea-okra cropping system was significantly greater (p < 0.05) with an imposition of 2 kg B ha–1 as compared to control (Figure 2B). The extent of growth in equivalent yield was 25.7, 21.3, and 22.5% for CEY of cowpea, CEY of okra and total CEY, respectively over control. The increase was always higher with 2 kg B ha–1 than the rest of the B levels including control.
Figure 2. Effect of graded B levels on crop yield (Mg ha–1) in cauliflower-cowpea–okra cropping system (A); effect of graded B levels on CEY (Mg ha–1) in cauliflower-cowpea–okra cropping system (B).
Soil physico-chemical properties
Imposition of differential B levels could not significantly (p ≤ 0.05) affect the soil BD under cauliflower-cowpea-okra cropping system (Table 1), however, a decrease in BD was noticed in B applied plots as compared to control (B0). In general, the experimental soils were mostly acidic in reaction with pH 4.83–4.89 (Table 1), with no significant difference among the B treatments. Continuous two years of experimentation decreased the BD in soil by 2.4% as compared to the initial status on addition of a maximum B level of 2 kg ha–1. However, improvements in soil pH and SOC were noticed to the tune of 1.5 and 30% over the initial status when 2 kg B ha–1 was applied. The available N and P (Table 1), were also significantly (p ≤ 0.05) influenced by the addition of graded B levels wherein the greatest value was always with 2 kg B ha–1 application.
Table 1. Effect of graded B levels on soil physico-chemical properties at different crop stages in cauliflower-cowpea-okra cropping system (Pooled data 2015–017).
Implication of B on temporal dynamics in microbiological pools of soil organic matter
The soil microbial biomass carbon (MBC), Microbial quotient (MBC: SOC) microbial biomass-N, (MBN), Microbial biomass-P (MBP), and soil respiration (SR) improved significantly (p ≤ 0.05) on the application of different B doses at all the stages of crop growth over the control with few exceptions in case of MBN (Table 2). The highest values of MBC, MBC: SOC MBN, MBP, and SR in all the crop growth stages in different crops were with 2.0 kg B ha–1 addition. After 2-years of experimentation, there was a gain of 24.5, 12.1, 54.2, 34.4, and 36.1% in MBC, MBC: SOC MBN, MBP, and SR over the initial soil status. The order of the B treatments in respect of MBC, MBN, and SR at different crop growth stages was 2.0 kg B ha–1 > 1.5 B kg ha–1 > 1.0 B kg ha–1 > 0.5 B kg ha–1> 0 kg B ha–1 (Control), respectively.
Table 2. Effect of graded B levels on microbiological pools of soil organic matter at different crop stages in cauliflower-cowpea-okra cropping system (Pooled data 2015–2017).
Potentially mineralizable-C and potentially mineralizable-N
A significant (p ≤ 0.05) increasing trend in potentially mineralizable-C (PMC) and PMN with the B application rate was noted in the soil of CCOCS (Table 3). Comparatively, the higher recoveries of PMC and PMN were always higher with 2.0 kg B ha–1 direct and residual impact of B-fertilization in CCOCS at all the crop growth stages. The PMC and PMN across the crops and growth phases in the sequence were ranked as 2 kg B ha–1> 1.5 B kg ha–1> 1.0 B kg ha–1> 0.5 B kg ha–1> 0.0 kg B ha–1 (Control). Compared with the initial soil status, there was an improvement of 26.3 and 52.5% in respect of PMC and PMN content in soil due to the imposition of graded levels of B.
Table 3. Effect of graded B levels on PMN and PMC content in soil at different crop stages of cauliflower-cowpea-okra cropping system.
Microbial populations
There was significant (p ≤ 0.05) improvement in microbial populations (actinomycetes, bacteria and fungi) across B treatments in all crops in the sequence (Table 4). The addition of escalated B level of 2.0 kg B ha–1 led to a significant improvement in the status of microbial populations in the soil at different crop growth stages in cauliflower (direct application of B); and that in cowpea and okra (residual effect of B). Interestingly, this escalated B treatment led to an augmentation of 54.2, 55.3 and 53.7% of actinomycetes, bacteria and fungi population in comparison to the initial soil status.
Table 4. Effect of graded B levels on microbial populations in soil under cauliflower-cowpea-okra cropping system (Pooled data 2015–2017).
Soil enzymes
Soil enzymes viz., AS, DH, FDA, and PMA showed a significant (p ≤ 0.05) response to graded B levels in soil (Table 5). However, the urease enzyme (UE) showed an antagonistic effect, thereby, exhibiting a reciprocal response to the appliance of the graded levels of B. Across all crop growth phases in the sequence, higher enzyme activities were noticed under plots receiving higher B levels of 2.0 kg ha–1. The status of the content of soil enzymes in diverse crops and growth stages varied as: 2.0 kg B ha–1> 1.5 B kg ha–1> 1.0 B kg ha–1> 0.5 B kg ha–1> 0.0 kg B ha–1. Soil enzymes’ activity registered an increment to the tune of 44.7 (AS), 45.1 DH), 38.6 (FDA) and 46.7% (PMA), respectively over the initial soil status. On contrary to that, UE enzyme activity exhibited a decrement by 20.7% over the initial soil status.
Table 5. Effect of graded B levels on soil enzymes in cauliflower-cowpea-okra cropping system (Pooled data 2015–2017).
Soil quality index and functional diversity index
The GMEA was significantly (p ≤ 0.05) higher under plots receiving the highest B levels of 2 kg ha–1 at all the crop growth stages under CCOCS (Table 6). Irrespective of crop growth stages, direct application of 2 kg B ha–1 in cauliflower recorded comparatively higher GMEA than that under cowpea and okra (residual effect), respectively (Table 6). Similarly, the functional diversity indexes (H and SYI) exhibited the similar trend wherein their higher values were observed under plots receiving 2 kg B ha–1 (Table 6).
Table 6. Geometric mean enzyme activity (GMEA), Shannon diversity index (H) and Simpson diversity index (SYI) as affected by graded B levels under cauliflower-cowpea-okra cropping system (Pooled data 2015–2017).
Correlation between the soil properties and mean yield of cropping sequence
The univariate correlation coefficients (r) in between the soil properties (physico-chemical and biological) and mean yield of cropping sequence (MYCS) are illustrated by the correlogram [Auto correlation function (ACF) plot] (Figure 3). In general, the results exhibited an existence of a significant positive correlation (p < 0.01 and p < 0.01) between the soil properties and MYCS barring MBN, at different crop growth stages, thereby signaling a synergistic relationship between them. However, BD and UE activity exceptionally showed a negative non significant correlation with the rest of the parameters, while both are positively and strongly correlated to each other (p < 0.01). Selectively, pH is highly correlated (p < 0.01) with SOC, available N and P, MBC, MBP, AP, FDA, and MYCS but positively correlated (p < 0.05) with PMN. Similarly, SOC showed a highly significant positive correlation (p < 0.01) with available N and P, FDA but exhibited a positive correlation (p < 0.05) with MBP, AP, and MYCS, respectively.
Figure 3. Correlograms of Pearson correlation coefficients (r) matrix between the soil physico-chemical and biological entities under cauliflower–cowpea–okra cropping system. The correlation coefficient (r) values are significantly positive at p < 0.01 (**) and p < 0.05 (*) levels of probability (2-tailed); the color assigned to a point in the correlograms grid indicates the strength of a correlation between the soil biological entities, and r values correspond directly to the color codes ranging from red to blue, respectively. Right and left tilted ellipse in the correlograms grid indicate positive and negative correlation, respectively. where, Bd, bulk density; SOC, soil organic carbon; MBC, microbial biomass carbon; MINN, mineralizable nitrogen; AP, actinomycetes population; BP, bacterial population; FP, fungal population; AS, arylsulphatase activity; DHA, dehydrogenase activity; FDA, fluorescein di-acetate hydrolysis activity; PMA, phosphomonoesterase activity; UE, urease activity.
Clustered analysis
Hierarchical clustering (Figure 4) in respect of soil biological entities under different crops in CCOCS identified the distinct clusters based on similarity in function and other relevant biological attributes. In respect of cauliflower, three distinct clusters were formed viz., Cluster-I (PMC and PMN), Cluster-II (Microbiological pools of soil organic matter: MBN, MBC, and SR) and Cluster-III (Microbial population: AP, FP, and BP; Soil enzymes: AS, PMA, FDA, and DH). Similarly, in the case of cowpea, three distinct clusters were generated viz. Cluster-I (MBN, MBC, and SR), Cluster-II (PMC, PMN, and PMA), and Cluster-III (Microbial population: FP, AP, and BP; Soil enzymes: FDA, DH, and AS). Likewise in okra, similar clusters were formed viz. Cluster-I (PMC and PMN), Cluster-II (Microbiological pools of soil organic matter: MBN, MBC, and SR) and Cluster-III (Microbial population: FP, AP, and BP; Soil enzymes: FDA, DH, AS, and PMA). However, the UE enzyme formed a discrete outlier as this enzyme had a reciprocal relationship with the examined parameters (Figure 4).
Figure 4. Hierarchical clustering of the soil biological entities indicating similarities in different soil physico-chemical and biological entities in cauliflower-cowpea–okra cropping system. MBC, microbial biomass carbon; SR, soil respiration; MINC, mineralizable carbon; MINN, mineralizable nitrogen; AP, actinomycetes population; BP, bacterial population; FP, fungal population; AS, arylsulphatase activity; DHA, dehydrogenase activity; FDA, fluorescein di-acetate hydrolysis activity; PMA, phosphomonoesterase activity; UE, urease activity; VS, vegetative stage; CIS, curd initiation stage; MS, maturity stages are the Blevels, respectively.
Principal component analysis
Principal component analysis executed in respect of soil physico-chemical and biological properties in CCOCS extracted three principal components with eigenvalues equal or greater than unity (Supplementary Table 3), accounting cumulatively up to 95.56% of the total variance since they possessed eigenvalues >1.0 and explained >5% of the variance in the total dataset of the available data (Supplementary Table 3). The loading plot (Figure 5), generated three PCs with eigenvalues equal or greater than unity viz. PC1 (68.7%), PC2 (14.8%), and PC3 (8.2%), respectively. Barring the soil physiochemical properties, the loading plot (Figure 5) (denoted by blue lines), elucidated that PC1 had large positive loadings on BP, MBC, MBN, MBP, PMC, and SR and subsequently followed by soil enzymes, and they were highly correlated to each other. Similarly, PC 2 exerted higher loadings on PMC and PMN, whereas PC 3 had heavy loadings on MBN and UE (Figure 5), respectively. Contrarily, UE activity was negatively correlated with PC1 which is attributed to the reciprocal relation with added corresponding B levels. In case of PC2, it showed heavy loading on PMA and PMC. The respective score plots (denoted by red colored dots) of the crops in CCOCS were divided into four quadrants (I, II, III, and IV) based on component (1 and 2) scores (Figure 5) to allow for better visual discrimination of B levels on soil physico-chemical and biological properties in CCOCS. The scoreplot showed that the first quadrant identifies that the B levels of 1 and 1.5 kg ha–1 displayed positive heavy loading on some PC1 components viz. soil chemical properties (SOC and pH), soil enzymes (FDA, PMA, DH, and AS) and microbial population (AP and FP), respectively. Similarly, the 2nd quadrant, indicated that B level of 0.5 kg ha–1 had greater loadings on BD and UE wherein both these variables had reciprocal relation with the rest of the studied parameters. The 3rd quadrant was occupied by control (B level of 0 kg ha–1) which did not influence any soil parameters. Interestingly, the 4th quadrant harbored the most important parameters influencing the crop yield in the sequence namely BP, MBC, MBN, MBP, PMN, PMC, SR, respectively in the biplot (Figure 5), which in turn, was affected by higher B level of 2 B kg ha–1.
Figure 5. Three-dimensional graphical biplot showing the loading and score plot formed by principal components 1, 2, and 3 with different soil physicochemical and biological entities in cauliflower-cowpea–okra cropping system. Percentage values on PC1, PC2, and PC3 indicate the respective variance explained by the first three PCA axes; where, MBC, microbial biomass carbon; MINN, mineralizable nitrogen; AP, actinomycetes population; BP, bacterial population; FP, fungal population; AS, arylsulphatase activity; DHA, dehydrogenase activity; FDA, fluorescein di-acetate hydrolysis activity; PMA, phosphomonoesterase activity; UE, urease activity.
Stepwise multiple linear regressions for predicting the best model for crop yield
The stepwise multiple linear regressions (SMLR)exercised on MYCS showed the best fitting model that may produce maximum yield is enlisted in Table 7. The results of SMLR implied that the retained biological entities BP, MBC, MBN, MBP and PMC (Table 8) were the best predictors contributing 44.8, 22.17, 18.67, and 14.18%, respectively to the MYCS. Fitting MYCS as a dependent attribute (response variable) and biological properties as the independent attributes (predictor variables) (Eq. 1), a best-fitting regression model was generated below as:
Table 8. Stepwise regression variances analysis of different soil biological properties in cauliflower-cowpea-okra cropping system.
Path analysis
For the execution of path analysis, all the biological entities were assigned as predictor variables excluding soil physiochemical properties, while the MYCS under CCOS was assigned as a response variable to avoid the redundancy of data. In path analysis, the magnitude of the contribution of all the biological entities to MYCS was quantified by its corresponding path coefficient values. Results of the path analysis (Figure 6) showed that BP with path coefficient values of 2.09 had the highest and most significant direct effect on MYCS and had twelve numbers of indirect effects emanating from the rest of the twelve biological parameters under investigation. Barring, urease enzymes, these twelve biological parameters indirectly contributed to MYCS by largely linking to the BP and their indirect path coefficients through BP. The contributions and impact of the biological entities to MYCS can be ranked in decreasing order as BP > MBC > MBN > MBP > PMC > PMN > SR > DH > AP > FP > AS > FDA > UE > PMA, respectively.
Figure 6. Path diagram depicting the contribution of soil biological entities to the MYCS in cauliflower-cowpea–okra cropping system. Single-headed arrows, double-headed arrows and connectors signify the path coefficient (β) (direct effect), simple correlation coefficients between variables and mutual association, respectively; where, MYCS, mean yield of the cropping system; MBC, microbial biomass carbon; MBN, microbial biomass nitrogen; SR, soil respiration; PMC, potentially mineralizable carbon; PMN, potentially mineralizable nitrogen; AP, total actinomycetes population; BP, total bacterial population; FP, total fungal population; AS, arylsulphatase activity; DH, dehydrogenase activity; FDA, fluorescein di-acetate hydrolysis; PMA, phosphomonoesterase activity; UE, urease activity.
Discussion
Boron is one of the indispensable nutrients for the ideal growth, development, produce, and quality of crops (Shireen et al., 2018). In general, B being a vital nutrient plays a role in plant growth, phenols, lignification, tissue expansion, membrane-related reactions, ribose nucleic acid (RNA) metabolism, hydrocarbon metabolism, pollen germination and seed development which are directly implicated in increasing crop yield (Goldbach and Monika, 2007). The enhancement in crop yield as a result of B-fertilization could be ascribed to the improved availability and accessibility of nutrients to plants (Kumar et al., 2016, 2017), hence producing and mobilizing surplus carbohydrates and proteins along with its role in enhancing their translocation from the site of synthesis to the storage organs (Takkar and Randhawa, 1978; Verma et al., 2012). Moreover, B acts as a key role in many metabolic processes such as cell wall differentiation, cell development, N-metabolism, fertilization, fat metabolism, hormone metabolism, active salt absorption, and photosynthesis (Nason and McElory, 1963), which in turn contributed to higher fresh and dry matter yield of cauliflower. Similar findings in okra were also reported earlier (Saha et al., 2010; Rahman et al., 2017). Kumar and Sen (2004) reported that the application of B and Zn improved the yield and quality of okra seed. The beneficial impacts of B on curd quality and yield of cauliflower were acknowledged by Gupta (1993).
The decrease in BD under CCOCS might be due to an improvement in soil structure and porosity due to the addition of FYM. A slight increment in soil pH could be possibly due to the ligand exchange between OH–groups of soil Al and Fe(OH)2 and organic compounds, and the disintegration and binding of organic complexes of the applied FYM (Xu et al., 2006). Similarly, SOC in soil remained significantly (p ≤ 0.05) (Table 1) unaffected after two years of B application, however, an increment in SOC was observed possibly due to the SOC build-up through regular addition of FYM and desirable changes in biochemical and physical properties of soil (Ghosh et al., 2012; Bhupenchandra et al., 2022; Harish et al., 2022; Kumar et al., 2022). Also, another reason for the improvement in SOC could be due to the formation of a strong diol complex of B with organic matter in soils and the capacity of organic matter to improve CEC of soils (Bhupenchandra et al., 2021b). Increase in available N could be due the release of mineralized N by the addition of organic matter along with the concurrent release of N via symbiotic biological N fixation by cowpea roots, since B plays a vital role in biological N fixation and upsurges the number of effective nodules (Bolaños et al., 2001). The improvement in the status of available P could be explicated by the existence of positive interaction between P and B in the soil as both are in anionic forms and might have been involved in anion exchange (Bhupenchandra et al., 2021b).
Microbial biomass carbon is the measure of the C present within the living constituent of soil organic matter. Soil respiration (SR) is the CO2 released by the biological activities of soil organisms, involving plant roots, microbes, and soil animals are usually calculated as a flux of CO2 from the soil surface. Escalating B appliance quickly altered soil MBC content and soil respiration (Bilen et al., 2011). Improvement in MBC and CO2-C production in soil with the application of graded B could also be due to the continuous application of FYM in all the crops for two years in the sequence, which in turn, might have created a conducive atmosphere for intensified microbial activities in soil (Singh et al., 2020, 2021). Soil MBC, being an active and labile component of SOC dependent on the SOM (Chen et al., 2005; Kumar et al., 2022). There is no direct evidence of the effect of B application on microbial properties of soil. But, the enhancement in these biological properties with the incremental doses of B indicated a possible relationship between the B-fertilization and microbial activity of the soils. Subsequently, microbes are implicated in the assimilation of SOM, which further led to the release of B from organic complexes in soil (Kumar et al., 2016, 2017). Also, total B present in the soil fraction is closely related to SOM and was unleashed via microorganism action (Berger and Pratt, 1963). Upon intensifying the B appliance, the soil microbial biomass load in the soil quickly transformed and escalated vibrantly. Moreover, the accessibility of readily mineralized C and N, and improvement in soil physico-chemical properties might have enhanced the microbial population in soil (Bhardwaj and Datt, 1995; Kumar et al., 2022). Application of FYM improves the SOC pool by supplying organic matter in greater mineralizable form, thereby, delivering substrate for microbial utilization and this could be the cause of higher MBC and SR in the current study (Goldberg, 1997; Liu, 2000; Rajpoot et al., 2021).
The PMC, also known as biodegradable C, is the entirety of organic matter which can be decomposed through microbial action (Guo et al., 2019). Potentially mineralizable-N (PMN) is a measure of the active fraction of soil organic-N, predominantly accountable for the discharge of mineral-N via microbial accomplishment (Campbell and Curtin, 2007). The PMN is availed to plants and microorganisms in the form of NO3– by aerobic mineralization. It is a fraction of N linked to the microbial biomass which is positively related to MBC. With crop growth, PMN content in soil is augmented as it is the quantum of N that mineralize with time at the most favorable temperature and moisture. It comprises a diverse group of organic complexes which encompass microbial biomass, crop residues and humus. The increment of PMC content might be attributed to the soil application of well-decomposed FYM which acts as a substrate for microbial entities. The enhanced N-mineralization was observed during the symbiotic biological N-fixation by cowpea roots since B acts a vital role in the biological fixation of N and augments the number of effectual nodules (Bonilla et al., 1997; Bolaños et al., 2001), and hence, might have created signaling compound through the rhizobia infection on roots of legume crop (Dénarié and Cullimore, 1993; Spaink, 2000).
It is established that B is vital for symbiont/plant signaling, namely nod-gene activation by root plant exudates and nodule invasion (Redondo-Nieto et al., 2001). Moreover, B is essential for infectivity thread advancement and nodule initiation (Bolaños et al., 1996) due to the function of B as a modulator of the interactions amidst plants derived infection thread matrix glycoproteins and the bacteria cell surface. Boron stabilizes membranes, which aids the relationship between bacterial cell surfaces and the peri-bacteroid membranes, helping them in regulating symbiotic setup (Bolaños et al., 1996). Specifically, B is indispensable for the target of nodule-specific plant-derived glycoproteins that are critical for signaling bacteroid differentiation into a N2-fixing form (Bolaños et al., 2004b). Thus, all these factors created a conducive atmosphere for augmenting the microbial population in the soil. Several B tolerant bacterial strains belonging to the genus Bacillus, Chimaereicella, Pseudomonas, Microbacterium, Shewanella, Mycobacterium, and Rhodococcus have been reported with the ability to accumulate B from soil (Ahmed et al., 2007; Raja and Omine, 2013). There are reports on increased rhizosphere microbial populations by B in soybean (Sun et al., 2013). There were reports of B improving the population of diverse bacterial orders (Burkholderiales, Nitrosospherales, and Rhodospirales) (Vera et al., 2019). Boron aids in the enhancement of endomycorrhizae in roots owing to the action of indole 3-acetic acid (IAA) oxidase activity that activates IAA intensities eventually augmenting the translocation of carbohydrates to roots thereby improving fungi–mycorrhizal interaction and its fungal population (Lambert et al., 1980; Kumar et al., 2017). Related findings were reported by Bilen et al. (2011), where the highest population of bacteria and fungi production were observed with 2 kg B ha–1 in altered growth periods of the plant and diverse soil depths.
Even though the direct role of the effect of B on soil enzymes could not be established, however, the improvement in the status of enzyme activities (AS, DH, FDA, and PMA) was observed during the two years of experimentation. Possibly it was speculated that the enhanced enzyme activities could be attributed to improved soil condition due to the continuous addition of organic matter in the form of FYM in all the crops for two years which enhanced greater microbial activities. Since mineralised C and N from FYM enhanced the soil physico-chemical properties and the quantum of applied-B, the microbial populace and soil enzyme activities increased (Bhardwaj and Datt, 1995; Kumar et al., 2022). The improvement in the soil enzymatic activities might be attributed to the readily degradable organic matter added to the soil, which increases soil microorganisms and soil enzyme activities (Perucci, 1992). Also, DH enzyme activities in soil improved under graded B-fertilization. A similar finding was reported by Bilen et al. (2011) who observed a significant (p < 0.01) positive correlation with B and DH enzyme activity. DHA is a key indicator of microbial activity and organic matter stability since it is directly implicated in microbial respiration (Nikaeen et al., 2015). Improvement in the rhizosphere soil enzyme activity of soybean on the appliance of B was also reported (Goldberg, 1997; Liu, 2000). It was reported that surface residues encouraged the conservation of mineralizable C via residues that bettered the activities of soil enzymes ensuing in higher soil microbial biomass carbon and enhanced soil quality (Mohammadi et al., 2012). The increment in the soil microbial populations enhanced the rhizosphere metabolisms and bettered the soil enzyme activities on the appliance of B (Sun et al., 2013). Urease (urea amidohydrolase, EC 3.5.1.5) is N-related extracellular enzyme, the enzyme implicated in the degradation of urea cleaving urea to NH3 and CO2 (Kappaun et al., 2018). Evidently, UE enzyme activity de-escalation in all the crops at different growth stages was observed in B applied to soil as compared to control. The most likely reason may be due to the fact that B containing acids acts as a UE inhibitor in soil (Vera et al., 2019). It was reported that UE activity was inhibited and minimized by the application of higher B levels. Furthermore, this profound effect of B on the UE enzyme activity might be attributed to its action on water-soluble N, apart from the structural impairment to the enzyme and improved availability of N, as amino-boranes, which might have inhibited the urease activity via feedback mechanism (Kappaun et al., 2018). It was reported that the enzyme UE was immobilized on a membrane of microbes as inhibited B (Zaborska, 1995).
Soil functional diversity is governed by substrate amount, quality, and microbial accessibility. (Bending et al., 2002). Consequently, the plots which received higher C sources could alter microbial load of organic matter and the functional diversity of the microbial community in soils (Sall et al., 2006). Higher H values observed under the plots receiving escalated 2 kg B ha–1 might be due to the higher receipt of B dose which in turn improves SOC contents in the soil as there existed a positive correlation between B and SOC had been earlier documented (Bilen et al., 2011). Similarly, in the current investigation also, higher SYI values hold true in those plots where the highest B levels were applied and this is attributed to the greatest C availability by catabolic diversity (Lagomarsino et al., 2011) ensuing in better soil functional diversity.
The forming of clusters (Figure 4) with MBC, MBN, and SR in the soil in CCOCS may be possibly due to the fact that MBC is on the whole readily decomposable pool of organic material owing to the simplest structure and high quality of C and nutrients where soil bacteria used to nourish (Singh et al., 1992). The primary activation of microbial activity possibly ensued from swift catabolism of simple soluble C compounds, thereby, augmenting microbial biomass load in soil (Singh and Singh, 1991). Due to the possibility of similarity in their function, the soil enzymes were found to accumulate in similar clusters in the soil.
Normally, the first component explains most of the variability contained in the data set (Johnson and Wichern, 2002). The loading plot (Figure 5) exhibited that BP, MBN, PMC and MBC, PMN, actinomycetes, soil enzymes barring urease and subsequently followed by FP, and are highly correlated to each other in CCOCS. Because these parameters are highly correlated to each other as the angle within the variables of 0 or 180° reveals a correlation of 1 or −1, respectively (Kohler and Luniak, 2005). However, UE activity was negatively correlated with PC1 which is ascribed to the inverse link between levels (Vera et al., 2019). Results of the path analysis (Figure 6) also showed that barring urease enzymes, all twelve biological parameters indirectly contributed to MYCS by largely linking to the BP and their indirect path coefficients through BP. The contributions of the biological entities to mean yield of the cropping sequence followed the trend of BP > MBN > PMC > AP > PMC > PMN > AP > BP > FP > AS > DH > FDA > PMA > UE, respectively. The execution of SMLR regression models is considered best fitting (R2 > 0.9), as, R2 is in the range of 0.90 and 1 (Ogwueleka and Ogwueleka, 2010). In general, an R2 value ranging between 0.8 and 0.9 implies a good fitting and values lesser than 0.8 signify a poor model.
Conclusion
The results of the foregoing study revealed a tangible and significant impact of the graded levels of B-fertilization on soil biological entities under cauliflower-cowpea-okra cropping system in an acid Inceptisol. The key biological properties like BP, MBC, MBP, PMC, microbial population, and soil respiration were enhanced significantly with the incremental dose of B-fertilization. The order of the B treatments in respect of MBC, MBN and soil respiration at different crop growth stages was 2.0 kg B ha–1> 1.5 B kg ha–1> 1.0 B kg ha–1> 0.5 B kg ha–1> 0 kg B ha–1, respectively. Higher recoveries of PMC and PMN were noticed under 2 kg B ha–1 in cauliflower, cowpea and okra, respectively, at all the crop growth stages over control. Barring urease, the activities of all other important soil enzymes (AS, DH, FDA, and PMA) were increased significantly up to the application of 2 kg B ha–1. The positive impact of B-fertilization on these biological properties was observed at different growth stages of all three crops in the sequence which ultimately led to higher and sustainable crop production. A significant and positive relationship between these properties and crop yield greatly supported this observation. Multivariate analysis also confirmed the role of B-fertilization in the augmentation of the soil’s biological properties and yield enhancement. Overall, it was concluded that different soil physico-chemical and biological properties under the cauliflower-cowpea-okra cropping sequence can be invariably improved by the application of graded levels of B up to 2 kg B ha–1 in an acid Inceptisol. Future research entails more advanced research between B with soil microbial pools, microbial populations and soil enzymes to explore the precise mechanism of their interaction in soil. Comprehension of the mechanisms underlying established functions of B may explicate the significance of B and, in the end, lead to an advanced perception of its biological function, which has vital pragmatic implications in agriculture. The continuity of the residual impact as well as the beneficial effect of B-fertilization in such cropping sequences is therefore an important subject for future empirical research to elucidate its sustainability.
Data availability statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author contributions
IB: conceptualization, methodology formulation and implementation, data analysis, project administration, resource, and review. AB and DS: analysis and graphical works, editing, and results validation. AKC: data curation, reviewing, and editing. AdK and SC: editing of original and revised versions of the manuscript. AS, ED, GG, and LO: data collection and processing and original draft preparation. SB and MDS: review of literature and results compilation. MS and BAG: review of literature and basic analysis. AmK and SHD: results compilation and draft preparation. BG and MH: graphs and maps preparation. YD, KS, and ST: review of literature, results compilation, and English language correction. MR: review and editing of the final draft of manuscript, validation of statistical analysis, and formatting the manuscript. All authors contributed to the article and approved the submitted version.
Acknowledgments
We are thankful to the Department of Soil Science, Assam Agricultural University, Jorhat, Assam, India for supporting in carrying out the analysis of the present research work.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.996220/full#supplementary-material
References
Adam, G., and Duncan, H. (2001). Development of a sensitive and rapid method for the measurement of total microbial activity using fluorescein diacetate (FDA) in a range of soils. Soil Biol. Biochem. 33, 943–951. doi: 10.1016/S0038-0717(00)00244-3
Ahmad, W., Niaz, A., Kanwal, S., and Khalid Rasheed, M. (2009). Role of boron in plant growth: A review. J. Agric. Res. 47, 329–338.
Ahmed, I., Yokota, A., and Fujiwara, T. (2007). A novel highly boron tolerant bacterium, Bacillus boroniphilus sp. nov., isolated from soil, that requires boron for its growth. Extremophiles 11, 217–224. doi: 10.1007/s00792-006-0027-0
Anderson, J. P. E., and Domsch, K. H. (1989). Ratios of microbial biomass carbon to total organiccarbon in arable soils. Soil Biol. Biochem. 21, 471–479. doi: 10.1016/0038-0717(89)90117-X
Andrews, S. S., Karlen, D. L., and Mitchell, J. P. (2002). A comparison of soil quality indexing methods for vegetable production systems in Northern California. Agric. Ecosyst. Environ. 90, 25–45. doi: 10.1016/S0167-8809(01)00174-8
Andrewsi, S. S., and Carroll, C. R. (2001). Designing a Soil quality assessment tool for sustainable agroecosystem management. Ecol. Soc. Am. 11, 1573–1585. doi: 10.1890/1051-0761(2001)011[1573:DASQAT]2.0.CO;2
Aneja, K. R. (2003). Experiments in microbiology, plant pathology and biotechnology, 4th Edn. New Delhi: New age international publication.
Awale, R., Emeson, M. A., and Machado, S. (2017). Soil organic carbon pools as early indicators for soil organic matter stock changes under different tillage practices in Inland Pacific Northwest. Front. Ecol. Evol. 5:96. doi: 10.3389/fevo.2017.00096
Balota, E. L., Colozzi-Filho, A., Andrade, D. S., and Dick, R. P. (2003). Microbial biomass in soils under different tillage and crop rotation systems. Biol. Fertil. Soils 38, 15–20. doi: 10.1007/s00374-003-0590-9
Bandick, A. K., and Dick, R. P. (1999). Field management effects on soil enzyme activities. Soil Biol. Biochem. 31, 1471–1479. doi: 10.1016/S0038-0717(99)00051-6
Behera, B., Kancheti, M., Raza, M. B., Shiv, A., Vikas, M., Rathod, G., et al. (2022). Mechanistic insight on boron-mediated toxicity in plant vis-a-vis its mitigation strategies: A review. Int. J. Phytoremediation 17, 1–18. doi: 10.1080/15226514.2022.2049694
Bellaloui, N., Mengistu, A., and Kassem, M. A. (2014). “Role of boron nutrient in nodules growth and nitrogen fixation in soybean genotypes under water stress conditions,” in Advances in biology and ecology of nitrogen fixation, ed. T. Ohyama (London: IntechOpen Limited), 237–257. doi: 10.5772/56994
Bending, G. D., Turner, M. K., and Jones, J. E. (2002). Interactions between crop residue and soil organic matter quality and the functional diversity of soil microbial communities. Soil Biol. Biochem. 34, 1073–1082. doi: 10.1016/S0038-0717(02)00040-8
Berger, K. C., and Pratt, P. J. (1963). “Advances in secondary and micro-nutrient fertilization,” in Fertilizer technology and use, eds M. McVickar, G. Bridger, and L. Nelson (Madison, WI: Soil Science Society of America), 281–340.
Bhardwaj, K. K. R., and Datt, N. (1995). Effects of legume green-manuring on nitrogen mineralization and some microbiological properties in an acid rice soil. Biol. Fertil. Soils 19, 19–21. doi: 10.1007/BF00336341
Bhupenchandra, I., Basumatary, A., Dutta, S., Singh, A. H., Singh, L. K., Bora, S. S., et al. (2021b). Effect of boron fertilization on soil chemical properties, nutrients status in the soil and yield of crops under cauliflower-cowpea okra sequence in North East India. Commun. Soil Sci. Plant Anal. 52, 1301–1326. doi: 10.1080/00103624.2021.1879127
Bhupenchandra, I., Basumatary, A., Dutta, S., Singh, L. N., Das, A., Singh, L. K., et al. (2021a). Direct and residual impact of boron fertilization improves the crop yield, nutrient contents, nutrient uptake, and nutrient use efficiencies in cauliflower–cowpea–okra sequence in an acidic Inceptisol of north-east India. J. Plant Nutr. 45, 1–21. doi: 10.1080/01904167.2021.1994591
Bhupenchandra, I., Chongtham, S. K., Basumatary, A., Singh, A. H., Das, A., Choudhary, A. K., et al. (2022). Changes in soil properties, productivity and profitability as influenced by the adoption of site-specific integrated crop management technology in turmeric (Curcuma longa L.) in Eastern Himalayan acidic Inceptisol. Ind. Crops Prod. 180, 1–15. doi: 10.1016/j.indcrop.2022.114745
Bilen, S., Bilen, M., and Bardhan, S. (2011). The effects of boron management on soil microbial population and enzyme activities. Afr. J. Biotechnol. 10, 5311–5319.
Bolaños, L., Brewin, N. J., and Bonilla, I. (1996). Effects of boron on Rhizobium-legume cell-surface interactions and nodule development. Plant Physiol. 110, 1249–1256. doi: 10.1104/pp.110.4.1249
Bolaños, L., Cebrián, A., and Redondo-Nieto, M. (2001). Lectin-like glycoprotein PsNLEC-1 is not correctly glycosylated and targeted in boron-deficient pea nodules. Mol. Plant Microbe Interact. 14, 663–670. doi: 10.1094/MPMI.2001.14.5.663
Bolaños, L., Lukaszewski, K., Bonilla, I., and Blevins, D. (2004a). Why boron? Plant Physiol. Biochem. 42, 907–912. doi: 10.1016/j.plaphy.2004.11.002
Bolaños, L., Redondo-Nieto, M., and Rivilla, R. (2004b). Cell surface interactions of Rhizobium bacteroids and other bacterial strains with symbiosomal and peribacteroid membrane components from pea nodules. Mol. Plant Microbe Interact. 17, 216–223. doi: 10.1094/MPMI.2004.17.2.216
Bollag, J.-M., Smith, J. L., and Paul, E. A. (2017). “The significance of soil microbial biomass estimations,” in Soil biochemistry, 1st Edn, ed. J.-M. Bollag (New York, NY: Routledge), 357–398. doi: 10.1201/9780203739389-7
Bonde, T. A., Schnürer, J., and Rosswall, T. (1988). Microbial biomass as a fraction of potentially mineralizable nitrogen in soils from long-term field experiments. Soil Biol. Biochem. 20, 447–452. doi: 10.1016/0038-0717(88)90056-9
Bonilla, I., Mergold-Villaseñor, C., and Campos, M. E. (1997). The aberrant cell walls of boron-deficient bean root nodules have no covalently bound hydroxyproline-/proline-rich proteins. Plant Physiol. 115, 1329–1340. doi: 10.1104/pp.115.4.1329
Bowerman, B., O‘Connell, R. T., and Koehler, A. B. (2005). Forecasting, Time series, and regression: An applied approach, 4th Edn. Belmont, CA: Thomson Brooks/Cole.
Brdar-Jokanović, M. (2020). Boron toxicity and deficiency in agricultural plants. Int. J. Mol. Sci. 21:1424. doi: 10.3390/ijms21041424
Campbell, C. A., Ellert, B. H., and Jame, Y. W. (1993). “Nitrogen mineralization potential in soils,” in Soil sampling and methods of analysis, ed. M. R. Carter (Boca Raton, FL: Lewis), 341–349.
Campbell, C., and Curtin, D. (2007). “Mineralizable nitrogen,” in Soil sampling and methods of analysis, 2nd Edn, ed. Taylor Francis (Boca Raton, FL: CRC Press), 206. doi: 10.1201/9781420005271.ch46
Casida, L. E., Klein, D. A., and Santoro, T. (1964). Soil dehydrogenase activity. Soil Sci. 98, 371–376. doi: 10.1097/00010694-196412000-00004
Cattell, R. B. (1966). The scree test for the number of factors. Multivariate Behav. Res. 1, 245–276. doi: 10.1207/s15327906mbr0102_10
Chen, T. H., Chiu, C. Y., and Tian, G. L. (2005). Seasonal dynamics of soil microbial biomass incoastal sand dune forest. Pedobiologia 49, 645–653. doi: 10.1016/j.pedobi.2005.06.005
Choudhary, A. K., Rana, K. S., Dass, A., and Srivastav, M. (2014). Advances in vegetable agronomy. New Delhi: Post Graduate School, IARI.
Deluca, T. H., Pingree, M. R. A., and Gao, S. (2019). Assessing soil biological health in forest soils. Dev. Soil Sci. 36, 397–426. doi: 10.1016/B978-0-444-63998-1.00016-1
Dénarié, J., and Cullimore, J. (1993). Lipo-oligosaccharide nodulation factors: A new class of signaling molecules mediating recognition and morphogenesis. Cell 74, 951–954. doi: 10.1016/0092-8674(93)90717-5
Dick, R. P. (1994). “Soil enzyme activities as indicators of soil quality,” in Defining soil quality for a sustainable environment, eds J. Doran, D. Coleman, D. Bezdicek, and B. Stwart (Madison, WI: Soil Science Society of America, Special publication), 107–124. doi: 10.2136/sssaspecpub35.c7
Friedel, J. K., Munch, J. C., and Fischer, W. R. (1996). Soil microbial properties and the assessment of available soil organic matter in a haplic luvisol after several years of different cultivation and crop rotation. Soil Biol. Biochem. 28, 479–488. doi: 10.1016/0038-0717(95)00188-3
Furtak, K., and Gajda, A. M. (2018). “Activity and variety of soil microorganisms depending on the diversity of the soil tillage system,” in Sustainability of agroecosystems, ed. IntechOpen (London: Intech Open), 45–61. doi: 10.5772/intechopen.72966
Ghosh, S., Wilson, B. R., Ghoshal, S. K., Senapati, N., and Mandal, B. (2012). Organic amendmentsinfluence soil quality and carbon sequestration in the Indo-Gangetic plains of India. Agric. Ecosys. Environ. 156, 134–141. doi: 10.1016/j.agee.2012.05.009
Goldbach, H., and Monika, W. (2007). Boron in plants and animals: Is there a role beyond cell-wall structure? J. Plant Nutr. Soil Sci. 170, 39–48. doi: 10.1002/jpln.200625161
Goldberg, S. (1997). Reactions of boron with soils. Plant Soil 193, 35–48. doi: 10.1007/978-94-011-5580-9_3
Goswami, G., Deka, P., and Das, P. (2017). Diversity and functional properties of acid-tolerant bacteria isolated from tea plantation soil of Assam. 3 Biotech 7:229. doi: 10.1007/s13205-017-0864-9
Grattan, S. R., Díaz, F. J., Pedrero, F., and Vivaldi, G. A. (2015). Assessing the suitability of saline wastewaters for irrigation of Citrus spp.: Emphasis on boron and specific-ion interactions. Agric. Water Manage. 157, 48–58. doi: 10.1016/j.agwat.2015.01.002
GreenFacts (2002). Levels 2 - details on boron. Available online at: https://www.greenfacts.org/en/boron/boron-greenfacts-level2.pdf (assessed April 4, 2022).
Guo, Z., Han, J., and Xu, Y. (2019). The mineralization characteristics of organic carbon and particle composition analysis in reconstructed soil with different proportions of soft rock and sand. PeerJ 7:e7707. doi: 10.7717/peerj.7707
Gupta, U. C. (1993). “Factors affecting boron uptake by plants,” in Boron and its role in crop production, ed. U. C. Gupta (Boca Raton, FL: CRC Press Inc), 87–104.
Harish, M. N., Choudhary, A. K., Kumar, S., Dass, A., Singh, V. K., Sharma, V. K., et al. (2022). Double zero-tillage and foliar phosphorus fertilization coupled with microbial inoculants enhance maize productivity and quality in a maize–wheat rotation. Sci. Rep. 12, 1–23. doi: 10.1038/s41598-022-07148-w
Hinojosa, M. B., García-Ruíz, R., Viñegla, B., and Carreira, J. A. (2004). Microbiological Rates and enzyme activities as indicators of functionality in soils affected by the aznalcóllar toxic spill. Soil Biol. Biochem. 36, 1637–1644. doi: 10.1016/j.soilbio.2004.07.006
Jenkinson, D. S., and Powlson, D. S. (1976). The effects of biocidal treatments on metabolism in soil-V. A method for measuring soil biomass. Soil Biol. Biochem. 8, 209–213. doi: 10.1016/0038-0717(76)90005-5
Jjemba, P. K., and Alexander, M. (1999). Possible determinants of rhizosphere competence of bacteria. Soil Biol. Biochem. 31, 623–632. doi: 10.1016/S0038-0717(98)00168-0
Johnson, R. A., and Wichern, D. W. (2002). Applied multivariate statistical analysis. New Jersy, NJ: Prentice-Hall.
Kaiser, H. F. (1960). The application of electronic computers to factor analysis. Educ. Psychol. Meas. 20, 141–151. doi: 10.1177/001316446002000116
Kappaun, K., Piovesan, A. R., Carlini, C. R., and Ligabue-Braun, R. (2018). Ureases: Historical aspects, catalytic, and non-catalytic properties – A review. J. Adv. Res. 13, 3–17. doi: 10.1016/j.jare.2018.05.010
Kohler, U., and Luniak, M. (2005). Data inspection using biplots. Stata J. 5, 208–223. doi: 10.1177/1536867X0500500206
Kumar, A., Choudhary, A. K., and Suri, V. K. (2017). Agronomic bio–fortification and quality enhancement in okra–pea cropping system through arbuscular mycorrhizal fungi at varying phosphorus and irrigation regimes in Himalayan acid Alfisol. J. Plant Nutr. 40, 1213–1229. doi: 10.1080/01904167.2016.1267208
Kumar, A., Choudhary, A. K., Pooniya, V., Suri, V. K., and Singh, U. (2016). “Soil factors associated with micronutrient acquisition in crops- biofortification perspective,” in Biofortification of food crops, eds U. Singh, C. S. Praharaj, S. S. Singh, and N. P. Singh (New York, NY: Springer Publishers), 159–176. doi: 10.1007/978-81-322-2716-8_13
Kumar, A., Rana, K. S., Choudhary, A. K., Bana, R. S., Sharma, V. K., Gupta, G., et al. (2022). Sole- or dual-crop basis residue-mulching and Zn-fertilization lead to improved productivity, rhizo-modulation and soil health in zero-tilled pigeonpea–wheat cropping system. J. Soil Sci. Plant Nutr. 22, 1193–1214. doi: 10.1007/s42729-021-00723-6
Kumar, M., and Sen, N. L. (2004). Interaction effect of zinc and boron on okra (Abelmoschus esculentus L. Moench) cv. Prabhani kranti. Agnc. Sci. Digest 24, 307–308.
Lagomarsino, A., Benedetti, A., Marinari, S., Pompili, L., Moscatelli, M. C., Roggero, P. P., et al. (2011). Soil organic C variability and microbial functions in a Mediterranean agro-forest ecosystem. Biol. Fertil. Soils 47, 283–291. doi: 10.1007/s00374-010-0530-4
Lambert, D. H., Cole, H., and Baker, D. E. (1980). The role of boron in plant response to mycorrhizal infection. Plant Soil 57, 431–438. doi: 10.1007/BF02211700
Liu, P. (2000). The effect of molybdenum and boron on nutritional and physiological mechanism of yield and quality in soybean. Hangzhou: Zhejiang University.
Mohammadi, K., Heidari, G., and Karimi, M. T. (2012). Contrasting soil microbial responses to fertilization and tillage systems in canola rhizosphere. Saudi J. Biol. Sci. 19, 377–383. doi: 10.1016/j.sjbs.2012.05.001
Mooshammer, M., Wanek, W., and Hämmerle, I. (2014). Adjustment of microbial nitrogen use efficiency to carbon: Nitrogen imbalances regulates soil nitrogen cycling. Nat. Commun. 5, 1–7. doi: 10.1038/ncomms4694
Nadeem, F., Farooq, M., Nawaz, A., and Ahmad, R. (2019). Boron improves productivity and profitability of bread wheat under zero and plough tillage on alkaline calcareous soil. Field Crop Res. 239, 1–9. doi: 10.1016/j.fcr.2019.05.010
Nason, A., and McElory, W. D. (1963). “Modes of action of essential mineral elements,” in Plant physiology, ed. F. C. Steward (New York, NY: Academic Press), 451–536. doi: 10.1016/B978-0-12-395600-2.50017-6
Nikaeen, M., Nafez, A. H., Bina, B., Nabavi, B. F., and Hassanzadeh, A. (2015). Respiration and enzymatic activities as indicators of stabilization of sewage sludge composting. Waste Manag. 39, 104–110. doi: 10.1016/j.wasman.2015.01.028
Ogwueleka, T. C., and Ogwueleka, F. N. (2010). Modelling energy content of municipal solid waste using artificial neural network. Iran. J. Environ. Health Sci. Eng. 7, 259–266.
Panse, V. G., and Sukhatme, P. V. (1985). Statistical methods for agricultural workers. New Delhi: Indian Council of Agricultural Research Publication.
Pavan, F. S. A., Bettiol, W., and Cerri, C. C. (2005). Effect of sewage sludge on microbial biomass, basal respiration, metabolic quotient and soil enzymatic activity. Appl. Soil Ecol. 30, 65–77. doi: 10.1016/j.apsoil.2004.03.008
Paz-Ferreiro, J., Gascó, G., Gutiérrez, B., and Méndez, A. (2012). Soil biochemical activities and the geometric mean of enzyme activities after application of sewage sludge and sewage sludge biochar to soil. Biol. Fertil. Soils 48, 511–517. doi: 10.1007/s00374-011-0644-3
Perucci, P. (1992). Enzyme activity and microbial biomass in a field soil amended with municipal refuse. Biol. Fertil. Soils 14, 54–60. doi: 10.1007/BF00336303
Rahman, M., Hossain, I., and Ahmad, M. (2017). Effects of boron and zinc on yield and quality of okra seed. Adv. Biores. 8, 202–211.
Raja, C. E., and Omine, K. (2013). Characterization of boron tolerant bacteria isolated from a fly ash dumping site for bacterial boron remediation. Environ. Geochem. Health 35, 431–438. doi: 10.1007/s10653-012-9505-8
Rajpoot, S. K., Rana, D. S., and Choudhary, A. K. (2021). Crop and water productivity, energy auditing, carbon footprints and soil health indicators of Bt-cotton transplanting led system intensification. J. Environ. Manag. 300:113732. doi: 10.1016/j.jenvman.2021.113732
Rana, K. S., Choudhary, A. K., Sepat, S., Bana, R. S., and Dass, A. (2014). Methodological and analytical agronomy. New Delhi: Post Graduate School, IARI.
Redondo-Nieto, M., Rivilla, R., and El-Hamdaoui, A. (2001). Boron deficiency affects early infection events in the pea-Rhizobium symbiotic interaction. Aust. J. Plant Physiol. 28, 819–823. doi: 10.1071/PP01020
Saha, P., Chatterjee, R., and Das, N. R. (2010). Effect of foliar application of boron and molybdenum in sprouting broccoli (Brassica oleracea var. italica P.) under trai region of West Bengal. Res. J. Agric. Sci. 1, 335–337.
Sall, S. N., Masse, D., Ndour, N. Y. B., and Chotte, J. L. (2006). Does cropping modify the decomposition function and the diversity of the soil microbial community of tropical fallow soil? Appl. Soil Ecol. 31, 211–219. doi: 10.1016/j.apsoil.2005.05.007
Shelp, B. (1993). “Physiology and biochemistry of boron in plants,” in Boron and its Role in crop protection, ed. U. Gupta (Boca Raton, FL: CRC Press), 53–85.
Shireen, F., Nawaz, M., and Chen, C. (2018). Boron: Functions and Approaches to enhance its availability in plants for sustainable agriculture. Int. J. Mol. Sci. 19:1856. doi: 10.3390/ijms19071856
Shorrocks, V. M. (1997). The occurrence and correction of boron deficiency. Plant Soil 193, 121–148. doi: 10.1007/978-94-011-5580-9_9
Singh, B., Singh, Y., Sadana, U. S., and Meelu, O. P. (1992). Effect of organic amendments and moisture regimes on the kinetics of micronutrients in a calcareous sandy loam soil. J. Indian Soc. Soil Sci. 40, 114–118.
Singh, M., and Singh, T. B. (1991). Influence of organic and inorganic amendments, modified urea, and application methods on ammonia volatilization in saturated calcareous soil. Int. Rice Res. Newsl. 16, 17–18.
Singh, U., Choudhary, A. K., and Sharma, S. (2020). Comparative performance of conservation agriculture vis-a-vis organic and conventional farming in enhancing plant attributes and rhizospheric bacterial diversity in Cajanus cajan: A field study. Eur. J. Soil Biol. 99:103197. doi: 10.1016/j.ejsobi.2020.103197
Singh, U., Choudhary, A. K., and Sharma, S. (2021). Agricultural practices modulate the bacterial communities, and nitrogen cycling bacterial guild in rhizosphere: Field experiment with soybean. J. Sci. Food Agric. 101, 2687–2695. doi: 10.1002/jsfa.10893
Sinsabaugh, R. L., and Follstad Shah, J. J. (2012). Ecoenzymatic Stoichiometry and ecological theory. Annu. Rev. Ecol. Evol. Syst. 43, 313–343. doi: 10.1146/annurev-ecolsys-071112-124414
Spaink, H. P. (2000). Root nodulation and infection factors produced by Rhizobial bacteria. Annu. Rev. Microbiol. 54, 257–288. doi: 10.1146/annurev.micro.54.1.257
Sun, T., Wang, Y. P., and Wang, Z. Y. (2013). The effects of molybdenum and boron on the rhizosphere microorganisms and soil enzyme activities of soybean. Acta Physiol. Plant 35, 763–770. doi: 10.1007/s11738-012-1116-6
Tabatabai, M. A. (1994). “Soil enzymes,” in Methods of soil analysis: Part 2—Microbiological and biochemical properties, eds P. S. Bottomley, J. S. Angle, and R. W. Weaver (Madison, WI: Soil Science Society of America), 775–833.
Tabatabai, M. A., and Bremner, J. M. (1970). Arylsulfatase activity of soils. Soil Sci. Soc. Am. Proc. 34, 225–229. doi: 10.2136/sssaj1970.03615995003400020016x
Takkar, P. N., and Randhawa, N. (1978). Micronutrients in Indian agriculture. Fertil. News 23, 3–26.
Vera, A., Moreno, J. L., and García, C. (2019). Boron in soil: The impacts on the biomass, composition and activity of the soil microbial community. Sci. Total Environ. 685, 564–573. doi: 10.1016/j.scitotenv.2019.05.375
Verma, C. K., Prasad, K., and Yadav, D. D. (2012). Studies on the response of sulphur, zinc and boron levels on yield, economics and nutrients uptake of mustard [Brassica juncea (L.) Czern & Coss.]. Crop Res. 44, 75–78.
Walkley, A., and Black, I. A. (1934). An examination of the Degtjareff method for determiningsoil organic matter and a proposed modification of the chromic acid titration method. Soil Sci. 37, 29–38. doi: 10.1097/00010694-193401000-00003
Wold, S., Esbensen, K., and Geladi, P. (1987). Principal component analysis. Chemom. Intell. Lab. Syst. 2, 37–52. doi: 10.1016/0169-7439(87)80084-9
Xu, J. M., Tang, C., and Chen, Z. L. (2006). The role of plant residues in pH change of acid soils differing in initial pH. Soil Biol. Biochem. 38, 709–719. doi: 10.1016/j.soilbio.2005.06.022
Zaborska, W. (1995). Competitive inhibitors of free and chitosan-immobilized urease. Acta Biochim. Pol. 42, 115–118. doi: 10.18388/abp.1995_4677
Keywords: boron, Inceptisol, microbiological pools, microbial populations, microbial biomass carbon, soil enzymatic activities
Citation: Bhupenchandra I, Basumatary A, Choudhary AK, Kumar A, Sarkar D, Chongtham SK, Singh AH, Devi EL, Bora SS, Salam MD, Sahoo MR, Gudade BA, Kumar A, Devi SH, Gogoi B, Harish MN, Gupta G, Olivia LC, Devi YP, Sarika K, Thapa S and Rajawat MVS (2022) Elucidating the impact of boron fertilization on soil physico-chemical and biological entities under cauliflower-cowpea-okra cropping system in an Eastern Himalayan acidic Inceptisol. Front. Microbiol. 13:996220. doi: 10.3389/fmicb.2022.996220
Received: 17 July 2022; Accepted: 30 August 2022;
Published: 07 November 2022.
Edited by:
Anukool Vaishnav, Agroscope, SwitzerlandReviewed by:
Vishal Tripathi, GLA University, IndiaBilal Ahmed, Yeungnam University, South Korea
Rajesh Kumar Singh, Guangxi Academy of Agricultural Science, China
Sarvesh Pratap Kashyap, Banaras Hindu University, India
Copyright © 2022 Bhupenchandra, Basumatary, Choudhary, Kumar, Sarkar, Chongtham, Singh, Devi, Bora, Salam, Sahoo, Gudade, Kumar, Devi, Gogoi, Harish, Gupta, Olivia, Devi, Sarika, Thapa and Rajawat. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Mahendra Vikram Singh Rajawat, rajawat.mvs@gmail.com; Anil K. Choudhary, anilpau2010@gmail.com