Assessing the relevance of DNA metabarcoding compared to morphological identification for lake phytoplankton monitoring.

Phytoplankton is a key biological group used to assess the ecological status of lakes. The classical monitoring approach relies on microscopic identification and counting of phytoplankton species, which is time-consuming and requires high taxonomic expertise. High-throughput sequencing, combined with metabarcoding, has recently demonstrated its potential as an alternative approach for plankton surveys. Several studies have confirmed the relevance of the diatom metabarcoding approach to calculate biotic indices based on species ecology. However, phytoplankton communities have not yet benefited from such validation. Here, by comparing the results obtained with the two methods (molecular and microscopic counting), we evaluated the relevance of metabarcoding approach for phytoplankton monitoring by considering different metrics: alpha diversity, taxonomic composition, community structure and a phytoplankton biotic index used to assess the trophic level of lakes. For this purpose, 55 samples were collected in four large alpine lakes (Aiguebelette, Annecy, Bourget, Geneva) during the year 2021. For each sample, a metabarcoding analysis based on two genetic markers (16S and 23S rRNA) was performed, in addition to the microscopic count. Regarding the trophic level of lakes, significant differences were found between index values obtained with the two approaches. The main hypothesis to explain these differences comes from the incompleteness, particularly at the species level, of the barcode reference library for the two genetic markers. It is therefore necessary to complete reference libraries for using such species-based biotic indices with metabarcoding data. Besides this, species richness and diversity were higher in the molecular inventories than in the microscopic ones. Moreover, despite differences in taxonomic composition of the floristic lists obtained by the two approaches, their community structures were similar. These results support the possibility of using metabarcoding for phytoplankton monitoring but in a different way. We suggest exploring alternative approaches to index development, such as a taxonomy-free approach.

Differential global distribution of marine picocyanobacteria gene clusters reveals distinct niche-related adaptive strategies.

The ever-increasing number of available microbial genomes and metagenomes provides new opportunities to investigate the links between niche partitioning and genome evolution in the ocean, especially for the abundant and ubiquitous marine picocyanobacteria Prochlorococcus and Synechococcus. Here, by combining metagenome analyses of the Tara Oceans dataset with comparative genomics, including phyletic patterns and genomic context of individual genes from 256 reference genomes, we show that picocyanobacterial communities thriving in different niches possess distinct gene repertoires. We also identify clusters of adjacent genes that display specific distribution patterns in the field (eCAGs) and are thus potentially involved in the same metabolic pathway and may have a key role in niche adaptation. Several eCAGs are likely involved in the uptake or incorporation of complex organic forms of nutrients, such as guanidine, cyanate, cyanide, pyrimidine, or phosphonates, which might be either directly used by cells, for example for the biosynthesis of proteins or DNA, or degraded to inorganic nitrogen and/or phosphorus forms. We also highlight the enrichment of eCAGs involved in polysaccharide capsule biosynthesis in Synechococcus populations thriving in both nitrogen- and phosphorus-depleted areas vs. low-iron (Fe) regions, suggesting that the complexes they encode may be too energy-consuming for picocyanobacteria thriving in the latter areas. In contrast, Prochlorococcus populations thriving in Fe-depleted areas specifically possess an alternative respiratory terminal oxidase, potentially involved in the reduction of Fe(III) to Fe(II). Altogether, this study provides insights into how phytoplankton communities populate oceanic ecosystems, which is relevant to understanding their capacity to respond to ongoing climate change.

Metal bioavailable contamination engages richness decline, species turnover but unchanged functional diversity of stream macroinvertebrates at the scale of a French region.

