Exploring the Synergistic Impacts of Cover Crops and Fertilization on Soil Microbial Metabolic Diversity in Dryland Soybean Production Systems Using Biolog EcoPlates

: The metabolic diversity of soil microbiota embodies diverse functional capabilities that support ecosystem resilience, driving essential biogeochemical processes and facilitating the optimization of sustainable agricultural systems. Integrating cover crops into agricultural systems cultivates a diverse array of metabolic activities among soil microbes, synergistically enhancing ecosystem services and bolstering soil health for sustainable and productive farming practices. In an effort to gain deeper insights and expand our knowledge, we conducted a study examining the effects of cover crops and fertilizer sources, thereby shedding light on their combined impacts on the metabolic activity dynamics of soil microbial communities. In this investigation, we employed a split-plot design with two factors: (a) cover crop with three solo cover crop species—Cereal rye ( Secale cereale ), wheat ( Triticum aestivum ), hairy vetch ( Vicia villosa ) , and one mixture of mustard ( Brassica rapa ) and cereal rye ( Secale cereale ) (CC-mix), (b) Fertilizer source includes poultry litter, chemical fertilizer, and no-fertilizer treatments. We assessed the metabolic potential of soil microbiota by using carbon substrates utilizing Biolog EcoPlates. The ﬁndings revealed that the plots with CC-mix treatment exhibited greater metabolic diversity compared to the other treatments, while among the fertilizer sources, poultry litter demonstrated higher metabolic activity. Furthermore, both treatment factors predominantly metabolized carbohydrates and polymers compared to other carbon substrate categories. The principal component analysis accounted for 46.4% of the variance, collectively represented by PC1 and PC2, emphasizing the substantial contributions of carbohydrates, amino acids, and carboxylic acids to the observed metabolic diversity. Canonical correspondence analysis revealed that pH had positively correlated with microbial functional diversity, whereas total carbon (TC), total nitrogen (TN), and water-stable aggregates (WSA) showed a negative correlation. In conclusion, cover cropping and type of fertilizer source had a notable impact on soil microbial functional diversity, with the cover crop mixture exhibiting a more pronounced inﬂuence than the individual cover crop treatments.


