Next Article in Journal
Forage Yield, Quality, and Impact on Subsequent Cash Crop of Cover Crops in an Integrated Forage/Row Crop System
Next Article in Special Issue
Mitigation of CaCO3 Influence on Ipomoea batatas Plants Using Bacillus megaterium DSM 2894
Previous Article in Journal
Improvements in Sugarcane (Saccharum spp.) Varieties and Parent Traceability Analysis in Yunnan, China
Previous Article in Special Issue
Influence of Organomineral Fertiliser from Sewage Sludge on Soil Microbiome and Physiological Parameters of Maize (Zea mays L.)
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Agronomy
Volume 12
Issue 5
10.3390/agronomy12051213
agronomy-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Combined Organic and Inorganic Fertilization Can Enhance Dry Direct-Seeded Rice Yield by Improving Soil Fungal Community and Structure

by
Xiaohong Guo
1,
Jiajun Liu
1,
Lingqi Xu
1,
Fujing Sun
2,
Yuehan Ma
1,
Dawei Yin
1,
Qiang Gao
1,
Guiping Zheng
1 and
Yandong Lv
1,*
1
Heilongjiang Provincial Key Laboratory of Modern Agricultural Cultivation and Crop Germplasm, College of Agriculture, Heilongjiang Bayi Agricultural University, Daqing 163319, China
2
Qingdao Pony Test Co., Ltd., Qingdao 266104, China
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(5), 1213; https://doi.org/10.3390/agronomy12051213
Submission received: 14 April 2022 / Revised: 13 May 2022 / Accepted: 16 May 2022 / Published: 18 May 2022
(This article belongs to the Special Issue A Circular Economy: Chemical, Microbiological and Environmental Implications of Mineral and Organic Fertilizers Use in Soils)
Browse Figures
Review Reports Versions Notes

Abstract

:
Direct seeding of rice has emerged as a strategy for sustainable rice (Oryza sativa L.) production because of advantages, such as fewer production links, labor and farmland water-saving, easy mechanization, and high economic benefits. However, few studies have investigated the effects of different organic fertilizers on soil fungal community and rice yield in dry direct-seeded paddy fields. In order to select the best combination of organic fertilizer and chemical fertilizer, field experiments were used to evaluate the role of no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer; and OF4, attapulgite organic fertilizer + conventional NPK fertilizer on microbial structure and diversity and rice yield. Under Jishiwang bioorganic fertilizer + conventional NPK fertilization, the number of fungal OTUs was 365 and ranged from 1 to 9. The Ascomycota relative abundance was increased by 28.25% under Jishiwang bioorganic fertilizer application compared with CF, but the Basidiomycota decreased. Sordariomycetes and Leotiomycetes relative abundances were increased under organic fertilization. The relative abundance of dung saprotrophs, fungal parasites, and leaf saprotrophs was increased under organic fertilizer compared to CF, and animal pathogens decreased, but organic fertilizers also increased plant pathogens. Rice yield was increased under Jishiwang bioorganic fertilizer + conventional NPK fertilizer and was positively correlated with Ascomycota and Sordariomycetes relative abundances. The use of Jishiwang bioorganic fertilizer + conventional NPK fertilizer improves fungal community diversity and rice yield.

1. Introduction

Rice (Oryza sativa L.) is one of the most important food crops in the world, which can meet the food needs of more than 3 billion people in more than 100 countries across the world [1]. However, due to the increasing shortage of freshwater for crop irrigation and the reduction of the labor force in the rice-planting industry, traditional Chinese rice cultivation is under threat [2,3]. In this context, dry direct-seeded rice has emerged as a simplified and sustainable cultivation technology of great significance, but yields may be limited [4]. Due to phenological and physiological differences [5], the traditional fertilization technology applied to transplanted rice is not suitable for dry direct-seeded rice [6]. Till now, most studies have focused on the effect of chemical fertilizers on dry direct-seeded rice [7,8]. The excessive use of chemical fertilizers is the main cause of resource waste and environmental damage [9]. Organic fertilizers contain multiple nutrients required for plant growth, such as trace elements, sugars, and fats. They can improve the physicochemical properties of the soil, enhance microbial abundance and activities, and improve crop yield and quality [10,11]. Therefore, it is important to identify the best combination of organic and chemical fertilizers with respect to dry direct-seeded rice.
The soil microbial community includes two important components, soil fungi and bacteria. The soil fungal community plays diverse functions in the ecosystem and is involved in many ecological processes, such as improving plant nutrient utilization through biological nitrogen fixation and phosphorus dissolution in the soil and rhizosphere and enhancing plant growth and development [12,13,14], thereby playing an important role in various types of terrestrial ecosystems. Fungi are more tolerant to water stress compared with bacteria [15] because hyphae can explore the soil more profusely, gaining greater access to soil water [16]. Better access to soil pores through fungal hyphae allows a more efficient nutrient acquisition [17], thus alleviating the pressure of water stress. On the other hand, soil fungi also have important effects on soil organic compounds [18], thus affecting soil nutrient conversion [19]. Some studies also show that soil fungi are more sensitive to fertilization management measures than bacteria [20]. It has been reported that the long-term application of inorganic fertilizers reduces fungal diversity in soil under wheat cultivation [21]. Compared with chemical fertilizers, organic fertilizers improve soil fungal diversity and community structure [22]. Likewise, the combination of organic and inorganic fertilizers may increase the diversity of soil bacteria and fungi, making its distribution more homogeneous than that found in soils under inorganic fertilizers and in unfertilized soils [23]. It has also been reported that even in the short term, replacing part of the chemical fertilizer with organic fertilizers can optimize soil fungal community structure, promote soil beneficial fungi multiplication, and inhibit harmful soil bacteria [24]. Therefore, exploring the changes in soil fungal communities under different fertilization regimes is important to improve our understanding of the rice-microbial ecosystem.
At present, studies on the impact of organic fertilizers on soil fungal diversity are limited [25,26], and the consequences of such practices in dry direct-seeded rice fields were rarely reported. In this study, Illumina Miseq high-throughput sequencing was used to analyze the influence of different organic fertilizers on the soil fungal community and rice yield in a dry direct-seeded rice field, and to explore the best combination of organic and chemical fertilizers, in order to provide a theoretical basis for optimizing the soil microflora and selecting appropriate types of organic fertilizer in dry direct-seeded rice field.

2. Materials and Methods

2.1. Experimental Site

The fertilization regimes described in the following section were applied for three consecutive dry direct-seeded rice campaigns (2018–2020) under field conditions; the plots were located in Youyi farm, Heilongjiang Province (46.45° N, 131.49° E). The climate in this region corresponds to middle temperate continental monsoon, with annual mean temperature and precipitation of 3.1 °C and 501.2 mm, respectively (Figure 1). The soil type in this area is meadow black soil. The physicochemical properties of the soil at a depth of 0–20 cm are shown in Table 1.

2.2. Experimental Design