Freshwater ecosystems are the main source of water for sustaining life on earth, and the biodiversity they support is the main source of valuable goods and services for human populations. Despite growing recognition of the impairment of freshwater ecosystems by micropollutant contamination, different conceptual and methodological considerations can newly be addressed to improve our understanding of the ecological impact into these ecosystems. Here, we originally combined in situ ecotoxicology and community ecology concepts to unveil the mechanisms structuring macroinvertebrate communities along a regional contamination gradient. The novelty of our study lies in the use of an innovative biomonitoring approach (measurement of metal contents in caged crustaceans) allowing to quantify and compare on a regional scale the levels of bioavailable metal contamination to which stream communities are exposed. We were hence able to identify 23 streams presenting a significant gradient of bioavailable metal contamination within the same catchment area in the South West of France, from which we also obtained data on the composition of resident macroinvertebrate communities. Analyses of structural and functional integrity of communities revealed an unexpected decoupling between taxonomic and functional diversity of communities in response to bioavailable metal contamination. We show that despite the negative impact of bioavailable metal contamination exposure on taxonomic diversity (with an average species loss of 17% in contaminated streams), functional diversity is maintained through a process of non-random species replacement by functional redundant species at the regional scale. Such unanticipated findings call for a deeper characterization of metal-tolerant communities' ability to cope with environmental variability in multi-stressed ecosystems.

Presence of male mitochondria in somatic tissues and their functional importance at the whole animal level in the marine bivalve Arctica islandica.

Metazoans normally possess a single lineage of mitochondria inherited from the mother (♀-type mitochondria) while paternal mitochondria are absent or eliminated in fertilized eggs. In doubly uniparental inheritance (DUI), which is specific to the bivalve clade including the ocean quahog, Arctica islandica, ♂-type mitochondria are retained in male gonads and, in a few species, small proportions of ♂-type mitochondria co-exist with ♀-type in somatic tissues. To the best of our knowledge, we report, for the first time in metazoan, the natural occurrence of male and female individuals with exclusively ♂-type mitochondria in somatic tissues of the bivalve A. islandica. Mitochondrial genomes differ by ~5.5% at DNA sequence level. Exclusive presence of ♂-type mitochondria affects mitochondrial complexes partially encoded by mitochondrial genes and leads to a sharp drop in respiratory capacity. Through a combination of whole mitochondrial genome sequencing and molecular assays (gene presence and expression), we demonstrate that 1) 11% of individuals of an Icelandic population appear homoplasmic for ♂-type mitochondria in somatic tissues, 2) ♂-type mitochondrial genes are transcribed and 3) individuals with ♂-type mitochondria in somatic cells lose 30% of their wild-type respiratory capacity. This mitochondrial pattern in A. islandica is a special case of DUI, highlighted in individuals from both sexes with functional consequences at cellular and conceivably whole animal level.

Small freshwater ecosystems with dissimilar microbial communities exhibit similar temporal patterns.

Despite small freshwater ecosystems being biodiversity reservoirs and contributing significantly to greenhouse fluxes, their microbial communities remain largely understudied. Yet, microorganisms intervene in biogeochemical cycling and impact water quality. Because of their small size, these ecosystems are in principle more sensitive to disturbances, seasonal variation and pluri-annual climate change. However, how microbial community composition varies over space and time, and whether archaeal, bacterial and microbial eukaryote communities behave similarly remain unanswered. Here, we aim to unravel the composition and intra/interannual temporal dynamic patterns for archaea, bacteria and microbial eukaryotes in a set of small freshwater ecosystems. We monitored archaeal and bacterial community composition during 24 consecutive months in four ponds and one brook from northwestern France by 16S rRNA gene amplicon sequencing (microbial eukaryotes were previously investigated for the same systems). Unexpectedly for oxic environments, bacterial Candidate Phyla Radiation (CPR) were highly diverse and locally abundant. Our results suggest that microbial community structure is mainly driven by environmental conditions acting over space (ecosystems) and time (seasons). A low proportion of operational taxonomic units (OTUs) (<1%) was shared by the five ecosystems despite their geographical proximity (2-9 km away), making microbial communities almost unique in each ecosystem and highlighting the strong selective influence of local environmental conditions. Marked and similar seasonality patterns were observed for archaea, bacteria and microbial eukaryotes in all ecosystems despite strong turnovers of rare OTUs. Over the 2-year survey, microbial community composition varied despite relatively stable environmental parameters. This suggests that biotic associations play an important role in interannual community assembly.

