Single metabarcoding multiplex captures community‐level freshwater biodiversity and beyond

Abstract Cost‐effective and accurate quantification of biodiversity is important for biodiversity conservation, resource management, and forecasting. Traditional monitoring approaches have relied on direct observations, remote sensing, and mark‐recapture techniques, providing insights into species ecology and the impact of pollution and climate change on indicator species. However, these techniques are typically low throughput, expensive, and can be invasive. In addition, they cannot detect cryptic diversity and are biased toward species that leave identifiable remains. DNA‐based methods, such as metabarcoding or single marker gene assays, have enabled high throughput screening of a wide range of taxonomic groups, including ones without well‐preserved remains. When compared with traditional techniques, these approaches have high throughput, can resolve cryptic diversity, do not require taxonomic specialist skills, and are non‐invasive. However, although they are comparatively cheaper than traditional approaches, they are expensive when applied at the community‐level as single marker assays are amplified and sequenced independently. Multilocus approaches in which multiple gene markers are amplified in a single reaction are desirable to deliver community‐level assessments in a cost‐effective manner. Yet, they are uncommon because of technical challenges that may lead to biases in downstream analyses, such as index hopping and unbalanced representation of taxonomic groups. Here, we developed a highly multiplexed protocol that combines the early pooling of marker genes that target broad taxonomic groups and taxon‐specific markers in a single tube reaction. This step is followed by the pooling of up to 384 samples per locus (N = 15,636 samples) with unique dual‐indexed sequencing adapters in a single sequencing run. This approach dramatically reduces the costs of community‐level biodiversity quantification and lowers the need for input DNA without compromising output quality. We optimized the multiplex assay on lake freshwater sediment samples and benchmarked the assay on samples from other environmental matrices, demonstrating its direct application to the river and marine communities.