- 1Gansu Provincial Key Laboratory of Aridland Crop Science, Lanzhou, China
- 2College of Agronomy, Gansu Agricultural University, Lanzhou, China
- 3College of Resource and Environment, Gansu Agricultural University, Lanzhou, China
- 4State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing, China
Soil autotrophic bacteria play a crucial role in regulating CO2 fixation and crop productivity. However, the information is limited to how fertilization amendments alter soil autotrophic bacterial community, crop yield, and carbon emission efficiency (CEE). Here, we estimated the impact of the structure and co-occurrence network of soil autotrophic bacterial community on maize yield and CEE. A long-term field experiment was conducted with five fertilization treatments in semiarid Loess Plateau, including no amendment (NA), chemical fertilizer (CF), chemical fertilizer plus commercial organic fertilizer (SC), commercial organic fertilizer (SM), and maize straw (MS). The results showed that fertilization amendments impacted the structure and network of soil Calvin–Benson–Bassham (CBB) (cbbL) gene-carrying bacterial community via changing soil pH and NO3–N. Compared with no amendment, the cbbL-carrying bacterial diversity was increased under the SC, SM, and MS treatments but decreased under the CF treatment. Soil autotrophic bacterial network contained distinct microbial modules that consisted of closely associated microbial species. We detected the higher abundances of soil cbbL-carrying bacterial genus Xanthobacter, Bradyrhizobium, and Nitrosospira. Structural equation modeling further suggested that the diversity, composition, and network of autotrophic bacterial community had strongly positive relationships with CEE and maize yield. Taken together, our results suggest that soil autotrophic bacterial community may drive crop productivity and CEE, and mitigate the atmospheric greenhouse effect.
Introduction
The growing emissions of greenhouse gasses, carbon dioxide (CO2), nitrous oxide (N2O), and methane (CH4), contributing to global warming are of great concern worldwide (Mora et al., 2018). Soil respiration is the second largest carbon flux between the terrestrial ecosystem and the atmosphere. Carbon dioxide released into the atmosphere through soil respiration has a significant effect on the global C budget (Peylin et al., 2013) and global warming (Cherubini et al., 2011). Annually, CO2 emission from soil accounts for 25% of the global total C exchange between the atmosphere and terrestrial sources (Schlesinger and Andrews, 2000), and agroecosystems represent a significant source of CO2 emission. Agronomic management practices such as tillage and fertilization have a marked effect on soil ecosystem and soil CO2 emission (Ding et al., 2007; Dong et al., 2009). Plants assimilate atmospheric CO2 through photosynthesis, increase carbon use efficiency, and reduce the atmospheric CO2 concentration (Eric et al., 2018). Therefore, reducing CO2 emissions and increasing carbon assimilation under agricultural practices are effective strategies for mitigating atmospheric greenhouse effect in the context of global climate change.
Soil microorganisms play key roles in mediating soil organic matter and nutrient cycling dynamics and profoundly impact soil productivity and sustainability (Theuerl and Buscot, 2010; Nimmo et al., 2013). Soil autotrophic bacteria play an essential role in soil C fixation through the Calvin–Benson–Bassham (CBB) cycle to utilize atmosphere CO2 (Selesi et al., 2005). Ribulose-1,5-bisphosphate carboxylase/oxygenase (RubisCO) is responsible for catalyzing the first rate-limiting step in the CBB cycle for autotrophic CO2 fixation. To date, four RubisCO forms (forms I–IV) have been found that differ in structure, catalytic property, and O2 sensitivity. Form I of RubisCO is the most abundant among the four forms. The cbbL gene, which encodes a large subunit of RubisCO I, has been often used as a phylogenetic marker to investigate the autotrophic bacterial community (Kovaleva et al., 2011). Previous studies have reported that fertilization treatments significantly affect the cbbL-carrying bacterial community structure and, consequently, facilitate carbon sequestration (Yuan et al., 2012b). So far, the question remains as to whether and how fertilization treatments alter the structure and network of the soil cbbL-carrying bacterial community.
Fertilization treatments can significantly affect nutrient-supplying (Yousaf et al., 2017), soil quality, and grain yield (Steiner et al., 2007; Alburquerque et al., 2013). Long-term reliance on chemical fertilizers leads to a decline in soil quality (Liu et al., 2010), nutrient imbalance (Mandal et al., 2007), and reduced fertilizer use efficiency (Duan et al., 2014). Organic fertilizer plus chemical fertilizers can increase soil organic C and enhance soil fertility (Mandal et al., 2007; Diacono and Montemurro, 2011). Furthermore, fertilization enhances CO2 fixation, carbon emission efficiency (CEE), and crop yield via promoting soil autotrophic bacterial activity (Selesi et al., 2005; Yuan et al., 2012b; Souza et al., 2015). However, information on the direction and extent of autotrophic bacteria that govern soil CEE and crop yield under different fertilization treatments is limited.
On the Chinese semiarid Loess Plateau, maize is an important crop, where low rainfall (Xin and Lu, 2011) and low heat units (Gan et al., 2013) challenge crop yield. The fully plastic mulched ridge-furrow system has become a widely adopted practice for producing high crop yields and water use efficiency in recent years (Rong et al., 2013).
The intent of our study was to investigate the impact of the structure and network of soil autotrophic bacterial community on CEE and maize yield. We performed a long-term field experiment with five fertilization treatments on the semiarid Loess Plateau. Specifically, the objectives of this study were to (i) investigate the diversity, composition, and co-occurrence network of soil autotrophic bacterial community under fertilization treatments; and (ii) explore the underlying mechanisms of soil properties and the autotrophic bacterial community in regulating C emission, CEE, and maize yield. We hypothesized that organic manure plus chemical fertilizers increased soil autotrophic bacterial diversity and changed composition and network of soil autotrophic bacterial community. Furthermore, we expected that soil properties and autotrophic bacterial community jointly regulated carbon emission, CEE, and crop yield.
Materials and Methods
Field Experiment Description
The long-term fertilization experiment commenced at the Rainfed Agricultural Experimental Station (35°28′N, 104°44′E, elevation 1,971 m above sea level) of Gansu Agricultural University in Gansu province, China. The experimental site is at the temperature semiarid zone in the western Loess Plateau, with a long-term average solar radiation of 593 kJ cm–2 and an annual average sunshine of 2,477 h. The mean air temperature is 6.4°C, with accumulated temperature > 0°C of 2,934°C and > 10°C of 2,239°C. The aeolian soil at the experimental site is locally known as Huangmian (Gong, 2001), which is a Calcaric Cambisol according to the Food and Agriculture Organization (FAO) soil classification (FAO, 1990).
The long-term field experiment followed a completely randomized design with three replications. The experiment was conducted in 2012, which consist of 15 plots with the following treatments: no amendment (NA) and four fertilization treatments with the same nitrogen (N) (200 kg N ha–1) and phosphate (150 kg P2O5 ha–1) input: (1) chemical fertilizer (CF), applied 200 kg N ha–1 as urea (46% N) plus 150 kg P2O5 ha–1 triple superphosphate (16% P2O5); (2) 50% chemical N fertilizer plus 50% organic N fertilizer (SC), applied 3.03 t ha–1 of commercial organic fertilizer plus 100 kg N ha–1 as urea and 120 kg P2O5 ha–1 as triple superphosphate; (3) organic fertilizer (SM), applied 6.06 t ha–1 of commercial organic fertilizer plus 90 kg P2O5 ha–1 as triple superphosphate; and (4) maize straw (MS), applied at 28.5 t ha–1 plus 36 kg ha–1 as triple superphosphate. The commercial organic fertilizer used contained 3.3% N, 1.0% P2O5, and 0.7% K2O, with cow manure as ripened organic fertilizer (Gansu Daxing Agricultural Technology Co., Gansu, China). The maize straw obtained from the previous harvesting season was air-dried, shredded to ≤ 5 cm, weighed, and then applied to the field plots. The maize straw contained 0.7% N, 0.4% P2O5, and 0.5% K2O. Representative samples of maize straw were collected at the time of application and analyzed for nutrient concentration using an Elementar Vario MACRO Cube (Elementar, Hanau, Germany). All fertilizers were evenly applied on the soil surface and incorporated by moldboard plowing to a depth of 25 cm in the spring.
The experimental plots were 42.9 m2 (13 m long and 3.3 m wide) and consisted of narrow ridges (15 cm high × 40 cm wide) alternated with wide ridges (10 cm high × 70 cm wide) (Supplementary Figure 1). Plastic film (0.01 mm thick and 140 cm wide; Lanzhou Green Garden Corp., Lanzhou, China) covered all ridges to increase soil temperature and reduce evaporative losses. After the plastic film is placed, a handheld device made holes through the film’s furrows to channel rainfall from the ridges into the root zone. Maize (cv. Pioneer 335) was sown in late April and harvested in mid-October with a plant spacing in each row at 35 cm with a density of 52,500 plants ha–1. Weeds were manually removed during the maize growing season. All agronomic practices, except for the fertilization amendments, were kept uniform in this study.
Soil Sampling and Soil Property Assays
At silking stage of maize in early August 2019, five randomly selected soil cores (5 cm in diameter) to 20-cm depth from each plot were passed through a 2-mm sieve to remove residues and mixed as one composite sample. A subsample from each composite sample was immediately placed with dry ice for transport to the laboratory and stored at –80°C for DNA analysis, while the remaining soil was air-dried for soil chemical analysis.
Soil pH was measured in distilled water (soil:water ratio of 1:2.5 m/v) using a pH meter (Mettler Toledo FE20, Shanghai, China). A modified Walkley–Black wet oxidation method and the Kjeldahl method were used to determine soil organic carbon (SOC) and total nitrogen (TN), respectively (Bao, 2000). Soil NO3–N and NH4–N were extracted with 2 M of KCl and measured using a continuous flow analyzer (Skalar, Breda, Netherlands). Available phosphorus (AP) was extracted with sodium bicarbonate and measured using the molybdenum-blue method (Olsen, 1954).
Soil DNA Extraction and Amplification
Total genome DNA from soil samples was extracted using an OMEGA soil DNA Kit (Omega Bio-Tek, Doraville, GA, United States). The quantity and quality of DNA were determined by agarose gel (1%) electrophoresis and a NanoDrop 2000 (Thermo Fisher Scientific, Waltham, MA, United States). According to the concentration, DNA was diluted to 1 ng μl–1 using sterile water. The specific primer K2f/V2r, ACCAYCAAGCCSAAGCTSGG/GCCTTCSAGCTTGCCSACCRC was used to amplify the autotrophic bacterial community (Tolli and King, 2005). All PCR were carried out in 30-μl reactions with 15 μl of PCR Master Mix (New England Biolabs, Ipswich, MA, United States), 0.2 μM of forward and reverse primers, and 10 ng of template DNA. Thermal cycling consisted of initial denaturation at 98°C for 1 min, followed by 30 cycles of denaturation at 98°C for 10 s, annealing at 50°C for 30 s, and elongation at 72°C for 60 s. The PCR products remained at 72°C for 5 min, and then the same volume of 1 × loading buffer was mixed with the PCR products and analyzed by electrophoresis on 2% agarose gel. The PCR products were then mixed in equivalent ratios and purified with a GeneJET Gel Extraction Kit (Thermo Fisher Scientific, Waltham, MA, United States).
Illumina MiSeq Sequencing
The library was prepared using the NEB Next Ultra II DNA library prep kit for Illumina (New England Biolabs, Ipswich, MA, United States) according to manufacturer’s instructions. The library quality was assessed using a Qubit@ 2.0 Fluorometer (Thermo Fisher Scientific, Ipswich, MA, United States) and the Agilent Bioanalyzer 2100 system (Agilent Tec., Palo Alto, CA, United States). The library was sequenced on an Illumina MiSeq platform (Shanghai Biozeron Biotechnology Co., Ltd, Pudong, Shanghai, China), generating 250-bp paired-end reads. Paired-end reads from the original bacterial DNA fragments were merged using FLASH (version 1.2.7)1. Sequence analyses were performed with the UPARSE (version 7.1)2 software package using the UPARSE-OTU and UPARSE-OTUref algorithms (Edgar, 2013), and the operational taxonomic unit (OTU) was clustered at a distance of ≤ 0.03 (approximately 97% of the sequence similarity). Alpha-diversity of the autotrophic bacterial community was calculated after rarefying all samples to an equal number of 24,138 reads. Alpha-diversity analysis was performed using MOTHUR software (Schloss et al., 2009), and diversity indices were determined, including the Chao 1, Shannon, and Simpson indices (Caporaso et al., 2010). We have deposited sequencing reads in the Sequence Read Archive at National Center for Biotechnology Information (NCBI) under the accession number PRJNA700109.
Soil Respiration, Total Carbon Emission, Grain Yield, and Carbon Emission Efficiency
In order to minimize the root effect, the plastic tube was pressed in the big ridge with a distance of 25 cm far from maize plants. Polyvinyl chloride collars that were 20 cm in diameter and 11.5 cm in height were pressed into the soil post-crop emergence to a depth of 9.5 cm, leaving 2 cm exposed above the soil surface and under the plastic mulch. These collars were used to locate the LI-8100 soil chamber (LI-COR Inc., Lincoln, NE, United States) and provided an airtight seal. Soil respiration was measured at 09:00, 11 times per year from May 3 to October 16 in 2018 and from May 5 to October 17 in 2019. Twelve hours before the soil measurements, the plastic mulch was removed from the collars to release emissions under the mulch. After the sampling, the plastic mulch was replaced.
Carbon emission (kg ha–1) was calculated based on soil respiration (μmol m–2 s–1) using the following equation (Zhai et al., 2011):
where Rs is soil respiration (μmol CO2 m–2 s–1) measured at 15-day intervals during the growing season; i + 1 and i are the current and previous sampling dates, respectively; t is days after sowing; 0.1584 converts μmol CO2 m–2 s–1 to g CO2 m–2 h–1; 0.2727 converts g CO2 m–2 h–1 to g C m–2 h–1; and 24 and 10 convert g C m–2 h–1 to kg ha–1 for the growing season.
Plots were harvested at physiological maturity of maize, and yield was determined on an oven-dry-weight basis by drying the grain at 105°C for 45 min and subsequently to constant weight at 85°C (Lamptey et al., 2017).
CEE (kg kg–1) was calculated from crop yield per unit of C emission using the following equation (Qin et al., 2013):
Data Analysis
Analysis of variance was performed by Fisher’s protected least significant (P < 0.05) difference test using SPSS software (v. 20.0, IBM Corp., Chicago, IL, United States). Pearson’s correlation coefficients were used to assess linear associations among soil properties, communities of soil cbbL-carrying bacteria, crop yield, total carbon emission (TCE), and CEE across soil amendment. All statistical analyses were conducted based on 15 samples (5 treatments × 3 replicates). Principal coordinate analysis (PCoA) was used to evaluate the Bray–Curtis distances of soil cbbL-carrying bacterial community compositions using the “vegan” package in R (version 4.0.2) software.
To describe the complex co-occurrence patterns in soil cbbL-carrying bacterial networks, co-occurrence networks were constructed using Spearman’s correlation and the Kullback–Leibler dissimilarity (KLD) measures (Faust and Raes, 2012). The OTUs detected in more than two-thirds of the soil samples were kept for network construction. A valid co-occurrence was considered a statistically robust correlation between species when the correlation coefficient (r) was > 0.6 or less than -0.6 and the P value was < 0.05. The co-occurrence network visualization was conducted using Gephi software (Bastian et al., 2009), and modules were defined as clusters of closely interconnected nodes (i.e., groups of co-occurring microbes) (Layeghifard et al., 2017).
Random forest tool was used to quantitatively evaluate the important predictors of crop yield and CEE. Random forest modeling was conducted using the RandomForest package (Bento et al., 2002), and the rfPermute package (Archer, 2016) was used to determine the model significance and predictor importance. The significant predictors were further chosen to perform a structural equation modeling (SEM) analysis. SEM analysis was applied to determine the contributions of soil properties and soil cbbL-carrying bacterial community to the maize yield and CEE. SEM analysis was conducted via the robust maximum likelihood evaluation method using AMOS 22.0 (SPSS, Chicago, IL, United States). Latent variables were used to integrate the effects of multiple conceptually related observed variables into a single-composite effect, aiding interpretation of model results. The SEM fitness was examined on the basis of a nonsignificant chi-square test (P > 0.05), the goodness-of-fit index, and the root mean square error of approximation (Hooper et al., 2008).
Results
Soil Properties, C Emission, Crop Yield, and Carbon Emission Efficiency
Soil chemical properties were significantly affected by soil amendment (P < 0.05, Table 1). Soil TN, NO3–N, and AP were significantly (P < 0.05) increased under four fertilization treatments compared with the NA treatment, while soil pH exhibited a significant decreasing trend (P < 0.05). SOC under the SC, SM, and MS treatments was significantly greater than that under NA and CF treatments (P < 0.01). However, there was no significant difference in NH4–N among the five treatments (P = 0.454).
During the maize growing seasons, soil respiration exhibited a similar pattern across five treatments in 2018 and 2019 (Figure 1). Soil respiration gradually increased from sowing in May and reached a maximum in July and then decreased to a minimum at harvest in October. Compared with the NA treatment, the MS treatment significantly increased in mean soil respiration and C emission by 63.6% and 56.7% in 2018 and by 69.5% and 68.8% in 2019, followed by the CF, SC, and SM treatments (Table 2). Crop yield and CEE significantly increased under CF and SC treatments compared with SM, MS, and NA treatments in 2018 and 2019 (P < 0.01). However, TCE was not significantly different under CF, SC, and SM treatments (Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-g001.jpg)
Figure 1. Soil respiration under five treatments during the 2018 and 2019 growing seasons. Vertical bars indicate the least significant difference value at P < 0.05. NA, no amendment; CF, chemical fertilizer; SC, 50% chemical N fertilizer plus 50% organic N fertilizer; SM, organic fertilizer; MS, maize straw.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-t002.jpg)
Table 2. Mean soil respiration (Rs), total carbon emission (TCE), crop yield, and carbon emission efficiency (CEE) in 2018 and 2019 under five treatments.
Community Structure of Soil Autotrophic Bacteria
With the use of Illumina sequencing of cbbL-carrying bacteria amplicons, a total of 362,070 sequences were obtained after quality control. The alpha-diversity of the autotrophic bacterial community indicated by Shannon index and Chao 1 richness was significantly (P < 0.05) higher under the MS treatment than under the CF, NA, SC, and SM treatments (Supplementary Table 1). Across all samples, the autotrophic bacterial community consisted mainly of Alphaproteobacteria (26.7%), Betaproteobacteria (11.2%), and Actinobacteria (9.6%) (Figure 2). The dominant community closely matched facultative chemoautotrophic hydrogen oxidizers within the genera Xanthobacter (11.1%), Bradyrhizobium (6.6%), and Nitrosospira (6.4%) (Supplementary Figure 2). The results of PCoA, based on the Bray–Curtis distance, showed that soil autotrophic bacterial community structure significantly differed among five treatments (Figure 3). The PCoA1 and PCoA2 explained 40.8% and 8.4%, respectively, of the variation in soil autotrophic bacterial community composition. The composition of soil autotrophic bacterial community under SM and SC treatments was clustered together with the NA treatment, whereas those communities under the CF and MS treatments were clearly separated from the NA treatment. The diversity of soil autotrophic bacteria was significantly correlated with pH (r = 0.56, P < 0.05), TN (r = 0.66, P < 0.01), and SOC (r = 0.80, P < 0.01), while the community composition was significantly correlated with SOC (r = 0.74, P < 0.01) (Table 3). The diversity and composition of soil autotrophic bacterial community had significantly positive relationships with TCE (r = 0.75, P < 0.01 and r = 0.72, P < 0.01) but negative correlations with CEE (r = -0.66, P < 0.01 and r = -0.62, P < 0.01) and crop yield (r = -0.61, P < 0.05 and r = -0.53, P < 0.05) (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-g002.jpg)
Figure 2. Relative abundance of the autotrophic bacterial phyla under five treatments. NA, no amendment; CF, chemical fertilizer; SC, 50% chemical N fertilizer plus 50% organic N fertilizer; SM, organic fertilizer; MS, maize straw. The numbers before the treatment name indicate the sampling replications. For example, 1-NA, 2-NA, and 3-NA mean the sampling was taken from replicates 1, 2, and 3 of the field plots, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-g003.jpg)
Figure 3. The structure of soil autotrophic bacterial community by principal coordinate analysis (PCoA), which is constrained by five treatments based on Bray–Curtis distances. NA, no amendment; CF, chemical fertilizer; SC, 50% chemical N fertilizer plus 50% organic N fertilizer; SM, organic fertilizer; MS, maize straw.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-t003.jpg)
Table 3. Correlation coefficients between soil properties, the diversity and composition of soil autotrophic bacterial community, crop yield, total carbon emission (TCE), and carbon emission efficiency (CEE).
The Autotrophic Bacterial Co-occurrence Networks
To examine the co-occurrence pattern of soil cbbL-carrying bacterial communities in five treatments, the network was constructed from the 15 samples (Supplementary Figure 3). In total, there were 533 nodes and 4,526 links and three dominant modules (modules I, II, and III) in the autotrophic bacterial network (Supplementary Figure 3). The autotrophic bacterial network displayed more positive correlations (4,472 edges) than negative correlations (54 edges). Modules I, II, and III of the bacterial networks consisted of 122, 205, and 206 nodes, respectively. At the phylum level, the relative abundances of Alphaproteobacteria and Gemmatimonadetes were significantly greater in module I than in modules II and III; those of Betaproteobacteria, Deltaproteobacteria, Planctomycetes, Verrucomicrobia, and Acidobacteria were significantly greater in module III than in modules I and II; and those of Actinobacteria and Gammaproteobacteria were significantly greater in modules II and III than in module I (Figure 4A, P < 0.05). At the genus level, the relative abundances of Xanthobacter (Alphaproteobacteria), Gemmatirosa (Gemmatimonadetes), Nitrosospira (Betaproteobacteria), Mycobacterium (Actinobacteria), and Pirellula (Planctomycetes) followed a similar trend of individual phyla (Figure 4B, P < 0.05). Module I was positively correlated with NO3–N (r = 0.75, P < 0.01), CEE (r = 0.95, P < 0.001), and crop yield (r = 0.79, P < 0.01) (Table 3). Modules II and III were positively associated with TN (r = 0.59, P < 0.05 and r = 0.59, P < 0.05), SOC (r = 0.64, P < 0.05 and r = 0.78, P < 0.01), and TCE (r = 0.74, P < 0.01 and r = 0.83, P < 0.01) (Table 3).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-g004.jpg)
Figure 4. Relative abundance of different dominant modules in soil autotrophic bacterial networks at phylum (A) and genera (B) levels. Different small letters indicate the significant difference among modules at P < 0.05. Bars represent standard errors (n = 15).
Soil Properties and Soil Autotrophic Bacterial Community Affected Crop Yield and Carbon Emission Efficiency
Random forest modeling showed that soil pH (8.6% and 7.9%, P < 0.05) and NO3–N (9.4% and 8.8%, P < 0.01) were the two primary predictors of crop yield and CEE (Supplementary Figure 4). Module I (8.2%, P < 0.01; and 9.7%, P < 0.01) of the bacterial networks exhibited more pronounced effects on maize yield and CEE than diversity (2.5%, P > 0.05; and 6.5%, P < 0.05) and composition (1.2%, P > 0.05; and 6.6%, P < 0.05) of the autotrophic bacterial community (Supplementary Figure 4). SEM further suggested that soil properties were significantly correlated with the autotrophic bacterial community, CEE, and maize yield. The autotrophic bacterial community had strong positive relationships with CEE and crop yield according to diversity, composition, and module I of the bacterial network (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/664120/fmicb-12-664120-HTML/image_m/fmicb-12-664120-g005.jpg)
Figure 5. Impacts of soil properties and soil autotrophic bacterial community on crop yield and carbon emission efficiency (CEE) as estimated using structural equation modeling (SEM) analysis. Based on random forest analyses, the significant predictors were chosen to perform the SEM analysis. The latent variables inside the ellipse were used to integrate the effects of multiple conceptually related observed variables into a single-composite effect. Soil properties are represented by pH and nitrate nitrogen (NO3–N). Soil autotrophic bacterial community (cbbL) is represented by the diversity (Shannon index), composition [first principal coordinates (PCoA1)], and two modules (module II and module III) of the soil cbbL-carrying bacteria. *P < 0.05; **P < 0.01; ***P < 0.001.
Discussion
Fertilization Treatments Affected Soil Autotrophic Bacterial Community
Based on the cbbL gene-carrying bacteria analysis, our results showed that organic fertilization treatments increased the diversity of soil autotrophic bacteria compared with the NA and CF treatments. Observations of large autotrophic bacterial diversity were presumably derived from changes in soil properties as affected by organic treatments. Organic fertilizer with abundant SOC and available nutrients may provide spatially adaptive living micro-habitat for the autotrophic bacterial community and thus promoted the growth of different microbial taxa and increased soil microbial diversity (Yuan et al., 2012b; Francioli et al., 2016). In addition, highest autotrophic bacterial diversity with MS suggested that maize straw are easily degradable organic materials, which provided sufficient substrate for microorganisms and increased soil microbial diversity (Yuan et al., 2012b). Soil pH has been reported as a good predictor that plays crucial roles in driving the microbial community structure and functions (Fierer and Jackson, 2006; da Jesus et al., 2009; Jiang et al., 2016). The decomposition of SOC can provide an adequate supply of nutrients and energy and significantly changed the autotrophic bacterial communities (Liao et al., 2020).
Fertilization treatments differentially modulated the autotrophic bacterial community structure with distinct shifts in bacterial diversity. The cbbL-carrying bacterial communities were mainly assigned to Xanthobacter, Bradyrhizobium, Nitrosospira, Aminobacter, Mycobacterium, Nocardia, Gemmatirosa, and Pirellula, etc. Most of these bacteria are previously reported as facultative autotrophs (Selesi et al., 2005; Tolli and King, 2005; Yuan et al., 2015), which play significant roles in CO2 fixation and degrade several carbon-containing organic compounds (Guo et al., 2015; Liao et al., 2020). However, we also found that some genera may not be autotrophs, such as the genus Pirellula and Gemmatirosa (DeBruyn et al., 2013; Glöckner et al., 2003). We further revealed that the autotrophic bacterial network contained distinct microbial modules that consisted of closely associated microbial species. The high ratio of positive edges to negative edges in the modules may suggest the ecological interactions (i.e., commensalism and mutualism) among soil autotrophic microorganisms (Zhang et al., 2018), and these associations in the microbial networks may also represent niche sharing among microorganisms (Berry and Widder, 2014). The modules in co-occurrence network generally consist of highly interconnected microorganisms and carry out various functions in maintaining ecosystem processes of nutrient exchange and resource availability (de Menezes et al., 2015; Jiang et al., 2016). Modules II and III significantly related to SOC and TN suggested that the interconnected taxa in two modules may share the same habitat preferences and common ecological niche. These modules in the bacterial network can organize a higher-order structured meta-module to reflect modules’ functional relevant interactions (Jiang et al., 2015). However, it is worth noting that caution is warranted in referring to bacterial interactions from network patterns because correlations do not necessarily represent causal relationship.
The Mechanism of Autotrophic Bacterial Community Driving Carbon Emission Efficiency and Crop Yield
The variations in soil autotrophic bacterial community offer significant potentials for CO2 fixation. A study based on 14CO2 labeling has recently estimated that microbial CO2 assimilation processes were annually responsible for roughly 4% of the total CO2 fixation in terrestrial ecosystems (Yuan et al., 2012a). The present study suggested that autotrophic bacterial diversity was positively correlated with TCE and negatively correlated with CEE and plant productivity. Long-term fertilization study showed that soil autotrophic CO2 fixation processes were more important for the nutrient-poor environment; chemical fertilization-driving changes in soil autotrophic microbial community depressed soil CO2 fixation (Liao et al., 2020). In the present study, long-term organic fertilizer increased organic matter and soil autotrophic bacterial diversity but decreased CEE and crop yield as compared with chemical fertilizer treatment. These results suggested that nutrient-rich organic fertilizer treatments showed inhibitory effects on the autotrophic metabolic pathway and may consequently decrease CEE and crop yield. Theoretical and empirical evidence supported that microbial diversity plays a crucial role in predicting multiple ecological functions in terrestrial ecosystems, including nutrient cycling and climate regulation (Bardgett and Van Der Putten, 2014; Delgado-Baquerizo et al., 2016; Domeignoz-Horta et al., 2020). However, it is still a matter of debate how and to what extent soil autotrophic bacterial diversity drives microbial carbon metabolism. We speculated that the autotrophic bacterial diversity reduced soil RubisCO and CO2 fixation activities of bacterial populations, thereby declining CEE and crop production. High diversity suppresses microbial community function primarily due to the potential relationships among taxa in the modules of the autotrophic bacterial networks (Chen et al., 2019; Saleem et al., 2019). Highly similar microbial species shared with overlapping microhabitats (niches) may induce minimal stimulation in microbial C metabolic functioning as diversity increases (Maynard et al., 2017). In this case, reduced diversity of the autotrophic bacterial community is probably beneficial to CEE and crop yield improvement. As such, our findings highlight that soil autotrophic microorganisms may play key roles in driving carbon assimilation and contributing to soil C pool and crop productivity. However, many autotrophic microorganisms are facultative autotrophic in nature, and they have both heterotrophic and autotrophic metabolic pathways (Miltner et al., 2005). Therefore, soil autotrophic microorganisms and CO2 fixation only have a certain degree of correlation. Their causal relationships in the soil environment need further study. Despite that the structure and network of autotrophic bacterial community have been considered, we cannot rule out the possible impact of the cbbL copy numbers. Future research by stable isotope-based evidence could help gain novel insights regarding how active autotrophs interact with each other in the autotrophic bacterial community.
Conclusion
In the present study, long-term fertilization significantly changed soil autotrophic bacterial diversity, composition, and co-occurrence network through soil pH and NO3–N. Soil autotrophic bacterial community had significant relationships with crop yield and CEE. Taken together, our results suggest that soil autotrophic bacterial community may drive crop productivity and CEE and mitigate greenhouse effect. Further research should continue to take advantage of pot incubation to verify the causal relationships between soil autotrophic bacterial community, crop yield, and CEE. If confirmed, we will provide novel insights into the importance of soil autotrophic bacterial community on CEE and crop yield in agricultural systems.
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 at: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA700109.
Author Contributions
LL, JX, and YJ designed the experiments. JW, LW, and YJ curated the literature data. JW, RZ, and YJ elaborated the data, ran the statistical analysis, and arranged the tables and figures. JW wrote the first draft of the manuscript. LL, JX, YJ, LW, and ZL critically revised the manuscript. All authors read and approved the final manuscript.
Funding
This work was supported by the National Natural Science Foundation of China (31761143004, 31660373, and 41922048), Fostering Foundation for the Excellent Ph.D. Dissertation of Gansu Agricultural University (YB2020002), Department of Science and Technology of Gansu Province (GSPT-2018-56), and Outstanding Youth Scholar Program of Jiangsu Province (BK20180049).
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.
Acknowledgments
We appreciate the assistance in the field and laboratory by students of the Rainfed Agricultural Experimental Station of Gansu Agricultural University.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2021.664120/full#supplementary-material
Footnotes
References
Alburquerque, J. A., Salazar, P. F., Barron, V., Torrent, J., Campillo, M. C. D., Gallardo, A., et al. (2013). Enhanced wheat yield by biochar addition under different mineral fertilization levels. Agron. Sustain. Dev. 33, 475–484.
Archer, E. (2016). rfPermute: Estimate Permutation p-Values for Random Forest Importance Metrics R package version 2.18.
Bardgett, R. D., and Van Der Putten, W. H. (2014). Belowground biodiversity and ecosystem functioning. Nature 515, 505–511. doi: 10.1038/nature13855
Bastian, M., Heymann, S., and Jacomy, M. (2009). “Gephi: An Open Source Software for Exploring and Manipulating Networks,” in Proceedings of the Third International Conference on Weblogs and Social Media, ICWSM 2009, (San Jose, CA), 17–20.
Bento, A. P., Gaulton, A., Hersey, A., Bellis, L. J., Chambers, J., Davies, M., et al. (2002). Classification and Regression by randomForest. R News 23.
Berry, D., and Widder, S. (2014). Deciphering microbial interactions and detecting keystone species with co-occurrence networks. Front. Microbiol. 5:219.
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335. doi: 10.1038/NMETH.F.303
Chen, L., Jiang, Y., Liang, C., Luo, Y., Xu, Q., Han, C., et al. (2019). Competitive interaction with keystone taxa induced negative priming under biochar amendments. Microbiome 7, 1–18. doi: 10.1186/s40168-019-0693-7
Cherubini, F., Peters, G. P., Berntsen, T., Strømman, A. H., and Hertwich, E. (2011). CO2 emissions from biomass combustion for bioenergy: atmospheric decay and contribution to global warming. GCB Bioenergy 3, 413–426. doi: 10.1111/j.1757-1707.2011.01102.x
da Jesus, E., Marsh, T. L., Tiedje, J. M., and de Moreira, F. M. (2009). Changes in land use alter the structure of bacterial communities in Western Amazon soils. ISME J. 3, 1004–1011. doi: 10.1038/ismej.2009.47
de Menezes, A. B., Prendergast-Miller, M. T., Richardson, A. E., Toscas, P., Farrell, M., Macdonald, L. M., et al. (2015). Network analysis reveals that bacteria and fungi form modules that correlate independently with soil parameters. Environ. Microbiol 17, 2677–2689. doi: 10.1111/1462-2920.12559
DeBruyn, J. M., Fawaz, M. N., Peacock, A. D., Dunlap, J. R., Nixon, L. T., Cooper, K. E., et al. (2013). Gemmatirosa kalamazoonesis gen. nov., sp. nov., a member of the rarely-cultivated bacterial phylum Gemmatimonadetes. J. Gen. Appl. Microbiol. 59, 305–312. doi: 10.2323/jgam.59.305
Delgado-Baquerizo, M., Maestre, F. T., Reich, P. B., Jeffries, T. C., Gaitan, J. J., Encinar, D., et al. (2016). Microbial diversity drives multifunctionality in terrestrial ecosystems. Nat. Commun 7, 1–8.
Diacono, M., and Montemurro, F. (2011). Long-term effects of organic amendments on soil fertility. Sustainable Agriculture 2, 761–786. doi: 10.1007/978-94-007-0394-0_34
Ding, W., Meng, L., Yin, Y., Cai, Z., and Zheng, X. (2007). CO2 emission in an intensively cultivated loam as affected by long-term application of organic manure and nitrogen fertilizer. Soil Boil. Biochem. 39, 669–679. doi: 10.1016/j.soilbio.2006.09.024
Domeignoz-Horta, L. A., Pold, G., Liu, X.-J. A., Frey, S. D., Melillo, J. M., and DeAngelis, K. M. (2020). Microbial diversity drives carbon use efficiency in a model soil. Nat. Commun. 11, 1–10. doi: 10.1038/s41467-020-17502-z
Dong, W., Hu, C., Chen, S., and Zhang, Y. (2009). Tillage and residue management effects on soil carbon and CO2 emission in a wheat–corn double-cropping system. Nutr. Cycl. Agroecosyst. 83, 27–37. doi: 10.1007/s10705-008-9195-x
Duan, Y., Xu, M., Gao, S., Yang, X., Huang, S., Liu, H., et al. (2014). Nitrogen use efficiency in a wheat–corn cropping system from 15 years of manure and fertilizer applications. Field Crop Res. 157, 47–56. doi: 10.1016/j.fcr.2013.12.012
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996. doi: 10.1038/nmeth.2604
Eric, D. M., Galvao, D. A., and Way, D. A. (2018). Plant carbon metabolism and climate change: elevated CO2 and temperature impacts on photosynthesis, photorespiration and respiration. New Phytol. 221, 32–49. doi: 10.1111/nph.15283
FAO (1990). Soil map of the world. Revised legend. World Resources Report 60. Italy: Food and Agriculture Organization of the United Nations Rome.
Faust, K., and Raes, J. (2012). Microbial interactions: from networks to models[J]. Nat Rev Microbiol. 10, 538–550.. doi: 10.1038/nrmicro2832
Fierer, N., and Jackson, R. B. (2006). The diversity and biogeography of soil bacterial communities. Proceedings of the National Academy of Sciences 103, 626–631. doi: 10.1073/pnas.0507535103
Francioli, D., Schulz, E., Lentendu, G., Wubet, T., Buscot, F., and Reitz, T. (2016). Mineral vs. organic amendments: microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilization strategies. Front. Microbiol. 7:1446. doi: 10.3389/fmicb.2016.01446
Gan, Y., Siddique, K. H. M., Turner, N. C., Li, X. G., Niu, J. Y., Yang, C., et al. (2013). Chapter Seven–Ridge-Furrow Mulching Systems—An Innovative Technique for Boosting Crop Productivity in Semiarid Rain-Fed Environments. Adv. Agron. 118, 429–476. doi: 10.1016/B978-0-12-405942-9.00007-4
Glöckner, F. O., Kube, M., Bauer, M., Teeling, H., Lombardot, T., Ludwig, W., et al. (2003). Complete genome sequence of the marine planctomycete Pirellula sp. strain 1. Proceedings of the National Academy of Sciences 100, 8298–8303.
Guo, G., Kong, W., Liu, J., Zhao, J., Du, H., Zhang, X., et al. (2015). Diversity and distribution of autotrophic microbial community along environmental gradients in grassland soils on the Tibetan Plateau. Appl. Microniol. Biot. 99, 8765–8776. doi: 10.1007/s00253-015-6723-x
Hooper, D., Coughlan, J., and Mullen, M. R. (2008). Structural Equation Modelling: Guidelines for Determining Model Fit. Elec. J. Buss. Res. Methods 6, 141–146.
Jiang, Y., Liang, Y., Li, C., Wang, F., Sui, Y., Suvannang, N., et al. (2016). Crop rotations alter bacterial and fungal diversity in paddy soils across East Asia. Soil Boil. Biochem. 95, 250–261. doi: 10.1016/j.soilbio.2016.01.007
Jiang, Y., Sun, B., Li, H., Liu, M., Chen, L., and Zhou, S. (2015). Aggregate-related changes in network patterns of nematodes and ammonia oxidizers in an acidic soil. Soil Boil. Biochem. 88, 101–109. doi: 10.1016/j.soilbio.2015.05.013
Kovaleva, O. L., Tourova, T. P., Muyzer, G., Kolganova, T. V., and Sorokin, D. Y. (2011). Diversity of RuBisCO and ATP citrate lyase genes in soda lake sediments. FEMS Microbiol. Ecol. 75, 37–47. doi: 10.1111/j.1574-6941.2010.00996.x
Lamptey, S., Li, L., Xie, J., Zhang, R., Yeboah, S., and Antille, D. L. (2017). Photosynthetic response of maize to nitrogen fertilization in the semiarid Western Loess Plateau of China. Crop Sci. 57, 2739–2752. doi: 10.2135/cropsci2016.12.1021
Layeghifard, M., Hwang, D. M., and Guttman, D. S. (2017). Disentangling Interactions in the Microbiome: A Network Perspective. Trends Microbiol 25, 217–228. doi: 10.1016/j.tim.2016.11.008
Liao, H., Qin, F., Wang, K., Zhang, Y., Hao, X., Chen, W., et al. (2020). Long-term chemical fertilization-driving changes in soil autotrophic microbial community depresses soil CO2 fixation in a Mollisol. Sci. Total. Environ 748, 141317. doi: 10.1016/j.scitotenv.2020.141317
Liu, E., Yan, C., Mei, X., He, W., Bing, S. H., Ding, L., et al. (2010). Long-term effect of chemical fertilizer, straw, and manure on soil chemical and biological properties in northwest China. Geoderma 158, 173–180. doi: 10.1016/j.geoderma.2010.04.029
Mandal, A., Patra, A. K., Singh, D., Swarup, A., and Masto, R. E. (2007). Effect of long-term application of manure and fertilizer on biological and biochemical activities in soil during crop development stages. Bioresour Technol 98, 3585–3592. doi: 10.1016/j.biortech.2006.11.027
Maynard, D. S., Crowther, T. W., and Bradford, M. A. (2017). Competitive network determines the direction of the diversity–function relationship. Proceedings of the National Academy of Sciences 114, 11464–11469. doi: 10.1073/pnas.1712211114
Miltner, A., Kopinke, F. D., Kindler, R., Selesi, D., Hartmann, A., and Kästner, M. (2005). Non-phototrophic CO2 fixation by soil microorganisms. Plant and Soil 269, 193–203. doi: 10.1007/s11104-004-0483-1
Mora, C., Spirandelli, D., Franklin, E. C., Lynham, J., Kantar, M. B., Miles, W., et al. (2018). Broad threat to humanity from cumulative climate hazards intensified by greenhouse gas emissions. Nat. Clim. Chang. 8, 1062–1075. doi: 10.1038/s41558-018-0315-6
Nimmo, J., Lynch, D. H., and Owen, J. (2013). Quantification of nitrogen inputs from biological nitrogen fixation to whole farm nitrogen budgets of two dairy farms in Atlantic Canada. Nutr. Cycl. Agroecosyst. 96, 93–105. doi: 10.1007/s10705-013-9579-4
Olsen, S. R. (1954). Estimation of available phosphorus in soils by extraction with sodium bicarbonate. Washington: United States Department of Agriculture.
Peylin, P., Law, R., Gurney, K., Chevallier, F., Jacobson, A., Maki, T., et al. (2013). Global atmospheric carbon budget: results from an ensemble of atmospheric CO2 inversions. Biogeosciences 10, 6699–6720. doi: 10.5194/bg-10-6699-2013
Qin, A., Huang, G., Chai, Q., Yu, A., and Huang, P. (2013). Grain yield and soil respiratory response to intercropping systems on arid land. Field Crop. Res 144, 1–10. doi: 10.1016/j.fcr.2012.12.005
Rong, L., Hou, X., Jia, Z., Han, Q., Ren, X., and Yang, B. (2013). Effects on soil temperature, moisture, and maize yield of cultivation with ridge and furrow mulching in the rainfed area of the Loess Plateau, China. Agric. Water Manage 116, 101–109. doi: 10.1016/j.agwat.2012.10.001
Saleem, M., Hu, J., and Jousset, A. (2019). More than the sum of its parts: microbiome biodiversity as a driver of plant growth and soil health. Annu. Rev. Ecol. Evol. Syst. 50, 145–168. doi: 10.1146/annurev-ecolsys-110617-062605
Schlesinger, W. H., and Andrews, J. A. (2000). Soil respiration and the global carbon cycle. Biogeochemistry 48, 7–20. doi: 10.1023/A:1006247623877
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
Selesi, D., Schmid, M., and Hartmann, A. (2005). Diversity of green-like and red-like ribulose-1, 5-bisphosphate carboxylase/oxygenase large-subunit genes (cbbL) in differently managed agricultural soils. Appl. Environ. Microbiol. 71, 175–184. doi: 10.1128/AEM.71.1.175-184.2005
Souza, R. D., Ambrosini, A., and Passaglia, L. M. (2015). Plant growth-promoting bacteria as inoculants in agricultural soils. Genet. Mol. Biol 38, 401–419. doi: 10.1590/S1415-475738420150053
Steiner, C., Teixeira, W. G., Lehmann, J., Nehls, T., Macêdo, J. L. V., Blum, W. E. H., et al. (2007).). Long term effects of manure, charcoal and mineral fertilization on crop production and fertility on a highly weathered Central Amazonian upland soil. Plant Soil 291, 275–290. doi: 10.1007/s11104-007-9193-9
Theuerl, S., and Buscot, F. (2010). Laccases: toward disentangling their diversity and functions in relation to soil organic matter cycling. Biol. Fertil. Soils 46, 215–225. doi: 10.1007/s00374-010-0440-5
Tolli, J. D., and King, G. M. (2005). Diversity and structure of bacterial chemolithotrophic communities in pine forest and agroecosystem soils. Appl. Environ. Microbiol. 71, 8411–8418. doi: 10.1128/AEM.71.12.8411-8418.2005
Xin, Z., and Lu, X. X. (2011). Spatiotemporal variation in rainfall erosivity on the Chinese Loess Plateau during the period 1956–2008. Reg. Envir. Chang 11, 149–159. doi: 10.1007/s10113-010-0127-3
Yousaf, M., Li, J., Lu, J., Ren, T., Cong, R., Fahad, S., et al. (2017). Effects of fertilization on crop production and nutrient-supplying capacity under rice-oilseed rape rotation system. Sci. Rep. 7, 1270. doi: 10.1038/s41598-017-01412-0
Yuan, H., Ge, T., Chen, C., O’Donnell, A. G., and Wu, J. (2012a). Significant role for microbial autotrophy in the sequestration of soil carbon. Appl. Environ. Microbiol. 78, 2328–2336. doi: 10.1128/AEM.06881-11
Yuan, H., Ge, T., Chen, X., Liu, S., Zhu, Z., Wu, X., et al. (2015). Abundance and diversity of CO2-assimilating bacteria and algae within red agricultural soils are modulated by changing management practice. Microb. Eco. 70, 971–980. doi: 10.1007/s00248-015-0621-8
Yuan, H., Ge, T., Wu, X., Liu, S., Tong, C., Qin, H., et al. (2012b). Long-term field fertilization alters the diversity of autotrophic bacteria based on the ribulose-1,5-biphosphate carboxylase/oxygenase (RubisCO) large-subunit genes in paddy soil. Appl. Microbiol. Biotechnol. 95, 1061–1071. doi: 10.1007/s00253-011-3760-y
Zhai, L., Liu, H., Zhang, J., Huang, J., and Wang, B. (2011). Long-Term Application of Organic Manure and Mineral Fertilizer on N2O and CO2 Emissions in a red soil from cultivated maize-wheat rotation in China. Agr. Sci. China 10, 1748–1757. doi: 10.1016/s1671-2927(11)60174-0
Keywords: fertilization amendment, soil autotrophic bacteria, carbon emission efficiency, crop yield, semiarid Loess Plateau
Citation: Wang J, Xie J, Li L, Luo Z, Zhang R, Wang L and Jiang Y (2021) The Impact of Fertilizer Amendments on Soil Autotrophic Bacteria and Carbon Emissions in Maize Field on the Semiarid Loess Plateau. Front. Microbiol. 12:664120. doi: 10.3389/fmicb.2021.664120
Received: 04 February 2021; Accepted: 10 May 2021;
Published: 18 June 2021.
Edited by:
Yurong Liu, Huazhong Agricultural University, ChinaReviewed by:
Jiajun Hu, Shanghai University, ChinaHua Qin, Zhejiang Agriculture and Forestry University, China
Hangwei Hu, The University of Melbourne, Australia
Copyright © 2021 Wang, Xie, Li, Luo, Zhang, Wang and Jiang. 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: Lingling Li, lill@gsau.edu.cn; Yuji Jiang, yjjiang@issas.ac.cn