Role of functionally dominant species in varying environmental regimes: evidence for the performance-enhancing effect of biodiversity

Synthetic phylogenetically diverse communities promote denitrification and stability

Denitrification is an important process of the global nitrogen cycle as some of its intermediates are environmentally important or related to global warming. However, how the phylogenetic diversity of denitrifying communities affects their denitrification rates and temporal stability remains unclear. Here we selected denitrifiers based on their phylogenetic distance to construct two groups of synthetic denitrifying communities: one closely related (CR) group with all strains from the genus Shewanella and the other distantly related (DR) group with all constituents from different genera. All synthetic denitrifying communities (SDCs) were experimentally evolved for 200 generations. The results showed that high phylogenetic diversity followed by experimental evolution promoted the function and stability of synthetic denitrifying communities. Specifically, the productivity and denitrification rates were significantly (P < 0.05) higher with Paracocus denitrificans as the dominant species (since the 50th generation) in the DR community than those in the CR community. The DR community also showed significantly (t = 7.119, df = 10, P < 0.001) higher stability through overyielding and asynchrony of species fluctuations, and showed more complementarity than the CR group during the experimental evolution. This study has important implications for applying synthetic communities to remediate environmental problems and mitigate greenhouse gas emissions.

Synthetic Denitrifying Communities Reveal a Positive and Dynamic Biodiversity-Ecosystem Functioning Relationship during Experimental Evolution

Biodiversity is vital for ecosystem functions and services, and many studies have reported positive, negative, or neutral biodiversity-ecosystem functioning (BEF) relationships in plant and animal systems. However, if the BEF relationship exists and how it evolves remains elusive in microbial systems. Here, we selected 12 Shewanella denitrifiers to construct synthetic denitrifying communities (SDCs) with a richness gradient spanning 1 to 12 species, which were subjected to approximately 180 days (with 60 transfers) of experimental evolution with generational changes in community functions continuously tracked. A significant positive correlation was observed between community richness and functions, represented by productivity (biomass) and denitrification rate, however, such a positive correlation was transient, only significant in earlier days (0 to 60) during the evolution experiment (180 days). Also, we found that community functions generally increased throughout the evolution experiment. Furthermore, microbial community functions with lower richness exhibited greater increases than those with higher richness. Biodiversity effect analysis revealed positive BEF relationships largely attributable to complementary effects, which were more pronounced in communities with lower richness than those with higher richness. This study is one of the first studies that advances our understanding of BEF relationships and their evolutionary mechanisms in microbial systems, highlighting the crucial role of evolution in predicting the BEF relationship in microbial systems. IMPORTANCE Despite the consensus that biodiversity supports ecosystem functioning, not all experimental models of macro-organisms support this notion with positive, negative, or neutral biodiversity-ecosystem functioning (BEF) relationships reported. The fast-growing, metabolically versatile, and easy manipulation nature of microbial communities allows us to explore well the BEF relationship and further interrogate if the BEF relationship remains constant during long-term community evolution. Here, we constructed multiple synthetic denitrifying communities (SDCs) by randomly selecting species from a candidate pool of 12 Shewanella denitrifiers. These SDCs differ in species richness, spanning 1 to 12 species, and were monitored continuously for community functional shifts during approximately 180-day parallel cultivation. We demonstrated that the BEF relationship was dynamic with initially (day 0 to 60) greater productivity and denitrification among SDCs of higher richness. However, such pattern was reversed thereafter with greater productivity and denitrification increments in lower-richness SDCs, likely due to a greater accumulation of beneficial mutations during the experimental evolution.

Independent Effects of Species Removal and Asynchrony on Invariability of an Intertidal Rocky Shore Community

