Phylogeny and historical biogeography of the subtribe Aporiina (Lepidoptera: Pieridae): implications for the origin of Australian butterflies

COI Barcodes combined with multilocus data for representative Aporia taxa shed light on speciation in the high altitude Irano-Turanian mountain plateaus (Lepidoptera: Pieridae)

Even though the high plateaus of Qinghai-Tibet and Iran share many faunal elements, the historical biogeography of the species present in this area are not very well understood. We present a complete COI barcode library for Aporia Hübner and a first comprehensive phylogeny for the genus including all known species and majority of subspecies using ten available genes (COI-COII, ND1, ND5, Cytb, EF-1a, Wg, 16S, 28S-D2/D3 and 28S-D8). We then focus on A. leucodice (Eversmann, 1843) and related taxa in order to resolve some long-standing taxonomic issues in this species-group. Based on DNA sequence data as well as morphology, we raise Aporia illumina (Grum-Grshimailo 1890) stat. nov. (= pseudoillumina Tshikolovets 2021 syn. nov.) as a distinct species and designate a lectotype; synonymize Aporia leucodice leucodice Eversmann, 1843 (= A. l. morosevitshae Sheljuzhko, 1908 syn. nov.); and describe a new species, Aporia ahura sp. nov., from the Central Alborz Mountains in northern Iran.

To and fro in the archipelago: Repeated inter-island dispersal and New Guinea's orogeny affect diversification of Delias, the world's largest butterfly genus

The world's largest butterfly genus Delias, commonly known as Jezebels, comprises ca. 251 species found throughout Asia, Australia, and Melanesia. Most species are endemic to islands in the Indo-Australian Archipelago or to New Guinea and nearby islands in Melanesia, and many species are restricted to montane habitats over 1200 m. We inferred an extensively sampled and well-supported molecular phylogeny of the group to better understand the spatial and temporal dimensions of its diversification. The remarkable diversity of Delias evolved in just ca. 15-16 Myr (crown age). The most recent common ancestor of a clade with most of the species dispersed out of New Guinea ca. 14 Mya, but at least six subsequently diverging lineages dispersed back to the island. Diversification was associated with frequent dispersal of lineages among the islands of the Indo-Australian Archipelago, and the divergence of sister taxa on a single landmass was rare and occurred only on the largest islands, most notably on New Guinea. We conclude that frequent inter-island dispersal during the Neogene-likely facilitated by frequent sea level change-sparked much diversification during that period. Many extant New Guinea lineages started diversifying 5 Mya, suggesting that orogeny facilitated their diversification. Our results largely agree with the most recently proposed species group classification system, and we use our large taxon sample to extend this system to all described species. Finally, we summarize recent insights to speculate how wing pattern evolution, mimicry, and sexual selection might also contribute to these butterflies' rapid speciation and diversification.

Comprehensive phylogeny of Pieridae butterflies reveals strong correlation between diversification and temperature

Temperature is thought to be a key factor influencing global species richness patterns. We investigate the link between temperature and diversification in the butterfly family Pieridae by combining next generation DNA sequences and published molecular data with fine-grained distribution data. We sampled nearly 600 pierid butterfly species to infer the most comprehensive molecular phylogeny of the family and curated a distribution dataset of more than 800,000 occurrences. We found strong evidence that species in environments with more stable daily temperatures or cooler maximum temperatures in the warm seasons have higher speciation rates. Furthermore, speciation and extinction rates decreased in tandem with global temperatures through geological time, resulting in a constant net diversification.

Spatio-Temporal Evolutionary Patterns of the Pieridae Butterflies (Lepidoptera: Papilionoidea) Inferred from Mitogenomic Data

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.

Molecular phylogeny, classification, biogeography and diversification patterns of a diverse group of moths (Geometridae: Boarmiini)

