Effects of antidepressant exposure on aquatic communities assessed by a combination of morphological identification, functional measurements, environmental DNA metabarcoding and bioassays.

The antidepressant fluoxetine is frequently detected in aquatic ecosystems, yet the effects on aquatic communities and ecosystems are still largely unknown. Therefore the aim of this study is to assess the effects of the long-term application of fluoxetine on key components of aquatic ecosystems including macroinvertebrate-, zooplankton-, phytoplankton- and microbial communities and organic matter decomposition by using traditional and non-traditional assessment methods. For this, we exposed 18 outdoor mesocosms (water volume of 1530 L and 10 cm of sediment) to five different concentrations of fluoxetine (0.2, 2, 20 and 200 µg/L) for eight weeks, followed by an eight-week recovery period. We quantified population and community effects by morphological identification, environmental DNA metabarcoding, in vitro and in vivo bioassays and measured organic matter decomposition as a measure of ecosystem functioning. We found effects of fluoxetine on bacterial, algal, zooplankton and macroinvertebrate communities and decomposition rates, mainly for the highest (200 µg/L) treatment. Treatment-related decreases in abundances were found for damselfly larvae (NOEC of 0.2 µg/L) and Sphaeriidae bivalves (NOEC of 20 µg/L), whereas Asellus aquaticus increased in abundance (NOEC <0.2 µg/L). Fluoxetine decreased photosynthetic activity and primary production of the suspended algae community. eDNA assessment provided additional insights by revealing that the algae belonging to the class Cryptophyceae and certain cyanobacteria taxa were the most negatively responding taxa to fluoxetine. Our results, together with results of others, suggest that fluoxetine can alter community structure and ecosystem functioning and that some impacts of fluoxetine on certain taxa can already be observed at environmentally realistic concentrations.

Assessing ecological responses to exposure to the antibiotic sulfamethoxazole in freshwater mesocosms.

Antibiotics are a contaminant class of worldwide concern as they are frequently detected in aquatic ecosystems. To better understand the impacts of antibiotics on aquatic ecosystems, we conducted an outdoor mesocosm experiment in which aquatic communities were exposed to different concentrations of the antibiotic sulfamethoxazole (0, 0.15, 1.5, 15 and 150 µg/L). These concentrations include mean (0.15 µg/L) and maximum detected concentrations (15 and 150 µg/L) in aquatic ecosystems worldwide. Sulfamethoxazole was applied once a week for eight consecutive weeks to 1530 L outdoor mesocosms in the Netherlands, followed by an eight-week recovery period. We evaluated phytoplankton-, bacterial- and invertebrate responses during and after sulfamethoxazole exposure and assessed impacts on organic matter decomposition. Contrary to our expectations, consistent treatment-related effects on algal and bacterial communities could not be demonstrated. In addition, sulfamethoxazole did not significantly affect zooplankton and macroinvertebrate communities. However, some effects on specific taxa were observed, with an increase in Mesostoma flatworm abundance (NOEC of <0.15 µg/L). In addition, eDNA analyses indicated negative impacts on the insects Odonata at a sulfamethoxazole concentration of 15 µg/L. Overall, environmentally relevant sulfamethoxazole concentration did not result in direct or indirect impairment of entire aquatic communities and ecological processes in our mesocosms. However, several specific macroinvertebrate taxa demonstrated significant (in)direct effects from sulfamethoxazole. Comparison of the results with the literature showed inconsistent results between studies using comparable, environmentally relevant, concentrations. Therefore, our study highlights the importance of testing the ecological impacts of pharmaceuticals (such as sulfamethoxazole) across multiple trophic levels spanning multiple aquatic communities, to fully understand its potential ecological threats.

Environmental DNA metabarcoding reveals comparable responses to agricultural stressors on different trophic levels of a freshwater community.