Ecological stability depends on interactions between different levels of biological organization. The insurance effects occur when increasing species diversity leads to more temporally invariable (i.e., more stable) community-level properties, due in part to asynchronous population-level fluctuations. While the study of insurance effects has received considerable attention, the role of dominant species that contribute with particular functional traits across different level of organizations is less understood. Using a field-based manipulative experiment, we investigated how species richness and different types of parameters at the population level, such as the invariability of dominants, population invariability, and population asynchrony, influence the community invariability. The experiment involved the repetitive removal of the canopy forming alga Mazzaella laminarioides (hereafter “Mazzaella”) during 32 months in two rocky intertidal sites of northern-central Chile. We predicted that the invariability of dominants enhances community invariability, that the effect of multispecies population-level parameters on community invariability are dependent on species richness, and that subdominant algae are unable to fully compensate the loss of canopies of the dominant species. Biomass of algae and mobile invertebrates was quantified over time. We observed independent effects of Mazzaella removal and community-wide asynchrony on community invariability. While canopy removal reduced community invariability, population asynchrony boosted community invariability regardless of the presence of canopies. In addition, filamentous and foliose algae were unable to compensate the loss of biomass triggered by the experimental removal of Mazzaella. Canopy removal led to a severe decrement in the biomass of macrograzers, while, at the same time, increased the biomass of mesograzers. Asynchrony stemmed from compensatory trophic responses of mesograzers to increased abundances of opportunistic algae. Thus, further work on consumer-resource interactions will improve our understanding of the links between population- and community-level aspects of stability.

Formation, characterization and modeling of emergent synthetic microbial communities

Microbial communities colonize plant tissues and contribute to host function. How these communities form and how individual members contribute to shaping the microbial community are not well understood. Synthetic microbial communities, where defined individual isolates are combined, can serve as valuable model systems for uncovering the organizational principles of communities. Using genome-defined organisms, systematic analysis by computationally-based network reconstruction can lead to mechanistic insights and the metabolic interactions between species. In this study, 10 bacterial strains isolated from the Populus deltoides rhizosphere were combined and passaged in two different media environments to form stable microbial communities. The membership and relative abundances of the strains stabilized after around 5 growth cycles and resulted in just a few dominant strains that depended on the medium. To unravel the underlying metabolic interactions, flux balance analysis was used to model microbial growth and identify potential metabolic exchanges involved in shaping the microbial communities. These analyses were complemented by growth curves of the individual isolates, pairwise interaction screens, and metaproteomics of the community. A fast growth rate is identified as one factor that can provide an advantage for maintaining presence in the community. Final community selection can also depend on selective antagonistic relationships and metabolic exchanges. Revealing the mechanisms of interaction among plant-associated microorganisms provides insights into strategies for engineering microbial communities that can potentially increase plant growth and disease resistance. Further, deciphering the membership and metabolic potentials of a bacterial community will enable the design of synthetic communities with desired biological functions.

Diversity Manipulation of Psychrophilic Bacterial Consortia for Improved Biological Treatment of Medium-Strength Wastewater at Low Temperature

Psychrophilic bacteria are valuable biocatalysts to develop robust bioaugmentation formulations for enhanced wastewater treatment at low temperatures or fluctuating temperature conditions. Here, using different biodiversity indices [based on species richness (SR), phylogenetic diversity (PD) and functional diversity (FD)], we studied the effects of microbial diversity of artificial bacterial consortia on the biomass gross yields (measured through OD₆₀₀) and removal efficiency of soluble chemical oxygen demand (mg sCOD removed/mg sCOD introduced) in synthetic, medium-strength wastewater. We built artificial consortia out of one to six bacterial strains isolated at 4°C through combinatorial biodiversity experiments. Increasing species richness resulted in improved sCOD removal efficiency (i.e., 0.266 ± 0.146, 0.542 ± 0.155, 0.742 ± 0.136, 0.822 ± 0.019 for mono-, tri-, penta-and hexacultures, respectively) and higher biomass gross yields (i.e., 0.065 ± 0.052, 0.132 ± 0.046, 0.173 ± 0.049, 0.216 ± 0.019 for mono-, tri-, penta,- and hexacultures, respectively). This positive relationship between biodiversity, sCOD removal and biomass gross yield was also observed when considering metabolic profiling (functional diversity) or evolutionary relationships (phylogenetic diversity). The positive effect of biodiversity on sCOD removal efficiency could be attributed to the selection of a particular, best-performing species (i.e., Pedobacter sp.) as well as complementary use of carbon resources among consortia members (i.e., complementarity effects). Among the biodiversity indices, PD diversity metrics explained higher variation in sCOD removal than SR and FD diversity metrics. For a more effective bioaugmentation, our results stress the importance of using phylogenetically diverse consortia, with an increased degradation ability, instead of single pure cultures. Moreover, PD could be used as an assembly rule to guide the composition of mixed cultures for wastewater bioaugmentation under psychrophilic conditions.