How diatom-, invertebrate- and fish-based diagnostic tools can support the ecological assessment of rivers in a multi-pressure context: Temporal trends over the past two decades in France.

The degradation of aquatic ecosystems, induced by worldwide intensification in the use of both land and aquatic resources, has highlighted the critical need for innovative methods allowing an objective quantification and ranking of anthropogenic pressure effects on aquatic organisms. Such diagnostic tools have a great potential for defining robust management responses to anthropogenic pressures. Our objective was to explore how the outputs of three diagnostic tools (based on benthic diatoms, macroinvertebrates and fishes) could be combined to (i) disentangle the temporal effects of multiple pressures over two decades and (ii) provide policy-relevant information for stream managers and decision makers. The diagnostic tools estimated, using taxonomy- and trait-based metrics, the impairment probabilities of biotic assemblages over time by different pressure categories, describing the alteration of water quality, hydromorphology and land use related to anthropogenic activities, in French streams (number of sites = 312). The main result shows that a large proportion of the time series exhibited no significant temporal patterns over the two decades (61.5% to 87.8%, depending on the used tests). Among time series exhibiting significant change, positive trends in impairment probabilities (i.e., degradation) were less frequent than negative ones, indicating a modest improvement in water quality at national scale over the study period. However, trends can be substantially different according to hydroecoregion and pressure category. The three biological compartments displayed convergent temporal responses according to the pressure category and regional context (e.g., lowland plains vs. mountains, pristine vs. agricultural regions). Altogether, this study proposes a unifying approach to integrate a vast amount of information in a single ecological diagnosis using an unparalleled database on natural and anthropized environments. Strengthening the synthesis of biological information provided by various biological compartments should be a priority before implementing evidence-based sustainable conservation and restoration actions.

How to quantify the links between bioavailable contamination in watercourses and pressures of anthropogenic land cover, contamination sources and hydromorphology at multiple scales?

Active biomonitoring permits the quantification of biological exposure to chemicals through measurements of bioavailable concentrations in biota and biological markers of toxicity in organisms. It enables respective comparison of the levels of contamination between sites and sampling campaigns. Caged gammarids are recently proposed as relevant probes for measuring bioavailable contamination in freshwater systems. The purpose of the present study was to develop a multi-pressure and multiscale approach, considering metallic contamination levels (from data based on active biomonitoring) as a response to pressures (combination of individual stressors). These pressures were anthropogenic land cover, industry density, wastewater treatment plant density, pressures on stream hydromorphological functioning, riverside vegetation and bioavailability factors. A dataset combining active biomonitoring and potentially related pressures was established at the French national scale, with 196 samplings from 2009 to 2016. The links between pressures and metallic contamination were defined and modelled via structural equation modeling (more specifically partial least squares - path modeling). The model enabled the understanding of the respective influences of pressures on metallic bioconcentration in caged sentinel organisms. Beyond validating the local influence of industries and wastewater treatment plants on metallic contamination, this model showed a complementary effect of driving forces of anthropogenic land cover (leading to human activities). It also quantified a significant influence of pressures on stream hydromorphological functioning, presence of vegetation and physico-chemical parameters on metal bioconcentration. This hierarchical multi-pressure approach could serve as a concept on how pressures and contamination (assessed by active biomonitoring) can be connected. Its future application will enable better understanding of environmental pressures leading to contamination in freshwater ecosystems.

Investigating microbial associations from sequencing survey data with co-correspondence analysis.