Understanding how and why some groups have become more species-rich than others, and how past biogeography may have shaped their current distribution, are questions that evolutionary biologists have long attempted to answer. We investigated diversification patterns and historical biogeography of a hyperdiverse lineage of Lepidoptera, the geometrid moths, by studying its most species-rich tribe Boarmiini, which comprises ca. 200 genera and ca. known 3000 species. We inferred the evolutionary relationships of Boarmiini based on a dataset of 346 taxa, with up to eight genetic markers under a maximum likelihood approach. The monophyly of Boarmiini is strongly supported. However, the phylogenetic position of many taxa does not agree with current taxonomy, although the monophyly of most major genera within the tribe is supported after minor adjustments. Three genera are synonymized, one new combination is proposed, and four species are placed in incertae sedis within Boarmiini. Our results support the idea of a rapid initial diversification of Boarmiini, which also implies that no major taxonomic subdivisions of the group can currently be proposed. A time-calibrated tree and biogeographical analyses suggest that boarmiines appeared in Laurasia ca. 52 Mya, followed by dispersal events throughout the Australasian, African and Neotropical regions. Most of the transcontinental dispersal events occurred in the Eocene, a period of intense geological activity and rapid climate change. Diversification analyses showed a relatively constant diversification rate for all Boarmiini, except in one clade containing the species-rich genus Cleora. The present work represents a substantial contribution towards understanding the evolutionary origin of Boarmiini moths. Our results, inevitably biased by taxon sampling, highlight the difficulties with working on species-rich groups that have not received much attention outside of Europe. Specifically, poor knowledge of the natural history of geometrids (particularly in tropical clades) limits our ability to identify key innovations underlying the diversification of boarmiines.

Molecular phylogeny and historical biogeography of Parnara butterflies (Lepidoptera: Hesperiidae)

The butterfly genus Parnara (Hesperiinae: Baorini), of which some are major pests of economic crops (e.g., rice, wild rice stems and sugarcane), currently consists of 10 species and several subspecies and has a highly disjunct distribution in Australia, Africa, and Asia. We determined the systematic relationships and biogeographical history of the genus by reconstructing the phylogeny based on eight genes and 101 specimens representing all 10 recognized species. Four species delimitation methods (ABGD, bPTP, GMYC and BPP) were also employed to assess the taxonomic status of each species. Based on these results and analyses, we recognize 11 extant species in the genus. The status of the taxon P. naso poutieri (Boisduval, 1833) from Madagascar is revised as a distinct species, Parnara poutieri (Boisduval, 1833) stat. rev. The subspecies P. guttata mangala (Moore, 1866) syn. nov. is synonymized with P. guttata guttata (Bremer & Grey, 1853), while P. bada (Moore, 1878) is provisionally treated as a complex of two species, namely P. bada and P. apostata (Snellen, 1886). The monophyly of Parnara is strongly supported, with the following relationships: P. amalia + ((P. monasi + (P. poutieri + P. naso)) + ((P. kawazoei + P. bada complex) + (P. ganga + (P. ogasawarensis + (P. guttata + P. batta))))). Divergence time and ancestral range estimates indicate that the common ancestor of Parnara originated in an implausible area of Australia, Africa, and Oriental region in the mid-Oligocene and then differentiated in the late Miocene-late Pliocene. Dispersal and range expansion have played an important role in diversification of the genus in Asia and Afica. Relatively stable geotectonic plates at the time when most extant lineages appeared during the late Miocene-early Pliocene might have been the factor responsible for the relatively constant low dynamic rate of diversification within the group.

Evaluating an Alleged Mimic of the Monarch Butterfly: Neophasia (Lepidoptera: Pieridae) Butterflies are Palatable to Avian Predators

Some taxa have adopted the strategy of mimicry to protect themselves from predation. Butterflies are some of the best representatives used to study mimicry, with the monarch butterfly, Danaus plexippus (Lepidoptera: Nymphalidae) a well-known model. We are the first to empirically investigate a proposed mimic of the monarch butterfly: Neophasia terlooii, the Mexican pine white butterfly (Lepidoptera: Pieridae). We used captive birds to assess the palatability of N. terlooii and its sister species, N. menapia, to determine the mimicry category that would best fit this system. The birds readily consumed both species of Neophasia and a palatable control species but refused to eat unpalatable butterflies such as D. plexippus and Heliconius charithonia (Lepidoptera: Nymphalidae). Given some evidence for mild unpalatability of Neophasia, we discuss the results considering modifications to classic mimicry theory, i.e., a palatability-based continuum between Batesian and Müllerian mimicry, with a quasi-Batesian intermediate. Understanding the ecology of Neophasia in light of contemporary and historical sympatry with D. plexippus could shed light on the biogeography of, evolution of, and predation pressure on the monarch butterfly, whose migration event has become a conservation priority.

Behavioral Evidence for Host Transitions in Plant, Plant Parasite, and Insect Interactions

Specialized herbivorous insects have the ability to transition between host plant taxa, and considering the co-evolutionary history between plants and the organisms utilizing them is important to understanding plant insect interactions. We investigated the role of a pine tree parasite, dwarf mistletoe (Arceuthobium spp.) M. Bieb. Santalales: Viscaceae, in mediating interactions between Neophasia (Lepidoptera: Pieridae) butterflies and pine trees, the butterflies' larval hosts. Mistletoe is considered the butterflies' ancestral host, and the evolutionary transition to pine may have occurred recently. In Arizona, United States, we studied six sites in pine forest habitats: three in Neophasia menapia (Felder and R. Felder, 1859) habitat and three in Neophasia terlooii Behr, 1869 habitat. Each site contained six stands of trees that varied in mistletoe infection severity. Butterfly behavior was observed and ranked at each stand. Volatile compounds were collected from trees at each site and analyzed using gas chromatography-mass spectroscopy. Female butterflies landed on or patrolled around pine trees (i.e., interacted) more than males, and N. terlooii interacted more with pine trees than N. menapia. Both butterfly species interacted more with tree stands harboring greater mistletoe infection, and N. terlooii interacted more with heavily infected tree stands than did N. menapia. The influence of mistletoe on Neophasia behavior may be mediated by differences in tree volatiles resulting from mistletoe infection. Volatile profiles significantly differed between infected and uninfected pine trees. The role of mistletoe in mediating butterfly interactions with pines has implications for conservation biology and forest management, and highlights the importance of understanding an organism's niche in an evolutionary context.

Differences in winter cold hardiness reflect the geographic range disjunction of Neophasia menapia and Neophasia terlooii (Lepidoptera: Pieridae)

Predicting how rapid climate change will affect terrestrial biota depends on a thorough understanding of an organism's biology and evolutionary history. Organisms at their range boundaries are particularly sensitive to climate change. As predominantly terrestrial poikilotherms, insects are often geographically limited by extremes in ambient temperatures. We compared the cold hardiness strategies of two geographically widespread butterflies, the pine white, Neophasia menapia, and the Mexican pine white, N. terlooii (Lepidoptera: Pieridae), at the near-contact zone of their range boundaries. Eggs are laid on pine needles and are exposed to harsh winter conditions. Eggs were collected from wild-caught butterflies, and we determined the supercooling point (SCP) and lower lethal temperature (LLT₅₀) of overwintering eggs. The SCP of Neophasia menapia eggs (-29.0 ± 0.6 °C) was significantly lower than that of N. terlooii eggs (-21.8 ± 0.7 °C). Both species were freeze-intolerant and capable of surviving down to their respective SCPs (LLT₅₀ of N. menapia between -30 and -31 °C, N. terlooii between -20 and -21 °C). Cold exposure time did not affect the survival of N. menapia, but N. terlooii experienced somewhat greater mortality at sub-freezing temperatures during longer exposures. Our results, coupled with an analysis of microclimate data, indicate that colder winters in northern Arizona may contribute to the northern range limit for N. terlooii. Furthermore, careful analysis of historical weather data indicates that mortality from freezing is unlikely in southern Arizona but possible in northern Arizona. Movements of Neophasia range boundaries could be monitored as potential biological responses to climate change.

Systematics and biogeography of the Gondwanan Orthocladiinae (Diptera: Chironomidae)

Restrictions to effective dispersal and gene flow caused by the fragmentation of ancient supercontinents are considered to have driven diversification and speciation on disjunct landmasses globally. Investigating the role that these processes have played in the development of diversity within and among taxa is crucial to understanding the origins and evolution of regional biotas. Within the chironomid (non-biting midge) subfamily Orthocladiinae (Diptera: Chironomidae), a group of genera that are distributed across the austral continents (Australia, New Zealand, South America) have been proposed to represent a relict Gondwanan clade. We used a molecular approach to resolve relationships among taxa with the aim to determine the relative roles that vicariance and dispersal may have played in the evolution of this group. Continental biotas did not form monophyletic groups, in accordance with expectations given existing morphological evidence. Patterns of phylogenetic relationships among taxa did not accord with expected patterns based on the geological sequence of break-up of the Gondwanan supercontinent. Likewise, divergence time estimates, particularly for New Zealand taxa, largely post-dated continental fragmentation and implied instead that several transoceanic dispersal events may have occurred post-vicariance. Passive dispersal of gravid female chironomid adults is the most likely mechanism for transoceanic movement, potentially facilitated by West Wind Drift or anti-cyclone fronts. Estimated timings of divergence among Australian and South American Botryocladius, on the other hand, were congruent with the proposed ages of separation of the two continents from Antarctica. Taken together, these data suggest that a complex relationship between both vicariance and dispersal may explain the evolution of this group. The sampling regime we implemented here was the most intensive yet performed for austral members of the Orthocladiinae and unsurprisingly revealed several novel taxa that will require formal description.