Dry direct-seeded rice was subjected to 6 experimental treatments under a complete randomized design, with the fertilization regime as the sole variable factor, including no fertilizer (F0), conventional NPK fertilizer (CF), and 4 combined organic-inorganic fertilizers: biochar (Shenyang Longtai Biological Engineering Co., Ltd., Shenyang, Liaoning Province, China) + conventional NPK fertilizer (OF1), seaweed bioorganic fertilizer (Qingdao Pulebaiwo Biotechnology Co., Ltd., Qingdao, Shandong Province, China) + conventional NPK fertilizer (OF2), Jishiwang bioorganic fertilizer (Jiamusi Sanxing Agricultural Technical Service Co., Ltd., Jiamusi, Heilongjiang Province, China) + conventional NPK fertilizer (OF3), and attapulgite organic fertilizer (Dingfengyuan concave High-tech Development Co., Ltd., Zhangye, Gansu Province, China) + conventional NPK fertilizer (OF4). Conventional NPK fertilizer (CF) was used as the control, and no fertilizer use (F0) was also included.
In the conventional fertilization regime, a total of 204.05 kg N ha−1 (urea, 46% N) was applied, distributed as basal fertilizer (88.82 kg N ha−1), tillering fertilizer (69.00 kg ha−1), and panicle fertilizer 46.23 kg N ha−1. The total amount of phosphate fertilizer (calcium triple superphosphate, 46% P2O5) was P2O5 80.04 kg ha−1, all as base fertilizer, and the total amount of potassium fertilizer (potassium sulfate, 50% K2O) was K2O 115.00 kg ha−1, divided into basal (63.75 kg K2O ha−1) and panicle (51.25 kg K2O ha−1) fertilizers. Organic fertilizers comprised biochar (organic matter 26.5%); seaweed bioorganic fertilizer (N + P2O5 + K2O ≥ 6.0%, organic matter ≥ 50.0%, living bacteria ≥ 20 million g−1); Jishiwang bioorganic fertilizer (N +P2O5 +K2O ≥ 6.0%, organic matter content ≥ 40.0%, living bacteria ≥ 20 million g−1), and attapulgite organic fertilizer (N + P2O5 + K2O ≥ 6.0%, organic matter ≥ 41.5%); these four compounds were applied before tillage at the doses of 15 t ha−1, 300 kg ha−1, 300 kg ha−1, and 750 kg ha−1, respectively.
The rice variety Longjing 31 (11 leaves on the main stem) was used in this experiment; this genotype is characterized by about 130 days to maturity and an active accumulated temperature ≥ 10 °C of about 2350 °C.
Before sowing, the seeds were coated with the Liangdun seed coating agent, and the dry seeds were mechanically sown on 30 April 2018, 18 April 2019, and 19 April 2020. In all cases, the seeding rate was 210 kg ha−1, and the seeding depth was 2 cm, with a row spacing of 20 cm. Each experimental plot covered 300 m2 (length 20 m, width 15 m), and 3 replicates per treatment were set, making a total of 18 plots. Basal fertilizer was applied at 4 cm aside from the seedling belt, 5 cm depth; a soil covering chain was used, and the soil was compacted twice after sowing. Harvest was performed on 20 September, 23 September, and 27 September in the rice production cycles of 2018, 2019, and 2020, respectively.

2.3. Soil Sampling and Rice Yield and Yield Components Measurements

Four points in each experimental plot were randomly selected to obtain four soil samples from the 0–10 and 10–20 cm soil layers after the rice harvest in September 2020; subsamples from each layer were collected and mixed (same layer). Fresh soil samples were transported in ice bags to the laboratory and kept at −80 °C until analysis. In addition, 3 points in each plot were randomly selected at rice maturity to collect 30 representative panicles to calculate the theoretical yield based on the measurements of yield components. The number of panicles per hill was determined and converted to the number of panicles per square meter. Filled and unfilled spikelets were separated by submerging them in tap water, counted, and weighed. Spikelet number per panicle, seed-setting rate, and 1000-grain weight were determined. The actual yield was assessed by collecting rice plants from a 5 m2-area in each plot; final data were obtained after adjusting the moisture content to 14%.

2.4. Soil DNA Extraction and Illumina HiSeq Sequencing

Total fungal DNA was extracted from the soil samples collected using Power Soil DNA Isolation Kit (MO BIO Laboratories), according to the manufacturer’s protocol. DNA quality and quantity were assessed by measuring the absorbances and calculating 260 nm/280 nm and 260 nm/230 nm ratios. The DNA was stored at −80 °C until further processing.
The fungal ITS1 region was amplified using primers ITS1-F (5′-CTTGGTCATTTAGAGGAAGTAA-3′) and ITS2-R (5′-GCTGCGTTCTTCATCGATGC-3′) combined with Illumina adapter sequences and barcodes. The PCR was performed in a total reaction volume of 20 µL: H2O 13.25 µL, 10× PCR ExTaq Buffer 2.0 µL, DNA template (100 ng/mL) 0.5 µL, primer 1 (10 mmol/L) 1.0 µL, primer 2 (10 mmol/L) 1.0 µL, dNTP 2.0 µL, and ExTaq polymerase (5 U/mL) 0.25 µL. After an initial denaturation at 95 °C for 5 min, amplification was achieved by 30 cycles of 30 s at 95 °C, 20 s at 58 °C, and 6 s at 72 °C, followed by a final extension step at 72 °C for 7 min. Then, the amplified products were purified and recovered using 1.0% agarose gel electrophoresis. The library construction and sequencing steps were performed by Beijing Biomarker Technologies Co., Ltd.

2.5. Bioinformatics and Statistical Analyses