Freshwater habitats are under stress from agricultural land use, most notably the influx of neonicotinoid pesticides and increased nutrient pressure from fertilizer. Traditional studies investigating the effects of stressors on freshwater systems are often limited to a narrow range of taxa, depending heavily on morphological expertise. Additionally, disentanglement of multiple simultaneous stressors can be difficult in field studies, whereas controlled laboratory conditions do not accurately reflect natural conditions and food webs. To overcome these drawbacks, we investigated the impacts of two agricultural stressors (the neonicotinoid insecticide thiacloprid and fertilizer) in full-factorial design in a semi-natural research site, using environmental DNA sampling to study three different taxonomic groups representing three trophic levels: bacteria (decomposers), phytoplankton (primary producers), and chironomids (consumers). The results show considerable impact of both stressors across trophic levels, with an additive effect of fertilizer and thiacloprid on community composition at all levels. These findings suggest that agricultural stressors affect the entire food web, either directly or through cascade reactions. They are also consistent with morphological assessments that were performed in the same study site, even at a lower number of replicates. The study presented shows that the use of multimarker environmental DNA provides a more comprehensive assessment of stressor impacts across multiple trophic levels, at a higher taxonomic resolution than traditional surveys. Additionally, many putative novel bioindicators for both agricultural stressors were discovered. We encourage further investigations into stressors impacts at different trophic levels, which will lead to more effective monitoring and management of freshwater systems.

Organic matter reduces the amount of detectable environmental DNA in freshwater.

Environmental DNA (eDNA) is used for monitoring the occurrence of freshwater organisms. Various studies show a relation between the amount of eDNA detected and target organism abundance, thus providing a potential proxy for reconstructing population densities. However, environmental factors such as water temperature and microbial activity are known to affect the amount of eDNA present as well. In this study, we use controlled aquarium experiments using Gammarus pulex L. (Amphipoda) to investigate the relationship between the amount of detectable eDNA through time, pH, and levels of organic material. We found eDNA to degrade faster when organic material was added to the aquarium water, but that pH had no significant effect. We infer that eDNA contained inside cells and mitochondria is extra resilient against degradation, though this may not reflect actual presence of target species. These results indicate that, although estimation of population density might be possible using eDNA, measured eDNA concentration could, in the future, be corrected for local environmental conditions in order to ensure accurate comparisons.

Increased performance of DNA metabarcoding of macroinvertebrates by taxonomic sorting.

DNA-based identification through the use of metabarcoding has been proposed as the next step in the monitoring of biological communities, such as those assessed under the Water Framework Directive (WFD). Advances have been made in the field of metabarcoding, but challenges remain when using complex samples. Uneven biomass distributions, preferential amplification and reference database deficiencies can all lead to discrepancies between morphological and DNA-based taxa lists. The effects of different taxonomic groups on these issues remain understudied. By metabarcoding WFD monitoring samples, we analyzed six different taxonomic groups of freshwater organisms, both separately and combined. Identifications based on metabarcoding data were compared directly to morphological assessments performed under the WFD. The diversity of taxa for both morphological and DNA-based assessments was similar, although large differences were observed in some samples. The overlap between the two taxon lists was 56.8% on average across all taxa, and was highest for Crustacea, Heteroptera, and Coleoptera, and lowest for Annelida and Mollusca. Taxonomic sorting in six basic groups before DNA extraction and amplification improved taxon recovery by 46.5%. The impact on ecological quality ratio (EQR) scoring was considerable when replacing morphology with DNA-based identifications, but there was a high correlation when only replacing a single taxonomic group with molecular data. Different taxonomic groups provide their own challenges and benefits. Some groups might benefit from a more consistent and robust method of identification. Others present difficulties in molecular processing, due to uneven biomass distributions, large genetic diversity or shortcomings of the reference database. Sorting samples into basic taxonomic groups that require little taxonomic knowledge greatly improves the recovery of taxa with metabarcoding. Current standards for EQR monitoring may not be easily replaced completely with molecular strategies, but the effectiveness of molecular methods opens up the way for a paradigm shift in biomonitoring.

