Genomic basis for an informed conservation management of Pelophylax water frogs in Luxembourg.

Genetic identification methods have become increasingly important for species that are difficult to identify in the field. A case in point is Pelophylax water frogs. While their morphological determination is highly complex, they include species protected under EU law and some that are classified as invasive. Additionally, genetic data can provide insights into their complex breeding systems, which may or may not involve the reproductive dependency of one species on another. Here, we generate baseline data for water frog monitoring in Luxembourg. We applied a countrywide sampling approach and used SNPs generated by ddRAD sequencing to identify individuals and infer the breeding systems present in the country. We found Pelophylax lessonae and P. kl. esculentus throughout Luxembourg, mostly living in syntopy. In general, a reproductive dependency of P. kl. esculentus on P. lessonae (L-E system) was revealed. Besides this general system, we detected triploid P. kl. esculentus in six ponds. This indicates a modified L-E system with reproductive dependency of the triploids on the diploid P. kl. esculentus. The invasive P. cf. bedriagae was detected in three ponds in southern Luxembourg, with evidence for hybridization with native water frogs. In addition to the ddRAD data, we tested a simple genetic method for future monitoring based on the MND1 marker. It showed in almost all cases, an identical species identification as the ddRAD data and was successfully applied to DNA extracts from mouth swabs. Combining this method with our baseline data will enable informed choices for the protection of native water frog species in Luxembourg.

Individual small in-stream barriers contribute little to strong local population genetic structure five strictly aquatic macroinvertebrate taxa.

Water flow in river networks is frequently regulated by man-made in-stream barriers. These obstacles can hinder dispersal of aquatic organisms and isolate populations leading to the loss of genetic diversity. Although millions of small in-stream barriers exist worldwide, their impact on dispersal of macroinvertebrates remains unclear. Therefore, we, therefore, assessed the effects of such barriers on the population structure and effective dispersal of five macroinvertebrate species with strictly aquatic life cycles: the amphipod crustacean Gammarus fossarum (clade 11), three snail species of the Ancylus fluviatilis species complex and the flatworm Dugesia gonocephala. We studied populations at nine weirs and eight culverts (3 pipes, 5 tunnels), built 33-109 years ago, mainly in the heavily fragmented catchment of the river Ruhr (Sauerland, Germany). To assess fragmentation and barrier effects, we generated genome-wide SNP data using ddRAD sequencing and evaluated clustering, differentiation between populations up- and downstream of each barrier and effective migration rates among sites and across barriers. Additionally, we applied population genomic simulations to assess expected differentiation patterns under different gene flow scenarios. Our data show that populations of all species are highly isolated at regional and local scales within few kilometers. While the regional population structure likely results from historical processes, the strong local differentiation suggests that contemporary dispersal barriers exist. However, we identified significant barrier effects only for pipes (for A. fluviatilis II and III) and few larger weirs (>1.3 m; for D. gonocephala). Therefore, our data suggest that most small in-stream barriers can probably be overcome by all studied taxa frequently enough to prevent fragmentation. However, it remains to be tested if the strong local differentiation is a result of a cumulative effect of small barriers, or if larger in-stream barriers, land use, chemical pollution, urbanization, or a combination of these factors impede gene flow.

Introgressive hybridisation between domestic pigs (Sus scrofa domesticus) and endemic Corsican wild boars (S. s. meridionalis): effects of human-mediated interventions.