Introduction
Soil is a highly diverse, intricate, and dynamic ecosphere that comprises numerous living and non-living elements persistently communicating with each other.Soil health can undergo changes over time as a result of both natural occurrences and human activities.Soil health is improved through effective management and land-use choices that consider the diverse functions of soil [1].Unsustainable soil management practices worsened soil erosion, a deterioration in soil microbiome activity, and a consequent reduction in nutrient availability for plant uptake [2].Nature-based solutions (NBS) strategies, such as reduced soil tillage, soil covering with crop residuals or cover crops, soil amendments, and crop rotation to achieve microbial diversity, have been suggested as a substitution for conventional agriculture for improving soil health characteristics [3].The association between soil biology and agricultural sustainability is paramount, as microbes play a pivotal role in facilitating the decomposition of intricate organic compounds into plant-available forms [4].
Over the past few years, our comprehension of soil microbial communities in dryland ecosystems has broadened, as these organisms play a vital role in maintaining soil health and ecosystem functions [5].For instance, the soil microbiome mediates various essential ecosystem services, such as nitrogen incorporation [6], organic carbon decomposition, carbon capture [7], bio-remediation of environmental pollutants [3], and resistance to biotic and abiotic stresses [8,9].Nonetheless, such ecosystem services are compromised by intensive farming systems and tillage practices [10] due to excessive synthetic fertilizers and monoculture cropping systems [10].This lowers soil microbial functional diversity, which significantly influences crop performance [11].Expanding diversified cropping systems promotes sustainable agriculture, potentially improving ecosystem functioning by enhancing soil microbial diversity as observed in natural ecosystems [12,13].The objective is to enhance the functional diversity of soil microbiota by promoting microbial diversity, thereby boosting its capacity to metabolize a broad range of organic substances.Introducing diversity in organic matter residues in agricultural systems, e.g., cover crop mixtures, is anticipated to enhance functional diversity [14,15].Integration of cover crops and organic amendments alters the abundance, diversity, and metabolic activities of microoganisms present in the agricultural soils, and these altered microbial groups have various effects on their soil ecosystem services [16], which require further exploration.
A large body of literature specifies that, different cover crops with varying crop rotations and periods have diverse impacts and roles on soil health properties and microorganisms [17,18].The legume cover crops are potent nitrogen fixers (i.e., symbiotic association with Rhizobium) and can transform atmospheric nitrogen into available forms (NH 3 and NH 4 ) [19].Legumes generate residue with a comparatively higher nitrogen content that is readily degradable and accessible to crops and soil microbes.Grass species with their root system have strong nutrient scavenging capability, specifically for N, and hence have been used to reduce N leaching [20].Brassicaceous species cover crops have a deep root system, can provide an allelopathic effect to control weeds, and have bio fumigation properties for the suppression of soil pathogens [21,22].Cover crops are characterized by their rapid growth rate and exhibit a range of C:N ratios, leading to the production of residues that possess a stoichiometry that differs from microbial biomass [23].Winter cover crops such as cereal rye and hairy vetch are widely utilized due to their ability to withstand low temperatures, incorporate substantial biomass [24], N fixing [25], and stimulate soil biological activity [26].Three-year research conducted by Buyer et al. (2010) reported that incorporating vetch and rye into crop rotations significantly enhanced soil microbial biomass and metabolic activity [27].In addition, cover crops, compost, and manure are recognized as vital factors in sustainable farming as they supply vital nutrients to plants, enhance soil structure and quality, and add organic matter [28].
Single-cover crop species cannot deliver all the benefits for soil health [29].Therefore, multifunctional cover crops, such as combining legumes and grasses or integrating multispecies mixture into crop production systems, can be more advantageous for improving soil health since they can expand the variety of substrates available to the ecosystem and provide a number of ecosystem services [29,30].The inclusion of plant mixtures with varying C:N ratios can provide diverse substrate options for the soil system, which might expand the niche breadth of the soil microbiota, thus enhancing its metabolic potential [31].For instance, a combination of oats, radish, and vetch showed significant enhancement in the abundance of gram-positive bacteria, and phosphatase function using phospholipid fatty acid analysis more effectively than a combination of oats and radish [32].Given that different cover crop species have chemically different root exudates, Housman et al. (2021) reported that higher microbial biomass and improved enzymatic activity in multiple cover crops compared to solo species also differed effects between species [33].Moreover, elevated soil microbial functional activity was observed in cover crop mixture plots compared to solo treatments [30].
The biological variation of microbial community composition is being studied from three perspectives, viz., species diversity, genetic diversity, and functional diversity [34,35].The functional diversity of soil microbiome is a component of biodiversity that comprises a broad range of metabolic behaviors that can affect many elements of ecosystem functions, such as ecological stability, nutrient cycling, ecosystem dynamics, and so on [36].Identifying the functional microbial diversity in soils entails selecting the most appropriate approaches for examining the maximum diversity present in the samples.
Cultivation-dependent approaches require the preparation of cell suspensions, dilution, inoculation on solid media, and phenotypic examination of microbial species.The majority of research on the structure of heterotrophic microbial communities has used isolate-based techniques [37].The isolation of microbial samples varies among species, and cell type (cell, spore, or mycelium) generates bias.Furthermore, designing nutrient media and appropriate conditions for all members of microbial communities presents a significant challenge, which are major disadvantages of these approaches [36].Microbial lipid-based techniques without actual culturing of the microorganisms, namely, PLFA and fatty acid methyl ester (FAME), have been extensively used for the examination of microbiome metabolic and compositional changes due to agricultural practices [38].These techniques, however, unquestionably have drawbacks, including complicated operational protocols, time-consuming assays, poor precision, and inconsistent repeatability [39,40].
Community-Level Physiological Profiling (CLPP) offers an alternative approach, comparing and categorizing microbial communities based on their utilization of carbon substrates.The Biolog ® EcoPlates TM (Biolog Inc., Hayward, CA, USA) has been specifically designed to examine the functional diversity of microbial communities.This plate has 31 separate carbon source wells, each with three replications, representing six different types of carbon compound families: carbohydrates, carboxylic acids, amino acids, polymers, amines, and phenolic compounds, with a blank well as a control [41].The principle behind this assay involves the reduction reaction of water-soluble colorless triphenyl tetrazolium chloride to purple triphenyl formazan [41].This technique is straightforward to implement, as it does not necessitate the use of isolated cultures and preserves the metabolic traits of microbial communities at optimal levels [42].
Soybean (Glycine max L.) is a legume oil seed crop known to fix atmospheric nitrogen through symbiotic interactions with soil bacteria [43].It has numerous applications, including the production of soy milk, soy sauces, tofu, edible oil, economically recoverable phytohormones, and biodiesel.Traditional soybean monoculture has been coupled with numerous ecological concerns, especially regarding soil health parameters, such as loss of organic matter [44] and available phosphorus [45].Soil conservation and sustainable agriculture approaches such as cover cropping, crop rotations, and organic fertilizers have widely been recognized as substitutes for soybean monocropping with positive influences on soil organic carbon, aggregate soil stability, nutrient cycling, N fixation, and high-water holding capacity [46][47][48].
The ecosystem services of different cover crops with fertilizer sources vary, which allows us to evaluate cover crop-mediated soil health promotion processes specific to microbial functional diversity and to correlate soil physicochemical characteristics such as pH, organic carbon, active carbon (POXC), total nitrogen, and glomalin in dryland soybean production systems.Overall, the objective of this experiment was to examine the effect of cover crops and different types of fertilizer sources (organic and inorganic) on microbial community composition, functional diversity, and soil health parameters.
In this experiment, we executed of two-year field study in which cover crop solo and mixture species were integrated into dryland soybean production systems to test the following hypotheses:

•
Soil microbial functional diversity and soil health characteristics were positively affected by cover crop mixtures and organic amendments in dryland soybean production systems.

•
Cover crop mixture leads to higher microbial functional diversity compared to solo cover crop treatments.

