Refocusing the microbial rare biosphere concept through a functional lens.

The influential concept of the rare biosphere in microbial ecology has underscored the importance of taxa occurring at low abundances yet potentially playing key roles in communities and ecosystems. Here, we refocus the concept of rare biosphere through a functional trait-based lens and provide a framework to characterize microbial functional rarity, a combination of numerical scarcity across space or time and trait distinctiveness. We demonstrate how this novel interpretation of the rare biosphere, rooted in microbial functions, can enhance our mechanistic understanding of microbial community structure. It also sheds light on functionally distinct microbes, directing conservation efforts towards taxa harboring rare yet ecologically crucial functions.

Warming effects on lizard gut microbiome depend on habitat connectivity.

Climate warming and landscape fragmentation are both factors well known to threaten biodiversity and to generate species responses and adaptation. However, the impact of warming and fragmentation interplay on organismal responses remains largely under-explored, especially when it comes to gut symbionts, which may play a key role in essential host functions and traits by extending its functional and genetic repertoire. Here, we experimentally examined the combined effects of climate warming and habitat connectivity on the gut bacterial communities of the common lizard (Zootoca vivipara) over three years. While the strength of effects varied over the years, we found that a 2°C warmer climate decreases lizard gut microbiome diversity in isolated habitats. However, enabling connectivity among habitats with warmer and cooler climates offset or even reversed warming effects. The warming effects and the association between host dispersal behaviour and microbiome diversity appear to be a potential driver of this interplay. This study suggests that preserving habitat connectivity will play a key role in mitigating climate change impacts, including the diversity of the gut microbiome, and calls for more studies combining multiple anthropogenic stressors when predicting the persistence of species and communities through global changes.

Beyond reaction norms: the temporal dynamics of phenotypic plasticity.

Phenotypic plasticity can allow organisms to cope with environmental changes. Although reaction norms are commonly used to quantify plasticity along gradients of environmental conditions, they often miss the temporal dynamics of phenotypic change, especially the speed at which it occurs. Here, we argue that studying the rate of phenotypic plasticity is a crucial step to quantify and understand its adaptiveness. Iteratively measuring plastic traits allows us to describe the actual dynamics of phenotypic changes and avoid quantifying reaction norms at times that do not truly reflect the organism's capacity for plasticity. Integrating the temporal component in how we describe, quantify, and conceptualise phenotypic plasticity can change our understanding of its diversity, evolution, and consequences.

Molecular ecology of microbiomes in the wild: Common pitfalls, methodological advances and future directions.

The study of microbiomes across organisms and environments has become a prominent focus in molecular ecology. This perspective article explores common challenges, methodological advancements, and future directions in the field. Key research areas include understanding the drivers of microbiome community assembly, linking microbiome composition to host genetics, exploring microbial functions, transience and spatial partitioning, and disentangling non-bacterial components of the microbiome. Methodological advancements, such as quantifying absolute abundances, sequencing complete genomes, and utilizing novel statistical approaches, are also useful tools for understanding complex microbial diversity patterns. Our aims are to encourage robust practices in microbiome studies and inspire researchers to explore the next frontier of this rapidly changing field.

A robust approach to estimate relative phytoplankton cell abundances from metagenomes.

Phytoplankton account for >45% of global primary production, and have an enormous impact on aquatic food webs and on the entire Earth System. Their members are found among prokaryotes (cyanobacteria) and multiple eukaryotic lineages containing chloroplasts. Genetic surveys of phytoplankton communities generally consist of PCR amplification of bacterial (16S), nuclear (18S) and/or chloroplastic (16S) rRNA marker genes from DNA extracted from environmental samples. However, our appreciation of phytoplankton abundance or biomass is limited by PCR-amplification biases, rRNA gene copy number variations across taxa, and the fact that rRNA genes do not provide insights into metabolic traits such as photosynthesis. Here, we targeted the photosynthetic gene psbO from metagenomes to circumvent these limitations: the method is PCR-free, and the gene is universally and exclusively present in photosynthetic prokaryotes and eukaryotes, mainly in one copy per genome. We applied and validated this new strategy with the size-fractionated marine samples collected by Tara Oceans, and showed improved correlations with flow cytometry and microscopy than when based on rRNA genes. Furthermore, we revealed unexpected features of the ecology of these ecosystems, such as the high abundance of picocyanobacterial aggregates and symbionts in the ocean, and the decrease in relative abundance of phototrophs towards the larger size classes of marine dinoflagellates. To facilitate the incorporation of psbO in molecular-based surveys, we compiled a curated database of >18,000 unique sequences. Overall, psbO appears to be a promising new gene marker for molecular-based evaluations of entire phytoplankton communities.

