Zooplankton community structure along a pollution gradient at fine geographical scales in river ecosystems: The importance of species sorting over dispersal

The release of anthropogenic pollution into freshwater ecosystems has largely transformed biodiversity and its geographical distribution patterns globally. However, for many communities including ecologically crucial ones such as zooplankton, it is largely unknown how different communities respond to environmental pollution. Collectively, dispersal and species sorting are two competing processes in determining the structure and geographical distribution of zooplankton communities in running water ecosystems such as rivers. At fine geographical scales, dispersal is usually considered as the dominant factor; however, the relative role of species sorting has not been evaluated well, mainly because significant environmental gradients rarely exist along continuously flowing rivers. The Chaobai River in northern China represents a rare system, where a significant environmental gradient exists at fine scales. Here, we employed high‐throughput sequencing to characterize complex zooplankton communities collected from the Chaobai River, and tested the relative roles of dispersal and species sorting in determining zooplankton community structure along the pollution gradient. Our results showed distinct patterns of zooplankton communities along the environmental gradient, and chemical pollutant‐related factors such as total phosphorus and chlorophyll‐a were identified as the major drivers for the observed patterns. Further partial redundancy analyses showed that species sorting overrode the effect of dispersal to shape local zooplankton community structure. Thus, our results reject the dispersal hypothesis and support the concept that species sorting caused by local pollution can largely determine the zooplankton community structure when significant environmental gradients exist at fine geographical scales in highly polluted running water ecosystems.     

Not all are free‐living: high‐throughput DNA metabarcoding reveals a diverse community of protists parasitizing soil metazoa

Protists, the most diverse eukaryotes, are largely considered to be free‐living bacterivores, but vast numbers of taxa are known to parasitize plants or animals. High‐throughput sequencing (HTS) approaches now commonly replace cultivation‐based approaches in studying soil protists, but insights into common biases associated with this method are limited to aquatic taxa and samples. We created a mock community of common free‐living soil protists (amoebae, flagellates, ciliates), extracted DNA and amplified it in the presence of metazoan DNA using 454 HTS. We aimed at evaluating whether HTS quantitatively reveals true relative abundances of soil protists and at investigating whether the expected protist community structure is altered by the co‐amplification of metazoan‐associated protist taxa. Indeed, HTS revealed fundamentally different protist communities from those expected. Ciliate sequences were highly over‐represented, while those of most amoebae and flagellates were under‐represented or totally absent. These results underpin the biases introduced by HTS that prevent reliable quantitative estimations of free‐living protist communities. Furthermore, we detected a wide range of nonadded protist taxa probably introduced along with metazoan DNA, which altered the protist community structure. Among those, 20 taxa most closely resembled parasitic, often pathogenic taxa. Therewith, we provide the first HTS data in support of classical observational studies that showed that potential protist parasites are hosted by soil metazoa. Taken together, profound differences in amplification success between protist taxa and an inevitable co‐extraction of protist taxa parasitizing soil metazoa obscure the true diversity of free‐living soil protist communities.     

The risks of using molecular biodiversity data for incidental detection of species of concern

Incidental detection of species of concern (e.g., invasive species, pathogens, threatened and endangered species) during biodiversity assessments based on high‐throughput DNA sequencing holds significant risks in the absence of rigorous, fit‐for‐purpose data quality and reporting standards. Molecular biodiversity data are predominantly collected for ecological studies and thus are generated to common quality assurance standards. However, the detection of certain species of concern in these data would likely elicit interest from end users working in biosecurity or other surveillance contexts (e.g., pathogen detection in health‐related fields), for which more stringent quality control standards are essential to ensure that data are suitable for informing decision‐making and can withstand legal or political challenges. We suggest here that data quality and reporting criteria are urgently needed to enable clear identification of those studies that may be appropriately applied to surveillance contexts. In the interim, more pointed disclaimers on uncertainties associated with the detection and identification of species of concern may be warranted in published studies. This is not only to ensure the utility of molecular biodiversity data for consumers, but also to protect data generators from uncritical and potentially ill‐advised application of their science in decision‐making.     

Advancing biodiversity assessments with environmental DNA: Long‐read technologies help reveal the drivers of Amazonian fungal diversity

