Dynamics of bacterial community succession in a salt marsh chronosequence: evidences for temporal niche partitioning

The mechanisms underlying community assembly and promoting temporal succession are often overlooked in microbial ecology. Here, we studied an undisturbed salt marsh chronosequence, spanning over a century of ecosystem development, to understand bacterial succession in soil. We used 16S rRNA gene-based quantitative PCR to determine bacterial abundance and multitag 454 pyrosequencing for community composition and diversity analyses. Despite 10-fold lower 16S rRNA gene abundances, the initial stages of soil development held higher phylogenetic diversities than the soil at late succession. Temporal variations in phylogenetic β-diversity were greater at initial stages of soil development, possibly as a result of the great dynamism imposed by the daily influence of the tide, promoting high immigration rates. Allogenic succession of bacterial communities was mostly driven by shifts in the soil physical structure, as well as variations in pH and salinity, which collectively explained 84.5% of the variation concerning community assemblage. The community assembly data for each successional stage were integrated into a network co-occurrence analysis, revealing higher complexity at initial stages, coinciding with great dynamism in turnover and environmental variability. Contrary to a spatial niche-based perspective of bacterial community assembly, we suggest temporal niche partitioning as the dominant mechanism of assembly (promoting more phylotype co-occurrence) in the initial stages of succession, where continuous environmental change results in the existence of multiple niches over short periods of time.     

Temporal variability in soil microbial communities across land-use types

Although numerous studies have investigated changes in soil microbial communities across space, questions about the temporal variability in these communities and how this variability compares across soils have received far less attention. We collected soils on a monthly basis (May to November) from replicated plots representing three land-use types (conventional and reduced-input row crop agricultural plots and early successional grasslands) maintained at a research site in Michigan, USA. Using barcoded pyrosequencing of the 16S rRNA gene, we found that the agricultural and early successional land uses harbored unique soil bacterial communities that exhibited distinct temporal patterns. α-Diversity, the numbers of taxa or lineages, was significantly influenced by the sampling month with the temporal variability in α-diversity exceeding the variability between land-use types. In contrast, differences in community composition across land-use types were reasonably constant across the 7-month period, suggesting that the time of sampling is less important when assessing β-diversity patterns. Communities in the agricultural soils were most variable over time and the changes were significantly correlated with soil moisture and temperature. Temporal shifts in bacterial community composition within the successional grassland plots were less predictable and are likely a product of complex interactions between the soil environment and the more diverse plant community. Temporal variability needs to be carefully assessed when comparing microbial diversity across soil types and the temporal patterns in microbial community structure can not necessarily be generalized across land uses, even if those soils are exposed to the same climatic conditions.     

Landscape Position Influences Microbial Composition and Function via Redistribution of Soil Water across a Watershed

Subalpine forest ecosystems influence global carbon cycling. However, little is known about the compositions of their soil microbial communities and how these may vary with soil environmental conditions. The goal of this study was to characterize the soil microbial communities in a subalpine forest watershed in central Montana (Stringer Creek Watershed within the Tenderfoot Creek Experimental Forest) and to investigate their relationships with environmental conditions and soil carbonaceous gases. As assessed by tagged Illumina sequencing of the 16S rRNA gene, community composition and structure differed significantly among three landscape positions: high upland zones (HUZ), low upland zones (LUZ), and riparian zones (RZ). Soil depth effects on phylogenetic diversity and β-diversity varied across landscape positions, being more evident in RZ than in HUZ. Mantel tests revealed significant correlations between microbial community assembly patterns and the soil environmental factors tested (water content, temperature, oxygen, and pH) and soil carbonaceous gases (carbon dioxide concentration and efflux and methane concentration). With one exception, methanogens were detected only in RZ soils. In contrast, methanotrophs were detected in all three landscape positions. Type I methanotrophs dominated RZ soils, while type II methanotrophs dominated LUZ and HUZ soils. The relative abundances of methanotroph populations correlated positively with soil water content (R = 0.72, P < 0.001) and negatively with soil oxygen (R = −0.53, P = 0.008). Our results suggest the coherence of soil microbial communities within and differences in communities between landscape positions in a subalpine forested watershed that reflect historical and contemporary environmental conditions.     

Temporal variation of β-diversity and assembly mechanisms in a bacterial metacommunity

The turnover of community composition across space, β-diversity, is influenced by different assembly mechanisms, which place varying weight on local habitat factors, such as environmental conditions and species interactions, and regional factors such as dispersal and history. Several assembly mechanisms may function simultaneously; however, little is known about how their importance changes over time and why. Here, we implemented a field survey where we sampled a bacterial metacommunity consisting of 17 rock pools located at the Swedish Baltic Sea coast at 11 occasions during 1 year. We determined to which extent communities were structured by different assembly mechanisms using variation partitioning and studied changes in β-diversity across environmental gradients over time. β-Diversity was highest at times of high overall productivity and environmental heterogeneity in the metacommunity, at least partly due to species sorting, that is, selection of taxa by the prevailing environmental conditions. In contrast, dispersal-driven assembly mechanisms were primarily detected at times when β-diversity was relatively low. There were no indications for strong and persistent differences in community composition or β-diversity between permanent and temporary pools, indicating that the physical disturbance regime is of relatively minor importance. In summary, our study clearly suggests that there are temporal differences in the relative importance of different assembly mechanisms related to abiotic factors and shows that the temporal variability of those factors is important for a more complete understanding of bacterial metacommunity dynamics.     

Niche and metabolic principles explain patterns of diversity and distribution: theory and a case study with soil bacterial communities

The causes of biodiversity patterns are controversial and elusive due to complex environmental variation, covarying changes in communities, and lack of baseline and null theories to differentiate straightforward causes from more complex mechanisms. To address these limitations, we developed general diversity theory integrating metabolic principles with niche-based community assembly. We evaluated this theory by investigating patterns in the diversity and distribution of soil bacteria taxa across four orders of magnitude variation in spatial scale on an Antarctic mountainside in low complexity, highly oligotrophic soils. Our theory predicts that lower temperatures should reduce taxon niche widths along environmental gradients due to decreasing growth rates, and the changing niche widths should lead to contrasting α- and β-diversity patterns. In accord with the predictions, α-diversity, niche widths and occupancies decreased while β-diversity increased with increasing elevation and decreasing temperature. The theory also successfully predicts a hump-shaped relationship between α-diversity and pH and a negative relationship between α-diversity and salinity. Thus, a few simple principles explained systematic microbial diversity variation along multiple gradients. Such general theory can be used to disentangle baseline effects from more complex effects of temperature and other variables on biodiversity patterns in a variety of ecosystems and organisms.     

Impact of Logging and Forest Conversion to Oil Palm Plantations on Soil Bacterial Communities in Borneo

Tropical forests are being rapidly altered by logging and cleared for agriculture. Understanding the effects of these land use changes on soil bacteria, which constitute a large proportion of total biodiversity and perform important ecosystem functions, is a major conservation frontier. Here we studied the effects of logging history and forest conversion to oil palm plantations in Sabah, Borneo, on the soil bacterial community. We used paired-end Illumina sequencing of the 16S rRNA gene, V3 region, to compare the bacterial communities in primary, once-logged, and twice-logged forest and land converted to oil palm plantations. Bacteria were grouped into operational taxonomic units (OTUs) at the 97% similarity level, and OTU richness and local-scale α-diversity showed no difference between the various forest types and oil palm plantations. Focusing on the turnover of bacteria across space, true β-diversity was higher in oil palm plantation soil than in forest soil, whereas community dissimilarity-based metrics of β-diversity were only marginally different between habitats, suggesting that at large scales, oil palm plantation soil could have higher overall γ-diversity than forest soil, driven by a slightly more heterogeneous community across space. Clearance of primary and logged forest for oil palm plantations did, however, significantly impact the composition of soil bacterial communities, reflecting in part the loss of some forest bacteria, whereas primary and logged forests did not differ in composition. Overall, our results suggest that the soil bacteria of tropical forest are to some extent resilient or resistant to logging but that the impacts of forest conversion to oil palm plantations are more severe.     

Spatial and Temporal Variation of Archaeal,Bacterial and Fungal Communities in Agricultural Soils

Background

Soil microbial communities are in constant change at many different temporal and spatial scales. However, the importance of these changes to the turnover of the soil microbial communities has been rarely studied simultaneously in space and time.

Methodology/Principal Findings

In this study, we explored the temporal and spatial responses of soil bacterial, archaeal and fungal β-diversities to abiotic parameters. Taking into account data from a 3-year sampling period, we analyzed the abundances and community structures of Archaea, Bacteria and Fungi along with key soil chemical parameters. We questioned how these abiotic variables influence the turnover of bacterial, archaeal and fungal communities and how they impact the long-term patterns of changes of the aforementioned soil communities. Interestingly, we found that the bacterial and fungal β-diversities are quite stable over time, whereas archaeal diversity showed significantly higher fluctuations. These fluctuations were reflected in temporal turnover caused by soil management through addition of N-fertilizers.

Conclusions

Our study showed that management practices applied to agricultural soils might not significantly affect the bacterial and fungal communities, but cause slow and long-term changes in the abundance and structure of the archaeal community. Moreover, the results suggest that, to different extents, abiotic and biotic factors determine the community assembly of archaeal, bacterial and fungal communities.     

青藏高原北部戈壁植物群落物种、功能与系统发育β多样性分布格局及其影响因素

戈壁荒漠广泛分布于全球干旱和极旱区域, 是我国陆地生态系统的重要组成部分。由于自然环境恶劣和交通条件限制, 目前有关戈壁植物群落物种、功能和系统发育等多维度β多样性形成机制的系统研究还很缺乏, 严重制约着对戈壁植物多样性维持机制的认知。本文以青藏高原北部61个典型戈壁生境植物群落为研究对象, 通过构建系统发育树和测量8个关键功能性状, 获取戈壁生境的物种、功能和系统发育β多样性, 比较3个维度β多样性格局与零模型的差异, 同时量化环境距离和地理距离对其的相对影响, 以探讨戈壁植物多样性的形成机制。结果显示: (1)戈壁植物的物种、功能和系统发育β多样性均表现出显著的距离衰减效应; (2)戈壁植物的物种、功能和系统发育β多样性均表现为非随机的格局; (3)由于功能性状趋同进化, 植物功能和系统发育β多样性变化趋势并不一致; (4)环境差异对植物3个维度β多样性均有着比空间距离更为重要的影响, 且土壤含水量、地表砾石盖度等局域生境因素的影响比气候更为强烈。以上结果表明, 戈壁植物的β多样性可能主要由局域生境过滤作用控制, 且不同维度的β多样性分布格局并不一致。     

Elevational Gradients in β-Diversity Reflect Variation in the Strength of Local Community Assembly Mechanisms across Spatial Scales

Despite long-standing interest in elevational-diversity gradients, little is known about the processes that cause changes in the compositional variation of communities (β-diversity) across elevations. Recent studies have suggested that β-diversity gradients are driven by variation in species pools, rather than by variation in the strength of local community assembly mechanisms such as dispersal limitation, environmental filtering, or local biotic interactions. However, tests of this hypothesis have been limited to very small spatial scales that limit inferences about how the relative importance of assembly mechanisms may change across spatial scales. Here, we test the hypothesis that scale-dependent community assembly mechanisms shape biogeographic β-diversity gradients using one of the most well-characterized elevational gradients of tropical plant diversity. Using an extensive dataset on woody plant distributions along a 4,000-m elevational gradient in the Bolivian Andes, we compared observed patterns of β-diversity to null-model expectations. β-deviations (standardized differences from null values) were used to measure the relative effects of local community assembly mechanisms after removing sampling effects caused by variation in species pools. To test for scale-dependency, we compared elevational gradients at two contrasting spatial scales that differed in the size of local assemblages and regions by at least an order of magnitude. Elevational gradients in β-diversity persisted after accounting for regional variation in species pools. Moreover, the elevational gradient in β-deviations changed with spatial scale. At small scales, local assembly mechanisms were detectable, but variation in species pools accounted for most of the elevational gradient in β-diversity. At large spatial scales, in contrast, local assembly mechanisms were a dominant force driving changes in β-diversity. In contrast to the hypothesis that variation in species pools alone drives β-diversity gradients, we show that local community assembly mechanisms contribute strongly to systematic changes in β-diversity across elevations. We conclude that scale-dependent variation in community assembly mechanisms underlies these iconic gradients in global biodiversity.     

Comparison of Diversities and Compositions of Bacterial Populations Inhabiting Natural Forest Soils

The diversity and composition of soil bacterial communities were compared among six Austrian natural forests, including oak-hornbeam, spruce-fir-beech, and Austrian pine forests, using terminal restriction fragment length polymorphism (T-RFLP, or TRF) analysis and sequence analysis of 16S rRNA genes. The forests studied differ greatly in soil chemical characteristics, microbial biomass, and nutrient turnover rates. The aim of this study was to relate these differences to the composition of the bacterial communities inhabiting the individual forest soils. Both TRF profiling and clone sequence analysis revealed that the bacterial communities in soils under Austrian pine forests, representing azonal forest types, were distinct from those in soils under zonal oak-hornbeam and spruce-fir-beech forests, which were more similar in community composition. Clones derived from an Austrian pine forest soil were mostly affiliated with high-G+C gram-positive bacteria (49%), followed by members of the α-Proteobacteria (20%) and the Holophaga/Acidobacterium group (12%). Clones in libraries from oak-hornbeam and spruce-fir-beech forest soils were mainly related to the Holophaga/Acidobacterium group (28 and 35%), followed by members of the Verrucomicrobia (24%) and the α-Proteobacteria (27%), respectively. The soil bacterial communities in forests with distinct vegetational and soil chemical properties appeared to be well differentiated based on 16S rRNA gene phylogeny. In particular, the outstanding position of the Austrian pine forests, which are determined by specific soil conditions, was reflected in the bacterial community composition.     

Discriminating ecological processes affecting different dimensions of α‐ and β‐diversity in desert plant communities

Understanding the spatial distribution of plant diversity and its drivers are major challenges in biogeography and conservation biology. Integrating multiple facets of biodiversity (e.g., taxonomic, phylogenetic, and functional biodiversity) may advance our understanding on how community assembly processes drive the distribution of biodiversity. In this study, plant communities in 60 sampling plots in desert ecosystems were investigated. The effects of local environment and spatial factors on the species, functional, and phylogenetic α‐ and β‐diversity (including turnover and nestedness components) of desert plant communities were investigated. The results showed that functional and phylogenetic α‐diversity were negatively correlated with species richness, and were significantly positively correlated with each other. Environmental filtering mainly influenced species richness and Rao quadratic entropy; phylogenetic α‐diversity was mainly influenced by dispersal limitation. Species and phylogenetic β‐diversity were mainly consisted of turnover component. The functional β‐diversity and its turnover component were mainly influenced by environmental factors, while dispersal limitation dominantly effected species and phylogenetic β‐diversity and their turnover component of species and phylogenetic β‐diversity. Soil organic carbon and soil pH significantly influenced different dimensions of α‐diversity, and soil moisture, salinity, organic carbon, and total nitrogen significantly influenced different dimensions of α‐ and β‐diversity and their components. Overall, it appeared that the relative influence of environmental and spatial factors on taxonomic, functional, and phylogenetic diversity differed at the α and β scales. Quantifying α‐ and β‐diversity at different biodiversity dimensions can help researchers to more accurately assess patterns of diversity and community assembly.     

Reproducibility and quantitation of amplicon sequencing-based detection

To determine the reproducibility and quantitation of the amplicon sequencing-based detection approach for analyzing microbial community structure, a total of 24 microbial communities from a long-term global change experimental site were examined. Genomic DNA obtained from each community was used to amplify 16S rRNA genes with two or three barcode tags as technical replicates in the presence of a small quantity (0.1% wt/wt) of genomic DNA from Shewanella oneidensis MR-1 as the control. The technical reproducibility of the amplicon sequencing-based detection approach is quite low, with an average operational taxonomic unit (OTU) overlap of 17.2%±2.3% between two technical replicates, and 8.2%±2.3% among three technical replicates, which is most likely due to problems associated with random sampling processes. Such variations in technical replicates could have substantial effects on estimating β-diversity but less on α-diversity. A high variation was also observed in the control across different samples (for example, 66.7-fold for the forward primer), suggesting that the amplicon sequencing-based detection approach could not be quantitative. In addition, various strategies were examined to improve the comparability of amplicon sequencing data, such as increasing biological replicates, and removing singleton sequences and less-representative OTUs across biological replicates. Finally, as expected, various statistical analyses with preprocessed experimental data revealed clear differences in the composition and structure of microbial communities between warming and non-warming, or between clipping and non-clipping. Taken together, these results suggest that amplicon sequencing-based detection is useful in analyzing microbial community structure even though it is not reproducible and quantitative. However, great caution should be taken in experimental design and data interpretation when the amplicon sequencing-based detection approach is used for quantitative analysis of the β-diversity of microbial communities.     

Indoor-Air Microbiome in an Urban Subway Network: Diversity and Dynamics

Subway systems are indispensable for urban societies, but microbiological characteristics of subway aerosols are relatively unknown. Previous studies investigating microbial compositions in subways employed methodologies that underestimated the diversity of microbial exposure for commuters, with little focus on factors governing subway air microbiology, which may have public health implications. Here, a culture-independent approach unraveling the bacterial diversity within the urban subway network in Hong Kong is presented. Aerosol samples from multiple subway lines and outdoor locations were collected. Targeting the 16S rRNA gene V4 region, extensive taxonomic diversity was found, with the most common bacterial genera in the subway environment among those associated with skin. Overall, subway lines harbored different phylogenetic communities based on α- and β-diversity comparisons, and closer inspection suggests that each community within a line is dependent on architectural characteristics, nearby outdoor microbiomes, and connectedness with other lines. Microbial diversities and assemblages also varied depending on the day sampled, as well as the time of day, and changes in microbial communities between peak and nonpeak commuting hours were attributed largely to increases in skin-associated genera in peak samples. Microbial diversities within the subway were influenced by temperature and relative humidity, while carbon dioxide levels showed a positive correlation with abundances of commuter-associated genera. This Hong Kong data set and communities from previous studies conducted in the United States formed distinct community clusters, indicating that additional work is required to unravel the mechanisms that shape subway microbiomes around the globe.     

How Does Conversion of Natural Tropical Rainforest Ecosystems Affect Soil Bacterial and Fungal Communities in the Nile River Watershed of Uganda?

Uganda''s forests are globally important for their conservation values but are under pressure from increasing human population and consumption. In this study, we examine how conversion of natural forest affects soil bacterial and fungal communities. Comparisons in paired natural forest and human-converted sites among four locations indicated that natural forest soils consistently had higher pH, organic carbon, nitrogen, and calcium, although variation among sites was large. Despite these differences, no effect on the diversity of dominant taxa for either bacterial or fungal communities was detected, using polymerase chain reaction-denaturing gradient gel electrophoresis (PCR-DGGE). Composition of fungal communities did generally appear different in converted sites, but surprisingly, we did not observe a consistent pattern among sites. The spatial distribution of some taxa and community composition was associated with soil pH, organic carbon, phosphorus and sodium, suggesting that changes in soil communities were nuanced and require more robust metagenomic methods to understand the various components of the community. Given the close geographic proximity of the paired sampling sites, the similarity between natural and converted sites might be due to continued dispersal between treatments. Fungal communities showed greater environmental differentiation than bacterial communities, particularly according to soil pH. We detected biotic homogenization in converted ecosystems and substantial contribution of β-diversity to total diversity, indicating considerable geographic structure in soil biota in these forest communities. Overall, our results suggest that soil microbial communities are relatively resilient to forest conversion and despite a substantial and consistent change in the soil environment, the effects of conversion differed widely among sites. The substantial difference in soil chemistry, with generally lower nutrient quantity in converted sites, does bring into question, how long this resilience will last.     

Estimating Bacterial Diversity for Ecological Studies: Methods,Metrics, and Assumptions

Methods to estimate microbial diversity have developed rapidly in an effort to understand the distribution and diversity of microorganisms in natural environments. For bacterial communities, the 16S rRNA gene is the phylogenetic marker gene of choice, but most studies select only a specific region of the 16S rRNA to estimate bacterial diversity. Whereas biases derived from from DNA extraction, primer choice and PCR amplification are well documented, we here address how the choice of variable region can influence a wide range of standard ecological metrics, such as species richness, phylogenetic diversity, β-diversity and rank-abundance distributions. We have used Illumina paired-end sequencing to estimate the bacterial diversity of 20 natural lakes across Switzerland derived from three trimmed variable 16S rRNA regions (V3, V4, V5). Species richness, phylogenetic diversity, community composition, β-diversity, and rank-abundance distributions differed significantly between 16S rRNA regions. Overall, patterns of diversity quantified by the V3 and V5 regions were more similar to one another than those assessed by the V4 region. Similar results were obtained when analyzing the datasets with different sequence similarity thresholds used during sequences clustering and when the same analysis was used on a reference dataset of sequences from the Greengenes database. In addition we also measured species richness from the same lake samples using ARISA Fingerprinting, but did not find a strong relationship between species richness estimated by Illumina and ARISA. We conclude that the selection of 16S rRNA region significantly influences the estimation of bacterial diversity and species distributions and that caution is warranted when comparing data from different variable regions as well as when using different sequencing techniques.     

Phyllosphere Bacterial Community of Floating Macrophytes in Paddy Soil Environments as Revealed by Illumina High-Throughput Sequencing

The phyllosphere of floating macrophytes in paddy soil ecosystems, a unique habitat, may support large microbial communities but remains largely unknown. We took Wolffia australiana as a representative floating plant and investigated its phyllosphere bacterial community and the underlying driving forces of community modulation in paddy soil ecosystems using Illumina HiSeq 2000 platform-based 16S rRNA gene sequence analysis. The results showed that the phyllosphere of W. australiana harbored considerably rich communities of bacteria, with Proteobacteria and Bacteroidetes as the predominant phyla. The core microbiome in the phyllosphere contained genera such as Acidovorax, Asticcacaulis, Methylibium, and Methylophilus. Complexity of the phyllosphere bacterial communities in terms of class number and α-diversity was reduced compared to those in corresponding water and soil. Furthermore, the bacterial communities exhibited structures significantly different from those in water and soil. These findings and the following redundancy analysis (RDA) suggest that species sorting played an important role in the recruitment of bacterial species in the phyllosphere. The compositional structures of the phyllosphere bacterial communities were modulated predominantly by water physicochemical properties, while the initial soil bacterial communities had limited impact. Taken together, the findings from this study reveal the diversity and uniqueness of the phyllosphere bacterial communities associated with the floating macrophytes in paddy soil environments.     

Microdiversity characterizes prevalent phylogenetic clades in the glacier-fed stream microbiome

Glacier-fed streams (GFSs) are extreme and rapidly vanishing ecosystems, and yet they harbor diverse microbial communities. Although our understanding of the GFS microbiome has recently increased, we do not know which microbial clades are ecologically successful in these ecosystems, nor do we understand potentially underlying mechanisms. Ecologically successful clades should be more prevalent across GFSs compared to other clades, which should be reflected as clade-wise distinctly low phylogenetic turnover. However, methods to assess such patterns are currently missing. Here we developed and applied a novel analytical framework, “phyloscore analysis”, to identify clades with lower spatial phylogenetic turnover than other clades in the sediment microbiome across twenty GFSs in New Zealand. These clades constituted up to 44% and 64% of community α-diversity and abundance, respectively. Furthermore, both their α-diversity and abundance increased as sediment chlorophyll a decreased, corroborating their ecological success in GFS habitats largely devoid of primary production. These clades also contained elevated levels of putative microdiversity than others, which could potentially explain their high prevalence in GFSs. This hitherto unknown microdiversity may be threatened as glaciers shrink, urging towards further genomic and functional exploration of the GFS microbiome.Subject terms: Climate-change ecology, Microbial ecology, Community ecology, Theoretical ecology     

Microbial Biomass and Activity in Lead-Contaminated Soil

Microbial community diversity, potential microbial activity, and metal resistance were determined in three soils whose lead contents ranged from 0.00039 to 48 mmol of Pb kg of soil⁻¹. Biomass levels were directly related to lead content. A molecular analysis of 16S rRNAs suggested that each soil contained a complex, diverse microbial community. A statistical analysis of the phospholipid fatty acids indicated that the community in the soil having the highest lead content was not related to the communities in the other soils. All of the soils contained active microbial populations that mineralized [¹⁴C]glucose. In all samples, 10 to 15% of the total culturable bacteria were Pb resistant and had MIC of Pb for growth of 100 to 150 μM.