The effects of spatial and temporal replicate sampling on eDNA metabarcoding.

BACKGROUND: The heterogeneous nature of environmental DNA (eDNA) and its effects on species detection and community composition estimates has been highlighted in several studies in the past decades. Mostly in the context of spatial distribution over large areas, in fewer occasions looking at spatial distribution within a single body of water. Temporal variation of eDNA, similarly, has mostly been studied as seasonality, observing changes over large periods of time, and often only for small groups of organisms such as fish and amphibians. METHODS: We analyzed and compared small-scale spatial and temporal variation by sampling eDNA from two small, isolated dune lakes for 20 consecutive weeks. Metabarcoding was performed on the samples using generic COI primers. Molecular operational taxonomic unit (MOTUs) were used to assess dissimilarities between spatial and temporal replicates. RESULTS: Our results show large differences between samples taken within one lake at one point in time, but also expose the large differences between temporal replicates, even those taken only 1 week apart. Furthermore, between-site dissimilarities showed a linear correlation with time frame, indicating that between-site differences will be inflated when samples are taken over a period of time. We also assessed the effects of PCR replicates and processing strategies on general patterns of dissimilarity between samples. While more inclusive PCR replicate strategies lead to higher richness estimations, dissimilarity patterns between samples did not significantly change. CONCLUSIONS: We conclude that the dissimilarity of temporal replicates at a 1 week interval is comparable to that of spatial replicate samples. It increases, however, for larger time intervals, which suggests that population turnover effects can be stronger than community heterogeneity. Spatial replicates alone may not be enough for optimal recovery of taxonomic diversity, and cross-comparisons of different locations are susceptible to inflated dissimilarities when performed over larger time intervals. Many of the observed MOTUs could be classified as either phyto- or zooplankton, two groups that have gained traction in recent years as potential novel bio-indicator species. Our results, however, indicate that these groups might be susceptible to large community shifts in relatively short periods of time, highlighting the need to take temporal variations into consideration when assessing their usability as water quality indicators.

Incorporation of an invasive plant into a native insect herbivore food web.

The integration of invasive species into native food webs represent multifarious dynamics of ecological and evolutionary processes. We document incorporation of Prunus serotina (black cherry) into native insect food webs. We find that P. serotina harbours a herbivore community less dense but more diverse than its native relative, P. padus (bird cherry), with similar proportions of specialists and generalists. While herbivory on P. padus remained stable over the past century, that on P. serotina gradually doubled. We show that P. serotina may have evolved changes in investment in cyanogenic glycosides compared with its native range. In the leaf beetle Gonioctena quinquepunctata, recently shifted from native Sorbus aucuparia to P. serotina, we find divergent host preferences on Sorbus- versus Prunus-derived populations, and weak host-specific differentiation among 380 individuals genotyped for 119 SNP loci. We conclude that evolutionary processes may generate a specialized herbivore community on an invasive plant, allowing prognoses of reduced invasiveness over time. On the basis of the results presented here, we would like to caution that manual control might have the adverse effect of a slowing down of processes of adaptation, and a delay in the decline of the invasive character of P. serotina.

Evolution of endemism on a young tropical mountain.

Tropical mountains are hot spots of biodiversity and endemism, but the evolutionary origins of their unique biotas are poorly understood. In varying degrees, local and regional extinction, long-distance colonization, and local recruitment may all contribute to the exceptional character of these communities. Also, it is debated whether mountain endemics mostly originate from local lowland taxa, or from lineages that reach the mountain by long-range dispersal from cool localities elsewhere. Here we investigate the evolutionary routes to endemism by sampling an entire tropical mountain biota on the 4,095-metre-high Mount Kinabalu in Sabah, East Malaysia. We discover that most of its unique biodiversity is younger than the mountain itself (6 million years), and comprises a mix of immigrant pre-adapted lineages and descendants from local lowland ancestors, although substantial shifts from lower to higher vegetation zones in this latter group were rare. These insights could improve forecasts of the likelihood of extinction and 'evolutionary rescue' in montane biodiversity hot spots under climate change scenarios.

