Jump dispersal drives the relationship between micro- and macroevolutionary dynamics in the Sicydiinae (Gobiiformes: Oxudercidae) of Sundaland and Wallacea.

Insular biodiversity hotspots of Southeast Asia are remarkable for their biodiverse faunas. With a marine larval phase lasting up to several months, the freshwater fish subfamily Sicydiinae has colonized most islands of these hotspots. However, Sicydiinae diversity is still poorly understood in Southeast Asia. With the objective to estimate intraspecific genetic diversity and infer past demography, we conducted the molecular inventory of Sicydiinae species in Sundaland and Wallacea using 652 bp of the mitochondrial cytochrome oxidase I gene, species delimitation methods and Bayesian Skyline plot reconstructions. In total, 24 Molecular Operational Taxonomic Units are delimited among the 603 sequences belonging to 27 species and five genera. Two cases of discordance between morphology and mitochondrial sequence are observed suggesting ongoing speciation and/or introgression in two genera. Multiple new occurrences are reported, either for a single biodiversity hotspot or both, some of which corresponding to observations of a few individuals far from the range distribution of their conspecifics. Among the ten species or species group whose intraspecific diversity was examined, high levels of genetic diversity and past population expansion are revealed by Tajima's D tests and Bayesian Skyline Plot reconstructions. Together these results indicate that long-distance dispersal is common and suggest that most endemic species originated through founder events followed by population expansion. Patterns of sexual dimorphism and males' coloration among diverging species pair seem to point to sexual selection as an important mechanism contributing to speciation in the Sicydiinae of Sundaland and Wallacea.

Predicting environmental stressor levels with machine learning: a comparison between amplicon sequencing, metagenomics, and total RNA sequencing based on taxonomically assigned data.

Introduction: Microbes are increasingly (re)considered for environmental assessments because they are powerful indicators for the health of ecosystems. The complexity of microbial communities necessitates powerful novel tools to derive conclusions for environmental decision-makers, and machine learning is a promising option in that context. While amplicon sequencing is typically applied to assess microbial communities, metagenomics and total RNA sequencing (herein summarized as omics-based methods) can provide a more holistic picture of microbial biodiversity at sufficient sequencing depths. Despite this advantage, amplicon sequencing and omics-based methods have not yet been compared for taxonomy-based environmental assessments with machine learning. Methods: In this study, we applied 16S and ITS-2 sequencing, metagenomics, and total RNA sequencing to samples from a stream mesocosm experiment that investigated the impacts of two aquatic stressors, insecticide and increased fine sediment deposition, on stream biodiversity. We processed the data using similarity clustering and denoising (only applicable to amplicon sequencing) as well as multiple taxonomic levels, data types, feature selection, and machine learning algorithms and evaluated the stressor prediction performance of each generated model for a total of 1,536 evaluated combinations of taxonomic datasets and data-processing methods. Results: Sequencing and data-processing methods had a substantial impact on stressor prediction. While omics-based methods detected a higher diversity of taxa than amplicon sequencing, 16S sequencing outperformed all other sequencing methods in terms of stressor prediction based on the Matthews Correlation Coefficient. However, even the highest observed performance for 16S sequencing was still only moderate. Omics-based methods performed poorly overall, but this was likely due to insufficient sequencing depth. Data types had no impact on performance while feature selection significantly improved performance for omics-based methods but not for amplicon sequencing. Discussion: We conclude that amplicon sequencing might be a better candidate for machine-learning-based environmental stressor prediction than omics-based methods, but the latter require further research at higher sequencing depths to confirm this conclusion. More sampling could improve stressor prediction performance, and while this was not possible in the context of our study, thousands of sampling sites are monitored for routine environmental assessments, providing an ideal framework to further refine the approach for possible implementation in environmental diagnostics.

Metagenomics versus total RNA sequencing: most accurate data-processing tools, microbial identification accuracy and perspectives for ecological assessments.

Metagenomics and total RNA sequencing (total RNA-Seq) have the potential to improve the taxonomic identification of diverse microbial communities, which could allow for the incorporation of microbes into routine ecological assessments. However, these target-PCR-free techniques require more testing and optimization. In this study, we processed metagenomics and total RNA-Seq data from a commercially available microbial mock community using 672 data-processing workflows, identified the most accurate data-processing tools, and compared their microbial identification accuracy at equal and increasing sequencing depths. The accuracy of data-processing tools substantially varied among replicates. Total RNA-Seq was more accurate than metagenomics at equal sequencing depths and even at sequencing depths almost one order of magnitude lower than those of metagenomics. We show that while data-processing tools require further exploration, total RNA-Seq might be a favorable alternative to metagenomics for target-PCR-free taxonomic identifications of microbial communities and might enable a substantial reduction in sequencing costs while maintaining accuracy. This could be particularly an advantage for routine ecological assessments, which require cost-effective yet accurate methods, and might allow for the incorporation of microbes into ecological assessments.