Description of the Experimental Sites
The study area was situated at the Pontotoc Ridge-Flatwoods Branch Experimental Station in Pontotoc, MS, USA (34 • 09 N, 88 • 58 W).The soil type at the research location belongs to the Atwood silt loam series, characterized as a semi-active, mixed soil with fine-silty texture, thermic Typic Paleudalf, situated on a moderate 3% slope.The soil here consists of 13.9% clay, 17.6% sand, and 68.5% silt (Soil Survey, Natural Resources Conservation Service (NRCS)).The investigation was carried out under rainfed and zero tillage conditions, with initial soil characterization at the experimental site revealing 1.57% organic matter and a pH of 6.67.Prior to this study, the experimental site had been planted with corn and soybeans for 2016 and 2017, respectively.Cover crops were established in 2016 and maintained until 2019, while soybeans were cultivated from 2017-2019.The study site has an average January temperature of 5.4 • C, an average July temperature of 23 • C, and an annual average rainfall of 1483 mm [49].

Experimental Design and Field Methods
The study was arranged in a split-plot arrangement with a randomized block design and three replications each.In this experiment, two factors-cover crops as a primary factor and fertilizer source as a secondary factor were taken into consideration.The experimental design comprised 15 whole plots of 167.2 m 2 each, with 45 sub-plots of 74.3 m 2 designated for the different cover crop and fertilizer source treatments.

Soil Sampling and Analysis of Physicochemical Properties
At each plot, soil samples were collected at rhizosphere region 1-15 cm depth after the termination of cover crops.Roughly 500 g of soil was collected into a Ziplock bag and preserved at a −20 • C freezer used for estimation of the metabolic diversity of microbial species.Soil physicochemical characters such as total C, active carbon (POXC), total nitrogen, glomalin (EE-GRSP), water stable aggregate (WSA), and soil pH were estimated from samples that were air-dried overnight, grounded, and passed through a 2.0 mm sieve.Detailed protocols for estimating these soil parameters were described in our previous studies [16,50].

Biolog TM EcoPlates
The community-level physiological profiling of soil microbiota was analyzed using Biolog TM EcoPlates.Each plate consists of 31 wells with different carbon sources and one control well.All these are replicated three times to control variation in the sample inoculum.The consumption of carbon substrates in wells by microbial communities results in the color development of their respective wells.This color development is quantified spectrophotometrically using a plate reader.Firstly, 10 g of soil was measured and placed in a 250 mL conical flask, then added 100 mL of sterilized NaCl solution (0.85%), later agitated at 150 rpm for 30 min, and then left undisturbed in a refrigerator for 30 min.The clear suspension was transferred to a petri dish, and with the help of a multi-channel pipette, samples were loaded into a Biolog EcoPlate.The inoculated sample plates were placed in the dark during incubation at 25 • C, and absorbance was measured using a BioTek 800 TS (Agilent TM , Santa Clara, CA, USA) microplate reader every 24 h for five days.

Soil Microbial Communities Catabolic Profiling 2.5.1. Quantification of Average Well Color Development in Biolog EcoPlates Wells
The metabolic process of soil microbial populations in individual wells of Biolog EcoPlate leads to the formation of formazan, thus change in color of the tetrazolium dye is induced [42].The absorbance values of individual wells at 590 nm wavelength represent metabolic activity and are expressed as average well color development (AWCD) [37].
Here, C i represents the absorbance value (OD) of each well containing a carbon source; R represents the absorbance value of a blank well (without carbon substrate); and n represents the total number of wells with carbon substrates.

Determination of Diversity Indices of Microbial Populations:
Average well-color development estimated from Biolog EcoPlates can be utilized to estimate microbial metabolic diversity indices proposed by Yan (2011) [51].Absorbance measurements and optical density (OD) at 96 h were employed to estimate the average well-color development and metabolic diversity parameters.This time point was selected as it exhibited an ideal range of absorbance values for the analysis.The following functional diversity parameters were calculated: (1) Shannon Diversity Index (H ): Here, P i represents the ratio between ( C i − R) and the total absorbance of the total plate wells.ln denotes the natural logarithm with respect to P i .The Shannon diversity index (H ) is used to measure the microbial communities' heterogeneousness based on the concept of uncertainty.Higher uncertainty refers to greater diversity present in communities, and it gives microbial species diversity and evenness within the community.
(2) Shannon Evenness Index (E): This method centers on the uniformity of absorbance values at each well for all utilized carbon substrates [41].In this equation, S represents the total absorbance value of 31 wells.The Shannon evenness index (E) identifies the evenness of microbial-type abundance in the communities.
(3) Inverse Simpson Diversity Index (1/D): The Simpson Diversity Index accounts for both the overall number of species and the proportional representation of each individual species [52].In this formula, n i denotes the relative absorbance value in each ith well; N is the sum of the absorbance values of 31 wells [23].The inverse Simpson diversity index (1/D) defines the species diversity by estimating the probability of one species encountering another [53].
(5) Substrate richness (SR): Substrate richness (R) refers to the count of carbon substrates that are utilized by soil microorganisms.It was estimated as the sum of oxidized wells, which had to be at least 0.5 after 96 h incubation [56].
(6) Statistical analysis: Diversity parameters were calculated using Microsoft Excel Version 16.74 (Microsoft Corporation 2023, Redmond, WA, USA).Analysis of variance (ANOVA) and LSD (least significant differences) for two factors was performed for the cover crop and fertilizer source factors using JMP ® , Version 17.0.SAS Institute Inc., Cary, NC, 1989-2021.The plots were generated using the 'ggplot2' package Wickham (2016) within the R programming language (R Core Team, 2023) [57] and Sigma Plot, Version 15 (Systat Software Inc., San Jose, CA, USA).
To assess the functional composition of the bacterial community, we conducted a Principal Component Analysis (PCA) on the corrected values of Biolog absorbance.The Biolog absorbance values were adjusted by subtracting the absorbance values from the control well of the Biolog EcoPlate.PCA is a widely used statistical technique that simplifies complex multidimensional data into a smaller set of interpretable variables called principal components [37].This analysis enabled us to examine whether the measured variables could effectively differentiate the soil management treatments based on their principal components.PERMANOVA is a multivariate statistical method that allows for the analysis of dissimilarities in microbial community composition based on categorical or continuous variables.It helps to determine whether there are significant differences in community structure between different treatments or groups [58].In addition to PCA, Canonical Correspondence Analysis (CCA) was applied to explore and visualize the relationships between microbial communities and environmental variables.CCA helped us identify the key environmental factors that influenced the composition of microbial communities, providing valuable insights into the ecological processes that shape the dynamics of the microbial community [59].
The PERMANOVA (Permutational Multivariate Analysis of Variance) test, CCA (Canonical Correspondence Analysis), and PCA (Principal Component Analysis) were performed using a vegan package in R programming language [60].

