2.7 | Statistical analysis
Paired-end reads of the 16S rRNA gene were assembled using Flash (v1.2.11) (https://ccb.jhu.edu/software/FLASH/index.shtml) to obtain the raw tags. The relative abundance taxonomic summaries, beta diversity, and rarefactions were examined with QIIME (v1.9.1) (http://qiime.org/install/index.html). The OTU clustering of sequences was performed using UPARSE (v7.0.1090) (http://www.drive5.com/uparse/). Taxonomic classification was conducted using the ribosomal database project classifier (v2.11) (https://sourceforge.net/projects/rdp-classifier/). Concatenated sequences were detected using USEARCH (v7.0) (http://www.drive5.com/usearch/). Alpha and beta diversity estimates were calculated using MOTHUR (v1.30.2) (https://www.mothur.org/wiki/Download_mothur). Linear discriminant analysis (LDA) effect size (LEfSe) was performed to evaluate the different taxa among the soil depth, slope position, and forest types (http://huttenhower.sph.harvard.edu/galaxy/root?tool_id=lefse_upload). SEM was used to assess the effects of the G. jasminoidesintercropping.