Next Article in Journal
Economic Impact and Risk Assessment of Sand and Dust Storms (SDS) on the Oil and Gas Industry in Kuwait
Next Article in Special Issue
Organic Matter Composition of Manure and Its Potential Impact on Plant Growth
Previous Article in Journal
Resignification Practices of Youth in Zona da Mata, Brazil in the Transition Toward Agroecology
Previous Article in Special Issue
The Evaluation of Hazards to Man and the Environment during the Composting of Sewage Sludge
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Short-Term Effects of Different Organic Amendments on Soil Fungal Composition

1
Key Laboratory of Sugarcane Biology and Genetic Breeding, Ministry of Agriculture, Fujian Agriculture and Forestry University, Fuzhou 350002, China
2
College of Crop Science, Fujian Agriculture and Forestry University, Fuzhou 350002, China
3
College of Geographical Sciences, Fujian Normal University, Fuzhou 350007, China
4
Graduate School of Bio-Applications and Systems Engineering, Tokyo University of Agriculture and Technology, Koganei, Tokyo 184-8588, Japan
5
Japan Collection of Microorganisms, RIKEN BioResource Center, Tsukuba, Ibaraki 305–0074, Japan
6
Fujian Provincial Key Laboratory of Agro-ecological Processing and Safety Monitoring, College of Life Sciences, Fujian Agriculture and Forestry University, Fuzhou 350002, China
7
Key Laboratory of Crop Ecology and Molecular Physiology, Fujian Agriculture and Forestry University, Fuzhou 35002, China
*
Authors to whom correspondence should be addressed.
Sustainability 2019, 11(1), 198; https://doi.org/10.3390/su11010198
Submission received: 19 October 2018 / Revised: 26 December 2018 / Accepted: 27 December 2018 / Published: 3 January 2019
(This article belongs to the Collection Organic Waste Management)

Abstract

:
Fungi play an essential role in recovering the quality and fertility of soil. There is a limited understating of the complex response of fungal diversity to different organic materials in clay loam soil. Here, we report the response of soil fungi toward the short-term application of manure (M), sugarcane straw (S), and sugarcane straw plus manure (MS), including no organic material control (CK) at two different time points (50 and 100 days after application). Illumina sequencing was used to examine the fungal communities. Our results reveal a significant shift among the soil fungal community structure associated with each organic material application. After both time points, amendments—especially M and MS—decreased the fungal richness and stimulated the copiotrophic fungal group (Ascomycota) compared to the control soil (CK) and S-amended soil. On the contrary, as compared to the M and MS-amended soils, the CK and S-amended soils at both time points increased the fungal richness and stimulated the oligotrophic fungal groups. Organic material use, especially M and MS, showed variable results regarding pathogenic fungi enhancing the abundance of Lophodermium and Cercophora and decreasing Fusarium. Concerning the abundance of plant-beneficial fungi, Mortierella was reduced, and Podospora was increased by M and MS input. FUNGuild showed that the amendment of organic materials efficiently declined the abundance of endophytes and plant pathogens, but also enhanced the animal pathogens in terms of abundance with respect to CK at two time points. This study could be useful to provide a novel understanding of the management of soil-borne pathogens by organic amendments for the sustainable production of short-term crops.

1. Introduction

Fertilization is a well-known agricultural strategy in which the mineral and organic amendments are commonly practiced to enhance soil fertility and crop yield [1]. One of the leading concerns is the growing demand for crop production to feed human beings around the globe [2]. Conventional agriculture farming based on mineral fertilizers has made a significant contribution to meeting global food demands [3]. On the contrary, soil degradation and acidification are the significant obstacles attributed to the application of mineral fertilizers; the nitrogen especially contributes to a drastic reduction in soil fertility and crop yield [4,5]. In view of that, organic material application represents an alternative approach to mineral fertilizers in order to ensure human health, agricultural sustainability, food security, and crop yield [6,7]. Furthermore, organic materials frequently improve soil fertility and structure by enhancing soil organic matter contents and nutrient status, thereby increasing soil microbial biomass and activity [8,9,10]. Additionally, different types, amounts, and frequencies of soil amendments improve the soil physiochemical attributes, affecting soil fertility, and thus ensuring crop yields [1,8].
Fungi are the primary decomposers and carbon sequesters in soil and agroecosystems [11,12,13]. The shift in fungal community composition is not only due to the input of substrates but also the existence of fungi in the organic materials [14,15,16], which may be pathogenic to human and plant health [17,18]. However, most of the researchers have studied the response of fungal communities to organic materials in term of different types in soil ecosystems under long-term experiments [19,20]. Only a few studies were carried out in order to investigate the impacts of organic materials on soil fungal community structure and soil fertility under short-term experiments [21,22], which are supposed to be very different from long-term ones. In our previous study, we explored the impacts of sugarcane straw (S), goat manure (M), and goat manure plus sugarcane straw (MS) on soil fertility, soil enzymes activities, and bacterial composition at two different time points. The results revealed that M and MS amendments could effectively improve soil fertility and enzyme activities as well as enhance the plant beneficial bacteria [8]. However, the response of soil fungal communities to M, MS, and S-treated soil was not explored. Therefore, in this study, Illumina sequencing of the fungal ITS (Internal Transcribed Spacer) region was used to investigate the response of soil fungal communities to various organic materials (M, S, and MS) at two different time points (50 and 100 days after application). We assumed that variations in the fungal communities due to the application of different organic materials may change the abundance of various fungal taxa.

2. Material and Methods

2.1. Experimental Design and Soil Sample Scheme

