Species may co-occur due to responses to similar environmental conditions, biological associations, or simply because of coincident geographical distributions. Disentangling patterns of co-occurrence and potential biotic and abiotic interactions is crucial to understand ecosystem function. Here we used DNA metabarcoding data from litter and mineral soils collected from a longitudinal transect in Amazonia to explore patterns of co-occurrence. We compared data from different Amazonian habitat types, each with a characteristic biota and environmental conditions. These included non-flooded rainforests (terra-firme), forests seasonally flooded by fertile white waters (várzeas) or by unfertile black waters (igapós), and open areas associated with white sand soil (campinas). We ran co-occurrence network analyses based on null models and Spearman correlation for all samples and for each habitat separately. We found that one third of all operational taxonomic units (OTUs) were bacteria and two thirds were eukaryotes. The resulting networks were nevertheless mostly composed of bacteria, with fewer fungi, protists, and metazoans. Considering the functional traits of the OTUs, there is a combination of metabolism modes including respiration and fermentation for bacteria, and a high frequency of saprotrophic fungi (those that feed on dead organic matter), indicating a high turnover of organic material. The organic carbon and base saturation indices were important in the co-occurrences in Amazonian networks, whereas several other soil properties were important for the co-exclusion. Different habitats had similar network properties with some variation in terms of modularity, probably associated with flooding pulse. We show that Amazonian micro-organism communities form highly interconnected co-occurrence and co-exclusion networks, which highlights the importance of complex biotic and abiotic interactions in explaining the outstanding biodiversity of the region.
Original Article
Assessing biotic and abiotic interactions of micro-organisms in Amazonia through co-occurrence networks and DNA metabarcoding
https://doi.org/10.21203/rs.3.rs-184472/v1
This work is licensed under a CC BY 4.0 License
Journal Publication
published 18 Feb, 2021
You are reading this latest preprint version
Bacteria
Biodiversity
Fungi
Metabarcoding
Protists
Tropics
Patterns of co-occurrence among species have traditionally been explained by the species’ similar environmental preferences or tolerances, such as pH range, temperature, habitat type, and similar biogeographic histories generating coincident distributions [1, 2, 3]. Species co-occurrence can also be related to biotic interactions, such as predation and parasitism [4]. Such patterns can be described with the use of ecological network analyses, an important tool to investigate community structure and the function and stability of ecosystems [5]. Identifying co-occurrence patterns is essential to grasp the potential ecological interactions and their impact on biodiversity and ecosystem functioning [6, 7].
In the tropics, patterns of co-occurrence can be extremely complex due to the high diversity of species inhabiting different kinds of habitats. The study of species interactions in tropical environments has been largely neglected (but see [8, 9, 10, 11, 12]). Most studies of co-occurrences have focused on groups of symbiotic organisms, such as trees and mycorrhizas [13, 14, 15] and trees and seed dispersers [16, 17], and are largely centered around mutualism [18] and parasitism [19, 20, 21] . However, many other interactions among organisms may play a crucial role in biotic and/or abiotic interactions in the tropics, such as those related to nutrient cycling and organic matter decomposition. Unveiling these potential interactions is important to shed light on the ecological dynamics and ecosystem functions of biodiverse regions.
Micro-organisms account for the majority of biodiversity in any environment, including the tropics [22], being crucial for ecosystem dynamics such as nutrient recycling [23, 24] and much of organic soil decomposition [25, 26]. A negative correlation between Amazonian soil organic carbon and soil biodiversity has been found, which could be explained by high soil biodiversity accelerating organic decomposition and subsequently locking carbon in above-ground biomass [27].
Amazonia is the world’s largest tropical forest, covering ca. 3.6% of the global surface while comprising approximately 40% of all rainforest. It harbors 10% of the world’s known biodiversity [28] and potentially a quarter of the world's terrestrial species [29]. Amazonia is home to a large number of coexisting and potentially interacting species, which probably make up one of the most complex webs of life on Earth. As such, Amazonia provides a wide range of ecosystem services through its high above- and below-ground biodiversity [30], including water cycling and carbon storage [31, 32, 33], which are mediated through biotic and abiotic interactions. Amazonia can be divided into four main phytophysiognomies (hereafter habitats): non-flooded rainforests (terra-firme), forests that are seasonally flooded either by fertile white waters (várzeas) or by unfertile black waters (igapós), and naturally open areas associated with white sand soil (campinas). The geographical location, the degree of connectivity and isolation, and the soil characteristics of those four main habitat types are important determinants of their biota [34, 35, 36], which may consequently determine species co-occurrence patterns.
Recent studies have shown that Amazonian habitat types influence the composition of microbial communities [36]. Micro-organism diversity is the richest in campina habitats [36, 37, in contrast to known patterns of macro-organisms that show low diversity in these habitats [38], 39]. Campinas have a highly specialized biota, probably explained by multiple stressors that affect the habitat, such as poor soils and insular distribution [38]. Likewise, the flooded habitats have a highly specialized biota due to stress induced by seasonal flooding [40]. Furthermore, the soil characteristics of all these habitats are different and explain in part the composition of the biota, with pH and organic carbon being the most important soil properties for micro-organism richness and turnover [27]. Most micro-organisms seem to occur in just one locality and habitat type [36], indicating a potentially high degree of local specialization or low detectability. The diversity of microbial communities and their composition in Amazonia have already been investigated thanks to advances in high-throughput sequencing (HTS) methods [41, 42, 43, 44, 36, 27, 45]. However, most of these studies focused on single properties of the communities, such as alpha and beta diversity, and less on the potential interactions between micro-organisms as revealed through ecological networks. The integration of graph theory and network analysis into metabarcoding studies allows the exploration of complex community interactions beyond descriptive diversity patterns [46]. Furthermore, graph properties such as density (the rate of connection between nodes), transitivity (probability for the network to have adjacent nodes interconnected), and modularity (the degree of segregation between different modules) from networks may allow an assessment of an ecosystem’s integrity and resilience [47, 48, 49, 50].
Here we explore the potential interactions of micro-organisms (bacteria, fungi, and protists) and some plants and metazoans in soil using co-occurrence network analysis in four localities in Amazonia. These localities are separated by wide rivers and long distances (Fig. 1) and have previously been examined for diversity using both richness and the effective number of operational taxonomic units (OTUs [51]) [36, 27]. We also test for the importance of environmental filters based on the soil physico-chemical properties quantified in each sampled plot [27]. Based on previously identified patterns of OTU richness and distribution in each habitat [36], we expect that the co-occurrence networks will show different structures depending on habitat type. We formulated the following specific hypotheses, based on the current literature: H1: The Amazonia-wide network will be mostly composed of organisms associated with organic decomposition; H2: Environmental soil properties, especially pH and organic carbon [27, 52], will be the most important factors, acting as key nodes, to explain co-occurrence in all the networks; and H3: The habitat-specific networks will show different network structures, with the presence of flood pulse (the periodic inundation of floodplains along certain rivers) as the main factor explaining network properties. We expect that the high environmental stress in seasonal flooding forests will act as an abiotic filter, keeping the same set of specialized organisms that co-occur for long periods of time, thus, resulting in a dense and highly connected network (high transitivity, low modularity).
Sampling and sequencing – We sampled four main localities along a longitudinal transect (Fig. 1) in November 2015, a month of the dry season without inundation of seasonally flooded forests (around 5 months after inundation peak). The data was previous published in Ritter et al. [27, 36, 37]. In each locality, we sampled all of the four habitats (terra-firme, campinas, várzea, and igapó) whenever they were present (Fig. 1). We used data from soil and litter metabarcoding of the V3–V4 region of the ribosomal small subunit 16S rRNA, targeting bacteria with primers described in Klindworth et al. [53]; the V7 region of the ribosomal small subunit 18S rRNA with primers described in Guardiola et al. [54], as well as the mitochondrial cytochrome c oxidase subunit I (COI) gene, using the forward primer described in Leray et al. [55] and the reverse primer described in Meyer [56], both targeting eukaryotes. Details of sampling and sequencing can be found in Ritter et al. [36, 37]. Briefly, we filtered the low-quality sequences, discarded the chimeras, and clustered sequences into OTUs at 97% similarity using a ‘greedy’ algorithm that performs chimera filtering and OTU clustering simultaneously with the USEARCH / UPARSE v. 9.0.2132 Illumina paired reads pipeline [57].
Classification of OTUs – We classified OTUs into two categories by habitat occurrence: “specialist” if they occurred in only one habitat (e.g., only in campinas) and as “generalist” if they occurred in at least three habitats. These definitions follow earlier studies on the distribution of organisms in Amazonia as well as previous analyses of our data [36]. To classify the OTUs and produce the figures, we used the tidyverse v.1.2.1 [58] and ggplot2 v.3.2.1 [59] R packages in R v.3.6.0 [60].
We assigned OTUs according to the taxonomic classification published in Ritter et al. [36]. As far as possible, we used the class or order for the classification of the OTU. We assessed the bacterial OTUs for the functional traits using the FAPROTAX database [61], which we applied to all bacteria in our networks. Since several bacteria were assigned to multiple functional traits, we reviewed all classifications and assigned them to the most probable or the dominant functional trait based on the literature and our professional experience. For protists, we used the study of Adl et al. [62] as a reference. For fungi, we used the functional classification published for the same data as in Ritter et al. [63].
Co-occurrence network analysis – Network analyses were performed using abundance matrices as recommended [3] and the null model strategy described in Connor et al. [64], which is regarded as the most reliable strategy for inferring both co-occurrence and co-exclusion relationships in HTS datasets of microbial organisms [65]. The analysis was implemented using the NetworkNullHPC script (https://github.com/lentendu/NetworkNullHPC). Since our study included data from three independent sequencing runs, a normalization of the OTU-to-sample matrices was done to ensure comparability between the data. Drawing conclusions from data of independent sequencing runs without normalization is prone to misinterpretations because of technical biases, such as excessive false-positive correlations among sequencing data of the same runs. To account for this problem, NetworkNullHPC provides several options for normalizing data from which we chose the default one. This option considers the expected sequencing depth of the datasets (by default set to half of the median read abundance determined through all samples) and then performs a normalization by scaling the sum of the read counts in each sample to the expected sequencing depth. After normalization, NetworkNullHPC conducted the null model strategy by Connor et al. [64], which consists of i) a first class of null models, where noise is added to every entry of the OTU matrix; ii) the noise-added matrix is permuted, and the distribution of similarity scores in the permuted matrix is used to set the lower bound for the threshold; iii) the threshold is applied to derive the observed network. This network is used to construct the second class of null models, Erdős-Rényi, based on the average degree, and the Chung-Lu model, based on the average degree distribution of the network. This method calculates a consensus network, which contains every pairwise interaction that is present in at least 90% of the Monte Carlo samples to define the Spearman correlation threshold [64]. Furthermore, only OTUs whose read abundances were higher than 10% of the median read abundance determined through all samples were kept for the network analyses. Low-abundant OTUs (including singleton OTUs with a read abundance of 1) were discarded at this point. In the network analyses, OTUs are represented as nodes, and a statistically significant Spearman rank correlation, calculated through the null models described above, between two OTUs is represented by an edge between the respective OTUs. The networks contain only OTUs that have a significant co-occurrence or co-exclusion with at least one other OTU. We further included soil physico-chemical properties (soil texture, exchangeable bases, pH, organic carbon, phosphorus, and aluminum) of each sample from Ritter et al. [27] as environmental parameters into the network analysis for considering significant correlations between OTUs and abiotic variables.
We first produced a network with all OTUs recorded in all samples. We then compared the patterns of co-occurrence in each habitat type and produced networks for each habitat separately using all OTUs recorded in a determined habitat. Due to possible regional variation and differences in the number of replicates among habitats, we also constructed networks for each habitat within each locality. Finally, to obtain a better understanding of significant co-occurrences within different taxonomic groups, we calculated networks using OTUs of all samples for bacteria, all eukaryotes, fungi, protists, and metazoans separately. We then used the tidyverse and igraph v.1.2.4.1 [66] R packages to combine the resulting classification matrices to visualize and explore the networks with the interactive platform Gephi v.0.9.2, using the Yifan Hu layout [67].
For comparing co-occurrence patterns within micro-organism communities of the different habitats, we analyzed properties of the network structure by applying several measures of graph theory [68] (for a schematic representation of some network properties see Fig. S1): 1) the maximum component of the network, which is the subgraph that comprises the largest number of nodes (OTUs); 2) the diameter, which is the number of edges connected in the maximal shortest path length of the network and describes the size of the maximum component; 3) the average path length of the network, where length values closer to 1 indicate a more direct association of OTUs with each other; 4) the transitivity, which is the probability of two random nodes (OTUs) being directly or indirectly connected. This measure provides an indication of the clustering in the whole network, representing the presence of tripartite relationships (nodes connected by more than one path); high transitivity means a tightly connected community and might be indicative of degradation pathways or niche filtering [69]; 5) two modularity scores that assess the community structure in the network and show if the community consists of smaller groups of highly associated OTUs that are poorly associated with other OTUs. Modules may indicate different niches and have been used to study habitat preferences [69]. For this metric, we first calculated the modularity of the network and in a second step also identified Louvain communities in the networks, which are groups of highly associated OTUs. Inside these groups we calculated the modularity (LCmodularity) [70].
In addition to those measures, we calculated the assortativity, which is the tendency of nodes that share a specific attribute to be connected within the network, for each main taxonomic group [71]. In other words, we inferred the chance of OTUs of each taxonomic group being more connected among themselves than with OTUs of other taxonomic groups. Furthermore, we also calculated the degree of each node, that is, the number of direct connections it has to other nodes (Fig. S2 and S3). Finally, we calculated the key nodes, which are articulation point nodes that have a high betweenness centrality (the extent to which a node lies on paths between other nodes), and identified the OTUs that maintain the network structure and, potentially, mediate several of the processes in the community. Key nodes may represent keystone species, and their removal can cause the ecosystem to collapse [69]. These structural properties of the network allow comparison among complex networks, making it possible to compare the networks among habitats and therefore test our hypotheses. All network analyses were calculated using the igraph R package.
We obtained a total of 39,351 OTUs, of which 9,943 (25%) were classified as bacteria, 6,750 (17%) as fungi, 5,568 (14%) as protists, 5,107 (13%) as metazoans, and 438 (1%) as chloroplastida. Many OTUs (10,443; 26%) could not be classified to any resolved taxonomic level. We decided to keep them because even though some of them could be technically compromised, after all clean steps we consider the vast majority should be true OTUs not yet represented in available databases, considering the limited reference sequences for Amazonia. The bacteria were composed mostly of generalist OTUs (~50%), in contrast to a lower proportion of generalists among eukaryotes (~25%; Table 1; Fig. 2).
The network analysis considering all samples resulted in 72,577 edges (statistically significant correlations between pairs of OTUs) involving 2,350 OTUs (6% of the total), of which 61,145 edges involved 2,340 OTUs for the co-occurrence network, and 11,432 edges involved 1,343 OTUs for the co-exclusion network. The Spearman's rank correlation threshold, calculated based on the null models, was set to 0.43 for the co-occurrence and -0.42 for the co-exclusion networks. The number of OTUs, nodes, edges, and the Spearman's rank correlation threshold for habitat-specific networks are shown in Table S1.
Overall, our Amazonian network comprised highly connected OTUs. The maximum component comprised 2,337 OTUs (co-occurrence = 2,327, co-exclusion = 1,339) with a diameter (the maximal shortest path) of eight edges (both co-occurrence and co-exclusion). The transitivity – the measure of the clustering in the network – was 0.386 (co-occurrence = 0.432, co-exclusion = 0). The modularity score (range of -1 to 1) was low, -0.000002 (co-occurrence = -0.000008, co-exclusion = -0.00003), and the LC modularity was 0.469 (co-occurrence = 0.574, co-exclusion = 0.444).
The assortativity of the network showed that bacterial OTUs are mainly connected with other bacterial OTUs; similarly, eukaryotic OTUs are mainly connected with each other (Table S2). When we divided the bacteria by phyla and eukaryotes into fungi and protist realms, and into the metazoan kingdom, this tendency was still observed with few exceptions, that is, OTUs of one phylum or group will be more connected with OTUs of the same phyla (Table S2). Due to the limited number of OTUs classified as Chloroplastida, we did not calculate assortativity for this group.
For the degree attribute, the co-occurrence network showed that several bacterial groups were highly connected and served as hubs (Fig. S2A). The highest degree in the co-exclusion network, that is the nodes that were most connected (with more edges) with other nodes, was observed for several soil properties (pH, sodium [Na], potassium [K], phosphorus (P), calcium [Ca], magnesium [Mg], the sum of all exchangeable bases [SB], exchangeable aluminum [Al and H+Al], aluminium saturation index [m], Base Saturation Index [V], the effective cation exchange capacity [t], the cation exchange capacity [T], and soil texture), and OTUs classified as Acidobacteria, Actinobacteria, and Proteobacteria (Fig. S2B). On top of having a high degree, several soil properties (base saturation index, cation exchange capacity, exchangeable bases, exchangeable aluminum, aluminum saturation index, pH, and soil texture) were also identified as key nodes (Table S3) in the co-exclusion network (Fig. S2B). The values of each score by habitat are provided in Table 2. We found that a total of 124 OTUs were classified as key nodes, mostly bacteria (61%, Table S3). For networks for each taxonomic group, see Figures S4 – S7.
Although bacteria accounted for approximately half of the number of OTUs compared to eukaryotes in our data, they dominated the networks, stressing the importance of prokaryotes in soil communities. The networks are composed mostly of Proteobacteria (Alphaproteobacteria) and Acidobacteria (Fig. 3A). Other groups such as Chloroplastida, fungi, protists, and metazoans are also represented. Most of the OTUs classified as Chloroplastida belong to Magnoliophyta, mainly from the orders Fabales, Liliopsida, and Malpighiales. Most fungi are Dikarya, the majority belongs to Ascomycota and the remainder to Basidiomycota (Fig. S5). The metazoans are mostly comprised of Hexapoda (Fig. S6). Most protists belong to the supergroup SAR (mostly Rhizaria, followed by Alveolata, and Stramenopiles) followed by Amebozoa (Fig. S7). The co-exclusion network is smaller but still features all taxonomic groups (Fig. 3B). For networks scaled to indicate the degree of each node, see Figure S2.
Our first hypothesis (H1) - that most OTUs in the network will be associated with organic decomposition - was supported. Most bacteria were classified as chemoheterotrophic. For bacterial metabolism, there is a combination of respiration and fermentation, which are complementary in the context of organic decomposition. Also, most bacteria were classified as acidophilic or acid tolerant (Fig. 3 and S8). Most fungi were saprotrophic, and most protists were bacterivores (Fig. 3).
Our second hypothesis (H2) - that some soil properties (mainly pH and organic carbon) will be important factors in explaining the co-occurrence networks - was partially supported. Soil properties were important for the co-exclusion network and less important for the co-occurrence network (Fig. S2, Table S3). For the network considering all OTUs, only organic carbon (C) and base saturation index (v) were key nodes in explaining co-occurrence, whereas several other soil properties were key nodes for the co-exclusion, notably soil texture, exchangeable bases, aluminum, and pH.
Our third hypothesis (H3) - that habitat networks will show different network structures related to the presence of flood pulse as the main factor in explaining the network properties - was also partially supported. Campinas and terra-firme networks were more similar in structure and properties than igapós and várzeas (Fig. 4, Table 2). Várzeas showed by far the least complex co-occurrence network structure, but it was the habitat with the lowest number of replicates (12 samples, Fig. 4, Table 2). For the habitat networks within each locality, in general, there was a high modularity and transitivity, and a low average path length, diameter, and maximum components, with some regional variation (Table 2).
Due to the high diversity and turnover in Amazonia we could expected a neutral community assembly due to dispersion limitation [72, 73] and a high redundancy of functional traits in the microbial community [74], which would result in less connected co-occurrence networks. However, here we show that Amazonian micro-organism communities form highly interconnected co-occurrence and co-exclusion networks when compared to those of other regions [75, 12]. For instance, global scale samples of the Tara Ocean project, with 313 samples and more than 1 million OTUs, including viruses, showed a network composed of 29,912 OTUs and 127,000 edges with an average degree (number of edges per node) of 4.26 [76]. Our results, with a smaller sampling of 39 Amazonian plots (78 samples, 25% of the number of Tara Ocean samples) still resulted in a network of 72,577 edges involving 2,350 OTUs, with an average degree of 30.88. In a study surveying soil bacteria in France, which included 30 soil samples for different land use categories (forest, grassland, crop system, and vineyards), researchers found a range from 16,430 edges involving 1,083 nodes in forests to 2,046 edges involving 1,342 nodes in vineyards [12]. Our study comprised 24 samples for terra-firme, 24 for igapós, 18 for campinas, and 12 for várzeas, and ranged from 6,223 edges involving 1,047 nodes in terrra-firme to 321 edges involving 273 nodes in várzeas. Although these studies use different methods to calculate the networks, our results demonstrate the relatively high complexity of Amazonian networks and highlight the importance of biotic and/or abiotic interactions in this region.
H1: The Amazonian-wide network will be mostly composed by organisms associated with organic decomposition: As expected in our H1, considering bacterial metabolism, there is a combination of bacteria that rely on respiration (e.g., Acetobacteraceae and Planctomycetaceae) and those that rely on fermentation (e.g., Fusobacteriales and Lactobacillales). These two groups complement each other in ecological functionality and indicate a high turnover of organic material underpinned by both metabolisms. Some groups present in the networks are active in degrading cellulose (e.g., Acidothermaceae and Polyangiaceae), which is in agreement with our litter samples that had a high amount of plant debris. There are also a few nitrogen-fixating groups, but they do not seem to be very common. Yet, one third of the bacteria recorded in our networks are not assigned to any functional group or preference, stressing the need for further investigation of Amazonian bacterial communities. Furthermore, the majority of fungi found here were saprotrophic, that is, decomposers of organic matter and, therefore, important agents in carbon cycling [77]. Unfortunately, archaea were heavily underrepresented in the results (just 131 OTUs) of our sequencing approach. We do not consider this a biologically meaningful result but rather an artefact of the amplification process, since archaea are commonly found in soil samples, where they contribute to, e.g., ammonia oxidation [78, 79]. We therefore decided to exclude archaea from our analyses and refrain from making any statements about this group of organisms.
Deforestation alters patterns of co-occurrence, impacting ecosystem functions [10]. The edges of Amazonian forest fragments have an increase of soil organic carbon probably due to high tree mortality [80], but potentially also due to an alteration of the microbial community [81], which may have an impact on the rate of organic decomposition and the ecosystem resilience capability. With the increasing anthropogenic pressures in Amazonia [82], it is crucial to understand how biological interactions are linked to decomposition in Amazonia and the potential effects of microbial community alteration in the world’s most diverse forest.
H2: Environmental soil properties, especially pH and organic carbon (Ritter et al. 2018), will be the most important factors to explain co-occurrence in all networks: Our hypothesis H2 was partially supported. The organic carbon and base saturation indices were key nodes in explaining co-occurrence in the general network, whereas several other soil properties were key nodes for the co-exclusion, such as soil texture, exchangeable bases, aluminum, and pH. The composition of micro-organisms in the soil is usually determined by physico-chemical conditions such as phosphorous [83, 84, 85], pH [86, 87], and soil moisture [88]. Using the same data, Ritter et al. [27] showed that pH was the most important factor explaining the community turnover, and organic carbon the OTU richness [27]. This is the reason why we expected the pH to be important for the co-occurrence network, but instead we found it to be important just for co-exclusion.
Our results showed a predominance of bacteria that grow in low pH regimes, especially Acidobacteria [89], which are very abundant in soil. Amazonian soils are acidic, with pH varying from 3.65 to 5.14 [27], which explains the large presence of acidophilic organisms in both co-occurrence and co-exclusion networks. However, in less acidic soils, other groups could co-occur with acidophilic organisms, but the acidity tends to lead to increased competition. Phyla associated with low pH, such as Acidobacteria and Proteobacteria [85, 90], may be better competitors in more acidic soils. Furthermore, temperature appears to select for groups that are tolerant to high temperatures up to 40 °C (Fig. S8), which can occur in Amazonia all year. The soil properties were less important in the habitat-specific networks. As these soil properties are associated with the habitat types, and the same habitat presents similar soil physico-chemical conditions [27], the network analysis of habitats may already be controlled by the variation in soil properties.
H3: The habitat-specific networks will show different network structures, with the presence of flood pulse (the periodic inundation of floodplains along certain rivers) as the main factor explaining network properties: Our hypotheses H3 was also partially supported. We expected that the high environmental stress in igapós and várzeas due to the seasonal flooding would act as an abiotic filter, keeping the same set of specialized organisms that co-occur for long periods of time, which would result in a dense and highly connected network. However, both habitat-specific networks were sparsely connected considering analyses within localities (Fig. S9-S12). Both habitats remain submerged during most of the year, up to 240 days [91]. In our data, the seasonally flooded forests – the várzeas and igapós – were more similar to each other in terms of community composition [36], which could be related to similar environmental filters linked to stress by flooding [92, 93], although the transitivity, that may indicate niche filtering, of the habitat networks was similar. Other factors, such as the random colonization due to smaller area size than in terra-firme, fragmented habitat distribution [94, 95], or the seasonally introduced species with the flood pulse that brings organism from river curse, could to some extent randomize the presence of OTUs. Várzeas differ from igapós by being a more fertile habitat, as their waters carry sediments from the Andes [92]. Therefore, a higher rate of colonization by micro-organisms is to be expected in várzeas than in igapós, which are bathed by acidic, low-fertile waters [96, 97]. Várzeas may therefore have more micro-organisms that could survive this particular stress condition, not just specialists. Indeed, várzeas had fewer specialist OTUs than igapós [36]. The difference between the sampling size of igapós (24 replicates) and that of várzeas (12 replicates), though, makes a more general comparison between these habitats difficult.
Communities in campinas (Fig. 4A) resulted in a smaller network (fewer nodes and fewer edges) than terra-firme (Fig. 4C) and igapós (Fig. 4E, Table 2), but larger than várzeas (Fig. 4G). However, terra-firme and igapós had more replicates (24 each) than campinas (18 replicates). Regarding network properties, campinas and terra-firme are similar (Table 2). However, the habitat network of campinas shows two distinct modules (Fig. 4A) that either represent two different metabolic pathways, as the dominant bacterial phyla are different in each module, or indicate that there is a geographical clustering [69]. The geographical clustering could be related to the natural fragmentation of campina habitats, which makes colonization more random than in more continuous habitats such as terra-firme [94, 95]. Yet, comparing the campinas networks within localities, the localities bathed by acidic, sediment-poor black waters (Jaú and Cuieras) had more modular networks (Fig. S9-S12, Table2). That is also true for terra-firme, that had the low modularity (0.47 in co-occurrence network; Table 2), in other words a more connected network in Caxiuanã, a locality bathed by rich sediments from a white-water river (Fig. S9-S12). Stressful localities may have a more specialized co-occurrence community due to the environmental filters such as low pH, while less stressful localities may allow more organisms to co-occur.
General networks composition: Overall, the key nodes of our networks are mainly affiliated to Proteobacteria (29 nodes, of which 17 belong to Alphaproteobacteria), Acidobacteria (19), Planctomycetes (12), and Actinobacteria (10, Table S3). These are common groups of bacteria in soils in general [52] and also in our samples [36, 37], being the most frequent groups in our networks. Alphaproteobacteria are a highly diverse clade of Gram-negative bacteria with several biological functions, including metabolizing C1 compounds [98], fixating nitrogen [99], endosymbiosis such as the widespread and important genus Wolbachia [100], and also intracellular pathogenicity [101]. Actinobacteria are important for the decomposition of organic matter from soils [102], making a range of nutrients available to other organisms, which probably explains the importance of this group in the co-occurrence network. Planctomycetes, a phylum of mostly aquatic free-living bacteria, but also found in soils [103], were also rich in our samples, with their distribution linked to soil properties such as soil organic matter, Ca2+, and pH [103]. Interestingly, even though Planctomycetes are associated with aquatic environments, they were not the richest in seasonally flooded forests (Fig. 4).
The co-occurrence network for all samples is densely connected, with low modularity and without differentiation explained by the habitat types, since most OTUs were classified as habitat generalists. If habitat type was the strongest variable to explain co-occurrence, the network should present modules, since groups belonging to a particular module should have similar environmental preferences [69, 104]. For instance, through analysis of deforestation along an Amazonian transect, it was possible to differentiate modules more associated with different degrees of deforestation [10]. The lack of modularity in the network for all samples may be due to the majority of OTUs (> 98% of nodes in all our networks) being classified as habitat generalists (present in three or four habitat types). The specialist OTUs (present in just one habitat) were usually local specialists, i.e., recorded in one or very few plots (Fig. 2), making the co-occurrence and co-exclusion patterns for specialist organisms hard to detect. Finally, even though networks for habitats are different from each other, the resulting pattern when including all samples may not show any apparent modularity, as patterns for each habitat are overlaid, and OTUs may present additional connections with OTUs from other habitats, smoothing or hiding patterns from individual habitats.
The high number of local specialist OTUs is in agreement with general microbial patterns elsewhere [105, 106] and also with distribution patterns of Amazonian tree species [95]. It is also in agreement with the competitive exclusion theory, which postulates that species with low competitive skills should be excluded from the community by highly competitive species [107]. Yet, these “poor” competitors remain in the communities as rare species [108]. In a global-scale study, rare species were found to have the highest positive association with each other, whereas common species had more negative associations [11]. However, that study relied on the number of observations of macro-organism species (plants and animals) and defined rare species in terms of low abundance. Abundance is difficult to quantify in metabarcoding studies due to PCR biases, such as the false negatives [109]. These biases will affect the detection of OTUs with low abundance beyond the intrinsic low detection probability and the stochastic distribution in their habitats due to their small abundance [110]. Yet, we had a relatively highly populated co-exclusion networks, which may be related to a high competition in soil communities, and further supports the exclusion theory in Amazonian microbial soil communities.
As microbial metabolic functions can differ or overlap, indicating metabolic complementarity or redundancy, they can explain either co-occurrence due to mutualism or co-exclusion due to competition [69]. For instance, several bacterial phyla are syntrophic, with complementary ecosystem functionality, such as microbial interactions for anaerobic methane oxidation, thermodynamic degradation, and nutritional exchange [111]. These associations potentially reinforce co-occurrence patterns among these groups and can explain the bacterial dominance in the co-occurrence networks. On the other hand, bacteria are also dominant in the co-exclusion networks, probably due to competition between other bacteria and also with fungi and protists. Examples of antagonist relationships between bacteria and fungi are well established [112, 113]. Most protists recovered in the networks feed on bacteria and some on fungi (or both), with few eukaryotic parasites recorded. Although this pattern is very different from what was found in other Neotropical soils [43], our results can be explained due to the prevalence of bacteria in our networks.
We have shown that Amazonian soil communities have a densely connected network when compared with temperate forests or oceans. Bacteria dominate both the co-occurrence and co-exclusion networks. As some bacteria interact in ways that affect their growth and metabolism, these cooperative metabolic interactions can lead to increased growth of interacting bacteria and ultimately to positive co-occurrence patterns. On the other hand, competition for the same resources may lead to an inverse pattern and explain the dominance of bacteria in co-exclusion networks. Co-occurrence patterns may also reflect the response of different species to a common environmental factor rather than their direct interactions. All habitat networks in Amazonia were similar in structure and properties. Although co-occurrence networks cannot map out direct interactions in complex microbial communities, microbial network studies provide a way forward to understand potential intertaxon relationships in microbial communities.
Competing Interests
The authors declare no competing interests.
Funding
We acknowledge funding from the Alexander von Humboldt Foundation to C.D.R., and the Deutsche Forschungsgemeinschaft (#DU1319/5-1) to M.D. the Swedish Research Council, Swedish Foundation for Strategic Research, the Knut and Alice Wallenberg Foundation, and the Royal Botanic Gardens, Kew to A.A.
Author contributions
C.D.R., D.F, and M.D. designed the study; C.D.R. and D.F. led all analyses; M.E.T. reviewed the classification of the bacteria functional groups. M.D. reviewed the classification of the protists functional groups. C.D.R. led the writing of the manuscript with contributions from D.M and J.A. A.A., R.H.N., M.E.T., and M.D. reviewed the manuscript. All authors approved the submission.
Data Availability
All data is available in Ritter et al. 2019a & 2019b.
Ethics Approval
No ethical approval was necessary to conduct this research.
- Myers JA, Harms KE (2009) Seed arrival, ecological filters, and plant species richness: a meta‐analysis. Ecol Lett 12:1250–1260
- Stegen JC, Lin X, Konopka AE, Fredrickson JK (2012) Stochastic and deterministic assembly processes in subsurface microbial communities. ISME J 6:1653
- Blanchet FG, Cazelles K, Gravel D (2020) Co‐occurrence is not evidence of ecological interactions. Ecol Lett 23:1050–1063. https://doi.org/10.1111/ele.13525
- Hutchinson GE (1978) An introduction to population ecology. Hew Haven Connecticut Yale University Press 1978., New Haven
- Williams RJ, Howe A, Hofmockel KS (2014) Demonstrating microbial co-occurrence pattern analyses within and between ecosystems. Front Microbiol 5:1–10. https://doi.org/10.3389/fmicb.2014.00358
- Yang P, Yu S, Cheng L, Ning K (2019) Meta-network: optimized species-species network analysis for microbial communities. BMC Genomics 20:187
- Geary WL, Nimmo DG, Doherty TS, et al (2019) Threat webs: Reframing the co‐occurrence and interactions of threats to biodiversity. J Appl Ecol 56:1992–1997
- Swenson NG, Jones FA (2017) Community transcriptomics, genomics and the problem of species co‐occurrence. J Ecol 105:563–568
- Wood SA, Gilbert JA, Leff JW, et al (2017) Consequences of tropical forest conversion to oil palm on soil bacterial community and network structure. Soil Biol Biochem 112:258–268
- Khan MAW, Bohannan BJM, Nüsslein K, et al (2019) Deforestation impacts network co-occurrence patterns of microbial communities in Amazon soils. FEMS Microbiol Ecol 95:fiy230
- Calatayud J, Andivia E, Escudero A, et al (2019) Positive associations among rare species and their persistence in ecological assemblages. Nat Ecol Evol 1–6
- Karimi B, Dequiedt S, Terrat S, et al (2019) Biogeography of soil bacterial networks along a gradient of cropping intensity. Sci Rep 9:1–10
- Brearley FQ, Elliott DR, Iribar A, Sen R (2016) Arbuscular mycorrhizal community structure on co-existing tropical legume trees in French Guiana. Plant Soil 403:253–265
- Demenois J, Carriconde F, Bonaventure P, et al (2018) Impact of plant root functional traits and associated mycorrhizas on the aggregate stability of a tropical Ferralsol. Geoderma 312:6–16
- Pagano MC, Falcão N, Weber OB, et al (2019) Mycorrhizas in South American Anthropic Environments. In: Mycorrhizal Fungi in South America. Springer, pp 343–365
- Escribano-Avila G, Lara-Romero C, Heleno R, Traveset A (2018) Tropical seed dispersal networks: emerging patterns, biases, and keystone species traits. In: Ecological networks in the tropics. Springer, pp 93–110
- Machado-de-Souza T, Campos RP, Devoto M, Varassin IG (2019) Local drivers of the structure of a tropical bird-seed dispersal network. Oecologia 189:421–433
- Chomicki G, Weber M, Antonelli A, et al (2019) The Impact of Mutualisms on Species Richness. Trends Ecol Evol 1:698–711
- Staab M, Bruelheide H, Durka W, et al (2016) Tree phylogenetic diversity promotes host–parasitoid interactions. Proc R Soc B Biol Sci 283:20160275
- Galiana N, Hawkins BA, Montoya JM (2019) The geographical variation of network structure is scale dependent: understanding the biotic specialization of host–parasitoid networks. Ecography (Cop) 42:1175–1187
- Doliwa A, Dunthorn M, Rassoshanska E, et al (2020) Identifying potential hosts of short-branch Microsporidia. Microb Ecol 1–5
- Mora C, Tittensor DP, Adl S, et al (2011) How many species are there on earth and in the ocean? PLoS Biol 9:1–8. https://doi.org/10.1371/journal.pbio.1001127
- Dominati E, Patterson M, Mackay A (2010) A framework for classifying and quantifying the natural capital and ecosystem services of soils. Ecol Econ 69:1858–1868
- Cavicchioli R, Ripple WJ, Timmis KN, et al (2019) Scientists’ warning to humanity: microorganisms and climate change. Nat Rev Microbiol 17:569–586
- Krishna MP, Mohan M (2017) Litter decomposition in forest ecosystems: a review. Energy, Ecol Environ 2:236–249
- Four B, Cárdenas RE, Dangles O (2019) Traits or habitat? Disentangling predictors of leaf‐litter decomposition in Amazonian soils and streams. Ecosphere 10:e02691
- Ritter CD, Zizka A, Roger F, et al (2018) High-throughput metabarcoding reveals the effect of physicochemical soil properties on soil and litter biodiversity and community turnover across Amazonia. PeerJ 2018:e5661. https://doi.org/10.7717/peerj.5661
- Maretti CC (2014) Amazon: There is Hope! If we all do ‘the right thing’…; Deforestation, Protected Areas and Indigenous Territories: Past, evolution and… Which future?
- Dirzo R, Young HS, Galetti M, et al (2014) Defaunation in the Anthropocene. Science (80- ) 345:401–406
- Wardle DA, Bardgett RD, Klironomos JN, et al (2004) Ecological linkages between aboveground and belowground biota. Science (80- ) 304:1629–1633
- Hansen MC, Potapov P V, Moore R, et al (2013) High-resolution global maps of 21st-century forest cover change. Science (80- ) 342:850–853
- Ojea E, Martin-Ortega J, Chiabai A (2012) Defining and classifying ecosystem services for economic valuation: the case of forest water services. Environ Sci Policy 19:1–15
- Fearnside PM (2008) Amazon forest maintenance as a source of environmental services. An Acad Bras Cienc 80:101–114
- Ribas CC, Aleixo A, Nogueira ACR, et al (2012) A palaeobiogeographic model for biotic diversification within Amazonia over the past three million years. Proc R Soc B Biol Sci 279:681–689
- Cohn-Haft M, Naka LN, Fernandes AM (2007) Padrões de distribuição da avifauna da várzea dos rios Solimões e Amazonas. Conserv da várzea identificação e Caracter regiões biogeográficas Manaus, Ibama, ProVárzea, 356p 287–323
- Ritter CD, Zizka A, Barnes C, et al (2019) Locality or habitat? Exploring predictors of biodiversity in Amazonia. Ecography (Cop) 42:321–333. https://doi.org/10.1111/ecog.03833
- Ritter CD, Faurby S, Bennett DJ, et al (2019) The pitfalls of biodiversity proxies: Differences in richness patterns of birds, trees and understudied diversity across Amazonia. Sci Rep 9:1–13
- Adeney JM, Christensen NL, Vicentini A, Cohn‐Haft M (2016) White‐sand ecosystems in Amazonia. Biotropica 48:7–23
- Capurucho JMG, Borges SH, Cornelius C, et al (2020) Patterns and processes of diversification in Amazonian white sand ecosystems: insights from birds and plants. In: Neotropical Diversification: Patterns and Processes. Springer, pp 245–270
- Junk WJ, Wittmann F, Schöngart J, et al (2020) Large rivers and their floodplains: Structures, functions, evolutionary traits and management with special reference to the Brazilian rivers. Introd Large Rivers 55–101
- Dunthorn M, Kauserud H, Bass D, et al (2017) Yeasts dominate soil fungal communities in three lowland Neotropical rainforests. Environ Microbiol Rep 9:668–675. https://doi.org/10.1111/1758-2229.12575
- Kim J-S, Sparovek G, Longo RM, et al (2007) Bacterial diversity of terra preta and pristine forest soil from the Western Amazon. Soil Biol Biochem 39:684–690
- Mahé F, De Vargas C, Bass D, et al (2017) Parasites dominate hyperdiverse soil protist communities in Neotropical rainforests. Nat Ecol Evol 1:1–8. https://doi.org/10.1038/s41559-017-0091
- Rodrigues JLM, Pellizari VH, Mueller R, et al (2013) Conversion of the Amazon rainforest to agriculture results in biotic homogenization of soil bacterial communities. Proc Natl Acad Sci 110:988–993
- Zinger L, Taberlet P, Schimann H, et al (2019) Body size determines soil community assembly in a tropical forest. Mol Ecol 28:528–543. https://doi.org/10.1111/mec.14919
- Forster D, Bittner L, Karkar S, et al (2015) Testing ecological theories with sequence similarity networks: marine ciliates exhibit similar geographic dispersal patterns as multicellular organisms. BMC Biol 13:16
- Araújo MB, Rozenfeld A, Rahbek C, Marquet PA (2011) Using species co‐occurrence networks to assess the impacts of climate change. Ecography (Cop) 34:897–908
- Gao J, Barzel B, Barabási A-L (2016) Universal resilience patterns in complex networks. Nature 530:307–312
- de Araujo ASF, Mendes LW, Lemos LN, et al (2018) Protist species richness and soil microbiome complexity increase towards climax vegetation in the Brazilian Cerrado. Commun Biol 1:135. https://doi.org/10.1038/s42003-018-0129-0
- Delmas E, Besson M, Brice M, et al (2019) Analysing ecological networks of species interactions. Biol Rev 94:16–36
- Blaxter M, Mann J, Chapman T, et al (2005) Defining operational taxonomic units using DNA barcode data. 1935–1943. https://doi.org/10.1098/rstb.2005.1725
- Delgado-Baquerizo M, Oliverio AM, Brewer TE, et al (2018) A global atlas of the dominant bacteria found in soil. Science (80- ) 359:320–325. https://doi.org/10.1126/science.aap9516
- Klindworth A, Pruesse E, Schweer T, et al (2013) Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res 41:e1–e1
- Guardiola M, Uriz MJ, Taberlet P, et al (2015) Deep-sea, deep-sequencing: metabarcoding extracellular DNA from sediments of marine canyons. PLoS One 10:e0139633
- Leray M, Yang JY, Meyer CP, et al (2013) A new versatile primer set targeting a short fragment of the mitochondrial COI region for metabarcoding metazoan diversity: application for characterizing coral reef fish gut contents. Front Zool 10:34
- Meyer CP (2003) Molecular systematics of cowries (Gastropoda: Cypraeidae) and diversification patterns in the tropics. Biol J Linn Soc 79:401–459
- Edgar RC (2013) UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods 10:996–998
- Wickham H (2017) tidyverse: Easily Install and Load “Tidyverse” Packages (Version R package version 1.1. 1)
- Wickham H (2016) ggplot2: elegant graphics for data analysis. Springer
- R Core Team (2003) The R development core team. R A Lang. Environ. Stat. Comput. 1
- Louca S, Parfrey LW, Doebeli M (2016) Decoupling function and taxonomy in the global ocean microbiome. Science (80- ) 353:1272–1277
- Adl SM, Bass D, Lane CE, et al (2019) Revisions to the Classification, Nomenclature, and Diversity of Eukaryotes. J Eukaryot Microbiol 66:4–119. https://doi.org/10.1111/jeu.12691
- Ritter CD, Dunthorn M, Anslan S, et al (2020) Advancing biodiversity assessments with environmental DNA : Long-read technologies help reveal the drivers of Amazonian fungal diversity. 00:1–16. https://doi.org/10.1002/ece3.6477
- Connor N, Barberán A, Clauset A (2017) Using null models to infer microbial cooccurrence networks. PLoS One 12:1–23. https://doi.org/10.1371/journal.pone.0176751
- Carr A, Diener C, Baliga NS, Gibbons SM (2019) Use and abuse of correlation analyses in microbial ecology. ISME J 13:2647–2655
- Csardi G, Nepusz T (2006) The igraph software package for complex network research. InterJournal, Complex Syst 1695:1–9
- Bastian M, Heymann S, Jacomy M (2009) Gephi: an open source software for exploring and manipulating networks. In: Third international AAAI conference on weblogs and social media
- Karimi B, Maron PA, Boure NC-P, et al (2017) Microbial diversity and ecological networks as indicators of environmental quality. Environ Chem Lett 15:265–281
- Röttjers L, Faust K (2018) From hairballs to hypotheses–biological insights from microbial networks. FEMS Microbiol Rev 42:761–780
- Blondel VD, Guillaume J-L, Lambiotte R, Lefebvre E (2008) Fast unfolding of communities in large networks. J Stat Mech theory Exp 2008:P10008
- Kurtz ZD, Müller CL, Miraldi ER, et al (2015) Sparse and compositionally robust inference of microbial ecological networks. PLoS Comput Biol 11:e1004226
- Hubbell SP (2005) Neutral theory in community ecology and the hypothesis of functional equivalence. Funct Ecol 19:166–172
- Dambros C, Zuquim G, Moulatlet GM, et al (2020) The role of environmental filtering, geographic distance and dispersal barriers in shaping the turnover of plant and animal species in Amazonia. Biodivers Conserv 29:3609–3634
- Louca S, Polz MF, Mazel F, et al (2018) Function and functional redundancy in microbial systems. Nat Ecol Evol 2:936–943
- Barberán A, Bates ST, Casamayor EO, Fierer N (2012) Using network analysis to explore co-occurrence patterns in soil microbial communities. ISME J 6:343–351. https://doi.org/10.1038/ismej.2011.119
- Lima-Mendez G, Faust K, Henry N, et al (2015) Determinants of community structure in the global plankton interactome. Science (80- ) 348:1262073
- Swift MJ (1982) Basidiomycetes as components of forest ecosystems. In: Decomposer basidiomycetes: their biology and ecology. Cambridge University Press
- Leininger S, Urich T, Schloter M, et al (2006) Archaea predominate among ammonia-oxidizing prokaryotes in soils. Nature 442:806–809
- Bates ST, Berg-Lyons D, Caporaso JG, et al (2011) Examining the global distribution of dominant archaeal populations in soil. ISME J 5:908–917
- Barros HS, Fearnside PM (2016) Soil carbon stock changes due to edge effects in central Amazon forest fragments. For Ecol Manage 379:30–36
- Kroeger ME, Delmont TO, Eren AM, et al (2018) New biological insights into how deforestation in Amazonia affects soil microbial communities using metagenomics and metagenome-assembled genomes. Front Microbiol 9:1635
- Pereira EJ de AL, Ferreira PJS, de Santana Ribeiro LC, et al (2019) Policy in Brazil (2016–2019) threaten conservation of the Amazon rainforest. Environ Sci Policy 100:8–12
- Buckley DH, Schmidt TM (2001) Environmental factors influencing the distribution of rRNA from Verrucomicrobia in soil. FEMS Microbiol Ecol 35:105–112
- Faoro H, Alves AC, Souza EM, et al (2010) Influence of soil characteristics on the diversity of bacteria in the Southern Brazilian Atlantic Forest. Appl Environ Microbiol 76:4744–4749
- Navarrete AA, Kuramae EE, de Hollander M, et al (2013) Acidobacterial community responses to agricultural management of soybean in Amazon forest soils. FEMS Microbiol Ecol 83:607–621. https://doi.org/10.1111/1574-6941.12018
- Rousk J, Brookes PC, Bååth E (2010) Investigating the mechanisms for the opposing pH relationships of fungal and bacterial growth in soil. Soil Biol Biochem 42:926–934
- Wang JT, Zheng YM, Hu HW, et al (2015) Soil pH determines the alpha diversity but not beta diversity of soil fungal community along altitude in a typical Tibetan forest ecosystem. J Soils Sediments 15:1224–1232. https://doi.org/10.1007/s11368-015-1070-1
- Bates ST, Clemente JC, Flores GE, et al (2013) Global biogeography of highly diverse protistan communities in soil. ISME J 7:652–659. https://doi.org/10.1038/ismej.2012.147
- Barns SM, Takala SL, Kuske CR (1999) Wide distribution and diversity of members of the bacterial kingdom Acidobacterium in the environment. Appl Environ Microbiol 65:1731–1737
- Dai Z, Su W, Chen H, et al (2018) Long‐term nitrogen fertilization decreases bacterial diversity and favors the growth of Actinobacteria and Proteobacteria in agro‐ecosystems across the globe. Glob Chang Biol 24:3452–3461
- Junk WJ (1997) General aspects of floodplain ecology with special reference to Amazonian floodplains. In: The Central Amazon Floodplain. Springer, pp 3–20
- Myster RW (2016) The physical structure of forests in the Amazon Basin: a review. Bot Rev 82:407–427
- Parolin P de, De Simone O, Haase K, et al (2004) Central Amazonian floodplain forests: tree adaptations in a pulsing system. Bot Rev 70:357–380
- ter Steege H, Sabatier D, Castellanos H, et al (2000) An analysis of the floristic composition and diversity of Amazonian forests including those of the Guiana Shield. J Trop Ecol 16:801–828
- ter Steege H, Henkel TW, Helal N, et al (2019) Rarity of monodominance in hyperdiverse Amazonian forests. Sci Rep 9:1–15
- Scudeller VV, Vegas‐Vilarrúbia T (2018) Distribution and β‐diversity of tree species in igapó forests (Negro River basin, Brazilian Amazon). J Veg Sci 29:1052–1064
- Singer R, Aguiar IA (1986) Litter decomposing and ectomycorrhizal Basidiomycetes in an igapó forest. Plant Syst Evol 153:107–117
- Green PN (2006) Methylobacterium. Prokaryotes Vol 5 Proteobacteria Alpha Beta Subclasses 257–265
- Fisher RF, Long SR (1992) Rhizobium–plant signal exchange. Nature 357:655
- Hedges LM, Brownlie JC, O’Neill SL, Johnson KN (2008) Wolbachia and virus protection in insects. Science (80- ) 322:702
- Andersson SGE, Zomorodipour A, Andersson JO, et al (1998) The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature 396:133
- Eisenlord SD, Zak DR (2010) Simulated atmospheric nitrogen deposition alters actinobacterial community composition in forest soils. Soil Sci Soc Am J 74:1157–1166
- Buckley DH, Huangyutitham V, Nelson TA, et al (2006) Diversity of Planctomycetes in soil in relation to soil history and environmental heterogeneity. Appl Environ Microbiol 72:4522–4531
- Olesen JM, Bascompte J, Dupont YL, Jordano P (2007) The modularity of pollination networks. Proc Natl Acad Sci 104:19891–19896
- Ashby MN, Rine J, Mongodin EF, et al (2007) Serial analysis of rRNA genes and the unexpected dominance of rare members of microbial communities. Appl Environ Microbiol 73:4532–4542
- Nemergut DR, Costello EK, Hamady M, et al (2011) Global patterns in the biogeography of bacterial taxa. Environ Microbiol 13:135–144
- Hardin G (1960) The competitive exclusion principle. Science (80- ) 131:1292–1297
- Hutchinson GE (1961) The paradox of the plankton. Am Nat 95:137–145
- Santoferrara LF (2019) Current practice in plankton metabarcoding: optimization and error management. J Plankton Res 41:571–582
- MacArthur RH (1984) Geographical ecology: patterns in the distribution of species. Princeton University Press
- Morris BEL, Henneberger R, Huber H, Moissl-Eichinger C (2013) Microbial syntrophy: interaction for the common good. FEMS Microbiol Rev 37:384–406
- Mille‐Lindblom C, Fischer H, J. Tranvik L (2006) Antagonism between bacteria and fungi: substrate competition and a possible tradeoff between fungal growth and tolerance towards bacteria. Oikos 113:233–242
- Rousk J, Bååth E, Brookes PC, et al (2010) Soil bacterial and fungal communities across a pH gradient in an arable soil. ISME J 4:1340
Table 1 – Number and proportion (in brackets) of OTUs and their classification by habitat (generalist ≥ 3 and specialist = 1). More bacterial OTUs were classified as generalists (~50%) than was the case for eukaryotes (~25%), while eukaryotes OTUs were more specialist (~50%) than prokaryotes (~25%).
|
Group |
Total |
Generalist |
Specialist |
|
All |
39351 |
10875 (28%) |
20495 (52%) |
|
Bacteria |
9943 |
5007 (50%) |
2761 (28%) |
|
Protists |
6568 |
1659 (25%) |
3395 (52%) |
|
Fungi |
6750 |
1394 (21%) |
3843 (57%) |
|
Metazoa |
5107 |
1242 (24%) |
2763 (54%) |
|
Chloroplastida |
438 |
107 (24%) |
231 (53%) |
|
Unknown |
10443 |
1422 (14%) |
7477 (72%) |
Table 2 – Properties of habitat-specific networks for all samples in the same habitat (all) and for habitat within each locality (in west-to-east order: BC = Benjamin Constant; JAU = Jaú; CUI = Cuieras; and CXN = Caxiuanã).