Realistic scenarios of environmental disturbance lead to functionally important changes in benthic species-environment interactions

Changes in community structure concurrent with environmental forcing often form a precursor to changes in species diversity, and can have substantive consequences for ecosystem functioning. Here, we assess the effects of altered levels of evenness that are representative of different levels of eutrophication and changes in salinity associated with altered precipitation patterns, on the mediation of nutrient release by sediment-dwelling invertebrate communities. We find that an adjustment towards a more even distribution of species corresponds with an increase in sediment particle reworking that, in general, translates to increased levels of nutrient release. This response, however, is dependent on the functional role of each species in the community and is influenced by concomitant changes in salinity, especially when salinity extends beyond the range typically experienced by the community. Overall, our findings highlight the dynamic nature of species contributions to functioning and reinforce the importance of understanding when, and how, the mechanistic basis of species-environment interactions are modified as the influence of abiotic and biotic factors flex under periods of directional forcing.

Distinct Responses of the Nitrogen-Fixing Marine Cyanobacterium Trichodesmium to a Thermally Variable Environment as a Function of Phosphorus Availability

Surface temperature in the ocean is projected to be elevated and more variable in the future, which will interact with other environmental changes like reduced nutrient supplies. To explore these multiple stressor relationships, we tested the influence of thermal variation on the key marine diazotrophic cyanobacterium Trichodesmium erythraeum GBRTRLI101 as a function of the limiting nutrient phosphorus (P). Two constant temperature treatments represented current winter (22°C) and summer (30°C) mean values. Three variable temperature treatments fluctuated around the constant control values: Mean 22°C, either ± 2°C or ± 4°C; and mean 30°C ± 2°C. Each thermal treatment was grown under both P-replete (10 μmol/L) and P-limiting conditions (0.2 μmol/L). Effects of thermal variability on Trichodesmium were mainly found in the two winter variable temperature treatments (22°C ± 2°C or ± 4°C). P availability affected growth and physiology in all treatments and had significant interactions with temperature. P-replete cultures had higher growth and nitrogen and carbon fixation rates in the 22°C constant control, than in the corresponding variable treatments. However, physiological rates were not different in the P-replete constant and variable treatments at 30°C. In contrast, in P-limited cultures an advantage of constant temperature over variable temperature was not apparent. Phosphorus use efficiencies (PUE, mol N or C fixed h-1 mol cellular P-1) for nitrogen and carbon fixation were significantly elevated under P-limited conditions, and increased with temperature from 22 to 30°C, implying a potential advantage in a future warmer, P-limited environment. Taken together, these results imply that future increasing temperature and greater thermal variability could have significant feedback interactions with the projected intensification of P-limitation of marine N2-fixing cyanobacteria.

Predator coevolution and prey trait variability determine species coexistence

Predation is one of the key ecological mechanisms allowing species coexistence and influencing biological diversity. However, ecological processes are subject to contemporary evolutionary change, and the degree to which predation affects diversity ultimately depends on the interplay between evolution and ecology. Furthermore, ecological interactions that influence species coexistence can be altered by reciprocal coevolution especially in the case of antagonistic interactions such as predation or parasitism. Here we used an experimental evolution approach to test for the role of initial trait variation in the prey population and coevolutionary history of the predator in the ecological dynamics of a two-species bacterial community predated by a ciliate. We found that initial trait variation both at the bacterial and ciliate level enhanced species coexistence, and that subsequent trait evolutionary trajectories depended on the initial genetic diversity present in the population. Our findings provide further support to the notion that the ecology-centric view of diversity maintenance must be reinvestigated in light of recent findings in the field of eco-evolutionary dynamics.