A pot experiment was carried out from March to June in 2017 in Fujian Agriculture and Forestry, Fuzhou, Fujian, China (latitude: 26°05′9.60″ N; longitude: 119°14′3.60″ E). In this experiment, four different treatments were performed by amending S, M, and MS to soil in the pots (red PVC) with 120-mm diameter and 180-mm height. Soil without amendment was considered as control (CK). The collection of S was carried out from the university sugarcane field, chopped into small pieces, air dried for 20 days, and kept at room temperature. M was purchased from the market. The soil had a clay loam texture (43.1% silt, 36.6% sand, and 20.3% clay) with 0.97%, 0.07%, and 0.10% of total carbon (TC), total phosphorus (TP), and total nitrogen (TN), respectively. The physiochemical characteristics of S and M that were utilized in this experiment were as follows: (i) S:TN (0.45%), TC (44.15%), total potassium (TK) (0.5%), TP (0.04%), and C:N ratio (98.11); and (ii) M:TN (0.8%), TC (17%), TK (0.41%), TP (0.52%), and C:N ratio (21.25) [8]. For MS treatment, the 7.5% M, 7.5% S, and 85% soil were combined. For S and M amendments, the 15% of S and M respectively, was combined to 85% of soil [23]. During the experiment, the relative humidity of 75–80% and temperature of 18–30 °C were kept. Three repetitions were considered for each pot. On a weekly basis, one liter of distilled water was provided to each pot [24]. Soil sampling was carried out at two time points (50 and 100 days) after the amendment of organic materials [8,25], and then samples were sieved using two mm of mesh and kept at −80 °C until DNA extraction and molecular analyses were performed. Further details and information concerning this experiment have been described in our previous work [8].

2.2. DNA Extraction, Purification, and Fungal ITS Sequencing

For each sample, total genomic DNA was extracted by using a Fast DNATM Spin kit (MP Biomedical, Santa Ana, CA, USA) and purified through DNA purification kit (Tiangen Biotech Co., Ltd., Beijing, China), according to the manufacturer’s recommendations. DNA was assessed for quantity and quality by using Nanodrop 2000c (Thermo Fisher Scientific, Middletown, VA, USA) and stored at −20 °C until sequencing. The primers ITS5-1737F and ITS2-110 2043R were used to amplify the fungal ITS1 region [26]. PCR reactions were conducted in 30-μL mixtures with each primer (0.2 μM), Phusion® High-Fidelity PCR Master Mix (15 μL) (New England BioLabs, Ipswich, MA, USA), and DNA templates (10 ng). PCRs conditions were (98 °C for one minute, followed by 30 cycles of 98 °C for 10 seconds, 50 °C for 30 seconds, and 72 °C for 60 seconds with a final extension at 72 °C for five minutes). The PCR products were purified by using QIAquik Gel Extraction Kit (QIAGEN, Düsseldorf, Germany). A TruSeq® DNA PCR-Free Sample Preparation Kit (Illumina, San Diego, CA, USA) was used to generate the sequencing libraries, and their qualities were further evaluated on the Qubit@ 2.0 Fluorometer (Thermo Scientific, Waltham, MA, USA) and Agilent Bioanalyzer 2100 system. Finally, the DNA libraries were sequenced on an Illumina HISeq2500 platform by external service (Novogen, Beijing, China).

2.3. Bioinformatic and Statistical Analysis

The paired-end reads from the original DNA fragments were merged using FLASH [27,28] and according to the unique barcodes that were assigned to each sample. Based on 97% similarity, sequences were assigned to the same operational taxonomic unit (OTU). The representative sequences were selected for each OTU and in order to annotate the taxonomic information for each representative sequence, the Ribosomal Database Project (RDP) classifier was used [29]. The diversity (Simpson and Shannon) [30] and richness (Chao1 and ACE) indexes [31,32] were used to analyze alpha diversity by Quantitative Insights into Microbial Ecology (QIIME) and displayed using R software (version 2.15.3) [33]. Based on the observed species richness, the rarefaction curves were constructed, and the Venn diagram was used in order to display the unique and common OTUs between the eight soil samples. For beta diversity analysis, unweighted UniFrac principal coordinate analysis (PCoA) and weighted Unifrac UPGMA (Unweighted Pair Group Method with Arithmetic Mean) was employed to examine the differences in species complexity between samples. In order to examine the interaction between the fungal community structure and environmental parameters, redundancy analysis (RDA) was performed on the soil’s physiochemical characteristics that were analyzed in our previous work [8]. The RDA of soil properties and the relative abundance of dominant fungal phyla was performed by using the R software (version 2.15.3). Significant differences in the relative abundance of each fungal taxon among the amendments and within each sampling time were examined through one-way analysis of variance (ANOVA) using the DPS software (version 7.05, www.dpssoftware.co.uk).

2.4. Functional Assignment

FUNGuild v1.0 [34] was used in order to explore the fungal functional groups (guilds) in organically-amended soils at the two time points. In more detail, the OTUs were generally grouped based on the three main trophic modes as saprotrophs, pathotrophs, and symbiotrophs, while all of the OTUs that were not matching with any taxa in the database were categorized as “unassigned”.

3. Results

3.1. Fungal Species Richness, Alpha and Beta Diversity Analysis

