This study aimed to compare monoculture and mixed rubber plantations in terms of their soil bacterial and fungal composition. An Illumina MiSeq sequencing analysis was performed to investigate the composition and diversity of the soil bacterial and fungal communities among three different rubber (Hevea brasiliensis) plantations: monoculture, Mixed 1 (Hevea brasiliensis and Mytilaria laosensis), and Mixed 2 (Hevea brasiliensis and Michelia macclurei) in Hainan. The results showed that the bacterial composition of the three rubber plantations was basically similar. However, there was a significant difference in fungal communities among the three rubber plantations at both the phylum and operational taxonomic unit level. The species richness, Chao, and Shannon diversity of bacterial communities of monoculture rubber plantations were higher than the Mixed 1 and Mixed 2 rubber plantations, whereas all diversity indexes of fungal communities were relatively equal for the monoculture and mixed rubber plantations. Soil nutrition (such as total nitrogen and total potassium) and soil pH are the main drivers of the bacterial composition (p <.001). However, soil pH and water content are the main drivers of the fungal composition (p <.001), and to some extent, soil pH can increase soil bacteria diversity. We suggest that alkaline fertilizers should be applied in mixed rubber plantations to improve the soil pH and, consequently, to increase the total diversity of the rubber plantation.
Introduction
Rubber (Hevea brasiliensis Muell. Arg) is an important economic crop in the world, and its cultivation is one of the main sources of income in tropical region of Asia. During the last two decades, rubber plantations (RPs) seem to have extended quickly throughout Asia (Li et al., 2015), and a large area of natural forests have been changed into RPs (Chen et al., 2017; Yan et al., 2015). Deforestation influences local temperature, light, dampness, and litter conditions, bringing about changes of microhabitats and loss of forests (Chaudhary, Burivalova, Koh, & Hellweg, 2016). The conservation of tropical forests has been found to have global importance (Cole, Bhagwat, & Willis, 2014; Naeem, Duffy, & Zavaleta, 2012). Currently, southeast Asian croplands are under pressure due to increasing population and demands of crop products in the world (Koh, Miettinen, Liew, & Ghazoul, 2011; Sodhi, Koh, Brook, & Ng, 2004; Wilcove, Giam, Edwards, Fisher, & Koh, 2013). This has aroused the interest of biodiversity scholars around the word. Hainan is considered to be the largest island in the Indo-Burma biodiversity hotspot and harbors many endangered species (Myers, Mittermeier, Mittermeier, Fonseca, & Kent, 2000). Now, RPs compose one fourth of the vegetation on Hainan Island (Wang, Lan, Wu, & Xie, 2012), and they are one of the most important vegetation compositions in Hainan Province (Lan, Wang, Wu, & Xie, 2013). Studies have shown that cropping systems involving the removal of an agricultural byproduct can reduce soil organic carbon (Singh, Rai, & Singh, 2016). Moreover, forest conversion also affects the composition of soil microorganisms (Cai, Zhang, Yang, & Wang, 2018; Rodrigues et al., 2013).
Microbial activity plays a key role in nutrient cycling and amelioration of plant stresses and is responsible for a broad range of the ecological functions of soil (Singh et al., 2016). The soil microbial diversity is a critical marker of soil quality (He, Zheng, Chen, He, & Zhang, 2008; Kirk et al., 2004). It was reported that soil microbial diversity is affected by a variety of environmental factors (Balser, Kinzig, & Firestone, 2002; Schimel, 1995). Previous findings revealed that abiotic soil nutrients drive the soil bacterial communities on RPs in Hainan (Lan, Li, Wu, & Xie, 2017a). Soil pH is another important factor that can affect the soil bacterial community structure (Krashevska, Klarner, Widyastuti, Maraun, & Scheu, 2015). Seasonal changes are also main drivers of the bacterial composition of RPs and other tropical forests (Lan, Li, Lesueur, Wu, & Xie, 2018). In addition, other factors controlling the distribution of soil microorganisms are being studied (Rousk et al., 2010). For instance, it is important to know whether different planting patterns of RP (monoculture and mixed plantations) affect the composition and diversity of soil microbes.
Monoculture has often been criticized for reducing biodiversity, causing widespread outbreaks of pests and diseases (Erskine, Lamb, & Bristow, 2006; Lamb, Erskine, & Parrotta, 2005) and not supporting numerous conventional ecosystem services (Singh, Rai, Banyal, Chauhan, & Singh, 2018). Thus, there is an emphasis on growing mixed plantations rather than monoculture. This includes an expanded yield, better natural assurance, biodiversity protection, and remodeling (Montagnini, Gonzales, & Porras, 1995). We aimed to test if there were any differences in the composition and diversity of soil microorganisms between monoculture and mixed RPs. In this study, we examined the bacterial and fungal composition and diversity of monoculture and mixed RPs by utilizing an extensive throughput Illumina MiSeq sequencing analysis. The objectives of this study were (a) to identify the soil microbial composition and changes in community structure among monoculture, Mixed 1 and Mixed 2 plantations and (b) to understand the main factors that affect the microbial community. We aimed to test whether soil microbial diversity in mixed RPs was higher than in a pure RP.
Methods
Study Site
Hainan Province (18°10′ −20°10′N and 108°37′ −111°03′E), with a territory of 33,920 km2, is located in the south of China. Hainan Island is one of the most important hot spots of biodiversity conservation in the Indo-Burma region. Hainan Island’s tropical rainforests are situated at the northern edge of tropical Asia (Zhu & Zhou, 2002) and are recognized for their high biodiversity (Zheng, Cour, & Zou, 2002). This island has a tropical rainstorm environment with a wet period from May to October and a dry period from November to April. The annual precipitation was recorded more than 1,600 mm and the precipitation was irregularly distributed between wet and dry seasons (Francisco-Ortega et al., 2010). The lowest mean monthly temperature was recorded in January (16°C) and the highest in July 27.5°C. Precipitation frequently occurs on the eastern part of Hainan Island because of storms in the Pacific Ocean, while the western part of the island is dry (Xu, Tan, Wang, & Gai, 2016).
Sample Collection
The study site was located in Wenchang area, northeast of Hainan Island. We selected three plots (20 m × 20 m) on three RPs, that is, a monoculture RP, a Mixed 1 RP (Hevea brasiliensis and Mytilaria laosensis), and a Mixed 2 RP (Michelia macclurei). Soil sampling design in three rubber plantations were shown in Figure A1. Monoculture RPs have understory herbaceous species such as Parabaena sagittata, Cyrtococcum patens, Eupatorium odoratum, and Mimosa pudica. Because of the high canopy density, the coverage of the herbaceous plants in Mixed 1 is relatively lower than in the monoculture. Mixed 1 comprised Cyrtococcum patens and Ageratum conyzoides. The coverage of Mixed 2 comprised Cyrtococcum patens and Eupatorium odoratum.
Soil samples were randomly collected from each RP. Six soil samples were randomly gathered within each plot by using a 5 cm diameter stainless steel cylinder. The soil samples were then separated into three sections; the first was sieved throughout a 2 mm net instantly and stored at 4°C until analysis of MBC; the second one was air-dried and passed through a 0.25 mm sieve for SOM, soil pH, and TN; and the third one remained stored at −80°C for DNA extraction. The soil pH was estimated in a 1:1 soil:water mixture. The soil moisture was determined gravimetrically. Total phosphorus (TP) was digested by molybdenum-antimony, and measurement was performed by spectrophotometric method (Soil Science Society of China, 2000). Total potassium (TK) content was measured by flame photometry. MBC was examined by chloroform fumigation and an extraction technique utilizing the correction factor 0.35 (Shen et al., 2014). The forests and soil properties are shown in Table 1.
Table 1.
Soils Characteristics of the Three Rubber Plantations in Wenchang, Hainan China.
DNA Extraction, PCR Amplification, and Illumina MiSeq Sequencing
Microbial DNA was extracted from 5.0 g of soil using the E.Z.N.A.® Soil DNA Kit (Omega Biotek, Norcross, GA, USA) according to the manufacturer’s protocols. For bacteria, the V3–V4 hypervariable regions of the bacterial 16S rRNA gene were amplified with primers 338 F (5′-ACTCCTACGGGAGGCAGCAG-3′) and 806 R (5′-GGACTACHVGGGTWTCTAAT-3′) by a thermocycler polymerase chain reaction system (Xu et al., 2016; GeneAmp 9700, ABI, USA). Fungal Internal Transcribed Spacer (ITS) Region 1 was amplified using the ITS1F (5′-CTTGGTCATTTAGAGG AAGTAA-3′) and ITS2 (5′-GCTGCGTTCTTCATCG ATGC-3′) primer pair (Kerfahi, Tripathi, Dong, Go, & Adams, 2016). Polymerase chain reactions (PCRs) were conducted in triplicate and the tests were conducted using a 20 µL blend consisting of 4 µL of 5 × FastPfu Buffer, 2 µL of 2.5 m dNTPs, 0.8 µL of the sample each primer (5 µM), 0.4 µL of FastPfu Polymerase, and 10 ng of template DNA. Amplicons were separated on 2% agarose gels and purified utilizing the AxyPrep DNA gel extraction kit (Axygen Biosciences, Union City, CA, USA) as indicated by the manufacturer’s instructions and measured utilizing QuantiFluor™—ST (Promega, USA). Purged amplicons were mutually equimolar and combined end-sequenced (2 × 250) on an Illumina MiSeq stage was conducted as indicated by standard conventions.
Bioinformatics and Statistical Analyses
Basic FASTQ records were demultiplexed and separated utilizing QIIME (version 1.17). The adjusted groupings were collected into operational taxonomic units (OTUs) characterized through 97% comparability (Stackebrandt & Goebel, 1994) utilizing the CD-HIT-OTU driver (Wu, Zhu, Fu, Niu, & Li, 2011). The phylogenetic attachments of all 16S and 18S rRNA gene arrangements were examined by the Ribosomal Database Project classifier of the SILVA (SSU117/119) 16S and 18S rRNA catalog utilizing a confidence threshold of 70%. We calculated the exposure proportion using Good’s technique. The Shannon diversity index and nonparametric measures of genus richness (Chao1) were designed for each illustration using Mothur (Schloss et al., 2009). Principal coordinate analysis (PCoA) based on the Bray-Curtis distance matrix of family composition was achieved. To check the microbial community’s similarity index, the similarities were statistically analyzed to test whether there was any significant variance among the RPs.
The seven ecological variables included in the analysis were the type of RP (1: monoculture RP; 2: Mixed 1 RP; 3: Mixed 2 RP), soil organic matter (SOM), TP, total nitrogen (TN), TK, water content (WC), and soil pH. To identify relationship between microbial communities and ecological factors, a redundancy analysis was conducted. Using the Vegan package within R, seven variables and taxonomic composition (family level) data in 18 communities were analyzed. The numerical consequences were evaluated utilizing Monte Carlo permutation’s technique on 999 permutations. The statistical analyses were conducted using SPSS 16.0. Duncan’s technique was utilized to test the significant differences in compositions and diversity among the mixed and monoculture RPs.
Results
Taxonomic Compositions
There were a total of 1,896 OTUs belonging to 170 families and 7 phyla in the soil fungal communities from 18 soil samples (6 monoculture RPs, 6 Mixed 1 RPs, and 6 Mixed 2 RPs), among which there were a total of 1,758 OTUs belonging to 196 families and 23 phyla in the soil bacterial communities. At the phylum level, there was no difference in soil bacterial compositions between the monoculture and mixed RPs (Figure 1(a)). For fungal structure, the most abundant phylum in soils of monoculture, Mixed 1, and Mixed 2 RPs was Ascomycota. The relative abundance of Ascomycota was higher in monoculture soil than in the mixed RPs; however, the relative abundance of Basidiomycota was lower in the monoculture (Figure 1(b)). Mixed 1 and Mixed 2 RPs were fundamentally similar.
Bacterial Diversity
Diversity estimations at the OTU level of the 16S and 18S rRNA gene banks of three communities from Illumina MiSeq sequencing analysis are shown in Figure 2. In terms of soil bacterial community, the number of OTUs in the monoculture RPs was significantly higher than in the Mixed 1 and Mixed 2 RPs, whereas Mixed 2 was higher than Mixed 1. The Chao1 diversity and Shannon diversity indexes showed that monoculture was significantly higher than the Mixed 1 RP. The number of species, Chao diversity and Shannon diversity indexes of fungal communities had no significant differences among the monoculture, Mixed 1, and Mixed 2 RPs.
PCoA of the bacterial and fungal communities at the OTU level of the three RPs showed some spatial heterogeneity of bacterial and fungal taxonomic compositions between the soils of the monoculture and mixed RPs (Figure 3). For bacteria, the analysis of similarity (ANOSIM) results showed that there was no significant difference in taxonomic composition (p = .0569, at the OTU level) among the monoculture, Mixed 1, and Mixed 2 RPs. However, the ANOSIM results showed that monoculture, Mixed 1, and Mixed 2 RPs were different from each other in fungal taxonomic composition at the OTU level (p < .001).
Affecting Factors
Compared to Mixed 1 and Mixed 2, the WC was higher in the monoculture. Furthermore, the soil pH in the monoculture was significantly higher than in the Mixed 1 and Mixed 2 (Table 1). In addition, SOM, TN, and TP did not differ between the monoculture and mixed plantations. The highest TK was measured in Mixed 1 RPs, followed by monoculture and Mixed 2 RPs. Seven variables and the abundance data of the bacterial and fungal dominant family were used for redundancy analysis. For bacteria, this variable combination explains 44.95% of the total variance of family abundance (Figure 4(a)). TN, TK, pH, and WC described 22.47%, 18.94%, 15.23%, and 15.18% of the total variance, respectively. Additional factors including Rubber Forest, TP, and SOM on the bacterial community composition had relative little effect, accounting for 6.53%, 9.74%, and 11.19% of the variance, respectively. For fungi, this variables combination explained 27.77% of the total variance of the family abundance (Figure 4(b)). The pH and WC explained 15.14% and 8.96% of the total variance, respectively. Additional factors, that is, TN, Rubber Forest, TP, TK, and SOM for fungal community composition had little effect, describing only 5.50%, 5.49%, 3.61%, 3.15%, and 2.96% of the variance, respectively.
Discussion
Our results showed that the bacterial composition of the soil from the three RPs was basically similar. However, there was a significant difference in fungal communities among the three RPs at both the phylum and OTU level. The mixed RPs, that is, Mixed 1 and Mixed 2 RPs, were very similar in terms of phylum composition (ANOSIM results, p = .65) and OTU composition (ANOSIM results, p = .07). This indicated that intercropping other tree species into a RP has only a limited effect on soil bacterial composition. A previous study conducted in Hainan and Xishuangbanna showed that rubber forests and other tropical forests have different soil bacterial (Lan et al., 2017a; Lan, Li, Wu, & Xie, 2017b) and fungal compositions (Lan et al., 2017b). A new viewpoint has been drawn from our research, that is, monoculture and mixed RPs are similar in soil bacterial composition but different in fungal composition. Previous studies also showed mature rubber-based agroforestry systems can alleviate the adverse effects of rubber monoculture on bacterial biodiversity (Liu et al., 2019). Therefore, the effect of mixed RPs on the soil microbial community still needs long-term investigation.
In our investigation, the species richness, Chao, and Shannon diversity of the bacterial communities of the monoculture were higher than the Mixed 1 and Mixed 2 RPs, whereas all diversity indexes of the fungal communities were comparable between the monoculture and mixed RPs. These results showed that mixed RPs did not result in a higher diversity of bacteria. A study conducted in the Amazon showed that the local taxonomic and phylogenetic diversity of the soil bacteria was higher in agroecosystems compared to natural ecosystems (Rodrigues et al., 2013). A previous study observed that soil bacteria in a bamboo forest and a fir forest showed a prominent diversity, but with a lower total sequencing number than from native forest and mixed forest. A high variety of bacterial groups in monoculture rubber forests of agricultural practices would result in a considerable variety of soil microscopic organisms (Kerfahi et al., 2016; Lan et al., 2017, 2017a, 2017b). Previous studies have shown that the percentage of Gram-negative bacterial phospholipid fatty acids in jungle RPs was lower than in monoculture RPs (Krashevska et al., 2015). Our results further confirmed that monoculture RPs have a higher bacterial diversity than mixed forests. In addition, our results also revealed that there were no significant differences in fungal diversity of soils between monoculture and mixed RPs. Our findings are similar to previous study conducted in the Jambi Province of southwest Sumatra, Indonesia (Krashevska et al., 2015), which showed that the total fungal phospholipid fatty acids of jungle RPs and monoculture RPs were comparable.
Our results showed that soil nutrition (such as TN and TK) and soil pH are the main drivers of the bacterial composition. However, soil pH and WC are the main drivers of the fungal composition. Our results are in line with a previous study in Danzhou, Hainan Island, which revealed that soil nutrients explained 43.05% of the total variance of soil bacterial composition in RPs (Lan et al., 2017a). In fact, numerous studies have shown that soil chemical properties, including soil pH, nitrogen, and carbon proportions together correspond with the bacterial community’s composition and diversity (Lee, Barbier, Bottos, Mcdonald, & Cary, 2012; Smith, Barrett, Tusnady, Rejto, & Cary, 2010; Van Horn et al., 2013; Zeglin et al., 2011). Among the soil properties, pH is an imperative variable that impacts both soil bacterial communities and their alpha diversity (Fierer & Jackson, 2006). In our study, soil pH in the monoculture was relatively higher than in the mixed RP, which resulted in higher diversity among bacterial community in the monoculture than mixed RP. Very similar results were found in Hainan (Lan et al., 2018; Lan et al., 2017a) and Xishuangbanna (Lan et al., 2017b). Increased pH presumably is a key factor favoring nutrient capture in rubber production systems (Krashevska, Klarner, Widyastuti, Maraun, & Scheu, 2015). A previous study also confirmed that soil pH was the main driver of soil microbial composition in RPs (Krashevska et al., 2015).
Implications for Conservation
Our study revealed that there was a significant difference in fungal composition between monoculture and mixed plantations. The diversity of soil bacterial communities was higher in monoculture than in mixed RPs due to the higher soil pH in the monoculture. Our results also confirmed that soil pH is the main driver for both bacterial and fungal soil structure. From the point of view of soil microbial diversity, it seems that monoculture RPs are better than mixed RPs. However, from the perspective of plant diversity above ground, the mixed RP is better than the monoculture RP. Here, we suggest that more alkaline fertilizers should be applied in mixed RPs to improve the soil pH and soil microbial diversity. In this way, the total diversity would be enhanced to some extent. At the same time, this treatment can also increase the rubber latex yield.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research work was financially sponsored by the National Natural Science Foundation of China (Grant No. 31770661), the Finance Science and Technology Project of Hainan Province (ZDYF2019145), the Lancang-Mekong River International Cooperation Project (081720203994192003), and the Earmarked Fund for China Agriculture Research System (CARS-34- ZP3).