Environment drives high phylogenetic turnover among oceanic bacterial communities

Selection, drift and community interactions shape microbial biogeographic patterns in the Pacific Ocean

Despite accumulating data on microbial biogeographic patterns in terrestrial and aquatic environments, we still lack a comprehensive understanding of how these patterns establish, in particular in ocean basins. Here we show the relative significance of the ecological mechanisms selection, dispersal and drift for shaping the composition of microbial communities in the Pacific Ocean over a transect of 12,400 km between subantarctic and subarctic regions. In the epipelagic, homogeneous selection contributes 50-60% and drift least to the three mechanism for the assembly of prokaryotic communities whereas in the upper mesopelagic, drift is relatively most important for the particle-associated subcommunities. Temperature is important for the relative significance of homogeneous selection and dispersal limitation for community assembly. The relative significance of both mechanisms was inverted with increasing temperature difference along the transect. For eukaryotes >8 µm, homogeneous selection is also the most important mechanisms at two epipelagic depths whereas at all other depths drift is predominant. As species interactions are essential for structuring microbial communities we further analyzed co-occurrence-based community metrics to assess biogeographic patterns over the transect. These interaction-adjusted indices explained much better variations in microbial community composition as a function of abiotic and biotic variables than compositional or phylogenetic distance measures like Bray-Curtis or UniFrac. Our analyses are important to better understand assembly processes of microbial communities in the upper layers of the largest ocean and how they adapt to effectively perform in global biogeochemical processes. Similar principles presumably act upon microbial community assembly in other ocean basins.

Structure and function of the Arctic and Antarctic marine microbiota as revealed by metagenomics

BACKGROUND

The Arctic and Antarctic are the two most geographically distant bioregions on earth. Recent sampling efforts and following metagenomics have shed light on the global ocean microbial diversity and function, yet the microbiota of polar regions has not been included in such global analyses.

RESULTS

Here a metagenomic study of seawater samples (n = 60) collected from different depths at 28 locations in the Arctic and Antarctic zones was performed, together with metagenomes from the Tara Oceans. More than 7500 (19%) polar seawater-derived operational taxonomic units could not be identified in the Tara Oceans datasets, and more than 3,900,000 protein-coding gene orthologs had no hits in the Ocean Microbial Reference Gene Catalog. Analysis of 214 metagenome assembled genomes (MAGs) recovered from the polar seawater microbiomes, revealed strains that are prevalent in the polar regions while nearly undetectable in temperate seawater. Metabolic pathway reconstruction for these microbes suggested versatility for saccharide and lipids biosynthesis, nitrate and sulfate reduction, and CO₂ fixation. Comparison between the Arctic and Antarctic microbiomes revealed that antibiotic resistance genes were enriched in the Arctic while functions like DNA recombination were enriched in the Antarctic.

CONCLUSIONS

Our data highlight the occurrence of dominant and locally enriched microbes in the Arctic and Antarctic seawater with unique functional traits for environmental adaption, and provide a foundation for analyzing the global ocean microbiome in a more complete perspective. Video abstract.

The Role of Ocean Currents in the Temperature Selection of Plankton: Insights from an Individual-Based Model

Biogeography studies that correlate the observed distribution of organisms to environmental variables are typically based on local conditions. However, in cases with substantial translocation, like planktonic organisms carried by ocean currents, selection may happen upstream and local environmental factors may not be representative of those that shaped the local population. Here we use an individual-based model of microbes in the global surface ocean to explore this effect for temperature. We simulate up to 25 million individual cells belonging to up to 50 species with different temperature optima. Microbes are moved around the globe based on a hydrodynamic model, and grow and die based on local temperature. We quantify the role of currents using the “advective temperature differential” metric, which is the optimum temperature of the most abundant species from the model with advection minus that from the model without advection. This differential depends on the location and can be up to 4°C. Poleward-flowing currents, like the Gulf Stream, generally experience cooling and the differential is positive. We apply our results to three global datasets. For observations of optimum growth temperature of phytoplankton, accounting for the effect of currents leads to a slightly better agreement with observations, but there is large variability and the improvement is not statistically significant. For observed Prochlorococcus ecotype ratios and metagenome nucleotide divergence, accounting for advection improves the correlation significantly, especially in areas with relatively strong poleward or equatorward currents.

Glycoside Hydrolases across Environmental Microbial Communities

Across many environments microbial glycoside hydrolases support the enzymatic processing of carbohydrates, a critical function in many ecosystems. Little is known about how the microbial composition of a community and the potential for carbohydrate processing relate to each other. Here, using 1,934 metagenomic datasets, we linked changes in community composition to variation of potential for carbohydrate processing across environments. We were able to show that each ecosystem-type displays a specific potential for carbohydrate utilization. Most of this potential was associated with just 77 bacterial genera. The GH content in bacterial genera is best described by their taxonomic affiliation. Across metagenomes, fluctuations of the microbial community structure and GH potential for carbohydrate utilization were correlated. Our analysis reveals that both deterministic and stochastic processes contribute to the assembly of complex microbial communities.

Relative Roles of Deterministic and Stochastic Processes in Driving the Vertical Distribution of Bacterial Communities in a Permafrost Core from the Qinghai-Tibet Plateau, China

Understanding the processes that influence the structure of biotic communities is one of the major ecological topics, and both stochastic and deterministic processes are expected to be at work simultaneously in most communities. Here, we investigated the vertical distribution patterns of bacterial communities in a 10-m-long soil core taken within permafrost of the Qinghai-Tibet Plateau. To get a better understanding of the forces that govern these patterns, we examined the diversity and structure of bacterial communities, and the change in community composition along the vertical distance (spatial turnover) from both taxonomic and phylogenetic perspectives. Measures of taxonomic and phylogenetic beta diversity revealed that bacterial community composition changed continuously along the soil core, and showed a vertical distance-decay relationship. Multiple stepwise regression analysis suggested that bacterial alpha diversity and phylogenetic structure were strongly correlated with soil conductivity and pH but weakly correlated with depth. There was evidence that deterministic and stochastic processes collectively drived bacterial vertically-structured pattern. Bacterial communities in five soil horizons (two originated from the active layer and three from permafrost) of the permafrost core were phylogenetically random, indicator of stochastic processes. However, we found a stronger effect of deterministic processes related to soil pH, conductivity, and organic carbon content that were structuring the bacterial communities. We therefore conclude that the vertical distribution of bacterial communities was governed primarily by deterministic ecological selection, although stochastic processes were also at work. Furthermore, the strong impact of environmental conditions (for example, soil physicochemical parameters and seasonal freeze-thaw cycles) on these communities underlines the sensitivity of permafrost microorganisms to climate change and potentially subsequent permafrost thaw.

Biodegradation of dichlorodiphenyltrichloroethanes (DDTs) and hexachlorocyclohexanes (HCHs) with plant and nutrients and their effects on the microbial ecological kinetics

Four pilot-scale test mesocosms were conducted for the remediation of organochlorine pesticides (OCPs)-contaminated aged soil. The results indicate that the effects on degradation of hexachlorocyclohexanes (HCHs) and dichlorodiphenyltrichloroethanes (DDTs) were in the following order: nutrients/plant bioaugmentation (81.18 % for HCHs; 85.4 % for DDTs) > nutrients bioaugmentation > plant bioaugmentation > only adding water > control, and nutrients/plant bioaugmentation greatly enhanced the degradation of HCHs (81.18 %) and DDTs (85.4 %). The bacterial community structure, diversity and composition were assessed by 454-pyrosequencing of 16S recombinant RNA (rRNA), whereas the abundance of linA gene was determined by quantitative polymerase chain reaction. Distinct differences in bacterial community composition, structure, and diversity were a function of remediation procedure. Predictability of HCH/DDT degradation in soils was also investigated. A positive correlation between linA gene abundance and the removal ratio of HCHs was indicated by correlation analyses. A similar relationship was also confirmed between the degradation of HCHs/DDTs and the abundance of some assemblages (Gammaproteobacteria and Flavobacteria). Our results offer microbial ecological insight into the degradation of HCHs and DDTs in aged contaminated soil, which is helpful for the intensification of bioremediation through modifying plant-microbe patterns, and cessation of costly and time-consuming assays.

Abundance of broad bacterial taxa in the sargasso sea explained by environmental conditions but not water mass