Toward global integration of biodiversity big data: a harmonized metabarcode data generation module for terrestrial arthropods.

Metazoan metabarcoding is emerging as an essential strategy for inventorying biodiversity, with diverse projects currently generating massive quantities of community-level data. The potential for integrating across such data sets offers new opportunities to better understand biodiversity and how it might respond to global change. However, large-scale syntheses may be compromised if metabarcoding workflows differ from each other. There are ongoing efforts to improve standardization for the reporting of inventory data. However, harmonization at the stage of generating metabarcode data has yet to be addressed. A modular framework for harmonized data generation offers a pathway to navigate the complex structure of terrestrial metazoan biodiversity. Here, through our collective expertise as practitioners, method developers, and researchers leading metabarcoding initiatives to inventory terrestrial biodiversity, we seek to initiate a harmonized framework for metabarcode data generation, with a terrestrial arthropod module. We develop an initial set of submodules covering the 5 main steps of metabarcode data generation: (i) sample acquisition; (ii) sample processing; (iii) DNA extraction; (iv) polymerase chain reaction amplification, library preparation, and sequencing; and (v) DNA sequence and metadata deposition, providing a backbone for a terrestrial arthropod module. To achieve this, we (i) identified key points for harmonization, (ii) reviewed the current state of the art, and (iii) distilled existing knowledge within submodules, thus promoting best practice by providing guidelines and recommendations to reduce the universe of methodological options. We advocate the adoption and further development of the terrestrial arthropod module. We further encourage the development of modules for other biodiversity fractions as an essential step toward large-scale biodiversity synthesis through harmonization.

Best practices in metabarcoding of fungi: From experimental design to results.

The development of high-throughput sequencing (HTS) technologies has greatly improved our capacity to identify fungi and unveil their ecological roles across a variety of ecosystems. Here we provide an overview of current best practices in metabarcoding analysis of fungal communities, from experimental design through molecular and computational analyses. By reanalysing published data sets, we demonstrate that operational taxonomic units (OTUs) outperform amplified sequence variants (ASVs) in recovering fungal diversity, a finding that is particularly evident for long markers. Additionally, analysis of the full-length ITS region allows more accurate taxonomic placement of fungi and other eukaryotes compared to the ITS2 subregion. Finally, we show that specific methods for compositional data analyses provide more reliable estimates of shifts in community structure. We conclude that metabarcoding analyses of fungi are especially promising for integrating fungi into the full microbiome and broader ecosystem functioning context, recovery of novel fungal lineages and ancient organisms as well as barcoding of old specimens including type material.

Decay of similarity across tropical forest communities: integrating spatial distance with soil nutrients.

Understanding the mechanisms that drive the change of biotic assemblages over space and time is the main quest of community ecology. Assessing the relative importance of dispersal and environmental species selection in a range of organismic sizes and motilities has been a fruitful strategy. A consensus for whether spatial and environmental distances operate similarly across spatial scales and taxa, however, has yet to emerge. We used censuses of four major groups of organisms (soil bacteria, fungi, ground insects, and trees) at two observation scales (1-m2 sampling point vs. 2,500-m2 plots) in a topographically standardized sampling design replicated in two tropical rainforests with contrasting relationships between spatial distance and nutrient availability. We modeled the decay of assemblage similarity for each taxon set and site to assess the relative contributions of spatial distance and nutrient availability distance. Then, we evaluated the potentially structuring effect of tree composition over all other taxa. The similarity of nutrient content in the litter and topsoil had a stronger and more consistent selective effect than did dispersal limitation, particularly for bacteria, fungi, and trees at the plot level. Ground insects, the only group assessed with the capacity of active dispersal, had the highest species turnover and the flattest nonsignificant distance-decay relationship, suggesting that neither dispersal limitation nor nutrient availability were fundamental drivers of their community assembly at this scale of analysis. Only the fungal communities at one of our study sites were clearly coordinated with tree composition. The spatial distance at the smallest scale was more important than nutrient selection for the bacteria, fungi, and insects. The lower initial similarity and the moderate variation in composition identified by these distance-decay models, however, suggested that the effects of stochastic sampling were important at this smaller spatial scale. Our results highlight the importance of nutrients as one of the main environmental drivers of rainforest communities irrespective of organismic or propagule size and how the overriding effect of the analytical scale influences the interpretation, leading to the perception of greater importance of dispersal limitation and ecological drift over selection associated with environmental niches at decreasing observation scales.