At both time points, a total of 1,897,648 (average 79,069) reads were obtained from all of the samples. According to the rarefaction analysis curve, the OTU numbers for ITS were plateaued at 97% similarity after the 45,000 sequences (Figure 1; Table S1). This confirmed that the sequence depth was sufficient to capture the richness and diversity of soil samples derived from CK, M, MS, and S at both time points. A total of 1383, 940, 1087, and 1191 OTUs were found from the samples derived from CK, M, MS, and S after the first sampling (50 days) (Figure S1). In contrast, a total of 1238, 1133, 967, and 1054 OTUs were observed from CK, M, MS, and S after the second sampling (100 days) (Figure S1). The Chao1 and ACE indexes of S, M, and MS amendments decreased compared to CK at both time points (Figure 2a,b). The Shannon and Simpson alpha diversity indexes with M, MS, and S amendments increased at the first time point as compared with CK, while at the second time point, they increased by M, decreased by MS, and appeared equal in S, with respect to CK (Figure 2c,d). The results indicate that the species richness decreased and the diversity of fungi increased by the organic material input at both sampling times. According to cluster analysis, CK and S samples at both time points were clustered in the same groups (I and III, respectively), while M and MS samples of both time points were clustered in group II (Figure 3a). PCoA confirmed what was observed with the cluster analysis, demonstrating that the composition of fungal communities might be directly associated with different organic types (Figure 3b).

3.2. Difference of Fungal Abundance by Different Organic Amendments

The most dominant phyla observed in all of the samples were Ascomycota, Basidiomycota, Chytridiomycota, Zygomycota, Glomeromycota, and Neocallimastigomycota representing the 90–96% of the fungal sequences (Figure 4; Table S2). At the first time point, in CK, the Ascomycota and Zygomycota were significantly represented as lower and higher with respect to the M, MS, and S-amended soils. At the second time point, the Ascomycota and Zygomycota were the highest and lowest in MS treatment, as compared with other treatments. At the same time point, Basidiomycota abundance in S treatment was significantly higher than the other treatments. At the first time point, Podospora abundance was significantly higher in M and MS amendments than CK, and S. Cercophora was significantly higher in M and MS with respect to S and CK amendments (Figure S2; Table S3), and the abundance of Zopfiella was higher in MS and S-amended soils rather than MS and CK. Lophodermium abundance was significantly higher in M-amended soil than that of CK, MS, and S. The Mortierella abundance was lower in M, MS, and S amendments as compared to CK. A lower abundance of Fusarium was observed in M, MS, and S amendments with respect to CK. Concerning the second time point, Podospora abundance was significantly higher in M and MS amendments than CK and S. The Cercophora abundance was higher in M and MS amendments as opposed to S and CK. A significantly higher abundance of Zopfiella was observed in S-amended soil than M, MS, and CK. As compared to CK, MS, and S amendments, Lophodermium showed higher abundance in M. Besides, a significantly higher abundance of Mortierella was observed in M-amended soil as compared to CK, MS, and S, while a lower Fusarium abundance was observed in M, MS, and S-amended soils with respect to CK.

3.3. Functional Groups

The OTUs obtained from CK, M, MS, and S treatments during both time points were detected as trophic modes, symbiotrophs, pathotrophs, and saprotrophs, while all of the OTUs that were not matching with any taxa in the database were categorized as “unassigned” (Figure S3, Table S4). Concerning the first time point, the abundance of saprotrophs in CK (54.89%) and S (61.72%) was higher than M and MS amendments (44.33% and 50.46%, respectively). At the same time, significantly lower symbiotrophs were observed in M, MS, and S (15.87%, 16.97%, and 5.07%) amendments as opposed to CK (19.47%). Besides, a significantly higher proportion of pathotrophs was observed in M (11.14%) and S (14.31%) amendments as compared to CK (8.00%) and MS (7.24%). Concerning the second time point, the proportion of saprotrophs in S (78.54%) was significantly higher than that in the CK, M, and MS treatments (52.20%, 49.21%, and 58.41%, respectively). The proportion of symbiotrophs in M and MS-amended soils (13.87% and 19.69%, respectively) was higher than that of CK (12.19%). Finally, a significantly lower proportion of pathotrophs was observed in M (8.57%), S (7.51%), and MS (4.03%) amendments compared to CK (14.87%) at the second time point.
FUNGuild gave more detailed information on the trophic modes of fungal populations obtained from the selected soil (Figure 5; Table S5). Regarding the first time point, the relative abundance of endophytes in M, MS, and S treatments (15.29%, 15.85%, and 4.31%, respectively) was reduced compared to CK (18.43%). Animal pathogens increased in the M, MS, and S treatments (1.87%, 3.32%, and 10.11%, respectively) rather than CK (0.70%). The relative abundance of plant pathogens was enhanced in M (9.26%), while the MS and S treatments (3.89% and 4.09%, respectively) reduced compared to CK (7.29%). Control treatment (CK) was dominated by unidentified saprotrophs (27.78%), which was significantly higher than M, MS, and S (8.68%, 11.15%, and 15.32%, respectively). The proportion of wood saprotrophs in CK (21.94%) was higher than that of M, MS, and S (16.23%, 15.55%, and 7.21%, respectively). Compared with CK (1.78%), the dung saprotrophs was significantly higher in M, MS, and S treatments (18.33%, 22.03%, and 34.71%, respectively).
Regarding the second time point, the relative abundance of endophytes was reduced in M-amended soils (11.13%) and more strongly reduced in the S-amended soils (4.50%), while it increased in the MS-amended soil (17.24%) compared to the CK-amended soils (11.45%). Animal pathogens increased in the M, MS, and S-amended soils (2.43%, 2.09%, and 5.49%, respectively) compared with CK (1.72%). An opposite trend was observed for plant pathogens in M, MS, and S-amended soils (respectively 6.01%, 1.89%, and 1.99%) compared with CK (13.13%). In the CK-amended soils, a high abundance of unidentified saprotrophs (29.39%) was observed with respect to the M, MS, and S-amended soils (15.45%, 12.29% and 15.63%). Wood saprotroph abundance was higher (15.96%, 24.67% and 32.20% respectively) in the M, MS, and S-amended soils as compared to CK (15.69%). Finally, a greater abundance of dung saprotrophs in the M, MS, and S-amended soils (14.51%, 19.90%, and 16.49%, respectively) was found compared to CK (1.96%).