Mitochondrial genome sequencing, mapping, and assembly benchmarking for Culicoides species (Diptera: Ceratopogonidae).

BACKGROUND: Mitochondrial genomes are the most sequenced genomes after bacterial and fungal genomic DNA. However, little information on mitogenomes is available for multiple metazoan taxa, such as Culicoides, a globally distributed, megadiverse genus containing 1,347 species. AIM:  Generating novel mitogenomic information from single Culicoides sonorensis and C. biguttatus specimens, comparing available mitogenome mapping and de novo assembly tools, and identifying the best performing strategy and tools for Culicoides species. RESULTS: We present two novel and fully annotated mitochondrial haplotypes for two Culicoides species, C. sonorensis and C. biguttatus. We also annotated or re-annotated the only available reference mitogenome for C. sonorensis and C. arakawae. All species present a high similarity in mitogenome organization. The general gene arrangement for all Culicoides species was identical to the ancestral insect mitochondrial genome. Only short spacers were found in C. sonorensis (up to 30 bp), contrary to C. biguttatus (up to 114 bp). The mitochondrial genes ATP8, NAD2, NAD6, and LSU rRNA exhibited the highest nucleotide diversity and pairwise interspecific p genetic distance, suggesting that these genes might be suitable and complementary molecular barcodes for Culicoides identification in addition to the commonly utilized COI gene. We observed performance differences between the compared mitogenome generation strategies. The mapping strategy outperformed the de novo assembly strategy, but mapping results were partially biased in the absence of species-specific reference mitogenome. Among the utilized tools, BWA performed best for C. sonorensis while SPAdes, MEGAHIT, and MitoZ were among the best for C. biguttatus. The best-performing mitogenome annotator was MITOS2. Additionally, we were able to recover exogenous mitochondrial DNA from Bos taurus (biting midges host) from a C. biguttatus blood meal sample. CONCLUSIONS: Two novel annotated mitogenome haplotypes for C. sonorensis and C. biguttatus using High-Throughput Sequencing are presented. Current results are useful as the baseline for mitogenome reconstruction of the remaining Culicoides species from single specimens to HTS and genome annotation. Mapping to a species-specific reference mitogenome generated better results for Culicoides mitochondrial genome reconstruction than de novo assembly, while de novo assembly resulted better in the absence of a closely related reference mitogenome. These results have direct implications for molecular-based identification of these vectors of human and zoonotic diseases, setting the basis for using the whole mitochondrial genome as a marker in Culicoides identification.

A bright idea-metabarcoding arthropods from light fixtures.

Arthropod communities in buildings have not been extensively studied, although humans have always shared their homes with them. In this study we explored if arthropod DNA can be retrieved and metabarcoded from indoor environments through the collection of dead specimens in light fixtures to better understand what shapes arthropod diversity in our homes. Insects were collected from 45 light fixtures at the Centre for Biodiversity Genomics (CBG, Guelph, Canada), and by community scientists at 12 different residential homes in Southern Ontario. The CBG ground floor of the CBG showed the greatest arthropod diversity, especially in light fixtures that were continuously illuminated. The community scientist samples varied strongly by light fixture type, lightbulb used, time passed since lamp was last cleaned, and specimen size. In all cases, the majority of OTUs was not shared between samples even within the same building. This study demonstrates that light fixtures might be a useful resource to determine arthropod diversity in our homes, but individual samples are likely not representative of the full diversity.

Molecular Taxonomy and Diversification of Atlantic Skates (Chondrichthyes, Rajiformes): Adding More Pieces to the Puzzle of Their Evolutionary History.