Compendium of 530 metagenome-assembled bacterial and archaeal genomes from the polar Arctic Ocean.

The role of the Arctic Ocean ecosystem in climate regulation may depend on the responses of marine microorganisms to environmental change. We applied genome-resolved metagenomics to 41 Arctic seawater samples, collected at various depths in different seasons during the Tara Oceans Polar Circle expedition, to evaluate the ecology, metabolic potential and activity of resident bacteria and archaea. We assembled 530 metagenome-assembled genomes (MAGs) to form the Arctic MAGs catalogue comprising 526 species. A total of 441 MAGs belonged to species that have not previously been reported and 299 genomes showed an exclusively polar distribution. Most Arctic MAGs have large genomes and the potential for fast generation times, both of which may enable adaptation to a copiotrophic lifestyle in nutrient-rich waters. We identified 38 habitat generalists and 111 specialists in the Arctic Ocean. We also found a general prevalence of 14 mixotrophs, while chemolithoautotrophs were mostly present in the mesopelagic layer during spring and autumn. We revealed 62 MAGs classified as key Arctic species, found only in the Arctic Ocean, showing the highest gene expression values and predicted to have habitat-specific traits. The Artic MAGs catalogue will inform our understanding of polar microorganisms that drive global biogeochemical cycles.

Cascading effects of moth outbreaks on subarctic soil food webs.

The increasing severity and frequency of natural disturbances requires a better understanding of their effects on all compartments of biodiversity. In Northern Fennoscandia, recent large-scale moth outbreaks have led to an abrupt change in plant communities from birch forests dominated by dwarf shrubs to grass-dominated systems. However, the indirect effects on the belowground compartment remained unclear. Here, we combined eDNA surveys of multiple trophic groups with network analyses to demonstrate that moth defoliation has far-reaching consequences on soil food webs. Following this disturbance, diversity and relative abundance of certain trophic groups declined (e.g., ectomycorrhizal fungi), while many others expanded (e.g., bacterivores and omnivores) making soil food webs more diverse and structurally different. Overall, the direct and indirect consequences of moth outbreaks increased belowground diversity at different trophic levels. Our results highlight that a holistic view of ecosystems improves our understanding of cascading effects of major disturbances on soil food webs.

Connecting high-throughput biodiversity inventories: Opportunities for a site-based genomic framework for global integration and synthesis.

High-throughput sequencing (HTS) is increasingly being used for the characterization and monitoring of biodiversity. If applied in a structured way, across broad geographical scales, it offers the potential for a much deeper understanding of global biodiversity through the integration of massive quantities of molecular inventory data generated independently at local, regional and global scales. The universality, reliability and efficiency of HTS data can potentially facilitate the seamless linking of data among species assemblages from different sites, at different hierarchical levels of diversity, for any taxonomic group and regardless of prior taxonomic knowledge. However, collective international efforts are required to optimally exploit the potential of site-based HTS data for global integration and synthesis, efforts that at present are limited to the microbial domain. To contribute to the development of an analogous strategy for the nonmicrobial terrestrial domain, an international symposium entitled "Next Generation Biodiversity Monitoring" was held in November 2019 in Nicosia (Cyprus). The symposium brought together evolutionary geneticists, ecologists and biodiversity scientists involved in diverse regional and global initiatives using HTS as a core tool for biodiversity assessment. In this review, we summarize the consensus that emerged from the 3-day symposium. We converged on the opinion that an effective terrestrial Genomic Observatories network for global biodiversity integration and synthesis should be spatially led and strategically united under the umbrella of the metabarcoding approach. Subsequently, we outline an HTS-based strategy to collectively build an integrative framework for site-based biodiversity data generation.

Disentangling the assembly mechanisms of ant cuticular bacterial communities of two Amazonian ant species sharing a common arboreal nest.