To explore the potential linkage between distribution of marine bacterioplankton groups, environmental conditions, and water mass, we investigated the factors determining the abundance of bacterial taxa across the hydrographically complex Subtropical Convergence Zone in the Sargasso Sea. Based on information from 16S rRNA gene clone libraries from various locations and two depths, abundances of the predominant taxa (eubacteria, Archaea, Alphaproteobacteria, Gammaproteobacteria, Bacteroidetes, and the Roseobacter, SAR11, and SAR86 clades) were quantified by real-time PCR. In addition, the abundances of Synechococcus, Prochlorococcus, and picoalgae were determined by flow cytometry. Linear multiple-regression models determining the relative effects of eight environmental variables and of water mass explained 35 to 86% of the variation in abundance of the quantified taxa, even though only one to three variables were significantly related to any particular taxon's abundance. Most of the variation in abundance was explained by depth and chlorophyll a. The predominant phototrophs, Prochlorococcus and picoalgae, were negatively correlated with phosphate, whereas eubacteria, heterotrophic bacteria, and SAR86 were negatively correlated with nitrite. Water mass showed limited importance for explaining the abundance of the taxonomical groups (significant only for Roseobacter, explaining 14% of the variation). The results suggest the potential for predicting the abundance of broad bacterioplankton groups throughout the Sargasso Sea using only a few environmental parameters.

Phylogenetic differences in attached and free-living bacterial communities in a temperate coastal lagoon during summer, revealed via high-throughput 16S rRNA gene sequencing

Most of what is known about coastal free-living and attached bacterial diversity is based on open coasts, with high particulate and nutrient riverine supply, terrestrial runoffs, and anthropogenic activities. The Magdalen Islands in the Gulf of St. Lawrence (Canada) are dominated by shallow lagoons with small, relatively pristine catchments and no freshwater input apart from rain. Such conditions provided an opportunity to investigate coastal free-living and attached marine bacterial diversity in the absence of confounding effects of steep freshwater gradients. We found significant differences between the two communities and marked temporal patterns in both. Taxonomic richness and diversity were greater in the attached than in the free-living community, increasing over summer, especially within the least abundant bacterial phyla. The highest number of reads fell within the SAR 11 clade (Pelagibacter, Alphaproteobacteria), which dominated free-living communities. The attached communities had deeper phylum-level diversity than the free-living fraction. Distance-based redundancy analysis indicated that the particulate organic matter (POM) concentration was the main variable separating early and late summer samples with salinity and temperature changes also significantly correlated to bacterial community structure. Our approach using high-throughput sequencing detected differences in free-living versus attached bacteria in the absence of riverine input, in keeping with the concept that marine attached communities are distinct from cooccurring free-living taxa. This diversity likely reflects the diverse microhabitats of available particles, implying that the total bacterial diversity in coastal systems is linked to particle supply and variability, with implications for understanding microbial biodiversity in marine systems.

A unifying quantitative framework for exploring the multiple facets of microbial biodiversity across diverse scales

Recent developments of molecular tools have revolutionized our knowledge of microbial biodiversity by allowing detailed exploration of its different facets and generating unprecedented amount of data. One key issue with such large datasets is the development of diversity measures that cope with different data outputs and allow comparison of biodiversity across different scales. Diversity has indeed three components: local (α), regional (γ) and the overall difference between local communities (β). Current measures of microbial diversity, derived from several approaches, provide complementary but different views. They only capture the β component of diversity, compare communities in a pairwise way, consider all species as equivalent or lack a mathematically explicit relationship among the α, β and γ components. We propose a unified quantitative framework based on the Rao quadratic entropy, to obtain an additive decomposition of diversity (γ = α + β), so the three components can be compared, and that integrate the relationship (phylogenetic or functional) among Microbial Diversity Units that compose a microbial community. We show how this framework is adapted to all types of molecular data, and we highlight crucial issues in microbial ecology that would benefit from this framework and propose ready-to-use R-functions to easily set up our approach.

Understanding bacteriophage specificity in natural microbial communities

Studying the coevolutionary dynamics between bacteria and the bacteriophage viruses that infect them is critical to understanding both microbial diversity and ecosystem functioning. Phages can play a key role in shaping bacterial population dynamics and can significantly alter both intra- and inter-specific competition among bacterial hosts. Predicting how phages might influence community stability and apparent competition, however, requires an understanding of how bacteria-phage interaction networks evolve as a function of host diversity and community dynamics. Here, we first review the progress that has been made in understanding phage specificity, including the use of experimental evolution, we then introduce a new dataset on natural bacteriophages collected from the phyllosphere of horse chestnut trees, and finally we highlight that bacterial sensitivity to phage is rarely a binary trait and that this variation should be taken into account and reported. We emphasize that there is currently insufficient evidence to make broad generalizations about phage host range in natural populations, the limits of phage adaptation to novel hosts, or the implications of phage specificity in shaping microbial communities. However, the combination of experimental and genomic approaches with the study of natural communities will allow new insight to the evolution and impact of phage specificity within complex bacterial communities.

