Jun Chen1,2,3
Yasir Arafat1,2,4
Israr Ud Din5
Bo Yang1,2,3
Liuting Zhou1,2
Juanying Wang1,2
Puleng Letuma2,3
Hongmiao Wu1,2,3
Xianjin Qin2,3
Linkun Wu1,2,3
Sheng Lin1,2
Zhixing Zhang1,2,3*
Wenxiong Lin1,2,3*- 1College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou, China
- 2Key Laboratory of Crop Ecology and Molecular Physiology, Fujian Agriculture and Forestry University, Fuzhou, China
- 3Fujian Provincial Key Laboratory of Agroecological Processing and Safety Monitoring, Fujian Agriculture and Forestry University, Fuzhou, China
- 4Key Laboratory of Crop Genetic Breeding and Comprehensive Utilization of the Ministry of Education, Fujian Agriculture and Forestry University, Fuzhou, China
- 5Institute of Biotechnology and Genetic Engineering, The University of Agriculture Peshawar, Peshawar, Pakistan
Availability of nitrogen (N) in soil changes the composition and activities of microbial community, which is critical for the processing of soil organic matter and health of crop plants. Inappropriate application of N fertilizer can alter the rhizosphere microbial community and disturb the soil N homeostasis. The goal of this study was to assess the effect of different ratio of N fertilizer at various early to late growth stages of rice, while keeping the total N supply constant on rice growth performance, microbial community structure, and soil protein expression in rice rhizosphere. Two different N regimes were applied, i.e., traditional N application (NT) consists of three sessions including 60, 30 and 10% at pre-transplanting, tillering and panicle initiation stages, respectively, while efficient N application (NF) comprises of four sessions, i.e., 30, 30, 30, and 10%), where the fourth session was extended to anthesis stage. Soil metaproteomics combined with Terminal Restriction Fragment Length Polymorphism (T-RFLP) were used to determine the rhizosphere biological process. Under NF application, soil enzymes, nitrogen utilization efficiency and rice yield were significantly higher compared to NT application. T-RFLP and qPCR analysis revealed differences in rice rhizosphere bacterial diversity and structure. NF significantly decreased the specific microbes related to denitrification, but opposite result was observed for bacteria associated with nitrification. Furthermore, soil metaproteomics analysis showed that 88.28% of the soil proteins were derived from microbes, 5.74% from plants, and 6.25% from fauna. Specifically, most of the identified microbial proteins were involved in carbohydrate, amino acid and protein metabolisms. Our experiments revealed that NF positively regulates the functioning of the rhizosphere ecosystem and further enabled us to put new insight into microbial communities and soil protein expression in rice rhizosphere.
Introduction
Rice is the staple food for more than half of the world’s population (Xing and Zhang, 2010). Global demand for food is predicted to grow by 70% in next few decades (FAO, 2009), which may put enormous pressure on already threatened agricultural system. In China, the scope for expansion of the irrigated rice area is limited. Researchers have succeeded to increase the crop yield by developing improved rice varieties through molecular breeding, but the efforts are continued to further enhance the yield traits of rice genotypes in order to meet the growing demand for food (Kato et al., 2007; Xing and Zhang, 2010). However, modern cultivation practices such as sustainable fertilizers management may help to boost the potential of existing varieties to get higher yield.
Nitrogen, which is typically applied in the form of urea, is the most essential element that determines the yield of rice (Witt et al., 1999). Therefore, for improving Nitrogen Use Efficiency (NUE) of crops and reducing off-field losses, it is pertinent to understand the molecular mechanisms accountable for N homeostasis in intensively N fertilized agriculture system. Denitrifying bacteria expressing narG or napA, where the most representative one is narG (Smith et al., 2007; Reyna et al., 2010), encode for nitrate reductase (NR) and convert to , which is then reduced to NO by another bacteria harboring nirK or nirS gene (Priemé et al., 2002; Yan et al., 2003). Nitric oxide reductase encoded by the norB gene reduces NO to N2O, and in the final step nosZ encoding nitrous oxide reductase reduces N2O to N2 (Chon et al., 2011). Estimation and monitoring of the N2O or N2 emission in rice paddy field is technically difficult, however, expression analysis of the catabolic genes regulating biological denitrification and production of N2O and N2 could be useful biomarkers for the determination of N loss in the system (Morales et al., 2010; Bouskill et al., 2012).
In search of minimizing the environmental pollution and enhancing NUE, researchers have regulated the timing and dosage of N application according to plant demand to facilitate maximum cycling of plant available N and to ensure sustainable rice production (Iwasaki et al., 1992; Ramasamy et al., 1997; Zhang et al., 2014). Our previous study showed that appropriately extending N application has the advantages such as delaying leaf senescence, enhancing photosynthetic rate at middle-to-late stage of grain filling, and stable synthesis and transportation of assimilates (Zhang et al., 2011). However, there is no reports hitherto regarding proteomics focusing on the role of rice fertilizer regimes inflicting upon the soil microbial communities.
Microbes residing in the rhizosphere play key roles in nutrient acquisition, nitrogen cycling, and carbon cycling (Nannipieri et al., 2003; Cardon and Whitbeck, 2007). The diverse soil nutrients and plant secondary metabolites present in the exudates enhance the enrichment of specific functional or taxonomic bacterial groups residing in the rhizosphere (Bais et al., 2006; Zahar et al., 2008). Currently, various DNA-dependent strategies, such as denaturing gradient gel electrophoresis (DGGE) (Muyzer and Smalla, 1998; Chen et al., 2017), terminal restriction fragment length polymorphism (T-RFLP) (Osborn et al., 2000; Chen et al., 2018), and 16S rRNA gene sequence have been utilized to determine the biological diversity and composition of microbial communities in soil ecosystems (Kim et al., 2012; Arafat et al., 2017). However, comprehensive characterization of biological processes needs integrated approaches (Griffin et al., 2002). Metaproteomics is a recognized tool that is used to determine soil ecological processes and the environmental factors that affect the functioning of the rhizosphere soil ecosystem (Bastida et al., 2009). This approach has been used to elucidate the functioning of the microbial community, particularly, their role in ecosystem services, such as agriculture (Lin et al., 2013; Zampieri et al., 2016; Bastida et al., 2017; Starke et al., 2017), bioremediation (Williams et al., 2010; Guazzaroni et al., 2013; Bastida et al., 2016a), C cycling (Schneider et al., 2012; Chourey et al., 2013; Bastida et al., 2015a, b, 2016b; Starke et al., 2016), and climate change factors (i.e., temperature) (Hultman et al., 2015; Bastida et al., 2017; Liu et al., 2017). Both nucleic acids and proteomics studies have shown good correlation, at least at the phylogenetic level (Hultman et al., 2015; Bastida et al., 2016b, 2017). Therefore, a metaproteomics approach, to investigate ecosystem functioning in rhizosphere is significant because it provides both phylogenetic and functional information (Bastida et al., 2016b; 2017; Starke et al., 2016, 2017).
Studying the effect of N amendment on soil properties, rhizosphere microbial activities and their correlation with crop performance is therefore crucial for understanding the impact of sustainable N fertilization. In this study, we investigated the effects of N (in the form of urea) application in four sessions (NF) on (i) soil N availability, uptake and percent nitrogen use efficiency (% NUE), enzymatic activity and crop productivity; (ii) the associated microbial diversity using T-RFLP; and (iii) the protein profile of the NT and NF rhizospheric soils at late growing stage of rice through metaproteomics.
Materials and Methods
Plant Material and Cultivation
The large-panicle rice cultivar Jinhui 809 (indica) was used in this study. The experiment was carried out at the experimental station of Fujian Agriculture and Forestry University, Fuzhou, China in 2012 and 2013 (119.280 E, 26.080 N). Seedlings were transplanted at the 5-leaf stage (one seedling per hill with a spacing of 0.15 m × 0.15 m). Each plot measured 4 m × 4 m and received N (225 kg ha–1), P2O5 (112.5 kg ha–1), and K2O (180 kg ha–1) fertilizer. The experiment was designed as a complete randomized block (CRB design) with three replicates. Phosphorus was applied as a basal dressing, and the potassium was used as the top dressing. Two different application ratios of nitrogen were used as treatments. NT treatment was applied in three sessions, i.e., 60% before transplanting, 30% at tillering, and 10% at panicle initiation. NF treatment was applied in four sessions, i.e., 30% before transplanting, 30% at tillering, 30% at panicle initiation, and 10% at anthesis (Supplementary Figure S1). The physical and chemical properties of the soil were investigated before the experiment was initiated. Soil was sandy loam, pH value of 6.2, and total nitrogen of 2.20 g kg–1, available nitrogen of 40.6 mg kg–1, total phosphorus of 0.65 g kg–1, available phosphorus of 26.6 mg kg–1, total K of 1.05 g kg–1 and available K of 30.16 mg kg–1. The annual temperature and average yearly precipitation were 25–32°C and 900–1362 mm, respectively.
Plant Physiological Data and Nitrogen Utilization
Chlorophyll content was determined at seven growth stages of rice using SPAD 502 (Minolta Camera, Co., Osaka, Japan) chlorophyll meter, i.e., the early booting stage, the late booting stage, and five other periods that were measured every 7 days from the full panicle stage to maturity. At booting stage, the three uppermost fully expanded leaves were selected from each treatments. Three chlorophyll meter readings were taken around the midpoint of each leaf blade on one side of the midrib. At full panicle stage, SPAD measurements were made on the flag leaves of each treatments.
Rice plants were harvested, subjected to 105°C for 30 min to deactivate enzymatic activity and then dried at 80°C for 48 h until a constant weight was reached. The plant samples were ground and digested with H2SO4. Grain production efficiency (NUEg), Nitrogen dry matter production efficiency (NUEb) and the nitrogen physiological efficiency (NPE) were calculated as the following equation: NUEg = YR/TN, NUEb = TD/TN, NPE = Δ (YR-Y0)/Δ (TN-TN0), respectively. Among them, YR represents yield of rice with nitrogen; TN represents the sum of the nitrogen contents in each plant; TD represents total plant dry matter; Y0 represents yield of rice without nitrogen; TN0 represents the sum of the nitrogen contents in each plant with nitrogen (Bremner and Mulvaney, 1982).
Grain and Yield Parameters of Rice
In 2012 and 2013, number of grains per panicle, seed setting rate, 1000 grain weight, and grain yield of each treatment were recorded.
Soil Sampling of the Rice Rhizosphere
Soil samples were collected (three replications) 14 days after anthesis (in which the strong grain filling reaches peak while weak grain filling just begun, Supplementary Figure S1), from NT and NF treated soil for nutrient (NPK), enzymatic, microbial community, and metaproteomic analyses. Soil samples were collected after digging the plant samples, removing the loosely attached soils and then scraping the soil that was still attached to the rice. Fresh samples were used for soil enzyme and nutrition analysis. Other samples were stored at −80°C for soil microbial community and metaproteomics analyses (Yang and Zhang, 2009; Zhang et al., 2012).
Determination of Soil Nutrients (NPK) and Soil Enzymatic Activities
Soil urease [EC 3.5.1.5] activity was measured by incubating 5 g fresh soil with 30 ml of extracting solution at 37°C for 24 h. The formation of ammonium was determined spectrophotometrically at 578 nm (Wang et al., 2009). Soil invertase [EC 3.2.1.26] activity was measured by incubating 5 g fresh soil with 15 ml of 8% sucrose solution at 37°C for 24 h. The suspension reacted with 3, 5-dinitrosalicylic acid and absorbance was measured at 508 nm (Wang et al., 2009). Nitrate reductase [EC 1.7.1.3] activity of the soil was determined according to the method already described (Abdelmagid and Tabatabai, 1987).
The available and total amounts of principle nutrient components (NPK) of rice plants were measured using the methods described by Jackson and Barak (Jackson, 2005). All determinations were performed in triplicate, and means were assessed for significant differences (P ≤ 0.05) using SPSS 20.0.
Soil DNA Extraction and T-RFLP Analysis
Soil (1 g) from each treatment was used to extract DNA by using the BioFast Soil Genomic DNA Extraction Kit (Hangzhou, China). Extracted DNA was subsequently stored at −20°C until further use. The forward primer 27F-FAM (5′-AGAGTTTGATCCTGGCTCAG-3′) and reverse primer 1492R (5′-GGTTACCTTGTTACGACTT-3′) were used for 16S rRNA amplification. All PCR products were visualized with 1.2% agarose gel electrophoresis and were purified using the Gel Extraction Kit (OMEGA Bio-Tek, United States) according to the manufacturer’s instructions. Purified PCR products were digested separately with four enzymes (HaeIII, MspI, AluI, and AfaI), and the digested products were sent to Shanghai Sheng Gong for analysis (ABI automated sequencer analyzer; Model 3130 Applied Bio systems).
T-RFLP profiles were analyzed using Gene Marker software (Version 1.2). The terminal fragments between 30 and 600 bp were selected for further analysis. Affiliations of the fragments were determined via online T-RFLP analysis of the Ribosomal Database Project II (RDP II1). The relative abundance (Pi) of terminal restriction fragments (TRFs), Shannon’s diversity index (H), the Shannon-Wiener index (H’), Simpson’s diversity Index (D) and Pielou’ index were calculated according to the following formulas using EXCEL2013:Pi = n/N, H = −Σ Pi log (Pi), H′ = −Σ Pi ln (Pi), D = 1–Σ Pi2, E = H′/ln(S).
Where n is the peak area that can be recognized from the TRFs segments; N is the total peak area of the TRFs segments; and S is the number of peak areas that can be recognized from the TRFs segments.
Quantification of Functional Communities Involved in Denitrification
Real-time PCR quantification of genes encoding the key enzymes of nitrate reduction (narG encoding membrane-bound nitrate reductase) and denitrification (nirK encoding cd1 and copper nitrite reductase) were used to estimate the density of functional communities involved in nitrogen cycling by using primers and conditions previously described (Hallin et al., 2009). In brief, PCR was performed in 15 μl reaction containing 7.5 μl 2 × SYBR green I SuperReal Premix (TransGen Biotech, Beijing, China), 0.5 μl of each primer (10 μM) and template DNA (20 ng of total soil DNA or a serial dilution of plasmid DNA for standard curves).
Protein Extraction and 2D-PAGE
The soil proteins from two samples were extracted and purified using the following protocol developed in our laboratory (Wang et al., 2010). The protein concentration was determined using the Bradford assay. For the 2D-polyacrylamide gel electrophoresis, immobilized pH 4–7 gradient (IPG) precast gels (24 cm in length) were purchased from GE Healthcare (Bjellqvist et al., 1982). For protein separation, a 1300-μg soil protein sample was loaded onto each IPG strip. The samples were separated by IEF in the first dimension, and the second dimension SDS-PAGE was performed on 26 cm × 20 cm, 12% (v/v) polyacrylamide gels using an Ettan Dalt six multiple apparatus (GE Healthcare) at 16°C. After electrophoresis, the gels were stained with Coomassie Brilliant Blue G250, scanned with Imagescan, and analyzed with the ImageMaster software 5.0. Protein spots with more than 1.5- or 0.667-fold change from the normalized volume were considered differentially expressed. Differentially expressed protein spots were extracted from the gels and sequenced using tandem mass spectrometry (MS/MS) with a fuzzy logic feedback control and a Reflex III MALDI-TOF system (Bruker) equipped with delayed ion extraction.
Results
Plant Physiological Data and Nitrogen Utilization
By varying the ratio of N input with the rice growth stages, a significant effect on chlorophyll content was observed (Table 1, P ≤ 0.05). High chlorophyll content was observed at early full panicle stage (7–14 days) under both NF and NT, however NF significantly enhanced chlorophyll content at late booting stage (30 days) and late full panicle stage (7–21 days) compared to NT. A significantly higher N accumulation (281.11 ± 5.43) for NF than the NT (265.52 ± 325) was observed (Table 2, P ≤ 0.05). Similarly, nitrogen physiological efficiency and grain production efficiency (47.30 ± 207 and 56.75 ± 101) respectively were significantly higher under NF than NT (38.56 ± 1.63 and 53.57 ± 146) treatment. However, no significant difference was seen in nitrogen dry matter production efficiency with NT and NF treatments.
Table 1. The SPAD values of the leaves of rice plants subjected to different nitrogen treatments.
Table 2. The nitrogen utilization efficiencies of rice under different nitrogen treatments.
Grain and Yield Parameters of Rice
The measurement of yield and yield components is an important aspect to investigate the better cropping model based on two nitrogen rates. In 2012, NF treatment significantly increased number of grains per panicle by 5.02%, seed setting rate by 6.38%, 1000 grain weight by 4.53% and grain yield by 17.07%. However, for the number of productive panicles, no significant difference has been shown between NF and NT treatments. In 2013, the seed setting rate and grain yield (10.08 and 15.25%) respectively were significantly higher under NF than the NT (Table 3, P ≤ 0.05).
Table 3. The grain yield and its components of rice in the treatments with different nitrogen rates.
Soil Enzymatic Activity and Nutritional Status
Soil enzymes activity showed significant variability with N amendment. With NF significantly higher activity of urease and invertase (8.4 ± 0.32 and 34.11 ± 3.34) respectively was observed compared to urease (7.15 ± 0.21) and invertase (25.01 ± 0.16) under NT treatment (Table 4, P ≤ 0.05). However, under NT, nitrate reductase activity (2.35 ± 0.16) was higher than the NF, i.e., (1.44 ± 0.08). Similarly, high N and potassium (K) availability (98.28 ± 174 and 76.54 ± 301) respectively has been exhibited by NF compared to N (88.01 ± 1.59) and K (66.08 ± 2.91) in NT. Activities of soil enzymes (urease and invertase) was higher in the NF than that in the NT treatment, while no significant difference in Phosphorus (P) availability was observed for NT and NF treatments.
Table 4. The biochemical properties of the rhizosphere soils of rice under different nitrogen treatments.
Effect of Different Nitrogen Treatments on the Diversity and Evenness of Soil Microbes
The 16S rRNA fragments were digested by four restriction enzymes, i.e., MspI, HaeIII, AfaI, and AluI. These enzymes were combined for diversity and evenness analysis using the T-RFLP technique. The results showed that Shannon’s index, Pielou’s and Simpson’s indices were significantly higher in the NF than those in the NT treatment (Table 5).
Table 5. Diversity and evenness analysis of microbial community in the rhizospheric soil based on TRFLP data.
Effect of Different Nitrogen Rates on Soil Microbial Community
The presence and relative abundance of the main bacterial groups in the soil samples were estimated using NCBI databases from the T-RFLP profiles obtained in the MspI, HaeIII, AfaI, and AluI digests as mentioned above. A total of 135 bacteria were identified in the two treated soil samples, 67 and 68 types of bacteria were identified in the NF and NT treatments, respectively. These bacteria were grouped into Actinobacteria, Firmicutes, Chlorobi, Bacteroidetes, Tenericutes, Clostridium, Nitrospira, and some others (Figure 1). NF treatment increased the relative abundance of Proteobacteria by 2.59%, Actinobacteria by 0.19%, Tenericutes by 1.56% and some unknown bacteria by 9.03%, however, the relative abundance of Firmicutes, Chlorobi, and Bacteroidetes was decreased by 8.27, 0.24, and 8.42% respectively in NF treatment.
Figure 1. T-RFLP results showing Bacterial phyla percentage in the rice rhizosphere under different nitrogen rates. NT, traditional nitrogen application; NF, efficient nitrogen application. Vertical bars show standard deviations.
Dominant phyla with a relative abundance greater than 1% were further used for functional analysis. The results indicated that Proteobacteria, Actinobacteria, Firmicutes, and Bacteroidetes were the dominant phyla in the two soil samples. Some bacteria in the soil of the NF treatment have been studied to perform multiple ecological functions including bacteria participate in cellulose degradation, e.g., MspI 147 (Aneurinibacillus aneurinilyticus) (Raj et al., 2007), AluI 241 (Cytophaga spp.) (Li et al., 1997) and AluI 250 (Paenibacillus curdlanolyticus) (Schoina et al., 2011), bacteria involved in N cycle including AluI 228 (Moorella thermoacetica) (Silaghi-Dumitrescu et al., 2003), in nitrification, HaeII 331 (Nitrospira) (Dionisi et al., 2002), bacteria linked to pathogenicity, MspI 538 (Spiroplasma) (Tully et al., 1977) and HaeIII 209 (Staphylococcus kloosii) (Peer et al., 2011), and those of unknown function including AluI 282 (Haemophilus), AluI 229 (Aquifex pyrophilus), MspI 538 (Spiroplasma) (Supplementary Table S1). The bacteria identified in NT soil included the S cycle bacteria, MspI 74 (Desulfohalobium retbaense) (Ollivier et al., 1991), denitrification bacteria, MspI 88 (Neisseria denitrificans) (Stupak et al., 2015), mosquito eradication bacteria, HaeIII 308 (Bacillus sphaericus) (Myers and Yousten, 1978), the C cycle bacteria including HaeIII 338 (Weissella) (Choi et al., 2002), AluI 232 (Lactobacillus) (Cheng et al., 1991) and bacteria associated with pathogenicity are HaeIII 310 (Staphylococcus) (Lowy, 1998), MspI 148 (Aerococcus urinae) (Aguirre and Collins, 1992), AfaI 100 (Mycoplasma corogypsi) (Wettere et al., 2013), AfaI 106 (Mycoplasma neurolyticum) (Aleu and Thomas, 1966), and AluI 243 (Mobiluncus curtisii) (Meltzer et al., 2008) (Supplementary Table S2).
Effect of NF on Gene Copies of Enzymes Involved in Denitrification
The T-RFLP results indicated clear bacterial responses to the different nitrogen regimes with respect to nitrogen cycling. We performed qPCR analysis to confirm and quantify the genes involved in denitrification (narG encoding membrane-bound nitrate reductase and nirK encoding cd1 and copper nitrite reductase) in the soil samples (Figure 2). As expected, the qPCR results showed the amounts of the narG and nirK gene copies were increased by 272.00 and 99.92% respectively in the NT treatment, compared with those in the NF treatment.
Figure 2. Relative changes in the genes encoding the key enzymes of nitrate reduction (narG encoding membrane-bound nitrate reductase) and denitrification (nirK encoding cd1 and copper nitrite reductase) based on real-time PCR quantification involved in N cycling. Bars with different letters indicate significant differences at P ≤ 0.05. Different letters show significant differences determined by LSD’s test (P < 0.05, n = 4).
Metaproteomics Analysis of the Rice Rhizosphere Soil in Response to Different Nitrogen Treatments
Approximately 1021 and 970 protein spots were detected on the gel of the proteins extracted from the NF and NT soil samples respectively (Figure 3). At the same time, highly reproducible 2-DE maps were obtained from the two different soil samples with significant correlations among scatter plots. The correlation index between the NF and NT soil was 0.674 (Supplementary Figure S2). To obtain a metaproteomic profile for the rice soil, 167 protein spots with high resolution and repeatability, including all 67 differentially expressed proteins and 100 constitutively expressed proteins (Figure 3), were selected for identification and 128 protein spots were successfully analyzed using MALDI TOF-TOF MS (Supplementary Table S3).
Figure 3. (A) 2-D gel of proteins extracted from the NF soil. Arrows in figure point at proteins with differential expression compared to the NT. (B) 2-D gel of proteins extracted from the NT soil. Black circles in figure represent the same expression level.
According to Gene Ontology (GO) annotations, the identified proteins were classified into 14 Cellular Component (CC), 11 Molecular Function (MF), and 20 Biological Process (BP) categories (Supplementary Figure S3). Highly represented categories of the GO annotated proteins in MF were associated with cell (47.6%), cell part (47.6%), and catalytic activity (56.6%) and binding (48.8%), while metabolic process (71.4%) and cellular process (63.3%) in BP.
According to the putative physiological functions assigned using the KEGG database, these soil proteins were categorized into 11 groups, of which 88.28% were derived from microbes, 5.74% from plants, and 6.25% from fauna (Figure 4A). Most of these identified proteins were associated with protein metabolism constituting (19.53%), energy metabolism (11.72%), cell development and motility (13.28%), carbohydrate metabolism (9.38%), transcription (7.81%), signal transduction (7.81%), and defense response (11.72%). Based on the metaproteomic data, a tentative metabolic model for the rhizosphere soil proteins was proposed as shown in Figure 5. These soil proteins are collectively involved in carbohydrate/energy, protein and amino acid metabolism, transcription, signal transduction, defense response, genetic information processing, etc. Most of the microbe proteins were specifically related to carbohydrate/energy metabolism, and microbes might use the necessary energy and organic materials for their growth and material exchange. At the same time, some microbial proteins related to cell development and motility (including the surface layer protein, outer membrane protein A, flagellin, outer membrane porin OmpC family, ABC transporter, periplasmic oligopeptide-binding protein, OppA, etc.), signal transduction (including TonB-dependent receptor, UspA domain-containing protein, signal peptide protein, etc.), defense responses (including anti-oxidant AhpC-TSA family protein, heat shock protein Hsp70, Hsp20, superoxide dismutase, alkyl hydroperoxide reductase, etc.) and transcription (including DNA-directed RNA polymerase subunit alpha, transcriptional regulator, phage shock protein A, uroporphyrinogen decarboxylase, ribosomal RNA large subunit methyltransferase, spermidine synthase, etc.) were identified in the soil from the rice rhizosphere. However, some plant and fauna proteins were also identified in rhizosphere soil, and these may play an important role in the root colonization by microbes.
Figure 4. (A) Functional classification of all of the identified proteins. (B) The functional category distribution of all of the differentially expressed proteins originating from the microbes.
Figure 5. Proposed metabolic model for rhizosphere soil proteins as inferred by metaproteomic data. Identification numbers [E.C.-.-.-.-.] refer to the identified proteins. Red Upward arrows indicate the up-regulated proteins and green downward arrows show the down-regulated proteins. EMP, Embden-Meyerhof pathway; TCA, tricarboxylic acid cycle; GAC, glyoxylic acid cycle.
Differentially Expressed Proteins and Their Roles in the Rhizosphere Soils
The differentially expressed proteins have been shown to be originated from microbes (constituting 92.16%), plants (constituting 3.92%), and fauna (constituting 3.92%) (Table 6). Among these differentially expressed proteins of microbial origin, the largest functional group comprising 26.09% were the proteins involved in protein metabolism, 19.57% were related to cell development and motility, followed by 17.40 and 17.39% were of proteins associated with carbohydrate/energy metabolism and transcription respectively. The remainder of the proteins that originated from the microbes were related to defense response (13.04%), signal transduction (4.35%), and amino acid metabolism (2.17%) (Figure 4B). Furthermore, most of the microbial proteins related to protein metabolism (including spot 4, peptidase, M42 family protein; spot 8, elongation factor Tu; spot 19, co-chaperonin GroES; spot 42, elongation factor Tu and spot 67, chaperonin GroEL), cell development and motility (including spot 13, 20, 27 surface layer protein; spot 24, outer membrane protein A and spot 29, flagellin), energy metabolism (spot 28, nitroreductase family protein; spot 34, mitochondrial F-ATPase beta sub unit and spot 54, F0F1 ATP synthase subunit beta) and amino acid metabolism (spot17, dihydroorotase) were up-regulated in the NF soil compared with the NT soil (Table 6). These up-regulated microbe proteins involved in protein, energy, motility, and amino acid metabolism provide the necessary energy and organic materials for microbial growth, microbial mobility, and material exchange. It is likely that the microbes move to the rhizosphere to obtain carbohydrates originating from rice rhizosphere secretions. Most of the proteins related to plant stress defense (including spot 38, anti-oxidant AhpCTSA family protein; spot 37 and 50, heat shock protein Hsp20 and spot 64, heat shock protein 70) were up-regulated in the NF soil compared with the NT soil (Table 6). Several plant proteins related to amino acid metabolism (including spot 14, ethylene-responsive methionine synthase) and carbohydrate metabolism (including spot 61, glucanase) were up-regulated in the NT soil. In addition, one genetic information processing protein (spot 18, proteasome subunit beta type 3) and one protein with an unknown function (spot 51, female sterile) that originated from fauna were identified.
Table 6. Differentially expressed proteins identified by MALDI TOF-TOF MS.
Discussion
To achieve the sustainable development goals (SDGs), it is vital to investigate the impact of inorganic fertilizer applications on the structure and dynamics of the plant rhizosphere microbial community. The synchronization of N supply and crop need is essential for maximum uptake and utilization. In rice cultivation, appropriate reduction of N supply on the basal dressing accompanied by an increase during the late growth stage could significantly promote plant growth, especially at filling stage (Iwasaki et al., 1992; Ramasamy et al., 1997). Rice yield is significantly correlated with soil chemical properties and fertility (Lv et al., 2011). Some studies have confirmed the importance of available N and available K for the rice yield and soil fertility (Mussgnug et al., 2006; Boling et al., 2010). The present results clearly showed that, extending N application to the late growth stage of rice (NF treatment), highly synchronized N supply with the rice nutritional need, this in turn led to significant (P ≤ 0.05) increase in the chlorophyll content at late growth stage of rice. Moreover, utilization efficiency, seed-setting percentage and grain yields enhanced significantly compared to the traditional method of nitrogen application (NT). Other studies (Iwasaki et al., 1992; Ramasamy et al., 1997) also explained that an appropriate reduction of N supply on the basal dressing accompanied by an increase during the late growth stage could significantly promote rice growth, especially at filling stage.
Soil enzymes activities can reflect the status of soil microbial community and physico-chemical conditions of the soil (Puglisi et al., 2006; Burke et al., 2011). Soil nitrogen mineralization is a vital process, which supplies adequate amount of N needed for the plant growth and development (Chang et al., 2007). For N uptake by plants, an essential step in the process of N mineralization is the hydrolysis of urea by the activity of urease in the soil, which releases NH4 + -N (Sardans et al., 2008; Tao et al., 2009). Invertase is also crucial in releasing low molecular weight and simple forms of carbohydrates, which act as energy resources for microbes (Eliasson et al., 2005). Invertase catalyzes the hydrolysis of sucrose into glucose and fructose, and its activity determines the structure and abundance of the soil microbial community (Tao et al., 2009). The activities of the soil enzymes, such as urease and invertase, were significantly higher (P ≤ 0.05) in the NF than in the NT soil (Table 5), however, an opposite trend was true in the case of nitrate reductase. Low expression of narG and nirK genes, which encode the key enzymes of nitrate reduction under NF treatment, is consistent with the nitrate reductase activity. Knowles (1982) has indicated that nitrate reductase is a key enzyme in the process of denitrification under anaerobic conditions, which led to N2O emission to the atmosphere. A lower activity of nitrate reductase in NF soil can be linked to the shift in relative abundance of denitrifying bacteria, such as Neisseria denitrificans was lower in NF soil, while nitrite-oxidizing bacteria (NOB) were enhanced in NF soil, which catalyzes the conversion of NO2 to NO3, a major step in nitrification (Vanparys et al., 2007; Erguder et al., 2009), and hence restricting the N losses in the form of N2O or N2 gas (Knowles, 1982). Previously, it is reported that much of the added N may be lost to the system by denitrification (Thomas et al., 2006).
Soil microorganism diversity is crucial to determine soil organic matter decomposition, nutrient cycling, soil degradation, and bioremediation of soil contamination (Li et al., 2012). Shifts in the structure and composition of the microbial community are strong indicators of soil biological activity, soil quality and crop productivity of terrestrial agro-ecosystems (Edmeades, 2003). Based on results of diversity, the NF led to a significant increase (P ≤ 0.05) in the Shannon, Pielou, and Simpson diversity indices compared with the NT soil (Table 5), which led to more stable soil environment for plants (Chen et al., 2018). As we know, soil contains the largest pool of carbon on earth (Batjes, 1996). It was also found that the NF treatment increased the relative abundance of soil bacteria involved in cellulose degradation and the C cycle (Supplementary Table S1). Cellulose is the most abundantly produced biopolymer within this large carbon pool in the terrestrial environment and each year, photosynthetic fixation of CO2 yields more than 1011 tons of dry plant material worldwide (Li et al., 2009). As a large component of the plant structural carbon (30 to 50% of plant dry weight) (Schmidt, 2006), cellulose is one of the major constituents of soil carbon. The degradation of plant cellulose in soil is an important part of the terrestrial carbon cycle.
It is therefore suggested that the introduction of N cycle and C cycle bacteria have been shown to change NUE, the amount and composition of available nitrogen, and other essential nutrients in the rhizosphere soil. Shifts in the patterns of nitrogen and carbon utilization and the fixation potential of the microbial community in response to NF can have long-term effects on rice crop productivity.
In order to further unravel the rhizospheric biological process and its underlying mechanisms mediated by different nitrogen treatments (NF and NT) at late growing stage of rice, comparative metaproteomics was utilized to detect intricate interactions between nitrogen treatments, soil microflora and the plant root response in the form of root exudates in the rhizosphere ecosystem. Based on the metaproteomic data, most of the proteins (88.28%) were derived from microbes, with only 5.74% from plants, and 6.25% from fauna (Supplementary Table S3). In this study, the nitroreductase family protein [EC: 1.13.11.79] (spot 28) had a higher expression level in the NF compared with the NT soil. This family of proteins is commonly described as oxidoreductases (De Oliveira et al., 2007) and comprises a group of flavin adenine dinucleotide (FAD) or flavin mononucleotide (FMN) dependent enzymes that are able to metabolize nitroheterocyclic and nitroaromatic derivatives through the reducing power of nicotinamide adenine dinucleotide (NAD(P)H). These enzymes are found in bacterial species and to a lesser extent in eukaryotes (De Oliveira et al., 2010a). There is a little information on the biochemical functions of nitroreductases. Some studies suggest their possible involvement in the oxidative stress response (De Oliveira et al., 2010b).
The result of T-RFLP has shown that C cycles were accelerated in the NF treatment. Moreover, the analysis of differential soil metaproteomics showed that the mitochondrial F-ATPase beta subunit (spot 34), which had a higher expression level in the NF compared with the NT soil, undergoes a sequence of conformational changes leading to the formation of ATP from ADP (Leyva et al., 2003). The F0F1 ATP synthase subunit beta [EC: 3.6.3.14] (spot 54) also had a higher expression level in the NF soil; bacterial enzyme membrane-bound ATP synthases (F0F1-ATPases) catalyze the synthesis of ATP from ADP and inorganic phosphate through the power of an electrochemical ion gradient. On the other hand, ATP synthases function as ATPases under conditions of low driving force and ultimately generate a transmembrane ion gradient due to ATP hydrolysis (Deckers-Hebestreit and Altendorf, 1996). The terminal enzyme of the oxidative phosphorylation pathway, F0F1ATP synthase, is responsible for the majority of ATP synthesis in all living cells (Boyer, 1997). It was also found that TonB-dependent receptor (spot 62), a signal protein, had a higher expression level in the NF compared with the NT soil. From the outside of bacterial cells, TonB complexes sense signals and transmit them into the cytoplasm, which leads to transcriptional activation of target genes. The TonB protein interacts with outer membrane receptor proteins in Escherichia coli and carries out energy-dependent uptake of specific substrates and high-affinity binding in the periplasmic space (Chimento et al., 2003). TonB-dependent transporters (TBDTs) bind and transport ferric chelates called siderophores, as well as vitamin B12, nickel complexes, and carbohydrates (Noinaj et al., 2010).
Differential metaproteomics analysis also indicated that superoxide dismutase (SOD) [EC: 1.15.1.1] (spot 44) and alkyl hydroperoxide reductase/Thiol specific antioxidant/Mal allergen [EC: 1.11.1.15] (spot 56) related to defense responses were up-regulated in the NF compared with the NT soil. Superoxide dismutases (SOD) [EC: 1.15.1.1] catalyze the dismutation (or partitioning) of the superoxide () radical into either hydrogen peroxide (H2O2) or ordinary molecular oxygen (O2). Superoxide production is a by-product of oxygen metabolism and if not regulated properly, creates many kinds of cell damage. Hydrogen peroxide is also damaging, but less and is degraded by other enzymes such as catalase. Thus, SOD is an important antioxidant defense in almost all living cells exposed to oxygen. A ubiquitous family of antioxidant enzymes, i.e., peroxiredoxins (Prxs) [EC 1.11.1.15], also control cytokine-induced peroxide levels and therefore mediate signal transduction in mammalian cells (Rhee et al., 2005). Thiol-specific antioxidant (TSA), also referred to as TPx, is part of a novel family of antioxidant enzymes and was isolated initially from yeast and was later found in mammalian tissues. In the presence of a thiol reductant such as dithiothreitol (DTT) or thioredoxin, it has the capability to protect biomolecules from oxidative damage (Kim et al., 1988). One of the few proteins such as glucanase derived from plants was up-regulated. Glucanases are enzymes that break down glucans, i.e., polysaccharides comprised of several glucose sub-units, and play an important role in carbohydrate cycling. The up-regulated expression of the network of proteins involved in nitrification, cellulose degradation, energy, defense, and signaling mechanisms indicate the efficiency of the soil community to maintain the agricultural soil ecosystem to ensure healthy plant growth.
Although the number of the identified soil protein was limited, combined with soil enzyme assays and T-RFLP analysis, our metaproteomic results afford us a solid foundation to understand the interactions between the soil microorganisms and plants in the soil ecosystem.
Conclusion
Our study provides evidence that the application of nitrogen in four sessions (NF) was able to improve the physiological status of rice, especially at late growing stage, indicating higher leaf chlorophyll contents, increased NUE along with rice productivity in the NF treatment. Which resulted from the interaction of rice plants with the nitrogen nutrient and microbial flora in the rhizosphere soil at late growing stage of rice. It was found that the NF could elicit beneficial shifts in the composition and diversity of soil microbial communities at late growing stage of rice, enhance N availability, enzymatic activity, the abundance of bacteria functioning in cellulose degradation, C and N cycling as well as nitrification in contrast to the NT. In which there existed the higher abundance of undesired denitrification and S cycle bacteria detected by T-RFLP and qPCR approaches. The metaproteomic study indicated the abundance of microbe-derived proteins related to nitrogen fixation, energy, signaling, defense, carbohydrate metabolism, and the candidate proteins involved in these functions were up-regulated in NF soil, indicating a healthy network for nutrient cycling through the enhancement of the beneficial microbial community. Our results could open new avenues for modulating the root microbiome to enhance crop production and sustainability.
Data Availability Statement
The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.
Author Contributions
WL, JC, YA, and SL conceived the study. WL, JC, and IU wrote the manuscript. JC, BY, and JW performed the experiments. JC, LW, LZ, PL, and HW performed the statistical analyses. XQ and ZZ involved in the field management and soil sampling. All authors discussed the results and commented on the manuscript.
Funding
This work was supported by grants from the National Natural Science Foundation of China (Nos. 81573530, 31271670, and 31401306), Natural Science Foundation of Fujian Province (2017J01803), Project of Key Laboratory of Ministry of Education (GBMUC-2018-006), Science and Technology Development Fund of Fujian Agriculture and Forestry University (KF2015043), and Scientific Research Foundation of Graduate School of Fujian Agriculture and Forestry University. We are thankful to Major Agricultural Extension Services, Fujian Province, China (KNJ-153015), Weng Hongwu Academic Innovation Research Fund of Peking University and Original Research Fund of Peking University.
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2019.02623/full#supplementary-material
Footnotes
References
Abdelmagid, H., and Tabatabai, M. (1987). Nitrate reductase activity of soils. Soil Biol. Biochem. 19, 421–427. doi: 10.1016/0038-0717(87)90033-2
Aguirre, M., and Collins, M. D. (1992). Phylogenetic analysis of some Aerococcus-like organisms from urinary tract infections: description of Aerococcus urinae sp. nov. Microbiology 138, 401–405. doi: 10.1099/00221287-138-2-401
Aleu, F., and Thomas, L. (1966). Studies of pplo infection: iii. Electron microscopic study of brain lesions caused by Mycoplasma neurolyticum toxin. J. Exp. Med. 124, 1083–1088. doi: 10.1084/jem.124.6.1083
Arafat, Y., Wei, X. Y., Jiang, Y. H., Chen, T., Saqib, H., Lin, S., et al. (2017). Spatial distribution patterns of root-associated bacterial communities mediated by root exudates in different aged ratooning tea monoculture systems. Int. J. Mol. Sci. 18, E1727. doi: 10.3390/ijms18081727
Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., and Vivanco, J. M. (2006). The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol. 57, 233–266. doi: 10.1146/annurev.arplant.57.032905.105159
Bastida, F., García, C., von Bergen, M., Moreno, J. L., Richnow, H. H., and Jehmlich, N. (2015a). Deforestation fosters bacterial diversity and the cyanobacterial community responsible for carbon fixation processes under semiarid climate: a metaproteomics study. Appl. Soil Ecol. 93, 65–67. doi: 10.1016/j.apsoil.2015.04.006
Bastida, F., Selevsek, N., Torres, I. F., Hernández, T., and García, C. (2015b). Soil restoration with organic amendments: linking cellular functionality and ecosystem processes. Sci. Rep. 5:15550. doi: 10.1038/srep15550
Bastida, F., Jehmlich, N., Lima, K., Morris, B., Richnow, H., Hernández, T., et al. (2016a). The ecological and physiological responses of the microbial community from a semiarid soil to hydrocarbon contamination and its bioremediation using compost amendment. J. Proteomics 135, 162–169. doi: 10.1016/j.jprot.2015.07.023
Bastida, F., Torres, I. F., Moreno, J. L., Baldrian, P., Ondoño, S., Ruiz-Navarro, A., et al. (2016b). The active microbial diversity drives ecosystem multifunctionality and is physiologically related to carbon availability in Mediterranean semiarid soils. Mol. Ecol. 25, 4660–4673. doi: 10.1111/mec.13783
Bastida, F., Moreno, J., Nicolas, C., Hernandez, T., and Garcia, C. (2009). Soil metaproteomics: a review of an emerging environmental science. Significance, methodology and perspectives. Eur. J. Soil Sci. 60, 845–859. doi: 10.1111/j.1365-2389.2009.01184.x
Bastida, F., Torres, I. F., Andrés-Abellán, M., Baldrian, P., López-Mondéjar, R., Větrovský, T., et al. (2017). Differential sensitivity of total and active soil microbial communities to drought and forest management. Global Change Biol. 23, 4185–4203. doi: 10.1111/gcb.13790
Batjes, N. H. (1996). Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sc. 47, 151–163. doi: 10.1111/ejss.12115
Bjellqvist, B., Ek, K., Righetti, P. G., Gianazza, E., Görg, A., Westermeier, R., et al. (1982). Isoelectric focusing in immobilized pH gradients: principle, methodology and some applications. J. Biochem. Bioph. Meth. 6, 317–339. doi: 10.1016/0165-022X(82)90013-6
Boling, A., Tuong, T., Van Keulen, H., Bouman, B., Suganda, H., and Spiertz, J. (2010). Yield gap of rainfed rice in farmers’ fields in Central Java. Indonesia. Agr. Syst. 103, 307–315. doi: 10.1016/j.agsy.2010.02.003
Bouskill, N., Tang, J., Riley, W. J., and Brodie, E. L. (2012). Trait-based representation of biological nitrification: model development, testing, and predicted community composition. Front. Microbiol. 3:364. doi: 10.3389/fmicb.2012.00364
Boyer, P. D. (1997). The ATP synthase—a splendid molecular machine. Annu. Rev. Biochem. 66, 717–749. doi: 10.1146/annurev.biochem.66.1.717
Bremner, J. M., and Mulvaney, C. (1982). Nitrogen—Total 1. Methods of soil analysis. Part 2. Chem. Microbiol. Properties 72, 595–624. doi: 10.1016/B978-0-12-386454-3.00522-4
Burke, D. J., Weintraub, M. N., Hewins, C. R., and Kalisz, S. (2011). Relationship between soil enzyme activities, nutrient cycling and soil fungal communities in a northern hardwood forest. Soil Biol. Biochem. 43, 795–803. doi: 10.1016/j.soilbio.2010.12.014
Cardon, Z. G., and Whitbeck, J. L. (2007). The Rhizosphere: An Ecological Perspective. Amsterdam: Eleisver.
Chang, E. H., Chung, R. S., and Tsai, Y. H. (2007). Effect of different application rates of organic fertilizer on soil enzyme activity and microbial population. Soil Sci. Plant Nutr. 53, 132–140. doi: 10.1111/j.1747-0765.2007.00122.x
Chen, J., Arafat, Y., Wu, L. K., Xiao, Z. G., Li, Q. S., Muhammad, A. K., et al. (2018). Shifts in soil microbial community, soil enzymes and crop yield under peanut/maize intercropping with reduced nitrogen levels. Appl. Soil Ecol. 124, 327–334. doi: 10.1016/j.apsoil.2017.11.010
Chen, J., Wu, L. K., Xiao, Z. G., Wu, Y. H., Wu, H. M., Qin, X. J., et al. (2017). Assessment of the diversity of Pseudomonas spp. and Fusarium spp. in radix pseudostellariae Rhizosphere under monoculture by combining DGGE and quantitative PCR. Front. Microbiol. 8:1748. doi: 10.3389/fmicb.2017.01748
Cheng, P., Mueller, R. E., Jaeger, S., Bajpai, R., and Iannotti, E. (1991). Lactic acid production from enzyme-thinned corn starch using Lactobacillus amylovorus. J. Ind. Microbiol. Biot. 7, 27–34. doi: 10.1007/bf01575599
Chimento, D. P., Kadner, R. J., and Wiener, M. C. (2003). The Escherichia coli outer membrane cobalamin transporter BtuB: structural analysis of calcium and substrate binding, and identification of orthologous transporters by sequence/structure conservation. J. Mol. Biol. 332, 999–1014. doi: 10.1016/j.jmb.2003.07.005
Choi, H.-J., Cheigh, C.-I., Kim, S.-B., Lee, J.-C., Lee, D.-W., Choi, S.-W., et al. (2002). Weissella kimchii sp. nov., a novel lactic acid bacterium from kimchi. Int. J. Syst. Evol. Micr. 52, 507–511. doi: 10.1016/S0149-2063(02)00191-5
Chon, K., Chang, J.-S., Lee, E., Lee, J., Ryu, J., and Cho, J. (2011). Abundance of denitrifying genes coding for nitrate (narG), nitrite (nirS), and nitrous oxide (nosZ) reductases in estuarine versus wastewater effluent-fed constructed wetlands. Ecol. Eng. 37, 64–69. doi: 10.1016/j.ecoleng.2009.04.005
Chourey, K., Nissen, S., Vishnivetskaya, T., Shah, M., Pfiffner, S., Hettich, R. L., et al. (2013). Environmental proteomics reveals early microbial community responses to biostimulation at a uranium-and nitrate-contaminated site. Proteomics 13, 2921–2930. doi: 10.1002/pmic.201300155
De Oliveira, I. M., Bonatto, D., and Henriques, J. A. P. (2010a). “Nitroreductases: enzymes with environmental, biotechnological and clinical importance,” in Current Research, Technology and Education in Applied Microbiology and Microbial Biotechnology, ed. A. Mendez-Vilas, (Copenhagen: Formatex), 1008–1019.
De Oliveira, I. M., Zanotto-Filho, A., Moreira, J. C. F., Bonatto, D., and Henriques, J. A. P. (2010b). The role of two putative nitroreductases, Frm2p and Hbn1p, in the oxidative stress response in Saccharomyces cerevisiae. Yeast 27, 89–102. doi: 10.1002/yea.1734
De Oliveira, I. M., Henriques, J. A. P., and Bonatto, D. (2007). In silico identification of a new group of specific bacterial and fungal nitroreductases-like proteins. Biochem. Bioph. Res. Co. 355, 919–925. doi: 10.1016/j.bbrc.2007.02.049
Deckers-Hebestreit, G., and Altendorf, K. (1996). The F0F1-type ATP synthases of bacteria: structure and function of the F0 complex. Annu. Rev. Microbiol. 50, 791–824. doi: 10.1146/annurev.micro.50.1.791
Dionisi, H. M., Layton, A. C., Harms, G., Gregory, I. R., Robinson, K. G., and Sayler, G. S. (2002). Quantification of Nitrosomonas oligotropha-like ammonia-oxidizing bacteria and Nitrospira spp. from full-scale wastewater treatment plants by competitive PCR. Appl. Environ. Microbiol. 68, 245–253. doi: 10.1128/AEM.68.1.245-253.2002
Edmeades, D. C. (2003). The long-term effects of manures and fertilisers on soil productivity and quality: a review. Nutr. Cycl. Agroecosys. 66, 165–180. doi: 10.1023/a:1023999816690
Eliasson, P. E., McMurtrie, R. E., Pepper, D. A., Stromgren, M., Linder, S., and Agren, G. I. (2005). The response of heterotrophic CO2 flux to soil warming. Glob. Chang. Biol. 11, 167–181. doi: 10.1111/j.1365-2486.2004.00878.x
Erguder, T. H., Boon, N., Wittebolle, L., Marzorati, M., and Verstraete, W. (2009). Environmental factors shaping the ecological niches of ammonia-oxidizing archaea. FEMS Microbiol. Rev. 33, 855–869. doi: 10.1111/j.1574-6976.2009.00179.x
FAO (2009). Global Agriculture Towards 2050. In: How to feed the World 2050: High-Level Expert Forum. Rome: Food and Agriculture Organization of the United Nations.
Griffin, T. J., Gygi, S. P., Ideker, T., Rist, B., Eng, J., Hood, L., et al. (2002). Complementary profiling of gene expression at the transcriptome and proteome levels in Saccharomyces cerevisiae. Mol. Cell. Proteomics 1, 323–333. doi: 10.1074/mcp.M200001-MCP200
Guazzaroni, M. E., Herbst, F. A., Lores, I., Tamames, J., Peláez, A. I., Lopez-Cortés, N., et al. (2013). Metaproteogenomic insights beyond bacterial response to naphthalene exposure and bio-stimulation. ISME J. 7, 122–136. doi: 10.1038/ismej.2012.82
Hallin, S., Jones, C. M., Schloter, M., and Philippot, L. (2009). Relationship between N-cycling communities and ecosystem functioning in a 50-year-old fertilization experiment. ISME J. 3, 597–605. doi: 10.1038/ismej.2008.128
Hultman, J., Waldrop, M. P., Mackelprang, R., David, M. M., McFarland, J., Blazewicz, S. J., et al. (2015). Multi-omics of permafrost, active layer and thermokarst bog soil microbiomes. Nature 521, 208. doi: 10.1038/nature14238
Iwasaki, Y., Mae, T., Makino, A., Ohira, K., and Ojima, K. (1992). Nitrogen accumulation in the inferior spikelet of rice ear during ripening. Soil Sci. Plant. Nutr. 38, 517–525. doi: 10.1080/00380768.1992.10415083
Jackson, M. L. (2005). Soil Chemical Analysis: Advanced Course. Madison, UW: Libraries Parallel Press, 432–433. doi: 10.1097/01.ss.0000223482.86664.14
Kato, T., Shinmura, D., and Taniguchi, A. (2007). Activities of enzymes for sucrose-starch conversion in developing endosperm of rice and their association with grain filling in extra-heavy panicle types. Plant Prod. Sci. 10, 442–450. doi: 10.1626/pps.10.442
Kim, K., Kim, I., Lee, K.-Y., Rhee, S., and Stadtman, E. (1988). The isolation and purification of a specific “protector” protein which inhibits enzyme inactivation by a thiol/Fe (III)/O2 mixed-function oxidation system. J. Biol. Chem. 263, 4704–4711. doi: 10.1117/1.482739
Kim, O. S., Cho, Y. J., Lee, K., Yoon, S. H., Kim, M., Na, H., et al. (2012). Introducing EzTaxon-e: a prokaryotic 16S rRNA gene sequence database with phylotypes that represent uncultured species. Int. J. Syst. Evol. Micr. 62, 716–721. doi: 10.1099/ijs.0.038075-0
Leyva, J. A., Bianchet, M. A., and Amzel, L. M. (2003). Understanding ATP synthesis: structure and mechanism of the F1-ATPase. Mol. Membr. Biol. 20, 27–33. doi: 10.1080/0968768031000066532
Li, X., Zhou, Z. X., Zhang, S., Jin, F., and Gao, P. (1997). Production and properties of a cellobiose-oxidizing enzyme from a newly isolatedcellulolytic bacterium Cytophaga sp. LX-7. World J. Microb. Biot. 13, 683–688. doi: 10.1023/a:1018527106547
Li, X. H., Yang, H. J., Roy, B., Wang, D., Yue, W. F., Jiang, L. J., et al. (2009). The most stirring technology in future: cellulase enzyme and biomass utilization. Afr. J. Biotechnol. 8, 2418–2422. doi: 10.1186/1475-2859-8-30
Li, Y., Chen, Y. L., Li, M., Lin, X. G., and Lin, R. J. (2012). Effects of arbuscular mycorrhizal fungi communities on soil quality and the growth of cucumber seedlings in a greenhouse soil of continuously planting cucumber. Pedosphere. 22, 79–87. doi: 10.1016/S1002-0160(11)60193-8
Lin, W. X., Wu, L. K., Lin, S., Zhang, A. J., Zhou, M. M., Lin, R., et al. (2013). Metaproteomic analysis of ratoon sugarcane rhizospheric soil. BMC Microbiol. 13:135. doi: 10.1186/1471-2180-13-135
Liu, D., Keiblinger, K. M., Schindlbacher, A., Wegner, U., Sun, H., Fuchs, S., et al. (2017). Microbial functionality as affected by experimental warming of a temperate mountain forest soil—a metaproteomics survey. Appl. Soil Ecol. 117, 196–202. doi: 10.1016/j.apsoil.2017.04.021
Lowy, F. D. (1998). Staphylococcus aureus infections. N. Engl. J. Med. 339:2026. doi: 10.1056/NEJM199808203390806
Lv, M., Li, Z., Che, Y., Han, F., and Liu, M. (2011). Soil organic C, nutrients, microbial biomass, and grain yield of rice (Oryza sativa L.) after 18 years of fertilizer application to an infertile paddy soil. Biol. Fert. Soils. 47, 777–783. doi: 10.1007/s00374-011-0584-y
Meltzer, M. C., Desmond, R. A., and Schwebke, J. R. (2008). Association of Mobiluncus curtisii with recurrence of bacterial vaginosis. Sex. Transm. Dis. 35, 611–613. doi: 10.1097/olq.0b013e318167b105
Morales, S. E., Cosart, T., and Holben, W. E. (2010). Bacterial gene abundances as indicators of greenhouse gas emission in soils. ISME J. 4, 799–808. doi: 10.1038/ismej.2010.8
Mussgnug, F., Becker, M., Son, T., Buresh, R., and Vlek, P. (2006). Yield gaps and nutrient balances in intensive, rice-based cropping systems on degraded soils in the Red River Delta of Vietnam. Field Crop Res. 98, 127–140. doi: 10.1016/j.fcr.2005.12.012
Muyzer, G., and Smalla, K. (1998). Application of denaturing gradient gel electrophoresis (DGGE) and temperature gradient gel electrophoresis (TGGE) in microbial ecology. Antonie van Leeuwenhoek 73, 127–141. doi: 10.1023/A:1000669317571
Myers, P., and Yousten, A. A. (1978). Toxic activity of Bacillus sphaericus SSII-1 for mosquito larvae. Infect. Immun. 19, 1047–1053. doi: 10.1007/BF01646065
Nannipieri, P., Ascher, J., Ceccherini, M., Landi, L., Pietramellara, G., and Renella, G. (2003). Microbial diversity and soil functions. Eur. J. Soil Sc. 54, 655–670. doi: 10.1046/j.1351-0754.2003.0556.x
Noinaj, N., Guillier, M., Barnard, T. J., and Buchanan, S. K. (2010). TonB-dependent transporters: regulation, structure, and function. Annu. Rev. Microbiol. 64, 43–60. doi: 10.1146/annurev.micro.112408.134247
Ollivier, B., Hatchikian, C., Prensier, G., Guezennec, J., and Garcia, J. L. (1991). Desulfohalobium retbaense gen. nov., sp. nov., a halophilic sulfate-reducing bacterium from sediments of a hypersaline lake in Senegal. Int. J. Syst. Evol. Micr. 41, 74–81. doi: 10.1099/00207713-41-1-74
Osborn, A. M., Moore, E. R., and Timmis, K. N. (2000). An evaluation of terminal-restriction fragment length polymorphism (T-RFLP) analysis for the study of microbial community structure and dynamics. Environ. Microbiol. 2, 39–50. doi: 10.1046/j.1462-2920.2000.00081.x
Peer, M., Nasir, R., Kakru, D., Fomda, B., Bashir, G., and Sheikh, I. (2011). Sepsis due to linezolid resistant Staphylococcus cohnii and Staphylococcus kloosii: first reports of linezolid resistance in coagulase negative staphylococci from India. Indian J. Med. Microbi. 29, 60–62. doi: 10.4103/0255-0857.76527
Priemé, A., Braker, G., and Tiedje, J. M. (2002). Diversity of nitrite reductase (nirK and nirS) gene fragments in forested upland and wetland soils. Appl. Environ. Microbiol. 68, 1893–1900. doi: 10.1128/AEM.68.4.1893-1900.2002
Puglisi, E., Del Re, A., Rao, M., and Gianfreda, L. (2006). Development and validation of numerical indexes integrating enzyme activities of soils. Soil Biol. Biochem. 38, 1673–1681. doi: 10.1016/j.soilbio.2005.11.021
Raj, A., Chandra, R., Reddy, M. M. K., Purohit, H. J., and Kapley, A. (2007). Biodegradation of kraft lignin by a newly isolated bacterial strain, Aneurinibacillus aneurinilyticus from the sludge of a pulp paper mill. World J. Microb. Biot. 23, 793–799. doi: 10.1007/s11274-006-9299-x
Ramasamy, S., Ten Berge, H., and Purushothaman, S. (1997). Yield formation in rice in response to drainage and nitrogen application. Field Crop. Res. 51, 65–82. doi: 10.1016/S0378-4290(96)01039-8
Reyna, L., Wunderlin, D. A., and Genti-Raimondi, S. (2010). Identification and quantification of a novel nitrate-reducing community in sediments of Suquía River basin along a nitrate gradient. Environ. Pollut. 158, 1608–1614. doi: 10.1016/j.envpol.2009.12.014
Rhee, S. G., Chae, H. Z., and Kim, K. (2005). Peroxiredoxins: a historical overview and speculative preview of novel mechanisms and emerging concepts in cell signaling. Free Radical. Bio. Med. 38, 1543–1552. doi: 10.1016/j.freeradbiomed.2005.02.026
Sardans, J., Penuelas, J., and Estiarte, M. (2008). Changes in soil enzymes related to C and N cycle and in soil C and N content under prolonged warming and drought in a Mediterranean shrubland. Appl. Soil Ecol. 39, 233–235. doi: 10.1016/j.apsoil.2007.12.011
Schneider, T., Keiblinger, K. M., Schmid, E., Sterflinger-Gleixner, K., Ellersdorfer, G., Roschitzki, B., et al. (2012). Who is who in litter decomposition? Metaproteomics reveals major microbial players and their biogeochemical functions. ISME J. 6, 1749–1762. doi: 10.1038/ismej.2012.11
Schoina, C., Stringlis, I. A., Pantelides, I. S., Tjamos, S. E., and Paplomatas, E. J. (2011). Evaluation of application methods and biocontrol efficacy of Paenibacillus alvei strain K-165, against the cotton black root rot pathogen Thielaviopsis basicola. Biol. Control 58, 68–73. doi: 10.1016/j.biocontrol.2011.04.002
Silaghi-Dumitrescu, R., Coulter, E. D., Das, A., Ljungdahl, L. G., Jameson, G. N., Huynh, B. H., et al. (2003). A flavodiiron protein and high molecular weight rubredoxin from Moorella thermoacetica with nitric oxide reductase activity. Biochemistry 42, 2806–2815. doi: 10.1021/bi027253k
Smith, C. J., Nedwell, D. B., Dong, L. F., and Osborn, A. M. (2007). Diversity and abundance of nitrate reductase genes (narG and napA), nitrite reductase genes (nirS and nrfA), and their transcripts in estuarine sediments. Appl. Environ. Microbiol. 73, 3612–3622. doi: 10.1002/gepi.10246
Starke, R., Bastida, F., Abadía, J., García, C., Nicolás, E., and Jehmlich, N. (2017). Ecological and functional adaptations to water management in a semiarid agroecosystem: a soil metaproteomics approach. Sci. Rep. 7:10221. doi: 10.1038/s41598-017-09973-w
Starke, R., Kermer, R., Ullmann-Zeunert, L., Baldwin, I. T., Seifert, J., Bastida, F., et al. (2016). Bacteria dominate the short-term assimilation of plant-derived N in soil. Soil Biol. Biochem. 96, 30–38. doi: 10.1016/j.soilbio.2016.01.009
Stupak, R., Makauskas, N., Radzevičius, K., and Valančius, Z. (2015). Optimization of intracellular product release from Neisseria denitrificans using microfluidizer. Prep. Biochem. Biotech. 45, 667–683. doi: 10.1080/10826068.2014.940539
Tao, W. H., Shen, Y. F., and Li, S. Q. (2009). Effect of N fertilization on the mineral nitrogen in soil in wind-water erosion crisscross region on the Loess Plateau. Journal of Northwest A & F University 10, 103–108.
Thomas, R. B., Van Bloem, S. J., and Schlesinger, W. H. (2006). “Climate change and symbiotic nitrogen fixation in agroecosystems,” in Agroecosystems in a Changing Climate, eds G. R. Edwards, P. C. D. Newton, and P. A. Niklaus, (Boca Raton, FL: CRC Press), 85–116. doi: 10.1201/9781420003826.ch4
Tully, J. G., Whitcomb, R. F., Clark, H. F., and Williamson, D. L. (1977). Pathogenic mycoplasmas: cultivation and vertebrate pathogenicity of a new spiroplasma. Science 195, 892–894. doi: 10.1126/science.841314
Vanparys, B., Spieck, E., Heylen, K., Wittebolle, L., Geets, J., Boon, N., et al. (2007). The phylogeny of the genus Nitrobacter based on comparative rep-PCR, 16S rRNA and nitrite oxidoreductase gene sequence analysis. Syst. Appl. Microbiol. 30, 297–308. doi: 10.1016/j.syapm.2006.11.006
Wang, H. B., Zhang, Z. X., Li, H., He, H. B., Fang, C. X., Zhang, A. J., et al. (2010). Characterization of metaproteomics in crop rhizospheric soil. J. Proteome Res. 10, 932–940. doi: 10.1021/pr100981r
Wang, Q. Y., Zhou, D. M., and Cang, L. (2009). Microbial and enzyme properties of apple orchard soil as affected by long-term application of copper fungicide. Soil Biol. Biochem. 41, 1504–1509. doi: 10.1016/j.soilbio.2009.04.010
Wettere, A. V., Ley, D., Scott, D., Buckanoff, H., and Degernes, L. (2013). Mycoplasma corogypsi–associated polyarthritis and tenosynovitis in black vultures (Coragyps atratus). Vet. Pathol. 50, 291–298. doi: 10.1177/0300985812457791
Williams, M. A., Taylor, E. B., and Mula, H. P. (2010). Metaproteomic characterization of a soil microbial community following carbon amendment. Soil Biol. Biochem. 42, 1148–1156. doi: 10.1016/j.soilbio.2010.03.021
Witt, C., Dobermann, A., Abdulrachman, S., Gines, H., Guanghuo, W., Nagarajan, R., et al. (1999). Internal nutrient efficiencies of irrigated lowland rice in tropical and subtropical Asia. Field Crop. Res. 63, 113–138. doi: 10.1016/S0378-4290(99)00031-3
Xing, Y. Z., and Zhang, Q. F. (2010). Genetic and molecular bases of rice yield. Annu. Rev. Plant Bio.l 61, 421–442. doi: 10.1146/annurev-arplant-042809-112209
Yan, T. F., Fields, M. W., Wu, L. Y., Zu, Y. G., Tiedje, J. M., and Zhou, J. Z. (2003). Molecular diversity and characterization of nitrite reductase gene fragments (nirK and nirS) from nitrate-and uranium-contaminated groundwater. Environ. Microbiol. 5, 13–24. doi: 10.1046/j.1462-2920.2003.00393.x
Yang, J., and Zhang, J. (2009). Grain-filling problem in ‘super’rice. J. Exp. Bot. 61, 1–5. doi: 10.1093/jxb/erp348
Zahar, H. F. E., Marol, C., Berge, O., Rangel-Castro, J. I., Prosser, J. I., Balesdent, J., et al. (2008). Plant host habitat and root exudates shape soil bacterial community structure. ISME J. 2, 1221–1230. doi: 10.1038/ismej.2008.80
Zampieri, E., Chiapello, M., Daghino, S., Bonfante, P., and Mello, A. (2016). Soil metaproteomics reveals an inter-kingdom stress response to the presence of black truffles. Sci. Rep. 6:25773. doi: 10.1038/srep25773
Zhang, A., Liu, R., Gao, J., Yang, S., and Chen, Z. (2014). Regulating N application for rice yield and sustainable eco-agro development in the upper reaches of Yellow River Basin, China. The Scientific World J. 2014, 239279. doi: 10.1155/2014/239279
Zhang, Z. X., Chen, J., Li, Z., Chen, H. F., Fang, C. X., Huang, J., et al. (2012). Differential proteomic expressions between superior and inferior spikelets of rice in response to varied nitrogen treatments. Aust. J. Crop. Sci. 6, 316–325.
Keywords: nitrogen utilization efficiency, T-RFLP, soil enzyme, soil metaproteomics, rice
Citation: Chen J, Arafat Y, Ud Din I, Yang B, Zhou L, Wang J, Letuma P, Wu H, Qin X, Wu L, Lin S, Zhang Z and Lin W (2019) Nitrogen Fertilizer Amendment Alter the Bacterial Community Structure in the Rhizosphere of Rice (Oryza sativa L.) and Improve Crop Yield. Front. Microbiol. 10:2623. doi: 10.3389/fmicb.2019.02623
Received: 18 March 2019; Accepted: 28 October 2019;
Published: 14 November 2019.
Edited by:
Helena Freitas, University of Coimbra, PortugalReviewed by:
Wu Xiong, Utrecht University, NetherlandsQiangsheng Wang, Nanjing Agricultural University, China
Copyright © 2019 Chen, Arafat, Ud Din, Yang, Zhou, Wang, Letuma, Wu, Qin, Wu, Lin, Zhang and Lin. 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: Zhixing Zhang, zinsser347@126.com; Wenxiong Lin, lwx@fafu.edu.cn