Mathew Seymour

Environmental DNA metabarcoding: transforming how we survey animal and plant communities

The genomic revolution has fundamentally changed how we survey biodiversity on earth. High‐throughput sequencing (“HTS”) platforms now enable the rapid sequencing of DNA from diverse kinds of environmental samples (termed “environmental DNA” or “eDNA”). Coupling HTS with our ability to associate sequences from eDNA with a taxonomic name is called “eDNA metabarcoding” and offers a powerful molecular tool capable of noninvasively surveying species richness from many ecosystems. Here, we review the use of eDNA metabarcoding for surveying animal and plant richness, and the challenges in using eDNA approaches to estimate relative abundance. We highlight eDNA applications in freshwater, marine and terrestrial environments, and in this broad context, we distill what is known about the ability of different eDNA sample types to approximate richness in space and across time. We provide guiding questions for study design and discuss the eDNA metabarcoding workflow with a focus on primers and library preparation methods. We additionally discuss important criteria for consideration of bioinformatic filtering of data sets, with recommendations for increasing transparency. Finally, looking to the future, we discuss emerging applications of eDNA metabarcoding in ecology, conservation, invasion biology, biomonitoring, and how eDNA metabarcoding can empower citizen science and biodiversity education.

Annual time-series analysis of aqueous eDNA reveals ecologically relevant dynamics of lake ecosystem biodiversity

The use of environmental DNA (eDNA) in biodiversity assessments offers a step-change in sensitivity, throughput and simultaneous measures of ecosystem diversity and function. There remains, however, a need to examine eDNA persistence in the wild through simultaneous temporal measures of eDNA and biota. Here, we use metabarcoding of two markers of different lengths, derived from an annual time series of aqueous lake eDNA to examine temporal shifts in ecosystem biodiversity and in an ecologically important group of macroinvertebrates (Diptera: Chironomidae). The analyses allow different levels of detection and validation of taxon richness and community composition (β-diversity) through time, with shorter eDNA fragments dominating the eDNA community. Comparisons between eDNA, community DNA, taxonomy and UK species abundance data further show significant relationships between diversity estimates derived across the disparate methodologies. Our results reveal the temporal dynamics of eDNA and validate the utility of eDNA metabarcoding for tracking seasonal diversity at the ecosystem scale.

Active colonization dynamics and diversity patterns are influenced by dendritic network connectivity and species interactions

Habitat network connectivity influences colonization dynamics, species invasions, and biodiversity patterns. Recent theoretical work suggests dendritic networks, such as those found in rivers, alter expectations regarding colonization and dispersal dynamics compared with other network types. As many native and non-native species are spreading along river networks, this may have important ecological implications. However, experimental studies testing the effects of network structure on colonization and diversity patterns are scarce. Up to now, experimental studies have only considered networks where sites are connected with small corridors, or dispersal was experimentally controlled, which eliminates possible effects of species interactions on colonization dynamics. Here, we tested the effect of network connectivity and species interactions on colonization dynamics using continuous linear and dendritic (i.e., river-like) networks, which allow for active dispersal. We used a set of six protist species and one rotifer species in linear and dendritic microcosm networks. At the start of the experiment, we introduced species, either singularly or as a community within the networks. Species subsequently actively colonized the networks. We periodically measured densities of species throughout the networks over 2 weeks to track community dynamics, colonization, and diversity patterns. We found that colonization of dendritic networks was faster compared with colonization of linear networks, which resulted in higher local mean species richness in dendritic networks. Initially, community similarity was also greater in dendritic networks compared with linear networks, but this effect vanished over time. The presence of species interactions increased community evenness over time, compared with extrapolations from single-species setups. Our experimental findings confirm previous theoretical work and show that network connectivity, species-specific dispersal ability, and species interactions greatly influence the dispersal and colonization of dendritic networks. We argue that these factors need to be considered in empirical studies, where effects of network connectivity on colonization patterns have been largely underestimated.

Experimental evidence for strong stabilizing forces at high functional diversity of aquatic microbial communities