Conservation and long-term management plans of marine species need to be based upon the universally recognized key-feature of species identity. This important assignment is particularly challenging in skates (Rajiformes) in which the phenotypic similarity between some taxa and the individual variability in others, hampers accurate species identification. Here, 432 individual skate samples collected from four major ocean areas of the Atlantic were barcoded and taxonomically analysed. A BOLD project ELASMO ATL was implemented with the aim of establishing a new fully available and well curated barcode library containing both biological and molecular information. The evolutionary histories of the 38 skate taxa were estimated with two concatenated mitochondrial markers (COI and NADH2) through Maximum Likelihood and Bayesian inference. New evolutionary lineages within the genus Raja were discovered off Angola, where paleogeographic history coupled with oceanographic discontinuities could have contributed to the establishment of isolated refugia, playing a fundamental role among skates' speciation events. These data successfully resolved many taxonomic ambiguities, identified cryptic diversity within valid species and demonstrated a highly cohesive monophyletic clustering among the order, laying the background for further inference of evolutionary patterns suitable for addressing management and conservation issues.

Exploring the vertebrate fauna of the Bird's Head Peninsula (Indonesia, West Papua) through DNA barcodes.

Biodiversity knowledge is widely heterogeneous across the Earth's biomes. Some areas, due to their remoteness and difficult access, present large taxonomic knowledge gaps. Mostly located in the tropics, these areas have frequently experienced a fast development of anthropogenic activities during the last decades and are therefore of high conservation concerns. The biodiversity hotspots of Southeast Asia exemplify the stakes faced by tropical countries. While the hotspots of Sundaland (Java, Sumatra, Borneo) and Wallacea (Sulawesi, Moluccas) have long attracted the attention of biologists and conservationists alike, extensive parts of the Sahul area, in particular the island of New Guinea, have been much less explored biologically. Here, we describe the results of a DNA-based inventory of aquatic and terrestrial vertebrate communities, which was the objective of a multidisciplinary expedition to the Bird's Head Peninsula (West Papua, Indonesia) conducted between 17 October and 20 November 2014. This expedition resulted in the assembly of 1005 vertebrate DNA barcodes. Based on the use of multiple species-delimitation methods (GMYC, PTP, RESL, ABGD), 264 molecular operational taxonomic units (MOTUs) were delineated, among which 75 were unidentified and an additional 48 were considered cryptic. This study suggests that the diversity of vertebrates of the Bird's Head is severely underestimated and considerations on the evolutionary origin and taxonomic knowledge of these biotas are discussed.

Assessment of current taxonomic assignment strategies for metabarcoding eukaryotes.

The effective use of metabarcoding in biodiversity science has brought important analytical challenges due to the need to generate accurate taxonomic assignments. The assignment of sequences to genus or species level is critical for biodiversity surveys and biomonitoring, but it is particularly challenging as researchers must select the approach that best recovers information on species composition. This study evaluates the performance and accuracy of seven methods in recovering the species composition of mock communities by using COI barcode fragments. The mock communities varied in species number and specimen abundance, while upstream molecular and bioinformatic variables were held constant, and using a set of COI fragments. We evaluated the impact of parameter optimization on the quality of the predictions. Our results indicate that BLAST top hit competes well with more complex approaches if optimized for the mock community under study. For example, the two machine learning methods that were benchmarked proved more sensitive to reference database heterogeneity and completeness than methods based on sequence similarity. The accuracy of assignments was impacted by both species and specimen counts (query compositional heterogeneity) which ultimately influence the selection of appropriate software. We urge researchers to: (i) use realistic mock communities to allow optimization of parameters, regardless of the taxonomic assignment method employed; (ii) carefully choose and curate the reference databases including completeness; and (iii) use QIIME, BLAST or LCA methods, in conjunction with parameter tuning to better assign taxonomy to diverse communities, especially when information on species diversity is lacking for the area under study.

Biodiversity inventory of the grey mullets (Actinopterygii: Mugilidae) of the Indo-Australian Archipelago through the iterative use of DNA-based species delimitation and specimen assignment methods.

DNA barcoding opens new perspectives on the way we document biodiversity. Initially proposed to circumvent the limits of morphological characters to assign unknown individuals to known species, DNA barcoding has been used in a wide array of studies where collecting species identity constitutes a crucial step. The assignment of unknowns to knowns assumes that species are already well identified and delineated, making the assignment performed reliable. Here, we used DNA-based species delimitation and specimen assignment methods iteratively to tackle the inventory of the Indo-Australian Archipelago grey mullets, a notorious case of taxonomic complexity that requires DNA-based identification methods considering that traditional morphological identifications are usually not repeatable and sequence mislabeling is common in international sequence repositories. We first revisited a DNA barcode reference library available at the global scale for Mugilidae through different DNA-based species delimitation methods to produce a robust consensus scheme of species delineation. We then used this curated library to assign unknown specimens collected throughout the Indo-Australian Archipelago to known species. A second iteration of OTU delimitation and specimen assignment was then performed. We show the benefits of using species delimitation and specimen assignment methods iteratively to improve the accuracy of specimen identification and propose a workflow to do so.