Compensatory responses can alter the form of the biodiversity-function relation curve

There is now strong evidence that ecosystem properties are influenced by alterations in biodiversity. The consensus that has emerged from over two decades of research is that the form of the biodiversity-functioning relationship follows a saturating curve. However, the foundation from which these conclusions are drawn mostly stems from empirical investigations that have not accounted for post-extinction changes in community composition and structure, or how surviving species respond to new circumstances and modify their contribution to functioning. Here, we use marine sediment-dwelling invertebrate communities to experimentally assess whether post-extinction compensatory mechanisms (simulated by increasing species biomass) have the potential to alter biodiversity-ecosystem function relations. Consistent with recent numerical simulations, we find that the form of the biodiversity-function curve is dependent on whether or not compensatory responses are present, the cause and extent of extinction, and species density. When species losses are combined with the compensatory responses of surviving species, both community composition, dominance structure, and the pool and relative expression of functionally important traits change and affect species interactions and behaviour. These observations emphasize the importance of post-extinction community composition in determining the stability of ecosystem functioning following extinction. Our results caution against the use of the generalized biodiversity-function curve when generating probabilistic estimates of post-extinction ecosystem properties for practical application.

Species contributions to ecosystem process and function can be population dependent and modified by biotic and abiotic setting

There is unequivocal evidence that altered biodiversity, through changes in the expression and distribution of functional traits, can have large impacts on ecosystem properties. However, trait-based summaries of how organisms affect ecosystem properties often assume that traits show constancy within and between populations and that species contributions to ecosystem functioning are not overly affected by the presence of other species or variations in abiotic conditions. Here, we evaluate the validity of these assumptions using an experiment in which three geographically distinct populations of intertidal sediment-dwelling invertebrates are reciprocally substituted. We find that the mediation of macronutrient generation by these species can vary between different populations and show that changes in biotic and/or abiotic conditions can further modify functionally important aspects of the behaviour of individuals within a population. Our results demonstrate the importance of knowing how, when, and why traits are expressed and suggest that these dimensions of species functionality are not sufficiently well-constrained to facilitate the accurate projection of the functional consequences of change. Information regarding the ecological role of key species and assumptions about the form of species–environment interactions needs urgent refinement.

Resolving Species Level Changes in a Representative Soil Bacterial Community Using Microfluidic Quantitative PCR

Rapid advances in genome sequencing technologies enable determination of relative bacterial abundances and community composition, yet, changes at the species level remain difficult to detect despite importance for certain ecological inferences. We present a method for extraction and direct quantification of species composition of a predefined multispecies bacterial community using microfluidic-based quantitative real-time PCR (qPCR). We employ a nested PCR approach based on universal 16S rRNA gene pre-amplification followed by detection and quantification of absolute abundance of bacterial species using microfluidic array of parallel singleplex qPCR reactions. Present microfluidic qPCR supports 2,304 simultaneous reactions on a single chip, while automatic distribution of samples and reactants minimizes pipetting errors and technical variations. The utility of the method is illustrated using a synthetic soil bacterial community grown in two contrasting environments – sand microcosms and batch cultures. The protocol entails extraction of total nucleic acid, preparation of genomic DNA, and steps for qPCR assessment of bacterial community composition. This method provides specific and sensitive quantification of bacterial species requiring only 2 ng of community DNA. Optimized extraction protocol and preamplification step allow for rapid, quantitative, and simultaneous detection of candidate species with high throughput. The proposed method offers a simple and accurate alternative to present sequencing methods especially when absolute values of species abundance are required. Quantification of changes at the species level contributes to the mechanistic understanding of the roles of particular species in a bacterial community functioning.

Defining the functional traits that drive bacterial decomposer community productivity