Unveiling the mechanisms that promote coexistence in biological communities is a fundamental problem in ecology. Stable coexistence of many species is commonly observed in natural communities. Most of these natural communities, however, are composed of species from multiple trophic and functional groups, while theory and experiments on coexistence have been focusing on functionally similar species. Here, we investigated how functional diversity affects the stability of species coexistence and productivity in multispecies communities by characterizing experimentally all pairwise species interactions in a pool of 11 species of eukaryotes (10 protists and one rotifer) belonging to three different functional groups. Species within the same functional group showed stronger competitive interactions compared to among-functional group interactions. This often led to competitive exclusion between species that had higher functional relatedness, but only at low levels of species richness. Communities with higher functional diversity resulted in increased species coexistence and community biomass production. Our experimental findings and the results of a stochastic model tailored to the experimental interaction matrix suggest the emergence of strong stabilizing forces when species from different functional groups interact in a homogeneous environment. By combining theoretical analysis with experiments we could also disentangle the relationship between species richness and functional diversity, showing that functional diversity per se is a crucial driver of productivity and stability in multispecies community.

Inferring species interactions in ecological communities: a comparison of methods at different levels of complexity

1. Natural communities commonly contain many different species and functional groups, and multiple types of species interactions act simultaneously, such as competition, predation, commensalism or mutualism. However, experimental and theoretical investigations have generally been limited by focusing on one type of interaction at a time or by a lack of a common methodological and conceptual approach to measure species interactions. 2. We compared four methods to measure and express species interactions. These approaches are, with increasing degree of model complexity, an extinction-based model, a relative yield model and two generalized Lotka-Volterra (LV) models. All four approaches have been individually applied in different fields of community ecology, but rarely integrated. We provide an overview of the definitions, assumptions and data needed for the specific methods and apply them to empirical data by experimentally deriving the interaction matrices among 11 protist and rotifer species, belonging to three functional groups. Furthermore, we compare their advantages and limitations to predict multispecies community dynamics and ecosystem functioning. 3. The relative yield method is, in terms of final biomass production, the best method in predicting the 11-species community dynamics from the pairwise competition experiments. The LV model, which is considering equilibrium among the species, suffers from experimental constraints given the strict equilibrium assumption, and this may be rarely satisfied in ecological communities. 4. We show how simulations of a LV stochastic community model, derived from an empirical interaction matrix, can be used to predict multispecies community dynamics across multiple functional groups. 5.Our work unites available tools to measure species interactions under one framework. This improves our ability to make management-oriented predictions of species coexistence/extinction and to compare ecosystem processes across study systems.

Connectivity in a pond system influences migration and genetic structure in threespine stickleback

Neutral genetic structure of natural populations is primarily influenced by migration (the movement of individuals and, subsequently, their genes) and drift (the statistical chance of losing genetic diversity over time). Migration between populations is influenced by several factors, including individual behavior, physical barriers, and environmental heterogeneity among populations. However, drift is expected to be stronger in populations with low immigration rate and small effective population size. With the technological advancement in geological information systems and spatial analysis tools, landscape genetics now allows the development of realistic migration models and increased insight to important processes influencing diversity of natural populations. In this study, we investigated the relationship between landscape connectivity and genetic distance of threespine stickleback (Gasterosteus aculeatus) inhabiting a pond complex in Belgjarskógur, Northeast Iceland. We used two landscape genetic approaches (i.e., least-cost-path and isolation-by-resistance) and asked whether gene flow, as measured by genetic distance, was more strongly associated with Euclidean distance (isolation-by-distance) or with landscape connectivity provided by areas prone to flooding (as indicated by Carex sp. cover)? We found substantial genetic structure across the study area, with pairwise genetic distances among populations (DPS) ranging from 0.118 to 0.488. Genetic distances among populations were more strongly correlated with least-cost-path and isolation-by-resistance than with Euclidean distance, whereas the relative contribution of isolation-by-resistance and Euclidian distance could not be disentangled. These results indicate that migration among stickleback populations occurs via periodically flooded areas. Overall, this study highlights the importance of transient landscape elements influencing migration and genetic structure of populations at small spatial scales.

Acidity promotes degradation of multi-species environmental DNA in lotic mesocosms

