Pieridae is one of the largest and almost cosmopolitan groups of butterflies, which plays an important role in natural ecosystems; however, to date, its phylogeny and evolutionary history have not been fully resolved. In this study, we obtained the complete or nearly complete mitochondrial genomes of 100 pierid taxa (six newly sequenced, sixty extracted from the whole-genome data, and thirty-four directly available from GenBank). At the same time, for the first time, we conducted comparative mitogenomic and phylogenetic analyses based on these mitogenomic data, to further clarify their spatio-temporal evolutionary patterns. Comparative mitogenomic analysis showed that, except for cox2, the GC content of each of the 13 protein-coding genes (PCGs) in the rapidly diverging subfamily Pierinae was higher than in its sister group Coliadinae. Moreover, the dN/dS values of nine genes (atp6, atp8, cox1, cox3, cob, nad1, nad3, nad5, and nad6) in Pierinae were also relatively higher than those in its sister group, Coliadinae. Phylogenetic analysis showed that all the resultant phylogenetic trees were generally in agreement with those of previous studies. The Pierinae family contained six clades in total with the relationship of (Leptosiaini + (((Nepheroniini + Arthocharidini) + Teracolini) + (Pierini + Elodini))). The Pieridae originated in the Palearctic region approximately 72.3 million years ago in the late Cretaceous, and the subfamily Pierinae diverged from this family around 57.9 million years ago in the Oriental region, shortly after the K–Pg mass extinction event; in addition, the spatio-temporal evolutionary patterns of Pierinae were closely correlated with geological events and environmental changes, as well as the host plant coevolutionary scenario in Earth’s history. However, some incongruencies were observed between our results and those of previous studies in terms of shallow phylogenies for a few taxa, and should be further investigated.

Mountain ecosystems harbor evolutionarily unique and exceptionally rich biodiversity, particularly in insects. In this study, we characterized the diversity, community stability, and assembly mechanisms of butterflies on a subtropical mountain in the Chebaling National Nature Reserve, Guangdong Province, China, using grid-based monitoring across the entire region for two years. The results showed that species richness, abundance, and Faith’s phylogenetic diversity decreased with increasing elevation; taxonomic diversity played a considerable role in mediating the effects of environmental changes on stability. Moreover, our results showed that stochastic processes are dominant in governing the assembly of butterfly communities across all elevational gradients, with habitats at an elevation of 416–580 m subjected to the strongest stochastic processes, whereas heterogeneous selection processes displayed stronger effects on the assembly of butterfly communities at 744–908 m, 580–744 m, and 908–1072 m, with abiotic factors inferred as the main driving forces. In addition, significant differences were detected between the barcode tree and the placement tree for the calculated β-NTI values at 416–580 m. Overall, this study provides new insights into the effects of environmental change on the stability and assembly of butterflies in Chebaling, which will be beneficial for biodiversity conservation and policy development.