Owing to the intensified domestication process with artificial trait selection, introgressive hybridisation between domestic and wild species poses a management problem. Traditional free-range livestock husbandry, as practiced in Corsica and Sardinia, is known to facilitate hybridisation between wild boars and domestic pigs (Sus scrofa). Here, we assessed the genetic distinctness and genome-wide domestic pig ancestry levels of the Corsican wild boar subspecies S. s. meridionalis, with reference to its Sardinian conspecifics, employing a genome-wide single nucleotide polymorphism (SNP) assay and mitochondrial control region (mtCR) haplotypes. We also assessed the reliance of morphological criteria and the melanocortin-1 receptor (MC1R) coat colour gene to identify individuals with domestic introgression. While Corsican wild boars showed closest affinity to Sardinian and Italian wild boars compared to other European populations based on principal component analysis, the observation of previously undescribed mtCR haplotypes and high levels of nuclear divergence (Weir's θ > 0.14) highlighted the genetic distinctness of Corsican S. s. meridionalis. Across three complementary analyses of mixed ancestry (i.e., STRUCTURE, PCADMIX, and ELAI), proportions of domestic pig ancestry were estimated at 9.5% in Corsican wild boars, which was significantly higher than in wild boars in Sardinia, where free-range pig keeping was banned in 2012. Comparison of morphologically pure- and hybrid-looking Corsican wild boars suggested a weak correlation between morphological criteria and genome-wide domestic pig ancestry. The study highlights the usefulness of molecular markers to assess the direct impacts of management practices on gene flow between domestic and wild species.

When barcoding fails: development of diagnostic nuclear markers for the sibling caddisfly species Sericostoma personatum (Spence in Kirby & Spence, 1826) and Sericostoma flavicorne Schneider, 1845.

The larval stages of the central European sibling caddisfly species Sericostoma personatum (Spence in Kirby and Spence, 1826) and S. flavicorne Schneider, 1845 are morphologically similar and can only be distinguished by differences in coloration in late larval instars. Identification using the mitochondrial barcoding gene, i.e., the Cytochrome c Oxidase 1, is impossible, as both species share the same highly differentiated haplotypes due to introgression. Nuclear gene markers obtained through double digest restriction site associate sequencing (ddRAD seq), however, can reliably distinguish both species, yet the method is expensive as well as time-consuming and therefore not practicable for species determination. To facilitate accurate species identification without sequencing genome-wide markers, we developed nine diagnostic nuclear RFLP markers based on ddRAD seq data. The markers were successfully tested on geographically distinct populations of the two Sericostoma species in western Germany, on known hybrids, and on another sericostomatid caddisfly species, Oecismus monedula (Hagen, 1859) that sometimes shares the habitat and can be morphologically confounded with Sericostoma. We describe a simple and fast protocol for reliable species identification of S. personatum and S. flavicorne independent of the life cycle stage of the specimens.

Company matters: The presence of other genotypes alters traits and intraspecific selection in an Arctic diatom under climate change.

Arctic phytoplankton and their response to future conditions shape one of the most rapidly changing ecosystems on the planet. We tested how much the phenotypic responses of strains from the same Arctic diatom population diverge and whether the physiology and intraspecific composition of multistrain populations differs from expectations based on single strain traits. To this end, we conducted incubation experiments with the diatom Thalassiosira hyalina under present-day and future temperature and pCO2 treatments. Six fresh isolates from the same Svalbard population were incubated as mono- and multistrain cultures. For the first time, we were able to closely follow intraspecific selection within an artificial population using microsatellites and allele-specific quantitative PCR. Our results showed not only that there is substantial variation in how strains of the same species cope with the tested environments but also that changes in genotype composition, production rates, and cellular quotas in the multistrain cultures are not predictable from monoculture performance. Nevertheless, the physiological responses as well as strain composition of the artificial populations were highly reproducible within each environment. Interestingly, we only detected significant strain sorting in those populations exposed to the future treatment. This study illustrates that the genetic composition of populations can change on very short timescales through selection from the intraspecific standing stock, indicating the potential for rapid population level adaptation to climate change. We further show that individuals adjust their phenotype not only in response to their physicochemical but also to their biological surroundings. Such intraspecific interactions need to be understood in order to realistically predict ecosystem responses to global change.

DNA barcode reference libraries for the monitoring of aquatic biota in Europe: Gap-analysis and recommendations for future work.

Effective identification of species using short DNA fragments (DNA barcoding and DNA metabarcoding) requires reliable sequence reference libraries of known taxa. Both taxonomically comprehensive coverage and content quality are important for sufficient accuracy. For aquatic ecosystems in Europe, reliable barcode reference libraries are particularly important if molecular identification tools are to be implemented in biomonitoring and reports in the context of the EU Water Framework Directive (WFD) and the Marine Strategy Framework Directive (MSFD). We analysed gaps in the two most important reference databases, Barcode of Life Data Systems (BOLD) and NCBI GenBank, with a focus on the taxa most frequently used in WFD and MSFD. Our analyses show that coverage varies strongly among taxonomic groups, and among geographic regions. In general, groups that were actively targeted in barcode projects (e.g. fish, true bugs, caddisflies and vascular plants) are well represented in the barcode libraries, while others have fewer records (e.g. marine molluscs, ascidians, and freshwater diatoms). We also found that species monitored in several countries often are represented by barcodes in reference libraries, while species monitored in a single country frequently lack sequence records. A large proportion of species (up to 50%) in several taxonomic groups are only represented by private data in BOLD. Our results have implications for the future strategy to fill existing gaps in barcode libraries, especially if DNA metabarcoding is to be used in the monitoring of European aquatic biota under the WFD and MSFD. For example, missing species relevant to monitoring in multiple countries should be prioritized for future collaborative programs. We also discuss why a strategy for quality control and quality assurance of barcode reference libraries is needed and recommend future steps to ensure full utilisation of metabarcoding in aquatic biomonitoring.

Fishing in troubled waters: Revealing genomic signatures of local adaptation in response to freshwater pollutants in two macroinvertebrates.

Local adaptation is of fundamental importance for populations to cope with fast, human-mediated environmental changes. In the past, analyses of local adaptation were restricted to few model species. Nowadays, due to the increased affordability of high-throughput sequencing, local adaptation can be studied much easier by searching for patterns of positive selection using genomic data. In the present study, we analysed effects of wastewater treatment plant and ore mining effluents on stream invertebrate populations. The two different anthropogenic stressors have impacted on stream ecosystems over different time scales and with different potencies. As target organisms we selected two macroinvertebrate species with different life histories and dispersal capacities: the caddisfly Glossosoma conformis and the flatworm Dugesia gonocephala. We applied a genome-wide genetic marker technique, termed ddRAD (double digest restriction site associated DNA) sequencing, to identify local adaptation. Ten and 18% of all loci were identified as candidate loci for local adaptation in D. gonocephala and G. conformis, respectively. However, after stringent re-evaluation of the genomic data, strong evidence for local adaptation remained only for one population of the flatworm D. gonocephala affected by high copper concentration from ore mining. One of the corresponding candidate loci is arnt, a gene associated with the response to xenobiotics and potentially involved in metal detoxification. Our results support the hypotheses that local adaptation is more likely to play a central role in environments impacted by a stronger stressor for a longer time and that it is more likely to occur in species with lower migration rates. However, these findings have to be interpreted cautiously, as several confounding factors may have limited the possibility to detect local adaptation. Our study highlights how genomic tools can be used to study the adaptability and thus resistance of natural populations to changing environments and we discuss prospects and limitations of the methods.

Genome-wide single-nucleotide polymorphism data reveal cryptic species within cryptic freshwater snail species-The case of the Ancylus fluviatilis species complex.

DNA barcoding utilizes short standardized DNA sequences to identify species and is increasingly used in biodiversity assessments. The technique has unveiled an unforeseeably high number of morphologically cryptic species. However, if speciation has occurred relatively recently and rapidly, the use of single gene markers, and especially the exclusive use of mitochondrial markers, will presumably fail in delimitating species. Therefore, the true number of biological species might be even higher. One mechanism that can result in rapid speciation is hybridization of different species in combination with polyploidization, that is, allopolyploid speciation. In this study, we analyzed the population genetic structure of the polyploid freshwater snail Ancylus fluviatilis, for which allopolyploidization was postulated as a speciation mechanism. DNA barcoding has already revealed four cryptic species within A. fluviatilis (i.e., A. fluviatilis s. str., Ancylus sp. A-C), but early allozyme data even hint at the presence of additional cryptic lineages in Central Europe. We combined COI sequencing with high-resolution genome-wide SNP data (ddRAD data) to analyze the genetic structure of A. fluviatilis populations in a Central German low mountain range (Sauerland). The ddRAD data results indicate the presence of three cryptic species within A. fluviatilis s. str. occurring in sympatry and even syntopy, whereas mitochondrial sequence data only support the existence of one species, with shared haplotypes between species. Our study hence points to the limitations of DNA barcoding when dealing with organismal groups where speciation is assumed to have occurred rapidly, for example, through the process of allopolyploidization. We therefore emphasize that single marker DNA barcoding can underestimate the true species diversity and argue in strong favor of using genome-wide data for species delimitation in such groups.

ddrage: A data set generator to evaluate ddRADseq analysis software.

High-throughput sequencing makes it possible to evaluate thousands of genetic markers across genomes and populations. Reduced-representation sequencing approaches, like double-digest restriction site-associated DNA sequencing (ddRADseq), are frequently applied to screen for genetic variation. In particular in nonmodel organisms where whole-genome sequencing is not yet feasible, ddRADseq has become popular as it allows genomewide assessment of variation patterns even in the absence of other genomic resources. However, while many tools are available for the analysis of ddRADseq data, few options exist to simulate ddRADseq data in order to evaluate the accuracy of downstream tools. The available tools either focus on the optimization of ddRAD experiment design or do not provide the information necessary for a detailed evaluation of different ddRAD analysis tools. For this task, a ground truth, that is, the underlying information of all effects in the data set, is required. Therefore, we here present ddrage, the ddRAD Data Set Generator, that allows both developers and users to evaluate their ddRAD analysis software. ddrage allows the user to adjust many parameters such as coverage and rates of mutations, sequencing errors or allelic dropouts, in order to generate a realistic simulated ddRADseq data set for given experimental scenarios and organisms. The simulated reads can be easily processed with available analysis software such as stacks or pyrad and evaluated against the underlying parameters used to generate the data to gauge the impact of different parameter values used during downstream data processing.

Deciphering the origin of mito-nuclear discordance in two sibling caddisfly species.

An increasing number of phylogenetic studies have reported discordances among nuclear and mitochondrial markers. These discrepancies are highly relevant to widely used biodiversity assessment approaches, such as DNA barcoding, that rely almost exclusively on mitochondrial markers. Although the theoretical causes of mito-nuclear discordances are well understood, it is often extremely challenging to determine the principal underlying factor in a given study system. In this study, we uncovered significant mito-nuclear discordances in a pair of sibling caddisfly species. Application of genome sequencing, ddRAD and DNA barcoding revealed ongoing hybridization, as well as historical hybridization in Pleistocene refugia, leading us to identify introgression as the ultimate cause of the observed discordance pattern. Our novel genomic data, the discovery of a European-wide hybrid zone and the availability of established techniques for laboratory breeding make this species pair an ideal model system for studying species boundaries with ongoing gene flow.

RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity.

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher Fst values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.