Differences in AWCD over Time in Soils with Cover Crops and Fertilizer Source Treatment
Generally, the amount of carbon substrate utilization is directly proportional to the metabolic ability of the respective microbial communities, as determined by AWCD [37].Our results showed that AWCD values increased over time, indicating that soil microbiome has a greater ability to utilize carbon substrates when there is a longer interaction with the substrate (Figure 1A,B).The results revealed that the AWCD values from 0 to 24 h for all treatments were small, indicating a noticeable lag phase in the first 24 h.After 24 h, the metabolic activity of microbial populations rapidly increased for all treatments, indicat-ing soil microbes were at their growth phase and capable of utilizing carbon substrates in plates.
Generally, the amount of carbon substrate utilization is directly proportional to the metabolic ability of the respective microbial communities, as determined by AWCD [37].Our results showed that AWCD values increased over time, indicating that soil microbiome has a greater ability to utilize carbon substrates when there is a longer interaction with the substrate (Figure 1A,B).The results revealed that the AWCD values from 0 to 24 h for all treatments were small, indicating a noticeable lag phase in the first 24 h.After 24 h, the metabolic activity of microbial populations rapidly increased for all treatments, indicating soil microbes were at their growth phase and capable of utilizing carbon substrates in plates.Cereal rye treatment consistently exhibited the lowest AWCD at all incubation periods.Native vegetation treatment had significantly higher AWCD after 24, 48h, and 72 h than other cover crop treatments followed by CC-mix (p < 0.05) (Figure 1A).After 96 h of incubation, both native vegetation and CC-mix had similar AWCD, which was higher than the remaining cover crop treatments.In regard to fertilizer source treatments, AWCD was highest with PL (Poultry litter) treatment, followed by chemical fertilizer.The fertilizer treatment exhibited the lowest AWCD after a 96-h incubation period (Figure 1B).During the entire incubation period, the AWCD values for each cover crop followed the following pattern NV > CC-mix > VT > WT > CR.These results suggest that the metabolic activity of soil microbial communities in utilizing carbon substrates differs among cover crops and fertilizer sources.