3.4. Correlation Analysis between Soil Attributes and Fungal Community

Soil features can affect the fungal communities under different organic matter treatments. Of the analyzed soil features, C/N, NH4+, AK, TC, AP, and TN were significant conditional influences, while pH, TP, AN, and NO3 concentrations had an indicative influence on Ascomycota, Glomeromycota, Chytridiomycota, and Neocalliumastigomycota. Finally, all of these soil features were negatively correlated with Zygomycota (Figure 6).

4. Discussion

The input of organic materials into the agricultural farming system has been recognized as a promising approach in order to enhance soil fertility. In our previous research, we have demonstrated that organic materials such as S, M, and MS were more effective at improving soil fertility in a short-term experiment (Table S6) [8]. Although S, M, and MS have a positive effect on soil fertility, their addition to soil could leave destructive or beneficial effects on soil health by increasing or suppressing pathogenic and beneficial fungi. Thus, with the current research, we studied the effects of S, M, and MS, focusing on the fungal communities in amended soil at two time points.
We observe a decrease in species richness in the M and MS-amended soils as compared to the CK and S-amended ones (Figure 2a,b), and the possible explanation could be due to the high availability of substrate incorporated in the soil affecting the oligotrophic fungal groups. Particularly, the higher oligotrophic population in CK and S-amended soils could be due to the lower nutrient availability, which in turn, enhanced the species richness in these treatments [35]. The destructive and beneficial impacts of fungal taxa later organic inputs have been widely reported [35,36,37,38]. It is likely that a higher availability of nutrients by M, MS, and S amendments may result in a shift of fungal taxa abundance that may be propitious or harmful to crop, human, and animals, thus influencing agroecosystem stability and performance.
Short-term organic fertilizers could have significant influences on soil fungal community structures. Clustering and PCoA analysis revealed visible changes in the fungal community compositions due to the different amendments (Figure 3a,b), which is consistent with previous results [21,22]. Ascomycota were significantly higher in term of relative abundance in soil amended with M and MS at both time points (Figure 4; Table S2). Ascomycota is known as the most ubiquitous and diverse phyla of eukaryotes as well as the decomposition of the organic substrates (e.g., leaf litter, wood, and manure), and has been found as the predominant fungal phyla in organic-amended soils [21]. At both time points, the relative abundance of Zygomycota in M, MS, and S-amended soils was significantly reduced, which is consistent with previous findings [21,39]. At first time point, the relative abundance of Basidiomycota in the soil amended with MS and S was lower, while an opposite pattern was observed in response to the amendment with S at the second time point. However, Basidiomycota usually degrades plant litter in forests with high lignin content [40]. It is likely that the high lignin content delivered by the S has led toward the increase in their abundance in the second time point.
Pathogenic fungi commonly acquire nutrients invading the host cells, so they are considered to have an adverse influence on other members of the fungal community [41]. The symbiotrophic fungi could lead to benefits for crop quality, health, and nutrition [42,43]. At the first time point, a lower abundance of symbiotrophs was found in organically-amended soils than CK, while the reverse pattern was observed in M and MS2-amended soils at the second time point. At the first time point, a considerably higher proportion of pathotrophs was observed on M and S amended soil as compared to CK and this trend was opposed in organically-amended soils at the second time point (Figure S3; Table S4). The results of FUNGuild showed that the amendment of organic materials efficiently decreased the abundance of endophytes and plant pathogens, but also enhanced the animal pathogens in terms of the abundance compared with CK at two time points. Additionally, goat manure input not only enhanced the proportion of endophytes at the second time point, but also enhanced the plant pathogens at the first time point (Figure 5; Table S5).
When investigating the data obtained through Illumina sequencing, we considered for unique fungal taxa (beneficial or pathogenic) that are suppressed or promoted by various management approaches. Thus, we can also report on the environmental role of fungal taxa regarding their behaviors in other systems [35]. In this contest, our study contributes to understanding how the input of organic materials may increase the proliferation of pathogenic and beneficial fungal taxa. Lophodermium spp. cause needle cast disease in a wide range of plants, and it was [44], enhanced significantly in M-amended soils at the first time point compared to the second time point (Figure S2; Table S3). Cercophora were highly abundant genus in M and MS-amended soils rather than S and CK soil at both time points. Plant diseases such as leaf spots in various crops are attributed to Cercophora [45,46,47,48], which is consistent with the previous study under long-term organic amendment [35]. Conversely, Fusarium is pathogenic fungi that were highly dominant in CK-amended soils instead of organically-amended soils at both time points. Several studies have shown that the amendment of organic materials has a significant impact on controlling the plant diseases, which were caused by Fusarium sp. [49,50]. Concerning the beneficial fungi, Mortierella was highly abundant in CK instead of the organically-amended soil at the first time point. Mortierella has been known as a potential biological control agent for various plant pathogens [37], as well as a biological control agent for insects and pests [38]. Additionally, the Mortierella genus is also known as a phosphate-solubilizing fungi; they could support arbuscular mycorrhizal fungi (AMF) colonization and mitigate the damaging effects of salt on plant growth, including soil enzyme activity [36]. Members affiliated to the Podospora genus are considered to be the most abundant members of healthy soil [51], and can suppress Verticillium wilt on tomato [52]. These were highly abundant in M and MS-amended soils at both time points. Furthermore, the Zopfiella genus was enhanced in S-amended soil at the first time point than a second time point. Zopfiella ssp. survives in healthy soils and produces antifungal compounds that control the disease of various crops caused by soil-borne pathogens [51,52,53]. These results demonstrated that fungal communities and functions render differently to the short-term amendment of different organic materials in soils at two different time points. These results are useful in providing a novel understanding in order to manage soil-borne pathogens by organic amendments for the sustainable production of short-term crops.