Fungi are a key component of tropical biodiversity. However, due to their inconspicuous and largely subterranean nature, they are usually neglected in biodiversity inventories. The goal of this study was to identify the key determinants of fungal richness, community composition, and turnover in tropical rainforests. We tested specifically for the effect of soil properties, habitat, and locality in Amazonia. For these analyses, we used high‐throughput sequencing data of short and long reads of fungal DNA present in soil and organic litter samples, combining existing and novel genomic data. Habitat type (phytophysiognomy) emerges as the strongest factor explaining fungal community composition. Naturally open areas—campinas—are the richest habitat overall. Soil properties have different effects depending on the soil layer (litter or mineral soil) and the choice of genetic marker. We suggest that campinas could be a neglected hotspot of fungal diversity. An underlying cause for their rich diversity may be the overall low soil fertility, which increases the reliance on biotic interactions essential for nutrient absorption in these environments, notably ectomycorrhizal fungi–plant associations. Our results highlight the advantages of using both short and long DNA reads produced through high‐throughput sequencing to characterize fungal diversity. While short reads can suffice for diversity and community comparison, long reads add taxonomic precision and have the potential to reveal population diversity.     

Changes in ectomycorrhizal fungal community composition and declining diversity along a 2‐million‐year soil chronosequence

Ectomycorrhizal (ECM) fungal communities covary with host plant communities along soil fertility gradients, yet it is unclear whether this reflects changes in host composition, fungal edaphic specialization or priority effects during fungal community establishment. We grew two co‐occurring ECM plant species (to control for host identity) in soils collected along a 2‐million‐year chronosequence representing a strong soil fertility gradient and used soil manipulations to disentangle the effects of edaphic properties from those due to fungal inoculum. Ectomycorrhizal fungal community composition changed and richness declined with increasing soil age; these changes were linked to pedogenesis‐driven shifts in edaphic properties, particularly pH and resin‐exchangeable and organic phosphorus. However, when differences in inoculum potential or soil abiotic properties among soil ages were removed while host identity was held constant, differences in ECM fungal communities and richness among chronosequence stages disappeared. Our results show that ECM fungal communities strongly vary during long‐term ecosystem development, even within the same hosts. However, these changes could not be attributed to short‐term fungal edaphic specialization or differences in fungal inoculum (i.e. density and composition) alone. Rather, they must reflect longer‐term ecosystem‐level feedback between soil, vegetation and ECM fungi during pedogenesis.     

Toward an ecoregion scale evaluation of eDNA metabarcoding primers: A case study for the freshwater fish biodiversity of the Murray–Darling Basin (Australia)

High‐throughput sequencing of environmental DNA (i.e., eDNA metabarcoding) has become an increasingly popular method for monitoring aquatic biodiversity. At present, such analyses require target‐specific primers to amplify DNA barcodes from co‐occurring species, and this initial amplification can introduce biases. Understanding the performance of different primers is thus recommended prior to undertaking any metabarcoding initiative. While multiple software programs are available to evaluate metabarcoding primers, all programs have their own strengths and weaknesses. Therefore, a robust in silico workflow for the evaluation of metabarcoding primers will benefit from the use of multiple programs. Furthermore, geographic differences in species biodiversity are likely to influence the performance of metabarcoding primers and further complicate the evaluation process. Here, an in silico workflow is presented that can be used to evaluate the performance of metabarcoding primers on an ecoregion scale. This workflow was used to evaluate the performance of published and newly developed eDNA metabarcoding primers for the freshwater fish biodiversity of the Murray–Darling Basin (Australia). To validate the in silico workflow, a subset of the primers, including one newly designed primer pair, were used in metabarcoding analyses of an artificial DNA community and eDNA samples. The results show that the in silico workflow allows for a robust evaluation of metabarcoding primers and can reveal important trade‐offs that need to be considered when selecting the most suitable primer. Additionally, a new primer pair was described and validated that allows for more robust taxonomic assignments and is less influenced by primer biases compared to commonly used fish metabarcoding primers.     

Terra incognita of soil biodiversity: unseen invasions under our feet