Bacteria living on the cuticle of ants are generally studied for their protective role against pathogens, especially in the clade of fungus-growing ants. However, little is known regarding the diversity of cuticular bacteria in other ant host species, as well as the mechanisms leading to the composition of these communities. Here, we used 16S rRNA gene amplicon sequencing to study the influence of host species, species interactions and the pool of bacteria from the environment on the assembly of cuticular bacterial communities on two phylogenetically distant Amazonian ant species that frequently nest together inside the roots system of epiphytic plants, Camponotus femoratus and Crematogaster levior. Our results show that (a) the vast majority of the bacterial community on the cuticle is shared with the nest, suggesting that most bacteria on the cuticle are acquired through environmental acquisition, (b) 5.2% and 2.0% of operational taxonomic units (OTUs) are respectively specific to Ca. femoratus and Cr. levior, probably representing their respective core cuticular bacterial community, and (c) 3.6% of OTUs are shared between the two ant species. Additionally, mass spectrometry metabolomics analysis of metabolites on the cuticle of ants, which excludes the detection of cuticular hydrocarbons produced by the host, were conducted to evaluate correlations among bacterial OTUs and m/z ion mass. Although some positive and negative correlations are found, the cuticular chemical composition was weakly species-specific, suggesting that cuticular bacterial communities are prominently environmentally acquired. Overall, our results suggest the environment is the dominant source of bacteria found on the cuticle of ants.

Latent Dirichlet Allocation reveals spatial and taxonomic structure in a DNA-based census of soil biodiversity from a tropical forest.

High-throughput sequencing of amplicons from environmental DNA samples permits rapid, standardized and comprehensive biodiversity assessments. However, retrieving and interpreting the structure of such data sets requires efficient methods for dimensionality reduction. Latent Dirichlet Allocation (LDA) can be used to decompose environmental DNA samples into overlapping assemblages of co-occurring taxa. It is a flexible model-based method adapted to uneven sample sizes and to large and sparse data sets. Here, we compare LDA performance on abundance and occurrence data, and we quantify the robustness of the LDA decomposition by measuring its stability with respect to the algorithm's initialization. We then apply LDA to a survey of 1,131 soil DNA samples that were collected in a 12-ha plot of primary tropical forest and amplified using standard primers for bacteria, protists, fungi and metazoans. The analysis reveals that bacteria, protists and fungi exhibit a strong spatial structure, which matches the topographical features of the plot, while metazoans do not, confirming that microbial diversity is primarily controlled by environmental variation at the studied scale. We conclude that LDA is a sensitive, robust and computationally efficient method to detect and interpret the structure of large DNA-based biodiversity data sets. We finally discuss the possible future applications of this approach for the study of biodiversity.

Global Trends in Marine Plankton Diversity across Kingdoms of Life.

The ocean is home to myriad small planktonic organisms that underpin the functioning of marine ecosystems. However, their spatial patterns of diversity and the underlying drivers remain poorly known, precluding projections of their responses to global changes. Here we investigate the latitudinal gradients and global predictors of plankton diversity across archaea, bacteria, eukaryotes, and major virus clades using both molecular and imaging data from Tara Oceans. We show a decline of diversity for most planktonic groups toward the poles, mainly driven by decreasing ocean temperatures. Projections into the future suggest that severe warming of the surface ocean by the end of the 21st century could lead to tropicalization of the diversity of most planktonic groups in temperate and polar regions. These changes may have multiple consequences for marine ecosystem functioning and services and are expected to be particularly significant in key areas for carbon sequestration, fisheries, and marine conservation. VIDEO ABSTRACT.

Altered trophic interactions in warming climates: consequences for predator diet breadth and fitness.

Species interactions are central in predicting the impairment of biodiversity with climate change. Trophic interactions may be altered through climate-dependent changes in either predator food preferences or prey communities. Yet, climate change impacts on predator diet remain surprisingly poorly understood. We experimentally studied the consequences of 2°C warmer climatic conditions on the trophic niche of a generalist lizard predator. We used a system of semi-natural mesocosms housing a variety of invertebrate species and in which climatic conditions were manipulated. Lizards in warmer climatic conditions ate at a greater predatory to phytophagous invertebrate ratio and had smaller individual dietary breadths. These shifts mainly arose from direct impacts of climate on lizard diets rather than from changes in prey communities. Dietary changes were associated with negative changes in fitness-related traits (body condition, gut microbiota) and survival. We demonstrate that climate change alters trophic interactions through top-predator dietary shifts, which might disrupt eco-evolutionary dynamics.