5. Conclusions

Overall, our short-term observations exposed that M and MS amendments increased soil pH and nutrient availability, especially C, N, and P as compared to CK and S amendments at both time points, which in turn reduced species’ richness and affected the proliferation of some genera that are plant pathogenic or beneficial. These changes induced by organic amendment may affect the stability and productivity of the agroecosystem. Soil fungi respond differently to the inputs of organic fertilizers under short-term experiments, which offers novel insights into the potential of managing soil-borne and soil-beneficial fungi for sustainable agricultural productivity. More research on field levels is needed in order to study the quick response of soil pathogenic and beneficial fungi to different organic amendments in the short-term cropping system.

Supplementary Materials

The following are available online at https://www.mdpi.com/2071-1050/11/1/198/s1, Figure S1: Venn diagram for soil samples at two time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days, Figure S2: Top 10 Fungi with high relative abundance at genus level in the soil samples. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days, Figure S3: Variation in fungal function at different time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days, Table S1: Effective tags from the soil samples derived after organic amendment in both time points, Table S2: Composition of different fungal phyla in soil samples derived after different organic amendments at both time points, Table S3: Top fungi with high relative abundance at genus level in the soil samples, Table S4: Variation in fungal function under organic-amended soils at different time points, Table S5: Shift in structure of fungal functional groups (guilds) inferred by FUNGuild at two different time points, Table S6: Physiochemical characteristics of the soil.

Author Contributions

All authors contributed to the intellectual input and provided assistance to this study and manuscript preparation. M.T. and W.I. designed the research and conducted the experiments. M.T., F.K., and Z.P. analyzed the data. M.T. wrote the manuscript. C.G.L., H.Z., W.L., and S.L. reviewed the manuscript. W.L. and H.Z. supervised the work and approved the manuscript for publication.

Funding

This research was funded by the Public Welfare Industrial Research Project of China grant number 20150311904-01 and the Earmarked Fund for the Modern Agriculture Technology of China via grant number CARS-20-4-3.

Acknowledgments