Assessing species diversity of Coral Triangle artisanal fisheries: A DNA barcode reference library for the shore fishes retailed at Ambon harbor (Indonesia).

The Coral Triangle (CT), a region spanning across Indonesia and Philippines, is home to about 4,350 marine fish species and is among the world's most emblematic regions in terms of conservation. Threatened by overfishing and oceans warming, the CT fisheries have faced drastic declines over the last decades. Usually monitored through a biomass-based approach, fisheries trends have rarely been characterized at the species level due to the high number of taxa involved and the difficulty to accurately and routinely identify individuals to the species level. Biomass, however, is a poor proxy of species richness, and automated methods of species identification are required to move beyond biomass-based approaches. Recent meta-analyses have demonstrated that species richness peaks at intermediary levels of biomass. Consequently, preserving biomass is not equal to preserving biodiversity. We present the results of a survey to estimate the shore fish diversity retailed at the harbor of Ambon Island, an island located at the center of the CT that display exceptionally high biomass despite high levels of threat, while building a DNA barcode reference library of CT shore fishes targeted by artisanal fisheries. We sampled 1,187 specimens and successfully barcoded 696 of the 760 selected specimens that represent 202 species. Our results show that DNA barcodes were effective in capturing species boundaries for 96% of the species examined, which opens new perspectives for the routine monitoring of the CT fisheries.

Disentangling the taxonomy of the subfamily Rasborinae (Cypriniformes, Danionidae) in Sundaland using DNA barcodes.

Sundaland constitutes one of the largest and most threatened biodiversity hotspots; however, our understanding of its biodiversity is afflicted by knowledge gaps in taxonomy and distribution patterns. The subfamily Rasborinae is the most diversified group of freshwater fishes in Sundaland. Uncertainties in their taxonomy and systematics have constrained its use as a model in evolutionary studies. Here, we established a DNA barcode reference library of the Rasborinae in Sundaland to examine species boundaries and range distributions through DNA-based species delimitation methods. A checklist of the Rasborinae of Sundaland was compiled based on online catalogs and used to estimate the taxonomic coverage of the present study. We generated a total of 991 DNA barcodes from 189 sampling sites in Sundaland. Together with 106 previously published sequences, we subsequently assembled a reference library of 1097 sequences that covers 65 taxa, including 61 of the 79 known Rasborinae species of Sundaland. Our library indicates that Rasborinae species are defined by distinct molecular lineages that are captured by species delimitation methods. A large overlap between intraspecific and interspecific genetic distance is observed that can be explained by the large amounts of cryptic diversity as evidenced by the 166 Operational Taxonomic Units detected. Implications for the evolutionary dynamics of species diversification are discussed.

A reference library for Canadian invertebrates with 1.5 million barcodes, voucher specimens, and DNA samples.

The reliable taxonomic identification of organisms through DNA sequence data requires a well parameterized library of curated reference sequences. However, it is estimated that just 15% of described animal species are represented in public sequence repositories. To begin to address this deficiency, we provide DNA barcodes for 1,500,003 animal specimens collected from 23 terrestrial and aquatic ecozones at sites across Canada, a nation that comprises 7% of the planet's land surface. In total, 14 phyla, 43 classes, 163 orders, 1123 families, 6186 genera, and 64,264 Barcode Index Numbers (BINs; a proxy for species) are represented. Species-level taxonomy was available for 38% of the specimens, but higher proportions were assigned to a genus (69.5%) and a family (99.9%). Voucher specimens and DNA extracts are archived at the Centre for Biodiversity Genomics where they are available for further research. The corresponding sequence and taxonomic data can be accessed through the Barcode of Life Data System, GenBank, the Global Biodiversity Information Facility, and the Global Genome Biodiversity Network Data Portal.

Validation of COI metabarcoding primers for terrestrial arthropods.