Accurate quantification of biodiversity is fundamental to understanding ecosystem function and for environmental assessment. Molecular methods using environmental DNA (eDNA) offer a non-invasive, rapid, and cost-effective alternative to traditional biodiversity assessments, which require high levels of expertise. While eDNA analyses are increasingly being utilized, there remains considerable uncertainty regarding the dynamics of multispecies eDNA, especially in variable systems such as rivers. Here, we utilize four sets of upland stream mesocosms, across an acid–base gradient, to assess the temporal and environmental degradation of multispecies eDNA. Sampling included water column and biofilm sampling over time with eDNA quantified using qPCR. Our findings show that the persistence of lotic multispecies eDNA, sampled from water and biofilm, decays to non-detectable levels within 2 days and that acidic environments accelerate the degradation process. Collectively, the results provide the basis for a predictive framework for the relationship between lotic eDNA degradation dynamics in spatio-temporally dynamic river ecosystems.Mathew Seymour et al. investigate the persistence of environmental DNA (eDNA) in river systems in environments of varying pH. Using four sets of upland stream mesocosms, they find that eDNA degrades to non-detectable levels within two days and this degradation is accelerated in acidic environments.

Performance of amplicon and shotgun sequencing for accurate biomass estimation in invertebrate community samples

New applications of DNA and RNA sequencing are expanding the field of biodiversity discovery and ecological monitoring, yet questions remain regarding precision and efficiency. Due to primer bias, the ability of metabarcoding to accurately depict biomass of different taxa from bulk communities remains unclear, while PCR‐free whole mitochondrial genome (mitogenome) sequencing may provide a more reliable alternative. Here, we used a set of documented mock communities comprising 13 species of freshwater macroinvertebrates of estimated individual biomass, to compare the detection efficiency of COI metabarcoding (three different amplicons) and shotgun mitogenome sequencing. Additionally, we used individual COI barcoding and de novo mitochondrial genome sequencing, to provide reference sequences for OTU assignment and metagenome mapping (mitogenome skimming), respectively. We found that, even though both methods occasionally failed to recover very low abundance species, metabarcoding was less consistent, by failing to recover some species with higher abundances, probably due to primer bias. Shotgun sequencing results provided highly significant correlations between read number and biomass in all but one species. Conversely, the read–biomass relationships obtained from metabarcoding varied across amplicons. Specifically, we found significant relationships for eight of 13 (amplicons B1FR‐450 bp, FF130R‐130 bp) or four of 13 (amplicon FFFR, 658 bp) species. Combining the results of all three COI amplicons (multiamplicon approach) improved the read–biomass correlations for some of the species. Overall, mitogenomic sequencing yielded more informative predictions of biomass content from bulk macroinvertebrate communities than metabarcoding. However, for large‐scale ecological studies, metabarcoding currently remains the most commonly used approach for diversity assessment.

Marine ecology: Genetics from a drop in the ocean

Analysis of environmental DNA (eDNA) extracted from just 30 litres of seawater from the Arabian Gulf provides genetic insights into populations of the largest fish in the world.

Dynamic species classification of microorganisms across time, abiotic and biotic environments—A sliding window approach

The development of video-based monitoring methods allows for rapid, dynamic and accurate monitoring of individuals or communities, compared to slower traditional methods, with far reaching ecological and evolutionary applications. Large amounts of data are generated using video-based methods, which can be effectively processed using machine learning (ML) algorithms into meaningful ecological information. ML uses user defined classes (e.g. species), derived from a subset (i.e. training data) of video-observed quantitative features (e.g. phenotypic variation), to infer classes in subsequent observations. However, phenotypic variation often changes due to environmental conditions, which may lead to poor classification, if environmentally induced variation in phenotypes is not accounted for. Here we describe a framework for classifying species under changing environmental conditions based on the random forest classification. A sliding window approach was developed that restricts temporal and environmentally conditions to improve the classification. We tested our approach by applying the classification framework to experimental data. The experiment used a set of six ciliate species to monitor changes in community structure and behavior over hundreds of generations, in dozens of species combinations and across a temperature gradient. Differences in biotic and abiotic conditions caused simplistic classification approaches to be unsuccessful. In contrast, the sliding window approach allowed classification to be highly successful, as phenotypic differences driven by environmental change, could be captured by the classifier. Importantly, classification using the random forest algorithm showed comparable success when validated against traditional, slower, manual identification. Our framework allows for reliable classification in dynamic environments, and may help to improve strategies for long-term monitoring of species in changing environments. Our classification pipeline can be applied in fields assessing species community dynamics, such as eco-toxicology, ecology and evolutionary ecology.