Delimiting species, revealing cryptic diversity in Molytinae (Coleoptera: Curculionidae) weevil through DNA barcoding

The subfamily Molytinae (Coleoptera: Curculionidae), being the second largest group within the family Curculionidae, exhibits a diverse range of hosts and poses a serious threat to agricultural and forestry industries. We used 1,290 cytochrome c oxidase subunit I (COI) barcodes to assess the efficiency of COI barcodes in species differentiation and uncover cryptic species diversity within weevils of Molytinae. The average Kimura 2-parameter distances within species, genus, and subfamily were 2.90%, 11.0%, and 22.26%, respectively, indicating significant genetic differentiation at both levels. Moreover, there exists a considerable degree of overlap between intraspecific (0%-27.50%) and interspecific genetic distances (GDs; 0%-39.30%). The application of Automatic barcode gap discovery, Assemble Species by Automatic Partitioning, Barcode Index Number, Poisson Tree Processes (PTP), Bayesian Poisson Tree Processes (bPTP), and jMOTU resulted in the identification of 279, 275, 494, 322, 320, and 279 molecular operational taxonomic units, respectively. The integration of 6 methods successfully delimited species of Molytinae in 86.6% of all examined morphospecies, surpassing a threshold value of 3% GD (73.0%). A total of 28 morphospecies exhibiting significant intraspecific divergences were assigned to multiple MOTUs, respectively, suggesting the presence of cryptic diversity or population divergence. The identification of cryptic species within certain morphological species in this study necessitates further investigation through comprehensive taxonomic practices in the future.