Predictable bacterial composition and hydrocarbon degradation in Arctic soils following diesel and nutrient disturbance

Increased exploration and exploitation of resources in the Arctic is leading to a higher risk of petroleum contamination. A number of Arctic microorganisms can use petroleum for growth-supporting carbon and energy, but traditional approaches for stimulating these microorganisms (for example, nutrient addition) have varied in effectiveness between sites. Consistent environmental controls on microbial community response to disturbance from petroleum contaminants and nutrient amendments across Arctic soils have not been identified, nor is it known whether specific taxa are universally associated with efficient bioremediation. In this study, we contaminated 18 Arctic soils with diesel and treated subsamples of each with monoammonium phosphate (MAP), which has successfully stimulated degradation in some contaminated Arctic soils. Bacterial community composition of uncontaminated, diesel-contaminated and diesel+MAP soils was assessed through multiplexed 16S (ribosomal RNA) rRNA gene sequencing on an Ion Torrent Personal Genome Machine, while hydrocarbon degradation was measured by gas chromatography analysis. Diversity of 16S rRNA gene sequences was reduced by diesel, and more so by the combination of diesel and MAP. Actinobacteria dominated uncontaminated soils with <10% organic matter, while Proteobacteria dominated higher-organic matter soils, and this pattern was exaggerated following disturbance. Degradation with and without MAP was predictable by initial bacterial diversity and the abundance of specific assemblages of Betaproteobacteria, respectively. High Betaproteobacteria abundance was positively correlated with high diesel degradation in MAP-treated soils, suggesting this may be an important group to stimulate. The predictability with which bacterial communities respond to these disturbances suggests that costly and time-consuming contaminated site assessments may not be necessary in the future.

Pole-to-pole biogeography of surface and deep marine bacterial communities

The Antarctic and Arctic regions offer a unique opportunity to test factors shaping biogeography of marine microbial communities because these regions are geographically far apart, yet share similar selection pressures. Here, we report a comprehensive comparison of bacterioplankton diversity between polar oceans, using standardized methods for pyrosequencing the V6 region of the small subunit ribosomal (SSU) rRNA gene. Bacterial communities from lower latitude oceans were included, providing a global perspective. A clear difference between Southern and Arctic Ocean surface communities was evident, with 78% of operational taxonomic units (OTUs) unique to the Southern Ocean and 70% unique to the Arctic Ocean. Although polar ocean bacterial communities were more similar to each other than to lower latitude pelagic communities, analyses of depths, seasons, and coastal vs. open waters, the Southern and Arctic Ocean bacterioplankton communities consistently clustered separately from each other. Coastal surface Southern and Arctic Ocean communities were more dissimilar from their respective open ocean communities. In contrast, deep ocean communities differed less between poles and lower latitude deep waters and displayed different diversity patterns compared with the surface. In addition, estimated diversity (Chao1) for surface and deep communities did not correlate significantly with latitude or temperature. Our results suggest differences in environmental conditions at the poles and different selection mechanisms controlling surface and deep ocean community structure and diversity. Surface bacterioplankton may be subjected to more short-term, variable conditions, whereas deep communities appear to be structured by longer water-mass residence time and connectivity through ocean circulation.

Quantifying phylogenetic beta diversity: distinguishing between 'true' turnover of lineages and phylogenetic diversity gradients

The evolutionary dissimilarity between communities (phylogenetic beta diversity PBD) has been increasingly explored by ecologists and biogeographers to assess the relative roles of ecological and evolutionary processes in structuring natural communities. Among PBD measures, the PhyloSor and UniFrac indices have been widely used to assess the level of turnover of lineages over geographical and environmental gradients. However, these indices can be considered as ‘broad-sense’ measures of phylogenetic turnover as they incorporate different aspects of differences in evolutionary history between communities that may be attributable to phylogenetic diversity gradients. In the present study, we extend an additive partitioning framework proposed for compositional beta diversity to PBD. Specifically, we decomposed the PhyloSor and UniFrac indices into two separate components accounting for ‘true’ phylogenetic turnover and phylogenetic diversity gradients, respectively. We illustrated the relevance of this framework using simple theoretical and archetypal examples, as well as an empirical study based on coral reef fish communities. Overall, our results suggest that using PhyloSor and UniFrac may greatly over-estimate the level of spatial turnover of lineages if the two compared communities show contrasting levels of phylogenetic diversity. We therefore recommend that future studies use the ‘true’ phylogenetic turnover component of these indices when the studied communities encompass a large phylogenetic diversity gradient.