Pollination of Specklinia by nectar-feeding Drosophila: the first reported case of a deceptive syndrome employing aggregation pheromones in Orchidaceae.

BACKGROUND AND AIMS: The first documented observation of pollination in Pleurothallidinae was that of Endrés, who noticed that the 'viscid sepals' of Specklinia endotrachys were visited by a 'small fly'. Chase would later identify the visiting flies as being members of the genus Drosophila. This study documents and describes how species of the S. endotrachys complex are pollinated by different Drosophila species. METHODS: Specimens of Specklinia and Drosophila were collected in the field in Costa Rica and preserved in the JBL and L herbaria. Flies were photographed, filmed and observed for several days during a 2-year period and were identified by a combination of non-invasive DNA barcoding and anatomical surveys. Tissue samples of the sepals, petals and labellum of Specklinia species were observed and documented by SEM, LM and TEM. Electroantennogram experiments were carried out on Drosophila hydei using the known aggregation pheromones ethyl tiglate, methyl tiglate and isopropyl tiglate. Floral compounds were analysed by gas chromatography-mass spectometry using those same pheromones as standards. KEY RESULTS: Flowers of S. endotrachys, S. pfavii, S. remotiflora and S. spectabilis are visited and pollinated by several different but closely related Drosophila species. The flies are arrested by aggregation pheromones, including ethyl tiglate, methyl tiglate and isopropyl tiglate, released by the flowers, and to which at least D. hydei is very sensitive. Visible nectar drops on the adaxial surface of sepals are secreted by nectar-secreting stomata, encouraging male and female Drosophila to linger on the flowers for several hours at a time. The flies frequently show courtship behaviour, occasionally copulating. Several different Drosophila species can be found on a single Specklinia species. CONCLUSIONS: Species of the S. endotrachys group share a similar pollination syndrome. There seem to be no species-specific relationships between the orchids and the flies. It is not expected that Specklinia species will hybridize naturally as their populations do not overlap geographically. The combination of pheromone attraction and nectar feeding is likely to be a generalized pollination syndrome in Pleurothallidinae.

Cyberdiversity: improving the informatic value of diverse tropical arthropod inventories.

In an era of biodiversity crisis, arthropods have great potential to inform conservation assessment and test hypotheses about community assembly. This is because their relatively narrow geographic distributions and high diversity offer high-resolution data on landscape-scale patterns of biodiversity. However, a major impediment to the more widespread application of arthropod data to a range of scientific and policy questions is the poor state of modern arthropod taxonomy, especially in the tropics. Inventories of spiders and other megadiverse arthropods from tropical forests are dominated by undescribed species. Such studies typically organize their data using morphospecies codes, which make it difficult for data from independent inventories to be compared and combined. To combat this shortcoming, we offer cyberdiversity, an online community-based approach for reconciling results of independent inventory studies where current taxonomic knowledge is incomplete. Participating scientists can upload images and DNA barcode sequences to dedicated databases and submit occurrence data and links to a web site (www.digitalSpiders.org). Taxonomic determinations can be shared with a crowdsourcing comments feature, and researchers can discover specimens of interest available for loan and request aliquots of genomic DNA extract. To demonstrate the value of the cyberdiversity framework, we reconcile data from three rapid structured inventories of spiders conducted in Vietnam with an independent inventory (Doi Inthanon, Thailand) using online image libraries. Species richness and inventory completeness were assessed using non-parametric estimators. Community similarity was evaluated using a novel index based on the Jaccard replacing observed with estimated values to correct for unobserved species. We use a distance-decay framework to demonstrate a rudimentary model of landscape-scale changes in community composition that will become increasingly informative as additional inventories participate. With broader adoption of the cyberdiversity approach, networks of information-sharing taxonomists can more efficiently and effectively address taxonomic impediments while elucidating landscape scale patterns of biodiversity.