We acknowledge the Public Welfare Industrial Research Project of China, (20150311904-01) and the Earmarked Fund for the Modern Agriculture Technology of China (CARS-20-4-3). We want to praise the donor agencies of China for supporting this project.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Saha, S.; Gopinath, K.A.; Mina, B.L.; Gupta, H.S. Influence of continuous application of inorganic nutrients to a Maize-Wheat rotation on soil enzyme activity and grain quality in a rainfed Indian soil. Eur. J. Soil Biol. 2008, 44, 521–531. [Google Scholar] [CrossRef]
  2. Foley, J.A.; Ramankutty, N.; Brauman, K.A.; Cassidy, E.S.; Gerber, J.S.; Johnston, M.; Mueller, N.D.; O’Connell, C.; Ray, D.K.; West, P.C.; et al. Solutions for a cultivated planet. Nature 2011, 478, 337–342. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Seufert, V.; Ramankutty, N.; Foley, J.A. Comparing the yields of organic and conventional agriculture. Nature 2012, 485, 229–232. [Google Scholar] [CrossRef] [PubMed]
  4. Barak, P.; Jobe, B.O.; Krueger, A.R.; Peterson, L.A.; Laird, D.A. Effects of long-term soil acidification due to nitrogen fertilizer inputs in Wisconsin. Plant Soil 1997, 197, 61–69. [Google Scholar] [CrossRef]
  5. Ju, X.-T.; Xing, G.-X.; Chen, X.-P.; Zhang, S.-L.; Zhang, L.-J.; Liu, X.-J.; Cui, Z.-L.; Yin, B.; Christie, P.; Zhu, Z.-L. Reducing environmental risk by improving N management in intensive Chinese agricultural systems. Proc. Natl. Acad. Sci. USA 2009, 106, 3041–3046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Francioli, D.; Schulz, E.; Lentendu, G.; Wubet, T.; Buscot, F.; Reitz, T. Mineral vs. organic amendments: Microbial community structure, activity and abundance of agriculturally relevant microbes are driven by long-term fertilization strategies. Front. Microbiol. 2016, 7. [Google Scholar] [CrossRef] [PubMed]
  7. Das, S.; Jeong, S.T.; Das, S.; Kim, P.J. Composted cattle manure increases microbial activity and soil fertility more than composted swine manure in a submerged rice paddy. Front. Microbiol. 2017, 8. [Google Scholar] [CrossRef]
  8. Tayyab, M.; Islam, W.; Arafat, Y.; Pang, Z.; Zhang, C.; Lin, Y.; Waqas, M.; Lin, S.; Lin, W.; Zhang, H. Effect of Sugarcane Straw and Goat Manure on Soil Nutrient Transformation and Bacterial Communities. Sustainability 2018, 10, 2361. [Google Scholar] [CrossRef]
  9. Santonja, M.; Rancon, A.; Fromin, N.; Baldy, V.; Hättenschwiler, S.; Fernandez, C.; Montès, N.; Mirleau, P. Plant litter diversity increases microbial abundance, fungal diversity, and carbon and nitrogen cycling in a Mediterranean shrubland. Soil Biol. Biochem. 2017, 111, 124–134. [Google Scholar] [CrossRef] [Green Version]
  10. Liang, Q.; Chen, H.; Gong, Y.; Fan, M.; Yang, H.; Lal, R.; Kuzyakov, Y. Effects of 15 years of manure and inorganic fertilizers on soil organic carbon fractions in a wheat-maize system in the North China Plain. Nutr. Cycl. Agroecosyst. 2012, 92, 21–33. [Google Scholar] [CrossRef]
  11. Hoorman, J.J. The role of soil fungus. In Agriculture and Natural Resources; SAG-14-11; Ohio State University Extension: Columbus, OH, USA, 2016. [Google Scholar]
  12. Jones, M.D.M.; Forn, I.; Gadelha, C.; Egan, M.J.; Bass, D.; Massana, R.; Richards, T.A. Discovery of novel intermediate forms redefines the fungal tree of life. Nature 2011, 474, 200–204. [Google Scholar] [CrossRef] [PubMed]
  13. Tedersoo, L.; Bahram, M.; Põlme, S.; Kõljalg, U.; Yorou, N.S.; Wijesundera, R.; Ruiz, L.V.; Vasco-Palacios, A.M.; Thu, P.Q.; Suija, A.; Smith, M.E.; Sharp, C.; et al. Global diversity and geography of soil fungi. Science 2014, 346. [Google Scholar] [CrossRef]
  14. Powell, K.I.; Chase, J.M.; Knight, T.M. Invasive plants have scale-dependent effects on diversity by altering species-area relationships. Science 2013, 339, 316–318. [Google Scholar] [CrossRef]
  15. Ziska, L.H.; Blumenthal, D.M.; Runion, G.B.; Hunt, E.R.; Diaz-Soltero, H. Invasive species and climate change: An agronomic perspective. Clim. Chang. 2011, 105, 13–42. [Google Scholar] [CrossRef]
  16. Ehrenfeld, J.G. Ecosystem Consequences of Biological Invasions. Annu. Rev. Ecol. Evol. Syst. 2010, 41, 59–80. [Google Scholar] [CrossRef]
  17. Bradford, S.A.; Morales, V.L.; Zhang, W.; Harvey, R.W.; Packman, A.I.; Mohanram, A.; Welty, C. Transport and fate of microbial pathogens in agricultural settings. Crit. Rev. Environ. Sci. Technol. 2013, 43, 775–893. [Google Scholar] [CrossRef]
  18. Mosaddeghi, M.R.; Sinegani, A.A.S.; Farhangi, M.B.; Mahboubi, A.A.; Unc, A. Saturated and unsaturated transport of cow manure-borne Escherichia coli through in situ clay loam lysimeters. Agric. Ecosyst. Environ. 2010, 137, 163–171. [Google Scholar] [CrossRef]
  19. Koivunen, E.E.; Tully, K.L.; Swett, C.L. Co-managing soil and plant pathogens: Effects of organic amendments on soil fertility and fungal pathogen survival. Plant Soil 2018, 432, 171–189. [Google Scholar] [CrossRef]
  20. Sun, R.; Dsouza, M.; Gilbert, J.A.; Guo, X.; Wang, D.; Guo, Z.; Ni, Y.; Chu, H. Fungal community composition in sotterils subjected to long-term chemical fertilization is most influenced by the type of organic ma. Environ. Microbiol. 2016, 18, 5137–5150. [Google Scholar] [CrossRef]
  21. Guo, J.; Liu, W.; Zhu, C.; Luo, G. Bacterial rather than fungal community composition is associated with microbial activities and nutrient-use efficiencies in a paddy soil with short-term organic amendments. Plant Soil 2018, 424, 335–349. [Google Scholar] [CrossRef]
  22. Leite, M.F.A.; Pan, Y.; Bloem, J.; Ten Berge, H.; Kuramae, E.E. Organic nitrogen rearranges both structure and activity of the soil-borne microbial seedbank. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [PubMed]
  23. Sun, J.; Zhang, Q.; Zhou, J.; Wei, Q. Pyrosequencing technology reveals the impact of different manure doses on the bacterial community in apple rhizosphere soil. Appl. Soil Ecol. 2014, 78, 28–36. [Google Scholar] [CrossRef]
  24. Sun, H.; Wang, Q.; Liu, N.; Li, L.; Zhang, C.; Liu, Z.; Zhang, Y. Effects of different leaf litters on the physicochemical properties and bacterial communities in Panax ginseng-growing soil. Appl. Soil Ecol. 2017, 111, 17–24. [Google Scholar] [CrossRef]
  25. Stark, C.; Condron, L.M.; Stewart, A.; Di, H.J.; O’Callaghan, M. Influence of organic and mineral amendments on microbial soil properties and processes. Appl. Soil Ecol. 2007, 35, 79–93. [Google Scholar] [CrossRef] [Green Version]
  26. Huang, Y.; Kuang, Z.; Wang, W.; Cao, L. Exploring potential bacterial and fungal biocontrol agents transmitted from seeds to sprouts of wheat. Biol. Control 2016, 98, 27–33. [Google Scholar] [CrossRef]
  27. Tan, Y.; Cui, Y.; Li, H.; Kuang, A.; Li, X.; Wei, Y.; Ji, X. Rhizospheric soil and root endogenous fungal diversity and composition in response to continuous Panax notoginseng cropping practices. Microbiol. Res. 2017, 194, 10–19. [Google Scholar] [CrossRef]
  28. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  29. Wang, Q.; Garrity, G.M.; Tiedje, J.M.; Cole, J.R. Naive Bayesian classifier for rapid assignment of rRNA\nsequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 2007, 73, 5261–5267. [Google Scholar] [CrossRef] [PubMed]
  30. Keylock, C.J. Simpson diversity and the Shannon-Wiener index as special cases of a generalized entropy. Oikos 2005, 109, 203–207. [Google Scholar] [CrossRef]
  31. Chao, A. Non-parametric estimation of the number of classes in a population. Ann. Math. Stat. 1984, 20, 265–270. [Google Scholar]
  32. Chao, A.; Lee, S.M. Estimating the number of classes via sample coverage. J. Am. Stat. Assoc. 1992, 87, 210–217. [Google Scholar] [CrossRef]
  33. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2014. [Google Scholar]
  34. Nguyen, N.H.; Song, Z.; Bates, S.T.; Branco, S.; Tedersoo, L.; Menke, J.; Schilling, J.S.; Kennedy, P.G. FUNGuild: An open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecol. 2016, 20, 241–248. [Google Scholar] [CrossRef]
  35. Hartmann, M.; Frey, B.; Mayer, J.; Mäder, P.; Widmer, F.; Ma, P. Distinct soil microbial diversity under long-term organic and conventional farming. ISME J. 2015, 9, 1177–1194. [Google Scholar] [CrossRef] [PubMed]
  36. Zhang, H.; Wu, X.; Li, G.; Qin, P. Interactions between arbuscular mycorrhizal fungi and phosphate-solubilizing fungus (Mortierella sp.) and their effects on Kostelelzkya virginica growth and enzyme activities of rhizosphere and bulk soils at different salinities. Biol. Fert. Soils 2011, 47, 543–554. [Google Scholar] [CrossRef]
  37. Tagawa, M.; Tamaki, H.; Manome, A.; Koyama, O.; Kamagata, Y. Isolation and characterization of antagonistic fungi against potato scab pathogens from potato field soils. FEMS Microbiol. Lett. 2010, 305, 136–142. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Edgington, S.; Thompson, E.; Moore, D.; Hughes, K.A.; Bridge, P. Investigating the insecticidal potential of Geomyces (Myxotrichaceae: Helotiales) and Mortierella (Mortierellacea: Mortierellales) isolated from Antarctica. SpringerPlus 2014, 3, 1–8. [Google Scholar] [CrossRef] [PubMed]
  39. Wang, J.; Song, Y.; Ma, T.; Raza, W.; Li, J.; Howland, J.G.; Huang, Q.; Shen, Q. Impacts of inorganic and organic fertilization treatments on bacterial and fungal communities in a paddy soil. Appl. Soil Ecol. 2017, 112, 42–50. [Google Scholar] [CrossRef]
  40. Blackwood, C.B.; Waldrop, M.P.; Zak, D.R.; Sinsabaugh, R.L. Molecular analysis of fungal communities and laccase genes in decomposing litter reveals differences among forest types but no impact of nitrogen deposition. Environ. Microbiol. 2007, 9, 1306–1316. [Google Scholar] [CrossRef]
  41. Anthony, M.A.; Frey, S.D.; Stinson, K.A. Fungal community homogenization, shift in dominant trophic guild, and appearance of novel taxa with biotic invasion. Ecosphere 2017, 8. [Google Scholar] [CrossRef]
  42. Igiehon, N.O.; Babalola, O.O. Biofertilizers and sustainable agriculture: Exploring arbuscular mycorrhizal fungi. Appl. Microbiol. Biotechnol. 2017, 101, 4871–4881. [Google Scholar] [CrossRef]
  43. Rouphael, Y.; Franken, P.; Schneider, C.; Schwarz, D.; Giovannetti, M.; Agnolucci, M.; De Pascale, S.; Bonini, P.; Colla, G. Arbuscular mycorrhizal fungi act as biostimulants in horticultural crops. Sci. Horticult. 2015, 196, 91–108. [Google Scholar] [CrossRef]
  44. Agarwal, G.P.; Singh, S.M. Microascus cinereus infection of human nail. Indian J. Med. Sci. 1980, 34, 263–265. [Google Scholar] [PubMed]
  45. Kimber, R.B.E.; Paull, J.G.; Scott, E.S.; Dyson, C.B.; Davidson, J.A. Temporal and Spatial Development of Cercospora Leaf Spot of Faba Bean Influenced by In Situ Inoculum. Plant Dis. 2016, 100, 1823–1830. [Google Scholar] [CrossRef]
  46. Smith, D.H. Management of Peanut Foliar Dis. Plant Dis. 1980, 64, 357. [Google Scholar] [CrossRef]
  47. Balau, A.M.; Faretra, F. Molecular method for detection of Cercospora beticola Sacc. Lucra Ştiinţifice Univ. Stiinte Agric. Med. Vet. 2010, 53, 181–183. [Google Scholar]
  48. Skaracis, G.; Pavli, O.; Biancardi, E. Cercospora Leaf Spot Disease of Sugar Beet. Sugar Tech. 2010, 12, 220–228. [Google Scholar] [CrossRef]
  49. Bonanomi, G.; Antignani, V.; Pane, C.; Scala, F. Suppression of soilborne fungal diseases with organic amendments. J. Plant Pathol. 2007, 89, 311–324. [Google Scholar]
  50. Xiong, W.; Guo, S.; Jousset, A.; Zhao, Q.; Wu, H. Soil Biology & Biochemistry Bio-fertilizer application induces soil suppressiveness against Fusarium wilt disease by reshaping the soil microbiome. Soil Biol. Biochem. 2017, 114, 238–247. [Google Scholar] [CrossRef]
  51. Xu, L.; Ravnskov, S.; Larsen, J.; Nilsson, R.H.; Nicolaisen, M. Soil fungal community structure along a soil health gradient in pea fields examined using deep amplicon sequencing. Soil Biol. Biochem. 2012, 46, 26–32. [Google Scholar] [CrossRef]
  52. Dutta, B.K. Studies on some fungi isolated from the rhizosphere of tomato plants and the consequent prospect for the control of Verticillium wilt. Plant Soil 1981, 63, 209–216. [Google Scholar] [CrossRef]
  53. Huang, X.; Liu, L.; Wen, T.; Zhu, R.; Zhang, J.; Cai, Z. Illumina MiSeq investigations on the changes of microbial community in the Fusarium oxysporum f.sp. cubense infected soil during and after reductive soil disinfestation. Microbiol. Res. 2015, 181, 33–42. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Rarefaction curves for organic-amended soil samples of two time points at an operational taxonomic unit (OTU) threshold of 97% sequence similarity. CK1.1, CK1.2, and CK1.3: control soil at 50 days; M1.1, M1.2, and M1.3: goat manure-treated soil at 50 days; MS1.1, MS1.2, and MS1.3: goat manure plus sugarcane straw-treated soil at 50 days; S1.1, S1.2, and S1.3: sugarcane straw-treated soil at 50 days. CK2.1, CK2.2, and CK2.3: control soil at 100 days; M2.1, M2.2, and M2.3: goat manure-treated soil at 100 days; MS2.1, MS2.2, and MS2.3: goat manure plus sugarcane straw-treated soil at 100 days; S1.1, S1.2, and S1.3: sugarcane straw-treated soil at 100 days.