Microbial communities, which drive major ecosystem functions, consist of a wide range of interacting species. Understanding how microbial communities are structured and the processes underlying this is crucial to interpreting ecosystem responses to global change but is challenging as microbial interactions cannot usually be directly observed. Multiple efforts are currently focused to combine next-generation sequencing (NGS) techniques with refined statistical analysis (e.g., network analysis, multivariate analysis) to characterize the structures of microbial communities. However, most of these approaches consider a single table of sequencing data measured for several samples. Technological advances now make it possible to collect NGS data on different taxonomic groups simultaneously for the same samples, allowing us to analyse a pair of tables. Here, an analytical framework based on co-correspondence analysis (CoCA) is proposed to study the distributions, assemblages and interactions between two microbial communities. We show the ability of this approach to highlight the relationships between two microbial communities, using two data sets exhibiting various types of interactions. CoCA identified strong association patterns between autotrophic and heterotrophic microbial eukaryote assemblages, on the one hand, and between microalgae and viruses, on the other. We demonstrate also how CoCA can be used, complementary to network analysis, to reorder co-occurrence networks and thus investigate the presence of patterns in ecological networks.

Multisubstance Indicators Based on Caged Gammarus Bioaccumulation Reveal the Influence of Chemical Contamination on Stream Macroinvertebrate Abundances across France.

Most anthropogenic stressors affecting freshwater systems are qualitatively known. However, the quantitative assessment of contaminant exposure and effects to aquatic communities is still difficult, limiting the understanding of consequences on aquatic ecosystem functioning and the implementation of effective management plans. Here, multisubstance indicators based on caged gammarid bioaccumulated contamination data are proposed (for metals and persistent organic pollutants, POPs) to map the bioavailable contamination level of freshwater ecosystems at a large spatial scale. We assessed the ability of these indicators to highlight the relationships between chemical exposure gradients and alteration in the abundance of macroinvertebrate populations on a data set of 218 watercourses distributed throughout France. We identified spatial regional heterogeneities in the levels of bioavailable contamination of metals (18 compounds) and POPs (43 compounds). Besides this, a degradation of Gammaridae, Ephemeridae, and Hydrobiidae densities with increasing levels of metal contamination are identified relative to Baetidae, Chironomidae, and Hydropsychidae. We show here that active biomonitoring allows the establishment of multisubstance indicators of bioavailable contamination, which reliably quantify chemical exposure gradients in freshwater ecosystems. Our ability to identify species-specific responses to chemical exposure gradients demonstrates the promising possibility to further decipher the effects of chemical contamination on macroinvertebrate assemblages through this type of indicator.

Local forcings affect lake zooplankton vulnerability and response to climate warming.

While considerable insights on the ecological consequences of climate change have been gained from studies conducted on remote lakes, little has been done on lakes under direct human exposure. Ecosystem vulnerability and responses to climate warming might yet largely depend on the ecological state and thus on local anthropogenic pressures. We tested this hypothesis through a paleolimnological approach on three temperate large lakes submitted to rather similar climate warming but varying intensities of analogous local forcings (changes in nutrient inputs and fisheries management practices). Changes in the structure of the cladoceran community were considered as revealing for alterations, over the time, of the pelagic food web. Trajectories of the cladoceran communities were compared among the three study lakes (Lakes Geneva, Bourget, and Annecy) over the last 70-150 years. Generalized additive models were used to develop a hierarchical understanding of the respective roles of local stressors and climate warming in structuring cladoceran communities. The cladoceran communities were not equally affected by climate warming between lakes. In Lake Annecy, which is the most nutrient-limited, the cladoceran community was essentially controlled by local stressors, with very limited impact of climate. In contrast, the more nutrient-loaded Lakes Geneva and Bourget were more sensitive to climate warming, although the magnitude of their responses and the pathways under which climate warming affected the communities varied between the two lakes. Finally, our results demonstrated that lake vulnerability and responses to climate warming are modulated by lake trophic status but can also be altered by fisheries management practices through changes in fish predation pressure.