DNA barcoding of Dutch birds.

The mitochondrial cytochrome c oxidase subunit I (COI) can serve as a fast and accurate marker for the identification of animal species, and has been applied in a number of studies on birds. We here sequenced the COI gene for 387 individuals of 147 species of birds from the Netherlands, with 83 species being represented by > 2 sequences. The Netherlands occupies a small geographic area and 95% of all samples were collected within a 50 km radius from one another. The intraspecific divergences averaged 0.29% among this assemblage, but most values were lower; the interspecific divergences averaged 9.54%. In all, 95% of species were represented by a unique barcode, with 6 species of gulls and skua (Larus and Stercorarius) having at least one shared barcode. This is best explained by these species representing recent radiations with ongoing hybridization. In contrast, one species, the Lesser Whitethroat Sylvia curruca showed deep divergences, averaging 5.76% and up to 8.68% between individuals. These possibly represent two distinct taxa, S. curruca and S. blythi, both clearly separated in a haplotype network analysis. Our study adds to a growing body of DNA barcodes that have become available for birds, and shows that a DNA barcoding approach enables to identify known Dutch bird species with a very high resolution. In addition some species were flagged up for further detailed taxonomic investigation, illustrating that even in ornithologically well-known areas such as the Netherlands, more is to be learned about the birds that are present.

Which specimens from a museum collection will yield DNA barcodes? A time series study of spiders in alcohol.

We report initial results from an ongoing effort to build a library of DNA barcode sequences for Dutch spiders and investigate the utility of museum collections as a source of specimens for barcoding spiders. Source material for the library comes from a combination of specimens freshly collected in the field specifically for this project and museum specimens collected in the past. For the museum specimens, we focus on 31 species that have been frequently collected over the past several decades. A series of progressively older specimens representing these 31 species were selected for DNA barcoding. Based on the pattern of sequencing successes and failures, we find that smaller-bodied species expire before larger-bodied species as tissue sources for single-PCR standard DNA barcoding. Body size and age of oldest successful DNA barcode are significantly correlated after factoring out phylogenetic effects using independent contrasts analysis. We found some evidence that extracted DNA concentration is correlated with body size and inversely correlated with time since collection, but these relationships are neither strong nor consistent. DNA was extracted from all specimens using standard destructive techniques involving the removal and grinding of tissue. A subset of specimens was selected to evaluate nondestructive extraction. Nondestructive extractions significantly extended the DNA barcoding shelf life of museum specimens, especially small-bodied species, and yielded higher DNA concentrations compared to destructive extractions. All primary data are publically available through a Dryad archive and the Barcode of Life database.

Spider hosts (Arachnida, Araneae) and wasp parasitoids (Insecta, Hymenoptera, Ichneumonidae, Ephialtini) matched using DNA barcodes.

The study of parasitoids and their hosts suffers from a lack of reliable taxonomic data. We use a combination of morphological characters and DNA sequences to produce taxonomic determinations that can be verified with reference to specimens in an accessible collection and DNA barcode sequences posted to the Barcode of Life database (BOLD). We demonstrate that DNA can be successfully extracted from consumed host spiders and the shed pupal case of a wasp using non-destructive methods. We found Acrodactylaquadrisculpta to be a parasitoid of Tetragnathamontana; Zatypotapercontatoria and Zatypotabohemani both are parasitoids of Neottiurabimaculata. Zatypotaanomala is a parasitoid of an as yet unidentified host in the family Dictynidae, but the host species may be possible to identify in the future as the library of reference sequences on BOLD continues to grow. The study of parasitoids and their hosts traditionally requires specialized knowledge and techniques, and accumulating data is a slow process. DNA barcoding could allow more professional and amateur naturalists to contribute data to this field of study. A publication venue dedicated to aggregating datasets of all sizes online is well suited to this model of distributed science.