Influence of Cover Cropping and Fertilizer Source Treatments on Microbial Metabolic Diversity Indices
Absorbance measured as OD values at 96 h were used for the estimation of functional diversity parameters and average well-color development since these indicated the optimal range of absorbance.The influence of cover crops and fertilization on functional diversity indices such as the Shannon diversity index (H'), Shannon evenness index (E), Cereal rye treatment consistently exhibited the lowest AWCD at all incubation periods.Native vegetation treatment had significantly higher AWCD after 24, 48 h, and 72 h than other cover crop treatments followed by CC-mix (p < 0.05) (Figure 1A).After 96 h of incubation, both native vegetation and CC-mix had similar AWCD, which was higher than the remaining cover crop treatments.In regard to fertilizer source treatments, AWCD was highest with PL (Poultry litter) treatment, followed by chemical fertilizer.The fertilizer treatment exhibited the lowest AWCD after a 96-h incubation period (Figure 1B).During the entire incubation period, the AWCD values for each cover crop followed the following pattern NV > CC-mix > VT > WT > CR.These results suggest that the metabolic activity of soil microbial communities in utilizing carbon substrates differs among cover crops and fertilizer sources.

Influence of Cover Cropping and Fertilizer Source Treatments on Microbial Metabolic Diversity Indices
Absorbance measured as OD values at 96 h were used for the estimation of functional diversity parameters and average well-color development since these indicated the optimal range of absorbance.The influence of cover crops and fertilization on functional diversity indices such as the Shannon diversity index (H ), Shannon evenness index (E), McIntosh index (U), substrate richness (R), and inverse Simpson index (1/D) of soil microbial communities is shown in Table 1.
All the functional diversity indices except 1/D demonstrated significant differences (p < 0.05).The H , E, R, and 1/D indices were most pronounced in the CC-mix treatment, succeeded by native vegetation and vetch treatments, whereas the values observed in the cereal rye treatment were considerably lower.Additionally, the native vegetation treatment demonstrated the highest AWCD and U values.
In the poultry litter treatment, the highest values were observed across all functional diversity indices, including AWCD, and were closely followed by the chemical fertilizer treatment (Table 1).Notably, significant differences were observed between fertilizer source treatments for diversity indices, with the exception of 1/D.The microbial functional diversity indices also exhibited significant variations, which were influenced by the interaction between cover crop and fertilizer source treatments.Notably, the CC-mix + PL treatment demonstrated the highest values for all the diversity indices, followed by CC-mix + CL treatment for AWCD, U, and R indices and NV + NO for H and E. The soil samples under CR + NO treatment, on the other hand, displayed lower values for all the diversity indices.These results generally align with the patterns observed for the main effects of cover crop, except for the WT + PL treatment, which displayed higher values compared to the vetch treatment in combination with other fertilizer sources (Table 2).

Classification of Carbon Substrate Utilization Categories in Biolog EcoPlate
The Biolog EcoPlate comprises 31 carbon substrates, which are classified into six distinct categories: carbohydrates, carboxylic acids, amino acids, polymers, amines, and phenolic acids (Table 3).The metabolic activity of soil microbial populations on these six different categories of carbon substrates distinctly varies among cover crops and fertilizer source treatments.Figure 2A,B showed that the relative utilization of different carbon group substrates significantly differs with treatments (p < 0.001).
In the CC-mix treatment, higher utilization of carbon substrates in the polymer and carboxylic acid groups was observed, while the microbial communities from native vegetation treatment showed the highest consumption of carbohydrates, amino acids, amines, and phenolic compounds.Among fertilizer source treatments, the chemical fertilizer treatment showed greater carbon substrate utilization across all carbon groups except polymers, while the highest consumption of polymer group carbon substrates was evident in the poultry litter treatment.Collectively, both treatments indicated that carbohydrates and polymers were the primary carbon substrate groups harnessed by soil microbial communities residing within their respective treatments.The metabolic activity of soil microbial populations on these six different categories of carbon substrates distinctly varies among cover crops and fertilizer source treatments.Figure 2A and B showed that the relative utilization of different carbon group substrates significantly differs with treatments (p < 0.001).In the CC-mix treatment, higher utilization of carbon substrates in the polymer and carboxylic acid groups was observed, while the microbial communities from native vegetation treatment showed the highest consumption of carbohydrates, amino acids, amines, and phenolic compounds.Among fertilizer source treatments, the chemical fertilizer treatment showed greater carbon substrate utilization across all carbon groups except polymers, while the highest consumption of polymer group carbon substrates was evident in the poultry litter treatment.Collectively, both treatments indicated that carbohydrates and polymers were the primary carbon substrate groups harnessed by soil microbial communities residing within their respective treatments.

Analysis of Carbon Substrate Utilization by Soil Microbial Communities Using Principal Component Analysis (PCA)
To establish the degree of separation of treatments in relation to carbon substrate utilization, each of the cover crop and fertilizer treatments was subjected to principal component analysis and Bray-Curtis ordination.Principal components-PC1 and PC2 scores for 31 carbon substrates are given in Table 3.The higher loading scores represent larger

Analysis of Carbon Substrate Utilization by Soil Microbial Communities Using Principal Component Analysis (PCA)
To establish the degree of separation of treatments in relation to carbon substrate utilization, each of the cover crop and fertilizer treatments was subjected to principal component analysis and Bray-Curtis ordination.Principal components-PC1 and PC2 scores for 31 carbon substrates are given in Table 3.The higher loading scores represent larger effects of carbon substrate on the principal components.The first and second principal components described 34.4% and 12% variation, respectively.As shown in Figure 3A, there were distinct differences in the utilization of soil microbial communities from different cover crop treatments (PERMANOVA, p = 0.009).
The results from the principal component analysis showed that PC1 was positively correlated to carbon sources, including L-Asparagine, Tween 40, D-Malic acid, D-Galacturonic acid, D-Mannitol, Putrescine, Gamma-Amino-Butyric acid, N-Acetyl-D-Glucosamine, L-Arginine, and D-Glucosaminic acid, as evidenced by their high loading scores exceeding 0.5 (Table 3).On the other hand, PC2 showed a positive correlation with Glucose-1-Phosphate, Alpha-D-Lactose, and I-Erythritol.Collectively, these results highlight the significant contribution of carbohydrates, amino acids, and carboxylic acids to the two principal components.

Association between Soil Microbial Metabolic Diversity and Soil Physicochemical Characteristics
The effects of cover crops and fertilizer source treatments on soil physicochemical characteristics are given in Table 4.
A canonical correspondence analysis (CCA) was carried out on the soil physicochemical properties and carbon substrate groups (Figure 4).The initial two CCA axes accounted for 79.5% and 12.1% of the variability in data, respectively.The cumulative amount of interpretation of the first 2 CCA axes reached 92.06 %, which could reflect the association between microbial community species based on carbon group utilization and soil parameters.The soil pH = 0.7367 indicated a substantial positive correlation with the major axis of the CCA, whereas total carbon TC = −0.28704,total nitrogen TN = −0.31267,and POXC = −0.25727displayed a negative association.
effects of carbon substrate on the principal components.The first and second principal components described 34.4% and 12% variation, respectively.As shown in Figure 3A, there were distinct differences in the utilization of soil microbial communities from different cover crop treatments (PERMANOVA, p = 0.009).The results from the principal component analysis showed that PC1 was positively correlated to carbon sources, including L-Asparagine, Tween 40, D-Malic acid, D-Galacturonic acid, D-Mannitol, Putrescine, Gamma-Amino-Butyric acid, N-Acetyl-D-Glucosamine, L-Arginine, and D-Glucosaminic acid, as evidenced by their high loading scores exceeding 0.5 (Table 3).On the other hand, PC2 showed a positive correlation with Glucose-1-Phosphate, Alpha-D-Lactose, and I-Erythritol.Collectively, these results highlight the significant contribution of carbohydrates, amino acids, and carboxylic acids to the two principal components.

Association between Soil Microbial Metabolic Diversity and Soil Physicochemical Characteristics
The effects of cover crops and fertilizer source treatments on soil physicochemical characteristics are given in Table 4.The CCA also revealed a strong positive correlation between carbon substrate groups and primary axes, such as carbohydrates (0.94), polymers (0.67), and a bit less to carboxylic acids (0.39) and amino acids (0.23).
ical properties and carbon substrate groups (Figure 4).The initial two CCA axes accounted for 79.5% and 12.1% of the variability in data, respectively.The cumulative amount of interpretation of the first 2 CCA axes reached 92.06 %, which could reflect the association between microbial community species based on carbon group utilization and soil parameters.The soil pH = 0.7367 indicated a substantial positive correlation with the major axis of the CCA, whereas total carbon TC = −0.28704,total nitrogen TN = −0.31267,and POXC = −0.25727displayed a negative association.

Discussion
This is one of the few experiments investigating the synergistic effects of cover cropping and fertilizer sources on both the physicochemical properties of the soil and the functional diversity of the microbial community.We hypothesized that the synergistic effects of cover crops and fertilizer sources would cause significant differences in the metabolic diversity and distribution of soil microbial populations.The results showed that different cover crop and fertilizer source treatments exerted different impacts on the metabolic activities of soil microorganisms in soybean fields.
AWCD is proportional to the number and species diversity of soil microbes that can metabolize carbon substrates, and it indicates the carbon source catabolic activity of microbial communities to utilize the carbon substrates [37,61].In the present investigation, AWCD of soil microbial populations from cover crop and fertilizer source treatments exhibited classic microbial growth over time as the duration of incubation increased (Figure 1A,B), and these results were in accordance with other findings [40,53].In this experiment, during the entire incubation period, native vegetation and cover crop mix (cereal rye + mustard) exhibited higher AWCD values than wheat, vetch, and cereal rye.Elevated values of average wellcolor development (AWCD) indicate a higher level of metabolic activity exhibited by soil microbial populations [62].Additionally, plots with cover crop mix treatment showed high Shannon diversity and substrate richness, whereas native vegetation showed high U and E. Drost et al. (2020) reported multiple cover crop species results in elevated amounts of substrate niches that increased the functional diversity of soil microbial populations [30].
Incorporating a mix of two or more species in cover cropping can provide additional benefits by promoting both the abundance of beneficial soil microbial populations and the diversity of the soil microbiota [63].The quality of the residue from cover crops, such as cereal rye, mustard, vetch, and wheat, differed significantly.Additionally, different types of fertilizer sources were used in our two-year study.Despite these variations, our findings showed that the community-level physiological profiling using the Biolog EcoPlates exhibited convergence in diversity indices between cover crop mix and native vegetation treatments.
Studies have shown that mixed cover crops favor the release of certain allelochemicals, primary metabolites, and secondary metabolites significantly alter the soil microbial abundance and their functions [53].However, cover crop biomass, quality, different soil types, C:N ratio of residues, environmental circumstances, and their interactions significantly influence the abundance and functional diversity of soil microbial populations [64].Consistently, the cereal rye treatment exhibited the lowest average well-color development (AWCD) values, as well as lower diversity index values, which is consistent with previous findings [28].The AWCD of the vetch cover crop was moderate relative to other treatments, as were the indexes of microbial functional diversity.Vetch, being a leguminous cover crop, has a low C: N ratio and decomposes rapidly, leading to the availability of more nitrogen in the soil.This is reflected in the significant differences observed in the soil's physicochemical properties compared to other treatments (p = 0.0018, Table 4).
The AWCD values of fertilizer source treatments were also significantly different.Plots with poultry litter treatment showed higher AWCD than chemical fertilizer and no fertilizer (Table 1); these findings are similar to the other discoveries.Many researchers reported application of organic amendments results in elevated soil microbiological activity as determined by dehydrogenase function, which reflects the complete range of oxidative activity of microbial populations [65,66].The utilization of organic amendments has been widely recognized for its beneficial impact on the structure and functioning of soil microbial communities [67][68][69][70].A study conducted by Marinari et al. (2000) [67], it was found that the application of organic fertilizers resulted in enhanced soil microbial activity.The researchers attributed this effect to a synergistic interaction between the microbial populations present in the soil and the organic materials used as fertilizers.
The combined effects of fertilizer source and cover crop demonstrated significant differences in functional diversity parameters.Specifically, the soil samples from the CCmix + poultry litter treatment plot exhibited higher values across all functional diversity indices.This outcome is likely attributed to the synergistic benefits derived from the combination of mixed cover crop species and the application of poultry litter as a fertilizer.Furthermore, it is interesting to note that the soil samples from the wheat + poultry litter treatment plot also displayed higher values compared to the main effects of wheat and even outperformed the main effects of the vetch treatment.This suggests that the addition of poultry litter to the wheat plot might have a greater contribution to the microbial communities in that specific plot.
Additionally, our observations revealed that cereal rye, as a cover crop in the main effects analysis, displayed lower values in the functional diversity indices.Interestingly, this pattern persisted when examining the combined effects of cereal rye with different fertilizer sources.The consistently lower values suggest that the influence of cover crops, as a collective factor, exerts a stronger impact on the functional diversity indices compared to individual cover crop effects alone.This finding highlights the significance of considering the combined effects of cover crops and fertilizer sources in shaping microbial functional diversity in the soil.
A total of 31 carbon substrates are classified into six categories, and the metabolism patterns of these six categories were significantly different among treatments (p < 0.001).Cover crop mix had a higher metabolism rate of polymers and carboxylic acids, whereas microbial communities under native vegetation treatment showed higher utilization of carbohydrates, amines, and phenolic compounds.Overall, both cover crop and fertilizer source treatments showed that carbohydrates and polymers were mainly used by the microbial populations in the soil.These results were in accordance with Lan et al. (2019) [53].
Nivelle et al. (2016) reported higher metabolism of phenolic compounds and carbohydrates under no-tillage with oats, phacelia, and flax [71].Many researchers have reported a high degree of carbohydrate utilization by soil microbial communities in the uppermost soil layer (0-10 cm) [28,72], owing to the presence of sugars and organic matter generated by the decomposition process and rhizodeposits.This abundance of carbohydrates in the topsoil may offer greater accessibility to soil microbial communities.The higher utilization of polymers (Glycogen, alpha-cyclodextrin, tween-80, and tween-40) was observed by CC-mix treatment.Generally, these carbon substrates are more stable or recalcitrant; specifically, tween-40 and tween-80 don't represent plant-derived products [73]; furthermore, they are often found in processed organic matter [74].The microbial communities also exhibited better utilization of amines which includes phenylethylamine and putrescine.Additionally, phenylethylamine was discovered in fermented soybean, whereas putrescine was produced in a variety of soybean plant material, including immature roots [75,76].These results infer that the major crop also has a role in structuring metabolic patterns of soil microbial populations.
Fertilizer source treatments also exhibited significant differences among the utilization of carbon groups.Poultry litter had higher consumption of polymers, whereas chemical fertilizer had higher utilization of carbohydrates, carboxylic acids, amino acids, amines, and phenolic compounds.Guanghua et al. (2008) reported chemical fertilizer treatment metabolized more polymers, whereas chemical fertilizer + FYM showed higher utilization rates in carboxylic acids, which contradicted our findings [77].
Through CCA analysis, it was determined that total carbon, total nitrogen, and POXC had a negative correlation with carbon substrate categories, indicating that these parameters have an effect on the metabolism of different carbon substrate categories in different cover crops and fertilizer sources.The soil characters' pH and EEGSP positively correlated with the primary axis.Similar results were reported by Lan et al. (2019) [53]; however, they conducted an investigation in forest soils of Eucalyptus trees.
The present investigation revealed varying utilization patterns of the six categories of carbon substrates among the soil microbial communities subjected to different treatments (Figure 2A,B); this may be attributed to the varied structures of soil microbial communities resulting from the different treatments involving cover crop and fertilizer source.However, carbohydrates were found to be the dominant carbon substrate utilized across all treatments, as indicated by the results of principal component analysis (PCA).The PCA results further indicated that carbohydrates made the largest contribution to both PC1 and PC2, followed by amino acids and carboxylic acids.Polymers are also highly metabolized by soil microbial populations in CC-mix treatment; however, PCA didn't reveal in the first two components.
Treatments strongly influenced microbial functional diversity distribution in both ordination plots; the differences in utilization of carbon substrates across cover crop and fertilizer source plots significantly influenced microbial functional diversity distribution (Figure 3A,B).The PCA separation suggests the significant combined effects of cover cropping and fertilizer sources on soil microbial metabolic activity and functional diversity (Figure 3B).The Bray-Curtis ordination plot clearly shows the grouping of wheat, cover crop mix, and native vegetation, whereas vetch and cereal rye formed separate groups.These results infer that grass species cover crops harboring soil microorganisms with similar catabolic profiles, whereas vetch, which is leguminous species distributed separately.The effect of fertilizer sources on vetch and cereal rye is also evident in the PCA ordination plot (Figure 3B).The use of different fertilizer sources with cereal rye resulted in the formation of three separate groups, indicating that the application of poultry litter, chemical fertilizer, or no fertilizer creates distinct microbial communities with unique catabolic profiles.Similarly, vetch with poultry litter and no fertilizer formed a single group, while vetch with chemical fertilizer formed a separate group.These results infer that combined impacts of cover cropping and fertilizer source significantly influence the metabolic activity and community structure of soil microbiota.To further confirm our findings, we conducted PERMANOVA, which revealed that both cover crops and fertilizer sources significantly influenced the metabolic patterns of soil microbial communities (p = 0.0012 (PCA), p = 0.009 (Bray-Curtis)).Our results are consistent with previous studies by Nair and Ngouajio (2012), Bucher and Lanyon (2005), and Gomez et al. (2005) that investigated the impact of organic manures such as compost, dairy manure, and FYM on the metabolic activity of soil microbiota [28,78,79].
Recent breakthroughs in plant-microbe interaction studies have demonstrated that plants may alter their rhizosphere microbial communities [80,81].In addition to this, several primary metabolites, including sugars and organic acids, are secreted by soybean into the rhizosphere [82,83].In addition to soybean root exudates, the root system of cover crops may supply alternate carbon pools that impact soil microorganisms' metabolic activity and community structure in soybean fields.

Conclusions
Our study demonstrated that the Biolog EcoPlates approach utilizing communitylevel physiological profiling methodology was capable of detecting short-term changes resulting from management techniques.Our results showed that organic amendment considerably enhanced microbial functional diversity.We also demonstrated how the implementation of cover cropping impacted the functional diversity and metabolic capacity of the soil microbiota.Soil microbes occur in the active/labile portion of soil organic matter participating in various soil ecological services; thus, soil microbial metabolic activities may respond more rapidly than physicochemical properties of soils to changes in management techniques such as the use of cover crops, organic amendments application, and tillage operations or environmental circumstances [78].Our results also demonstrated that cover cropping with multispecies can enhance and alter the metabolic diversity of soil microbial communities and the association between cover crop species and microbial communities' metabolic patterns.Our findings suggest that the incorporation of cover crops, especially muti-species, could potentially benefit dryland soybean production through enhanced soil microbial metabolic diversity.However, it is important to note that these results only reflect a portion of the entire microbial communities in the soil to a limited extent; some microbial species don't metabolize carbon substrates in the Biolog EcoPlates, and some might be in a dormant state and couldn't be accounted for the overall metabolic diversity of microbial communities.Novel approaches like soil metagenomics (amplicon sequencing and shotgun metagenomics) should be utilized to further explore how carbon substrate associated with management practices affect the abundance, species diversity, and in-depth metabolic patterns of soil microbial populations.

Figure 1 .
Figure 1.Changes in average well color development (AWCD) of the soil microbial communities in cover crop (A) and fertilizer source (B) treatments with incubation time.

Figure 1 .
Figure 1.Changes in average well color development (AWCD) of the soil microbial communities in cover crop (A) and fertilizer source (B) treatments with incubation time.

Figure 2 .
Figure 2. Utilization of carbon substrate categories by soil microbial communities from the cover crop (A) and fertilizer source (B) treatments.Bars with letters denotes significance at the level 0.05 level using Fisher's protected LSD.(CR = Cereal rye; CC-mix = Cereal rye and mustard; NV = Native vegetation; VT = Vetch; WT = Wheat, CF = Chemical fertilizer; PL = Poultry litter; and NO = No fertilizer).

Figure 2 .
Figure 2. Utilization of carbon substrate categories by soil microbial communities from the cover crop (A) and fertilizer source (B) treatments.Bars with letters denotes significance at the level 0.05 level using Fisher's protected LSD.(CR = Cereal rye; CC-mix = Cereal rye and mustard; NV = Native vegetation; VT = Vetch; WT = Wheat, CF = Chemical fertilizer; PL = Poultry litter; and NO = No fertilizer).

Figure 4 .
Figure 4. Canonical Correspondence Analysis (CCA) ordination plot revealing the relationship between cover crop treatments and soil physicochemical properties.

Table 1 .
Main effects of cover crop and fertilizer source on soil microbial community diversity indices.Variable in column with no letters is not significant at the 0.05 level using Fisher's protected LSD; Standard error in parenthesis.2CC-mix=Cereal rye and mustard; 3 CC*FT = Interaction between cover crop and fertilizer source (ns = non-significant, s* = significant.AWCD = Average color development; H = Shannon Diversity Index; 1 E = Shannon evenness index; U = Mcinthosh Index; R = Substrate richness; 1/D = Inverse Simpson Diversity, CL = Chemical fertilizer, PL = Poultry litter.

Table 2 .
Combined effects of cover crop and fertilizer source on soil microbial community diversity indices.
1Variable in column with no letters are not significant at the 0.05 level using Fisher's protected LSD; Standard Error in parenthesis.

Table 3 .
Carbon substrates present in Biolog EcoPlates related to PC1 and PC2.

Table 4 .
Main effects of cover crops and fertilizer source treatments on soil physicochemical characteristics.

Table 4 .
Main effects of cover crops and fertilizer source treatments on soil physicochemical characteristics.Variable in column with no letters are not significant at the 0.05 level using Fisher's protected LSD; Standard Error in parenthesis.2CC-mix=Cereal rye and mustard; 3 CC*FT = Interaction between cover crop and fertilizer source (ns = non-significant).TC = Total carbon; TN = Total nitrogen; WSA = Water-stable aggregate; EEGSP = Easily extractable soil glomalin; POXC = Permanganate oxidizable carbon.