The relative importance of ecological drivers of arbuscular mycorrhizal fungal distribution varies with taxon phylogenetic resolution.

The phylogenetic depth at which arbuscular mycorrhizal (AM) fungi harbor a coherent ecological niche is unknown, which has consequences for operational taxonomic unit (OTU) delineation from sequence data and the study of their biogeography. We tested how changes in AM fungi community composition across habitats (beta diversity) vary with OTU phylogenetic resolution. We inferred exact sequence variants (ESVs) to resolve phylotypes at resolutions finer than provided by traditional sequence clustering and analyzed beta diversity profiles up to order-level sequence clusters. At the ESV level, we detected the environmental predictors revealed with traditional OTUs or at higher genetic distances. However, the correlation between environmental predictors and community turnover steeply increased at a genetic distance of c. 0.03 substitutions per site. Furthermore, we observed a turnover of either closely or distantly related taxa (respectively at or above 0.03 substitutions per site) along different environmental gradients. This study suggests that different axes of AM fungal ecological niche are conserved at different phylogenetic depths. Delineating AM fungal phylotypes using DNA sequences should screen different phylogenetic resolutions to better elucidate the factors that shape communities and predict the fate of AM symbioses in a changing environment.

One-locus-several-primers: A strategy to improve the taxonomic and haplotypic coverage in diet metabarcoding studies.

In diet metabarcoding analyses, insufficient taxonomic coverage of PCR primer sets generates false negatives that may dramatically distort biodiversity estimates. In this paper, we investigated the taxonomic coverage and complementarity of three cytochrome c oxidase subunit I gene (COI) primer sets based on in silico analyses and we conducted an in vivo evaluation using fecal and spider web samples from different invertivores, environments, and geographic locations. Our results underline the lack of predictability of both the coverage and complementarity of individual primer sets: (a) sharp discrepancies exist observed between in silico and in vivo analyses (to the detriment of in silico analyses); (b) both coverage and complementarity depend greatly on the predator and on the taxonomic level at which preys are considered; (c) primer sets' complementarity is the greatest at fine taxonomic levels (molecular operational taxonomic units [MOTUs] and variants). We then formalized the "one-locus-several-primer-sets" (OLSP) strategy, that is, the use of several primer sets that target the same locus (here the first part of the COI gene) and the same group of taxa (here invertebrates). The proximal aim of the OLSP strategy is to minimize false negatives by increasing total coverage through multiple primer sets. We illustrate that the OLSP strategy is especially relevant from this perspective since distinct variants within the same MOTUs were not equally detected across all primer sets. Furthermore, the OLSP strategy produces largely overlapping and comparable sequences, which cannot be achieved when targeting different loci. This facilitates the use of haplotypic diversity information contained within metabarcoding datasets, for example, for phylogeography and finer analyses of prey-predator interactions.

Body size determines soil community assembly in a tropical forest.

Tropical forests shelter an unparalleled biological diversity. The relative influence of environmental selection (i.e., abiotic conditions, biotic interactions) and stochastic-distance-dependent neutral processes (i.e., demography, dispersal) in shaping communities has been extensively studied for various organisms, but has rarely been explored across a large range of body sizes, in particular in soil environments. We built a detailed census of the whole soil biota in a 12-ha tropical forest plot using soil DNA metabarcoding. We show that the distribution of 19 taxonomic groups (ranging from microbes to mesofauna) is primarily stochastic, suggesting that neutral processes are prominent drivers of the assembly of these communities at this scale. We also identify aluminium, topography and plant species identity as weak, yet significant drivers of soil richness and community composition of bacteria, protists and to a lesser extent fungi. Finally, we show that body size, which determines the scale at which an organism perceives its environment, predicted the community assembly across taxonomic groups, with soil mesofauna assemblages being more stochastic than microbial ones. These results suggest that the relative contribution of neutral processes and environmental selection to community assembly directly depends on body size. Body size is hence an important determinant of community assembly rules at the scale of the ecological community in tropical soils and should be accounted for in spatial models of tropical soil food webs.