Figure 1. Rarefaction curves for organic-amended soil samples of two time points at an operational taxonomic unit (OTU) threshold of 97% sequence similarity. CK1.1, CK1.2, and CK1.3: control soil at 50 days; M1.1, M1.2, and M1.3: goat manure-treated soil at 50 days; MS1.1, MS1.2, and MS1.3: goat manure plus sugarcane straw-treated soil at 50 days; S1.1, S1.2, and S1.3: sugarcane straw-treated soil at 50 days. CK2.1, CK2.2, and CK2.3: control soil at 100 days; M2.1, M2.2, and M2.3: goat manure-treated soil at 100 days; MS2.1, MS2.2, and MS2.3: goat manure plus sugarcane straw-treated soil at 100 days; S1.1, S1.2, and S1.3: sugarcane straw-treated soil at 100 days.
Sustainability 11 00198 g001
Figure 2. Box plots of the Chao1, ACE, Shannon and Simpson diversity indices of soil samples derived by various organic amendments at both time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Figure 2. Box plots of the Chao1, ACE, Shannon and Simpson diversity indices of soil samples derived by various organic amendments at both time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Sustainability 11 00198 g002
Figure 3. (a) Weighted Unifrac UPGMA clustering analysis of fungal communities associated with soil samples derived by various organic-amended soils at both time points. (b) Unweighted UniFrac (UUF) principal coordinate analysis (PCoA) plot based on the ITS sequencing genes from soil samples of both time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Figure 3. (a) Weighted Unifrac UPGMA clustering analysis of fungal communities associated with soil samples derived by various organic-amended soils at both time points. (b) Unweighted UniFrac (UUF) principal coordinate analysis (PCoA) plot based on the ITS sequencing genes from soil samples of both time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Sustainability 11 00198 g003
Figure 4. Fungi with high relative abundance at the phylum level in the soil samples derived from control and organic-amended soils at both time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Figure 4. Fungi with high relative abundance at the phylum level in the soil samples derived from control and organic-amended soils at both time points. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Sustainability 11 00198 g004
Figure 5. Variation in fungal function and the structure of the fungal functional groups (guilds) observed in this study. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Figure 5. Variation in fungal function and the structure of the fungal functional groups (guilds) observed in this study. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Sustainability 11 00198 g005
Figure 6. Redundancy analysis of soil properties and the relative abundance of dominant fungal phyla. These factors such as C/N, NH4+, AK, TC, AP, and TN, pH, AN, TP, and NO3 concentrations are shown in blue and red text. Soil properties’ arrows in terms of length are showing the strength of interaction among the overall fungal community and soil properties. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Figure 6. Redundancy analysis of soil properties and the relative abundance of dominant fungal phyla. These factors such as C/N, NH4+, AK, TC, AP, and TN, pH, AN, TP, and NO3 concentrations are shown in blue and red text. Soil properties’ arrows in terms of length are showing the strength of interaction among the overall fungal community and soil properties. CK1: control soil at 50 days; M1: goat manure-amended soil at 50 days; MS1: goat manure plus sugarcane straw-amended soil at 50 days; S1: sugarcane straw-amended soil at 50 days. CK2: control soil at 100 days; M2: goat manure-amended soil at 100 days; MS2: goat manure plus sugarcane straw-amended soil at 100 days; S2: sugarcane straw-amended soil at 100 days.
Sustainability 11 00198 g006