Microbial communities are essential to a wide range of ecologically and industrially important processes. To control or predict how these communities function, we require a better understanding of the factors which influence microbial community productivity. Here, we combine functional resource use assays with a biodiversity–ecosystem functioning (BEF) experiment to determine whether the functional traits of constituent species can be used to predict community productivity. We quantified the abilities of 12 bacterial species to metabolise components of lignocellulose and then assembled these species into communities of varying diversity and composition to measure their productivity growing on lignocellulose, a complex natural substrate. A positive relationship between diversity and community productivity was caused by a selection effect whereby more diverse communities were more likely to contain two species that significantly improved community productivity. Analysis of functional traits revealed that the observed selection effect was primarily driven by the abilities of these species to degrade β-glucan. Our results indicate that by identifying the key functional traits underlying microbial community productivity we could improve industrial bioprocessing of complex natural substrates.

Specific arrangements of species dominance can be more influential than evenness in maintaining ecosystem process and function

The ecological consequences of species loss are widely studied, but represent an end point of environmental forcing that is not always realised. Changes in species evenness and the rank order of dominant species are more widespread responses to directional forcing. However, despite the repercussions for ecosystem functioning such changes have received little attention. Here, we experimentally assess how the rearrangement of species dominance structure within specific levels of evenness, rather than changes in species richness and composition, affect invertebrate particle reworking and burrow ventilation behaviour - important moderators of microbial-mediated remineralisation processes in benthic environments - and associated levels of sediment nutrient release. We find that the most dominant species exert a disproportionate influence on functioning at low levels of evenness, but that changes in biomass distribution and a change in emphasis in species-environmental interactions become more important in governing system functionality as evenness increases. Our study highlights the need to consider the functional significance of alterations to community attributes, rather than to solely focus on the attainment of particular levels of diversity when safeguarding biodiversity and ecosystems that provide essential services to society.

Resource-dependent attenuation of species interactions during bacterial succession

Bacterial communities are vital for many economically and ecologically important processes. The role of bacterial community composition in determining ecosystem functioning depends critically on interactions among bacterial taxa. Several studies have shown that, despite a predominance of negative interactions in communities, bacteria are able to display positive interactions given the appropriate evolutionary or ecological conditions. We were interested in how interspecific interactions develop over time in a naturalistic setting of low resource supply rates. We assembled aquatic bacterial communities in microcosms and assayed the productivity (respiration and growth) and substrate degradation while tracking community composition. The results demonstrated that while bacterial communities displayed strongly negative interactions during the early phase of colonisation and acclimatisation to novel biotic and abiotic factors, this antagonism declined over time towards a more neutral state. This was associated with a shift from use of labile substrates in early succession to use of recalcitrant substrates later in succession, confirming a crucial role of resource dynamics in linking interspecific interactions with ecosystem functioning.

Anthropogenic sources of underwater sound can modify how sediment-dwelling invertebrates mediate ecosystem properties

Coastal and shelf environments support high levels of biodiversity that are vital in mediating ecosystem processes, but they are also subject to noise associated with mounting levels of offshore human activity. This has the potential to alter the way in which species interact with their environment, compromising the mediation of important ecosystem properties. Here, we show that exposure to underwater broadband sound fields that resemble offshore shipping and construction activity can alter sediment-dwelling invertebrate contributions to fluid and particle transport--key processes in mediating benthic nutrient cycling. Despite high levels of intra-specific variability in physiological response, we find that changes in the behaviour of some functionally important species can be dependent on the class of broadband sound (continuous or impulsive). Our study provides evidence that exposing coastal environments to anthropogenic sound fields is likely to have much wider ecosystem consequences than are presently acknowledged.

The extent of functional redundancy changes as species' roles shift in different environments