The bioinformatic analysis in this study was completed using the Biomarker Biocloud Platform (www.biocloud.org, accessed on 5 May 2021). To obtain the raw tags, paired-end reads were merged using FLASH (v1.2.7, http://ccb.jhu.edu/software/FLASH/, accessed on 16 May 2021) [27]. Then, raw tags were filtered and clustered in the next steps. The merged tags were compared with the primers, and the tags with more than six mismatches were discarded by the FASTX-Toolkit. Tags with an average quality score < 20 in a 50 bp sliding window were truncated using Trimmomatic (http://www.usadellab.org/cms/?page=trimmomatic, accessed on 25 May 2021) [28], and tags shorter than 350 bp were removed. We identified possible chimeras by employing UCHIME, a tool included in mothur (http://drive5.com/uchime, accessed on 28 May 2021). The denoised sequences were clustered using USEARCH (version 10.0), and tags with a similarity ≥ 97% were regarded as an operational taxonomic unit. Taxonomic classification was assigned to all units by searching against the Silva database (http://www.arb-silva.de, accessed on 8 June 2021) using the UCLUST classifier in QIIME.
Four fungal alpha diversity indices (Chao1 richness estimator, ACE richness estimator, Shannon–Wiener diversity index, and Simpson diversity index) were calculated by mothur v1.30 [29]. The variation in fungal communities among soil samples was visualized by a principal coordinates analysis (PCoA) based on the Bray–Curtis dissimilarity index using the Vegan package.
The data from the 3 years were pooled, and group means were assessed using one-way ANOVA followed by the least significant difference (LSD) test, performed using IBM SPSS v20.0 statistical software (SPSS, Chicago, IL, USA). A Pearson correlation analysis was conducted to assess the relationships between yield and dominant soil fungi at phylum and class levels. Figures were generated using GraphPad Prism 7 and Origin 2021 software.

3. Results

3.1. Paddy Soils Sequencing Throughput

The raw sequences obtained through Illumina sequencing from the soil samples collected were qualitatively controlled, assembled, and filtered for chimeras. A total of 49,855–62,521 raw reads and 46,126–58,375 effective reads were obtained for fungi (Table 2); the effective reads accounted for 89.80–93.37% of the raw reads. The number of raw reads and effective reads and the percentage of effective reads were lower under all combined treatments (OF1, OF2, OF3, and OF4) than CF. The sequences were similar in length across all treatments, ranging between 241 and 251 bp, and the coverage level was over 0.99 in all cases, indicating that the sequencing data were reliable and could truly reflect the community composition of soil fungi.

3.2. Effect of Organic-Inorganic Fertilizers on the Alpha Diversity of the Soil Fungal Community

Figure 2a–d shows the effect of the fertilization treatments applied to the dry direct-seeded rice on several fungi alpha diversity indices. Unfertilized plots (F0) exhibited smaller fungi species richness (ACE and Chao1) and diversity (Simpson) than fertilized plots. Moreover, species richness under combined organic-inorganic fertilizers tended to diminish compared with CF, but no significant differences were detected. Among fertilized treatments, OF1 was associated with the lowest ACE and Chao1 indices and the highest Simpson and Shannon indices, followed by those found for OF3. OF2 and OF4 led to minor changes in these parameters compared with CF. Briefly, treatment with organic fertilizers had different effects. The ACE, Chao1, and Simpson and Shannon indexes changed between different treatments, but there was no significant difference with CF. This shows that organic fertilizer type and composition are relevant for determining the alpha diversity of soil fungi communities in dry direct-seeded rice fields.

3.3. Effect of Combined Organic-Inorganic Fertilizers on Soil Fungal OTUs

As the Venn graph depicts, the number of fungal OTUs shared by all soil samples was 365 (Figure 3), while a reduced number of unique OTUs were recognized for treatments F0, CF, OF1, OF2, OF3, and OF4: 3, 7, 9, 1, 2, and 6, respectively. This result indicates that in the dry direct-seeded rice field, the application of biochar tended to increase the number of fungal species in the soil compared with conventional fertilization. In contrast, the application of the seaweed bioorganic fertilizer and the Jishiwang bioorganic fertilizer reduced these values, a finding that may be related to the competitive pressure exerted by the microorganisms contained in both organic fertilizers.
The differences in soil fungal community composition among soils subjected to different fertilization regimes during dry direct-seeded rice cultivation were subsequently analyzed by a principal coordinates analysis (PCoA) at the OTU level (Figure 4). This analysis revealed that fungal community variability under each fertilization treatment was explained by 37.69% and 11.74% in the axes PC1 and PC2, respectively; the total interpretation rate was 49.44%. It may be noticed that OTUs associated with the fungal community in non-fertilized plots (F0) were mainly distributed on the zero axis of PC2, scattered and away from those corresponding to fertilized treatments. OTUs representing the fungal community in CF, OF1, and OF3 were concentrated in the second quadrant and clustered, suggesting that the fungal community composition in these samples was similar. OTUs linked to the fungal community in OF2 and OF4 occupied mainly the third quadrant, close to each other and CF, OF1, and OF3, indicating that the application of organic-inorganic fertilizers influenced the fungal community structure in dry direct-seeded rice fields.

3.4. Effect of Organic-Inorganic Fertilizers on Soil Fungal Community at Phyla and Class Levels

In this study, eight fungal phyla were identified in the soil samples analyzed, and five phyla had a relative abundance > 1%: Ascomycota, Basidiomycota, Mortierellomycota, Rozellomycota, and Chytridiomycota (Figure 5a). Compared with CF, the relative abundance of Ascomycota in OF1 and OF3-treated soils increased by 2.49% and 28.25%, respectively, while it decreased by 20.05% and 17.99% in OF2 and OF4. Interestingly, Basidiomycota’s relative abundance showed the opposite trend of Ascomycota, with lower values under OF1 and OF3 compared with CF and higher values under OF2 (39.06%) and OF4 (63.51%). It is worth mentioning that the introduction of organic fertilizers in dry direct-seeded paddy soil resulted in an inhibitory effect on Mortierellomycota, reducing its relative abundance by 15.24% to 39.98% compared with CF. Rozellomycota’s relative abundance under OF1, OF2, and OF4 was greater than that found in CF; the highest increase corresponded to OF2 (58.21% over CF), while the highest decrease found corresponded to OF3 (69.49% below CF). Compared with CF, the relative abundance of soil Chytridiomycota increased by 104.23% and 31.54% in rice plots subjected to OF1 and OF4, respectively, and decreased by 6.15% and 24.23% in those under OF2 and OF3. Rice plots treated with OF1 had a notably greater proportion of Ascomycota, Rozellomycota, and Chytridiomycota (compared to CF), whereas Basidiomycota was dominant in plots treated with OF2 and OF4. Mortierellomycota only occupied a dominant position in CF.
At the class level, eight dominant fungal classes with relative abundance > 1% were identified, comprising, in decreasing order, Sordariomycetes (25.00–32.77%), Mortierellomycetes (7.63–12.24%), Dothideomycetes (5.09–10.43%), Leotiomycetes (2.45–13.75%), Tremellomycetes (1.76–6.44%), Agaricomycetes (1.53–4.89%), Eurotiomycetes (1.24–5.21%), and Pezizomycetes (0.84–3.92%) (Figure 5b). The relative abundance of Sordariomycetes under OF3 and OF4 increased compared to CF (7.09% and 3.86%, respectively), while it decreased in OF1 and OF2 treatments (8.14% and 19.61%, respectively). Furthermore, the incorporation of organic fertilizers considerably decreased the relative abundance of Mortierellomycetes and Tremellomycetes, with the lowest values for the former in OF1 (37.66% below CF) and for the latter in OF2 (72.67% below CF).
All organic-inorganic fertilizers increased the relative abundance of Pezizomycetes; increases reached 3.67 times (OF1), 1.60 times (OF2), 0.65 times (OF3), and 1.40 times (OF4) compared to CF. The relative abundance of Dothideomycetes in the soils sampled from OF1, OF2, and OF3 showed an increasing trend compared with CF, with the greatest change in OF1. The highest relative abundance of Leotiomycetes was detected under OF3, followed by that found under OF4 (3.15 and 0.83 times over CF, respectively). The relative abundance of Agaricomycetes and Eurotiomycetes was lowest in OF3. Sordariomycetes, Pezizomycetes, Dothideomycetes, Leotiomycetes, Agaricomycetes, and Eurotiomycetes occupied a relatively dominant position in the plots exposed to the combined organic and inorganic fertilization treatments, especially OF3.

3.5. Ecological Function of Soil Fungi

We explored the putative ecological function of the fungal community deduced from the sequencing analysis based on the FUNGuild software (Figure 6). The soil fungal community across the 6 treatments tested could be ascribed to 10 ecological functional groups: undefined saprotroph, dung saprotroph, wood saprotroph, animal pathogen, plant pathogen, plant saprotroph, fungal parasite, soil saprotroph, endophyte, and leaf saprotroph. Notably, the relative abundance of animal pathogens was greater than that of dung saprotrophs under CF, while the relative abundance of wood saprotrophs was greater than that of dung saprotrophs under OF3 treatment, which suggests that OF3 treatment had a particular impact on functional group distribution in the fungal community.
The relative abundances of dung saprotroph, plant pathogen, fungal parasite, and leaf saprotroph groups under organic-inorganic fertilizers were all higher than those found under conventional fertilizer treatment, while animal pathogens showed the opposite trend. The relative abundance of wood saprotroph and endophyte groups increased in OF2, OF3, and OF4 compared with conventional fertilization. Soil saprotrophs were at the highest relative abundance in the soil under attapulgite fertilizer. It can be seen that the long-term use of organic fertilizer increased the relative abundance of wood saprotroph and endophyte groups but decreased the relative abundance of animal pathogens.

3.6. Yield and Yield Components under Different Fertilizer Treatments in Dry Direct-Seeded Rice

As shown in Table 3, the number of panicles per square meter (PPSM) under OF2 significantly exceeded that of CF, reaching the highest value (668.3 panicles m−2); OF3 and OF4 treatments also resulted in higher PPSM than CF, while OF1 resulted in a lower value. The seed-setting rate (SSR) and 1000-grain weight (1000 GW) were highest under OF4 and OF1, respectively. Both the theoretical yield and the actual yield presented the following decreasing order, OF3 > OF2 > OF4 > CF > OF1 > F0. Significant differences for OF2 and OF3 compared to CF were detected for the theoretical yield, which increased by 19.0% and 23.4%, respectively, while there was no significant difference in the actual yield among fertilization treatments. Among the treatments, the theoretical yield and actual yield associated with OF3 were the highest; this was attributed to the increase in panicle number.

3.7. Yield and Dominant Fungal Taxa in Paddy Soil: Correlation Analysis

A correlation analysis was carried out between the actual yield and the dominant fungi at the phyla and class levels to further explore the relationship between these relevant parameters (Figure 7). Dry direct-seeded rice yield was significantly positively correlated with Ascomycota relative abundance (correlation coefficient, 0.42) and significantly negatively correlated with Chytridiomycota relative abundance (correlation coefficient: −0.49), suggesting that these conditions contribute to improving rice yield. It can also be seen in Figure 7a that Ascomycota relative abundance was significantly negatively associated with that of Basidiomycota and Mortierella, and Basidiomycota, Mortierellomycota, and Chytridiomycota relative abundances also showed negative correlations with each other. These findings suggest that Ascomycota is the fungal phylum with the greatest influence on yield, and complex, competitive interactions between fungal phyla occur in the soil.
Figure 7b shows that at the class level, only the relative abundance of Sordariomycetes was significantly and positively correlated with dry direct-seeded rice yield, with a correlation coefficient of 0.41; the remaining dominant fungi were positively correlated with rice yield but at non-significant levels. The relationships among dominant fungi in the soil at the class level differed from those detected at the phylum level: the relationships in the former tended to be mutually beneficial–competitive. Very significant or significant positive correlations were detected between Sordariomycetes and Mortierellomycetes, Leotiomycetes, and Tremellomycetes. Dothideomycetes were negatively correlated with Sordariomycetes, Tremellomycetes, and Agaricomycetes. Agaricomycetes, Eurotiomycetes, and Pezizomycetes were also negatively correlated with each other.

4. Discussion

4.1. Effect of Organic-Inorganic Fertilizers on Dry Direct-Seeded Paddy Soil Fungal Diversity

It was reported that organic fertilizers could regulate soil microbial community structure, and through the changes in the soil carbon utilization pattern they induce, organic fertilizers may increase soil microbial functional diversity and improve soil quality [30,31]. However, it was also shown that combined inorganic and organic fertilizers weakly improved soil fungal diversity [32].
Soil fungal ACE and Chao indices under organic-inorganic fertilizer treatments showed a slightly decreasing trend compared with CF, without significant differences. These findings suggest that combined fertilization schemes in dry direct-seeded rice fields might not be conducive to enhancing soil fungal richness. Ji et al. [33] showed that the abundance of fungal species in plots fertilized with organic fertilizers was lower than in plots exposed to inorganic fertilizers, considering the same soil layer. Chen et al. [34] showed that the application of organic fertilizers significantly reduced fungal diversity, while inorganic fertilizers’ addition reduced fungal richness. However, a study by Luo et al. [35] showed that both inorganic and organic fertilizers increased the taxonomic abundance of soil fungi, but the effect of the organic fertilizer was more obvious.
The findings of Hu et al. [36] highlighted that geographic isolation and fertilization systems adopted under agricultural management were important factors influencing the changes in fungal communities in response to contrasting fertilization regimes. Differences in soil type and analytical methods, among others, may account for variable results [34,37]. According to the comparative analysis of Simpson and Shannon indices, biochar application (OF1) led to the highest fungal diversity, followed by the Jishiwang bioorganic fertilizer (OF3). Li et al. [38] showed that biochar application increased fungal diversity in purple soil. It was proposed that biochar has many pores, which could provide shelter for soil microorganisms; on the other hand, pores might improve soil aeration, moisture, and nutrient availability [39,40]. However, the phenols and polyphenols, microelements, organic pollutants, and persistent free radicals present in biochar might also have negative effects on microorganisms [41,42,43].
We found that the application of seaweed bioorganic fertilizer (OF2) and attapulgite organic fertilizer (OF4) resulted in lower fungal diversity than that observed under conventional fertilization, indicating that the properties of the organic fertilizer itself might be relevant in terms of fungal diversity. In a 30-year-fertilization experiment, Sun et al. [44] found that the carbon composition in the organic fertilizer may influence soil fungal diversity.
Other studies, including that of Luo et al. [45] and Chen et al. [18], found that organic fertilizers promoted the growth and metabolism of certain specific microorganisms, stimulating their growth and inhibiting other microorganisms’ multiplication, resulting in reduced microbial diversity. In addition, it has been indicated that organic fertilizers containing microorganisms might lead to an exogenous biological invasion, representing a huge threat to biodiversity and ecosystems [46]. Summing up, the influence of organic fertilizers on the fungal diversity of soils under dry direct-seeded rice cultivation still needs intensive study.

4.2. Effect of Organic-Inorganic Fertilizers on Dry Direct-Seeded Paddy Soil Fungal Structure and Function

Across all treatments, Ascomycota and Basidiomycota were the most dominant phyla, whereas Sordariomycetes was the most dominant class. These findings are consistent with the results of Li et al. [47]. Nevertheless, the relative abundance of fungal phyla and classes varied under different fertilization regimes. Compared to CF, OF1 and OF3 increased the relative abundance of soil Ascomycota but decreased the relative abundance of Basidiomycota. Oppositely, OF2 and OF4 increased the relative abundance of soil Basidiomycota but decreased the relative abundance of Ascomycota. The relative abundance of Sordariomycetes decreased in OF1 and OF2 and increased in OF3 and OF4.
Fungal communities were regulated by resource type and availability [48], and although many fungal taxa may display a wide range of resource utilization capabilities, some taxa target specific organic compounds [49]. The organic fertilizers used in this study differed in their organic matter components, and this may explain, at least in part, the impact of these fertilizers on the fungal community structure. In addition, the correlation analysis demonstrated that Ascomycota and Sordariomycetes were significantly and positively correlated with rice yield. Therefore, the treatments that increased the relative abundance of these taxa also improved rice yield.
In this study, undefined saprotroph and dung saprotroph groups occupied the dominant positions in the fungal community, accounting for more than 55% in terms of relative abundance. Notably, the relative abundance of animal pathogens was greater than that of dung saprotrophs under CF treatment, while the relative abundance of wood saprotrophs was greater than that of dung saprotrophs under OF3 treatment. Thus, the application of organic fertilizers altered the ecological functional pattern of soil fungi. Compared with CF, the inclusion of organic fertilizers increased the relative abundance of dung saprotrophs, plant pathogens, fungal parasites, and leaf saprotrophs, while the relative abundance of animal pathogens showed a significantly decreasing trend. Although some studies have shown that organic fertilizers inhibited soil pathogens [25,36,50], we found a higher relative abundance of plant pathogens under the four organic-inorganic fertilizer treatments tested, compared to CF, while the highest relative abundance of animal pathogens was detected in CF, a finding inconsistent with those of some previous studies.
Fungi have a strong nitrogen and phosphorus absorption capacity [51]; fertilizers can also nourish pathogens while promoting the growth of beneficial fungi. Organic fertilizers need to reach the stage of decomposition to allow nutrient release, and because transformation is slow, these fertilizers can provide nutrients for plants for a long time, even in the mature period; this promoting effect may also exist and account for the higher relative abundance of plant pathogens in combined organic-inorganic fertilizer treatments than conventional fertilizer alone. Furthermore, Bonanomi et al. [52] reported that the ability of organic additions to inhibit plant diseases varied greatly according to the pathogens involved. Schönning et al. [53] suggested that pathogens or harmful microorganisms in organic fertilizers may induce deleterious effects. Further research on the relationship between organic fertilizers and plant pathogens thriving in dry direct-seeded paddy soil is necessary.

5. Conclusions

The influence of different organic fertilizers on soil fungal community and rice yield in dry direct-seeded rice was explored by using six fertilizer treatments. Compared with those associated with CF, the Simpson and Shannon indices associated with biochar organic fertilizer treatment increased by 1.44% and 2.79%, respectively. The effect of organic fertilizers on fungal diversity depends on the type of organic fertilizer used. Ascomycota, Basidiomycota, and Sordariomycetes were dominant in the fungal community of dry direct-seeded rice; the application of organic fertilizer improved the community structure of soil-dominant bacteria; compared with the conventionally used NPK fertilizer, organic fertilizers (seaweed biological organic fertilizer + conventional NPK fertilizer, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, attapulgite organic fertilizer + conventional NPK fertilizer) improved the ecological function of soil fungi, and organic fertilizers treatments increased the abundance of wood saprophytic bacteria with an average increase of 22.55–76.09%. The average reduction in animal pathogens was 28.61–73.84%. The yield of rice treated with Jishiwang bioorganic fertilizer+ conventional NPK fertilizer was the highest (15.4% higher than that obtained using the conventionally used NPK fertilizer). However, many points regarding the relationship between organic fertilizers and soil fungi remain unknown and deserve further research.

Author Contributions

Conceptualization, Q.G. and Y.M.; Data curation, F.S. and D.Y.; Formal analysis, J.L.; Funding acquisition, Y.L. and G.Z.; Investigation, L.X.; Project administration, G.Z.; Supervision, Y.L.; Writing—original draft, X.G.; Writing—review and editing, Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Key Research and Development Program of China, grant number 2017YFD0300502–6, and the Heilongjiang Bayi Agricultural University Support Program for San Heng San Zong, grant number TDJH201802.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available from the author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Birla, D.S.; Malik, K.; Sainger, M.; Chaudhary, D.; Jaiwal, R.; Jaiwal, P.K. Progress and challenges in improving the nutritional quality of rice (Oryza sativa L.). Crit. Rev. Food Sci. 2017, 57, 2455–2781. [Google Scholar] [CrossRef]
  2. Tao, Y.; Chen, Q.; Peng, S.B.; Wang, W.Q.; Nie, L.X. Lower global warming potential and higher yield of wet direct-seeded rice in Central China. Agron. Sustain. Dev. 2016, 36, 24. [Google Scholar] [CrossRef] [Green Version]
  3. Xie, H.M.; Wu, K.; Liu, W.Q.; Wei, G.L.; Lu, X.; Li, Z.L.; Wei, S.Q.; Liang, H.; Jiang, L.G. Effects of partial substitution of seaweed fertilizers and microbial inoculant for chemical fertilizer on rice yield and its components. Crops 2022, 1, 161–166. (In Chinese) [Google Scholar] [CrossRef]
  4. Farooq, M.; Siddique, K.H.M.; Rehman, H.; Aziz, T.; Lee, D.J.; Wahid, A. Rice direct seeding: Experiences, challenges and opportunities. Soil Till. Res. 2011, 111, 87–98. [Google Scholar] [CrossRef]
  5. Ya, Z.J.; Lu, S.C.; Hou, K. Development status, problems and application prospects of dry direct seeding rice. Crops 2020, 2, 9–15. (In Chinese) [Google Scholar] [CrossRef]
  6. Zhang, Q.C.; Wang, G.H. Optimal nitrogen application for direct-seeding early rice. Chin. J. Rice Sci. 2002, 16, 346–350. (In Chinese) [Google Scholar] [CrossRef]
  7. Liu, H.Y.; Won, P.L.P.; Banayo, N.P.M.; Nie, L.X.; Peng, S.B.; Kato, Y. Late-season nitrogen applications improve grain yield. and fertilizer-use efficiency of dry direct-seeded rice in the tropics. Field Crops Res. 2019, 233, 114–120. [Google Scholar] [CrossRef]
  8. Dhillon, B.S.; Kaur, R.; Khanna, R.; Dhillon, B.S.; Mangat, G.S.; Gewaily, E.E.; Ghoneim, A.M.; Osman, M.A.; Kumar, P.; Choudhary, S.K.; et al. Response of rice genotypes to nitrogen fertilizer under dry direct seeding in North-Western India. Agric. Res. J. 2021, 58, 532–534. [Google Scholar] [CrossRef]
  9. Yang, A.Z.; Tang, H.Y.; Nie, Y.C.; Zhang, X. Responses of cotton growth, yield, and biomass to nitrogen split application ratio. Eur. J. Agron. 2011, 35, 164–170. [Google Scholar] [CrossRef]
  10. Wan, L.J.; He, M.; Tian, Y.; Zheng, Y.Q.; Lv, Q.; Xie, R.J.; Ma, Y.Y.; Deng, L.; Yi, S.L.; Li, J. Effects of partial substitution of chemical fertilizer with organic fertilizer on Ponkan growth and yield and soil biological properties. J. Plant Nutr. Fertil. 2022, 28, 675–687. (In Chinese) [Google Scholar] [CrossRef]
  11. Chen, Y.M.; Xiao, H.J.; Zhao, H.; Qin, S.; Gou, J.L.; Zhao, T.F.; Hu, G.; Zhao, L.X.; He, C.X.; Zhang, M.; et al. Effects of commercial organic fertilizer as substitution of chemical fertilizer on growth, yield of spring maize and soil fertility. Southwest China J. Agric. Sci. 2022, 35, 148–152. (In Chinese) [Google Scholar] [CrossRef]
  12. Klaubauf, S.; Inselsbacher, E.; Zechmeister-Boltenstern, S.; Wanek, W.; Gottsberger, R.; Strauss, J.; Gorfer, M. Molecular diversity of fungal communities in agricultural soils from Lower Austria. Fungal Divers. 2010, 44, 65–75. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Sheteiwy, M.S.; Ali, D.F.S.; Xiong, Y.C.; Brestic, M.; Skalicky, M.; Hamoud, Y.A.; Ulhassan, Z.; Shaghaleh, H.; AbdElgawad, H.; Farooq, M.; et al. Physiological and biochemical responses of soybean plants inoculated with Arbuscular mycorrhizal fungi and Bradyrhizobium under drought stress. BMC Plant Biol. 2021, 21, 195. [Google Scholar] [CrossRef] [PubMed]
  14. Sheteiwy, M.S.; Elgawad, H.A.; Xiong, Y.C.; Macovei, A.; Brestic, M.; Skalicky, M.; Shaghaleh, H.; Hamoud, Y.A.; El-Sawah, A.M. Inoculation with Bacillus amyloliquefaciens and mycorrhiza confers tolerance to drought stress and improve seed yield and quality of soybean plant. Physiol. Plant. 2021, 172, 2153–2169. [Google Scholar] [CrossRef]
  15. Zhu, Y.Z.; Li, Y.Y.; Han, J.G.; Yao, H.Y. Effects of changes in water status on soil microbes and their response mechanism a review. Chin. J. Appl. Ecol. 2019, 30, 4323–4332. (In Chinese) [Google Scholar] [CrossRef]
  16. Vries, F.T.D.; Liiri, M.E.; Bjørnlund, L.; Bowker, M.A.; Christensen, S.; Setälä, H.M.; Bardgett, R.D. Land use alters the resistance and resilience of soil food webs to drought. Nat. Clim. Change 2012, 2, 276–280. [Google Scholar] [CrossRef]
  17. Barnard, R.L.; Osborne, C.A.; Firestone, M.K. Responses of soil bacterial and fungal communities to extreme desiccation and rewetting. ISME J. 2013, 7, 2229–2241. [Google Scholar] [CrossRef]
  18. Chen, D.M.; Chen, X.M.; Liang, Y.J.; Huo, X.J.; Zhang, C.H.; Duan, Y.Q.; Yuan, L. Influence of crop rotation on soil nutrients, microbial activities and bacterial community structures. Acta Pratacult. Sin. 2015, 24, 56–65. (In Chinese) [Google Scholar] [CrossRef]
  19. Babu, A.G.; Reddy, M.S. Diversity of arbuscular mycorrhizal fungi associated with plants growing in fly ash pond and their potential role in ecological restoration. Curr. Microbiol. 2011, 63, 273–280. [Google Scholar] [CrossRef]
  20. Thormann, M.N. Diversity and function of fungi in peatlands: A carbon cycling perspective. Can. J. Soil Sci. 2006, 86, 281–293. [Google Scholar] [CrossRef]
  21. Cai, Y.; Hao, M.D.; Zang, Y.F.; He, X.Y.; Zhang, L.Q. Effect of long-term fertilization on microbial diversity of black loessial soil based on 454 Sequencing technology. J. Nucl. Agric. Sci. 2015, 29, 344–350. (In Chinese) [Google Scholar] [CrossRef]
  22. Wang, Y.; Li, J.; Cao, Z.P.; Yang, H.F. Effects of long-term fertilization on soil fungi. Chin. J. Eco-Agric. 2014, 22, 1267–1273. (In Chinese) [Google Scholar] [CrossRef]
  23. Kamaa, M.; Mburu, H.; Blanchart, E.; Chibole, L.; Chotte, J.L.; Kibunja, C.; Lesueur, D. Effects of organic and inorganic fertilization on soil bacterial and fungal microbial diversity in the Kabete long-term trial, Kenya. Biol. Fertil. Soils 2011, 47, 315–321. [Google Scholar] [CrossRef]
  24. Yang, H.; Jin, F.R.; Guan, T.W.; Xu, H.X.; Hu, X.P.; Xie, Q. Short-term effect of partial substitution of inorganic fertilizer with organic fertilizer on soil fertility and fungal communities in greenhouse. Acta Agric. Boreali-Occident. Sin. 2021, 30, 422–430. (In Chinese) [Google Scholar] [CrossRef]
  25. Lang, J.J.; Hu, J.; Ran, W.; Xu, Y.C.; Shen, Q.R. Control of cotton Verticillium wilt and fungal diversity of rhizosphere soils by bio-organic fertilizer. Biol. Fertil. Soils 2012, 48, 191–203. [Google Scholar] [CrossRef]
  26. Tao, C.Y.; Li, R.; Xiong, W.; Shen, Z.Z.; Liu, S.S.; Wang, B.B.; Ruan, Y.Z.; Kowalchuk, G.A. Bio-organic fertilizers stimulate indigenous soil Pseudomonas populations to enhance plant disease suppression. Microbiome 2020, 8, 137. [Google Scholar] [CrossRef]
  27. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  28. Bolger, A.M.; Lohse, M.; Usadel, B. Trimmomatic: A flexible trimmer for Illumina sequence data. Bioinformatics 2014, 30, 2114–2120. [Google Scholar] [CrossRef] [Green Version]
  29. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; Sahl, J.W.; et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [Green Version]
  30. Zhong, W.H.; Gu, T.; Wang, W.; Zhang, B.; Lin, X.G.; Huang, Q.R.; Shen, W.S. The effects of mineral fertilizer and organic manure on soil microbial community and diversity. Plant Soil 2010, 326, 511–522. [Google Scholar] [CrossRef]
  31. Liu, J.; Liu, M.; Wu, M.; Jiang, C.Y.; Chen, X.F.; Cai, Z.J.; Li, Z.P. Soil pH rather than nutrients drive changes in microbial community following long-term fertilization in acidic Ultisols of southern China. J. Soil Sediments 2018, 18, 1853–1864. [Google Scholar] [CrossRef]
  32. Sun, R.; Chen, Y.; Han, W.X.; Dong, W.X.; Zhang, Y.M.; Hu, C.S.; Liu, B.B.; Wang, F.H. Different contribution of species sorting and exogenous species immigration from manure to soil fungal diversity and community assemblage under long-term fertilization. Soil Biol. Biochem. 2020, 151, 108049. [Google Scholar] [CrossRef]
  33. Ji, L.F.; Ni, K.; Ma, L.F.; Chen, Z.J.; Zhao, Y.Y.; Ruan, J.Y.; Guo, S.W. Effect of different fertilizer regimes on the fungal community of acidic tea-garden soil. Acta Ecol. Sin. 2018, 38, 8158–8166. (In Chinese) [Google Scholar] [CrossRef]
  34. Chen, C.; Zhang, J.N.; Lu, M.; Qin, C.; Chen, Y.H.; Yang, L.; Shen, Q.R. Microbial communities of an arable soil treated for 8 years with organic and inorganic fertilizers. Biol. Fertil. Soils 2016, 52, 455–467. [Google Scholar] [CrossRef]
  35. Luo, P.Y.; Han, X.R.; Wang, Y.; Han, M.; Shi, H.; Liu, N.; Bai, H.Z. Influence of long-term fertilization on soil microbial biomass, dehydrogenase activity, and bacterial and fungal community structure in a brown soil of northeast China. Ann. Microbiol. 2015, 65, 533–542. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Hu, X.J.; Liu, J.J.; Wei, D.; Zhu, P.; Cui, X.A.; Zhou, B.K.; Wang, G.H. Effects of over 30-year of different fertilization regimes on fungal community compositions in the black soils of northeast China. Agric. Ecosyst. Environ. 2017, 248, 113–122. [Google Scholar] [CrossRef]
  37. Ding, J.L.; Jiang, X.; Guan, D.W.; Zhao, B.S.; Ma, M.C.; Zhou, B.K.; Li, J. Influence of inorganic fertilizer and organic manure application on fungal communities in a long-term field experiment of Chinese Mollisols. Appl. Soil Ecol. 2017, 111, 114–122. [Google Scholar] [CrossRef]
  38. Li, Y.; Yang, Y.Q.; Shen, F.; Tian, D.; Zeng, Y.M.; Yang, G.; Deng, S.H. Partitioning biochar properties to elucidate their contributions to bacterial and fungal community composition of purple soil. Sci. Total Environ. 2019, 648, 1333–1341. [Google Scholar] [CrossRef]
  39. Quilliam, R.S.; Glanville, H.C.; Wade, S.C.; Jones, D.L. Life in the ‘charosphere’-Does biochar in agricultural soil provide a significant habitat for microorganisms? Soil Biol. Biochem. 2013, 65, 287–293. [Google Scholar] [CrossRef]
  40. Sizmur, T.; Quilliam, R.; Puga, A.P.; Moreno-Jiménez, E.; Beesley, L.; Gomez-Eyles, J.L. Application of biochar for soil remediation. In: Guo MX, He ZQ, Uchimiya SM. Agricultural and environmental applications of biochar: Advances and barriers. Soil Sci. Soc. Am. Inc. 2015, 63, 295–324. [Google Scholar] [CrossRef]
  41. Graber, E.R.; Harel, Y.M.; Kolton, M.; Cytryn, E.; Silber, A.; David, D.R.; Elad, Y. Biochar impact on development and productivity of pepper and tomato grown in fertigated soilless media. Plant Soil 2010, 337, 481–496. [Google Scholar] [CrossRef]
  42. Hieu, T.; Slawo, L.; Barry, D. Potential for misidentification of environmentally persistent free radicals as molecular pollutants in particulate matter. Environ. Sci. Technol. 2010, 44, 1933–1939. [Google Scholar] [CrossRef] [Green Version]
  43. Warnock, D.D.; Mummey, D.L.; McBride, B.; Major, J.; Lehmann, J.; Rillig, M.C. Influences of non-herbaceous biochar on arbuscular mycorrhizal fungal abundances in roots and soils: Results from growth-chamber and field experiments. Appl. Soil Ecol. 2010, 46, 450–456. [Google Scholar] [CrossRef]
  44. Sun, R.; Dsouza, M.; Gilbert, J.A.; Guo, X.S.; Wang, D.Z.; Guo, Z.B.; Chu, H.Y. Fungal community composition in soils subjected to long-term chemical fertilization is most influenced by the type of organic matter. Environ. Microbiol. 2016, 18, 5137–5150. [Google Scholar] [CrossRef] [PubMed]
  45. Luo, X.Q.; Hao, X.H.; Chen, T.; Deng, C.J.; Wu, J.S.; Hu, R.G. Effects of long-term different fertilization on microbial community functional diversity in paddy soil. Acta Ecol. Sin. 2009, 29, 740–748. (In Chinese) [Google Scholar]
  46. Charles, H.; Dukes, J.S. Impacts of invasive species on ecosystem services. In Biological Invasions; Nentwig, W., Ed.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 217–237. [Google Scholar]
  47. Li, J.; Yang, X.X.; Yang, W.P.; Li, W.G.; Jing, D.D.; Yang, Z.P.; Gao, Z.Q. Effect of Forage Rape Green Manure Returning to Field on Soil Fungal Community in Winter Wheat Field. Acta Microbiol. Sin. 2021, 61, 2869–2882. (In Chinese) [Google Scholar] [CrossRef]
  48. Moll, J.; Goldmann, K.; Kramer, S.; Stefan, H.; Ellen, K.; Sven, M.; Francois, B. Resource type and availability regulate fungal communities along arable soil profiles. Microb. Ecol. 2015, 70, 390–399. [Google Scholar] [CrossRef]
  49. Hanson, C.A.; Allison, S.D.; Bradford, M.A.; Wallenstein, M.D.; Treseder, K.K. Fungal taxa target different carbon sources in forest soil. Ecosystems 2008, 11, 1157–1167. [Google Scholar] [CrossRef]
  50. Cwalina-Ambroziak, B.; Bowszys, T. Changes in fungal communities in organically fertilized soil. Plant Soil Environ. 2009, 64, 25–32. [Google Scholar] [CrossRef] [Green Version]
  51. Wu, T.H. Can ectomycorrhizal fungi circumvent the nitrogen mineralization for plant nutrition in temperate forest ecosystems? Soil Biol. Biochem. 2011, 43, 1109–1117. [Google Scholar] [CrossRef]
  52. Bonanomi, G.; Antignani, V.; Pane, C.; Scala, F. Suppression of soilborne fungal diseases with organic amendments. J. Plant Pathol. 2007, 89, 325–340. [Google Scholar]
  53. Schönning, C.; Leeming, R.; Stenström, T.A. Faecal contamination of source-separated human urine based on the content of faecal sterols. Water Res. 2002, 36, 1965–1972. [Google Scholar] [CrossRef]
Agronomy 12 01213 g001 550
Figure 1. Temperature and rainfall during dry direct-seeded rice cultivation during 2018–2020.
Figure 1. Temperature and rainfall during dry direct-seeded rice cultivation during 2018–2020.
Agronomy 12 01213 g001
Agronomy 12 01213 g002 550
Figure 2. Alpha diversity (a, ACE index; b, Chao1 index; c, Simpson index; and d, Shannon index) of soil fungal community. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer. The same letters above the bars are not significantly different among treatments (ANOVA, LSD test, p < 0.05).
Figure 2. Alpha diversity (a, ACE index; b, Chao1 index; c, Simpson index; and d, Shannon index) of soil fungal community. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer. The same letters above the bars are not significantly different among treatments (ANOVA, LSD test, p < 0.05).
Agronomy 12 01213 g002
Agronomy 12 01213 g003 550
Figure 3. Venn graph showing the number of shared and unique fungal OTUs among fertilization treatments. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Figure 3. Venn graph showing the number of shared and unique fungal OTUs among fertilization treatments. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Agronomy 12 01213 g003
Agronomy 12 01213 g004 550
Figure 4. Principal coordinates analysis of soil fungal OTUs. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Figure 4. Principal coordinates analysis of soil fungal OTUs. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Agronomy 12 01213 g004
Agronomy 12 01213 g005 550
Figure 5. Dominant fungal phyla (a) and classes (b). F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Figure 5. Dominant fungal phyla (a) and classes (b). F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Agronomy 12 01213 g005
Agronomy 12 01213 g006 550
Figure 6. Dominant ecological functional groups in the fungal community. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Figure 6. Dominant ecological functional groups in the fungal community. F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF4, attapulgite organic fertilizer + conventional NPK fertilizer.
Agronomy 12 01213 g006
Agronomy 12 01213 g007a 550Agronomy 12 01213 g007b 550
Figure 7. Correlation analysis of yield and dominant soil fungi at phylum (a) and class (b) levels. *, ** Significant at P < 0.05 and P < 0.01, respectively.
Figure 7. Correlation analysis of yield and dominant soil fungi at phylum (a) and class (b) levels. *, ** Significant at P < 0.05 and P < 0.01, respectively.
Agronomy 12 01213 g007aAgronomy 12 01213 g007b
Table
Table 1. Physiochemical properties of the soil at the experiment site between 2018 and 2020.
Table 1. Physiochemical properties of the soil at the experiment site between 2018 and 2020.
YearHydrolyzable Nitrogen
(mg kg−1)		Available Phosphorus
(mg kg−1)		Available Potassium
(mg kg−1)		Organic Matter
(g kg−1)		pH
2018106.6729.55 171.9527.705.91
2019137.4024.40159.9039.006.70
2020128.4625.38163.4340.006.32
Table
Table 2. Sequencing throughput data for soil fungi according to fertilization treatments.
Table 2. Sequencing throughput data for soil fungi according to fertilization treatments.
Treatment 1Raw ReadsEffective ReadsAverage Sequence Length Effectiveness (%)Coverage
F049,85546,12622492.500.9989
CF62,52158,37524193.370.9991
OF155,99250,55725190.280.9990
OF256,73050,94724889.800.9990
OF357,70052,90024791.660.9991
OF458,25952,59824590.280.9990
1 F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF, attapulgite organic fertilizer + conventional NPK fertilizer.
Table
Table 3. Yield and yield components of dry direct-seeded rice 1.
Table 3. Yield and yield components of dry direct-seeded rice 1.
Treatment 2PPSMGNPPSSR (%)1000 GW(g)TY (Mg ha−1)AY (Mg ha−1)
F0436.3d44.0 c74.7 b20.88 b2992 d1924 b
CF590.0bc73.9 ab89.4 a21.37 ab8304 c7372 a
OF1557.5c71.8 ab88.4 a22.87 a8075 c7334 a
OF2668.3a74.4 ab92.0 a21.61 ab9883 ab7763 a
OF3663.8a78.9 a88.5 a22.12 ab10,244 a8512 a
OF4646.3ab70.6 b99.1 a21.63 ab8890 bc7570 a
1 PPSM, panicles per square meter; GNPP, grain number per panicle; SSR, seed setting rate;1000GW,1000-grain weight; TY, theoretical yield; AY, actual yield. 2 F0, no fertilizer (F0); CF, conventional NPK fertilizer, OF1, biochar + conventional NPK fertilizer; OF2, seaweed bioorganic fertilizer + conventional NPK fertilizer; OF3, Jishiwang bioorganic fertilizer + conventional NPK fertilizer, OF, attapulgite organic fertilizer + conventional NPK fertilizer; Different small letters in the same column indicate significant differences among treatments (p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, X.; Liu, J.; Xu, L.; Sun, F.; Ma, Y.; Yin, D.; Gao, Q.; Zheng, G.; Lv, Y. Combined Organic and Inorganic Fertilization Can Enhance Dry Direct-Seeded Rice Yield by Improving Soil Fungal Community and Structure. Agronomy 2022, 12, 1213. https://doi.org/10.3390/agronomy12051213

AMA Style

Guo X, Liu J, Xu L, Sun F, Ma Y, Yin D, Gao Q, Zheng G, Lv Y. Combined Organic and Inorganic Fertilization Can Enhance Dry Direct-Seeded Rice Yield by Improving Soil Fungal Community and Structure. Agronomy. 2022; 12(5):1213. https://doi.org/10.3390/agronomy12051213

Chicago/Turabian Style

Guo, Xiaohong, Jiajun Liu, Lingqi Xu, Fujing Sun, Yuehan Ma, Dawei Yin, Qiang Gao, Guiping Zheng, and Yandong Lv. 2022. "Combined Organic and Inorganic Fertilization Can Enhance Dry Direct-Seeded Rice Yield by Improving Soil Fungal Community and Structure" Agronomy 12, no. 5: 1213. https://doi.org/10.3390/agronomy12051213

APA Style

Guo, X., Liu, J., Xu, L., Sun, F., Ma, Y., Yin, D., Gao, Q., Zheng, G., & Lv, Y. (2022). Combined Organic and Inorganic Fertilization Can Enhance Dry Direct-Seeded Rice Yield by Improving Soil Fungal Community and Structure. Agronomy, 12(5), 1213. https://doi.org/10.3390/agronomy12051213

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop