The early evolutionary pattern and molecular adaptation mechanism of alpine Parnassius butterflies to high altitudes in Qinghai–Tibet Plateau are poorly understood up to now, due to difficulties in sampling, limited sequence data, and time calibration issues. Here, we present large-scale transcriptomic datasets of eight representative Parnassius species to reveal the phylogenetic timescale and potential genetic basis for high-altitude adaptation with multiple analytic strategies using 476 orthologous genes. Our phylogenetic results strongly supported that the subgenus Parnassius formed a well-resolved basal clade, and the subgenera Tadumia and Kailasius were closely related in the phylogenetic trees. In addition, molecular dating analyses showed that the Parnassius began to diverge at about 13.0 to 14.3 million years ago (middle Miocene), correlated with their hostplant’s spatiotemporal distributions, as well as geological and palaeoenvironmental changes of the Qinghai–Tibet Plateau. Moreover, the accelerated evolutionary rate, candidate positively selected genes and their potentially functional changes were detected, probably contributed to the high-altitude adaptation of Parnassius species. Overall, our study provided some new insights into the spatiotemporally evolutionary pattern and high altitude adaptation of Parnassius butterflies from the extrinsic and intrinsic view, which will help to address evolution, biodiversity, and conservation questions concerning Parnassius and other butterfly species.

Parnassius glacialis is a typical “Out of the QTP” alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for “Out of the QTP” of P. glacialis.

The origin and dispersal history of the large butterfly subfamily Nymphalinae are not fully understood, due to internal phylogenetic and time calibration issues. We conducted phylogenetic and dating analyses using mitochondrial and nuclear genes of biogeographically diverse groups of the Nymphalinae in order to resolve some controversial relationships and the paleobiogeographic pattern of the subfamily. Our results support the sister relationship of Vanessa (Tribe Nymphalini) and the Nymphalis-group, and the grouping of the three old-world genera (Rhinopalpa, Kallimoides and Vanessula) within Tribe Victorinini. Molecular dating analyses invoking two additional calibrations under the butterfly-host plant coevolutionary scenarios result in a relatively deeper divergence of the subfamily’s two major clades (Nymphalini and the Kallimoids), compatible with the Cretaceous floral turnover scenario during the so-called Cretaceous Terrestrial Revolution. Phylobiogeographic analyses reveal that the Oriental region is probably the center of early divergences for Nymphalinae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by repeated dispersals into the rest of the Old World and the New World during various periods beginning in Eocene. The biogeographic history indicates that temperature changes and host-plant diversification may have facilitated the dispersals of this butterfly subfamily, with accelerated global colonization during the middle to late Miocene.

We studied 239 imagoes of 12 Parnassius species collected from the mountains of the Qinghai–Tibet Plateau (QTP) and its neighbouring areas in China. We selected three mitochondrial gene (COI, ND1, and ND5) sequences, along with the homologous gene sequences of other Parnassius species from GenBank, to reconstruct the phylogenetic tree and biogeographic history of this genus. Our results show that Parnassius comprises eight monophyletic subgenera, with subgenus Parnassius at the basal position; the genus crown group originated during the Middle Miocene (ca. 16.99 Ma), and species diversification continued during sustained cooling phases after the Middle Miocene Climate Optimum (MMCO) when the QTP and its neighbouring regions experienced rapid uplift and extensive orogeny. A phylogenetic network analysis based on transcriptomes from GenBank suggests that ancient gene introgression might have contributed to the spread of the Parnassius genus to different altitudes. Ancestral area reconstruction indicates that Parnassius most likely originated in West China (QTP and Xinjiang) and then spread to America in two dispersal events as subgenera Driopa and Parnassius, along with their host plants Papaveraceae and Crassulaceae, respectively. Our study suggests that extensive mountain-building processes led to habitat fragmentation in the QTP, leading to the early diversification of Parnassius, and climate cooling after MMCO was the driving mechanism for the dispersal of Parnassius butterflies from West China to East Asia, Europe, and North America.

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.

Bees (Hymenoptera, Apoidea and Anthophila) are distributed worldwide and considered the primary pollinators of angiosperm. Megachilidae is one of the largest families of Anthophila. In this study, two complete mitogenomes of cuckoo bees in Megachilidae, namely Coelioxys fenestrata and Euaspis polynesia, were amplified and sequenced, with a length of 17,004 bp (C. fenestrata) and 17,682 bp (E. polynesia). The obtained results show that 37 mitogenomic genes and one putative control region were conserved within Hymenoptera. Truncated stop codon T was found in the cox3 gene of E. polynesia. The secondary structure of small (rrnS) and large (rrnL) rRNA subunits contained three domains (28 helices) and five domains (44 helices) conserved within Hymenoptera, respectively. Compared with ancestral gene order, gene rearrangement events included local inversion and gene shuffling. In order to reveal the phylogenetic position of cuckoo bees, we performed phylogenetic analysis. The results supported that all families of Anthophila were monophyletic, the tribe-level relationship of Megachilidae was Osmiini + (Anthidiini + Megachilini) and Coelioxys fenestrata was clustered to the Megachile genus, which was more closely related to Megachile sculpturalis and Megachile strupigera than Euaspis polynesia.