Molecular association and morphological characterisation of Himalopsyche larval types (Trichoptera, Rhyacophilidae)

DNA barcoding unveils a high diversity of caddisflies (Trichoptera) in the Mount Halimun Salak National Park (West Java; Indonesia)

Background

Trichoptera are one of the most diverse groups of freshwater insects worldwide and one of the main bioindicators for freshwater quality. However, in many areas, caddisflies remain understudied due to lack of taxonomic expertise. Meanwhile, globally increasing anthropogenic stress on freshwater streams also threatens Trichoptera diversity.

Methods

To assess the Trichoptera diversity of the area within and around the Mount Halimun Salak National Park (MHSNP or Taman Nasional Gunung Halimun Salak) in West Java (Indonesia), we conducted a molecular-morphological study on Trichoptera diversity using larvae from a benthic survey and adults from hand-netting. In addition to morphological identification, we applied four different molecular taxon delimitation approaches (Generalized Mixed Yule Coalescent, Bayesian Poisson Tree Processes, Automatic Barcode Gap Discovery and Assemble Species by Automatic Partitioning) based on DNA barcoding of Cytochrome-C-Oxidase I (COI).

Results

The molecular delimitation detected 72 to 81 Operational Taxonomic Units (OTU). Only five OTUs could be identified to species level by comparing sequences against the BOLD database using BLAST, and four more to the genus level. Adults and larvae could be successfully associated in 18 cases across six families. The high diversity of Trichoptera in this area highlights their potential as bioindicators for water quality assessment.

Conclusions

This study provides an example of how molecular approaches can benefit the exploration of hidden diversity in unexplored areas and can be a valuable tool to link life stages. However, our study also highlights the need to improve DNA barcode reference libraries of Trichoptera for the Oriental region.

The impact of sequencing depth and relatedness of the reference genome in population genomic studies: A case study with two caddisfly species (Trichoptera, Rhyacophilidae, Himalopsyche)

Whole genome sequencing for generating SNP data is increasingly used in population genetic studies. However, obtaining genomes for massive numbers of samples is still not within the budgets of many researchers. It is thus imperative to select an appropriate reference genome and sequencing depth to ensure the accuracy of the results for a specific research question, while balancing cost and feasibility. To evaluate the effect of the choice of the reference genome and sequencing depth on downstream analyses, we used five confamilial reference genomes of variable relatedness and three levels of sequencing depth (3.5×, 7.5× and 12×) in a population genomic study on two caddisfly species: Himalopsyche digitata and H. tibetana. Using these 30 datasets (five reference genomes × three depths × two target species), we estimated population genetic indices (inbreeding coefficient, nucleotide diversity, pairwise F ST, and genome-wide distribution of F ST) based on variants and population structure (PCA and admixture) based on genotype likelihood estimates. The results showed that both distantly related reference genomes and lower sequencing depth lead to degradation of resolution. In addition, choosing a more closely related reference genome may significantly remedy the defects caused by low depth. Therefore, we conclude that population genetic studies would benefit from closely related reference genomes, especially as the costs of obtaining a high-quality reference genome continue to decrease. However, to determine a cost-efficient strategy for a specific population genomic study, a trade-off between reference genome relatedness and sequencing depth can be considered.

Gene Flow and Diversification in Himalopsyche martynovi Species Complex (Trichoptera: Rhyacophilidae) in the Hengduan Mountains

The Hengduan Mountains are one of the most species-rich mountainous areas in the world. The origin and evolution of such a remarkable biodiversity are likely to be associated with geological or climatic dynamics, as well as taxon-specific biotic processes (e.g., hybridization, polyploidization, etc.). Here, we investigate the mechanisms fostering the diversification of the endemic Himalopsyche martynovi complex, a poorly known group of aquatic insects. We used multiple allelic datasets generated from 691 AHE loci to reconstruct species and RaxML phylogenetic trees. We selected the most reliable phylogenetic tree to perform network and gene flow analyses. The phylogenetic reconstructions and network analysis identified three clades, including H. epikur, H. martynovi sensu stricto and H. cf. martynovi. Himalopsyche martynovi sensu stricto and H. cf. martynovi present an intermediate morphology between H. epikur and H. viteceki, the closest known relative to the H. martynovi-complex. The gene flow analysis revealed extensive gene flow among these lineages. Our results suggest that H. viteceki and H. epikur are likely to have contributed to the evolution of H. martynovi sensu stricto and H. cf. martynovi via gene flow, and thus, our study provides insights in the diversification process of a lesser-known ecological group, and hints at the potential role of gene flow in the emergence of biological novelty in the Hengduan Mountains.