Genes Genet. Syst. (2017) 92, p. 1–20 Soil
bacterial community structure in the tropics
1
Geographical variation in soil bacterial community structure in tropical forests in Southeast Asia and temperate forests in Japan based on pyrosequencing analysis of 16S rRNA Natsumi Ito1, Hiroko Iwanaga1, Suliana Charles2, Bibian Diway2, John Sabang2, Lucy Chong2, Satoshi Nanami3, Koichi Kamiya4, Shawn Lum5, Ulfah J. Siregar6, Ko Harada4 and Naohiko T. Miyashita1* 1
Laboratory of Plant Genetics, Graduate School of Agriculture, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyoku, Kyoto, Kyoto Prefecture 606-8502, Japan 2 Sarawak Forestry Corporation, Lot 218, KCLD, Jalan Tapang, Kota Sentosa, 93250 Kuching, Sarawak, Malaysia 3 Laboratory of Plant Ecology, Graduate School of Science, Osaka City University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka, Osaka Prefecture 558-8585, Japan 4 Department of Forest Resources, Faculty of Agriculture, Ehime University, Tarumi, Matsuyama, Ehime Prefecture 790-8566, Japan 5 Department of Natural Sciences and Science Education, National Institute of Education, Nanyang Technological University, Nanyang Walk, Singapore 637616 6 Faculty of Forestry, Bogor Agricultural University and Southeast Asian Regional Centre for Tropical Biology, Jl. Raya Tajur Km. 6, Bogor 16134, Indonesia (Received 10 March 2016, accepted 19 October 2016; J-STAGE Advance published date: 21 December 2016)
Geographical variation in soil bacterial community structure in 26 tropical forests in Southeast Asia (Malaysia, Indonesia and Singapore) and two temperate forests in Japan was investigated to elucidate the environmental factors and mechanisms that influence biogeography of soil bacterial diversity and composition. Despite substantial environmental differences, bacterial phyla were represented in similar proportions, with Acidobacteria and Proteobacteria the dominant phyla in all forests except one mangrove forest in Sarawak, although highly significant heterogeneity in frequency of individual phyla was detected among forests. In contrast, species diversity (α-diversity) differed to a much greater extent, being nearly six-fold higher in the mangrove forest (Chao1 index = 6,862) than in forests in Singapore and Sarawak (~1,250). In addition, natural mixed dipterocarp forests had lower species diversity than acacia and oil palm plantations, indicating that aboveground tree composition does not influence soil bacterial diversity. Shannon and Chao1 indices were correlated positively, implying that skewed operational taxonomic unit (OTU) distribution was associated with the abundance of overall and rare (singleton) OTUs. No OTUs were represented in all 28 forests, and forest-specific OTUs accounted for over 70% of all detected OTUs. Forests that were geographically adjacent and/or of the same forest type had similar bacterial species composition, and a positive correlation was detected between species divergence (β-diversity) and direct distance between forests. Both α- and β-diversities were correlated with soil pH. These results suggest that soil bacterial communities in different forests evolve largely independently of each other and that soil bacterial communities adapt to their local environment, modulated by bacterial dispersal (distance effect) and forest type. Therefore, we conclude that the biogeography of soil bacteria communities described here is non-random, reflecting the influences of contemporary environmental factors and evolutionary history. Key words: biogeography, soil bacteria, tropical forest, Southeast Asia, 16S rRNA INTRODUCTION Edited by K. Ryo Takahashi * Corresponding author. E-mail:
[email protected] DOI: http://doi.org/10.1266/ggs.16-00013
Soil microbes are among the most diverse and abundant organisms on Earth (Whitman et al., 1998; Torsvik et al., 2002; Venter et al., 2004), and they play important roles
2
N. ITO et al.
in ecological systems, promoting decomposition of organic matter (i.e., plant litter and dead animals), supplying nutrients and stimulating growth and health of aboveground plant life, and facilitating the carbon and nitrogen cycles (Chapin, 1980; Högberg et al., 2001; Kowalchuk and Stephen, 2001; Hättenschwiler et al., 2005; van der Heijden et al., 2008). Soil microbes can exist in or on the surface of plant roots in a mutualistic or commensal manner, or they can remain physically apart from but near plant roots; in the latter arrangement, the bacteria modulate the supply or availability of plant nutrients in the soil, and the fractional rate of resources (van der Heijden et al., 2008). On the other hand, plant litter is also a growth substrate and a source of nutrients for soil microbes (Wu et al., 2011), and substances secreted from the roots of living plants also support the growth and survival of microbes (Bais et al., 2006). This suggests a mutual reciprocal relationship between the growth and productivity of the aboveground vegetation and those of soil microbes (Wardle et al., 2004; Harris, 2009). Molecular studies have recently examined the relationship between the above- and belowground biotas in forest environments. A 16S rRNA pyrosequencing analysis showed that the composition of the soil bacterial communities in natural forests in Australia varied with variation in the dominant aboveground vegetation (Hackl et al., 2004); for example, high-GC content Grampositive bacteria (Firmicutes and Actinomycetes) were most dominant (49% of all data) in pine forests, followed by Alphaproteobacteria. In contrast, in oak-hornbeam and spruce-fir-beech forests, Holophagae of Acidobacteria was most dominant, followed by Verrucomicrobia and Alphaproteobacteria, respectively. However, these forests also differed in soil chemical properties, microbe biomass and metabolic rate of nutrients, making it unclear which factor has the greatest influence on soil bacteria composition. A study using automated ribosomal internal transcribed spacer (ITS) analysis (ARISA) reported that the structure of bacterial and fungal communities in fynbos soils in South Africa is stable over time, but varies with the nature of the aboveground vegetation at different sampling points, suggesting that the growth of some bacteria and fungi may depend on the presence of particular plant species (Slabbert et al., 2010). In a study of forest soils in Oregon, USA, in which terminal restriction fragment length polymorphism (T-RFLP) analysis for eubacterial and archaeal 16S rRNA and length heterogeneity (LH)-PCR analysis for fungal ITS sequences were employed, it was reported that the fungal and archaeal community structures were different in forests dominated aboveground by Douglas-fir or red alder vegetation, but that the structure of the bacterial community did not vary (Yarwood et al., 2010). As noted above, these studies also suggest that the presence of certain plant species aboveground influences the growth and survival of spe-
cific microbes belowground. Nevertheless, many complex interacting factors, both biotic and abiotic, influence the distribution of biota in a single soil ecosystem, such that the exact factors that influence the structure of the soil microbial community are not completely understood (Wardle et al., 2004). As an example of biotic factors, both the developmental stage and the species of crops (alfalfa, clover and beans) were reported to influence the structure of the bacterial community in soil and in association with plant roots in cultivated land in Germany, based on the temperature gradient gel electrophoresis (TGGE) separation of 16S rRNA (Wieland et al., 2001). A change in land use management from forest to pasture also correlated with change in the structure of the soil bacterial community in Hawaii and Brazil, using GC content abundance analysis of SSU rRNA and T-RFLP of 16S rRNA, respectively (Nüsslein and Tiedje, 1999; Jesus et al., 2009). Additional physicochemical characteristics of soil, such as temperature (Fiere et al., 2005), salinity (Lozupone and Knight, 2007), pH (Fierer and Jackson, 2006; Tripathi et al., 2012), nutrient abundance (Fiere et al., 2012) and redox potential (DeAngelis et al., 2010), also affect growth, survival and function of soil microbes. A study of 98 forest soils in North and South America revealed a strong association between soil pH and the microbial community structure, based on T-RFLP and pyrosequencing analyses of 16S rRNA (Fierer and Jackson, 2006; Lauber et al., 2009). Even when two or more environmental conditions (e.g., temperature, moisture content, C/N ratio, texture, pH) were considered simultaneously, soil pH appeared to be most strongly correlated with bacterial diversity. Considering the relationship between above- and belowground biotas, in areas with high biodiversity for both animal and plant species, such as tropical forests, the soil microbial ecosystem is expected to be more complex than in areas with less biodiversity (Myers et al., 2000). One study on pasture and forest soils in tropical Brazil, using a culture-independent metagenomic approach based on 16S rRNA sequences, identified previously unreported sequences, and suggested that “immense” microbial diversity existed in the studied regions (Borneman and Triplett, 1997). On the other hand, subtropical soils in Brazil and Florida and temperate soils in Illinois and Canada were comparable with respect to the level of soil bacterial diversity (104 species per sample), based on pyrosequencing analyses of the 16S rRNA region (Roesch et al., 2007). Recently, by studying 25 temperate grassland sites from four continents (Africa, Australia, Europe and North America) by high-throughput sequencing of 16S rRNA (bacteria and archaea) and an ITS region (fungi), it was shown that plant α-diversity was poorly related to that for microbial groups (Prober et al., 2015). In our previous work on five tropical forests in Malaysia and one temperate forest in Japan, using pyrosequencing analy-
Soil bacterial community structure in the tropics ses of bacterial 16S rRNA (Miyashita et al., 2013), we also observed a similar composition of higher taxonomic taxa and similar bacterial species diversity in the six forests. The present study investigates bacterial diversity and composition in 28 forest soils based on pyrosequencing analyses of the 16S rRNA region. This extends our previous study (Miyashita et al., 2013) by including an additional 21 tropical forests in Malaysia, Indonesia and Singapore, and one temperate forest in Japan. These 28 forests had different aboveground vegetation, including dipterocarp forest, acacia plantation, oil palm plantation, limestone forest and mangrove. Thus, it should be possible to determine whether the conclusion of our previous study would hold for forest soil bacterial communities having more diverse geographic origins and aboveground vegetation. The aim of this study was to elucidate the environmental factors and mechanisms that influence the spatial biogeography of a soil bacterial community. For this purpose, we conducted the following analyses: 1) comparison of higher taxonomic composition and species diversity among soil bacterial communities, geographic regions and aboveground vegetation (forest type); 2) comparison of operational taxonomic unit (OTU) occurrence within and between forests; 3) examination of the distance effect by testing the correlation between β-diversity and direct distance between forests; 4) examination of the effect of geographical region and aboveground vegetation (forest type) on the soil bacterial community (composition); and 5) evaluation of the association between α- and β-diversities and physicochemical properties of soil. Primary conclusions were as follows: 1) among the investigated forests, bacterial phyla were represented in similar proportions, despite diverse environmental differences, but species diversity differed considerably; 2) OTU distribution was skewed, and forest-specific OTUs were predominant; 3) a significant and positive correlation was detected between β-diversity and direct distance between forests, suggesting local adaptation and bacterial dispersal among forests; 4) both geographical region and aboveground vegetation (forest type) influenced soil bacterial composition; 5) α- and β-diversities were correlated with soil pH. MATERIALS AND METHODS Soil sampling and DNA extraction Soil samples were collected from 26 tropical forests in Malaysia (Sarawak), Indonesia (Java, Sumatra) and Singapore and from two temperate forests in Japan in 2010–2011 (Supplementary Fig. S1). A detailed description of the environmental condition and forest type at each sampling location is given in Table 1. In each forest, ~10 g of soil located ~10 cm below the surface was sampled at 10 sampling sites 1–10 m apart from each other. The soil samples were mixed and transferred into a tube containing
3
extraction beads and buffer immediately after removal, and then frozen. Frozen soil samples were shipped to laboratories. Microbial DNA was extracted using a PowerSoil DNA Isolation Kit or a PowerMax Soil DNA Isolation Kit (MO BIO Laboratories, Carlsbad, CA, USA) according to the manufacturer’s instructions. Extracted DNA was purified using a PowerClean DNA Clean-Up Kit (MO BIO Laboratories). For Indonesian forests, 1 kg of soil was collected at two randomly chosen sampling sites to determine soil physicochemical properties (soil type, moisture content, organic C content, N content, and pH). The average value of the two sampling sites for the Indonesian samples was used for analyses (Supplementary Table S1). For the other forests, 1 kg of soil pooled from 10 sampling sites was used to determine physicochemical properties. However, some samples were not examined for one or more soil properties (Supplementary Table S1). All necessary permits were obtained for the described field studies. Research permits for the studied forests and export permits for the extracted DNA samples were issued by the Forest Department Sarawak, the National Parks Board (Singapore Botanical Gardens), the Ministry of Environment and Forestry and Indonesian Academy of Sciences through SEAMEO BIOTROP (Bogor, Indonesia), the Ashiu Experimental Forest (the Field Science Education Center, Kyoto University), and the Ibuki Property District Manager (Maibara, Shiga Prefecture). PCR amplification of bacterial 16S rRNA region Extracted microbial DNA was used as template for PCR amplification of part of the 16S rRNA region as described (Miyashita et al., 2013; Miyashita, 2015). Briefly, two PCR primers were used for reactions in this study: 9F 5’-GAGTTTGATCCTGGCTCAG-3’ and 533R 5’-TIACCGIIICTICTGGCAC-3’. These primers amplify a 510-bp fragment that includes three variable regions (V1, V2 and V3) of the 16S rRNA gene in bacteria, archaea and eukarya (Baker et al., 2003). To distinguish the 28 samples from each other, a unique tag sequence was attached to the 5’ end of the forward and reverse primers, as described (Miyashita et al., 2013). The PCR solution (50 μl) contained 0.2 mM dNTPs, 2.5 units of TaKaRa ExTaq (Takara Bio, Otsu, Japan), 10 × Taq buffer, 2 mM MgCl2, 0.2 μM each of the forward and reverse primers, and ~10 ng of template DNA. The PCR conditions used were, sequentially: 94 °C for 2 min; 30 cycles of 94 °C for 30 sec denaturation, 55 °C for 30 sec annealing and 72 °C for 1 min extension; and 72 °C for 5 min extension. For each sample, three replicate PCR products were combined and quantified by agarose gel electrophoresis, and ~200 ng of PCR amplicons from each of the 28 samples was prepared for pyrosequencing. The PCR amplicons were cleaned using the QIAquick PCR Purification Kit (Qiagen, Venlo, Netherlands) following the manufacturer’s instructions.
4
N. ITO et al. Table 1. List of sampling locations ID
Location
Country
Latitude/ Longitude
Altitude, m
Forest type
Date of soil collection
JASoil
Jasinga
Java, Indonesia
6.54S 106.53E
335
Oil palm plantation
2010.7.31
MARaca5
Maribaya, Parung Panjang
Java, Indonesia
6.39S 106.47E
105
Acacia mangium plantation (5 years)
2010.7.31
MARaca9
Maribaya, Parung Panjang
Java, Indonesia
6.39S 106.46E
97
Acacia mangium plantation (9 years)
2010.7.31
BTNmdf
Bukit Tigapuluh National Park Sumatra, Indonesia 0.84S 102.52E
289
Mixed dipterocarp forest
2010.8.2
BTNbur
Bukit Tigapuluh National Park Sumatra, Indonesia 0.79S 102.54E
115
Burned forest (1 month after)
2010.8.3
TIMrid
Bukit Timah National Park
Singapore
1.35N 103.78E
128
Mixed dipterocarp forest (ridge area)
2011.8.23
TIMval
Bukit Timah National Park
Singapore
1.35N 103.78E
100
Mixed dipterocarp forest (valley area)
2011.8.23
TIM2nd
Bukit Timah National Park
Singapore
1.35N 103.78E
89
Secondary forest (60 years)
2011.8.23
BAUlim
Bau, Fairy Cave
Sarawak, Malaysia
1.38N 110.12E
47
Limestone forest
2010.10.13
BAU2nd
Bau, Fairy Cave
Sarawak, Malaysia
1.38N 110.12E
37
Secondary forest
2010.10.13
SEMpin
Botanical Research Center, Semengoh
Sarawak, Malaysia
1.41N 110.31E
41
Dipterocarp plantation (Shorea pinanga)
2010.10.15
SEMspl
Botanical Research Center, Semengoh
Sarawak, Malaysia
1.41N 110.31E
47
Dipterocarp plantation (Shorea splendida)
2010.10.15
SEMmac
Botanical Research Center, Semengoh
Sarawak, Malaysia
1.42N 110.31E
36
Dipterocarp plantation (Shorea macrophylla)
2010.10.15
OPAoil
Opar
Sarawak, Malaysia
1.44N 110.10E
49
Oil palm plantation
2010.10.15
LUNaca
Lundu
Sarawak, Malaysia
1.56N 109.89E
23
Acacia mangium plantation (20 years)
2010.10.13
BKOmdf
Bako National Park
Sarawak, Malaysia
Missing
Mixed dipterocarp forest
2010.10.20
BKOpts
Bako National Park
Sarawak, Malaysia
1.71N 110.44E
15
Peat swamp forest
2010.10.20
BKOker
Bako National Park
Sarawak, Malaysia
1.72N 110.45E
157
Kerangas forest
2010.10.20
BKOman
Bako National Park
Sarawak, Malaysia
1.72N 110.45E
17
Mangrove
2010.10.20
NIAmdf
Niah National Park
Sarawak, Malaysia
3.82N 113.77E
35
Mixed dipterocarp forest
2010.10.26
NIAlim
Niah National Park
Sarawak, Malaysia
3.82N 113.77E
20
Limestone forest
2010.10.26
LHNrid
Lambir Hills National Park
Sarawak, Malaysia
4.19N 114.02E
162
Mixed dipterocarp forest (ridge area)
2011.4.18
LHNval
Lambir Hills National Park
Sarawak, Malaysia
4.18N 114.02E
76
Mixed dipterocarp forest (valley area)
2011.4.18
BER2nd
Bakam Experimental Reserve
Sarawak, Malaysia
4.26N 113.99E
76
Secondary forest
2011.4.19
BERbur
Bakam Experimental Reserve
Sarawak, Malaysia
4.26N 113.99E
84
Burned forest (20 months after)
2011.4.19
BERrem
Bakam Experimental Reserve
Sarawak, Malaysia
4.26N 113.99E
81
Mixed dipterocarp forest (remnant forest)
2011.4.19
Ashiu
Ashiu Forest Research Station
Kyoto, Japan
35.35N 135.73E
757
Temperate broadleaf and mixed forest
2011.8.2
IBUlim
Mt. Ibuki
Shiga, Japan
35.40N 136.39E
633
Limestone forest
2011.7.25
Missing
Soil bacterial community structure in the tropics Pyrosequencing and data analyses The mass of pooled PCR amplicons of the 16S rRNA region was single-end-sequenced using a 454 Life Science Genome Sequencer FLX (Roche, Basel, Switzerland) by Hokkaido System Science (Sapporo, Japan). After trimming 454 adapter sequences, and tag and primer sequences, PCR and pyrosequencing errors were corrected by a de-noising program in Quantitative Insight Into Microbial Ecology (QIIME, http://qiime.org) (Quince et al., 2009; Caporaso et al., 2010a). In this process, reads shorter than 300 bp and/or having a quality score lower than 50 were eliminated from the data set. A total of 267,073 sequences were used for further analyses. All of the sequences were clustered into OTUs based on their sequence identity using the uclust method (Caporaso et al., 2010b), and a cluster at 3% distance was defined as a “species” (OTU) following the conventional definitions (Bond et al., 1995; Schloss and Handelsman, 2005; Huber et al., 2007). For each OTU a representative sequence that was the most abundant sequence in an OTU cluster was selected, and a set of the representative sequences were aligned using the PyNAST method (Caporaso et al., 2010b) and taxonomically assigned based on the 16S rRNA database in the Ribosomal Database Project (RDP) at Michigan State University (Cole et al., 2009) using the uclust consensus taxonomy classifier. Chimeric sequences were identified using ChimeraSlayer (http://microbiomeutil.sourceforge. net/#A_CS), and OTUs represented by these chimeric sequences were excluded from all downstream analyses. The above-mentioned cluster analyses, alignment, taxonomic assignment and chimera identification were performed using QIIME tools. A total of 256,693 nonchimerical sequences were used for further analyses. To compare the composition of bacterial taxa among the forests, Spearman’s rank correlation was calculated for all the pairwise comparisons between forests for 29 phyla detected and for the top 30 ranked taxa of class and order, using Excel and a web tool (http://www.gen-info.osaka-u. ac.jp/testdocs/tomocom/spea.html). Heterogeneity in the taxon frequency among the forests was tested by chisquare tests, using marginal frequencies of taxon and location to obtain the expectation with Excel. Based on the results of cluster analyses, soil bacterial species diversity (α-diversity) in each forest was estimated; α-diversity was represented by two indices, the Chao1 index (Chao, 1984), which estimates the minimum number of OTUs based on the observed number of OTUs, and the Shannon H’ index (and equitability: Shannon index divided by loge (the number of OTUs)) (Shannon, 1948), which estimates the evenness of OTU distribution. Rarefaction curves (Sanders, 1968; Gotelli and Colwell, 2001) were created by a tool in QIIME. To normalize read number variation, 2,951 sequences (the minimum number of reads in all forest samples) were randomly resampled 100 times from each sample. For the resampled sequence sets, the
5
average values for the 100 resampled sequence sets were estimated for the Chao1 index and Shannon H’ index (and equitability). Pearson’s correlation coefficients were used to examine the relationship between α-diversity and environmental factors. To estimate β-diversity (differentiation) between samples, weighted UniFrac (Lozupone et al., 2007) and the 1 - abundant Jaccard index (Chao et al., 2005), which measures similarity in OTU composition between forests, were calculated for all pairwise comparisons. Based on the β-diversity distance matrices, Principal Coordinates Analysis (PCoA) plots (Anderson and Willis, 2003) and Unweighted Pair Group Method with Arithmetic Mean (UPGMA) dendrograms were constructed. Calculation of β-diversity indices and construction of dendrograms and PCoA plots were performed using QIIME tools. Pearson’s correlation coefficients were used to examine the relationship between coordinate values of PCoA and environmental factors. The latitude and longitude data of sampling locations (Table 1) were used to estimate the direct distance between all the pairs of sampling locations with a web tool (http://vldb.gsi.go.jp/ sokuchi/surveycalc/surveycalc/bl2stf.html) of the Geospatial Information Authority of Japan, Ministry of Land, Infrastructure, Transportation and Tourism. Since latitude and longitude data are missing for the mixed dipterocarp forest in Bako National Park (BKOmdf), the total number of pairs was 351 ( = 27 × 26 / 2). Statistical significance of the correlation coefficient between the pairwise direct distance and β-diversity (1 - abundant Jaccard index and weighted UniFrac) was examined by the Mantel test (Sokal and Rohlf, 2000). Mantel correlogram analysis (Sokal, 1986; Legendre and Legendre, 1998) was conducted to examine autocorrelation in the soil bacterial community structure. The number of distance classes was determined to be nine based on Sturge’s rule (the number of classes = 1 + log2(m), where m = the number of distances in the β-diversity triangular matrix; in this study, m = 351). Two different distance classes were used (Legendre and Legendre, 1998): equal distance and equal count (Supplementary Table S2). The Mantel test and Mantel correlogram analysis were performed using PASSaGE v2 (Pattern Analysis, Spatial Statistics and Geographic Exegesis) software (http://www.passagesoftware.net/, Rosenberg and Anderson, 2011). Statistical significance of the correlation coefficient for the Mantel test and Mantel statistics for the correlogram analysis was tested by 10,000 permutations. For the Mantel correlogram analysis, the progressive Bonferroni correction was used. Permutational ANOVA (PERMANOVA, Anderson, 2001) was applied to the estimated triangular β-diversity matrices to test the effect of soil sampling geographical region and aboveground forest type on β-diversity in terms of community structure (“location” and “dispersion”) within and between groups (region and forest type). One-way design for each
6
N. ITO et al.
of the groups, in which factors were fixed, was set up. To test the statistical significance of the pseudo-F value, the distribution of pseudo-F statistics obtained by permutation of raw data for 9,999 times by the unrestricted method was used. PERMANOVA was performed using the PRIMER v7 and PERMANOVA + add-on computer programs (http://www.primer-e.com/). The 16S rRNA sequences have been deposited under DDBJ accession nos. DRR014854 to DRR014881. RESULTS Bacterial composition at higher taxonomic levels in investigated forest soils A total of 29 of 55 phyla listed in the RDP database were detected, representing 89.2% of all sequence reads (Fig. 1 and Supplementary Fig. S2). The average number of phyla per forest was 16.4. A mangrove forest in the Bako National Park (BKOman) had the highest number of phyla (26) and a kerangas forest in the Bako National Park (BKOker) had the lowest number of phyla (13). The two temperate forests in Japan (Ashiu and IBUlim) had a moderate number of phyla (16 and 18, respectively). Nine phyla (Acidobacteria, Proteobacteria, Verrucomicrobia, Planctomycetes, Bacteroidetes, TM7, Actinobacteria, Firmicutes and OP10) were detected in all 28 forests, and three phyla (Bacteria_incertae_sedis, Gemmatimonadetes and
OP11) were found in 27 forests, the exceptional forests being BKOman, a peat swamp forest in Bako (BKOpts) and BKOpts, respectively. Several phyla were found uniquely in one forest, including Deferribacteres, Fusobacteria, Deinococcus-Thermus, and Thermotogae in the BKOman forest, Synergistetes in Bakam secondary forest (BER2nd), and Chlamydiae in Lambir Hills National Park Valley area (LHNval). A total of 10.8% of the 16S rRNA sequence reads were assigned to unclassified bacteria, a group belonging to the domain Bacteria but not assigned to any existing phylum, ranging from 4.8% in a mixed dipterocarp forest in Bako (BKOmdf) to 26.2% in a secondary forest in Bau (BAU2nd). As a whole, Acidobacteria were the most dominant phylum (41%) among the 29 phyla detected, and Proteobacteria followed (23%) (Fig. 1). The predominant five phyla (Acidobacteria, Proteobacteria, Verrucomicrobia, Planctomycetes and Bacteroidetes) accounted for more than 80% of the phyla detected. However, there was considerable variation in frequencies of bacterial phyla among forests. Exceptional composition of phyla was observed in the BKOman mangrove forest (Fig. 1), where Acidobacteria were represented at a lower frequency (3%) and Proteobacteria, Bacteroidetes, and Cyanobacteria were represented at higher frequencies (37%, 11% and 10%, respectively) than in other forests. To examine the similarity in phylum composition in the forests, Spearman’s
Fig. 1. Composition of phyla in the 28 soil samples.
Soil bacterial community structure in the tropics Table 2. Rank correlation of taxa at higher taxonomic levels (378 pairwise comparisons) Taxonomic level Phylum
Range
Average
Number of significant pairwise comparisons (%) 5%
1%
0.55–0.99
0.83
378 (100.0)
378 (100.0)
Class
−0.02–0.97
0.69
334 (88.4)
305 (80.7)
Order
−0.18–0.96
0.66
353 (93.4)
340 (89.9)
rank correlation was calculated for all pairwise comparisons (378 comparisons) of 28 forests (Table 2). Despite the variation in phylum frequency among the forests and above-mentioned environmental conditions, statistically significant rank correlations were detected at least at the 1% level in all comparisons, indicating similar taxonomic composition at the phylum level. The statistical significance of rank correlations probably reflects the dominance of the top five phyla in all forests studied. However, when heterogeneity in frequency of individual phyla was evaluated by chi-square tests, statistically significant heterogeneity was detected at least at the 1% level for all phyla detected. These results indicated that, despite the overall similarity in phylum composition, the representation of individual phyla in these forests was influenced by some forest-specific factor(s). At the class and order levels, about 87% and 77% of the total sequence reads were assigned to 61 classes and 67 orders, respectively (Supplementary Fig. S2). Compared with the phylum composition, the class and order compositions were more variable among forests (Supplementary Fig. S3). Some dominant taxa were common, and the top five or six dominant taxa accounted for more than 50% of assigned classes (Supplementary Fig. S3A). Acidobacteria subgroup Gp1 was the most dominant class (20%) of all classes detected, and Acidobacteria Gp2 (14%) and Alphaproteobacteria (10%) followed. In most forests, Acidobacteria Gp1 and Gp2 were the major classes of the phylum Acidobacteria, accounting for 80–90% of Acidobacteria (Supplementary Fig. S4A). However, instead of Acidobacteria Gp1 and Gp2, Gp4 and Gp6 were predominant in an oil palm forest in Jasinga (JASoil), a secondary forest in Bau (BAU2nd), and a limestone forest in Mt. Ibuki (IBUlim). Similarly, Acidobacteria Gp4 and Gp10 were predominant in BKOman, and Gp5 and Gp6 were predominant in limestone forests in Bau (BAUlim) and Niah (NIAlim). The class composition within Proteobacteria also differed among forests. In most forests, Alphaproteobacteria was the most dominant of five classes, but in some forests (BAU2nd, BAUlim, BKOman, NIAlim and IBUlim) mentioned above, Betaproteobacteria was the most dominant class (Supplementary Fig. S4B). At the class and order levels, Spearman’s rank correlations were lower than at the phylum level, and were not sig-
7
nificant for some pairs, but the percentage was low (11.6 and 6.6% for class and order, respectively, Table 2), still consistent with the generally similar taxonomic composition in the forests studied at the class and order levels. The non-significant rank correlations involved BAUlim and BKOman forests, ndicating that these two forests differed in class and order composition from the other forests. However, statistically significant heterogeneity was detected at the 1% level for all classes and orders as well as phyla, implying that forest-specific factors influence the frequency of taxa at the level of class and order. Soil bacterial community structure at the species level To examine bacterial composition and diversity at the species level, cluster analysis based on 97% sequence identity was carried out for the entire data set. Rarefaction curves at the 3% divergence level did not plateau for most of the forests, so that the observed number of OTUs was insufficiently represented for actual bacterial diversity (Supplementary Fig. S5). The total number of OTUs obtained from the cluster analysis was 21,223 for the 28 forests (Table 3). The average number of OTUs per forest was 1,414, ranging from 855 (BKOmdf, a mixed dipterocarp forest in Bako) to 2,992 (BKOman). The estimated Chao1 index for all 28 forests was 45,130.4 (Table 3). The average Chao1 index per forest was 2,569, ranging from 1,245.7 (SEMmac, a Shorea macrophylla plantation in Semengoh) to 6,861.8 (BKOman). These results imply a large variation in bacterial species diversity among investigated forests. The two Japanese forests had intermediate levels of species diversity, indicating that temperate forests did not necessarily have lower bacterial species diversity than tropical forests. All the Shannon indices estimated were high, indicating highly skewed OTU distribution in all the forests (Table 3). In fact, Pearson’s correlation between the Shannon index and Chao1 index or percentage of singleton OTUs was highly significant (r = 0.770, df = 26, P < 0.01 for Chao1 index and r = 0.689, df = 26, P < 0.01 for percentage of singleton OTUs), indicating that the bias of OTU distribution was associated with the abundance of overall and rare OTUs in the forest soils. The BKOman forest had the highest Shannon index and percentage of singleton OTUs, while forests with low Shannon indices (for example, TIMrid, a ridge area of mixed dipterocarp forest in the Bukit Timah National Park, BER2nd and SEMmac) had lower percentages of singleton OTUs. Although Pearson’s correlation between the equitability and percentage of singleton OTUs was also highly significant (r = 0.527, df = 26, P < 0.01), relatively high equitability values indicated the absence of highly predominant OTUs, consistent with the low maximum number of reads per OTU in all forests, except for a secondary forest in the Bakam Experimental Reserve (BER2nd, Table 3). Although Pearson’s correlation coefficients between
8
N. ITO et al. Table 3. Summary of bacterial diversity and abundance of the total data
ID
Number of reads
Number of OTUs
Chao1 index
Shannon Equitability index (H’) (EH’)
Number of singleton OTUs (%)
Number of forest-specific OTUs (%)
Average number of reads per OTU (%)
Max. number of reads per OTU (%)
JASoil
4680
1555
3172.8
6.72
0.914
872 (56.1)
758 (48.7)
3.0 (0.06)
71 (1.52)
MARaca5
9025
1607
2404.5
6.52
0.883
639 (39.8)
417 (25.9)
5.6 (0.06)
212 (2.35)
MARaca9
6589
1408
2204.5
6.43
0.886
600 (42.6)
329 (23.4)
4.7 (0.07)
208 (3.16)
BTNmdf
9883
1326
2272.7
5.86
0.815
564 (42.5)
296 (22.3)
7.5 (0.08)
559 (5.66)
BTNbur
10862
1751
2627.0
6.52
0.873
671 (38.3)
479 (27.4)
6.2 (0.06)
401 (3.69)
TIMrid
8153
868
1260.5
5.30
0.783
317 (36.5)
122 (14.1)
9.4 (0.12)
711 (8.72)
TIMval
15444
1361
2063.9
5.68
0.788
496 (36.4)
335 (24.6)
11.3 (0.07)
806 (5.22)
TIM2nd
17933
1375
1963.1
5.73
0.793
455 (33.1)
362 (26.3)
13.0 (0.07)
807 (4.50)
BAUlim
11086
2564
6051.8
6.80
0.866
1461 (57.0)
1324 (51.6)
4.3 (0.04)
355 (3.20)
BAU2nd
5094
1769
3702.4
6.90
0.923
993 (56.1)
834 (47.1)
2.9 (0.06)
57 (1.12)
SEMpin
13261
1248
1578.5
5.91
0.829
343 (27.5)
311 (24.9)
10.6 (0.08)
441 (3.33)
SEMspl
9118
1134
1478.5
5.96
0.847
357 (31.5)
205 (18.1)
8.0 (0.09)
345 (3.78)
SEMmac
7145
917
1245.7
5.72
0.838
299 (32.6)
111 (12.1)
7.8 (0.11)
257 (3.60)
OPAoil
7997
1581
2710.0
6.41
0.870
708 (44.8)
485 (30.7)
5.1 (0.06)
271 (3.39)
LUNaca
4784
896
1363.5
5.80
0.853
372 (41.5)
143 (16.0)
5.3 (0.11)
165 (3.45)
BKOmdf
5172
855
1303.2
5.76
0.854
344 (40.2)
173 (20.2)
6.0 (0.12)
129 (2.49)
BKOpts
2951
909
1616.8
5.99
0.879
492 (54.1)
376 (41.4)
3.2 (0.11)
115 (3.90)
BKOker
8985
1101
1989.6
5.74
0.819
456 (41.4)
313 (28.4)
8.2 (0.09)
616 (6.86)
BKOman
10045
2992
6861.8
7.10
0.887
1777 (59.4)
2826 (94.5)
3.4 (0.03)
322 (3.21)
NIAmdf
6187
980
1548.4
5.85
0.849
406 (41.4)
192 (19.6)
6.3 (0.10)
289 (4.67)
NIAlim
3859
1269
2830.9
6.40
0.896
754 (59.4)
401 (31.6)
3.0 (0.08)
85 (2.20)
LHNrid
12027
1323
2483.0
5.77
0.803
580 (43.8)
452 (34.2)
9.1 (0.08)
557 (4.63)
LHNval
14852
1543
2988.0
5.85
0.797
680 (44.1)
577 (37.4)
9.6 (0.06)
589 (3.97)
BER2nd
10748
1356
2056.7
5.46
0.758
536 (39.5)
355 (26.2)
7.9 (0.07)
1485 (13.82)
BERbur
5031
1108
1919.1
6.18
0.881
490 (44.2)
255 (23.0)
4.5 (0.09)
203 (4.03)
BERrem
15804
1543
2608.3
5.69
0.776
621 (40.2)
494 (32.0)
10.2 (0.06)
686 (4.34)
9036
998
1974.5
5.59
0.809
455 (45.6)
424 (42.5)
9.1 (0.10)
317 (3.51)
Ashiu IBUlim Total Average
10942
2266
5652.9
6.75
0.874
1259 (55.6)
1587 (70.0)
4.8 (0.04)
355 (3.24)
256693
21223
45130.4
7.73
0.776
10607 (50.0)
14936 (70.4)
12.1 (0.00)
7934 (3.09)
9168
1414
2569
6.08
0.839
643 (43.8)
533 (32.7)
6.8 (0.08)
408 (4.20)
the number of 16S rRNA sequence reads and species diversity estimates were not statistically significant (r = 0.296, df = 26 for the number of OTUs, and 0.155, df = 26 for Chao1 index), the positive correlation coefficients suggested some influence on the diversity estimates. To eliminate the effect of read number variation, the read number was normalized into the minimum number of reads per sample among the 28 forests (2,951), and then the number of OTUs, Chao1 index and Shannon index were recalculated for each sample (Table 4). There was still a large variation in species diversity among the for-
ests examined. The total number of OTUs was 12,660, with an average of 860 over the 28 forests, ranging from 576 (TIMrid) to 1,400 (BKOman). The Chao1 index for the total sample was 23,912.6, with an average of 1,567.0, and was lowest in TIMrid (908.2) and highest in BKOman (3,346.2). Even after adjusting the read number, corresponding with the positive correlation between Shannon and Chao1 indices detected for raw data, BKOman had the highest Shannon index (6.71), and TIMrid had the lowest (5.19). In fact, Pearson’s correlation coefficient between diversity measures for raw and normalized data
Soil bacterial community structure in the tropics
9
Table 4. Summary of α-diversity normalized to the minimum number of sequence reads (2951 reads) a) Forest JASoil MARaca5
Number of OTUs
Chao1 index
Shannon index (H’)
Equitability b) Area (EH’)
N
1191
2502.0
6.58
0.929
Tropical forest
973
1635.5
6.30
0.916
Indonesia
Number Chao1 Shannon Equitability of OTUs index index (H’) (EH’)
26
861.0
1562.9
5.91
0.875
5
977.6
1755.6
6.22
0.904
MARaca9
971
1642.9
6.26
0.910
Java
3
1045.0
1926.7
6.38
0.918
BTNmdf
768
1304.9
5.69
0.856
Sumatra
2
876.5
1499.0
5.99
0.883
BTNbur
985
1693.4
6.28
0.911
Singapore
3
645.0
1051.7
5.41
0.836
Sarawak
TIMrid
576
908.2
5.19
0.817
TIMval
681
1154.5
5.50
0.843
TIM2nd
678
1092.2
5.54
0.851
except BKOman Japan
BAUlim
1180
2645.7
6.46
0.913
BAU2nd
1290
2697.9
6.72
0.939
c) Forest type
SEMpin
724
1135.6
5.73
0.871
Acacia plantation
18
864.6
1594.6
5.91
0.874
17
833.1
1491.5
5.86
0.872
2
853.5
1620.0
5.96
0.882
3
889.7
1477.7
6.09
0.897
SEMspl
735
1135.9
5.81
0.880
Burned forest
2
928.0
1566.5
6.18
0.904
SEMmac
661
971.3
5.61
0.864
Dipterocarp plantation
3
706.7
1081.0
5.72
0.871
OPAoil
982
1728.3
6.21
0.901
Kerangas
1
689.0
1111.0
5.60
0.857
LUNaca
725
1154.7
5.72
0.868
Limestone
3
1122.3
2435.3
6.42
0.914
BKOmdf
674
1071.7
5.68
0.873
Mangrove
1
1400.0
3346.0
6.71
0.926
BKOpts
909
1611.2
5.99
0.879
Mixed dipterocarp
8
698.3
1168.5
5.56
0.850
BKOker
689
1111.2
5.60
0.857
Oil palm plantation
2
1086.5
2115.0
6.40
0.915
BKOman
1400
3346.2
6.71
0.926
Peat swamp forest
1
909.0
1611.0
5.99
0.879
NIAmdf
709
1170.8
5.73
0.873
Secondary forest
3
907.7
1700.7
5.85
0.858
Temperate forest
1
600.0
1012.0
7.87
0.853
NIAlim
1080
2432.5
6.33
0.906
LHNrid
706
1187.5
5.60
0.853
LHNval
741
1288.1
5.64
0.854
BER2nd
755
1312.1
5.28
0.797
BERbur
871
1440.4
6.07
0.897
BERrem
731
1260.8
5.48
0.831
Ashiu IBUlim Total Average
600
1012.3
5.45
0.853
1107
2228.2
6.46
0.922
12660
23912.6
7.77
0.822
860
1567.0
5.92
0.875
Average values of 100 resampled data are shown in a). Average values of a) are shown in b) and c).
was highly significant (r = 0.783, df = 26, P < 0.01 for the number of OTUs, and r = 0.854, df = 26, P < 0.01 for Chao1 index). This indicates that the overall conclusion with respect to soil bacterial species diversity was not influenced by sequence read variation, and that tropical forests did not necessarily have higher soil bacterial species diversity than temperate forests in Japan (Tables 3 and 4a). Soil bacterial species diversity also varied according to the geographical area of the soil sample. Forest soils in Indonesia, especially Java, had the highest bacterial spe-
cies diversity, while forest soils in Singapore had the lowest bacterial species diversity (Table 4b). As for the forest type, as noted above, the mangrove forest (BKOman) had exceptionally high bacterial species diversity and mixed dipteroparp forests, which may be considered natural forests with minimal human activity, had lower diversity than acacia and oil palm plantations, which are artificial forests with a single tree species (Table 4c), indicating that aboveground tree composition was not strongly related to forest soil bacterial diversity. It should be noted that there should be some caution when
10
N. ITO et al.
inferring these relationships among soil bacterial communities from these data, since almost all samples did not actually plateau in rarefaction curves (Supplementary Fig. S5). Occurrence of OTUs within and between investigated forests Based on the concerted evolution and conserved copy number variation of bacterial 16S rRNA cistrons, the number of 16S rRNA sequence reads per OTU may be regarded as the representative of the number of individual bacteria per species (OTU) in the forest soils investigated, although it is difficult to convert the number of sequence reads to the exact number of individual bacteria (Miyashita et al., 2013; Miyashita, 2015). For the total data, the most abundant OTU, which belonged to Acidobacteria subgroup Gp1, represented approximately 3% (7,934) of all reads (256,693), and the second most abundant OTU, which belonged to Acidobacteria subgroup Gp2, represented approximately 2% of reads (5,384) (Table 3). For individual forests, the average maximum number of reads per OTU was 4.2%; however, BER2nd had an exceptionally predominant OTU representing 13.8% of sequence reads. Besides these dominant OTUs with many sequence reads, the average number of reads per OTU was only 12.1 (0.00047%) for all 21,223 OTUs, ranging from 2.9 (0.06% in BAU2nd) to 13.0 (0.07% in TIM2nd) over the 28 forests. In addition, on average, 43.8% of OTUs in a forest were singletons, ranging from
27.5% (SEMpin, Shorea pinanga plantation in Semengo) to 59.4% (BKOman), and 50.0% of OTUs were singletons for all data. These results indicate that most OTUs (species) exist at very low density in each forest, and only a small number of OTUs reach high frequency. No OTUs were represented in all 28 forests. Five OTUs were represented in 27 forests, and four OTUs were represented in all forests except BKOman. On the other hand, forest-specific OTUs accounted for over 70% (14,936) of all OTUs in the total data (Table 3), of which 71.0% (10,607) were singletons. The forest-specific OTUs accounted for 32.7% of OTUs in each forest on average, and, surprisingly, 94.5% (2,826) of OTUs (2,992) in BKOman were uniquely detected in BKOman. This indicated that bacterial composition at the species (OTU) level was largely differentiated by forest, especially for BKOman. Relative divergence of soil bacterial communities in investigated forests To investigate the differentiation of soil bacterial communities, distance matrices calculated for all pairwise comparisons of 28 forests (378 pairwise comparisons), based on the abundant Jaccard index (Chao et al., 2005) and weighted UniFrac (Lozupone et al., 2007), were used to construct UPGMA trees and PCoA plots. In the UPGMA dendrogram based on the abundant Jaccard index (Fig. 2A), the BKOman forest, which had a distinct composition at the phylum level
Fig. 2. UPGMA dendrogram for the 28 forests constructing using (A) abundant Jaccard index and (B) weighted UniFrac distance.
Soil bacterial community structure in the tropics and an extremely high proportion of forest-specific OTUs, was located at the periphery of the tree. Although their external branches were the longest, two Japanese forests, Ashiu and IBUlim, formed a cluster and were included in a large cluster of all Southeast Asian forests, except BKOman. Forests in Singapore and Indonesia (two different islands, except for JASoil) also formed a distinct cluster. In general, forests adjacent to one another in Sarawak formed their own cluster, with the exception of BKOman, NIAlim and NIAmdf. These results indicate that geographic origin (i.e., physical distance) influenced the similarity of soil bacterial composition, rather than aboveground tree vegetation. However, in addition to the singularity of BKOman, the separation of NIAmdf and NIAlim, and the proximity of NIAlim to BAUlim and IBUlim, all of which had limestone soil, suggest that soil properties are important per se, in addition to geographical origin. In the UPGMA dendrogram based on weighted UniFrac (Fig. 2B), the 28 forests were separated into three large groups. BKOman formed the outer branch, representing a solo cluster. Three limestone forests (IBUlim, BAUlim and NIAlim), the oil palm plantation in Jasinga (JASoil) and the secondary forest in Bau (BAU2nd) formed another cluster separated from the third cluster of 22 forests. The five forests in the second cluster were also in close proximity in the tree based on the abundant Jaccard index (Fig. 2A), although IBUlim and JASoil were in different clusters. The two Japanese forests (Ashiu and IBUlim) were distantly located in different clusters based on weighted UniFrac, although they were close to each other in the tree based on the abundant Jaccard index (Fig. 2A). Furthermore, BTNmdf in Indonesia was
11
separated from the other three Indonesian forests, but was still in the third cluster. Despite these discrepancies, relative relationships were similar in the two trees based on different distance measures. In other words, geographic origin, rather than aboveground tree vegetation, also influenced the similarity, when based on weighted UniFrac, which considers DNA sequence variation (Lozupone et al., 2007). Two PCoA plots were constructed based on the two different distance measures, and both indicated differentiation of soil bacterial communities in the investigated forests (Fig. 3). In both plots, despite large variation in the aboveground environment, forests in Malaysia (except for limestone forests) and Singapore formed a distinct cluster, and limestone forests and the mangrove forest (BKOman) were relatively isolated. The two Japanese forests were well separated. However, there were also slight differences between the plots based on the different distance measures. In the plot based on the abundant Jaccard index (Fig. 3A), five Indonesian forests were separated from the rest and located on the negative side of the 2nd axis. On the other hand, in the plot based on weighted UniFrac (Fig. 3B), four Indonesian forests clustered with the forests in Malaysia and Singapore. In addition, in this plot BKOman was distantly located in the 2nd quadrant. Since the PCoA based on weighted UniFrac explained a larger fraction of the total variation (PC1: 43.2%) than that based on the abundant Jaccard index (19.1%), the former plot may be a more realistic representation of relationships among soil bacterial communities in these 28 forests.
Fig. 3. Principal coordinate analysis (PCoA) of the 28 soil samples using (A) abundant Jaccard index and (B) weighted UniFrac distance.
12
N. ITO et al.
A
B 0.8
1
0.6 weighted UniFrac
1 - Jaccard index
0.8
0.6
0.4
0.2
0.2
-2
0.4
0 -1 0 1 2 3 log10 (direct distance in kilometers)
4
-2
0 -1 0 1 2 3 log10 (direct distance in kilometers)
4
Fig. 4. Relationship between β-diversity ((A), 1 - abundant Jaccard Index, r = 0.45, P < 0.0001 and (B), weighted UniFrac, r = 0.11, P < 0.204) and direct distance in log10 (km) between forests.
1 0.8 0.6
Mantel rM
0.4 0.2 0 -0.2 -0.4 -0.6 -0.8 -1
-1
1
3
5
Distance class
7
9
Fig. 5. Mantel correlogram for the distance class with equal distance. Circles: 1 - abundant Jaccard index; squares: weighted UniFrac. Black symbols indicate significant Mantel statistics, after the progressive Bonferroni correction (α = 0.05); white, non-significant values.
Factors that influence soil bacterial composition (β-diversity) We also investigated the relationship between β-diversity (1 - abundant Jaccard index and weighted UniFrac) and direct distance between forests (Fig. 4). Pearson’s correlation coefficient was positive for both β-diversity estimates (r = 0.45, P < 0.0001, for 1 abundant Jaccard index and r = 0.11, P < 0.204 (NS), for weighted UniFrac), according to the Mantel test. This result indicated that a distance effect existed not only for adjacent forests, as noted in the UPGMA trees and PCoA plots above, but also for forests that are widely distributed over Southeast Asia and Japan. These positive correlations may reflect relatively lower β-diversities for pairs of forests less than 1 km apart and higher β-diversities for
pairs of forests close to 104 km apart (between Southeast Asia and Japanese forests), especially for 1 - abundant Jaccard index. Pairs of forests separated by 1 km to 1000 km had large variations in β-diversity, which may be due to large variation in above- and belowground environmental conditions. The Mantel correlogram analysis revealed that significant autocorrelation was detected for both β-diversities in the first distance classes for the equal-distance distance class data (Fig. 5); however, for the equal-count distance class data, none of the Mantel statistics was significant (Supplementary Fig. S6). The detected significant autocorrelation is consistent with the lower β-diversities between forests located relatively close to each other (less than 1 km), indicating that the autocorrelation is limited within a short distance. However, for the equal-distance distance class data, the number of counts (pairwise distances) was very low (3 – 6) in the distance classes with significant value (Supplementary Table S2), so that the significance in the permutation test might not be reliable due to low power of the test. Nevertheless, the positive correlation (Fig. 4) and autocorrelation (Fig. 5) indicated that biogeography of soil bacteria was not random and was influenced by factors, such as local adaptation and bacterial dispersal, which could reflect a distance effect. However, the exact factor(s) responsible for the distance effect and autocorrelation cannot be identified using the available data. Table 5 summarizes the average β-diversity within and between geographic regions and forest types. The average β-diversity within a geographic region was generally smaller than between regions (Table 5a), except for two Japanese forests, in which soil properties, especially pH, are distinct. This overall pattern was confirmed by the
Soil bacterial community structure in the tropics
13
Table 5. Average β-diversity within and between regions (a) and forest types (b) a) Region Java Sumatra Singapore West Sarawak East Sarawak Japan 1
Java
Sumatra Singapore West Sarawak East Sarawak Japan
0.5821 0.3672 0.550
0.355
0.346
0.241
0.820
0.626
0.198
0.396
0.245
0.167
0.867
0.780
0.772
0.788
0.415
0.365
0.372
0.414
0.830
0.725
0.676
0.764
0.656
0.387
0.301
0.280
0.373
0.319
0.911
0.868
0.880
0.908
0.897
0.811
0.420
0.380
0.391
0.439
0.423
0.563
1 - abundant Jaccard Index, and
2
weighted UniFrac.
Table 5. (continued) b) Forest type Acacia plantation Burned forest Dipterocarp plantation Kerangas Limestone Mangrove Mixed dipterocarp Oil palm plantation Peat swamp Secondary forest Temperate forest
Acacia Dipterocarp Mixed Oil palm Peat Secondary Temperate Burned Kerangas Limestone Mangrove plantation plantation dipterocarp plantation swamp forest forest 0.607 0.251 0.669
0.785
0.254
0.265
0.725
0.805
0.193
0.299
0.285
0.177
0.793
0.757
0.764
0.000
0.314
0.341
0.322
0.000
0.928
0.881
0.932
0.942
0.781
0.497
0.479
0.518
0.544
0.319
0.996
0.856
0.996
0.965
0.949
0.000
0.643
0.627
0.658
0.647
0.559
0.000
0.709
0.672
0.767
0.634
0.883
0.837
0.612
0.278
0.266
0.249
0.307
0.543
0.673
0.228
0.729
0.835
0.812
0.868
0.887
0.990
0.842
0.922
0.394
0.384
0.407
0.483
0.432
0.604
0.434
0.483
0.837
0.764
0.803
0.500
0.941
0.970
0.682
0.875
0.000
0.324
0.299
0.291
0.330
0.471
0.594
0.301
0.407
0.000
0.802
0.684
0.842
0.782
0.828
0.888
0.658
0.897
0.794
0.771
0.389
0.352
0.351
0.449
0.470
0.631
0.355
0.442
0.382
0.453
0.840
0.834
0.907
0.793
0.841
0.829
0.809
0.970
0.809
0.804
0.000
0.294
0.276
0.273
0.280
0.546
0.674
0.241
0.487
0.313
0.375
0.000
PERMANOVA test for 1 - abundant Jaccard index, but not for weighted Unifac, although the test result was fairly close to significance (Table 6a). This result was consistent with the distance effect detected above. On the other hand, it was not clear that average β-diversity
within a forest type was smaller than between forest types, except for dipterocarp plantations (Table 5b). These dipterocarp plantations were located close to one another in Semengo, so the lower value might not reflect forest type. However, for both β-diversity estimates, the
14
N. ITO et al. Table 6. PERMANOVA results of testing the effect of geographical region and forest type on β-diversity (1 - abundant Jaccard index and weighted UniFrac) β-diveristy: 1 - abundant Jaccard index
a) Region
β-diversity: weighted UniFrac
Source
df
SS
MS
pseudo-F
P
SS
MS
pseudo-F
Region
5
2.96
0.59
2.33
0.0001
0.48
0.10
1.26
Residual
22
5.5
0.25
1.69
0.08
Total
27
8.54
2.17
β-diveristy: 1 - abundant Jaccard index
b) Forest type
P 0.19
β-diversity: weighted UniFrac
Source
df
SS
MS
pseudo-F
P
SS
MS
pseudo-F
P
Forest type
10
4.71
0.47
2.09
0.0001
1.40
0.14
3.07
0.0004
Residual
17
3.83
0.23
0.77
0.05
Total
27
8.54
2.17
Classification of geographical region and forest type is the same as in Table 5. Table 7. Pearson’s correlation between α-diversity and environmental factors Raw data df
Chao1 index
Shannon index
Normalized data Chao1 index
Shannon index
Latitude
25
0.218
−0.085
−0.050
−0.112
Longitude
25
−0.138
−0.101
−0.128
−0.086
Altitude
25
0.156
−0.013
−0.038
−0.032
Moisture (%)
21
0.041
0.067
0.128
0.107
Organic C (%)
10
−0.131
−0.231
−0.300
−0.216
Organic N (%)
10
−0.211
−0.146
−0.279
−0.129
C/N ratio
10
0.270
−0.311
−0.323
−0.339
pH
26
0.846**
Sand %
26
−0.150
−0.310
−0.186
−0.311
Silt %
26
0.110
0.142
0.088
0.133
Clay %
26
0.102
0.299
0.176
0.311
0.761**
0.789**
0.703**
** Significant at 1% level.
Fig. 6. Relationship between α-diversity ((A), Chao1, r = 0.85, P < 0.01 and (B), Shannon indices, r = 0.79, P < 0.01) estimated for raw data and soil pH.
PERMANOVA test showed significant influence of forest type on soil bacterial community structure (Table 6b), even after eliminating singleton forest type from the test (Supplementary Table S3). These results indicate that forests, geographically adjacent and/or with the same forest type, generally tend to have similar belowground soil bacterial composition.
Relationship between α- and β-diversities and physicochemical properties in forest soil Although physicochemical properties were not measured for all the investigated forests (Supplementary Table S1), available data were used to examine the influence of the soil properties on both diversity estimates. Since Chao1 and Shannon indices were correlated as shown above, both measures of α-diversity were significantly and positively correlated with soil pH for the original and normalized data (Fig. 6, Table 7, and Supplementary Fig. S7). The forests with neutral and alkaline soil had higher α-diversity than forests with acidic soils. The other environmental factors (latitude, longitude, altitude, moisture content %, organic C content %, N content % and C/N ratio) did not show statistically significant correlation with α-diversity (Table 6). The principal coordinate analysis (PCoA) based on
Soil bacterial community structure in the tropics Table 8. Pearson correlation coefficients between the ordination score of PCoA and soil and environmental factors Factors
1 - abundant Jaccard weighted UniFrac PC1
PC2
PC1
PC2
Latitude
0.17
0.31
−0.10
0.06
Longitude
0.25
0.31
−0.19
0.02
Altitude, m
0.18
−0.20
−0.09
−0.03
Moisture content, %
0.25
0.05
−0.17
−0.18
Organic C, %
0.18
0.08
0.04
0.35
0.18
−0.13
0.15
0.24
0.08
0.79**
Organic N, % C/N ratio pH Sand, %
−0.06 0.78** −0.38
0.75** −0.04 0.48**
−0.79** 0.30
0.22 0.55**
**Significant at the 1% level.
weighted UniFrac and the 1 - abundant Jaccard index revealed a statistically significant association between soil pH and PCoA 1st axis value (Table 8), which explained 43.2% and 19.1%, respectively, of the variation in the original data (Fig. 3). These results indicated that soil pH strongly influenced not only the α-diversity of soil bacterial communities in individual forest soil but also the differentiation of soil bacterial communities among the investigated forest soils. A significant association was also detected between C/N ratio and soil texture (sand %) and PCoA 2nd axis value based on weighted UniFrac and the 1 - abundant Jaccard index (Table 8), which explained 10.8% and 10.2%, respectively, of the variation in the original data (Fig. 3), although the C/N ratio was not measured for forests in Singapore and many in Sarawak. No statistically significant relationships between the other environmental factors and axis 2 value were detected. DISCUSSION Taxonomic composition of soil bacterial communities in 28 forests in Southeast Asia and Japan Despite differences in the aboveground vegetation and environmental conditions, soil bacterial composition at the phylum level was largely similar in the 28 forests in Southeast Asia and Japan investigated in this study (Fig. 1). The top five phyla in all forests were Acidobacteria, Proteobacteria, Verrucomicrobia, Planctomycetes and Bacteroidetes, and these phyla have been reported to be major phyla in forest soil in Australia, North America and South America (Chow et al., 2002; Hackl et al., 2004; Roesch et al., 2007). Two phyla, Acidobacteria and Proteobacteria, were dominant in all forests except BKOman. These two phyla were reported to be dominant in tropical forests in the Malay Peninsula (Tripathi et al., 2012) and Sarawak (Miyashita et al., 2013) in Malaysia. The present study shows that Acidobacteria and Proteobacte-
15
ria are also predominant in the soil of tropical forests in Indonesia (Java and Sumatra) and Singapore in Southeast Asia. Acidobacteria and Proteobacteria are Gram-negative bacteria represented by many species that grow and survive under many environmental conditions (Stackebrandt et al., 1988; Quaiser et al., 2003); however, Acidobacteria are found primarily in soil (Sogins et al., 2006), while Proteobacteria are found in soil (Barns et al., 1999; Zinger et al., 2011; Wang et al., 2012), fresh water (Hugenholtz et al., 1998) and seawater (Hugenholtz et al., 1998; Venter et al., 2004; Zinger et al., 2011). In BKOman, an exceptional phylum composition was observed, characterized by high representation of Proteobacteria (≈40%) and low representation of Acidobacteria (4%). This may be related to the ability of Proteobacteria to thrive in seawater, which was detected in the mangrove forest BKOman but not in the other 27 forests studied. In addition to BKOman, in some forests, such as JASoil, BAU2nd, BAUlim, NIAlim and IBUlim, class compositions in Acidobacteria and Proteobacteria were different from those in the majority of the forests (Supplementary Figs. S3 and S4, A and B). In these six forests, the frequencies of Acidobacteria classes Gp1 and Gp2 were lower and other subgroups (Gp4, Gp5 and Gp6) were higher than in the other forests (Supplementary Fig. S4A). In BKOman, the frequency of Acidobacteria class Gp10 was distinctively high. The three limestone forest soils (BAUlim, NIAlim and IBUlim) had a relatively high frequency of Acidobacteria classes Gp5 and Gp6. On the other hand, in these six forests, except for JASoil, the frequency of Alphaproteobacteria was lower and that of Betaproteobacteria was higher than in the other forests (Supplementary Fig. S4B). Thus, the changes in class composition in Acidobacteria and Proteobacteria seemed to be linked, suggesting some kind of interaction between these dominant phyla, although there are many other possible explanations. This possibility needs to be further examined in future studies. Besides these six forests of unusual class composition of Acidobacteria and Proteobacteria, in the other forest soils in Southeast Asia and Japan, Acidobacteria Gp1 and Gp2 accounted for approximately 90% of Acidobacteria classes, and Alphaproteobacteria was most frequent among the Proteobacteria classes. These class compositions of Acidobacteria and Proteobacteria in Asia differ from those in other regions of the world. For example, in a study of 16S rRNA sequences in 32 soil samples from diverse environments in Europe and Australia in which Acidobacteria and Proteobacteria were the dominant phyla, class composition depended on the geographic origin of the soil sample, and Gp1, Gp4 and Gp6 of Acidobacteria and Alpha- and Betaproteobacteria of Proteobacteria were the dominant classes (Newton et al., 2011). A study on Acidobacteria by T-RFLP and DGGE methods on
16
N. ITO et al.
the 16S rRNA gene for 27 grassland and 30 forest plots in Germany showed that class composition of Acidobacteria differed among plots within an aboveground vegetation and between aboveground vegetations, and that Gp6 and Gp1 were dominant in grassland and forest soils, respectively (Janssen, 2006). Furthermore, pyrosequencing and clone library analyses of Acidobacteria 16S rRNA sequences from 87 soil samples from North and South America showed that Gp1 was dominant (Naether et al., 2012). These results indicate that Gp1 and Gp6 of Acidobacteria and Alpha- and Betaproteobacteria of Proteobacteria are dominant in soils from Europe and Australia. Although Alpha- and Betaproteobacteria of Proteobacteria were also dominant in Asian forest soils in this study, the predominance of Acidobacteria Gp2 and Gp5 may be a characteristic of Asian forest soils. Diversity of soil bacterial communities in investigated forests A cluster analysis showed that the number of OTUs and the number of species estimated (Chao1 index) were of the order of 103 in all forests (Tables 3 and 4), in agreement with the number of species estimated by PCR amplification of the 16S rRNA gene region (Curtis and Sloan, 2005; Gans et al., 2005). However, the actual diversity of the soil bacterial community in these samples is likely to be higher than 103, because, in all cases, the rarefaction curves for these data sets did not plateau. Many OTUs were singletons (43.8% on average) with few sequence reads per OTU per forest; this is consistent with the “rare biosphere” concept described in previous analyses of bacterial communities in deep-sea water (Sogins et al., 2006), agricultural soil (Ashby et al., 2007), prairie soil (Elshahed et al., 2008) and forest soil in Southeast Asia (Miyashita et al., 2013). Soil bacterial species diversity (the number of OTUs and Chao1 index) varied among forests. However, the relationship among forests with respect to the level of bacterial species diversity was not influenced by 16S rRNA sequence read variation, as indicated by a highly significant Pearson’s correlation between diversity estimates for raw and normalized data. These results suggested that factor(s) inherent to each forest affected the level of bacterial species diversity. The bacterial species diversity also varied by geographic area: Java samples had high diversity, and Singapore samples had low diversity. However, bacterial species diversity in Japanese forests was similar to that of the tropical forests. In fact, bacterial species diversity in Mt. Ibuki (IBUlim) was the third-highest value observed in this study (Table 4). Furthermore, forest type was not associated with bacterial species diversity. Thus, it can be concluded that aboveground condition did not directly influence soil bacterial species diversity in these forests. Bacterial diversity was reported to correlate with soil pH in American forests (Fierer and Jackson, 2006; Lauber
et al., 2009), such that diversity was lowest in acid soil, highest at nearly neutral pH, and then decreased with increasing pH. This suggests that extreme pH is stressful for unicellular bacteria. Although soil with pH > 8 was not detected in this study, Chao1 and Shannon indices correlated positively with soil pH (Fig. 6); the diversity increased as pH became alkaline (near pH 8) from acid (near pH 3), as reported for tropical forests in Malaysia (Tripathi et al., 2012; Miyashita et al., 2013). However, the influence of pH on soil bacterial diversity was not straightforward, since there was considerable variation in the Chao1 index among five forest soils with high pH (around 7), but slight variation in soil with lower pH (Fig. 6A). Although MARaca5 and BAUlim forest soils had almost the same pH (6.60 and 6.61, respectively, Supplementary Table S1), their Chao1 indices were 2404.5 and 6051.8, respectively. Furthermore, BAU2nd soil had a high pH (7.78), but its Chao1 index was not high (3702.4). As indicated by the positive correlation between Chao1 and Shannon indices, as pH increased, the Shannon index increased. This result suggests that as the number of species increases in a soil environment with pH, the density of species (the number of individuals per species) decreases. The effect of pH on the Shannon index was not the same as on the Chao1 index; there was little variation in the Shannon index among soils with higher pH, while considerable variation existed among soils with lower pH (Fig. 6B). This would indicate that the effect of pH on species density varied among forest soils at low pH. Taken together, these results suggest that the influence of pH differs between the two measures of α-diversity and that factor(s) other than pH influence the α-diversity of soil bacteria, although the effect of pH is the strongest. Relationship among soil bacterial community structures in investigated forests In the UPGMA dendrogram based on the abundant Jaccard index, which took account of the similarity of OTU composition, geographically adjacent forests formed clusters, except for BKOman, NIAlim, NIAmdf and JASoil (Fig. 2A). This indicated that soil bacterial composition at the species (OTU) level was related to geographical distance, but not to the type of aboveground vegetation. In contrast, the topology of the dendrogram based on weighted UniFrac distance (Fig. 2B) differed, in that the 28 forests were divided into three large clusters, not necessarily in agreement with the geographical relatedness among forests. Because weighted UniFrac considers phylogenetic relationships among OTUs, the difference in the topology between the two dendrograms indicated that OTUs in geographically adjacent forests did not necessarily share similar evolutionary history in terms of sequence evolution of 16S rRNA. The PCoA analysis based on the abundant Jaccard index divided the forests roughly
Soil bacterial community structure in the tropics into four quadrants (Fig. 3A), while that based on the weighted UniFrac distance showed three distantly separated groups (Fig. 3B), corresponding with the topology of the weighted UniFrac-based UPGMA dendrogram. The difference in the scattering of both PCoA plots also reflects the fact that they are based on different distance measures. It is worth noting that in the dendrograms based on the abundant Jaccard index or weighted UniFrac, and in the PCoA plot based on weighted UniFrac, the BKOman forest was clearly separated from the other forests. This result suggests that the structure of the bacterial community in BKOman is unique at the level of species. In particular, 94.5% of the OTUs in BKOman were unique to BKOman, in the context of this study. A mixture of fresh water and seawater irrigates the soil in BKOman, a mangrove forest, and it was reported that salinity influences the structure of bacterial communities (Lozupone and Knight, 2007). This could promote growth of salttolerant and anaerobic bacteria in BKOman soil, contributing to the unique structure of its bacterial community at the levels of phylum and species. Consistent with this idea, of the OTUs detected only in BKOman, the two most abundant were Cyanobacteria, which is also a major bacterial phylum of seawater (Gotelli and Colwell, 2001; Gans et al., 2005). However, many characteristics of BKOman other than the salinity of its soil, such as water temperature and oxygen level associated with tidal action, could contribute to the unique structure of its soil bacterial community. To examine the uniqueness of soil bacterial community structure in BKOman and the effect of salinity on the soil bacterial community structure, it would be necessary not only to study other mangrove forest soils, but also to measure salinity in other forest soils in future studies. Environmental properties and mechanisms that influence soil bacterial species composition To investigate the association between soil bacterial community structure, especially species composition, and soil environmental properties, Pearson’s correlation coefficient was examined for the 1st and 2nd ordination scores of PCoA (Table 8). A highly significant correlation was detected between the 1st ordinate values and soil pH, with a less strong correlation observed for eight other physicochemical properties of soil. We conclude that soil pH strongly influences α-diversity of soil bacteria as well as composition of the soil bacterial community at the level of species, most likely because the pH of the environment has a large effect on bacterial growth and survival. However, in the two PCoA plots based on different β-diversity, three limestone forests, IBUlim, BAUlim and NIAlim, in addition to Bau2nd, lie close together in a single cluster. Although limestone usually makes soil alkaline, NIAlim had acidic soil (Supplementary Table S1), indi-
17
cating that pH as well as other properties of soil influence soil bacterial composition. In addition to pH, a highly significant Pearson’s correlation was detected between the 2nd ordination values and C/N ratio and soil texture (sand %), indicating that these two variables also affect the composition of the soil bacterial community. Indeed, it has been reported that bacterial diversity and composition differed between Bornean rainforest soils with different C/N ratios and soil textures, suggesting the influence of the composition of carbon and nitrogen, and soil pH, in addition to oxygen availability (Russo et al., 2012). However, in this study, C/N ratio and soil texture were not associated with α-diversity, but instead were correlated with the 2nd ordinate values of PCoA, which have secondary importance. Since C/N ratio was not measured for all the forests studied, this discussion of the C/N ratio should be taken with caution. Another important observation of the present study was the statistically significant correlation between measures of β-diversity (1 - abundant Jaccard index) and physical distance between forests (Fig. 4), confirmed by the PERMANOVA test for the effect of geographical region on soil bacterial community structure (Table 6). This implies that the biogeography of soil bacterial composition is non-random, such that bacterial composition is influenced by local adaptation and limited dispersal (Hughes Martiny et al., 2006). In our previous studies, multiple samples were taken ca. 10 meters apart in each of six forests in Sarawak and Japan, with the result that bacterial species composition was more similar within a forest than between forests (Miyashita et al., 2013; Miyashita, 2015). These previous results are consistent with a very small spatial scale distance effect, probably reflecting local adaptation and bacterial movement within a forest. In this study, Mantel correlogram analysis revealed a significant autocorrelation in β-diversity between forests further separated (up to 28.8 km for the 1 - abundant Jaccard index and 0.55 km for weighted UniFrac, Fig. 5 and Supplementary Table S2), suggesting a relatively longer scale of distance effect. However, as mentioned above, the significance of this autocorrelation should be taken with some caution, because of the low power of the permutation test. On the other hand, the distance effect in this study is over a much longer spatial scale (up to 10,000 km), such that bacterial dispersal was probably limited, and bacteria or bacterial spores could be subject to long-distance passive dispersal by wind, rain or contact with birds (Hughes Martiny et al., 2006). In fact, some OTUs were shared between forests in Southeast Asia and Japan, and β-diversity estimates were less than the maximum value of one. Actually, in the Japanese forests Ashiu and IBUlim, 1- the percentage of forestspecific OTUs, which is the percentage of OTUs in a forest shared with other forests, was 0.575 and 0.3 (Table 3), respectively, indicating that a high percentage of OTUs in
18
N. ITO et al.
Japanese forests were shared with the Southeast Asian forests, given that the β-diversity estimates between the two Japanese forests were relatively high (Table 5). Although it was not clear from the average values of β-diversity (Table 5), PERMANOVA tests also indicated that aboveground forest type has a significant influence on soil bacterial species composition (β-diversity, Table 6), indicating that forests with the same forest type are likely to have similar soil bacterial composition. On the other hand, it was suggested in this study that soil bacterial species diversity (α-diversity) was not strongly influenced by aboveground forest type (Tables 3 and 4c). These results may be relevant to the finding from a study on temperate grassland sites that aboveground plant α-diversity can predict patterns of belowground soil bacterial composition (β-diversity) but not the pattern of α-diversity (Prober et al., 2015). However, since a diverse range of aboveground forest types was investigated in this study, it is reasonable to assume that many different factors other than plant α-diversity might have contributed to our results. Therefore, it would be very informative to take into account aboveground plant α-diversity for studies of soil bacterial community structure in future. Taken together, the results of this study suggest that soil bacterial community structure (species composition in particular) was influenced significantly by soil physicochemical properties, mainly pH (and salinity), which influence bacterial growth; by distance effects, which might be related to local adaptation, evolutionary history from common ancestors and migration of bacterial species among forests; and finally by aboveground forest type. However, one caveat is that the range of soil pH examined in this study was not wide, because high alkaline soil was not included. Also, the influence of salinity was inferred from the bacterial community structure in BKOman, from which the effect of mangrove itself was not eliminated. Moreover, salinity itself was not measured in this study, and should be quantified in future studies. In addition, the underlying causes of the distance effect observed here are not clear. In future studies on the environmental, ecological and historical factors that influence soil bacterial community structure in tropical and temperate forests, it will be important to analyze forests that are more ecologically diverse and more geographically distant than those in the present study. This work was supported by a grant (D-0901) from the Ministry of the Environment of Japan. This paper is contribution #614 from the Laboratory of Plant Genetics, Kyoto University, Japan.
REFERENCES Anderson, M. J. (2001) A new method for non-parametric multivariate analysis of variance. Austral Ecol. 26, 32–46.
Anderson, M. J., and Willis, T. J. (2003) Canonical analysis of principal coordinates: A useful method of constrained ordination for ecology. Ecology 84, 511–525. Ashby, M. N., Rine, J., Mongodin, E. F., Nelson, K. E., and Dimster-Denk, D. (2007) Serial analysis of rRNA genes and the unexpected dominance of rare members of microbial communities. Appl. Environ. Microbiol. 73, 4532–4542. Bais, H. P., Weir, T. L., Perry, L. G., Gilroy, S., and Vivanco, J. M. (2006) The role of root exudates in rhizosphere interactions with plants and other organisms. Annu. Rev. Plant Biol. 57, 233–266. Baker, G. C., Smith, J. J., and Cowan, D. A. (2003) Review and re-analysis of domain-specific 16S primers. J. Microbiol. Methods 55, 541–555. Barns, S. M., Takala, S. L., and Kuske, C. R. (1999) Wide distribution and diversity of members of the bacterial kingdom Acidobacterium in the environment. Appl. Environ. Microbiol. 65, 1731–1737. Bond, P. L., Hugenholtz, P., Keller, J., and Blackall, L. L. (1995) Bacterial community structures of phosphate-removing and non-phosphate-removing activated sludges from sequencing batch reactors. Appl. Environ. Microbiol. 61, 1910–1916. Borneman, J., and Triplett, E. W. (1997) Molecular microbial diversity in soils from Eastern Amazonia: Evidence for unusual microorganisms and microbial population shifts associated with deforestation. Appl. Environ. Microbiol. 63, 2647–2653. Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., Fierer, N., Peña, A. G., Goodrich, J. K., Gordon, J. I., et al. (2010a) QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. Caporaso, J. G., Bittinger, K., Bushman, F. D., DeSantis, T. Z., Andersen, G. L., and Knight, R. (2010b) PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26, 266–267. Chao, A. (1984) Nonparametric estimation of the number of classes in a population. Scand. J. Statist. 11, 265–270. Chao, A., Chazdone, R. L., Colwell, R. K., and Shen, T. J. (2005) A new statistical approach for assessing similarity of species composition with incidence and abundance data. Ecol. Lett. 8, 148–159. Chapin, F. S. (1980) The mineral nutrition of wild plants. Annu. Rev. Ecol. Syst. 11, 233–260. Chow, M. L., Radomski, C. C., McDermott, J. M., Davies, J., and Axelrood, P. E. (2002) Molecular characterization of bacterial diversity in Lodgepole pine (Pinus contorta) rhizosphere soils from British Columbia forest soils differing in disturbance and geographic source. FEMS Microbiol. Ecol. 42, 347–357. Cole, J. R., Wang, Q., Cardenas, E., Fish, J., Chai, B., Farris, R. J., Kulam-Syed-Mohideen, A. S., McGarrell, D. M., Marsh, T., Garrity, G. M., et al. (2009) The Ribosomal Database Project: improved alignments and new tools for rRNA analysis. Nucleic Acids Res. 37, D141–D145. Curtis, T. P., and Sloan, W. T. (2005) Exploring microbial diversity - a vast below. Science 309, 1331–1333. DeAngelis, K. M., Silver, W. L., Thompson, A. W., and Firestone, M. K. (2010) Microbial communities acclimate to recurring changes in soil redox potential status. Environ. Microbiol. 12, 3137–3149. Elshahed, M. S., Youssef, N. H., Spain, A. M., Sheik, C., Najar, F. Z., Sukharnikov, L. O., Roe, B. A., Davis, J. P., Schloss, P. D., Bailey, V. L., et al. (2008) Novelty and uniqueness patterns of rare members of the soil biosphere. Appl. Environ.
Soil bacterial community structure in the tropics Microbiol. 74, 5422–5428. Fierer, N., and Jackson, R. B. (2006) The diversity and biogeography of soil bacterial communities. Proc. Natl. Acad. Sci. USA 103, 626–631. Fiere, N., Craine, J. M., McLauchlan, K., and Schimel, J. P. (2005) Litter quality and the temperature sensitivity of decomposition. Ecology 86, 320–326. Fierer, N., Lauber, C. L., Ramirez, K. S., Zaneveld, J., Bradford, M. A., and Knight, R. (2012) Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients. ISME J. 6, 1007– 1017. Gans, J., Wolinsky, M., and Dunbar, J. (2005) Computational improvements reveal great bacterial diversity and high metal toxicity in soil. Science 309, 1387–1390. Gotelli, N. J., and Colwell, R. K. (2001) Quantifying biodiversity: procedures and pitfalls in the measurement and comparison of species richness. Ecol. Lett. 4, 379–391. Hackl, E., Zechmeister-Boltenstern, S., Bodrossy, L., and Sessitsch, A. (2004) Comparison of diversities and compositions of bacterial populations inhabiting natural forest soils. Appl. Environ. Microbiol. 70, 5057–5065. Harris, J. (2009) Soil microbial communities and restoration ecology: Facilitators or followers? Science 325, 573–574. Hättenschwiler, S., Tiunov, A. V., and Scheu, S. (2005) Biodiversity and litter decomposition in terrestrial ecosystems. Annu. Rev. Ecol. Evol. Syst. 36, 191–218. Högberg, P., Nordgren, A., Buchmann, N., Taylor, A. F., Ekblad, A., Högberg, M. N., Nyberg, G., Ottosson-Löfvenius, M., and Read, D. J. (2001) Large-scale forest girdling shows that current photosynthesis drives soil respiration. Nature 411, 789–792. Huber, J. A., Welch, D. B. M., Morrison, H. G., Huse, S. M., Neal, P. R., Butterfield, D. A., and Sogin, M. L. (2007) Microbial population structures in the deep marine biosphere. Science 318, 97–100. Hugenholtz, P., Pitulle, C., Hershberger, K. L., and Pace, N. R. (1998) Novel division level bacterial diversity in a Yellowstone hot spring. J. Bacteriol. 180, 366–376. Hughes Martiny, J. B., Bohannan, B. J. M., Brown, J. H., Colwell, R. K., Fuhrman, J. A., Green, J. L., Horner-Devine, M. C., Kane, M., Krumins, L. A., Kuske, C. R., et al. (2006) Microbial biogeography: putting microorganisms on the map. Nat. Rev. Microbiol. 4, 102–112. Janssen, P. H. (2006) Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 72, 1719–1728. Jesus, E. D. C., Marsh, T. L., Tiedje, J. M., and Moreira, F. M. D. S. (2009) Changes in land use alter the structure of bacterial communities in Western Amazon soils. ISME J. 3, 1004–1011. Kowalchuk, G. A., and Stephen, J. R. (2001) Ammonia-oxidizing bacteria: a model for molecular microbial ecology. Ann. Rev. Microbiol. 55, 485–529. Lauber, C. L., Hamady, M., Knight, R., and Fierer, N. (2009) Pyrosequencing-based assessment of soil pH as a predictor of soil bacterial community structure at the continental scale. Appl. Environ. Microbiol. 75, 5111–5120. Legendre, P., and Legendre, L. (1998) Numerical Ecology, Second English Edition, pp. 707–785. Elsevier, Paris. Lozupone, C. A., and Knight, R. (2007) Global patterns in bacterial diversity. Proc. Natl. Acad. Sci. USA 104, 11436–11440. Lozupone, C. A., Hamady, M., Kelley, S. T., and Knight, R. (2007) Quantitative and qualitative β diversity measures lead to different insights into factors that structure microbial com-
19
munities. Appl. Environ. Microbiol. 73, 1576–1585. Miyashita, N. T. (2015) Contrasting soil bacterial community structure between the phyla Acidobacteria and Proteobacteria in tropical Southeast Asian and temperate Japanese forests. Genes Genet. Syst. 90, 61–77. Miyashita, N. T., Iwanaga, H., Charles, S., Diway, B., Sabang, J., and Chong, L. (2013) Soil bacterial community structure in five tropical forests in Malaysia and one temprate forest in Japan revealed by pyrosequencing analyses of 16S rRNA gene sequence variation. Genes Genet. Syst. 88, 93–103. Myers, N., Mittermeier, R. A., Mittermeier, C. G., da Fonseca, G. A. B., and Kent, J. (2000) Biodiversity hotspots for conservation priorities. Nature 403, 853–858. Naether, A., Foesel, B. U., Naegele, V., Wüst, P. K., Weinert, J., Bonkowski, M., Alt, F., Oelmann, Y., Polle, A., Lohaus, G., et al. (2012) Environmental factors affect Acidobacterial communities below the subgroup level in grassland and forest soils. Appl. Environ. Microbiol. 78, 7398–7406. Newton, R. J., Jones, S. E., Eiler, A., McMahon, K. D., and Bertilsson, S. (2011) A guide to the natural history of freshwater lake bacteria. Microbiol. Mol. Biol. Rev. 75, 14–49. Nüsslein, K., and Tiedje, J. M. (1999) Soil bacterial community shift correlated with change from forest to pasture vegetation in a tropical soil. Appl. Environ. Microbiol. 65, 3622– 3626. Prober, S. M., Leff, J. W., Bates, S. T., Borer, E. T., Firn, J., Harpole, W. S., Lind, E. M., Seabloom, E. W., Adler, P. B., Bakker, J. D., et al. (2015) Plant diversity predicts beta but not alpha diversity of soil microbes across grasslands worldwide. Ecol. Lett. 18, 85–95. Quaiser, A., Ochsenreiter, T., Lanz, C., Schuster, S. C., Treusch, A. H., Eck, J., and Schleper, C. (2003) Acidobacteria form a coherent but highly diverse group within the bacterial domain: evidence from environmental genomics. Mol. Microbiol. 50, 563–575. Quince, C., Lanzén, A., Curtis, T. P., Davenport, R. J., Hall, N., Head, I. M., Read, L. F., and Sloan, W. (2009) Accurate determination of microbial diversity from 454 pyrosequencing data. Nat. Methods 6, 639–641. Roesch, L. F. W., Fulthorpe, R. R., Riva, A., Casella, G., Hadwin, A. K. M., Kent, A. D., Daroub, S. H., Camargo, F. A. O., Farmerie, W. G., and Triplett, E. W. (2007) Pyrosequencing enumerates and contrasts soil microbial diversity. ISME J. 1, 283–290. Rosenberg, M. S., and Anderson, C. D. (2011) PASSaGE: Pattern Analysis, Spatial Statistics and Geographic Exegesis. Version 2. Methods Ecol. Evol. 2, 229–232. Russo, S. E., Legge, R., Weber, K. A., Brodie, E. L., Goldfarb, K. C., Benson, A. K., and Tan, S. (2012) Bacterial community structure of contrasting soils underlying Bornean rain forests: Inferences from microarray and next-generation sequencing methods. Soil Biol. Biochem. 55, 48–59. Sanders, H. (1968) Marine benthic diversity: A comparative study. Am. Nat. 102, 243–282. Schloss, P. D., and Handelsman, J. (2005) Introducing DOTUR, a computer program for defining operational taxonomic units and estimating species richness. Appl. Environ. Microbiol. 71, 1501–1506. Shannon, C. E. (1948) A mathematical theory of communication. Bell Syst. Tech. J. 27, 379–423. Slabbert, E., Kongor, R. Y., Esler, K. J., and Jacobs, K. (2010) Microbial diversity and community structure in Fynbos soil. Mol. Ecol. 19, 1031–1041. Sogins, M. L., Morrison, H. G., Huber, J. A., Welch, D. M., Huse, S. M., Neal, P. R., Arrieta, J. M., and Herndl, G. J. (2006)
20
N. ITO et al.
Microbial diversity in the deep sea and the underexplored “rare biosphere”. Proc. Natl. Acad. Sci. USA 103, 12115– 12120. Sokal, R. R. (1986) Spatial data analysis and historical processes. In Data Analysis and Informatics, IV, (eds.: E. Diday et al.) pp. 29–43. North-Holland, Amsterdam. Sokal, R. R., and Rohlf, F. J. (2000) Biometry. Third Edition. W. H. Freeman and Company, New York, 803 p. Stackebrandt, E., Murray, R. G. E., and Trüper, H. G. (1988) Proteobacteria classis nov., a name for the phylogenetic taxon that includes the “purple bacteria and their relatives.” Int. J. Syst. Bacteriol. 38, 321–325. Torsvik, V., Øvreas, L., and Thingstad, T. F. (2002) Prokaryotic diversity—magnitude, dynamics, and controlling factors. Science 296, 1064–1066. Tripathi, B. M., Kim, M., Singh, D., Lee-Cruz, L., Lai-Hoe, A., Ainuddin, A. N., Go, R., Rahim, R. A., Husni, M. H. A., Chun, J., et al. (2012) Tropical soil bacterial communities in Malaysia: pH dominates in the equatorial tropics too. Microb. Ecol. 64, 474–484. van der Heijden, M. G. A., Bardgett, R. D., and van Straalen, N. M. (2008) The unseen majority: soil microbes as drivers of plant diversity and productivity in terrestrial ecosystems. Ecol. Lett. 11, 296–310. Venter, J. C., Remington, K., Heidelberg, J. F., Halpern, A. L., Rusch, D., Eisen, J. A., Wu, D., Paulsen, I., Nelson, K. E., Nelson, W., et al. (2004) Environmental genome shotgun sequencing of the Sargasso Sea. Science 304, 66–74. Wang, Y., Sheng, H., He, Y., Wu, J., Jiang, Y., Tam, N. F. Y., and Zhou, H. W. (2012) Comparison of the levels of bacterial
diversity in freshwater, intertidal wetland, and marine sediments by using millions of Illumina tags. Appl. Environ. Microbiol. 78, 8264–8271. Wardle, D. A., Bardgett, R. D., Klironomos, J. N., Setälä, H., van der Putton, W. H., and Wall, D. H. (2004) Ecological linkages between aboveground and belowground biota. Science 304, 1629–1633. Whitman, W., Coleman, D., and Wiebe, W. (1998) Prokaryotes: the unseen majority. Proc. Natl. Acad. Sci. USA 95, 6578– 6583. Wieland, G., Neumann, R., and Backhaus, H. (2001) Variation of microbial communities in soil, rhizosphere, and rhizoplane in response to crop species, soil type, and crop development. Appl. Environ. Microbiol. 67, 5849–5854. Wu, L., Feinstein, L. M., Valverde-Barrantes, O., Kershner, M. W., Leff, L. G., and Blackwood, C. B. (2011) Placing the effects of leaf litter diversity on saprotrophic microorganisms in the context of leaf type and habitat. Microb. Ecol. 61, 399–409. Yarwood, S. A., Bottomley, P. J., and Myrold, D. D. (2010) Soil microbial communities associated with Douglas-fir and red alder stands at high- and low- productivity forest sites in Oregon, USA. Microb. Ecol. 60, 606–617. Zinger, L., Amaral-Zettler, L. A., Fuhrman, J. A., HornerDevine, M. C., Huse, S. M., Welch, D. B. M., Martiny, J. B. H., Sogin, M., Boetius, A., and Ramette, A. (2011) Global patterns of bacterial beta-diversity in seafloor and seawater ecosystems. PLoS One 6, e24570, doi: 10.1371/journal. pone.0024570.