Whilst cartographers of the 19th century endeavoured to chart the last unknown lands, the great challenge for biologists in the 21st century is to fill the gaps on the biodiversity map of the Earth. And one of the largest gaps concerns the biodiversity of soils, a terra incognita right under our feet. The study of soil biodiversity, and particularly the complex communities of small invertebrates, has suffered from a severe ‘taxonomic impediment’ (Decaëns 2010 ) leading to great uncertainties about total species richness, phylogenetic diversity, geographical structure, temporal dynamics of soil organisms, and consequently about their role on ecosystem function (Bardgett & van der Putten 2014 ). However, the revolution in high‐throughput sequencing is now revealing the hidden biodiversity of the soil with unprecedented detail (e.g. Arribas et al. 2016 ). In a noteworthy from the Cover article in this issue of Molecular Ecology, Cicconardi et al. ( 2017 ) apply these new tools to study soil communities of Collembola in three distant oceanic islands of volcanic origin, obtaining a striking result: only 38 of 70 species (54%) are exclusively found in a single island, with the remaining shared among islands or with other distant regions, suggesting a massive recent introduction of soil species, whose impact is entirely unknown.     

Scaling up: A guide to high‐throughput genomic approaches for biodiversity analysis

The purpose of this review is to present the most common and emerging DNA‐based methods used to generate data for biodiversity and biomonitoring studies. As environmental assessment and monitoring programmes may require biodiversity information at multiple levels, we pay particular attention to the DNA metabarcoding method and discuss a number of bioinformatic tools and considerations for producing DNA‐based indicators using operational taxonomic units (OTUs), taxa at a variety of ranks and community composition. By developing the capacity to harness the advantages provided by the newest technologies, investigators can “scale up” by increasing the number of samples and replicates processed, the frequency of sampling over time and space, and even the depth of sampling such as by sequencing more reads per sample or more markers per sample. The ability to scale up is made possible by the reduced hands‐on time and cost per sample provided by the newest kits, platforms and software tools. Results gleaned from broad‐scale monitoring will provide opportunities to address key scientific questions linked to biodiversity and its dynamics across time and space as well as being more relevant for policymakers, enabling science‐based decision‐making, and provide a greater socio‐economic impact. As genomic approaches are continually evolving, we provide this guide to methods used in biodiversity genomics.     

Sequence capture using PCR‐generated probes: a cost‐effective method of targeted high‐throughput sequencing for nonmodel organisms

Recent advances in high‐throughput sequencing library preparation and subgenomic enrichment methods have opened new avenues for population genetics and phylogenetics of nonmodel organisms. To multiplex large numbers of indexed samples while sequencing predominantly orthologous, targeted regions of the genome, we propose modifications to an existing, in‐solution capture that utilizes PCR products as target probes to enrich library pools for the genomic subset of interest. The sequence capture using PCR‐generated probes (SCPP) protocol requires no specialized equipment, is highly flexible and significantly reduces experimental costs for projects where a modest scale of genetic data is optimal (25–100 genomic loci). Our alterations enable application of this method across a wider phylogenetic range of taxa and result in higher capture efficiencies and coverage at each locus. Efficient and consistent capture over multiple SCPP experiments and at various phylogenetic distances is demonstrated, extending the utility of this method to both phylogeographic and phylogenomic studies.     

PipeCraft: Flexible open‐source toolkit for bioinformatics analysis of custom high‐throughput amplicon sequencing data

High‐throughput sequencing methods have become a routine analysis tool in environmental sciences as well as in public and private sector. These methods provide vast amount of data, which need to be analysed in several steps. Although the bioinformatics may be applied using several public tools, many analytical pipelines allow too few options for the optimal analysis for more complicated or customized designs. Here, we introduce PipeCraft, a flexible and handy bioinformatics pipeline with a user‐friendly graphical interface that links several public tools for analysing amplicon sequencing data. Users are able to customize the pipeline by selecting the most suitable tools and options to process raw sequences from Illumina, Pacific Biosciences, Ion Torrent and Roche 454 sequencing platforms. We described the design and options of PipeCraft and evaluated its performance by analysing the data sets from three different sequencing platforms. We demonstrated that PipeCraft is able to process large data sets within 24 hr. The graphical user interface and the automated links between various bioinformatics tools enable easy customization of the workflow. All analytical steps and options are recorded in log files and are easily traceable.     