Assessing the ecological impacts of environmental change requires knowledge of the relationship between biodiversity and ecosystem functioning. The exact nature of this relationship can differ considerably between ecosystems, with consequences for the efficacy of species diversity as a buffer against environmental change. Using a microbial model system, we show that the relationship can vary depending on environmental conditions. Shapes suggesting functional redundancy in one environment can change, suggesting functional differences in another environment. We find that this change is due to shifting species roles and interactions. Species that are functionally redundant in one environment may become pivotal in another. Thus, caution is advised in drawing conclusions about functional redundancy based on a single environmental situation. It also implies that species richness is important because it provides a pool of species with potentially relevant traits. These species may turn out to be essential performers or partners in new interspecific interactions after environmental change. Therefore, our results challenge the generality of functional redundancy.

Deciphering microbial community robustness through synthetic ecology and molecular systems synecology

Microbial ecosystems exhibit specific robustness attributes arising from the assembly and interaction networks of diverse, heterogeneous communities challenged by fluctuating environmental conditions. Synthetic ecology provides new insights into key biodiversity-stability relationships and robustness determinants of host-associated or environmental microbiomes. Driven by the advances of meta-omics technologies and bioinformatics, community-centered approaches (defined as molecular systems synecology) combined with the development of dynamic and mechanistic mathematical models make it possible to decipher and predict the outcomes of microbial ecosystems under disturbances. Beyond discriminating the normal operating range and natural, intrinsic dynamics of microbial processes from systems-level responses to environmental forcing, predictive modeling is poised to be integrated within prescriptive analytical frameworks and thus provide guidance in decision-making and proactive microbial resource management.

Long-term effects of warming and ocean acidification are modified by seasonal variation in species responses and environmental conditions

Warming of sea surface temperatures and alteration of ocean chemistry associated with anthropogenic increases in atmospheric carbon dioxide will have profound consequences for a broad range of species, but the potential for seasonal variation to modify species and ecosystem responses to these stressors has received little attention. Here, using the longest experiment to date (542 days), we investigate how the interactive effects of warming and ocean acidification affect the growth, behaviour and associated levels of ecosystem functioning (nutrient release) for a functionally important non-calcifying intertidal polychaete (Alitta virens) under seasonally changing conditions. We find that the effects of warming, ocean acidification and their interactions are not detectable in the short term, but manifest over time through changes in growth, bioturbation and bioirrigation behaviour that, in turn, affect nutrient generation. These changes are intimately linked to species responses to seasonal variations in environmental conditions (temperature and photoperiod) that, depending upon timing, can either exacerbate or buffer the long-term directional effects of climatic forcing. Taken together, our observations caution against over emphasizing the conclusions from short-term experiments and highlight the necessity to consider the temporal expression of complex system dynamics established over appropriate timescales when forecasting the likely ecological consequences of climatic forcing.

Organism-sediment interactions govern post-hypoxia recovery of ecosystem functioning

Hypoxia represents one of the major causes of biodiversity and ecosystem functioning loss for coastal waters. Since eutrophication-induced hypoxic events are becoming increasingly frequent and intense, understanding the response of ecosystems to hypoxia is of primary importance to understand and predict the stability of ecosystem functioning. Such ecological stability may greatly depend on the recovery patterns of communities and the return time of the system properties associated to these patterns. Here, we have examined how the reassembly of a benthic community contributed to the recovery of ecosystem functioning following experimentally-induced hypoxia in a tidal flat. We demonstrate that organism-sediment interactions that depend on organism size and relate to mobility traits and sediment reworking capacities are generally more important than recovering species richness to set the return time of the measured sediment processes and properties. Specifically, increasing macrofauna bioturbation potential during community reassembly significantly contributed to the recovery of sediment processes and properties such as denitrification, bedload sediment transport, primary production and deep pore water ammonium concentration. Such bioturbation potential was due to the replacement of the small-sized organisms that recolonised at early stages by large-sized bioturbating organisms, which had a disproportionately stronger influence on sediment. This study suggests that the complete recovery of organism-sediment interactions is a necessary condition for ecosystem functioning recovery, and that such process requires long periods after disturbance due to the slow growth of juveniles into adult stages involved in these interactions. Consequently, repeated episodes of disturbance at intervals smaller than the time needed for the system to fully recover organism-sediment interactions may greatly impair the resilience of ecosystem functioning.