Yang Sean Xiao
Bo Zhou
Zhuangzhuang Han1
Shenzhou Liu
Wenzhi Zeng- 1College of Water Resources and Civil Engineering, China Agricultural University, Beijing, China
- 2Engineering Research Center for Agricultural Water-Saving and Water Resources, Ministry of Education, Beijing, China
- 3State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan, China
- 4Guangxi Hydraulic Research Institute, Nanning, China
- 5College of Water & Architectural Engineering, Shihezi University, Shihezi, China
Zinc is an essential minor element for rice growth and human health, which can also change the structure of the microorganisms. However, it remains unclear for the effects of zinc fertilizer on microbiome function in agricultural soils and crops. To solve this research gap, we investigated the relationship between improving rice (Oryza sativa L.) yield, Zn concentration, soil microbial community diversity, and function by the application of Zn fertilizer. The field trials included three rice varieties (Huanghuazhan, Nanjing9108, and Nuodao-9925) and two soil Zn levels (0 and 30 kg ha–1) in Jiangsu province, China. As a test, we studied the variety of soil bacterial composition, diversity, and function using 16S rRNA gene sequencing. The results showed that soil Zn application reduced the diversity of microbial community, but the bacterial network was more closely linked, and the metabolic function of bacterial community was improved, which increased the grain yield (17.34–19.52%) and enriched the Zn content of polished rice (1.40–20.05%). Specifically, redundancy analysis (RDA) and Mantel’s test results revealed soil total nitrogen (TN) was the primary driver that led to a community shift in the rice rhizosphere bacterial community, and soil organic carbon (SOC) was considered to have a strong influence on dominant phyla. Furthermore, network analysis indicated the most critical bacterial taxa were identified as Actinobacteria, Bacteroidetes, Proteobacteria, and Chloroflexi based on their topological roles of microorganisms. KEGG metabolic pathway prediction demonstrated that soil Zn application significantly (p < 0.05) improved lipid metabolism, amino acid metabolism, carbohydrate metabolism, and xenobiotic biodegradation. Overall, their positive effects were different among rice varieties, of which Nanjing-9108 (NJ9108) performed better. This study opens new avenues to deeply understand the plant and soil–microbe interactions by the application of fertilizer and further navigates the development of Zn-rich rice cultivation strategies.
Introduction
Rice (Oryza sativa L.), as the main source of zinc (Zn) intake for humans (Qin et al., 2012), is one of the most important food crops and staple food for more than one-third of the world’s population (Kasote et al., 2021) and is recognized as a major food crop in Asia due to its high demand and high planting rate (Nayar, 2014). However, Zn deficiency is a widespread micronutrient disorder in rice, resulting in reduced rice grain yield and poor nutrition quality (Fageria et al., 2002), continuous application of fertilizer and their low efficiency, especially N and P, and has caused environmental degradation (Bindraban et al., 2020). These problems have posed a major threat to sustainable rice production and food security (Rahman et al., 2022). Zn is one of the irreplaceable minor elements required for rice growth and other food crops (Broadley et al., 2007). It not only promotes plant metabolism, but also acts as an activator of cell division, protein synthesis, and gene expression (Yamaji et al., 2013). At the same time, as one of the essential micro-nutrients in the normal operation of various organs (Chasapis et al., 2012), Zn plays a vital role in human health (Ma et al., 2008). The average Zn content of rice in China is only 17 mg kg–1, which can be as low as 9.2 mg kg–1 in cooked rice (Chinese Center for Disease Control and Prevention). Considering the weight ratio, it was almost impossible to meet the recommended criteria for daily Zn intake (7.5–12.5 mg d–1) (Wang et al., 2016). The application of Zn fertilizer, such as zinc sulfate (ZnSO4), is a quite effective method to correct the deficiency of Zn and further improve in grain yield and increase grain Zn concentration (Khan et al., 2003; Rehman et al., 2018). Zn application methods mainly include soil application, foliar application, seed priming, and seed coating (Nadeem and Farooq, 2019; Nazir et al., 2021). Recent research has found that the application of 5–45 kg ZnSO4⋅7H2O ha–1 to the soil increased grain yield (13–60%) (Khan et al., 2002; Huang et al., 2019), while foliar application seems more effective in boosting grain Zn concentration (Fageria et al., 2009). In addition, other fertilizers (e.g., nitrogen fertilizer and bio-activated organic fertilizer), rice genotypes, and soil microorganisms are also influential factors in determining the Zn-rich capacity of rice (Ghasal et al., 2018; Wang et al., 2021). For instance, Hussain et al. (2020) used Bacillus sp. AZ6 as an inoculant of bio-activated organic fertilizers and found significant improvement in the crop growth, physiology, yield, and Zn content.
Soil microorganisms play important roles in plant performance by improving mineral nutrition, soil fertilizer, and soil quality (Weese et al., 2015; Gu et al., 2019). The nitrogen-fixing and nutrient mineralization processes carried by soil microorganisms can metabolize recalcitrant forms of N, P, and S to literate these elements for plants (Jacoby et al., 2017). In recent years, some scholars have proposed plant growth-promoting rhizobacteria (PGPR) as an alternative and eco-friendly technology that can enhance zinc solubilization and its availability to plants (Zeb et al., 2018; Bhatt and Maheshwari, 2020). Some bacteria, including Acinetobacter, Bacillus, and Pseudomonas, have been reported to solubilize zinc (Kumar et al., 2019). Thus, the production of biological fertilizers containing beneficial microorganisms may be an effective alternative to chemical fertilizers. In addition, potential interactions were found among bacterial communities, soil chemical properties, and crop growth status, which affected soil fertility and nutrient resorption characteristics of rice, leading to grain yield and Zn content variations in rice (Wang et al., 2021; Zhang et al., 2021). However, soil microorganisms are also shaped by many factors, such as farming manipulation, plant varieties, and soil nutrition (Berg and Smalla, 2009; Hartman et al., 2018). Changes of soil microbial community can be used to assess the quality of soil ecosystems (Pablo et al., 2014). Zn, as required nutrients for microorganisms, can be toxic to microorganisms by the displacement of essential metals from their native binding sites or through ligand interactions at high levels (Bruins et al., 2000). Co-occurrence network analysis can be used to explore changes in species interactions or microbial responses under environmental stress or human disturbance and to identify for key species (Deng et al., 2012; Bissett et al., 2013). However, soil microorganism and co-occurrence patterns to sustainable agricultural cultivation in rice are largely unknown.
In view of the above, we attempt to study the soil microorganisms, grain yield, and grain Zn concentration of rice affected by Zn fertilizer under sustainable agricultural practices (using microbial fertilizers instead of NPK fertilizer). The main objectives of this study are to test the following hypotheses: (i) the effect of Zn fertilizers differs among the common rice varieties in Jiangsu, China; (ii) the application of Zn fertilizers plays a significant role in shaping rhizosphere bacterial community and altering potential function. These findings obtained were beneficial to a more sustainable way of promoting Zn in rice through the regulation of the rhizosphere microbial community.
Materials and methods
Experimental background and zinc fertilizer input treatments
This study was conducted in Yancheng City, Jiangsu Province, China (33°11′ N, 119°52′ E). This area has a subtropical monsoon climate, the average annual air temperature is 15.45 ± 0.57°C, and the average annual precipitation is 1049.1 ± 304.1 mm (from 2011 to 2020; China National Knowledge Infrastructure, 2022). The value of available soil zinc in the field was detected as 1.07 ± 0.06 before the experiment. Three kinds of high-quality, high-yielding rice varieties (i.e., HHZ, Huanghuangzhan; NJ9108, Nanjing-9108; ND9925, Nuodao-9925) cultivated locally were chosen as the tested plants, and the critical limit of zinc (Zn) in the soil for rice ranges from 0.80 to 0.85 μg/g (Rahman et al., 2022). Two levels of basal Zn fertilizers (0 and 30 kg ha–1) were applied as ZnSO4⋅H2O ha–1 (effective Zn content 34.5%) before rice transplantation (25 May 2020). Therefore, six treatments were contained for this study. For each treatment, approximately 6000 m2 of field was divided into seven plots as replicates (see Supplementary Figure 1 for sampling diagram). Foliar Zn and microbial fertilizers were applied in all the treatments, where the foliar Zn application (3 kg ZnSO4⋅H2O ha–1 dissolved in 750 L of water) was sprayed at the booting stage, early grain-filling stage, and late grain-filling stage, respectively. The microbial fertilizers of 15,000 kg ha–1 [Manufacturer: Hebei Sanfeng Biofertilizer Co., Ltd; organic matter ≥ 70%, effective live microorganism (bacterial and fungal mix) ≥ 5 × 108 g–1], which contained a high content of composite microorganisms, mainly composed of Bacillus spp. and Trichoderma spp. were applied with basal fertilizer before rice transplantation. All the experimental treatments had the same seed planting and transplanting times on 21 April 2020 and 30 May 2020, respectively. In each hole, 3–5 seedlings were transplanted. The planting density was converted to 30 holes per m2 uniformly, and irrigation and tillage management were consistent across all treatments. The details of treatments in this study are shown in Table 1.
Table 1. Treatments, abbreviations, and fertilization status.
Testing of rice yield and zinc content in polished rice
A 0.25 m2 sample square of consistently growing rice plants were randomly collected at maturity for each treatment per plot. All unfilled spikelets were removed, and others were recorded. The grain yield was de-enzyme using an oven at 105°C for 30 min and then dried at 70°C to a constant weight. Subsequently, the grains were rolled out, and the dry matter mass of the rice grain yield was weighed. Following this treatment, the grain yield was adjusted to the standard moisture content of 0.14g H2O g–1 fresh weight and converted to yield per hectare by area ratio.
We determined the zinc content in polished rice according to the current Chinese National Food Safety Standard—Determination of Zinc in Food (GB 5009.14-2017). Briefly, 2 g of polished rice was put in a conical bottle with 10 mL of nitric acid (HNO3) and 0.5 mL of perchloric acid (HClO4) and then carried out the digestion on an adjustable electric heating furnace. After digestion of the sample, the absorbance was measured at 213.9 nm by flame atomization, which then was quantitatively compared with the standard series.
Rhizosphere sampling and chemical analysis
Samples of rhizosphere from each treatment group were collected from each replicate plot. For each plot, 5 sampling points were chosen in S-shaped pattern (Supplementary Figure 1), and a total of 10–12 rice plants from each sampling point were collected to ensure that the root system was as intact as possible. These five sub-samples were mixed well and used as one replicate. The sampled soil specimens included both loosely and tightly bound root soil. Loosely attached root soils from plots numbered 2, 4, and 6 were gently removed, passed through a 2 mm sieve, and then used to analyze the soil chemical properties (totally 18 samples). Tightly bound soils from plots 1–7 were separated using a sterile brush and then were re-suspended in 10 mL of a sterile 0.8% NaCl by vigorous shaking for 3 min (Angelo-Picard et al., 2004). These samples were used as rhizosphere soil (a total of 42 samples) and stored at –80°C for high-throughput sequencing. RDA and Mantel’s test were performed with bulk soils from plot numbered 2, 4, and 6 and rhizosphere soil from the same plots.
The soil chemical properties involved the soil pH, soil organic carbon (SOC), soil total nitrogen (TN), soil effective nitrogen (AN), soil total phosphorus (TP), soil total potassium (TK), cation exchange capacity (CEC), and soil available Zn (AZ) values. The Supplementary Part 1 provided the test methods, and the results are shown in Table 2.
Table 2. Zinc content in polished rice and grain yield of three rice cultivars.
Soil microbial and bioinformation analysis
DNA extraction, PCR amplification, high-throughput sequencing, and Real-time quantitative PCR
DNA was extracted from rhizosphere samples of six treatments (3 varieties × 2 Zn application levels; ca, 0.35 g), 42 samples in total (7 replicates per treatment), using a Fast DNA TM Spin Kit (MP Biomedicals LLC, United States) according to the manufacturer’s instructions. The V3-V4 region of the 16S rRNA gene fragment was amplified using a primer set 338F (5’ ACT CCT ACG GGA GGC AGC AG-3’) and 806R (5’-GGA CTA CHV GGG TWT CTA AT-3’) (Sangon Bioengineering, Co., Ltd, China). Paired-end format of purified amplicons was mixed well and sequenced using an Illumina Miseq sequencing platform according to the standard protocols by Shanghai Majorbio Bio-pharm Technology Co., Ltd (Shanghai, China). Flash (v 1.2.111) was used to remove the primer sequences and contig paired ends. QIIME, an open-source bioinformatics pipeline (v 1.9.12), was used to check the quality of the raw sequences. Clustering of high-quality sequences at 97% similarity to generate operational taxonomic units (OTUs) was conducted using UPARSE3. Valid sequences were matched against the SILVA database4 followed by a taxonomic comparison of representative OTU sequences at the 97% similarity level using the RDP classifier Bayesian algorithm (Stackebrandt and Goebel, 1994), with threshold of 70% was used to discriminate the phylum, order, family, and genus levels. A total of 9,500 OTUs were obtained by randomly sub sampling with the minimum number of sequences from the high-throughput test. Real-time quantitative PCR, targeting the bacterial 16S gene, was conducted with the primer sets Eub338 forward primer (5’-ACT CCT ACG GGA GGC AGC AG-3’) and Eub806 reverse primer (5’-GGA CTA CHV GGG TWT CTA AT-3’), pMD-18T (2692 bp, Takara Biotechnology Co., Ltd., Dalian, China) Easy vector was used to clone PCR products. The details are in Supplementary Part 2.
Co-occurrence networks construction
Co-occurrence network analysis of the rhizosphere bacterial communities under each treatment was conducted for six treatments (XRZS0, XRZS2, JRZS0, JRZS2, NRZS0, and NRZS2) with the same similarity threshold of 0.960 using the MENA platform on the Environmental Genomics Institute Website of the University of Oklahoma5 (Zhou et al., 2010). OTUs were screened and saved with the sequence numbers exceeding or equaling to 5 in at least three samples to filter for rare OTUs in each treatment. The results obtained from the MENA platform were visualized by Cytoscape (v 3.8.0) and Gephi software (v.0.9.2). The details can be found in the Supplementary Part 3.
Potential function prediction
The functional genes of rice rhizosphere bacteria were conducted by PICRUSt 2 (v 1.1.06), and KEGG functional annotation was used to obtain information on the annotation of OTUs at each KEGG function level and abundance (Langille et al., 2013; Wu et al., 2022). STAMP software (v 2.1.3) was used for statistical hypothesis tests and exploratory plots (Parks et al., 2014) using Welch’s t-test to compare the differences in two groups.
Statistical analysis
Soil chemical properties were recorded using Excel (v.2019). Two-way analysis of variance (ANOVA) was used to analyze whether basal Zn application, rice cultivars, and their interaction had a significant effect on soil chemical properties, grain yield, and Zn content in polished rice. The least significant difference (LSD) tests were performed at the p < 0.05 level. Alpha-diversity was calculated using Shannon, Evenness, and Sobs indices by MOTHUR (Schloss et al., 2009) (Wilcoxon rank-sum test. v 1.30.27). The differences in the microbial community structure were analyzed using principal coordinate analysis (PCoA) based on the Bray–Curtis distance (Wang et al., 2019). The linear discriminant analysis effect size (LEfSe, Afgan et al., 2018) was determined to identify the key rhizosphere bacteria (Kruskal-Wallis test p < 0.05 and logarithmic LDA score (log10) > 3.7). Redundant and canonical correlation analyses were conducted to illustrate the relationships between rhizosphere bacterial communities and the soil chemical properties by Canoco 5 (Legendre et al., 2011). Spearman’s correlations between key microbial communities (genus level) and soil chemical properties aimed to determine the relationship between key genus and environmental variables (at p < 0.05 and p < 0.01, respectively).
Results and analysis
Shift in rhizosphere microbial community structure and biomarkers after zinc treatment
As shown in Figure 1A, 16S rRNA gene abundance (per g of soil) for XRZS2 treatment (5.223 × 109) was 2.56% lower than XRZS0 (5.383 × 109), while JRZS2 and NRZS2 were both higher than the control treatment group (5.383 × 109 vs. 4.762 × 109 and 5.223 × 109 vs. 5.037 × 109), resulting in JRZS2 increasing by 13% and significantly differing from JRZS0 (p < 0.05, ANOVA Duncan’s test).
Figure 1. (A) Different soil Zn application and rice cultivars on the abundance of the rhizosphere soil bacterial gene copy numbers per gram sample; (B) microbial community composition of dominant bacterial phyla in rice rhizosphere soil. Phyla with a total relative abundance of < 0.5% are grouped in “Others”; (C) Sobs index; (D) Shannon index; (E–G) The volcano plot illustrates the up- and downregulated patterns of the rice rhizosphere bacterial microbiomes in different rice cultivars with Zn addition (ZS2), compared with control groups (ZS0). Each point represents an operational taxonomic unit (OTU), and the color represents up- and downregulated patterns of OTU. The x-axis represents the average abundance of OTU (as counts per million, CPM), and the y-axis represents the log2(fold change). (E) Huanghuazhan (HHZ), (F) Nanjing-9108 (NJ9108); (G) Nuodao9925 (ND9925). (H) Principal coordinate analysis (PCoA) of the rhizosphere soil bacterial community based on Bray–Curtis distance at the OTU level. The ellipses represent the 95% confidence intervals of each treatment; (I) bacterial biomarkers in different treatments and cultivars based on linear discriminant analysis (LDA) using LEfSe analysis. Different colors represent different treatments, and the circles from inside to outside correspond with phyla to genus. Colors coded within the cladogram indicate the taxa with different abundances in treatment by the Kruskal–Wallis test with p < 0.05 and logarithmic LDA score (log10) > 3.7. Genera with a relative abundance of less than 0.2% are not included.
A total of 2,412,863 (57,449 ± 7487/sample) microbial 16S V3-V4 valid reads were obtained from 42 samples through high-throughput sequencing analysis and classified into 9,500 operational taxonomic units (OTUs) at the 97% similarity level. The rarefaction curve indicated that the sequence depths of all tested samples were adequate for further analysis (Supplementary Figure 1). A total of 38 different bacterial phyla were found in the different treatment groups. Taxonomic analysis revealed that the rhizosphere of rice mainly included five phyla (Figure 1B), among which Actinobacteria (22.45–28.54%) was the most dominant phylum followed by Proteobacteria (22.44–31.67%), Chloroflexi (20.17–24.25%), and Acidobacteria (9.26–15.53%). After Zn treatment, Actinobacteria and Gemmatimonadetes significantly increased (5.22–27.09% and 6.36–33.40%, respectively), while Proteobacteria and Firmicutes decreased by 5.49–18.50% and 4.67–46.44%, respectively (p < 0.05, Wilcoxon rank-sum test).
The richness, diversity, and evenness of the bacterial community exhibited a similar decrease trend after soil Zn fertilizer application (ZS2). The sobs index of rhizosphere microorganisms of the three kinds of rice decreased by 0.4–4.92%, with HHZ and NJ9108 showing a significant difference (p < 0.05, Wilcoxon rank-sum test) (Figure 1C). While the Shannon diversity index decreased by 0.11–1.60%, only HHZ did not exhibit a discernible change (Figure 1D). The PCoA based on the Bray–Curtis distance showed that the structure of the bacterial communities differed significantly among the two kinds of soil Zn application treatments (Expect for ND9925) and the three tested cultivars (p = 0.001). Principal components explained 34.63% of all the differences in the data (Figure 1H).
The volcano plot illustrates the variance of bacterial abundance after Zn treatment (Figures 1E–G), large variation in bacterial community was observed in HHZ and NJ9108 with 11.70% (1058/9500) and 10.81% (1027/9500) OTUs changing significantly, while only 1.3% (121/9500) OTUs were observed in ND9925 treatment, which confirms the previous results. Differential OTUs mainly from Anaerolineaceae and Acidobacteria (family) were both significantly enriched and depleted in three treatments. In addition, the family Gemmatimonadaceae exhibits significant enriched and overlapped in three treatments.
To explore the importance of key species more accurately, LEfSe analysis was conducted to identify the biomarkers of bacterial communities in different treatments at the class level (Wilcoxon rank-sum, p < 0.05, LDA score > 3.7). A total of 45 bacterial classifications were enriched in the six treatments. At the class level, XRZS0 was notably enriched in Acidobacteria KD4_96 after applying Zn fertilizer (XRZS2). JRZS0 was significantly enriched in Gammaproteobacteria, Deltaproteobacteria, and Betaproteobacteria after applying soil Zn fertilizer (JRZS2). It was also significantly enriched in Actinobacteria and Gemmatimonadetes at the class level. NRZS0 was enriched in Alphaproteobacteria, while NRZS2 was not enriched. Zn fertilizer treatment enriched the phyla of Acidobacteria, Actinobacteria, and Gemmatimonadetes, families of Gaiellaies and Gemmatimonadaceae, and the genera of Gaiella (Figure 1I).
Dynamic changes of co-occurrence network association
We established six co-occurrence networks to further characterize the influence of Zn application and tested cultivars on bacterial networks (Figure 2 and Table 3). In the treatments without soil Zn application (ZS0), OTUs were mainly clustered into the top two modules, while the number of modules gradually increased with soil Zn application. For HHZ and NJ9108, the number of nodes in the group treated with soil Zn fertilizer decreased slightly (915 vs. 898 and 960 vs. 906), while the number of nodes in ND9925 remained similar (863 vs. 864). Compared to the controls, the number of total links was significantly higher (1546 vs. 1987 and 1382 vs. 1680) and the average path length was lower (11.905 vs. 8.656 and 11.263 vs. 9.485) in the HHZ and ND9925, amounting to no significant difference in any of the NJ9108 treatments. In addition, the higher avgK and avgCC were observed in ZS2. In summary, soil Zn fertilizer application had little effect on HHZ, but increased the complexity of all the bacterial networks.
Figure 2. Visualized networks of microbial co-occurrence in different treatments, (A–F) represent the networks of treatments XRZS0, JRZS0, NRZS0, XRZS2, JRZS2, and NRZS2, respectively. The nodes are colored by module. The module eigengenes from the top 10 large sub-modules explained 41.25–58.81% of the variations of the relative abundance across different replicates, suggesting that these eigengenes represented the module profiles relatively well. Nodes represent an operational taxonomic unit (97% sequence identifies a threshold, OTU). The size of each node is proportional to the node degree. The edges in green and red represent co-occurrence and mutual exclusion patterns among taxa.
Table 3. Topological features of the co-occurrence network of soil microbial communities of the rhizosphere of rice with different treatments and cultivars.
The topological roles of the OTUs are shown in Figure 3. The results demonstrated that most of the nodes (99.11%) were in the peripheral region, and their connections were mainly connected to the nodes in the module. All the nodes with Zi ≥ 2.5 or Pi ≥ 0.62 were determined as the keystone species. Therefore, nodes in the area of connectors (0.24%) and module hubs (0.65%) played a crucial role in the co-occurrence networks. See Supplementary Table 1 for details concerning the bacteria classified as nodes in the connector and module hub area. Among them, the connectors and module hubs of phyla Actinobacteria (e.g., g__Mycobacterium, g__Iamia, and g__Nocardioides), Proteobacteria (e.g., g__Nitrosospira, g__Steroidobacter, g__Shinella, and g__Geobacter), and Chloroflexum (e.g., g__Anaerolinea) were accounted for 15.2, 16.9, and 15.3%, respectively. Thus, these were the most important species in the rhizosphere. Furthermore, in the group ZS2, the number of Bacteroidetes (e.g., g__Chitinophaga and g__Pontibacter) in the connectors and module hubs decreased from 5% in the ZS0 to 0, whereby applying Zn fertilizer reduced the relative abundance of Bacteroidetes, considered as key bacteria, indicating that the species belonging to the phylum Bacteroidetes were more sensitive to Zn fertilizer. Intriguingly, as shown in Figure 3, Bacteroidetes were all derived from HHZ and NJ9108. Applying soil Zn fertilizer significantly changed the topological structure of the key bacteria communities in the rice rhizosphere.
Figure 3. Zi-Pi plot showing the distribution of OTUs based on their topological roles representing an OTU with (ZS2, triangle) or without (ZS0, circle) basal zinc fertilizer. The threshold values of Zi and Pi for categorizing OTUs are 2.5 and 0.62. The Venn diagram shows the similarities and differences between keystone hubs in networks between two treatments. (A) HHZ; (B) NJ9108; and (C) ND9925.
Potential functional prediction of bacterial community
The effect of basal zinc fertilizer application on the metabolic function of soil bacterial communities was investigated. The relative abundance of KEGG pathway (L1 level) in all samples is shown in Figure 4A. The main classifications include (ZS0 vs. ZS2) metabolism (52.24% vs. 52.74%), genetic information processing (15.70% vs. 15.60%), environment information processing (13.96% vs. 13.88%), unclassified (12.78% vs. 12.6%), cellular processes (3.61% vs. 3.5%), organismal systems (0.84% vs. 0.86%), and human diseases (0.85% vs. 0.82%). The results showed that all the above functions produced significant differences between ZS0 and ZS2, with mean proportion of metabolism and organismal systems being upregulated and others downregulated (Figure 4B).
Figure 4. (A) Relative percentage of predicted functional gene abundance (at L1 level by KEGG); (B) KEGG categories (L1 level) differing significantly between ZS0 and ZS2 based on Welch’s t-test, the mean proportion of majority pathway are overabundant within the rhizosphere soil of ZS0 compared to ZS2; (C) expression abundance multiples of bacterial function at the KEGG L2 level, most functions of ZS2 were upregulated compared to ZS0.
To investigate in more detail the functional changes of zinc-rich processes, we calculated expression abundance multiples based on L2 level (Figure 4C). Overall, almost all functions were upregulated in all classifications, for example, lipid metabolism, amino acid metabolism, carbohydrate metabolism, xenobiotic biodegradation and metabolism functions, and membrane transport, while only two functions were downregulated in the circulatory system of organismal systems and neurodegenerative diseases of unclassified, respectively. The results above indicated that the application of exogenous Zn fertilizer altered the potential functions of microorganisms, which may positively affect bacterial functions.
Variations in polished rice zinc content, grain yield, and soil chemical properties
The Zn content of polished rice and grain yield is shown in Table 2. Compared with the corresponding ZS0, the Zn content of polished rice in the XRZS2, JRZS2, and NRZS2 soil Zn application groups increased by 6.66, 20.05, and 1.40%, respectively, which was consistent with the variation trend of the soil available zinc content (p < 0.05). NJ9108 was the only one that differed significantly (p < 0.05). ZS2 had a relatively higher overall increase in yield. The yield of varieties HHZ, NJ9108, and ND9925 increased by 21.00, 17.34, and 19.52% after Zn treatment, respectively. Two-way ANOVA (Supplementary Table 1) suggested that the Zn content in polished rice and grain yield was significantly affected by cultivar differences and Zn application differences, respectively.
The shift in soil chemical properties among different treatments is provided in Table 4. The soil Zn application increased the soil total nitrogen (TN) by 14.58–34.20% and pH by 1.74–12.20%, but it dramatically decreased the available nitrogen (AN) by 30.56–41.48% and total phosphorus (TP) content (Figure 5A). Cultivar, treatment, and their interactions influenced almost all soil chemical characteristics. However, no effect was observed on total potassium (TK) concentration. Compared with the effect of cultivar, treat affected soil pH significantly, while cultivar had a significant effect on SOC. In addition, the interaction term cultivar × treat had no effect on AN.
Table 4. Two-way ANOVA for soil properties as affected by basal Zn application (treat), rice cultivar (cultivar), and the interaction (treat × cultivar).
Figure 5. (A) Soil properties at the top 15 cm layer in the experimental fields. SOC, soil organic carbon; TN, total nitrogen; AN, available nitrogen; TP, total phosphorus; TK, total potassium; CEC, cation exchange capacity; AZ, available zinc; ZS0, without basal zinc fertilization; ZS2, with basal zinc fertilization. (B) The relation among the bacterial community, environmental parameters, and varieties of rice and soil Zn application. Redundancy analysis (RDA) of the correlation between the bacterial community (OTU level) and soil chemical properties. Red lines represent significant factors, and the pink crosses represent OTUs; (C) correlations between soil chemical properties, grain yield, grain zinc content (GZC), and dominant bacterial phyla (with a total relative abundance > 3%). The relative abundance of dominant bacterial phyla based on Bray–Curtis distance is related to other factors by Mantel’s test. Line width indicates the Mantel’s r statistic, and color denotes the significance based on 999 permutations. The color gradient in heatmap corresponds to Pearson’s correlation coefficient.
The impact of environmental factors on bacterial community
Redundancy analysis (RDA) was conducted to investigate the drivers of change occurring in the microbial community structures (Figure 5B). All the edaphic variables explained up to 57.43% of the variance, with the first axis explaining 37.45% and the second axis explaining another 19.98%. The RDA plots indicate that TN (27.6%, F = 6.1, p = 0.002) resulted in the longest vector followed by SOC, TP, and AZ. Meanwhile, TN had a positive influence on the presence of Acidobacteria and Gemmatimonadetes during the Zn treatment process. In addition, Mantel’s analysis was conducted between dominant phyla, environmental factors, grain yield, and zinc content in polished rice (Figure 5C). The results showed that the phyla Proteobacteria, Acidobacteria, and Gemmatimonadetes were significantly affected by SOC; among them, Proteobacteria and Chloroflexi were also affected by TK and pH as well, which may have potential impact on grain yield. According to Mantel’s p and r statistics, soil SOC affects dominant phyla more than other environmental factors. Overall, TN appeared to be the critical chemical factor for assembling the microbial community structure of the rice rhizosphere, while SOC seemed to have strong influence on dominant phyla.
Discussion
Zn fertilizer applications changed soil properties and bacterial community characteristics
As we all know, the availability of Zn is significantly influenced by soil pH (Zeng et al., 2019), and Zn deficiency is more common in high pH soil (Welch and Graham, 2005). However, in this study, we found higher pH and lower TP after soil Zn application. Xu et al. (2003) reported that submerged conditions increase the pH of acidic soil and decrease the pH of alkaline soil. In addition, Hacisalihoglu (2020) reported the higher solubility of P than Zn, resulting in a shift of Zn from soluble to insoluble Zn (e.g., ZnS). These reports could be the reason for this result, which may also lead to a slight increase in soil available Zn and limit the uptake of zinc in rice. In view of Chen Y. L. et al.’s (2022) report, the available Zn accumulated with the number of years of ZnSO4 application, and the period of this study was also a limiting factor for effective Zn. Prasad et al. (2010) reported that higher Zn increases biomass production, resulting in the addition of large quantities of roots and stubbles, which perhaps build up SOC content. Our result is in line with Prasad. Soil microorganisms served as the major drivers in SOM and nutrient cycling (Rath and Rousk, 2015), whose composition and activity characteristics significantly influence soil nutrient transfer and accumulation (Zheng et al., 2018; Luo et al., 2019). The higher 16S gene copy number might in part explain why SOC content is higher in soil. In addition, SOM has a high CEC (Kleber et al., 2015) which can explain the CEC change trend as SOM. Intriguingly, the content of total nitrogen was opposite to the change of available nitrogen. This may be due to the fact that Zn fertilization reduces the ability of the soil to convert organically bound N to inorganic form.
Furthermore, plant grain nutrient status and grain yield are significantly impacted by nitrogen, including the transportation of Zn in plants (Xue et al., 2019). Gupta et al. (2016) identified the increased root zone translocators and other organic compounds as major drivers that can accelerate the transport and accumulation of Zn in the xylem through N-stimulated activities. The internal relationship between nitrogen and Zn has a remarkable synergistic effect in regulating the Zn concentration in rice grains. This also corresponded to the results that TN was the most critical soil factor during the process. The results of this study demonstrated a lower increase in yield and Zn content for polished rice than the results of previous studies (Khan et al., 2003; Phattarakul et al., 2012; Khampuang et al., 2021). This could be owing to the fact that no NPK fertilizer was used in this study, which has a positive effect on Zn enrichment (Xue et al., 2019).
Previous studies found applying Zn fertilizer did not significantly change the alpha diversity of bacterial community (Garcia-Gomez et al., 2018). However, the Shannon index of HHZ was observed not to change significantly, while that of NJ9108 and ND9925 decreased significantly. It is speculated that excess Zn suppresses or kills the Zn-sensitive species within a short time (e.g., Bacteroidetes), resulting in a decrease of the Sobs index (Zeng et al., 2011). As the accumulated time of Zn fertilizer increases, more microorganisms will become dominant species, thereby lowering the evenness index. After that, competition among the dominant microorganisms may further accelerate the extinction of some bacteria (Zhang et al., 2021), thus decreasing the richness index. The average size of networks (n) decreased after Zn application, in accordance with diversity, but the average links and average connectivity (avgK) increased. This was because microorganisms often form complex networks and change their responses to external disturbances (Jiao et al., 2019; Wu et al., 2022). In addition, species that are highly interconnected and interact frequently are grouped into modules. The modularity of ZS0 is higher than ZS2, one explanation is that Zn application may reduce the significance of the rhizosphere difference in rice (Ling et al., 2022).
Microbial mechanism of Zn fertilizer treatments and rice varieties promoted rice yield and zinc content
Chen Y. L. et al. (2022) reported that Actinobacteria thrived but Bacteroidetes were suppressed after 3 years of Zn application. Our results are exactly in line with it. A high percent abundance of Actinobacteria was observed in ZS2 treatments, which is associated with dissolving P and Zn in the soil, increasing the utilization ratio of P and Zn by plants (Rose et al., 2013). The families of Gaiellaies and the genera of Gaiella belonging to Actinobacteria are significantly enriched in this study, which are considered to have high Zn adaptability and storage potential (Pasti et al., 1990). In contrast, Bacteroidetes as a keystone species, selected by its topological roles, demonstrated suppression in ZS2 treatments, suggesting it is sensitive to Zn application. Based on the keystone species (genus level) selected, we used a correlation heat map based on Pearson’s coefficient to explore the interactions of soil chemical parameters and key genera (Supplementary Figure 5). The results indicated that the abundance of a single keystone had less correlation with the soil chemical parameters and grain yields, while Wang et al. (2021) also obtained similar results; it is speculated that these members of the bacterial community work together to affect the entire soil–plant system rather than acting independently.
The cluster analysis was conducted to explore the differences between rice varieties (Figure 6A). The results demonstrated that all treatments could be roughly divided into two groups. The first group contained XRZS0 and JRZS0 without soil Zn application, and the second group contained three Zn application treatments (ZS2) but NRZS0 without soil Zn application also merged. Overall, the structure of the bacterial community in the rice rhizosphere largely depended on whether soil Zn fertilizer was applied, followed by rice genotypes, which may have further influenced the grain yield and Zn concentration. Wissuwa et al. (2008) showed that rice yield and Zn content in polished rice varied with rice cultivars under the same growth conditions and were closely related to the different Zn absorption and utilization efficiency (White and Broadley, 2005; Song et al., 2015).
Figure 6. (A) Clustering analysis of bacterial communities under three rice cultivars and two soil Zn applications based on OTU abundance. Bray–Curtis similarity coefficients were used; (B) KEGG metabolism pathway (L2 level) differing significantly between ZS0 and ZS2 among three different rice cultivars (HHZ, NJ9108, and ND9925) with an effect size ≥ 10%.
Function predictions among the three tested cultivars (Figure 6B) showed that the metabolism functions of NJ9108 improved better, with significant improvement in lipid metabolism, amino acid metabolism, carbohydrate metabolism, and xenobiotic biodegradation metabolism. Lipid metabolism, as a resistance of microbial cells to environment, may be affected by alterations in gene expression during stress acclimation phase (Kosová et al., 2018). Amino acid metabolism is a primary metabolic process driving microbial metabolism and biosynthesis, which facilitates microbial utilization of amino acid. Amino acids are the main form of organic nitrogen, and its degradation is beneficial for microbial reproduction (Dai et al., 2015; Idrees et al., 2020). Carbohydrate metabolism is related to the biosynthesis of sugars and the decomposition of organic matter, which may be influenced by the different cellulose, xylan, and fructan in plant material of different rice species (Berlemont and Martiny, 2015). Furthermore, membrane transport and cell motility pathway may be related to the input or output of zinc ions and ensure the efficient accumulation and distribution of zinc ions in the cell (Navarrete et al., 2017). These variations in rice cultivars are in good agreement with Chen W. et al. (2022). The downregulated expression in HHZ and ND9925 treatments leads to low nutrient absorption (Kisand et al., 2012), consequently leading to low Zn absorption.
Limitations of the study
The findings of this study have to be seen in light of some limitations as follows:
(1) In order to reduce the application of chemical fertilizer, microbial fertilizers were applied as basal fertilizer. Inoculation of microorganisms into the soil may interfere with the interaction between native microorganisms and exogenous Zn fertilizers. Therefore, the shifts of microbial assemble and interactions among species are not exclusively a function of zinc application.
(2) Due to the outbreak of COVID-19 and time constraints, we only conducted a 1-year field trial. This study offered only a glimpse of the variation of soil microorganism and rice growth with soil Zn application. However, the environment, plants, and soil microorganisms all have variability and plasticity in response to the application of Zn or microbial fertilizer. It is necessary to conduct a long-term field trial to test the results.
(3) Based on the farm records of previous studies those conducted on these fields (blocks) involved in this study, there were no significant differences. However, as we did not apply randomized block design, the effects due to variations among different fields (blocks) could not be completely eliminated. Randomized block design would be a better choice in future studies, to minimize the effects of unpredictable variables on the results.
Conclusion
This study focuses on the effects of soil Zn fertilizers and rice varieties on the grain yield, grain Zn concentrations, bacterial community, and co-occurrence patterns under sustainable farming manipulation. We found that the Zn treatments (ZS2) decreased the soil bacterial diversity, strengthened the network interactions, and enhanced the yield and zinc concentration in grain. Soil total nitrogen (TN) and soil organic carbon (SOC) were identified as the primary drivers that led to shifts in community and dominant phyla, respectively. In addition, the potential functions of microorganisms were improved in ZS2 treatment. This is a novel insight to reveal the bacterial assemble and co-occurrence patterns under sustainable farming manipulation with different soil zinc application levels. Further studies should be conducted to gain a better understanding of key microbial communities and Zn enrichment mechanisms in rice by combining metagenomic and transcriptomic data and finally explore ways to guide these microorganisms in rice farming, which may contribute to Zn-rich rice cultivation strategies and achieve sustainable agriculture.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI – PRJNA845976.
Author contributions
YX: data analysis and writing original draft—revising and editing. BZ: conceptualization, resources, reviewing and editing, supervision, and funding acquisition. ZH, CD, and FJ: data analysis and curation, and reviewing and editing. SL: field experiment, and reviewing and editing. WZ: resources, and reviewing and editing. All authors contributed to the study conception and design, commented on previous versions of the manuscript and approved the final manuscript.
Funding
This research was supported by the National Natural Science Foundation of China (52079139 and 51709261) and the National Training Program of Innovation and Entrepreneurship for Undergraduates (202110019059).
Acknowledgments
We also appreciate the valuable suggestions from John P. Brooks from Genetics and Sustainable Agricultural Research Unit, United States Department of Agriculture.
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/fpls.2022.962246/full#supplementary-material
Abbreviations
TN, soil total nitrogen; AN, soil available nitrogen; TP, soil total phosphorus; TK, soil total potassium; CEC, soil cation exchange capacity; SOC, soil organic carbon; AZ, soil available zinc; PCoA, principal coordinate analysis; RDA, redundancy analysis; LEfSe, linear discriminant analysis effect size; LDA, linear discriminant analysis; OTU, operational taxonomic unit; HHZ, NJ908 and ND9925, Huanghuazhan (Japonica rice), Nanjing-9108 (Japonica rice), and Nuodao99-25 (glutinous rice); ZS0 and ZS2, treatments without added basal zinc fertilizer and with basal zinc fertilizer.
Footnotes
- ^ https://ccb.jhu.edu/software/FLASH/
- ^ http://qiime.org/install/index.html
- ^ http://www.drive5.com/uparse
- ^ https://www.arb-silva.de/
- ^ http://ieg2.ou.edu/MENA
- ^ http://picrust.github.io/picrust/
- ^ https://www.mothur.org/wiki/Download_mothur
References
Afgan, E., Baker, D., Batut, B., van den Beek, M., Bouvier, D., Cech, M., et al. (2018). The galaxy platform for accessible, reproducible and collaborative biomedical analyses. Nucleic Acids Res. 46, W537–W544. doi: 10.1093/nar/gky379
Angelo-Picard, C. D., Faure, D., Carlier, A., Uroz, S., Raffoux, A., Fray, R., et al. (2004). Bacterial populations in the rhizosphere of tobacco plants producing the quorum-sensing signals hexanoyl-homoserine lactone and 3-oxo-hexanoyl-homoserine lactone. FEMS Microbiol. Ecol. 51, 19–29. doi: 10.1016/j.femsec.2004.07.008
Berg, G., and Smalla, K. (2009). Plant species and soil type cooperatively shape the structure and function of microbial communities in the rhizosphere. FEMS Microbiol. Ecol. 68, 1–13. doi: 10.1111/j.1574-6941.2009.00654.x
Berlemont, R., and Martiny, A. C. (2015). Genomic potential for polysaccharide deconstruction in bacteria. Appl. Environ. Microbiol. 81, 1513–1519. doi: 10.1128/AEM.03718-14
Bhatt, K., and Maheshwari, D. K. (2020). Zinc solubilizing bacteria (Bacillus megaterium) with multifarious plant growth promoting activities alleviates growth in Capsicum annuum L. 3 Biotech 10:36. doi: 10.1007/s13205-019-2033-9
Bindraban, P. S., Dimkpa, C. O., and Pandey, R. (2020). Exploring phosphorus fertilizers and fertilization strategies for improved human and environmental health. Biol. Fertil. Soils 56, 299–317. doi: 10.1007/s00374-019-01430-2
Bissett, A., Brown, M. V., Siciliano, S. D., and Thrall, P. H. (2013). Microbial community responses to anthropogenically induced environmental change: Towards a systems approach. Ecol. Lett. 16, 128–139. doi: 10.1111/ele.12109
Broadley, M. R., White, P. J., Hammond, J. P., Zelko, I., and Lux, A. (2007). Zinc in plants. New Phytol. 173, 677–702. doi: 10.1111/j.1469-8137.2007.01996.x
Bruins, M. R., Kapil, S., and Oehme, F. W. (2000). Microbial resistance to metals in the environment. Ecotoxicol. Environ. Saf. 45, 198–207. doi: 10.1006/eesa.1999.1860
Chasapis, C. T., Loutsidou, A. C., Spiliopoulou, C. A., and Stefanidou, M. E. (2012). Zinc and human health: An update. Arch. Toxicol. 86, 521–534. doi: 10.1007/s00204-011-0775-1
Chen, W., Gao, Y., Yang, J., Fan, F. J., Zhang, W. G., Li, J. Y., et al. (2022). Taxonomical and functional bacterial community selection in the rhizosphere of the rice genotypes with different nitrogen use efficiencies. Plant Soil 470, 111–125. doi: 10.1007/s11104-021-05170-x
Chen, Y. L., Mi, H. Z., Zhang, Y. H., Zhang, G. Y., Li, C., Ye, Y., et al. (2022). Impact of ZnSO4 and ZnEDTA applications on wheat Zn biofortification, soil Zn fractions and bacterial community: Significance for public health and agroecological environment. Appl. Soil Ecol. 176:104484. doi: 10.1016/j.apsoil.2022.104484
China National Knowledge Infrastructure (2022). China Economic And Social Big Data Research Platform. Available Online at: https://data.cnki.net/ (accessed July 12, 2022).
Dai, Z. L., Wu, Z. L., Hang, S. Q., Zhu, W. Y., and Wu, G. Y. (2015). Amino acid metabolism in intestinal bacteria and its potential implications for mammalian reproduction. Mol. Hum. Reprod. 21, 389–409. doi: 10.1093/molehr/gav003
Deng, Y., Jiang, Y. H., Yang, Y., He, Z. L., and Luo, F. (2012). Molecular ecological network analyses. BMC Bioinform. 13:113. doi: 10.1186/1471-2105-13-113
Fageria, N. K., Baligar, V. C., and Clark, R. B. (2002). “Micronutrients in crop production,” in Advances In Agronomy, ed. D. L. Sparks (Cambridge: Academic Press), 185–268.
Fageria, N. K., Filho, M. P. B., Moreira, A., and Guimarães, C. M. (2009). Foliar fertilization of crop plants. J. Plant Nutr. 32, 1044–1064. doi: 10.1080/01904160902872826
Garcia-Gomez, G. C., Fernandez, M. D., Garcia, S., Obrador, A. F., Leton, M., and Babin, M. (2018). Soil pH effects on the toxicity of zinc oxide nanoparticles to soil microbial community. Environ. Sci. Pollut. Res. 25, 28140–28152. doi: 10.1007/s11356-018-2833-1
Ghasal, P. C., Shivay, Y. S., Pooniya, V., Choudhary, M., and Verma, R. K. (2018). Zinc partitioning in basmati rice varieties as influenced by Zn fertilization. Crop J. 6, 136–147. doi: 10.1016/j.cj.2017.09.001
Gu, S. S., Hu, Q. L., Cheng, Y. Q., Bai, L. Y., Liu, Z. H., Xiao, W. J., et al. (2019). Application of organic fertilizer improves microbial community diversity and alters microbial network structure in tea (Camellia sinensis) plantation soils. Soil And Tillage Res. 195:104356. doi: 10.1016/j.still.2019.104356
Gupta, N., Ram, H., and Kumar, B. (2016). Mechanism of zinc absorption in plants: Uptake, transport, translocation and accumulation. Rev. Environ. Sci. Biotechnol. 15, 89–109. doi: 10.1007/s11157-016-9390-1
Hacisalihoglu, G. (2020). Zinc (Zn): The Last Nutrient in the Alphabet and shedding light on Zn efficiency for the future of crop production under suboptimal Zn. Plant 9:1417. doi: 10.3390/plants9111471
Hartman, K., van der Heijden, M. G. A., Wittwer, R. A., Samiran, B., Jean-Claude, W., and Klaus, S. (2018). Cropping practices manipulate abundance patterns of root and soil microbiome members paving the way to smart farming. Microbiome 6:14. doi: 10.1186/s40168-017-0389-9
Huang, G. C., Ding, C. F., Zhou, Z. G., Zhang, T. L., and Wang, X. X. (2019). A tillering application of zinc fertilizer based on basal stabilization reduces Cd accumulation in rice (Oryza sativa L.). Ecotoxicol. Environ. Saf. 167, 338–344. doi: 10.1016/j.ecoenv.2018.10.044
Hussain, A., Zahir, Z. A., Ditta, A., Tahir, M. U., Ahmad, M., Mumtaz, M. Z., et al. (2020). Production and implication of bio-activated organic fertilizer enriched with zinc-solubilizing bacteria to boost up Maize (Zea mays L.) production and biofortification under two cropping seasons. Agronomy 10:39. doi: 10.3390/agronomy10010039
Idrees, M., Mohammad, A. R., Karodia, N., and Rahman, A. (2020). Multimodal role of amino acids in microbial control and drug development. Antibiotics 9:330. doi: 10.3390/antibiotics9060330
Jacoby, R., Peukert, M., Succurro, A., Koprivova, A., and Kopriva, S. (2017). The role of soil microorganisms in plant mineral nutrition-current knowledge and future directions. Front. Plant Sci. 8:1617. doi: 10.3389/fpls.2017.01617
Jiao, S., Wang, J. M., Wei, G. H., Chen, W. M., and Lu, Y. H. (2019). Dominant role of abundant rather than rare bacterial taxa in maintaining agro-soil microbiomes under environmental disturbances. Chemosphere 235, 248–259. doi: 10.1016/j.chemosphere.2019.06.174
Kasote, D., Sreenivasulu, N., Acuin, C., and Regina, A. (2021). Enhancing health benefits of milled rice: Current status and future perspectives. Crit. Rev. Food Sci. Nutr. 1–21. doi: 10.1080/10408398.2021.1925629 [Epub ahead of print].
Khampuang, K., Lordkaew, S., Dell, B., and Prom-u-thai, C. (2021). Foliar zinc application improved grain zinc accumulation and bioavailable zinc in unpolished and polished rice. Plant Prod. Sci. 24, 94–102. doi: 10.1080/1343943X.2020.1797512
Khan, M. U., Qasim, M., and Jamil, M. (2002). Response of rice to zinc fertilizer in calcareous soils of D. I. Khan. Asian J. Plant Sci. 1, 1–2. doi: 10.3923/ajps.2002.1.2
Khan, M. U., Qasim, M., Subhan, M., Jamil, M., and Ahmad, R. D. (2003). Response of rice to different methods of zinc application in calcareous soil. J. Appl. Sci. 3, 524–529. doi: 10.3923/jas.2003.524.529
Kisand, V., Valente, A., Lahm, A., Tanet, G., and Lettieri, T. (2012). Phylogenetic and functional metagenomic profiling for assessing microbial biodiversity in environmental monitoring. PLoS One 7:e43630. doi: 10.1371/journal.pone.0043630
Kleber, M., Eusterhues, K., Keiluweit, M., Mikutta, C., Mikutta, R., and Nico, P. (2015). “Chapter one - mineral–organic associations: formation, properties, and relevance in soil environments,” in advances In Agronomy, eds D. L. Sparks (Amsterdam: Elsevier), 1–140. doi: 10.1016/bs.agron.2014.10.005
Kosová, K., Vítámvás, P., Urban, M. O., Prášil, I. T., and Renaut, J. (2018). Plant abiotic stress proteomics: The major factors determining alterations in cellular proteome. Front. Plant Sci. 9:122. doi: 10.3389/fpls.2018.00122
Kumar, A., Dewangan, S., Lawate, P., Bahadur, I., and Prajapati, S. (2019). “Zinc-solubilizing bacteria: a boon for sustainable agriculture,” in Plant Growth Promoting Rhizobacteria for Sustainable Stress Management. Microorganisms for Sustainability, eds R. Sayyed, N. Arora, and M. Reddy (Berlin: Springer), doi: 10.1007/978-981-13-6536-2_8
Langille, M. G. I., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. doi: 10.1038/nbt.2676
Legendre, P., Oksanen, J., and Braak, C. J. F. (2011). Testing the significance of canonical axes in redundancy analysis. Methods Ecol. Evol. 2, 269–277. doi: 10.1111/j.2041-210X.2010.00078.x
Ling, N., Wang, T. T., and Kuzyakov, Y. (2022). Rhizosphere bacteriome structure and functions. Nat. Commun. 13:836. doi: 10.1038/s41467-022-28448-9
Luo, L. Y., Xie, L. L., Jin, D. C., Mi, B. B., Wang, D. H., Li, X. F., et al. (2019). Bacterial community response to cadmium contamination of agricultural paddy soil. Appl. Soil Ecol. 139, 100–106. doi: 10.1016/j.apsoil.2019.03.022
Ma, G., Jin, Y., Li, Y., Zhai, F., Kok, F., Jacobsen, E., et al. (2008). Iron and zinc deficiencies in China: What is a feasible and cost-effective strategy? Public Health Nutr. 11, 632–638. doi: 10.1017/S1368980007001085
Nadeem, F., and Farooq, M. (2019). Application of micronutrients in rice-wheat cropping system of south Asia. Rice Sci. 26, 356–371. doi: 10.1016/j.rsci.2019.02.002
Navarrete, A. A., Mellis, E. V., Escalas, A., Lavres, J., Lemos, L. N., Quaggio, J. A., et al. (2017). Zinc concentration affects the functional groups of microbial communities in sugarcane-cultivated soil. Agric. Ecosyst. Environ. 236, 187–197. doi: 10.1016/j.agee.2016.12.009
Nayar, N. M. (2014). “Chapter 6 - The Origin of Asian Rice,” in Origin and Phylogeny of Rices, ed. N. M. Nayar (San Diego: Academic Press), 169–253.
Nazir, Q., Wang, X., Hussain, A., Ditta, A., Aimen, A., Saleem, I., et al. (2021). Variation in Growth, Physiology, Yield, and Quality of Wheat under the Application of Different Zinc Coated Formulations. Appl. Sci. 11:4797. doi: 10.3390/app11114797
Pablo, S. A., Ana, N., and Luís, G. (2014). Soil biochemical alterations and microbial community responses under Acacia dealbata Link invasion. Soil Biol. Biochem. 79, 100–108. doi: 10.1016/j.soilbio.2014.09.008
Parks, D. H., Tyson, G. W., Hugenholtz, P., and Beiko, R. G. (2014). STAMP: Statistical analysis of taxonomic and functional profiles. Bioinformatics 30, 3123–3124. doi: 10.1093/bioinformatics/btu494
Pasti, M. B., Pometto, A. L., Nuti, M. P., and Crawford, D. L. (1990). Lignin-solubilizing ability of actinomycetes isolated from termite (termitidae) gut. Appl. Environ. Microbiol. 56, 2213–2218. doi: 10.1128/aem.56.7.2213-2218.1990
Phattarakul, N., Rerkasem, B., Li, L. J., Wu, L. H., Zou, C. Q., Ram, H., et al. (2012). Biofortification of rice grain with zinc through zinc fertilization in different countries. Plant Soil 361, 131–141. doi: 10.1007/s11104-012-1211-x
Prasad, R. K., Kumar, V., Prasad, B., and Singh, A. P. (2010). Long-term effect of crop residues and zinc fertilizer on crop yield, nutrient uptake and fertility build-up under rice-wheat cropping system in Calciorthents. J. Indian Soc. Soil Sci. 58, 205–211.
Qin, Y., Melse-Boonstra, A., Yuan, B., Pan, X., Dai, Y., Zhou, M., et al. (2012). Zinc biofortification of rice in China: A simulation of zinc intake with different dietary patterns. Nutrients 4, 517–528. doi: 10.3390/nu4060517
Rahman, M., Jahangir, M. M., Kibria, M. G., Hossain, M., Hosenuzzaman, M., Solaiman, Z. M., et al. (2022). Determination of critical limit of zinc for rice (Oryza sativa L.) and Potato (Solanum tuberosum L.) Cultivation In Floodplain Soils Of Bangladesh 14:167. doi: 10.3390/su14010167
Rath, K. M., and Rousk, J. (2015). Salt effects on the soil microbial decomposer community and their role in organic carbon cycling: A review. Soil Biol. Biochem. 81, 108–123. doi: 10.1016/j.soilbio.2014.11.001
Rehman, A., Farooq, M., Ozturk, L., Muhammad, A., and Kadambot, H. M. S. (2018). Zinc nutrition in wheat-based cropping systems. Plant Soil 422, 283–315. doi: 10.1007/s11104-017-3507-3
Rose, T. J., Impa, S. M., Rose, M. T., Pariasca-Tanaka, J., Mori, A., Heuer, S., et al. (2013). Enhancing phosphorus and zinc acquisition efficiency in rice: A critical review of root traits and their potential utility in rice breeding. Ann. Bot. 112, 331–345. doi: 10.1093/aob/mcs217
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537–7541. doi: 10.1128/AEM.01541-09
Song, W. Y., Lee, H. S., Jin, S. R., Ko, D., Martinoia, E., Lee, Y., et al. (2015). Rice PCR1 influences grain weight and Zn accumulation in grains. Plant Cell Environ. 38, 2327–2339. doi: 10.1111/pce.12553
Stackebrandt, E., and Goebel, B. M. (1994). Taxonomic note: a place for DNA-DNA reassociation and 16S rRNA sequence analysis in the present species definition in bacteriology. Microbiol. Soc. 44, 846–849. doi: 10.1099/00207713-44-4-846
Wang, J. L., Liu, K. L., Zhao, X. Q., Zhang, H. Q., Li, D., Li, J. J., et al. (2021). Balanced fertilization over four decades has sustained soil microbial communities and improved soil fertility and rice productivity in red paddy soil. Sci. Total Environ. 793:148664. doi: 10.1016/j.scitotenv.2021.148664
Wang, S. S., Lay, S., Yu, H. N., and Shen, S. R. (2016). Dietary guidelines for Chinese residents (2016): Comments and comparisons. J. Zhejiang Univ. Sci. 17, 649–656. doi: 10.1631/jzus.B1600341
Wang, Y., Randolph, T. W., Shojaie, A., and Ma, J. (2019). The generalized matrix decomposition biplot and its application to microbiome data. mSystems 4:e00504–19. doi: 10.1128/mSystems.00504-19
Weese, D. J., Heath, K. D., Dentinger, B. T. M., and Lau, J. A. (2015). Long-term nitrogen addition causes the evolution of less-cooperative mutualists. Evolution 69, 631–642. doi: 10.1111/evo.12594
Welch, R. M., and Graham, R. D. (2005). Agriculture: The real nexus for enhancing bioavailable micronutrients in food crops. J. Trace Elem. Med. Biol. 18, 299–307. doi: 10.1016/j.jtemb.2005.03.001
White, P. J., and Broadley, M. R. (2005). Biofortifying crops with essential mineral elements. Trends Plant Sci. 10, 586–593. doi: 10.1016/j.tplants.2005.10.001
Wissuwa, M., Ismail, A. M., and Graham, R. D. (2008). Rice grain zinc concentrations as affected by genotype, native soil-zinc availability, and zinc fertilization. Plant Soil 306, 37–48. doi: 10.1007/s11104-007-9368-4
Wu, C. C., Ma, Y. J., Wang, D., Shan, Y. P., Song, X. P., Hu, H. Y., et al. (2022). Integrated microbiology and metabolomics analysis reveal plastic mulch film residue affects soil microorganisms and their metabolic functions. J. Hazard. Mater. 432:127258. doi: 10.1016/j.jhazmat.2021.127258
Xu, R. K., Zhao, A. Z., Li, Q. M., Kong, X. L., and Ji, G. L. (2003). Acidity regime of the Red Soils in a subtropical region of southern China under field conditions. Geoderma 115, 75–84. doi: 10.1016/S0016-7061(03)00077-6
Xue, Y. F., Yue, S. C., Liu, D. Y., Zhang, W., Chen, X. P., and Zou, C. Q. (2019). Dynamic zinc accumulation and contributions of pre- and/or post-silking zinc uptake to grain zinc of maize as affected by nitrogen supply. Front. Plant Sci. 10:1203. doi: 10.3389/fpls.2019.01203
Yamaji, N., Xia, J. X., Mitani-Ueno, N., Yokosho, K., and Ma, J. F. (2013). Preferential delivery of zinc to developing tissues in rice is mediated by P-type heavy metal ATPase OsHMA2. Plant Physiol. 162, 927–933. doi: 10.1104/pp.113.216564
Zeb, H., Hussain, A., Naveed, M., Ditta, A., Ahmad, S., Jamshaid, M. U., et al. (2018). Compost enriched with ZnO and Zn-solubilising bacteria improves yield and Zn-fortification in flooded rice. Ital. J. Agron. 13, 310–316. doi: 10.4081/ija.2018.1295
Zeng, F. R., Ali, S., Zhang, H. T., Ouyang, Y. N., Qiu, B. Y., Wu, F. B., et al. (2011). The influence of pH and organic matter content in paddy soil on heavy metal availability and their uptake by rice plants. Environ. Pollut. 159, 84–91. doi: 10.1016/j.envpol.2010.09.019
Zeng, P., Guo, Z. H., Xiao, X. Y., and Peng, C. (2019). Dynamic response of enzymatic activity and microbial community structure in metal(loid)-contaminated soil with tree-herb intercropping. Geoderma 345, 5–16. doi: 10.1016/j.geoderma.2019.03.013
Zhang, X. J., Zhan, Y. B., Zhang, H., Wang, R. H., Tao, X. L., Zhang, L. P., et al. (2021). Inoculation of phosphate-solubilizing bacteria (Bacillus) regulates microbial interaction to improve phosphorus fractions mobilization during kitchen waste composting. Bioresour. Technol. 340:125714. doi: 10.1016/j.biortech.2021.125714
Zheng, H., Wang, X., Chen, L., Wang, Z. Y., Xia, Y., Zhang, Y. P., et al. (2018). Enhanced growth of halophyte plants in biochar-amended coastal soil: Roles of nutrient availability and rhizosphere microbial modulation. Plant Cell Environ. 41, 517–532. doi: 10.1111/pce.12944
Keywords: rice, zinc application, bacterial community, co-occurrence network, microbial function
Citation: Xiao YS, Zhou B, Han Z, Liu S, Ding C, Jia F and Zeng W (2022) Microbial mechanism of zinc fertilizer input on rice grain yield and zinc content of polished rice. Front. Plant Sci. 13:962246. doi: 10.3389/fpls.2022.962246
Received: 06 June 2022; Accepted: 29 July 2022;
Published: 25 August 2022.
Edited by:
Bo Zhang, Xinjiang Institute of Ecology and Geography (CAS), ChinaReviewed by:
Xiaofan Na, Lanzhou University, ChinaAllah Ditta, Shaheed Benazir Bhutto University, Pakistan
Copyright © 2022 Xiao, Zhou, Han, Liu, Ding, Jia and Zeng. 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: Bo Zhou, zhoubo89@cau.edu.cn; Wenzhi Zeng, zengwenzhi1989@whu.edu.cn