Metabarcoding can rapidly determine the species composition of bulk samples and thus aids biodiversity and ecosystem assessment. However, it is essential to use primer sets that minimize amplification bias among taxa to maximize species recovery. Despite this fact, the performance of primer sets employed for metabarcoding terrestrial arthropods has not been sufficiently evaluated. This study tests the performance of 36 primer sets on a mock community containing 374 insect species. Amplification success was assessed with gradient PCRs and the 21 most promising primer sets selected for metabarcoding. These 21 primer sets were also tested by metabarcoding a Malaise trap sample. We identified eight primer sets, mainly those including inosine and/or high degeneracy, that recovered more than 95% of the species in the mock community. Results from the Malaise trap sample were congruent with the mock community, but primer sets generating short amplicons produced potential false positives. Taxon recovery from both mock community and Malaise trap sample metabarcoding were used to select four primer sets for additional evaluation at different annealing temperatures (40-60 °C) using the mock community. The effect of temperature varied by primer pair but overall it only had a minor effect on taxon recovery. This study reveals the weak performance of some primer sets employed in past studies. It also demonstrates that certain primer sets can recover most taxa in a diverse species assemblage. Thus, based our experimental set up, there is no need to employ several primer sets targeting the same gene region. We identify several suitable primer sets for arthropod metabarcoding, and specifically recommend BF3 + BR2, as it is not affected by primer slippage and provides maximal taxonomic resolution. The fwhF2 + fwhR2n primer set amplifies a shorter fragment and is therefore ideal when targeting degraded DNA (e.g., from gut contents).

Survey of mislabelling across finfish supply chain reveals mislabelling both outside and within Canada.

Seafood has become one of the most heavily traded food commodities in the era of globalization. International seafood supply chains are complex and contend with many difficulties in bringing an enormous variety of products to market. A major challenge involves accurately labelling products such that they comply with a diverse set of regulatory frameworks, ranging from country-of-origin through to the final point of consumer sale. Thanks to DNA barcoding, seafood mislabelling is now recognized as a global problem, with potentially negative impacts on human health, economy and the environment. Mislabelling can result from species misidentification, use of inappropriate common names, incomplete and/or out-dated regulatory frameworks, or through market substitution. While prior studies have focused primarily on retail and food service establishments, this study used barcoding to assess rates of finfish mislabelling at multiple points in the supply chain within Ontario, Canada. A total of 203 specimens from 12 key targeted species were collected from varied importers, registered processing plants and retailers in Southern Ontario and identified using DNA barcoding. Species identity of samples was used to assess conformity of labelling against the Canadian Food Inspection Agency's (CFIA) Fish List, which revealed an overall mislabelling rate of 32.3% among targeted species. The mislabelling rate was significantly different between samples collected from importers and retailers. Among the mislabelled samples were seven samples that originated from US and were properly labelled according to US Food and Drug Administration (FDA) Seafood List. This study evaluated the integrity of chain of custody documents and identified discrepancies in 43 samples (21.4%). Implementing seafood traceability throughout the supply chain and harmonizing labelling regulations between countries can help to ensure industry compliance in a globalized market, while sampling at multiple points in the supply chain can help to reveal causes.

Metabarcoding a diverse arthropod mock community.

Although DNA metabarcoding is an attractive approach for monitoring biodiversity, it is often difficult to detect all the species present in a bulk sample. In particular, sequence recovery for a given species depends on its biomass and mitome copy number as well as the primer set employed for PCR. To examine these variables, we constructed a mock community of terrestrial arthropods comprised of 374 species. We used this community to examine how species recovery was impacted when amplicon pools were constructed in four ways. The first two protocols involved the construction of bulk DNA extracts from different body segments (Bulk Abdomen, Bulk Leg). The other protocols involved the production of DNA extracts from single legs which were then merged prior to PCR (Composite Leg) or PCR-amplified separately (Single Leg) and then pooled. The amplicons generated by these four treatments were then sequenced on three platforms (Illumina MiSeq, Ion Torrent PGM and Ion Torrent S5). The choice of sequencing platform did not substantially influence species recovery, although the Miseq delivered the highest sequence quality. As expected, species recovery was most efficient from the Single Leg treatment because amplicon abundance varied little among taxa. Among the three treatments where PCR occurred after pooling, the Bulk Abdomen treatment produced a more uniform read abundance than the Bulk Leg or Composite Leg treatment. Primer choice also influenced species recovery and evenness. Our results reveal how variation in protocols can have substantial impacts on perceived diversity unless sequencing coverage is sufficient to reach an asymptote.

When too much isn't enough: Does current food production meet global nutritional needs?

