Mitochondrial phylogeography of grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) endemic to the Qinghai–Tibetan plateau

Abstract

Abstract Grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) are the most damaging pests to alpine meadows in the Qinghai–Tibetan Plateau (QTP). Here, we conducted extensive sampling from 39 geographic populations covering almost the entire distribution of the eight QTP Gynaephora (Hübner) species to investigate phylogeographic patterns and speciation based on two mitochondrial genes (COI and ND5). A total of 40 haplotypes were detected in the 39 populations, with >70% of all haplotypes not shared between populations. The monophyletic QTP Gynaephora migrated from non‐QTP regions during the Pliocene, corresponding to the uplift of the QTP, suggesting a mode of transport into the QTP. Among the eight QTP Gynaephora species described by morphological characteristics, two species (G. alpherakii and G. menyuanensis) were recovered as monophyletic groups (Clades B and C), while the remaining six formed two monophyletic clades: Clade A (G. qinghaiensis, G. jiuzhiensis, and G. qumalaiensis) and Clade D (G. aureata, G. ruoergensis, and G. minora). These results suggested that the number of the QTP Gynaephora species may be overestimated and further studies based on both morphological and nuclear gene data are needed. Genetic differentiation and speciation of the QTP Gynaephora were likely driven by the QTP uplifts and associated climate fluctuations during the Pleistocene, indicated by divergence time estimation, suggesting that isolation and subsequent divergence was the dominant mode of speciation. The Sanjiangyuan region (i.e., Clade A, characterized by high genetic diversity) may have been a glacial refugium of the QTP Gynaephora, as supported by analyses of gene flow and biogeography. High levels of genetic diversity were found in QTP Gynaephora, without population expansion, which may explain the high‐altitude adaptation and outbreaks of grassland caterpillars in alpine meadows of the QTP. This study provides the largest phylogeographic analysis of QTP Gynaephora and improves our understanding of the diversity and speciation of QTP insects.

Keywords

• biogeography
• divergence time estimation
• gene flow
• insects
• moths
• population genetic structure