The study took place in the Hiji River on Shikoku Island, Japan, a catchment of 1,210 square kilometres that flows from the hills of Ehime Prefecture down to the Iyo Sea. Between 6 and 12 May 2021, the researchers sampled 60 sites spanning the river’s full range of environmental conditions, ultimately analysing 59 of them. At each site they used a kick-net to disturb the streambed and sweep up macroinvertebrates, collected in triplicate at upstream, middle, and downstream points. Alongside the animals themselves, they measured water temperature, pH, electrical conductivity, and dissolved oxygen on the spot, and carried back surface water for laboratory analysis of nitrogen and phosphorus compounds, the chemical signatures of nutrient pollution and productivity.

What makes this study remarkable is the technology behind it. Traditional population genetics requires capturing individual animals, identifying them under a microscope, and sequencing one specimen at a time, an approach so costly and laborious that most studies are limited to a single species or a handful of them. Haplotype-level bulk DNA metabarcoding changes the calculus entirely. The entire community sample, dozens of species ground up together, is processed as one. The team extracted DNA using the phenol-chloroform-isoamyl alcohol method, amplified a 421-base-pair region of the mitochondrial cytochrome oxidase subunit 1 gene, the workhorse marker of DNA barcoding, and sequenced the pool on an Illumina MiSeq with 300-base-pair paired-end reads.

The bioinformatics pipeline was correspondingly rigorous. Nearly seven million raw reads passed through the DADA2 denoising pipeline, which corrects sequencing errors and merges forward and reverse reads. The researchers then applied a gauntlet of filters that sets this study apart: sequences were trimmed to the expected amplicon length, checked with MACSE for internal frameshifts, stop codons, and deletions, and screened against an arthropod COI hidden Markov model to flag pseudogenes, the nuclear mitochondrial imposters that can contaminate metabarcoding datasets. After all filtering, 1,297 verified haplotypes remained. Each was assigned to a taxon by querying the Barcode of Life Database, an in-house database of morphologically identified crane flies, and NCBI, using a strict 97 percent identity threshold. The result: 889 haplotypes clustered into 198 operational taxonomic units, with a careful multimodality and phylogenetic screening confirming that only six OTUs, three percent, mixed haplotypes from more than one species, and those were set aside.

From this dataset the researchers distilled the classic currencies of population genetics. Fifty-nine OTUs, spanning seven orders from mayflies and stoneflies to caddisflies, beetles, true flies, dragonflies, and dobsonflies, were detected at five or more sites and subjected to analysis of molecular variance, yielding global GϕST values and pairwise ϕST measures of genetic differentiation. Haplotype richness per OTU ranged from one to 31 unique sequences. The team also compiled a trait proxy for dispersal: female wing length, scored in eight categories from the DISPERSE database, from wingless flies shorter than five millimetres to powerful fliers with wings exceeding fifty millimetres.

The headline finding is a striking asymmetry. Species detected at more sites across the catchment carried significantly more total genetic diversity, the γ-level diversity summed across the whole landscape, than species restricted to few sites. That fits classical theory: broader ranges expose populations to more heterogeneous habitats, more asynchronous local extinctions and recolonizations, and more founder effects and bottlenecks, all of which stack up unique haplotypes across the landscape. Consistent with this, wing length correlated positively with the number of occupied sites, meaning good dispersers really did live in more places. Yet the expected second half of the pattern never materialized: there was no correlation between species occurrence and β-genetic diversity, the genetic differentiation among populations. Wide-ranging species were not more genetically structured than rare ones.

The explanation the authors favour is gene flow. Species that disperse well enough to occupy many sites apparently homogenize their populations faster than extinction-colonization dynamics can differentiate them, erasing the β-diversity signal that habitat heterogeneity generates. Supporting this, neither river distance nor straight-line distance correlated with genetic differentiation in any of the OTUs examined, a pattern of absent isolation-by-distance that previous reviews have noted in aquatic insects, where gene flow routinely overwhelms genetic drift. Intriguingly, wing length itself did not predict genetic differentiation either, hinting that actual gene flow is governed less by wingspan than by species-specific behaviour, oviposition preferences, and demographic history. Some strong fliers, after all, preferentially lay eggs at nearby suitable sites rather than dispersing as far as their wings would allow.

The second major surprise concerns species-genetic diversity correlations, or SGDCs, the hypothesis that species diversity and intraspecific genetic diversity rise and fall together because both respond to the same environmental drivers. Across all 48 taxa tested, not a single one showed a significant α-SGDC, the correlation between local species richness and local haplotype richness. Only one species, a caddisfly of the genus Rhyacophila, displayed a significant β-SGDC, linking its population genetic differentiation to community compositional divergence across sites. The team attributes this near-universal absence to a fundamental decoupling: the mechanisms shaping the two levels of biodiversity differ. In the metapopulation, gene flow homogenizes genes; in the metacommunity, species sorting along environmental gradients, especially elevation, which explained variation in community composition for the majority of well-sampled OTUs, drives species turnover. One process pulls toward uniformity, the other toward divergence, and the correlation collapses.

Environmental analysis added texture to the story. As species occupied more sites, the range of conditions they experienced converged on the catchment-wide mean, and mean pairwise environmental differences correlated positively with haplotype richness. Individual chemicals left species-specific fingerprints: conductivity explained genetic structure in a gomphid dragonfly and a stonefly, phosphate in a limoniid crane fly, and pH in a caddisfly. None of this, however, added up to broad community-wide coupling of genes and species, and the single marginal correlation between pH and GϕST dissolved once OTUs with extreme negative values were excluded.

The implications reach well beyond one Japanese watershed. Haplotype-level metabarcoding opens the door to comparing population genetic structure across dozens of species simultaneously, transforming metapopulation and metacommunity research from single-species case studies into genuine community-scale science. For conservation, that matters enormously: as climate change and habitat degradation accelerate biodiversity loss, understanding whether wide-ranging species truly harbour reservoirs of adaptive genetic variation, or whether gene flow masks fragile structure, will shape which populations can weather rapid environmental change. The authors acknowledge the method’s limits, including biases in read-count-based allele frequency estimation, and point to future refinements such as quantitative MiSeq with internal standards and environmental DNA haplotyping, which could sample intraspecific variation across whole landscapes without disturbing a single animal. For now, the Hiji River study stands as a proof of concept that the genetic lives of entire communities can be read from a bucket of stream insects, and that the relationship between where species live and what their genes carry is far more intricate than theory alone predicted.

Subject of Research: Haplotype-level DNA metabarcoding of aquatic insects to link species occurrence with population genetic structure and species-genetic diversity correlations at the catchment scale

Article Title: Haplotype-level metabarcoding to study the relationship between species occurrence and population genetic structure at the catchment scale in aquatic insects

Article References: Logronio, D. J. C., Chiu, M.-C., Doloiras-Laraño, A. D., Kolcsár, L.-P., Serrana, J. M., Hamamoto, K., & Watanabe, K. (2026). Haplotype-level metabarcoding to study the relationship between species occurrence and population genetic structure at the catchment scale in aquatic insects. Metabarcoding and Metagenomics, 10, Article e185031. https://doi.org/10.3897/mbmg.10.185031

Image Credits: AI Generated

DOI: 10.3897/mbmg.10.185031

Keywords: haplotype-level metabarcoding, aquatic insects, population genetics, species-genetic diversity correlations, genetic diversity, gene flow, metapopulation, metacommunity, COI barcode, Hiji River, catchment ecology, biodiversity