Share and Cite

MDPI and ACS Style

Tayyab, M.; Islam, W.; Lee, C.G.; Pang, Z.; Khalil, F.; Lin, S.; Lin, W.; Zhang, H. Short-Term Effects of Different Organic Amendments on Soil Fungal Composition. Sustainability 2019, 11, 198. https://doi.org/10.3390/su11010198

AMA Style

Tayyab M, Islam W, Lee CG, Pang Z, Khalil F, Lin S, Lin W, Zhang H. Short-Term Effects of Different Organic Amendments on Soil Fungal Composition. Sustainability. 2019; 11(1):198. https://doi.org/10.3390/su11010198

Chicago/Turabian Style

Tayyab, Muhammad, Waqar Islam, Chol Gyu Lee, Ziqin Pang, Farghama Khalil, Sheng Lin, Wenxiong Lin, and Hua Zhang. 2019. "Short-Term Effects of Different Organic Amendments on Soil Fungal Composition" Sustainability 11, no. 1: 198. https://doi.org/10.3390/su11010198

APA Style

Tayyab, M., Islam, W., Lee, C. G., Pang, Z., Khalil, F., Lin, S., Lin, W., & Zhang, H. (2019). Short-Term Effects of Different Organic Amendments on Soil Fungal Composition. Sustainability, 11(1), 198. https://doi.org/10.3390/su11010198

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