Predicting fine‐scale distributions of peripheral aquatic species in headwater streams

Headwater species and peripheral populations that occupy habitat at the edge of a species range may hold an increased conservation value to managers due to their potential to maximize intraspecies diversity and species' adaptive capabilities in the context of rapid environmental change. The southern Appalachian Mountains are the southern extent of the geographic range of native Salvelinus fontinalis and naturalized Oncorhynchus mykiss and Salmo trutta in eastern North America. We predicted distributions of these peripheral, headwater wild trout populations at a fine scale to serve as a planning and management tool for resource managers to maximize resistance and resilience of these populations in the face of anthropogenic stressors. We developed correlative logistic regression models to predict occurrence of brook trout, rainbow trout, and brown trout for every interconfluence stream reach in the study area. A stream network was generated to capture a more consistent representation of headwater streams. Each of the final models had four significant metrics in common: stream order, fragmentation, precipitation, and land cover. Strahler stream order was found to be the most influential variable in two of the three final models and the second most influential variable in the other model. Greater than 70% presence accuracy was achieved for all three models. The underrepresentation of headwater streams in commonly used hydrography datasets is an important consideration that warrants close examination when forecasting headwater species distributions and range estimates. Additionally, it appears that a relative watershed position metric (e.g., stream order) is an important surrogate variable (even when elevation is included) for biotic interactions across the landscape in areas where headwater species distributions are influenced by topographical gradients.     

Emergence of nutrient co‐limitation through movement in stoichiometric meta‐ecosystems

Evidence that ecosystems and primary producers are limited in their productivity by multiple nutrients has caused the traditional nutrient limitation framework to include multiple limiting nutrients. The models built to mimic these responses have invoked local mechanisms at the level of the primary producers. In this paper, we explore an alternative explanation for the emergence of co‐limitation by developing a simple, stoichiometrically explicit meta‐ecosystem model with two limiting nutrients, autotrophs and herbivores. Our results show that differences in movement rates for the nutrients, autotrophs and herbivores can allow for nutrient co‐limitation in biomass response to emerge despite no local mechanisms of nutrient co‐limitation. Furthermore, our results provide an explanation to why autotrophs show positive growth responses to nutrients despite ‘nominal’ top‐down control by herbivores. These results suggest that spatial processes can be mechanisms for nutrient co‐limitation at local and regional scales, and can help explain anomalous results in the co‐limitation literature.     

Next‐generation sequencing to inventory taxonomic diversity in eukaryotic communities: a test for freshwater diatoms

The recent emergence of barcoding approaches coupled to those of next‐generation sequencing (NGS) has raised new perspectives for studying environmental communities. In this framework, we tested the possibility to derive accurate inventories of diatom communities from pyrosequencing outputs with an available DNA reference library. We used three molecular markers targeting the nuclear, chloroplast and mitochondrial genomes (SSU rDNA, rbcL and cox1) and three samples of a mock community composed of 30 known diatom strains belonging to 21 species. In the goal to detect methodological biases, one sample was constituted directly from pooled cultures, whereas the others consisted of pooled PCR products. The NGS reads obtained by pyrosequencing (Roche 454) were compared first to a DNA reference library including the sequences of all the species used to constitute the mock community, and second to a complete DNA reference library with a larger taxonomic coverage. A stringent taxonomic assignation gave inventories that were compared to the real one. We detected biases due to DNA extraction and PCR amplification that resulted in false‐negative detection. Conversely, pyrosequencing errors appeared to generate false positives, especially in case of closely allied species. The taxonomic coverage of DNA reference libraries appears to be the most crucial factor, together with marker polymorphism which is essential to identify taxa at the species level. RbcL offers a high resolving power together with a large DNA reference library. Although needing further optimization, pyrosequencing is suitable for identifying diatom assemblages and may find applications in the field of freshwater biomonitoring.     

Helping decision making for reliable and cost‐effective 2b‐RAD sequencing and genotyping analyses in non‐model species