Sustainably feeding the next generation is often described as one of the most pressing "grand challenges" facing the 21st century. Generally, scholars propose addressing this problem by increasing agricultural production, investing in technology to boost yields, changing diets, or reducing food waste. In this paper, we explore whether global food production is nutritionally balanced by comparing the diet that nutritionists recommend versus global agricultural production statistics. Results show that the global agricultural system currently overproduces grains, fats, and sugars while production of fruits and vegetables and protein is not sufficient to meet the nutritional needs of the current population. Correcting this imbalance could reduce the amount of arable land used by agriculture by 51 million ha globally but would increase total land used for agriculture by 407 million ha and increase greenhouse gas emissions. For a growing population, our calculations suggest that the only way to eat a nutritionally balanced diet, save land and reduce greenhouse gas emissions is to consume and produce more fruits and vegetables as well as transition to diets higher in plant-based protein. Such a move will help protect habitats and help meet the Sustainable Development Goals.

Slippage of degenerate primers can cause variation in amplicon length.

It is well understood that homopolymer regions should be avoided for primer binding to prevent off-target amplification. However, in metabarcoding, it is often difficult to avoid primer degeneracy in order to maximize taxa detection. We here investigate primer binding specificity using different primer sets from several invertebrate metabarcoding studies. Our results indicate that primers frequently bound 1-2 bp upstream in taxa where a homopolymer region was present in the amplification direction. Primer binding 1 bp downstream was observed less frequently. This primer slippage leads to taxon-specific length variation in amplicons and subsequent length variation in recovered sequences. Some widely used primer sets were severely affected by this bias, while others were not. While this variation will only have small impacts on the designation of Operational Taxonomic Units (OTUs) by clustering algorithms that ignore terminal gaps, primer sets employed in metabarcoding projects should be evaluated for their sensitivity to slippage. Moreover, steps should be taken to reduce slippage by improving protocols for primer design. For example, the flanking region adjacent to the 3' end of the primer is not considered by current primer development software although GC clamps in this position could mitigate slippage.

Estimating intraspecific genetic diversity from community DNA metabarcoding data.

BACKGROUND: DNA metabarcoding is used to generate species composition data for entire communities. However, sequencing errors in high-throughput sequencing instruments are fairly common, usually requiring reads to be clustered into operational taxonomic units (OTUs), losing information on intraspecific diversity in the process. While Cytochrome c oxidase subunit I (COI) haplotype information is limited in resolving intraspecific diversity it is nevertheless often useful e.g. in a phylogeographic context, helping to formulate hypotheses on taxon distribution and dispersal. METHODS: This study combines sequence denoising strategies, normally applied in microbial research, with additional abundance-based filtering to extract haplotype information from freshwater macroinvertebrate metabarcoding datasets. This novel approach was added to the R package "JAMP" and can be applied to COI amplicon datasets. We tested our haplotyping method by sequencing (i) a single-species mock community composed of 31 individuals with 15 different haplotypes spanning three orders of magnitude in biomass and (ii) 18 monitoring samples each amplified with four different primer sets and two PCR replicates. RESULTS: We detected all 15 haplotypes of the single specimens in the mock community with relaxed filtering and denoising settings. However, up to 480 additional unexpected haplotypes remained in both replicates. Rigorous filtering removes most unexpected haplotypes, but also can discard expected haplotypes mainly from the small specimens. In the monitoring samples, the different primer sets detected 177-200 OTUs, each containing an average of 2.40-3.30 haplotypes per OTU. The derived intraspecific diversity data showed population structures that were consistent between replicates and similar between primer pairs but resolution depended on the primer length. A closer look at abundant taxa in the dataset revealed various population genetic patterns, e.g. the stonefly Taeniopteryx nebulosa and the caddisfly Hydropsyche pellucidula showed a distinct north-south cline with respect to haplotype distribution, while the beetle Oulimnius tuberculatus and the isopod Asellus aquaticus displayed no clear population pattern but differed in genetic diversity. DISCUSSION: We developed a strategy to infer intraspecific genetic diversity from bulk invertebrate metabarcoding data. It needs to be stressed that at this point this metabarcoding-informed haplotyping is not capable of capturing the full diversity present in such samples, due to variation in specimen size, primer bias and loss of sequence variants with low abundance. Nevertheless, for a high number of species intraspecific diversity was recovered, identifying potentially isolated populations and taxa for further more detailed phylogeographic investigation. While we are currently lacking large-scale metabarcoding datasets to fully take advantage of our new approach, metabarcoding-informed haplotyping holds great promise for biomonitoring efforts that not only seek information about species diversity but also underlying genetic diversity.