A phylogeny of the oil bee tribe Ctenoplectrini (Hymenoptera: Anthophila) based on mitochondrial and nuclear data: evidence for early Eocene divergence and repeated out-of-Africa dispersal

The bee tribe Ctenoplectrini, with two genera, comprises nine species in tropical Africa and ten in Asia and Australia. Most of them collect floral oil, pollen, and nectar from Cucurbitaceae, but three species are thought to be cleptoparasites. The unusual morphology of Ctenoplectrini has made it difficult to infer their closest relatives, in turn preventing an understanding of these bees' geographic and temporal origin. We used two mitochondrial and two nuclear markers (4741 nucleotides) generated for most of the species to test the monophyly of the tribe, its relationships to other Apidae, and its biogeographic history. Ctenoplectrini are strongly supported as monophyletic and closest to the Long-horned bees, Eucerini. The presumably cleptoparasitic species form a clade (Ctenoplectrina) that is sister to the remaining species (Ctenoplectra), confirming the independent evolution of cleptoparasitism in this tribe. Tree topology and molecular dating together suggest that Ctenoplectrini originated in Africa in the Early Eocene and that Ctenoplectra dispersed twice from Africa to Asia, sometime in the Late Eocene, 30-40 my ago, from where one species reached the Australian continent via Indonesia and New Guinea in the mid-Miocene, c. 13 my ago. Dry and cool mid-Miocene climates also coincide with the divergence between Ctenoplectra bequaerti from West Africa and Ctenoplectra terminalis from East and South Africa, perhaps related to fragmentation of the equatorial African rainforest belt.

Biogeography of butterflies in the Australian monsoon tropics

The biogeography of butterflies within the monsoon tropical biome of northern Australia is reviewed in terms of patterns of species richness, endemism and area relationships. Available data indicate that the region supports a relatively rich fauna, comprising 265 species (~62% of the total Australian fauna), but endemism is low (6%). No genera are endemic to the monsoon tropics, but two (Neohesperilla, Nesolycaena) are characteristic components, embracing a total of seven species in the region, of which five are endemic. Three ecological specialists (Neohesperilla senta, Elodina walkeri, Candalides delospila), each associated with different vegetation types, appear to be characteristic elements of the monsoon tropics. Of 67 range-restricted species in the monsoon tropics, 15 (mostly associated with savanna) are endemic to the region, while 52 (mostly associated with rainforest) are non-endemic, occurring also in south-east Asia and/or mainland New Guinea. A pronounced attenuation in species richness from Cape York Peninsula across the Top End to the Kimberley is evident. Within the monsoon tropics, Cape York Peninsula stands out as an area of exceptional biodiversity, with 95% of the butterflies (251 species; 7 endemic species, 31 endemic subspecies/geographical forms) recorded from the entire region, compared with the Top End (123 species; 3 endemic species, 17 endemic subspecies/geographical forms). In contrast, the Kimberley has a comparatively depauperate fauna (85 species; 1 endemic species, 0 endemic subspecies) without strong Indonesian affinities, and contains only two range-restricted species. A sister-area relationship between Cape York Peninsula and the Top End–Kimberley is evident in one clade, Acrodipsas hirtipes (northern Cape York Peninsula) + A. decima (Top End), with a pairwise divergence of ~1% based on mtDNA, and is suspected in another, Nesolycaena medicea (southern Cape York Peninsula) and N. urumelia (Top End) + N. caesia (Kimberley); a further five species show similar sister-area relationships across the Carpentarian Gap but at the level of subspecies or geographical form. Three general and complementary hypotheses are proposed to explain patterns of geographical differentiation of butterflies in the monsoon tropics: (1) the Carpentarian Gap is a biogeographical filter, functioning as a barrier for some species but as a bridge for others; (2) divergence among taxa between Cape York Peninsula and the Top End–Kimberley has occurred fairly recently (Quaternary), probably through vicariance; and (3) the Bonaparte Gap, with the exception of Nesolycaena, is not a vicariant barrier for butterflies in the Top End and Kimberley.