High‐throughput sequencing has revolutionized population and conservation genetics. RAD sequencing methods, such as 2b‐RAD, can be used on species lacking a reference genome. However, transferring protocols across taxa can potentially lead to poor results. We tested two different IIB enzymes (AlfI and CspCI) on two species with different genome sizes (the loggerhead turtle Caretta caretta and the sharpsnout seabream Diplodus puntazzo) to build a set of guidelines to improve 2b‐RAD protocols on non‐model organisms while optimising costs. Good results were obtained even with degraded samples, showing the value of 2b‐RAD in studies with poor DNA quality. However, library quality was found to be a critical parameter on the number of reads and loci obtained for genotyping. Resampling analyses with different number of reads per individual showed a trade‐off between number of loci and number of reads per sample. The resulting accumulation curves can be used as a tool to calculate the number of sequences per individual needed to reach a mean depth ≥20 reads to acquire good genotyping results. Finally, we demonstrated that selective‐base ligation does not affect genomic differentiation between individuals, indicating that this technique can be used in species with large genome sizes to adjust the number of loci to the study scope, to reduce sequencing costs and to maintain suitable sequencing depth for a reliable genotyping without compromising the results. Here, we provide a set of guidelines to improve 2b‐RAD protocols on non‐model organisms with different genome sizes, helping decision‐making for a reliable and cost‐effective genotyping.     

Fluctuations in water level and the dynamics of zooplankton: a data‐driven modelling approach

1. The zooplankton in Lake Kinneret (Israel) have undergone large fluctuations in recent decades, which have been linked to both biotic and abiotic processes. 2. By applying a data‐driven modelling approach to a long‐term database, and focusing on key abiotic (lake‐level change) and biotic (prey abundance) variables, we attempted to identify the possible factors impacting the lake’s zooplankton community. 3. We hypothesised that changes in the predatory zooplankton (adult cyclopoids) assemblage are driven by changes in lake level during years of large changes. We further postulated that lake‐level changes would have a similar impact on the herbivorous zooplankton (cladocerans and cyclopoid copepodites) but an opposite effect on the microzooplankton. In the years of moderate changes to lake level, however, the abundance of predatory zooplankton would determine the size of the herbivore and microzooplankton populations rather than their food sources, that is, top‐down control. 4. The resulting decision trees supported the hypotheses stressing the importance of the annual rate of lake‐level change in shaping the zooplankton community in the lake. In addition, and in contrast to expectations, bottom‐up processes seem to play a role in determining zooplankton abundance.     

Local scale processes drive long‐term change in biodiversity of sandy beach ecosystems

Evaluating impacts to biodiversity requires ecologically informed comparisons over sufficient time spans. The vulnerability of coastal ecosystems to anthropogenic and climate change‐related impacts makes them potentially valuable indicators of biodiversity change. To evaluate multidecadal change in biodiversity, we compared results from intertidal surveys of 13 sandy beaches conducted in the 1970s and 2009–11 along 500 km of coast (California, USA). Using a novel extrapolation approach to adjust species richness for sampling effort allowed us to address data gaps and has promise for application to other data‐limited biodiversity comparisons. Long‐term changes in species richness varied in direction and magnitude among beaches and with human impacts but showed no regional patterns. Observed long‐term changes in richness differed markedly among functional groups of intertidal invertebrates. At the majority (77%) of beaches, changes in richness were most evident for wrack‐associated invertebrates suggesting they have disproportionate vulnerability to impacts. Reduced diversity of this group was consistent with long‐term habitat loss from erosion and sea level rise at one beach. Wrack‐associated species richness declined over time at impacted beaches (beach fill and grooming), despite observed increases in overall intertidal richness. In contrast richness of these taxa increased at more than half (53%) of the beaches including two beaches recovering from decades of off‐road vehicle impacts. Over more than three decades, our results suggest that local scale processes exerted a stronger influence on intertidal biodiversity on beaches than regional processes and highlight the role of human impacts for local spatial scales. Our results illustrate how comparisons of overall biodiversity may mask ecologically important changes and stress the value of evaluating biodiversity change in the context of functional groups. The long‐term loss of wrack‐associated species, a key component of sandy beach ecosystems, documented here represents a significant threat to the biodiversity and function of coastal ecosystems.